성능기반 내진설계와 증분동적해석(IDA): 붕괴 프라질리티와 붕괴여유비 CMR·ACMR (KDS 41 17 00 · FEMA P695)

요약 · 성능기반 내진설계(PBSD)의 최종 한계상태인 붕괴는 지진강도를 점증시키며 비선형 시간이력해석을 반복하는 증분동적해석(IDA)으로 포착한다. 이 글은 IDA의 IM–DM 정식화, 붕괴점(flatline·동적 불안정) 정의, 붕괴 프라질리티의 대수정규 적합(중앙값 Ŝ_CT·분산 β_RTR), 그리고 FEMA P695의 붕괴여유비 CMR=Ŝ_CT/S_MT·스펙트럼형상계수 SSF·보정여유비 ACMR을 유도한다. 22개 지반운동 예제로 Ŝ_CT=1.50 g, CMR=1.67, ACMR=2.18 ≥ ACMR_10%=1.96 판정과 MCE 붕괴확률 6.9 %까지 프로그램 재검산으로 제시한다(KDS 41 17 00 · ASCE 41-17 · FEMA P695).

STRUCTURAL ENGINEERING · 내진공학 / 성능기반설계

성능기반 내진설계와 증분동적해석(IDA)
— 붕괴 프라질리티와 붕괴여유비 CMR·ACMR (KDS 41 17 00 · FEMA P695)

성능기반설계(PBSD)증분동적해석(IDA)붕괴 프라질리티CMR·ACMR스펙트럼형상계수(SSF)FEMA P695KDS 41 17 00

현대 내진설계의 무게중심은 부재 단위의 힘·응력 검토에서 시스템 단위의 붕괴 안전성 정량화로 이동하고 있다. 성능기반 내진설계(PBSD)의 최종 한계상태는 붕괴방지(Collapse Prevention)이며, 이를 신뢰성 언어로 표현한 것이 붕괴 프라질리티(collapse fragility)다. 붕괴는 본질적으로 재료·기하 비선형이 지배하는 동적 현상이므로, 지반운동 강도를 점증시키며 비선형 시간이력해석(NLTHA)을 반복하는 증분동적해석(IDA, Incremental Dynamic Analysis)이 이를 포착하는 표준 도구가 된다.

1. 기호 정의 (Notation)

표 1. 주요 기호 및 단위

기호	정의	단위
IM = Sa(T1)	지진강도지표 — 1차모드 5% 감쇠 스펙트럼가속도	g
DM = θmax	손상지표 — 최대 층간변위각(interstory drift)	rad
ŜCT	붕괴 중앙값 강도(median collapse intensity)	g
SMT	T1에서의 MCE 스펙트럼가속도	g
βRTR	지반운동 간 분산(record-to-record dispersion)	–
βTOT	총 붕괴 불확실도(RTR+DR+TD+MDL 합성)	–
CMR	붕괴여유비 = ŜCT / SMT	–
SSF	스펙트럼형상계수(spectral shape factor, ≥1)	–
ACMR	보정 붕괴여유비 = SSF · CMR	–
μT	주기기반 연성도(period-based ductility)	–
Φ( )	표준정규 누적분포함수, Φ⁻¹은 그 역함수	–

2. 왜 IDA인가 — PBSD의 붕괴 한계상태

탄성 응답스펙트럼해석은 부재가 항복하기 전까지의 거동만 기술한다. 그러나 붕괴는 소성힌지의 연쇄 형성 → 강성·강도 저하 → P-Δ 증폭 → 동적 불안정(측방붕괴, sidesway collapse)이라는 비가역적 경로를 따르므로, 이를 판단하려면 강도가 점증하는 여러 지반운동에 대한 비선형 동적 응답을 직접 추적해야 한다. IDA는 하나의 지반운동을 스칼라 배수로 스케일링하여 IM을 0에서부터 붕괴까지 점증시키고, 각 강도수준에서 DM(예: θ_max)을 기록해 IM–DM 곡선(단일기록 IDA 곡선)을 얻는다. 여러 지반운동에 대해 이를 반복하면 곡선의 다발(bundle)이 형성되고, 그 통계가 시스템의 붕괴 능력을 정의한다.

붕괴점은 두 방식으로 정의한다. (i) 동적 불안정: 해석이 발산(수치적 비수렴)하여 DM이 무한대로 향하는 IM. (ii) flatline(수평화) 규준: IDA 곡선의 국부 기울기가 초기 탄성기울기의 20%까지 저하되는 점. 두 규준은 대개 근접한 IM을 준다. 본 글은 이렇게 얻은 기록별 붕괴강도 S_a,col,i의 집합을 출발점으로 삼는다.

3. 증분동적해석의 정식화

IDA는 지진강도지표 IM과 손상지표 DM 사이의 사상을 구성한다. IM은 구조물의 응답과 상관이 높고 분산이 작은 척도가 바람직하며, 통상 1차모드 스펙트럼가속도 S_a(T₁, 5%)를 사용한다. DM은 붕괴와 직결되는 최대 층간변위각 θ_max를 쓴다.

IM ≡ S_a(T₁, 5%)  [g]  ,  DM ≡ θ_max  [rad](3.1)

기록 i를 배수 λ로 스케일링하면 IM_i = λ·S_a,i(T₁)이며, 각 λ에서 NLTHA를 수행해 θ_max를 얻는다. 붕괴판정 규준을 만족하는 최소 IM이 그 기록의 붕괴강도 S_a,col,i이다. 초기에는 IM∝θ_max(탄성, equal-displacement 근사)이지만, 소성화가 진행되면 곡선이 우향으로 완만해지다가 flatline에 도달한다(그림 1).

4. 붕괴 프라질리티 — 대수정규 적합

N개 기록의 붕괴강도 {S_a,col,i}는 대수정규분포로 잘 근사된다. 붕괴 프라질리티는 "IM = S_a일 때 붕괴할 조건부 확률"이며 다음과 같다.

P[ C | S_a ] = Φ( ln( S_a / Ŝ_CT )β_RTR )(4.1)

여기서 중앙값 Ŝ_CT와 분산 β_RTR은 최우추정(대수공간의 평균·표준편차)으로 구한다.

Ŝ_CT = exp( 1N Σ_i ln S_a,col,i )  ,  β_RTR = √ 1N Σ_i ( ln S_a,col,i − ln Ŝ_CT )² (4.2)

Ŝ_CT는 지반운동의 절반이 붕괴하는 강도(50 % 분위)에 해당하고, β_RTR은 지반운동 특성의 무작위성만을 반영한 지반운동 간 분산이다. 설계·시험·모델링에서 오는 인식론적 불확실도는 5절에서 별도로 합성한다.

그림 1. 세 지반운동에 대한 IDA 곡선(모식도). 초기 탄성구간(IM∝θ_max)에서 소성화로 우향 완만화되다가, 기울기가 초기값의 20%로 저하되는 flatline에서 붕괴한다. 붕괴강도 S_a,col은 기록마다 달라(A≈1.18, B≈1.50, C≈1.95 g) 분산 β_RTR를 형성한다. 붉은 파선은 MCE 요구강도 S_MT=0.90 g. Figure 1 — 필자 직접 작도(inline SVG), 저작권 없음.

5. 붕괴여유비 CMR·SSF·ACMR

붕괴여유비(Collapse Margin Ratio)는 붕괴 중앙값 강도가 MCE 요구강도를 얼마나 초과하는가를 나타낸다.

CMR = Ŝ_CTS_MT(5.1)

그런데 IDA에 사용한 지반운동을 설계스펙트럼에 맞춰 스케일링하면, 실제 강진 기록이 갖는 스펙트럼 형상(peak/valley, ε 효과)보다 넓은 대역에서 응답이 커져 붕괴가 보수적으로 평가된다. FEMA P695는 이를 스펙트럼형상계수 SSF(≥1)로 보정한 보정 붕괴여유비 ACMR을 정의한다.

ACMR = SSF · CMR(5.2)

합격 판정에는 지반운동 무작위성 β_RTR에 더해 설계요구도(β_DR)·실험자료(β_TD)·수치모델(β_MDL)의 인식론적 불확실도를 제곱합근으로 합성한 총 불확실도를 쓴다.

β_TOT = √ β_RTR² + β_DR² + β_TD² + β_MDL² (5.3)

붕괴확률이 β_TOT의 대수정규를 따른다고 보면, 목표 붕괴확률 x에 대응하는 허용 ACMR은 아래 역함수로 산정된다(예: x=10 %면 Φ⁻¹(0.90)=1.28).

ACMR_x% = exp( β_TOT · Φ⁻¹(1−x) )(5.4)

판정 규칙: 개별 아키타입은 ACMR ≥ ACMR_20%, 성능그룹 평균은 ACMR̅ ≥ ACMR_10%를 만족해야 한다.

6. 수치예제 — 22개 지반운동 IDA

T₁=1.0 s, 주기기반 연성도 μ_T=6인 특수모멘트골조 아키타입을 FEMA P695 원거리(Far-Field) 세트 규모의 22개 지반운동으로 IDA 해석하여 얻은 기록별 붕괴강도(정렬)는 다음과 같다.

S_a,col,i (g) = 0.65 · 0.80 · 0.90 · 0.99 · 1.06 · 1.13 · 1.20 · 1.26 · 1.33 · 1.40 · 1.46 · 1.54 · 1.61 · 1.69 · 1.78 · 1.88 · 1.99 · 2.12 · 2.28 · 2.49 · 2.80 · 3.48  (N=22)

Step 1 — 프라질리티 적합 (식 4.2). 대수공간 평균·표준편차로부터

Ŝ_CT = exp( 122 Σ ln S_a,col,i ) = 1.50 g  ,  β_RTR(자료) = 0.41(6.1)

자료 분산 0.41은 FEMA P695가 연성 시스템(μ_T ≥ 3)에 규정한 β_RTR = 0.40과 사실상 일치한다. 이후 합성에는 규정값 0.40을 사용한다.

Step 2 — CMR (식 5.1). KDS 41 17 00 설계응답스펙트럼(설계지진)에 위험도계수를 반영해 산정한 MCE 스펙트럼가속도가 T₁=1.0 s에서 S_MT = 0.90 g(= 설계값 0.60 g × 1.5)라 하면

CMR = Ŝ_CTS_MT = 1.500.90 = 1.67(6.2)

Step 3 — SSF·ACMR (식 5.2). FEMA P695 표 7-1b(SDC D_max, T=1.0 s, μ_T=6 보간)에서 SSF = 1.31이므로

ACMR = SSF · CMR = 1.31 × 1.67 = 2.18(6.3)

Step 4 — 불확실도·허용여유비 (식 5.3, 5.4). 설계요구도·시험자료·모델링 품질을 모두 "양호(B, β=0.20)"로 평가하면

β_TOT = √0.40² + 0.20² + 0.20² + 0.20² = √0.28 = 0.53 (표 반올림 0.525)(6.4)

ACMR_10% = exp(0.525 × 1.28) = 1.96  ,  ACMR_20% = exp(0.525 × 0.842) = 1.56(6.5)

표 2. 붕괴여유비 산정 요약 (T₁=1.0 s, μ_T=6, SDC D_max)

항목	기호	값	비고
붕괴 중앙값 강도	ŜCT	1.50 g	식 4.2, 22개 기록
MCE 요구강도	SMT	0.90 g	KDS 41 17 00 스펙트럼×1.5
붕괴여유비	CMR	1.67	ŜCT/SMT
스펙트럼형상계수	SSF	1.31	FEMA P695 표 7-1b
보정 붕괴여유비	ACMR	2.18	SSF·CMR

표 3. 불확실도 합성과 합격 판정

βRTR	βDR	βTD	βMDL	βTOT	ACMR20%	ACMR10%	판정
0.40	0.20	0.20	0.20	0.525	1.56	1.96	합격

판정. ACMR = 2.18 ≥ ACMR_20% = 1.56 (개별 아키타입) 이자 ≥ ACMR_10% = 1.96 (성능그룹 평균) 이므로 붕괴 안전성 요구를 만족한다. 스펙트럼형상 보정 중앙값 Ŝ_CT′ = SSF·Ŝ_CT = 1.96 g에 대한 MCE 붕괴확률은 식 4.1 형식으로
P[C | MCE] = Φ( ln(0.90 / 1.96) / 0.525 ) = Φ(−1.49) ≈ 6.9 % (β_TOT 기준) 이다. 스펙트럼형상 보정 전(β_RTR만) 원자료 프라질리티로는 P ≈ 10.5 %로, SSF 보정이 붕괴여유를 실질적으로 증대시킴을 보인다.

그림 2. 22개 지반운동에서 적합한 붕괴 프라질리티 곡선(Ŝ_CT=1.50 g, β=0.41). MCE 요구강도 0.90 g에서 원자료 붕괴확률은 약 10.5 %, 스펙트럼형상 보정으로 중앙값이 1.96 g로 우향 이동해 유효 붕괴확률은 6.9 %로 감소한다. Figure 2 — 필자 직접 작도(inline SVG), 저작권 없음.

7. 기준 적용 — KDS 41 17 00 · ASCE 41-17 · FEMA P695

KDS 41 17 00(건축물 내진설계기준)은 성능기반 내진설계를 허용하며, 설계지진을 재현주기 약 2400년 지진(MCE)의 2/3로 정의한다. 따라서 MCE 스펙트럼가속도 = 설계값 × 1.5이고, 최상위 성능목표는 이 MCE에서의 붕괴방지(CP)이다. 다만 KDS는 붕괴여유비의 정량적 판정절차 자체는 규정하지 않으므로, 아래 표의 성능수준 체계는 KDS에서 취하고, 붕괴여유 정량화(CMR·ACMR)는 FEMA P695와 기존 건물 평가기준 ASCE 41-17의 방법을 준용한다.

표 4. 성능수준과 붕괴평가 지표의 대응 (개념 정리)

성능수준	영문	지진수준(예)	상태 / 지표
기능수행	OP / IO	빈번지진	즉시 사용, 미세 손상
인명안전	LS	설계지진	보수 가능, 국부 항복
붕괴방지	CP	MCE	붕괴 직전 — 본 글의 ŜCT, ACMR로 정량

주의. IDA의 신뢰도는 (i) 붕괴모드를 재현하는 열화 이력모델(강도·강성 저하, 핀칭), (ii) P-Δ 기하비선형의 포함, (iii) 목표 위험도에 부합하는 지반운동 선정·스케일링에 좌우된다. 세 요소 중 하나라도 부실하면 Ŝ_CT가 과대평가되어 ACMR이 비보수적이 된다. 지반운동 선정·스케일링(조건부평균스펙트럼 CMS 등)은 별도 주제로 다룰 가치가 있다.

8. 결론

붕괴는 개별 부재의 강도식이 아니라 시스템의 동적 불안정으로 정의되는 사건이며, 그 안전여유는 확률로 표현해야 한다. IDA는 지진강도 IM을 점증시켜 각 지반운동의 붕괴강도 S_a,col를 얻고, 이를 대수정규 프라질리티(Ŝ_CT, β_RTR)로 압축한다. FEMA P695의 CMR = Ŝ_CT/S_MT는 스펙트럼형상계수 SSF로 보정되어 ACMR이 되고, 총 불확실도 β_TOT로 정한 허용치와 비교해 합·부를 가른다. 본 예제의 아키타입은 ACMR = 2.18 ≥ ACMR_10% = 1.96, MCE 붕괴확률 ≈ 6.9 %로 붕괴 안전성을 확보했다. 핵심은 SSF·β_TOT가 각각 스펙트럼형상 편향과 인식론적 불확실도를 정량 반영한다는 점이며, 이 틀은 응답수정계수 R의 정당화(justification) 근거로도 확장된다. 모든 수치는 프로그램으로 재검산했다.

ⓘ 본 글은 이론·기준 해설이며 실제 설계·성능평가 시 프로젝트 조건과 현행 KDS·FEMA·ASCE 기준 원문으로 재검증해야 합니다 · 참고: KDS 41 17 00, ASCE/SEI 41-17, FEMA P695(2009), Vamvatsikos & Cornell(2002).