IDEA Member의 GMNIA 결과와 해석적 해의 비교
1. 목적
본 문서의 목적은 IDEA Member 애플리케이션의 GMNIA(기하학적 및 재료적 비선형 해석, 불완전성 포함) 모듈을 검증하는 것입니다. IDEA Member에서 산출된 저항력은 압축 기둥에 대한 EN 1993-1-1 [1] 해석적 해와 비교됩니다.
2. 모델 설명
GMNIA 모듈 검증을 위해 총 24개의 개별 케이스를 분석하였습니다. 모든 케이스는 동일한 단면 HEB 200과 동일한 강재 등급 S 355를 사용합니다. 네 가지 경계 조건(FF; PP; FP; FF)에 대해 기둥의 상대 세장비(0.5; 1.0; 1.5)를 변화시키며 검토하였습니다. 두 주축 방향에 대한 좌굴 저항력이 검증됩니다.
그림 1: 검증에 사용된 다양한 경계 조건
모든 케이스는 다음과 같은 방식으로 명명됩니다: "FR_0.5_Y", 여기서 "FR"은 경계 조건, "0.5"는 상대 세장비, "Y"는 좌굴 축을 나타냅니다.
3. 초기 불완전성
압축 기둥의 초기 불완전성을 산정하기 위해 세 가지 방법을 사용하였습니다. 이는 각각 A, B, C로 지정됩니다.
방법 A – EN 1993-1-1:2005, 표 5.1에 따름
표 1: 부재의 초기 활 불완전성 설계값 e0/L
방법 B – prEN 1993-1-1:2020, 2차 초안 [2], 조항 5.3.3.1에 따름
\[ \frac{e_0}{L}=\frac{α}{ε} \beta \]
여기서:
- e0 – 초기 불완전성
- α – 1993-1-1, 표 6.1 [1]에 따른 해당 좌굴 곡선에 의존하는 불완전성 계수
\[\varepsilon = \sqrt{\frac{235}{f_y}}\]
- fy – 기둥 항복 강도 [MPa]
- β – 표 2에 따른 기준 상대 활 불완전성
- L – 부재 길이
표 2: 기준 상대 활 불완전성
방법 C – EN 1993-1-1:2005 [1], 조항 5.3.2 (11)에 따른 EUGLI(등가 고유 전체 및 국부 초기 불완전성) 방법.
\[ e_0=\alpha (\bar \lambda - 0.2) \frac{M_{Rk}}{N_{Rk}} \]
여기서
- e0 – 초기 불완전성
- α – 1993-1-1, 표 6.1 [1]에 따른 해당 좌굴 곡선에 의존하는 불완전성 계수
- \( \bar \lambda \)– 부재 상대 세장비
- NRk – 단면의 특성 모멘트 저항력
- MRk – 단면의 축력에 대한 특성 저항력
이 초기 불완전성은 탄성 좌굴 모드 형상의 초기 불완전성에 대한 부재의 휨 모멘트 응답에 따라 조정됩니다.
표 3: 산출된 초기 불완전성 값 – y-y 축
표 4: 산출된 초기 불완전성 값 – z-z 축
4. 해석적 해
기둥의 좌굴 저항력 산정을 위해 EN 1993-1-1 [1], 조항 6.3에 따른 다음 방법을 사용합니다:
\[ N_{cr} = \frac{\pi ^2 E I}{L_{cr}^2} \]
\[ \bar \lambda = \sqrt{\frac{A f_y}{N_{cr}}} \]
\[ \phi = 0.5 \left [1 + \alpha \left (\bar \lambda - 0.2 \right ) + \bar \lambda ^2 \right] \]
\[ \chi = \frac{1}{\phi + \sqrt{\phi^2 - \bar \lambda ^2}} \]
\[ N_{b,Rd} = \frac{\chi A f_y}{\gamma_{M1}} \]
5. 결과
IDEA Member에서 산출된 극한 저항력(초기 불완전성 A, B, C에 대한)은 압연 단면(EN)의 해석적 값 및 웨브-플랜지 반경을 제외한 표현(Ew)과 비교됩니다.
5.1 강축 좌굴
강축 좌굴에 대한 결과는 아래 표에 요약되어 있습니다.
표 5: 산출된 저항력 값 – y-y 축
차트 1: 산출된 저항력 값 – y-y 축
차트 2: 산출된 저항력 비교 – y-y 축
GMNIA 결과는 유로코드 해에 비해 보수적입니다. 이는 부분적으로 IDEA Member의 단면 모델링 방식에 기인하며, 위 차트의 파란색 기둥 값에서 확인할 수 있듯이 이 영향은 2% 미만입니다.
초기 불완전성의 선택은 최종 저항력에 큰 영향을 미칩니다. 방법 C는 유로코드 해석적 해와 비교하여 약간 보수적인 결과(< 4%)를 나타내는 반면, 방법 A와 B는 10~16% 낮은 저항력을 산출합니다.
그림 2: 저항 한계에서의 기둥 PP_1.0_Y 및 플랜지 소성 변형률
5.2. 약축 좌굴
약축 좌굴에 대한 결과는 아래 표에 요약되어 있습니다.
표 6: 산출된 저항력 값 – z-z 축
차트 3: 산출된 저항력 값 – z-z 축
차트 4: 산출된 저항력 비교 – z-z 축
마찬가지로, 결과는 유로코드 해에 비해 보수적입니다. 위 차트의 파란색 기둥 값에서 확인할 수 있듯이 단면 모델링의 영향은 2% 미만입니다.
방법 C에 따라 선택된 초기 불완전성은 약간 보수적인 결과(< 6%)를 나타내는 반면, 방법 A와 B는 유로코드 해석적 해와 비교하여 10~26% 낮은 저항력을 산출합니다.
6. 참고 문헌
[1] EN 1993-1-1: Eurocode 3: 강구조 설계 – 제1-1부: 일반 규칙 및 건물 규칙, CEN, 2005.
[2] prEN 1993-1-1: Eurocode 3: 강구조 설계 – 제1-1부: 일반 규칙 및 건물 규칙, 2차 초안, CEN, 2017.