건축 및 토목 구조물의 안전성 확보는 단순히 재료의 강도에 의존하는 것이 아니라, 구조 전체의 응력 집중과 변형 거동을 정밀하게 파악하는 데 달려 있다. 이러한 목적을 달성하기 위해 건선 구조물의 비선형 변형률 시뮬레이션은 필수적인 해석 기법으로 자리 잡았다. 본 포스팅에서는 해당 논문 내용을 바탕으로, 비선형 변형률 시뮬레이션의 이론적 배경부터 수치 모델 구성, 다양한 구조 예제 실험까지 상세하게 소개한다. 특히 이 글의 메인 키워드인 “비선형 변형률 시뮬레이션”은 서론에서도 두 번 이상 반복해서 설명하며, SEO 최적화를 반영한 콘텐츠로 구성했다.
1. 비선형 변형률 시뮬레이션이 필요한 이유
구조물은 단순하지 않다
우리가 실제로 마주하는 구조물들은 단순한 기둥이나 슬래브의 조합이 아니다. 이음새, 구멍, 홈, 연결부와 같은 다양한 형상이 존재하고, 이들은 응력 집중의 주요 원인이 된다. 예를 들어, 고층 빌딩의 기둥에 케이블이나 덕트가 지나가기 위해 뚫리는 구멍 하나가 전체 구조물의 안정성을 저하시킬 수 있다. 따라서 실제 사용 조건을 반영한 비선형 모델링이 반드시 필요하다.
단순 선형 모델의 한계
기존의 구조해석은 대부분 선형 탄성 이론에 기반한다. 하지만 구조물은 실제 하중에 따라 소성 변형, 즉 영구 변형이 발생한다. 선형 모델로는 이를 예측할 수 없다. 이 한계를 극복하기 위해 비선형 변형률 시뮬레이션이 사용되며, 실제 구조물의 탄소성(elastoplastic) 거동을 반영해 보다 현실적인 해석이 가능하다.
2. 이론적 기반과 수치 모델 구성
수학적 기반 방정식
논문에서는 다음의 수학적 모델을 통해 3차원 구조물의 응력 상태를 계산한다.
- 평형 방정식 : 구조물 내부 응력이 외부 하중과 균형을 이루는 조건
- 기하학적 방정식 : 변형률과 변위 간의 관계
- 재료 방정식 : 탄소성 재료의 응력-변형률 관계
- 경계 조건 : 구조물의 외곽에서 적용되는 제약 조건
이 모든 식들은 선형 부분(탄성)과 비선형 부분(소성)을 함께 다루며, 특히 Ilyushin의 탄소성 함수는 항복 이후 거동을 예측하는 핵심 요소로 작용한다.
Ilyushin의 탄소성 모델
A.A. Ilyushin이 제시한 함수는 다음과 같다.
- 항복 이전 : 탄성 상태 (σ = Dε)
- 항복 이후 : 소성 변형 발생 (σ = Dε – Dω(ε))
이 함수는 재료가 항복점에 도달했을 때 발생하는 비선형 증가를 효과적으로 모델링한다. 실제 구조 해석에서 이 모델은 반복 수치 해석을 통해 응력 재분포 현상도 고려할 수 있다.
3. 다양한 구조 형상에 대한 수치 실험
[사례 1] 공동의 위치와 배치 방식이 구조에 미치는 영향
구조 내부에 구형 공동(cavity)을 단일로 배치하거나, 수직 혹은 수평 주기로 반복 배치했을 때 구조물의 응력 집중은 다음과 같이 달라진다.
| 구조 형태 | 탄성 상태의 최대 응력 σ/E | 탄소성 상태의 최대 응력 σ/E |
| 단일 공동 | 1.29×10⁻³ | 1.03×10⁻³ |
| 1차 주기 공동 | 1.13×10⁻³ | 1.29×10⁻³ |
| 2차 주기 공동 | 0.963×10⁻³ | 0.874×10⁻³ |
→ 2차 주기 배치는 구조 전체에 응력을 골고루 분산시켜 플라스틱 영역 확산을 억제함.
[사례 2] 공동 + 수직 균열 구조의 응력 완화 효과
두 개의 공동 사이에 수직 균열을 삽입했을 때, 응력 집중이 최대 14%까지 완화되는 현상이 관측되었다. 이는 설계자가 응력 분산 경로를 의도적으로 설계할 수 있다는 가능성을 시사한다.
- 균열 유무에 따라 응력 강도 지수(K)가 감소
- 이는 실제 철골구조물, 콘크리트 슬래브에 적용 시 비파괴 설계 전략으로 활용 가능
4. 반복 수치 해석과 응력 수렴 분석
수치 해석에서는 Ilyushin 모델 기반 반복법을 사용하여, 7회 반복에서 수렴하는 결과를 얻었다. 실제로는 3~4회 반복만으로도 실용적인 정확도에 도달할 수 있다.
| 반복 횟수 | 변형률 εi | 응력 강도 σi/𝐸 |
| 0회 | 0.11390 | 0.11390 |
| 5회 | 0.11174 | 0.09925 |
| 7회 | 0.11165 | 0.09922 |
→ 반복이 진행될수록 응력-변형 값은 안정화되며, 초기값과의 차이도 감소한다.
5. 홈(groove) 형상 변화와 구조 안정성
반지름에 따른 응력 변화
구조물 외벽에 홈을 추가하고 그 반경을 변경하면 응력 및 변형률도 달라진다.
| 반지름 r₃ | 변형률 εi | 응력 σi/𝐸 |
| 0.0125 | 0.2636 | 0.1744 |
| 0.2000 | 0.11165 | 0.0992 |
| 0.3500 | 0.0967 | 0.0915 |
→ 홈이 얕으면 응력 집중이 높고, 깊을수록 분산되는 패턴을 보인다.
→ 건축 실무에서는 철근 콘크리트 구조의 홈 깊이를 기준으로 설계에 반영 가능.
단면(Z축)별 응력 집중 분석
Z값에 따른 구조물 내부 응력 분포를 보면 다음과 같은 특징이 드러난다.
- Z = 0 : 중심부 응력 집중
- Z = 0.2 : 로젠지 모양의 응력 분산
- Z = 0.5 : 외벽 표면에서 응력 거의 없음
→ 이를 통해 공동이나 홈 위치 조절로 응력 흐름을 제어할 수 있음
6. 실무에서의 활용과 시사점
이 연구에서 제시한 비선형 변형률 시뮬레이션은 단순한 이론이 아닌, 실제 구조 설계 및 시공 단계에서도 다양한 방식으로 활용될 수 있다.
- 콘크리트 슬래브 설계 시, 공동 배치를 통해 하중 집중을 분산
- 프리캐스트 구조물 제작 시, 응력 흐름을 고려해 내부 공동 설계
- 리노베이션 구조물 분석 시, 응력 재분포를 고려한 보강 설계 가능
결국, 이런 수치 모델은 구조 해석가와 설계자에게 리스크를 줄이고 비용을 절감하는 도구가 되어준다.