기하강성행렬 기반 2차 해석: 유한요소 좌굴 고유치와 KDS·AISC 360-22·EN 1993-1-1 직접 2차해석 정식화의 비교

기하강성행렬 기반 2차 해석: 유한요소 좌굴 고유치와 KDS·AISC 360-22·EN 1993-1-1 직접 2차해석 정식화의 비교

Geometric Stiffness Matrix and Second-Order Analysis: FE Buckling Eigenvalue and the Formulation of Direct Second-Order Analysis in KDS, AISC 360-22, and EN 1993-1-1

초록. 매트릭스(행렬) 구조해석에서 기하강성행렬(geometric stiffness matrix, K_G)은 축력에 의한 강성 변화를 접선강성에 반영하여 좌굴과 2차효과(P-Δ·P-δ)를 일관되게 다루는 핵심 도구이다. 본 고는 퍼텐셜에너지 원리로부터 보-기둥 요소의 일관 기하강성행렬을 유도하고, 이를 이용한 선형 좌굴 고유치 문제 det(K_e+λK_G)=0의 구조를 분석한다. 이어 이 이론이 실무 설계기준의 2차해석 요구—KDS 14 31의 직접해석, AISC 360-22의 직접해석법(DAM), EN 1993-1-1의 2차해석·기하불완전성—로 어떻게 번역되는지 3개 코드를 병렬 비교한다. 단일 요소 캔틸레버 예제에서 K_G 고유치가 Euler 임계하중을 0.75% 오차로 재현함을 보이고, 강성감소·통보하중·증폭계수의 코드 간 철학 차이를 논한다.

주제어: 기하강성행렬(geometric stiffness matrix), 2차해석(second-order analysis), 직접해석법(direct analysis method, DAM), 좌굴 고유치(buckling eigenvalue), P-Δ 효과, 접선강성행렬(tangent stiffness matrix), 탄성임계하중배율(α_cr)

목차

1. 서론
2. 이론적 배경: 기하강성행렬의 유도와 좌굴 고유치
3. 설계기준 비교 (KR / US / EU)
4. 설계 예제: 단일 요소 좌굴과 층 증폭계수
5. 고찰
6. 결론
7. 기호 및 약어
8. 참고문헌

1. 서론

구조물의 안정성 설계는 지난 20년간 유효길이계수(effective length factor, K)에 의존하던 관행에서 벗어나, 해석 단계에서 2차효과와 기하불완전성을 직접 반영하는 방향으로 재편되었다. 이 전환의 수학적 토대가 바로 매트릭스 해석의 기하강성행렬이다. 부재에 작용하는 축력은 횡방향 강성을 변화시키며(압축은 감소, 인장은 증가), 이 효과를 요소 수준에서 정식화한 것이 K_G이다. 접선강성 K_t=K_e+K_G는 좌굴 고유치 해석, 반복 2차해석, 비선형 경로추적의 공통 출발점이 된다.

실무 설계자에게 이 주제가 중요한 이유는, 세 개의 주요 코드—KDS 14 31, AISC 360-22, EN 1993-1-1—가 모두 "직접 2차해석"을 표준 경로로 승격시켰으나, 그 정식화 철학이 서로 다르기 때문이다. AISC는 강성감소(0.8τ_b)와 통보하중(notional load)으로 불완전성과 비탄성을 해석에 주입하는 반면, EN은 명시적 기하불완전성 형상과 등가수평력을 사용하고, 강성감소는 별도 신뢰성 논리로 처리한다. 이 차이의 뿌리를 이해하려면 K_G의 유도 자체로 돌아가야 한다.

본 고는 부재·골조 수준의 탄성 2차효과와 좌굴 고유치에 한정한다. 재료비선형을 명시적으로 추적하는 소성힌지·화이버(fiber) 기반 비선형해석, 판좌굴·국부좌굴의 판이론적 처리, 동적 안정성(파라메트릭 공진)은 각각 별도 논고에서 다룬다. 즉 본 고의 K_G는 선형화된 기하강성(2차이론, second-order theory)에 국한하며, 대변형·후좌굴 거동(3차이론)은 범위 밖이다.

2. 이론적 배경: 기하강성행렬의 유도와 좌굴 고유치

2.1 퍼텐셜에너지로부터의 유도

축력 N(인장 양)을 받는 보-기둥 요소의 총 퍼텐셜에너지는 굽힘 변형에너지와 축력이 횡변위에 대해 하는 일(기하항)로 구성된다. 축력이 유발하는 2차 축방향 단축량은 횡변위 기울기의 제곱에 비례하므로, 기하항은 식 (1)로 표현된다.

$$\begin{array}{}\text{(1)}& \mathrm{\Pi }=\underset{U_b}{\underbrace{\frac{1}{2}{\int }_0^{L}EI(v^{″}{)}^{2}dx}}+\underset{U_g}{\underbrace{\frac{1}{2}{\int }_0^{L}N(v^{\prime }{)}^{2}dx}}\end{array}$$

여기서 v(x)는 횡변위, v′=dv/dx이다. 식 (1)의 두 번째 항의 물리적 기원은 축방향 변형률의 비선형(Green–Lagrange) 성분에 있다. 축방향 변형률을 유한변형 이론으로 전개하면 ε_xx=u′+½(v′)²이고, 이 중 ½(v′)² 항이 횡변위에 의한 부재 단축을 나타낸다. 축력 N이 이 단축에 대해 하는 일을 적분한 것이 식 (1)의 기하항 U_g이며, 따라서 K_G는 재료 물성이 아니라 현재 응력상태(축력)에 의존한다. 이 점이 탄성강성 K_e와의 근본적 차이로, 매 하중단계에서 축력이 갱신되면 K_G도 재계산되어야 하는 이유다. 압축(N<0)일 때 U_g는 음의 강성 기여를 주어 접선강성을 잠식한다.

변위장을 절점 자유도 d=[v₁, θ₁, v₂, θ₂]^T와 3차 Hermite 형상함수 N(x)로 보간하면 v=N d, v′=N′d이고, 두 에너지항은 각각 이차형식으로 정리된다. 정확히 이 지점에서 형상함수의 선택이 결과를 좌우한다. 굽힘의 정확해가 3차 다항식이므로 3차 Hermite 보간은 등단면·등축력 요소에 대해 굽힘 강성을 정확히 재현하지만, 기하항의 피적분함수 (N′)^T(N′)는 4차이므로 좌굴 형상에 대해 근사가 된다. 이 근사가 곧 2.4절에서 논하는 상측 수렴 오차의 근원이다.

$$\begin{array}{}\text{(2)}& U_b=\frac{1}{2}{\mathbf{d}}^T{\mathbf{k}}_e\mathbf{d},U_g=\frac{1}{2}{\mathbf{d}}^T{\mathbf{k}}_G\mathbf{d}\end{array}$$ $$\begin{array}{}\text{(3)}& {\mathbf{k}}_e={\int }_0^{L}EI({\mathbf{N}}^{″}{)}^T({\mathbf{N}}^{″})dx,{\mathbf{k}}_G=N{\int }_0^L({\mathbf{N}}^{\prime }{)}^T({\mathbf{N}}^{\prime })dx\end{array}$$

식 (3)의 적분을 수행하면 굽힘요소의 잘 알려진 탄성강성행렬과 함께, 일관 기하강성행렬(consistent geometric stiffness matrix)이 식 (4)로 얻어진다. "일관"이란 탄성강성과 동일한 3차 Hermite 형상함수로부터 유도되었음을 의미한다.

$$\begin{array}{}\text{(4)}& {\mathbf{k}}_G=\frac{N}{30L}\left[\begin{array}{cccc}36& 3L& -36& 3L\\ 3L& 4L^2& -3L& -L^2\\ -36& -3L& 36& -3L\\ 3L& -L^2& -3L& 4L^2\end{array}\right]\end{array}$$

그림 1은 이 요소의 자유도 규약과 축력에 의한 강성 변화를 도식화한다.

그림 1. 2차원 보-기둥 요소의 절점 자유도 [v₁, θ₁, v₂, θ₂]와 축력에 의한 접선강성 변화. 실선은 초기형상, 적색 곡선은 좌굴 변형형상을 나타낸다.

2.2 일관 강성 대 단순(P-Δ) 강성

회전 자유도를 포함하지 않고 절점 횡변위만으로 축력효과를 근사하면, 식 (4)는 이른바 단순(스트링) 기하강성 또는 P-Δ 기하강성으로 축약된다.

$$\begin{array}{}\text{(5)}& {\mathbf{k}}_G^{P\text{-}\mathrm{\Delta }}=\frac{N}{L}\left[\begin{array}{cccc}1& 0& -1& 0\\ 0& 0& 0& 0\\ -1& 0& 1& 0\\ 0& 0& 0& 0\end{array}\right]\end{array}$$

식 (5)는 부재 내 곡률에 의한 P-δ(member curvature)를 담지 못하고 절점 간 상대 변위에 의한 P-Δ(sidesway)만 포착한다. 상용 프로그램의 "P-Delta" 옵션 다수가 식 (5) 수준이며, 부재 좌굴이 지배적인 경우 반드시 부재를 다분할하거나 식 (4)의 일관 강성을 사용해야 한다. 그림 2는 두 효과의 기하학적 구분을 보인다.

그림 2. P-Δ(절점 상대 횡변위)와 P-δ(부재 내 곡률) 효과의 구분. 단순 기하강성 식 (5)는 좌측만, 일관 기하강성 식 (4) 또는 부재 다분할은 우측까지 포착한다.

2.3 접선강성과 좌굴 고유치 문제

요소 기하강성을 좌표변환·조립하여 전역 접선강성행렬을 구성한다. 하중을 기준하중 P_ref에 배율 λ를 곱한 것으로 두면, 기하강성은 축력에 선형이므로 K_G(λ)=λK_G(P_ref)이다.

$$\begin{array}{}\text{(6)}& {\mathbf{K}}_t(\lambda )={\mathbf{K}}_e+\lambda {\mathbf{K}}_G\end{array}$$

좌굴은 접선강성이 특이(singular)해져 인접평형이 존재하는 순간으로 정의되므로, 선형(고유치) 좌굴해석은 식 (7)의 일반화 고유치 문제로 귀결된다.

$$\begin{array}{}\text{(7)}& det({\mathbf{K}}_e+\lambda {\mathbf{K}}_G)=0⟺{\mathbf{K}}_e\mathit{\varphi }=-\lambda {\mathbf{K}}_G\mathit{\varphi }\end{array}$$

최소 양의 고유치 λ₁이 탄성임계하중배율(elastic critical load factor)이며, 대응 고유벡터 φ가 좌굴모드이다. 이 λ₁은 뒤에서 다룰 EN 1993-1-1의 α_cr 및 AISC의 층 전체 좌굴하중 P_e,story와 개념적으로 동일한 양이다. 다만 식 (7)은 선형화 좌굴로서 하중경로에 따른 K_G의 갱신을 무시하므로, 실제 분기하중을 상회 추정하는 경향이 있다. 이는 2.4절의 수렴 특성과 4장 예제에서 정량화한다.

식 (7)은 두 가지 실무적 함의를 갖는다. 첫째, 이는 대칭 양정치 K_e와 대칭(부정치일 수 있는) K_G의 일반화 고유치 문제이므로, 표준적으로 K_e를 Cholesky 분해한 뒤 표준 고유치 문제로 환원하여 부공간 반복법(subspace iteration)이나 Lanczos법으로 최소 몇 개의 λ만 효율적으로 추출한다. 대형 골조에서 전체 스펙트럼을 구할 필요는 없다. 둘째, 인장이 지배하는 부재(N>0)는 K_G에 양의 기여를 주어 실효 강성을 높이는데(그림 1의 stress stiffening), 이는 케이블·막구조의 형상강성과 회전기계 블레이드의 원심강화를 설명하는 동일 메커니즘이다. 즉 식 (4)의 K_G는 좌굴뿐 아니라 인장구조 해석의 공통 커널이다.

2.4 이론의 역사와 수렴 특성

기하강성행렬의 체계적 정식화는 Turner 외(1960)의 대변형 항 처리와 Martin, Gallagher 등의 1960년대 매트릭스 안정성 연구에서 정착되었고, Timoshenko & Gere(1961)의 고전 안정론이 폐형 해의 준거를 제공하였다. 교과서적 종합은 McGuire, Gallagher & Ziemian(2000)과 Bathe(2014), Cook 외(2002)에 정리되어 있다. 핵심 성질은 다음과 같다. 변위기반 요소는 강성을 과대평가하므로 좌굴하중을 상측에서 근사하고, 요소 분할을 늘리면 정해에 단조 수렴한다. 단일 요소로도 일관 기하강성은 오차 1% 내외의 실용 정확도를 보이는데, 이는 4장에서 캔틸레버에 대해 0.75% 오차로 확인된다.

여기서 실무자가 유의할 세부는 어떤 기하강성을 소프트웨어가 사용하는가이다. 일관 기하강성(식 4)은 굽힘 좌굴모드를 요소 내부에서도 근사하므로 요소당 정확도가 높은 반면, 스트링 강성(식 5)만 구현한 프로그램은 부재 좌굴을 절점 간 직선 변위로만 근사하여 한 부재를 하나의 요소로 두면 P-δ를 전혀 포착하지 못한다. 안정함수(stability functions)를 사용하는 정밀 정식화도 있는데, 이는 폐형 안정론 해를 요소 강성에 직접 삽입하여 축력에 대해 초월함수(삼각·쌍곡선) 형태의 강성계수를 쓰는 방식으로, 요소당 정확도는 최고이나 축력의 부호·크기에 따라 계수가 달라져 반복 재계산이 필요하다. 실무 상용해석에서는 계산 안정성과 자동화 편의 때문에 다항식 기반 일관 기하강성에 요소 다분할을 결합하는 절충이 표준으로 자리잡았다. 부재당 2~4분할이면 안정함수 수준의 정확도를 얻는다는 것이 통설이다.

3. 설계기준 비교 (KR / US / EU)

세 코드 모두 2차효과를 (a) 명시적 2차해석, (b) 근사 증폭, (c) 유효길이 기반의 경로로 허용하되, 어떤 조건에서 무엇을 강제하는가와 불완전성·강성감소의 주입 방식에서 갈라진다. 표 1이 그 핵심을 병렬한다.

표 1. 2차해석·안정성 설계 조항의 KR/US/EU 비교

항목	KDS 14 31 :2022 (한국)	AISC 360-22 (미국)	EN 1993-1-1:2005 (유럽)
표준 경로	직접해석(2차해석) 원칙, 유효길이법·1차해석 보완 허용	직접해석법(DAM), Chapter C 본문에 격상; 유효길이법은 Appendix 7로 이관	2차해석 원칙; 부재·골조 불완전성 반영(§5.2, §5.3)
2차효과 무시 한계	2차/1차 변위비 또는 αcr 기준으로 판정(AISC 계열 논리 준용)	Δ2nd/Δ1st ≤ 1.5 (또는 B2 ≤ 1.5)일 때 유효길이법 허용	αcr ≥ 10 (탄성해석) 또는 ≥ 15 (소성해석)이면 2차효과 무시 가능
기하불완전성	통보하중(중력하중의 약 0.2%) 또는 명시적 초기변형	통보하중 Ni = 0.002αYi (Chapter C)	등가수평력 φ = φ0αhαm, 부재 초기 굽힘 e0
강성감소	해석용 강성에 감소계수 적용(직접해석 채택 시)	0.8τb를 모든 강성에 적용 (τb는 축력비의 함수)	해석 강성 감소 대신 부재검토에 불완전성·재료계수 γM1로 반영
근사 증폭	증폭계수법 허용	B1(부재 P-δ), B2(층 P-Δ), Appendix 8	사이드스웨이 증폭 1/(1−1/αcr), §5.2.2
안전형식	하중저항계수설계(LRFD), 하중계수+저항계수 φ	LRFD(φ) 또는 ASD(Ω); DAM은 두 형식 모두 지원	부분계수형식: 작용 γF, 재료 γM (γM1 안정성)

차이의 기원. AISC의 DAM은 "해석에 불완전성과 비탄성을 넣어 부재는 K=1로 검토한다"는 단일 철학을 관철한다. 0.8 계수는 잔류응력·재료 불확실성을 포함한 계통적 강성 저하를, τ_b는 축력에 따른 비탄성 강성 저하를 대표하며, 통보하중 0.002는 시공 허용오차 L/500의 초기 기울기를 등가 수평력으로 환산한 값이다. 반면 EN 1993-1-1은 불완전성을 물리적 형상(전역 기울기 φ₀=1/200, 부재 초기 처짐 e₀)으로 명시하고, 비탄성은 해석 강성이 아니라 부재 좌굴검토의 감소계수 χ와 γ_M1로 흡수한다. 두 접근은 신뢰성 목표 β를 만족시키되, 보정 데이터(미국의 컬럼 시험 데이터베이스 대 유럽 ECCS 좌굴곡선)와 설계 문화가 달라 형식이 갈렸다. KDS 14 31은 AISC 계열의 직접해석 논리를 근간으로 채택하되 국내 하중·재료 조항과 정합시켜 운용된다. 세 코드의 공통 좌굴량은 모두 식 (7)의 λ₁—즉 α_cr 또는 P_e,story/ΣP—라는 점에서, 매트릭스 해석의 기하강성 고유치가 코드 간 공통 언어임을 알 수 있다.

보정 데이터의 계보를 좀 더 짚으면, 미국의 압축재 강도식은 SSRC(Structural Stability Research Council)가 정리한 다수 열간압연·용접 단면의 컬럼 시험과 잔류응력 측정에 근거한 단일 곡선(SSRC 곡선 2 계열)으로 정착된 반면, 유럽은 단면 형상·제작방법·좌굴축에 따라 a₀·a·b·c·d의 다섯 좌굴곡선(ECCS)으로 세분하였다. 이 세분화가 EN이 불완전성을 부재별 초기처짐 e₀과 감소계수 χ로 부재검토에 명시적으로 배분하는 형식을 택한 역사적 배경이다. 반대로 미국은 시스템 신뢰도 목표 β≈2.6~3.0(중력지배 대 조합)을 만족하도록 해석 강성감소와 통보하중을 계통 수준에서 보정하여, 부재검토를 K=1로 단순화하는 전략을 택했다. 어느 쪽도 물리적으로 "더 옳지" 않으며, 동일한 시험모집단을 서로 다른 방식으로 분해하여 유사한 목표 신뢰도에 도달한 두 갈래로 이해하는 것이 정확하다. 실무 검토자는 이 계보 차이를 알고 있어야, 두 형식의 안전여유를 이중 계상하는 과오를 피할 수 있다.

3.1 증폭계수의 정식화

AISC 360-22 Appendix 8의 층 증폭계수 B₂는 층 전체의 사이드스웨이 좌굴하중 P_e,story에 대한 층 중력하중 ΣP의 비로 정의된다.

$$\begin{array}{}\text{(8)}& B_2=\frac{1}{1-{\displaystyle \frac{\alpha \sum P_{\text{story}}}{P_{e,\text{story}}}}}\ge 1,\alpha =1.0(\text{LRFD})\end{array}$$

EN 1993-1-1 §5.2.2의 사이드스웨이 증폭은 탄성임계하중배율 α_cr=F_cr/F_Ed로 표현되어 식 (9)가 되며, α_cr=P_e,story/ΣP를 대입하면 식 (8)과 형태가 일치한다.

$$\begin{array}{}\text{(9)}& k_{\text{amp}}=\frac{1}{1-{\displaystyle \frac{1}{{\alpha }_{cr}}}},{\alpha }_{cr}\ge 10\Rightarrow \text{2차효과 무시}\end{array}$$

즉 두 코드의 증폭식은 대수적으로 동형이며, 차이는 "무시 한계"의 설정—AISC는 변위비 1.5(=α_cr≈3), EN은 α_cr≥10—에 있다. EN이 더 보수적인 한계를 두는 것은 소성해석 확장(≥15)과 국가부속서(National Annex) 변동을 함께 고려한 결과이다. AISC의 DAM 강성감소는 식 (10)으로 요약된다.

$$\begin{array}{}\text{(10)}& E{I}^{\ast }=0.8{\tau }_{b}EI,{\tau }_b=\{\begin{array}{ll}1.0,& \alpha P_r/P_y\le 0.5\\ 4{\displaystyle \frac{\alpha P_r}{P_y}}(1-{\displaystyle \frac{\alpha P_r}{P_y}}),& \alpha P_r/P_y>0.5\end{array}\end{array}$$

그림 3. 층 증폭계수 B₂=1/(1−ρ), ρ=ΣP/P_e,story=1/α_cr. AISC(변위비 1.5, ρ≈0.33)와 EN(α_cr≥10, ρ=0.1)의 2차효과 무시 한계를 함께 표시하였다. 두 코드의 증폭식은 동형이나 임계 판정선이 다르다.

3.2 불완전성·비탄성 주입 경로의 대조

세 코드의 실질적 차이는 "동일한 α_cr을 알고 난 뒤 무엇을 하는가"에 있다. 그림 4는 AISC 직접해석법(DAM)과 EN 2차해석 절차를 병렬 흐름으로 대비한다. AISC-DAM은 (i) 모든 부재 강성을 0.8τ_b로 감소시키고, (ii) 각 층에 통보하중 N_i=0.002αY_i를 가한 상태에서, (iii) P-Δ와 P-δ를 모두 포함한 2차해석을 수행한 뒤, (iv) 부재를 유효길이계수 K=1로 검토한다. 즉 안정성의 대부분을 해석이 흡수하고 부재검토는 단면 강도문제로 단순화된다.

EN 1993-1-1은 (i) 전역 불완전성 φ=φ₀α_hα_m을 등가수평력으로, 필요 시 부재 초기처짐 e₀을 형상으로 부여하고, (ii) α_cr≥10이면 1차해석, 그렇지 않으면 2차해석(또는 사이드스웨이 증폭)을 수행하며, (iii) 해석 강성은 감소시키지 않고 (iv) 부재는 좌굴감소계수 χ(좌굴곡선 a~d)와 γ_M1로 검토한다. 불완전성은 물리적 형상으로, 비탄성은 부재검토로 분리 배치된다. 두 경로가 유사한 신뢰도를 주도록 보정되어 있으나, 강성감소를 해석에 넣는 AISC와 넣지 않는 EN의 차이 때문에 동일 골조의 요구 단면이 갈릴 수 있다. 이 두 절차 중 어느 것이 "더 합리적"인가는 논쟁적이며, DAM은 소프트웨어 의존적·자동화 친화적인 반면 EN은 물리적 투명성이 높다는 것이 통상적 평가다.

그림 4. AISC 직접해석법(DAM)과 EN 1993-1-1 2차해석 절차의 병렬 비교. 진입 판정선인 α_cr은 두 코드가 공유하나, 불완전성과 비탄성을 해석(AISC)에 넣을지 부재검토(EN)에 넣을지가 갈린다.

4. 설계 예제: 단일 요소 좌굴과 층 증폭계수

4.1 캔틸레버 기둥의 좌굴하중 — 단일 요소 K_G 고유치

고정-자유(캔틸레버) 기둥을 하나의 보-기둥 요소로 모델링하여 식 (7)의 고유치로 임계하중을 구하고, 폐형 Euler 해와 비교한다. 절점 1은 고정(v₁=θ₁=0), 절점 2는 자유이므로 활성 자유도는 [v₂, θ₂]이다. 압축하중 P에 대해 N=−P를 식 (4)에 대입한다.

절점 2 블록만 추출하면(단위: EI, L):

$$\begin{array}{}\text{(11)}& {\mathbf{k}}_e=\frac{EI}{L^3}\left[\begin{array}{cc}12& -6L\\ -6L& 4L^2\end{array}\right],{\mathbf{k}}_G=-\frac{P}{30L}\left[\begin{array}{cc}36& -3L\\ -3L& 4L^2\end{array}\right]\end{array}$$

무차원 하중 p=PL²/EI를 도입하면 접선강성 블록은 EI/L³을 공통인수로 하여 식 (11a)가 된다.

$$\begin{array}{}\text{(11a)}& {\mathbf{K}}_t=\frac{EI}{L^3}\left[\begin{array}{cc}12-1.2p& (-6+0.1p)L\\ (-6+0.1p)L& (4-\frac{2}{15}p)L^2\end{array}\right]\end{array}$$

비자명 해 조건 det(K_t)=0에서 L²와 (EI/L³)²은 소거된다. 대각항 곱에서 비대각항 제곱을 빼면

$$\begin{array}{}\text{(12)}& (12-1.2p)(4-\frac{2}{15}p)-(-6+0.1p{)}^2=0.15p^2-5.2p+12=0\end{array}$$

가 되어 식 (12)의 특성방정식이 얻어진다. (검산: 48−6.4p+0.16p² − (36−1.2p+0.01p²) = 0.15p²−5.2p+12.)

근의 공식으로 p=(5.2±√(5.2²−4·0.15·12))/(2·0.15)=(5.2±√19.84)/0.3이며, 작은 근이 임계값이다.

최소 양의 근: p_cr = (5.2 − 4.4542)/0.30 = 2.486 →  P_cr,FE = 2.486 EI/L²

폐형 Euler 해(유효길이계수 K=2): P_cr,Euler = π²EI/(2L)² = 2.4674 EI/L²

상대오차 = (2.486 − 2.4674)/2.4674 = +0.75 % (상측 근사)

단 하나의 일관 기하강성 요소가 정해를 0.75% 오차로 재현하며, 부호(+)는 변위기반 요소의 강성 과대평가 성질(2.4절)과 정확히 일치한다. 요소를 2개로 분할하면 오차는 약 0.1% 이하로 떨어진다. 표 2는 요소 분할에 따른 수렴을 요약한다.

표 2. 캔틸레버 기둥 임계하중의 요소 분할별 수렴 (일관 기하강성)

요소 수 n	Pcr / (EI/L2)	Euler 대비 오차	비고
1	2.486	+0.75 %	식 (12) 해
2	2.4699	+0.10 %	단조 수렴(상측)
4	2.4676	+0.01 %	실무상 수렴
∞ (정해)	2.4674	0	π2/4

주: n=2 이상의 값은 동일한 일관 기하강성 정식(식 4)의 다요소 조립에 대한 수치 고유치해로, 단조 상측 수렴 경향을 예시한다.

4.2 층 2차효과 — 증폭계수 산정과 코드 간 대조

1개 층에 대해 총 중력하중 ΣP=4000 kN, 층 사이드스웨이 좌굴하중 P_e,story=20000 kN이라 하자. 그러면 ρ=ΣP/P_e,story=0.20, 등가 α_cr=1/ρ=5.0이다.

AISC 360-22 (식 8, LRFD α=1.0):

$$\begin{array}{}\text{(13)}& B_2=\frac{1}{1-1.0\cdot 0.20}=1.25\end{array}$$

EN 1993-1-1 (식 9):

$$\begin{array}{}\text{(14)}& k_{\text{amp}}=\frac{1}{1-{\displaystyle \frac{1}{5.0}}}=1.25\end{array}$$

ρ = 4000/20000 = 0.20,  α_cr = 5.0

AISC B₂ = 1.25  ≡  EN k_amp = 1.25 (동일)

판정: α_cr=5.0 < 10 이므로 EN 기준 2차해석 필수. AISC는 변위비(≈1.25 < 1.5)만으로는 유효길이법도 가능하나, DAM 채택 시 강성감소 0.8τ_b를 병행 적용.

두 코드가 동일한 증폭 1.25를 산출하는 것은 식 (8)·(9)의 동형성에서 필연이다. 그러나 실무 결과는 갈릴 수 있다. AISC-DAM은 해석 강성을 0.8배로 줄이므로 P_e,story가 감소해 ρ가 커지고 증폭이 상승하며, 부재는 K=1로 검토한다. EN은 해석 강성을 유지하되 부재 좌굴검토에서 χ·γ_M1로 안전여유를 확보한다. 따라서 동일 골조라도 요구 단면이 달라질 수 있으며, 이 차이는 4.1절에서 본 K_G 고유치(=P_e,story)를 어느 강성으로 계산하느냐에서 시작된다.

5. 고찰

선형 좌굴해석의 한계. 식 (7)의 고유치는 하중경로 중 K_G의 갱신을 무시한 분기 추정이다. 실제 골조는 불완전성으로 인해 분기 이전부터 휘어지는 극한점(limit-point) 거동을 보이므로, α_cr은 안정성 지표일 뿐 설계강도가 아니다. 세 코드가 α_cr(또는 변위비)을 스위치로만 쓰고 실제 검토는 증폭 2차해석+불완전성으로 넘기는 이유가 여기 있다.

단순 P-Δ의 함정. 식 (5)만 지원하는 해석은 부재 좌굴(P-δ)을 놓친다. 세장한 압축재가 지배하는 골조에서 부재를 분할하지 않고 단순 P-Δ만 켜면 임계하중을 위험측으로 과대평가할 수 있다. AISC의 B₁·B₂ 분리와 EN의 부재 불완전성 e₀는 모두 이 누락을 보완하는 장치이며, 직접해석 소프트웨어에서는 요소 다분할 또는 일관 기하강성 사용이 실무 표준이다.

강성감소의 신뢰성 논리. AISC의 0.8 계수는 특정 부재가 아니라 시스템의 계통적 강성 저하를 반영하도록 보정된 값으로, 이를 EN식 부재 감소계수와 이중 적용하면 과도한 보수성이 된다. 안전형식(φ 대 γ_M)을 혼용하지 않도록 signposting이 필요하며, 국내 KDS 적용 시에도 어느 강성으로 P_e,story를 산정했는지 계산서에 명시하는 것이 검토 상 필수이다.

모델링 상의 실무 쟁점. 층 좌굴하중 P_e,story의 산정 자체가 세심함을 요한다. AISC Appendix 8은 층 접근법(각 기둥의 강성 기여를 합산)을 제시하나, 강접골조와 가새골조가 혼재하거나 경사재·전이보가 있는 불규칙 골조에서는 층 개념이 모호해진다. 이 경우 식 (7)의 전역 고유치 좌굴해석이 층 공식보다 신뢰성 높은 α_cr을 준다. 다만 전역 좌굴모드가 국부(단일 부재) 모드로 나타나면 그 λ는 층 사이드스웨이 지표로 부적절하므로, 고유벡터 형상을 반드시 확인하여 사이드스웨이 모드를 선별해야 한다. 또한 중력기둥(leaning column)의 하중은 자체 횡강성이 없더라도 횡저항 시스템의 P_e,story를 잠식하므로 ΣP에 반드시 포함해야 하며, 이 누락은 실무에서 흔한 위험측 오류다.

최근 동향. 2세대 유로코드(EN 1993-1-1:2022)는 불완전성·2차해석 절차를 재정비하고 α_cr 기반 판정을 정교화하였으며, 국가부속서 의존 파라미터가 여전히 결과를 좌우한다. 연구 측면에서는 화이버 기반 비탄성 좌굴해석과 기하강성의 결합, 신뢰성기반(β 목표) 통보하중 재보정이 활발하다. 매트릭스 해석의 K_G는 이 모든 확장의 불변 핵심으로 남는다.

6. 결론

1. 기하강성행렬 K_G는 퍼텐셜에너지의 축력 일항(식 1)에서 유도되며, 일관 형태(식 4)는 탄성강성과 동일한 3차 Hermite 형상함수에 기반해 P-δ까지 포착한다.
2. 선형 좌굴은 일반화 고유치 문제 det(K_e+λK_G)=0으로 귀결되며, 최소 고유치 λ₁이 코드의 α_cr 및 P_e,story와 동일한 안정성 지표이다.
3. 단일 일관 요소로 캔틸레버 임계하중을 2.486 EI/L²로 산출하여 Euler 해를 +0.75% 오차(상측)로 재현하였고, 요소 분할로 단조 수렴함을 확인하였다.
4. AISC 360-22의 B₂와 EN 1993-1-1의 사이드스웨이 증폭은 대수적으로 동형(1/(1−1/α_cr))이나, 2차효과 무시 한계(변위비 1.5 대 α_cr≥10)가 달라 판정이 갈린다.
5. 불완전성·비탄성 주입 방식—AISC의 강성감소 0.8τ_b+통보하중 대 EN의 명시적 초기형상+χ·γ_M1—의 차이는 보정 데이터와 설계철학에서 비롯되며, 동일 골조에서도 요구 단면을 다르게 만들 수 있다.
6. 단순 P-Δ 강성(식 5)만 사용하는 해석은 부재 좌굴을 놓치므로, 세장 압축재 지배 골조에서는 요소 다분할 또는 일관 기하강성 사용이 필수이다.

핵심 요약 (Key Takeaways)

매트릭스 해석의 기하강성행렬은 좌굴 고유치, 2차 증폭, 직접해석법을 잇는 공통 언어이다. 세 코드의 증폭식은 α_cr로 통일되지만, 불완전성과 강성감소를 "해석에 넣을지(AISC/KDS-DAM) 부재검토에 넣을지(EN)"의 철학이 결과를 가른다. 계산서에는 어떤 강성으로 P_e,story=λ₁ΣP를 구했는지, 안전형식(φ 대 γ_M1)을 어떻게 분리했는지를 반드시 명시해야 한다.

기호 및 약어

표 3. 기호 및 단위

기호	의미	단위
E, I, EI	탄성계수, 단면2차모멘트, 휨강성	MPa, mm4, N·mm2
Ke, KG, Kt	탄성·기하·접선 강성행렬	N/mm 등
N	부재 축력(인장 양)	kN
P, Pcr	압축하중, 임계(좌굴)하중	kN
Pe,story	층 사이드스웨이 탄성좌굴하중	kN
p	무차원 하중 PL2/EI	—
λ, λ1	하중배율, 최소 좌굴 고유치	—
αcr	탄성임계하중배율 Fcr/FEd	—
B1, B2	부재(P-δ)·층(P-Δ) 증폭계수	—
τb	비탄성 강성감소계수(AISC)	—
φ0, e0	전역 초기기울기, 부재 초기처짐(EN)	rad, mm
φ, γM1	강도저항계수(LRFD), 안정성 부분계수(EN)	—
v, θ	절점 횡변위, 회전	mm, rad

참고문헌

1. 국가건설기준센터, KDS 14 31 :2022 강구조 설계기준(하중저항계수설계법), 국토교통부, 2022.
2. 국가건설기준센터, KDS 41 17 00 :2022 건축물 내진설계기준, 국토교통부, 2022.
3. American Institute of Steel Construction, ANSI/AISC 360-22, Specification for Structural Steel Buildings, AISC, Chicago, 2022 (Chapter C, Appendix 7, Appendix 8).
4. European Committee for Standardization, EN 1993-1-1:2005, Eurocode 3: Design of Steel Structures — Part 1-1: General Rules and Rules for Buildings, CEN, Brussels, 2005 (§5.2, §5.3); 2세대 개정 EN 1993-1-1:2022 참조.
5. American Society of Civil Engineers, ASCE/SEI 7-22, Minimum Design Loads and Associated Criteria for Buildings and Other Structures, ASCE, 2022.
6. McGuire, W., Gallagher, R. H., and Ziemian, R. D., Matrix Structural Analysis, 2nd ed., John Wiley & Sons, 2000.
7. Bathe, K.-J., Finite Element Procedures, 2nd ed., K.-J. Bathe (Prentice Hall), 2014.
8. Cook, R. D., Malkus, D. S., Plesha, M. E., and Witt, R. J., Concepts and Applications of Finite Element Analysis, 4th ed., John Wiley & Sons, 2002.
9. Timoshenko, S. P., and Gere, J. M., Theory of Elastic Stability, 2nd ed., McGraw-Hill, 1961.
10. Ziemian, R. D. (ed.), Guide to Stability Design Criteria for Metal Structures, 6th ed., John Wiley & Sons, 2010.
11. Nair, R. S., "A Model Specification for Stability Design by Direct Analysis," Engineering Journal, AISC, Vol. 46, No. 1, 2009, pp. 29–37.
12. American Institute of Steel Construction, Design Guide 28: Stability Design of Steel Buildings, AISC, 2013.