증분동적해석(IDA) — 원논문이 정의한 곡선, 붕괴 규칙, 그리고 알고리즘

Vamvatsikos와 Cornell의 2002년 IDA 논문을 읽어 IM·DM의 정의, 연화·경화·위빙과 구조적 부활, DM 규칙과 IM 규칙 및 FEMA 20% 접선 기울기, 16·50·84% 분위수 요약, 그리고 기록당 해석 횟수를 몇 회로 줄이는 hunt & fill 알고리즘을 정리한다.

증분동적해석(IDA)은 지진기록 하나를 여러 강도로 키워 가며 비선형 시간이력해석을 반복해, 강도와 응답의 관계를 하나의 곡선으로 만드는 방법이다. 쓰는 목적은 붕괴까지의 전 구간을 보는 것이다 — 한 번의 시간이력해석이 한 점을 주는 데 비해, IDA는 탄성에서 항복을 지나 동적 불안정까지 이어지는 연속적인 그림을 준다.

지금 성능기반 내진설계에서 쓰는 용어와 붕괴 판정 규칙 대부분이 이 논문에서 정리됐다. 이 글은 Vamvatsikos & Cornell, "Incremental Dynamic Analysis", Earthquake Engineering and Structural Dynamics 31(3): 491–514 (2002)의 본문과 그림을 인용한다. 인용한 수치와 규칙은 모두 논문에서 읽은 것이다.

1. 한 점이 아니라 곡선을 얻는다

논문은 방법을 정적 푸시오버와의 유비로 소개한다 — "정적해석 한 번에서 증분 정적 푸시오버로 넘어가는 것과 유비적으로, 한 번의 시간이력해석을 지진 하중을 스케일링하는 증분 해석으로 확장하는 데 이른다".

세 가지 용어가 이 방법의 뼈대다.

용어내용

IM (강도척도)기록을 얼마나 키웠는지 나타내는 값. 논문의 주 선택은 1차 모드 스펙트럼 가속도 Sa(T1, 5%)이다

DM (손상척도)구조 응답. 주로 최대 층간변형각 θmax를 쓴다

IDA 곡선한 기록에 대해 IM을 올리며 얻은 DM의 궤적

대상 구조물도 구체적이다 — 논문은 주기 4초의 20층 강재 모멘트골조, 2.2초의 9층, 1.3초의 3층(파단되는 접합부와 P–Δ 효과 포함), 그리고 1.8초의 5층 강재 가새골조를 예제로 쓴다.

2. 곡선은 단조롭게 눕지 않는다

IDA 곡선의 성질을 다루는 절이 이 논문에서 가장 실무적이다. 곡선에는 "연화(softening)와 경화(hardening)가 번갈아 나타나는 구간"이 있다. 강도를 올렸는데 응답이 줄어드는 구간이 실제로 생긴다.

같은 5층 가새골조에 네 개의 기록을 넣으면 "점진적으로 붕괴로 향하는 응답에서, 빠르고 비단조적인 앞뒤 뒤틀림 거동까지" 서로 다른 모양이 나온다. 층별로 곡선을 그리면 이유가 보인다 — "2층의 극단적인 연화가 퓨즈처럼 작동해 다른 층의 부담을 덜어 주는 복잡한 위빙(weaving) 상호작용"이다.

극단적인 경우에는 이름이 따로 붙는다.

이러한 현상들이 구조적 부활(structural resurrection)을 설명한다. 이는 경화의 극단적인 경우로, 계가 붕괴로 향하다가 더 높은 지진 강도에서 다시 회복되는 것이다. — Vamvatsikos & Cornell, §3

그래서 강도를 조금 키웠더니 응답이 좋아졌다는 결과는 계산 오류가 아니라 알려진 현상일 수 있다. 논문은 부활 이후 구간을 어떻게 다룰지도 규정한다 — 보수적으로 가려면 "첫 번째(IM 기준) flatline 위쪽의 곡선 부분은 폐기"한다.

3. 붕괴를 어디로 정의할 것인가

IDA 곡선 위에서 한계상태를 잡는 규칙은 두 계열이다.

규칙정의와 성질

DM 기반"DM ≥ C_DM이면 한계상태를 넘었다"는 형식. 예로 FEMA의 θmax = 2%는 1형 접합부 강재 모멘트골조의 즉시입주(IO) 수준을 표시한다. 구현이 간단하고 붕괴 이외의 성능수준에 자연스럽다

IM 기반"IM ≥ C_IM이면 한계상태를 넘었다"는 형식. 붕괴 영역을 하나로 명확히 가른다는 장점이 있으나, 그 값을 곡선마다 개별적으로 정해야 한다

FEMA 20% 접선 기울기사실상 IM 기반 규칙 — "접선 기울기가 탄성 기울기의 20%가 되는 곡선상의 마지막 점"을 용량점으로 삼는다

20% 규칙의 논리를 논문이 직접 설명한다 — 곡선이 평평해지는 것은 동적 불안정의 지표이고, DM이 "무한대"로 가는 것은 수치적으로 얻을 수 없으므로 "θmax의 증가율이 초기(탄성) 증가율의 5배가 되는 지점으로 물러서서" 용량을 표시한다는 것이다.

여기에 두 가지 주의가 붙는다.

위빙 때문에 그런 점이 여러 개 생긴다 — "구조물이 붕괴로 향하는 듯하다가 조금 더 높은 IM에서 회복되는" 지점들이며, "원칙적으로 낮은 쪽 점들은 용량 후보에서 폐기해야 한다".

곡선이 매끄럽지 않다 — IDA 곡선은 기껏해야 구간별로 매끄럽고, 접선 기울기는 매끄러운 보간을 통해 근사해야 한다.

FEMA의 전역 붕괴 정의는 두 규칙의 OR 결합이다 — "20% 기울기 IM 규칙과 C_DM = 10% DM 규칙의 OR 접속이며, 둘 중 하나라도 성립하면 용량을 정의한다". 여기서 "20% 강성 조건은 임박한 붕괴를 탐지하고, 10% 상한은 모델을 신뢰할 수 없는 영역의 과도한 θmax 값을 막는다".

그리고 논문은 DM 기반 붕괴 규칙에 대해 날카로운 단서를 단다 — 붕괴 용량의 경우 그런 규칙은 "사실 모델 결함의 신호일 수 있다. 모델이 충분히 현실적이라면 유한한 DM 출력이 아니라 비수렴으로 붕괴를 보여야 한다".

4. 여러 기록을 하나로 요약한다

기록 하나의 IDA 곡선은 그 기록의 성질이다. 논문은 5층 가새골조에 30개 기록을 적용한 다중기록 IDA를 보이고, 요약 방법을 정한다 — 응답 스펙트럼에 평균·중앙값·16%·84% 스펙트럼이 있듯, 16% · 50% · 84% 분위수 IDA 곡선을 정의한다.

평균이 아니라 분위수를 쓰는 이유가 명확하다.

각 IM 수준에서 평균을 계산하는 대신 표본 중앙값과 16%, 84% 분위수를 계산한다. 이 값들은 각각 기록의 50%, 84%, 16%에서 붕괴가 일어날 때에야 무한대가 된다. — Vamvatsikos & Cornell, §5

평균은 기록 하나가 붕괴하는 순간 무한대가 되지만, 중앙값은 절반이 붕괴할 때까지 유한하게 남는다. 여기에 유용한 성질이 하나 더 붙는다 — 연속성과 단조성을 가정하면 "IM이 주어졌을 때 DM의 x% 분위수를 잇는 선은, DM이 주어졌을 때 IM의 (100−x)% 분위수를 잇는 선과 같다". 수요 관점과 용량 관점이 같은 곡선을 가리킨다.

5. 해석 횟수를 줄이는 알고리즘

IDA의 비용은 기록 수 × 강도 수준 수다. 논문은 단순한 등간격 증가(stepping)의 문제를 먼저 지적한다.

스텝 크기가 기록마다 너무 크거나 너무 작다.

flatline의 높이가 기록과 IM에 따라 달라서 해석 횟수의 분배가 불균형해진다 — 낮은 IM에서 평평해지는 곡선은 점을 적게 받고, 높은 IM까지 가는 곡선은 많이 받는다.

수요와 용량의 해상도가 결합된다 — 둘 다 스텝 크기와 같아진다.

해법은 hunt & fill 세 단계다.

헌팅 — 스텝을 기하급수 또는 이차급수로 키우며 붕괴가 나올 때까지 올라간다. "몇 번의 해석만으로 flatline을 괄호에 넣는다".

이분법 — 붕괴(예: 비수렴)가 검출되면 붕괴하는 최저 IM과 붕괴하지 않는 최고 IM 사이를 좁혀, "수요 해상도와 무관하게 용량에 대해 규정된 정확도"를 얻는다.

빈틈 채우기 — 커진 스텝이 남긴 간격을 절반씩 나눠 채워 수요 추정의 해상도를 올린다.

논문은 이 알고리즘이 "기록당 비선형 해석 횟수를 몇 회 수준으로 줄일 수 있다"고 적는다. 특히 관심사가 개별 IDA 곡선의 세부가 아니라 요약된 결과일 때 그렇다.

6. 정적 푸시오버와의 관계

논문은 중앙값 IDA 곡선과 정적 푸시오버(SPO) 곡선을 같은 축에 겹쳐 놓는다 — 20층 모멘트골조와 5층 가새골조 두 경우다. 결론은 "두 곡선이 비슷한 DM 범위를 보인다"는 것이다. 1차 모드 하중 패턴으로 얻은 정적 곡선이 중앙값 동적 응답의 형태를 상당 부분 예고한다.

이 관계가 이후 SPO2IDA 같은 근사 도구의 출발점이 됐다. 다만 이 논문 안에서의 결론은 형태의 유사성이지 등가성이 아니다.

7. 정리

IDA는 기록 하나를 여러 강도로 스케일링해 IM–DM 곡선을 만든다. 논문의 표준 선택은 IM = Sa(T1, 5%), DM = θmax다.

곡선에는 연화와 경화가 번갈아 나타나고, 극단적인 경우 구조적 부활 — 붕괴로 향하다 더 높은 강도에서 회복 — 이 나타난다.

층별 IDA에서는 한 층의 극단적 연화가 퓨즈로 작동해 다른 층을 덜어 주는 위빙 상호작용이 보인다.

붕괴 판정은 DM 기반과 IM 기반으로 나뉘며, FEMA의 20% 접선 기울기 규칙은 IM 기반이다 — θmax 증가율이 초기의 5배가 되는 마지막 점.

FEMA 전역 붕괴는 20% 기울기 규칙과 θmax = 10% 규칙의 OR이고, 10% 상한은 모델을 신뢰할 수 없는 영역을 막는 장치다.

DM 기반 붕괴 규칙은 모델 결함의 신호일 수 있다 — 현실적인 모델이라면 비수렴으로 붕괴를 보여야 한다.

다중기록 IDA는 16% · 50% · 84% 분위수로 요약하며, 이 값들은 각각 기록의 84% · 50% · 16%가 붕괴할 때 무한대가 된다.

IM 주어진 DM의 x% 분위수 곡선은 DM 주어진 IM의 (100−x)% 분위수 곡선과 같다(연속·단조 가정).

hunt & fill 알고리즘은 헌팅(기하급수 증가) → 이분법(용량 정확도) → 빈틈 채우기(수요 해상도)로 기록당 해석 횟수를 몇 회로 줄인다.

이 논문이 20년 넘게 인용되는 이유는 새로운 해석 기법을 제안해서가 아니라 이미 여러 형태로 흩어져 있던 것에 이름과 규칙을 붙였기 때문이다. 강도를 올리며 반복한다는 발상 자체는 1977년 Bertero까지 거슬러 올라간다고 논문이 밝힌다. 실무에서 이 논문을 다시 볼 자리는 붕괴 판정이다 — 해석이 발산한 지점을 곧바로 용량으로 쓸지, 20% 기울기까지 물러설지, 그리고 여러 후보점 중 어느 것을 택할지가 결과를 바꾼다.