SNS 공유 미리보기
좌굴 고유값은 안전율일까? 초기 결함과 비선형 내하력 읽기
좌굴 고유값의 하중 배율을 실제 내하력으로 오해하지 않으려면 무엇을 확인해야 할까요? Euler 계산, 초기 결함의 진폭, 비선형 평형 경로를 구분해 정리합니다.
카카오톡, 페이스북 등 SNS 공유 시
위와 같은 형태로 노출됩니다.
좌굴 고유값은 안전율일까? 초기 결함과 비선형 내하력 읽기
좌굴 해석 화면에 첫 번째 고유값이 4.2로 나오면 ‘설계 하중의 4.2배까지 안전하다’고 읽고 싶어질 수 있어요. 하지만 그 숫자가 어떤 하중 패턴을 곱하는지, 이미 들어 있던 하중은 무엇인지, 완벽한 형상과 탄성 재료를 가정했는지부터 알아야 합니다. 보기 좋은 좌굴 모드 하나가 실제 붕괴 하중을 대신하지는 않죠.
이번 글은 2026년 10월 5~11일 신규 발표를 소개하는 뉴스가 아니라 공식 해석 문서로 정리한 실무 해설입니다. 앞선 FE 검증 글의 일반 절차를 반복하기보다, 좌굴 고유값에서 비선형 내하력 검토로 넘어갈 때 생기는 해석의 틈을 좁혀볼게요. 수치 예제는 직접 계산한 가상 사례이며, 실물 시험이나 저자의 설계 승인 결과가 아닙니다.
STEP 1고유값은 어떤 하중을 몇 배 하는 숫자인가요
선형 고유값 좌굴은 기준 상태 주변에서 구조의 강성이 특이해지는 하중 조합을 추정합니다. 하중을 조금 더했을 때 평형 형상이 달라질 수 있는 방향을 모드로 보여주는 거예요. 실제 구조를 하중 증가에 따라 끝까지 눌러본 결과가 아니므로, 소성 변형이나 접촉의 변화를 자동으로 포함한 붕괴 시험처럼 읽으면 안 됩니다.
우선 하중 장부를 나눠보세요. 기준 상태를 만든 선행 하중을 P, 좌굴 단계에서 배율을 찾는 추가 하중 패턴을 Q라고 두면, 해당 문서의 예측 하중 조합은 P + λQ로 표현됩니다. λ는 Q에 곱하는 고유값이에요. P까지 함께 배율을 곱하는지 여부는 해석 단계의 정의에 달렸으므로, 프로그램의 결과 이름만 보고 전체 하중의 배수라고 결론 내리지 않아야 합니다.
예측 좌굴 하중 조합 = P + λQ
예시: 같은 방향의 축압축 선행 하중 P = 20 kN,
추가 기준 하중 Q = 10 kN, 계산되었다고 가정한 λ = 4.2
예측 총 압축하중 = 20 + 4.2 × 10 = 62 kN
배율은 추가 하중 Q에만 곱해요
회색: 유지되는 P = 20 kN · 청록색: λQ = 4.2 × 10 = 42 kN
빗금: 선행 하중까지 4.2배로 늘린 잘못된 해석
여기서 4.2는 산술 설명을 위해 가정한 고유값이고, 실제 FE를 실행해 얻은 값은 아닙니다. 초기 하중 20 kN까지 묶어 4.2 × 30 = 126 kN으로 읽으면 전혀 다른 결과가 됩니다. 더구나 P와 Q가 압력, 축력, 온도 하중처럼 서로 다른 공간 분포를 갖는다면 숫자 하나로 더할 수 없어요. 각 하중 성분의 조합과 작용 위치를 유지해야 합니다.
고유벡터의 변형 크기도 실제 변위가 아닙니다. 모드 그림의 크게 휜 형상은 모드를 알아보기 위한 정규화와 화면 확대가 반영된 표현일 수 있어요. 그 그림에서 최대 변위가 특정 길이로 보인다고 실제 구조가 그만큼 휜 뒤 좌굴한다는 뜻은 아닙니다. 모드는 위험한 변형 패턴을 찾는 자료로 쓰고, 변위와 하중의 실제 관계는 별도로 구해야 해요.
이 구분만 해도 보고서의 첫 문장이 달라집니다. ‘안전율 4.2’ 대신 ‘정의한 기준 상태와 추가 하중 패턴에 대한 첫 고유값 4.2’라고 적고, 어떤 가정 아래의 값인지 붙일 수 있죠. 설계 기준에서 정한 저항 계수나 하중 계수까지 검토하지 않은 숫자에 안전율이라는 이름을 붙이지 않는 것이 출발점입니다.
STEP 2간단한 기둥 계산으로 단위와 경계조건을 확인해요
복잡한 모델에 앞서 이상적인 핀–핀 기둥을 생각해볼게요. 길이 L = 2,000 mm, 탄성계수 E = 200,000 N/mm², 직사각형 단면은 폭 b = 40 mm, 약축 휨에 쓰이는 두께 h = 10 mm라고 가정합니다. 재료는 선형 탄성이고, 기둥은 곧고, 압축력은 중심을 지나며, 양 끝은 이상적인 핀으로 지지된 조건이에요.
단면적 A = bh = 400 mm²
약축 단면2차모멘트 I = bh³/12 = 3,333.33 mm⁴
핀–핀 유효 길이 계수 K = 1
Euler 하중 Pcr = π²EI/(KL)² ≈ 1,644.93 N = 1.645 kN
평균 압축응력 Pcr/A ≈ 4.112 N/mm² = 4.112 MPa
단위는 N/mm² × mm⁴ ÷ mm²이므로 최종적으로 N이 됩니다. 이 예는 매우 세장한 얇은 기둥을 택해 이상적인 탄성 좌굴 계산을 보여주는 것입니다. 단면의 약축과 강축을 바꾸면 I가 달라지고, 같은 형상이라도 횡지지와 단부 회전 구속이 달라지면 유효 길이가 달라져요. 특히 두께 h가 세제곱으로 들어가므로 단면 방향을 잘못 읽은 오류가 크게 드러납니다.
같은 이상화에서 K가 1에서 2가 되면 식의 분모가 네 배가 되어 계산 하중은 약 0.411 kN이 됩니다. 이는 유효 길이 가정의 민감도를 보여주는 비교예요. 실제 접합부가 핀인지 고정인지 결정해주는 값은 아닙니다. 볼트 연결이나 용접부를 그림만 보고 완전 고정으로 놓기보다 회전 구속의 근거와 주변 구조의 유연성을 확인해야 합니다.
이 계산과 단순 FE 모델이 맞는지는 단위·단면·구속을 확인하는 작은 시험으로 쓸 수 있어요. 하지만 계산이 맞았다고 쉘 패널이나 보강판의 실제 내하력까지 검증된 것은 아닙니다. 국부좌굴, 초기 휨, 잔류응력, 편심과 재료 비선형이 들어오면 이 간단한 모델 밖으로 나갑니다. 확인한 범위를 작게 정확히 써두는 것이 더 유용하죠.
짧고 두꺼운 기둥에서는 Euler 식이 예측한 응력에 이르기 전에 재료가 비선형 영역에 들어갈 수 있습니다. 따라서 실제 재료의 응력–변형률 관계와 세장비 적용 범위를 확인해야 해요. 여기서는 항복강도나 허용응력을 임의로 넣어 합격 여부를 만들지 않았습니다. 손계산의 역할은 정답 도장을 찍는 것이 아니라 해석 가정을 드러내는 데 있어요.
STEP 3초기 결함은 보기 좋은 모드 한 개로 끝나지 않아요
완벽하게 곧고 대칭인 모델은 실제 제작물보다 높은 대칭성을 갖습니다. 좌굴 이후의 경로를 찾을 때는 측정된 형상 편차나 합리적으로 정의한 초기 결함을 모델에 넣어야 할 수 있어요. 고유값 해석에서 얻은 모드를 그 형상 후보로 쓸 수 있지만, 첫 모드가 언제나 실제 결함이거나 가장 불리한 조합이라고 단정할 수는 없습니다.
모드를 이용할 때는 진폭의 정의부터 분리하세요. 예를 들어 최대 법선 방향 변위가 1이 되도록 무차원 정규화한 모드 형상 φ를 쓰고, 결함 진폭 a를 mm로 정했다면 초기 형상 변화는 aφ입니다. 여러 모드를 더할 때는 각 모드의 계수뿐 아니라 최종 형상에서의 최대 편차와 위치도 확인해야 해요. 서로 다른 부위의 변형이 더해지거나 상쇄될 수 있기 때문입니다.
초기 형상: ximperfect = xideal + Σ aiφi
φi: 정규화 기준을 명시한 무차원 형상
ai: 그 정규화에 대응하는 길이 계수 [mm]
가령 두께 t = 4 mm인 판의 계산 연습에서 0.1t와 0.5t를 비교하면 진폭은 0.4 mm와 2.0 mm입니다. 이 값은 민감도 시험의 설명용 예시일 뿐 제작 공차나 규격상의 초기 결함값이 아니에요. 실제 적용값은 측정 형상, 해당 설계 기준, 제작 허용차와 검토 목적에서 근거를 찾아야 합니다. 임의로 작은 결함만 넣고 결과가 좋다는 이유로 채택하면 검토 범위가 숨겨집니다.
첫 몇 개 고유값이 가까우면 서로 다른 모드가 영향을 줄 수 있어요. 전체 휨과 국부 판좌굴이 경쟁하는 구조에서는 한 모드와 한 부호만 계산한 결과가 충분한지 따져봐야 합니다. 후보의 선택 이유, 진폭 범위, 부호와 조합을 기록하고, 어떤 결과를 설계 검토에 사용했는지 추적할 수 있게 남겨두세요.
형상 파일을 옮길 때 생기는 단순 실수도 큽니다. 노드 번호가 바뀌었거나 좌표계가 다르면 의도한 위치에 결함이 들어가지 않을 수 있어요. Abaqus 문서는 초기 결함을 가져올 때 노드 레이블을 사용하며 두 모델의 호환성을 자동으로 확인하지 않는다고 설명합니다. 입력 완료 메시지보다 변형 전후 좌표와 실제 편차 분포를 확인하는 일이 중요하죠.
STEP 4비선형 해석에서는 하중–변위 경로를 읽어요
비선형 검토의 질문은 ‘어떤 모드가 가능한가’에서 ‘정의한 형상과 재료가 하중을 따라 어떻게 변하는가’로 이동합니다. 큰 변형에 따른 기하학적 변화, 실제에 필요한 재료 비선형, 접촉과 경계조건의 변화를 포함할지 결정해야 해요. 모든 옵션을 켜는 것보다 내하력을 지배할 물리를 빠뜨리지 않는 것이 중요합니다.
일반적인 하중 제어는 한계점 부근에서 진행하기 어려울 수 있어요. 변위 제어가 적절한 경로도 있지만, 하중과 변위가 함께 되돌아가는 경우에는 다른 추적 방법이 필요할 수 있습니다. Riks 같은 호장법은 하중 계수와 변위를 함께 미지수로 두고 평형 경로를 추적하는 방법입니다. 곡선의 내려가는 부분을 구했다는 사실이 그 상태를 실제 장비가 안정적으로 유지한다는 뜻은 아니에요.
Riks 방법이 모든 실패의 해결책도 아닙니다. 하중을 하나의 매개변수에 따라 비례적으로 변화시킬 수 있는지, 접촉 분리나 국부 불안정이 방법의 가정과 맞는지 확인해야 해요. 충격이나 빠른 동적 좌굴을 준정적 평형 경로만으로 설명할 수도 없습니다. 문제의 시간 척도와 실제 하중 경로에 맞는 절차를 골라야 합니다.
출력은 변형 그림보다 하중–변위 관계를 먼저 준비해두면 좋아요. 예를 들어 하중 작용점의 축방향 변위, 대표 위치의 횡변위, 지점 반력 합, 소성 변형의 진행과 접촉 상태를 함께 추적합니다. 하중 계수만 저장하면 고정 하중이 있는 경우 총 하중과 혼동할 수 있으므로 실제 적용 하중 또는 반력과 대응시켜 확인하세요.
계산이 멈춘 마지막 하중을 곧바로 내하력이라고 적지는 마세요. 증분이 너무 크거나 접촉 설정이 불안정해서 수렴하지 않은 것일 수도 있습니다. 반대로 계산이 잘 수렴했다고 물리적 가정이 옳다는 보장도 없어요. 문제 구간의 잔차, 변형 모드, 에너지와 접촉 변화를 확인하고 필요한 경우 작은 증분이나 다른 제어 방법으로 경로를 점검해야 합니다.
안정화를 사용했다면 그 영향도 결과에 남겨야 합니다. 인공적인 도움 없이도 주요 경로가 유지되는지, 안정화 관련 에너지가 해석 결론을 바꾸지 않는지 확인할 필요가 있어요. 여기서 모든 구조에 통하는 한 가지 에너지 비율을 합격 기준으로 제시하지는 않습니다. 적용 해석 절차와 설계 검토의 요구에 맞는 허용 수준을 정하고 근거를 기록해야 해요.
STEP 5설계 판단으로 넘길 때 세 가지 결과를 나눠요
첫째는 이상화한 모델의 고유값과 모드입니다. 하중 패턴과 경계조건을 명시하고, 이 값으로 무엇을 선별했는지 적으세요. 둘째는 특정 초기 결함·재료·하중 경로를 넣은 비선형 응답입니다. 셋째는 적용 규격과 계수, 사용성 조건까지 포함해 정한 설계 판단이에요. 세 결과를 한 문장으로 합치면 계산과 책임의 경계가 흐려집니다.
비선형 곡선의 최고점만 중요한 것도 아닙니다. 설계 목적에 따라 과도한 횡변위, 국부 소성 변형, 접합부 한계 또는 사용성 문제가 그보다 먼저 나타날 수 있어요. 어떤 한계 상태를 내하력 판단에 쓰는지 먼저 정해두면 솔버의 마지막 프레임을 편의상 결론으로 삼는 일을 줄일 수 있습니다.
검토 결과에는 바꿔 본 조건을 함께 남겨주세요. 경계조건의 강성, 초기 결함 형상과 진폭, 재료 곡선, 하중 편심이 달라질 때 핵심 결론이 얼마나 움직였는지 확인하는 것입니다. 이때 모든 조합을 무작정 늘리기보다 물리적으로 가능한 범위와 근거를 정하고, 결론을 좌우하는 조건에 계산을 집중하는 편이 효율적이에요.
메시도 변위가 매끄럽게 보이는지만 확인하지 않습니다. 예상한 국부좌굴 파장을 표현하는지, 요소 선택과 적분 방식에 따라 모드나 경로가 크게 달라지지 않는지를 살펴야 해요. 좌굴 전 전체 강성이 비슷해도 짧은 파장의 모드를 표현하는 능력은 다를 수 있습니다. 선형 고유값의 수렴과 비선형 경로의 안정성을 같은 확인으로 대체하지 마세요.
간단한 인계 문장이라면 이렇게 쓸 수 있습니다. ‘이 값은 명시한 기준 하중에 대한 탄성 고유값 예측이다. 실제 내하력 검토에는 별도로 정한 초기 결함과 재료·기하 비선형을 사용했고, 판정은 해당 기준의 한계 상태에 따른다.’ 아직 비선형 검토를 하지 않았다면 그 부분을 완료한 것처럼 쓰지 말고 추가 확인 항목으로 남겨야 합니다.
결국 좌굴 결과를 잘 읽는다는 것은 숫자를 더 보수적으로 깎는 일이 아니에요. 어떤 가정이 어떤 값에 들어 있는지 분리하고, 실제 구조에서 중요한 가정을 근거 있게 확인하는 일입니다. 고유값은 훌륭한 출발점이지만, 실제 내하력과 설계 적합성의 최종 대답은 그 다음 검토까지 이어져야 합니다.
REFERENCES
2026-10-12 KST 확인. 아래 설명은 각 자료에서 참고한 범위이며, 본문의 가상 수치 예제는 직접 계산했습니다.
- 원문 보기
Abaqus · Eigenvalue Buckling Prediction기준 상태와 섭동 하중의 합 P + λQ, 모드의 임의 정규화, 초기 결함에 민감한 구조에서의 한계를 참고했습니다. 원문 보기 ↗
기준 상태와 섭동 하중의 합 P + λQ, 모드의 임의 정규화, 초기 결함에 민감한 구조에서의 한계를 참고했습니다.
- 원문 보기
Abaqus · Introducing a Geometric Imperfection into a Model좌굴 모드의 중첩 또는 직접 지정한 형상으로 초기 결함을 정의하는 방식과 노드 대응 확인을 참고했습니다. 원문 보기 ↗
좌굴 모드의 중첩 또는 직접 지정한 형상으로 초기 결함을 정의하는 방식과 노드 대응 확인을 참고했습니다.
- 원문 보기
Abaqus · Unstable Collapse and Postbuckling AnalysisRiks 방법의 평형 경로 추적, 비례 하중 조건과 한계를 참고했습니다. 원문 보기 ↗
Riks 방법의 평형 경로 추적, 비례 하중 조건과 한계를 참고했습니다.
- 원문 보기
Abaqus · Buckling of a Cylindrical Shell under Uniform Axial Pressure좌굴 형상에 맞는 메시와 초기 결함 민감도를 포함한 하중·변위 해석의 역할을 참고했습니다. 원문 보기 ↗
좌굴 형상에 맞는 메시와 초기 결함 민감도를 포함한 하중·변위 해석의 역할을 참고했습니다.
- 원문 보기
COMSOL · Buckling Analysis of a Truss Tower이상적인 기둥의 Euler 좌굴하중과 지지 조건에 따른 유효 길이 개념을 참고했습니다. 본문 단면과 수치는 별도 가정입니다. 원문 보기 ↗
이상적인 기둥의 Euler 좌굴하중과 지지 조건에 따른 유효 길이 개념을 참고했습니다. 본문 단면과 수치는 별도 가정입니다.
- 원문 보기
Abaqus · Eigenvalue Buckling Theory기준 상태에 더하는 탄성 섭동, 거의 선형인 좌굴 전 응답이라는 고유값 예측 전제를 참고했습니다. 원문 보기 ↗
기준 상태에 더하는 탄성 섭동, 거의 선형인 좌굴 전 응답이라는 고유값 예측 전제를 참고했습니다.
결과표의 배율은 어느 하중에 곱해지나요? 그 한 줄을 설명한 뒤, 초기 결함과 실제 하중–변위 경로가 설계 판단에 어떻게 이어지는지 확인해보세요.
