기둥은 재료가 버티는데도 왜 옆으로 무너지는가: 좌굴·불완전성·평형경로와 구조 안정성의 역학

압축시험에서 재료가 아직 탄성영역인데도 가느다란 기둥이 옆으로 크게 휘는 장면은 처음 보면 이상합니다. 단면의 평균압축응력 P/A는 낮고 재료도 부서지지 않았는데, 구조물은 더 이상 하중을 안정적으로 지탱하지 못합니다. 이 모순은 좌굴을 ‘압축강도 부족’으로 보면 풀리지 않습니다. 좌굴은 재료가 얼마나 센가보다 현재의 평형형상이 작은 교란 뒤에도 유지되는가를 묻는 문제입니다.

이 글은 기존의 공식 순서나 해석 프로그램 절차를 반복하지 않습니다. 하나의 이상적인 기둥을 두고 길이, 단면 방향, 끝단 구속, 초기 휨, 하중 편심을 한 번에 하나씩 바꾸는 사고실험을 진행합니다. 초급 독자는 그림의 ‘질문–관찰–판단’만 따라가도 되고, 숙련자는 같은 장면을 에너지의 2차 변분, 분기이론, 후좌굴 평형경로와 결함 민감도로 다시 읽을 수 있습니다.

이 글의 핵심 문장압축력은 단지 변형을 만드는 외력이 아니다. 변형된 형상에서는 압축력이 추가모멘트를 만들며 구조의 횡방향 복원강성을 낮춘다. 좌굴은 그 복원강성이 사라지거나 평형경로가 한계점에 도달하는 현상이다.

1. 응력 계산만으로는 답할 수 없는 질문

짧고 두꺼운 압축부재는 단면 전체가 압축항복하거나 재료의 압축한계에 접근하는 방식으로 내력을 잃을 수 있습니다. 이 경우 힘을 면적으로 나눈 응력이 직접적인 판단변수입니다. 그러나 길고 가는 부재에서는 작은 횡변위가 생기는 순간 축력의 작용선과 부재 중심선 사이에 거리가 생깁니다. 그 거리가 Δ라면 압축력은 추가모멘트 PΔ를 만들고, 추가모멘트는 다시 횡변위를 키웁니다. 원인과 결과가 서로 증폭하는 피드백이 형성됩니다.

그래서 좌굴에서 첫 질문은 “압축응력이 얼마인가?”가 아니라 “작게 옆으로 밀었다 놓았을 때 원래 형상으로 돌아오는가?”입니다. 돌아오면 안정평형, 새로운 위치에 머무르면 임계 또는 중립상태, 더 멀어지면 불안정평형입니다. 같은 축력과 같은 재료라도 길이, 굽힘강성, 지점조건이 바뀌면 답이 달라지는 이유가 여기에 있습니다.

MIT 16.001 좌굴 강의는 이 현상의 출발점을 변형 전 형상이 아니라 변형된 형상에서 평형을 세우는 것으로 설명합니다. 선형 정적해석에서 완전히 곧은 기둥의 횡변위가 0으로 나왔다는 사실은 안정성을 입증하지 않습니다. 그 해석은 직선 평형경로가 작은 교란에 대해 안정한지 묻지 않았을 수 있습니다.

2. 초급자의 구슬 그림과 전문가의 에너지 판정은 같은 말이다

그릇 바닥의 구슬은 조금 움직여도 다시 내려옵니다. 평평한 바닥의 구슬은 옮겨진 곳에 머물고, 산봉우리의 구슬은 작은 교란에도 멀어집니다. 이 비유를 수식으로 번역하면 전체 위치에너지 Π의 곡률을 조사하는 문제가 됩니다. 평형에서는 1차 변분이 0이고, 그 평형이 안정하려면 작은 허용변위 방향에 대한 2차 변분이 양수여야 합니다.

δΠ = 0은 “평형인가?”에 답하고, δ²Π > 0은 “그 평형이 안정한가?”에 답합니다. 임계점에서 특정 변형모드 방향의 2차 변분이 0이 되면 복원강성이 사라집니다. 임계 이후 음수가 되면 직선형상은 존재할 수 있어도 안정하지 않습니다. MIT 2.080 구조역학 Lecture 9는 이 에너지 판정을 이산 기둥모델에서 명시적으로 유도합니다.

두 층으로 읽기처음 배우는 독자는 “교란 뒤 돌아오는가”만 기억하면 충분합니다. 고급 해석에서는 같은 질문을 접선강성행렬의 고유값 또는 전체 위치에너지의 2차 변분으로 계산합니다.

3. Euler 기둥은 왜 고유치 문제가 되는가

길이 L, 일정한 굽힘강성 EI를 가진 완전한 핀–핀 기둥에 중심압축력 P가 작용한다고 합시다. 재료는 선형탄성이고 부재는 세장하며, 하중은 보존력이고 변형은 지배식을 선형화할 수 있을 만큼 작다고 가정합니다. 변형된 중심선을 v(x)라 두면 평형식은 다음 형태가 됩니다.

EI v⁗(x) + P v″(x) = 0핀 끝에서는 횡변위와 끝모멘트가 0입니다. 이 경계조건을 동시에 만족하는 비자명한 해는 임의의 P가 아니라 이산적인 하중에서만 존재합니다.

첫 번째 모드는 반 사인파이고 가장 낮은 임계하중은 Pcr = π²EI/L²입니다. 여기서 중요한 사실은 세 가지입니다. 첫째, 길이가 두 배면 기준 임계하중은 1/4이 됩니다. 둘째, 단면2차모멘트가 두 배면 임계하중은 두 배가 됩니다. 셋째, 재료의 항복강도는 이 완전탄성 기준식에 직접 등장하지 않습니다. 따라서 Euler 값은 모든 기둥의 실제 내력이 아니라, 가정이 명확한 비교 기준선입니다.

해석 화면의 모드가 50 mm 휘어 보이더라도 실제 기둥이 50 mm 변위했다는 뜻은 아닙니다. 고유벡터의 크기는 정규화 방식과 화면 배율에 따라 달라집니다. 물리적으로 읽어야 할 것은 변형의 상대적 패턴, 축, 파장과 결합된 자유도입니다.

4. 하나의 기둥을 바꾸는 변수격리 실험

복잡한 골조를 바로 해석하면 어느 변화가 결과를 만들었는지 분리하기 어렵습니다. 먼저 같은 재료와 같은 기준 기둥에서 한 변수만 바꿔 봅시다. 기준 임계하중을 1.00으로 정규화하면 인과관계가 선명해집니다.

바꾼 것직접 달라지는 항먼저 확인할 관찰그대로라고 가정한 것

길이 L임계하중은 L²에 반비례좌굴 반파장의 길이EI와 끝단 이상화

단면 방향축별 I와 회전반경 r약축 방향의 횡변형면적과 재료

끝단 회전강성유효한 모드 길이끝단 회전과 모멘트부재 자체의 EI

브레이싱 위치허용되는 모드와 파장브레이스 반력과 변위브레이스가 충분히 작동한다는 가정

초기 휨·편심분기 전 변위와 PΔ하중 초기에 이미 생기는 횡변위기준 기하와 재료모델

변수격리의 목적은 실제 구조를 단순하다고 믿는 것이 아닙니다. 복잡한 결과를 설명할 수 있는 최소한의 인과를 먼저 확보하는 것입니다. 실제 모델로 돌아간 뒤에는 바꾼 변수가 서로 결합하는지 확인해야 합니다. 예를 들어 약축 좌굴이 예상되어도 접합부의 비틀림 유연성이 크면 휨–비틀림 결합모드가 먼저 나타날 수 있습니다.

5. 유효길이는 표에서 고르는 숫자가 아니라 모드의 길이다

Pcr = π²EI/(KL)²에서 KL은 이상화된 좌굴모드의 유효길이입니다. 이를 단순한 끝단기호의 속성으로 보면 오류가 생깁니다. 실제 골조에서 기둥 끝은 완전한 핀이나 완전한 고정 사이 어딘가에 있고, 접합부 회전강성, 인접 보와 기둥의 강성, 횡이동 가능성, 브레이싱의 위치와 강성이 함께 모드를 결정합니다.

브레이스가 도면에 그려져 있다는 사실만으로 좌굴길이가 나뉘지 않습니다. 예상 모드를 막을 충분한 강성과 그 반력을 전달할 하중경로가 있어야 합니다. 반대로 접합부를 수치모델에서 완전고정하면 실제보다 끝단 회전을 억제해 임계하중을 과대평가할 수 있습니다. 유효길이는 먼저 모드형상과 변형구속을 읽고, 그 거동을 등가 Euler 기둥의 길이로 번역하는 순서가 안전합니다.

이 관점에서는 회전반경 r = √(I/A)도 단순한 단면표 값이 아닙니다. 가능한 모든 축에 대해 I와 r을 따로 보고, 실제로 허용되는 좌굴방향과 연결해야 합니다. 면적이 같아도 I가 작은 약축의 임계응력이 훨씬 낮을 수 있습니다.

6. 완전한 분기는 현실의 경로가 아니다

완전히 곧고 중심하중을 받으며 대칭인 이상모델에서는 직선 평형경로가 임계점에서 좌우 두 개의 휜 경로로 갈라지는 피치포크 분기를 보일 수 있습니다. 그러나 실제 부재에는 초기 휨, 단면 편차, 하중 편심, 잔류응력, 지점 비대칭이 존재합니다. 작은 결함은 좌우 대칭을 깨고, 임계점 전부터 한쪽으로 횡변위를 발생시킵니다.

Koiter의 1945년 박사논문 저장소 기록은 완전계의 임계상태와 초기결함에 민감한 실제 응답을 연결한 고전적 기반입니다. 중요한 결론은 “결함이 임계하중에서 일정한 값을 뺀다”가 아닙니다. 결함은 평형경로 전체를 펼치며, 결함의 크기뿐 아니라 형상, 방향과 여러 모드의 조합이 결과를 바꿀 수 있습니다.

NASA SP-8007 Rev.2가 얇은 원통 셸에서 기하학적 결함의 측정과 제작공정별 결함 서명을 강조하는 이유도 같습니다. 이 자료는 일반 기둥의 감쇠계수를 제공하는 문서가 아닙니다. 오히려 모든 구조에 하나의 임의 결함계수를 적용할 수 없고, 구조형식과 제조과정에 맞는 결함 정보가 필요하다는 좋은 반례입니다. 얇은 셸의 높은 결함 민감도를 보통의 압연 기둥에 그대로 이식해서도 안 됩니다.

7. 작은 편심은 임계점 전에 이미 2차효과를 만든다

하중이 중심에서 e만큼 벗어나면 처음부터 1차모멘트 Pe가 존재합니다. 기둥이 휘면서 생기는 추가편심 v까지 포함하면 모멘트는 대략 P(e+v)가 됩니다. 이상적인 핀–핀 탄성기둥의 중앙 최대모멘트에는 다음과 같은 secant 증폭이 나타납니다.

Mmax = Pe · sec[(π/2)√(P/Pcr)]

이 식은 실제 내력식이 아니라, 완전탄성 기둥에서 임계하중에 접근할수록 작은 편심의 효과가 얼마나 비선형적으로 커지는지 보여주는 사고도구입니다.

하중비 P/Pcr이상 secant 증폭읽어야 할 의미

0.25약 1.41작은 편심이 이미 모멘트를 키움

0.50약 2.251차모멘트의 두 배를 넘는 2차효과

0.75약 4.79강성저하가 응답을 지배하기 시작

0.90약 12.38완전계의 점근선 근처로 실제 결함·항복을 무시할 수 없음

점근적으로 무한대가 된다는 수학적 결과를 실제 응력이 무한대가 된다고 해석하면 안 됩니다. 실제 기둥은 그 전에 재료항복, 큰 변형, 접촉변화와 결함 때문에 다른 경로를 택합니다. 이 표의 역할은 내력을 정하는 것이 아니라 “Euler 하중 아래이므로 2차효과가 작다”는 추론이 틀릴 수 있음을 보여주는 것입니다.

8. Euler 기준선이 무너지는 세 영역

첫째, 짧은 부재에서는 재료의 압축항복이나 국부 손상이 전체 휨좌굴보다 먼저 일어날 수 있습니다. 둘째, 중간 세장비에서는 단면 일부가 항복하면서 유효 굽힘강성이 하중경로에 따라 변합니다. 셋째, 얇은 판으로 이루어진 단면에서는 전체 중심선이 휘기 전에 국부좌굴이 시작될 수 있습니다.

Shanley의 Inelastic Column Theory는 비탄성 좌굴에서 단순히 탄성계수 E를 임의의 작은 값으로 교체하는 방식이 왜 충분하지 않은지 보여준 고전적 연구입니다. 하중 증가 중 단면의 일부가 소성화되고 일부가 탄성제하될 수 있으므로, 접선강성과 평형경로를 일관되게 다뤄야 합니다. 중간 세장비의 실제 거동은 완전탄성 Euler 곡선과 단순 압축항복 중 작은 값을 고르는 것보다 복잡합니다.

따라서 손계산 기준선을 사용할 때는 예측 좌굴응력, 재료의 비선형 시작영역, 단면의 판세장성과 가능한 접합부 파괴를 함께 비교해야 합니다. 이 비교는 특정 안전계수를 정하려는 것이 아니라 어떤 물리모델이 다음 단계에 필요한지 선택하기 위한 것입니다.

9. 좌굴모드를 먼저 분류해야 요소를 고를 수 있다

전체 휨좌굴은 부재 중심선이 한 축을 중심으로 옆으로 이동하는 모드입니다. 비틀림좌굴은 단면 회전이 지배하고, 휨–비틀림좌굴은 이동과 회전이 결합됩니다. 국부좌굴은 플랜지나 웨브 같은 판 요소가 훨씬 짧은 파장으로 주름집니다. 보의 횡비틀림좌굴은 압축플랜지의 횡이동과 단면 비틀림이 결합된 별도의 문제입니다.

1차원 보 요소는 전체 휨좌굴과 골조 상호작용을 효율적으로 표현하지만, 판의 국부 주름을 직접 만들 수 없습니다. 셸 요소는 국부와 전체모드의 상호작용을 표현할 수 있지만 예상 반파장 안에 충분한 요소가 필요합니다. 솔리드 요소는 접촉과 두께방향 응력을 자세히 볼 수 있지만 얇은 부재의 좌굴에서는 계산비용이 급증합니다. “가장 정교해 보이는 요소”가 아니라 예상모드의 자유도와 파장을 담는 가장 단순한 모델을 선택해야 합니다.

10. 임계하중이 같아도 파괴의 급격함은 다르다

임계점 이후에도 변위와 함께 하중이 증가하면 안정 후좌굴입니다. 최대하중 뒤 하중이 감소하면 한계점과 연화가 나타납니다. 경로가 하중–변위 평면에서 급격히 돌아가면 스냅스루 또는 스냅백 거동이 생길 수 있습니다. 이 세 경로는 초기 임계하중만으로 구분할 수 없습니다.

Riks의 1979년 원 논문은 하중이나 변위 하나만을 진행변수로 삼지 않고 평형경로의 호길이를 따라 한계점과 분기점 부근을 추적하는 접근을 제시했습니다. Abaqus의 불안정 붕괴·후좌굴 문서도 하중크기를 미지수로 포함하는 호길이 방법이 불안정 정적응답을 추적하는 이유를 설명합니다.

그러나 호길이 방법은 모든 불안정을 자동으로 해결하는 버튼이 아닙니다. 갑작스러운 접촉상실, 국부 동적좌굴, 속도의존 재료, 비비례하중에서는 정적 평형경로 자체가 충분한 설명이 아닐 수 있습니다. 해법은 먼저 예상되는 물리현상과 제어가능한 실험변수를 결정한 뒤 선택해야 합니다.

11. 고유치 해석의 숫자를 안전여유로 읽으면 안 되는 이유

선형 고유치 좌굴은 기준상태의 접선강성이 특이해지는 하중배수와 모드를 계산합니다. Abaqus 고유치 좌굴 이론 문서의 표기처럼 기준하중을 PB, 섭동하중 패턴을 Q라 하면 예상 좌굴상태는 PB + λiQ입니다. 따라서 화면에 표시된 λ만 복사하면 기준하중과 섭동하중의 크기를 잃어버립니다.

증거직접 답하는 질문직접 답하지 못하는 질문

Euler·손계산 기준선단순 경계와 탄성가정에서 크기·축·길이 민감도실제 결함을 가진 최대하중

선형 고유치이상모델의 임계하중배수와 가능한 모드모드 진폭, 영구변형, 후좌굴 경로

결함 포함 비선형해석선택한 결함·재료·구속 아래 평형경로와 한계점입력하지 않은 제조편차와 모든 실제 불확실성

실험시험체와 시험경계가 보인 실제 경로다른 형상·공정·경계로의 자동 일반화

고유치 1.82는 곧바로 “안전율 1.82”가 아닙니다. 이상화된 기하, 탄성재료, 선택한 경계와 하중패턴에서 접선강성이 특이해지는 배수일 뿐입니다. 가까운 고유값이 여러 개라면 실제 결함은 한 모드가 아니라 여러 모드의 조합을 자극할 수 있으므로 낮은 모드를 묶어서 살펴봐야 합니다.

12. 비선형 모델에 정답을 몰래 넣지 않는 법

완전 대칭모델은 아무 교란이 없으면 직선경로에 머무를 수 있습니다. 그래서 비선형 후좌굴해석에는 초기결함이 필요합니다. 가장 좋은 근거는 측정된 초기형상입니다. 측정값이 없다면 제작공정에서 예상되는 결함형상, 여러 고유모드의 조합, 하중편심과 지점 유연성을 민감도 변수로 다뤄야 합니다. Abaqus 초기결함 문서는 고유모드 중첩, 이전 정적해석 변위와 직접 좌표입력이라는 서로 다른 방법을 구분합니다.

첫 번째 고유모드를 하나의 임의 크기로 넣은 결과는 결론이 아니라 한 시나리오입니다. 결함의 부호를 뒤집고, 가까운 모드를 조합하고, 크기를 단계적으로 바꾸며 최대하중뿐 아니라 초기강성, 모드 전환, 소성영역, 반력평형과 종료원인을 비교해야 합니다. 메시를 세분할수록 최대하중이 계속 떨어진다면 미해상 국부모드, 손상모델의 길이척도 또는 접촉정의를 의심해야 합니다.

지점과 하중도 같은 수준으로 의심해야 합니다. 끝단 모든 노드를 강체로 묶으면 단면변형과 회전을 과도하게 막을 수 있고, 잘못된 자유도를 풀면 강체운동이 첫 고유모드로 나타날 수 있습니다. 압축하중을 실제 도입판에 분포시킬지, 기준점 결합으로 줄지에 따라 국부모드가 달라질 수 있습니다.

13. 숫자 하나 대신 증거 사슬을 만든다

좋은 좌굴보고서는 가장 화려한 변형그림보다 가정의 추적성이 선명합니다. 다음 순서를 권합니다.

손계산 기준선: 단위, 길이, 축, I와 경계조건을 명시하고 길이 두 배 민감도를 확인합니다.

고유치·모드: 여러 낮은 모드의 하중간격과 형상, 강체모드 여부를 기록합니다.

결함 근거: 측정, 공정, 형상허용범위 또는 모드조합 중 무엇을 사용했는지 남깁니다.

비선형 경로: 하중–축단축, 하중–횡변위와 에너지·반력평형을 함께 봅니다.

민감도 포락선: 메시, 결함, 재료, 끝단과 브레이싱 강성을 바꿉니다.

실험 상관: 최대하중만 맞추지 말고 초기강성, 좌굴위치, 모드와 후좌굴경로를 함께 비교합니다.

모든 결과가 한 점에 모일 필요는 없습니다. 불확실성에 따라 결과가 움직이는 방향과 폭이 재현 가능하면 더 정직한 판단이 됩니다. 최종 산출물은 하나의 ‘좌굴계수’가 아니라 가정과 민감도를 포함한 판단범위여야 합니다.

14. 현장에서 휜 기둥을 볼 때 남겨야 할 증거

파괴 후 사진 한 장만으로 원인을 단정하기 어렵습니다. 전체 휨좌굴이 먼저였는지, 국부좌굴이나 접합부 변형이 먼저였는지, 하중제거 뒤 형상이 얼마나 회복됐는지에 따라 해석이 달라집니다. 가능하면 하중축과 부재 중심선의 정렬, 양 끝 회전, 브레이싱의 이동과 반력흔적, 국부주름의 위치와 파장, 잔류변형 방향을 같은 좌표계로 기록해야 합니다.

시험에서는 횡변위 한 점만 측정하지 말고 여러 높이의 변위와 단면회전을 기록해야 모드형상을 재구성할 수 있습니다. 하중–축단축 곡선만으로는 전체 휨과 국부모드 전환을 구분하기 어렵습니다. 초기형상 측정은 하중 전 상태에서 수행해야 하며, 이후 측정한 변형에서 초기결함을 빼야 실제 증가분을 얻을 수 있습니다.

분석자는 “마지막으로 크게 보인 모드”와 “처음 내력을 떨어뜨린 모드”를 구분해야 합니다. 후좌굴 과정에서 2차 모드가 시각적으로 지배적이어도 최초 한계점은 다른 국부모드나 접합부 유연성에서 시작했을 수 있습니다.

15. 초보자에서 전문가까지 가져갈 하나의 정신모형

처음 배우는 단계: 압축력은 옆으로 휜 기둥에 PΔ 모멘트를 만들며, 길고 약한 축일수록 복원하기 어렵다고 기억합니다.

계산을 시작하는 단계: Euler 값을 절대내력이 아니라 단위와 경계조건을 검산하는 기준선으로 사용합니다. 길이, I, 축과 모드가 예상대로 변하는지 확인합니다.

해석을 수행하는 단계: 고유치와 비선형해석이 답하는 질문을 분리하고, 결함의 근거·메시·지점·하중도입 민감도를 남깁니다.

전문가 단계: 분기점, 한계점, 안정·불안정 후좌굴, 모드 상호작용과 결함 민감도를 평형경로와 에너지 관점에서 통합합니다. 한 계산값보다 관찰 가능한 증거와 모델 가정 사이의 추적성을 우선합니다.

최종 판단좌굴은 “언제 옆으로 휘는가”라는 한 문장이 아닙니다. 어떤 평형경로가 존재하고, 그중 어느 경로가 안정하며, 실제 결함이 어느 경로를 선택하고, 그 판단을 어떤 증거로 재현할 수 있는가를 묻는 문제입니다.

16. 원문 자료와 이 글에서 사용한 범위

MIT 16.001 — Buckling of Beams: 변형된 형상의 평형, Euler 기둥, 경계조건, 초기결함과 편심응답.

MIT 2.080 — Structural Mechanics Lecture 9: 위치에너지의 2차 변분과 안정·불안정 판정.

W.T. Koiter — Over de stabiliteit van het elastisch evenwicht: 완전계의 안정성과 초기결함 민감도를 연결하는 고전적 기반. 저장소 메타데이터와 서지정보를 사용했습니다.

NASA SP-8007-2020 Rev.2: 얇은 원통 셸의 결함 측정, 제작공정별 결함 서명과 결함 민감도. 일반 기둥의 계수로 전용하지 않았습니다.

F.R. Shanley — Inelastic Column Theory: 부분 항복과 비탄성 기둥이 단순 Euler 치환보다 복잡한 이유.

E. Riks — An incremental approach to the solution of snapping and buckling problems: 한계점과 분기점을 포함하는 평형경로의 수치 추적.

Abaqus — Eigenvalue buckling prediction: 기준하중·섭동하중·고유값과 모드의 정확한 의미.

Abaqus — Introducing a geometric imperfection: 모드 중첩, 정적변위와 직접입력에 의한 초기결함 정의.

Abaqus — Unstable collapse and postbuckling analysis: 호길이 기반 불안정 평형경로 해석의 목적과 제한.

이 글의 식과 도해는 일반 역학을 설명하기 위한 것입니다. 실제 구조물의 내력판정에는 해당 구조형식의 재료, 제작오차, 접합부, 하중이력, 환경과 요구 신뢰도를 별도로 정의해야 합니다.