명시적 동해석

구조해석 고체역학 마지막 수정: 2026-07-06 04:51:33

명시적 동해석 Explicit Dynamic Analysis
분야	고체역학 × 구조해석
시간적분	중심차분법 (조건부 안정)
안정성 조건	CFL 조건
대표 소프트웨어	LS-DYNA, Abaqus/Explicit, PAM-CRASH
주 사용처	충돌, 성형, 낙하, 폭발 해석

1. 개요[편집]

암시적 해석이 “정답을 찾을 때까지 반복해서 조인다”면, 명시적 해석은 “일단 다음 순간으로 냅다 튀어나간다”.

명시적 동해석(Explicit Dynamic Analysis)은 구조물의 운동방정식을 시간에 대해 적분할 때, 다음 시각의 상태를 현재 시각의 정보만으로 직접(explicitly) 계산하는 구조해석 기법이다. 방정식을 반복적으로 풀어 평형을 맞추는 암시적(implicit) 방식과 달리, 명시적 방식은 각 시간 스텝에서 연립방정식을 풀지 않는다. 이 한 가지 차이가 명시적 동해석의 모든 성격을 결정한다. 짧은 시간에 매우 빠르고 격렬하게 벌어지는 현상 — 자동차 충돌, 금속 성형, 낙하 충격, 폭발 — 을 다루는 데 독보적이다.¹

지배 방정식은 반이산화된 운동방정식이다.

\mathbf{M}\ddot{\mathbf{u}} + \mathbf{C}\dot{\mathbf{u}} + \mathbf{K}\mathbf{u} = \mathbf{f}^{\text{ext}}

여기서 $\mathbf{M}$ 은 질량행렬, $\mathbf{C}$ 는 감쇠행렬, $\mathbf{K}$ 는 강성행렬이다. 명시적 기법의 마법은 이 방정식을 어떻게 시간적분하느냐에 있다.

2. 중심차분 시간적분[편집]

명시적 동해석의 심장은 중심차분법(central difference method)이다. 시각 $t_n$ 에서의 가속도로부터 속도와 변위를 반 스텝씩 어긋나게(leapfrog) 갱신한다.

\ddot{\mathbf{u}}_n = \mathbf{M}^{-1}\left( \mathbf{f}^{\text{ext}}_n - \mathbf{f}^{\text{int}}_n \right)

\dot{\mathbf{u}}_{n+1/2} = \dot{\mathbf{u}}_{n-1/2} + \Delta t \, \ddot{\mathbf{u}}_n, \qquad \mathbf{u}_{n+1} = \mathbf{u}_n + \Delta t \, \dot{\mathbf{u}}_{n+1/2}

여기서 결정적인 트릭이 하나 있다. 질량행렬 $\mathbf{M}$ 을 대각 행렬(집중질량, lumped mass)로 만들면 $\mathbf{M}^{-1}$ 이 그냥 각 성분의 역수라, 역행렬 계산이 사실상 나눗셈 몇 번으로 끝난다.² 강성행렬을 조립해 연립방정식을 푸는 과정 자체가 사라지는 것이다. 내부력 $\mathbf{f}^{\text{int}}$ 도 전역 행렬 없이 각 요소에서 응력을 적분해 절점력으로 모으는 식으로 요소 단위로 처리한다. 그 결과 한 스텝의 비용이 극도로 저렴하고, 병렬화도 자연스럽다. 이 구조는 입자 동역학에서 쓰는 베를레 적분과 수학적으로 형제지간이다.

3. 명시적 vs 암시적[편집]

명시적과 암시적의 차이는 구조해석자의 세계관을 가르는 근본적 분기점이다. 표로 정리하면 이렇다.

구분	명시적 (Explicit)	암시적 (Implicit)
스텝당 연립방정식	안 푼다	푼다 (뉴턴-랩슨법)
스텝당 비용	매우 저렴	비쌈
시간증분 크기	아주 작아야 함 (조건부 안정)	크게 가능 (무조건 안정)
강한 비선형·접촉	매우 강건	수렴 실패 잦음
적합한 시간 규모	밀리초 이하 초고속 현상	준정적·저속 현상

핵심 트레이드오프는 이거다. 명시적은 스텝 하나하나가 싸지만 스텝을 수십만 번 밟아야 하고, 암시적은 스텝이 비싸지만 큰 보폭으로 성큼성큼 간다. 그래서 충돌처럼 짧고 격렬한 현상은 명시적이, 몇 시간에 걸친 크리프 같은 준정적 현상은 암시적이 유리하다. 특히 명시적은 접촉·파단 같은 극악한 비선형에서도 방정식을 안 풀기 때문에 좀처럼 죽지 않는다는 게 최대 강점이다.

4. 임계 시간증분과 CFL 조건[편집]

명시적 해석의 아킬레스건은 조건부 안정성이다. 시간증분 $\Delta t$ 가 어떤 임계값을 넘으면 해가 즉시 발산한다. 이 임계값은 CFL 조건에 지배되며, 물리적으로는 “한 스텝 동안 응력파(stress wave)가 가장 작은 요소 하나를 건너뛰어서는 안 된다”는 뜻이다.

\Delta t \leq \Delta t_{\text{crit}} = \frac{L_{\min}}{c}, \qquad c = \sqrt{\frac{E}{\rho}}

여기서 $L_{\min}$ 은 가장 작은 요소의 특성 길이, $c$ 는 재료의 음속(탄성파 속도)이다. 이 식이 함의하는 바가 무섭다. 격자 안에 딱 하나 있는 아주 작은 요소가 전체 해석의 시간증분을 지배한다. 성형 해석에서 실수로 만든 찌그러진 요소 하나가 $\Delta t$ 를 100배 줄여 해석 시간을 며칠로 늘려버리는 참사가 흔하다.³

5. 질량 스케일링[편집]

CFL의 저주에서 벗어나는 실무적 꼼수가 질량 스케일링(mass scaling)이다. 임계 시간증분 $\Delta t_{\text{crit}} = L_{\min}/\sqrt{E/\rho}$ 을 보면, 밀도 $\rho$ 를 키우면 음속 $c$ 가 느려져 $\Delta t_{\text{crit}}$ 이 커진다. 그래서 문제가 되는 작은 요소에만 인위적으로 질량을 더해 시간증분을 늘리는 것이다. 공짜 점심은 없다. 질량을 늘리면 관성력이 커져 결과가 왜곡될 수 있으므로, 준정적 문제에서 운동에너지가 내부에너지의 5~10% 이하로 유지되는지 반드시 감시해야 한다.⁴ “질량 스케일링은 잘 쓰면 약, 못 쓰면 독”이라는 격언이 괜히 있는 게 아니다.

6. 활용 분야와 LS-DYNA[편집]

명시적 동해석이 진가를 발휘하는 무대는 짧고 격렬한 사건들이다.

충돌 해석: 자동차 정면·측면 충돌, 승객 안전(에어백, 더미) 시뮬레이션. NCAP 별점이 명시적 해석 위에서 결정된다.
금속 성형: 판재 프레스, 롤 성형, 단조. 큰 소성 변형과 접촉 해석, 마찰이 뒤엉킨 문제라 명시적이 제격이다.
낙하·충격: 스마트폰 낙하, 포장재 충격, 조류 충돌(bird strike).
폭발·관통: 방탄, 폭발 하중 해석. 극단적 변형과 파단, 자기 접촉(self-contact)이 난무하는 영역.

이 분야의 절대 강자가 LS-DYNA다. 1976년 로런스 리버모어 국립연구소의 존 할퀴스트(John Hallquist)가 개발한 DYNA3D에서 출발했으며, 원래는 핵무기 낙하 충격 해석용이었다.⁵ 이 밖에 Abaqus/Explicit, PAM-CRASH, RADIOSS 등이 경쟁한다. 명시적 코드가 다루는 격렬한 접촉·관통 문제는 개념적으로 충돌 감지와 맞닿아 있어, 게임 물리 엔진의 강체 동역학과도 사촌 관계다.

7. 관련 문서[편집]

“명시적”과 “암시적”은 정적/동적 여부가 아니라 시간적분 방식을 가리키는 말이다. 헷갈리는 초심자가 “정적 명시 해석”을 찾다가 혼란에 빠지는 일이 잦은데, 명시적은 본질적으로 관성이 있는 동적 해석이다. ↩
대각 질량행렬의 역행렬은 그냥 대각 성분의 역수다. 이 “역행렬 없는 역행렬” 트릭이 명시적 해석 속도의 90%를 책임진다. 반대로 정합질량(consistent mass)을 쓰면 이 마법이 깨지므로, 명시적 코드는 거의 항상 집중질량을 쓴다. ↩
그래서 명시적 해석자는 메시 품질에 병적으로 민감하다. “가장 작고 못생긴 요소가 나의 주말을 결정한다”는 건 과장이 아니라 CFL 조건의 직접적 귀결이다. ↩
LS-DYNA 후처리에서 GLSTAT의 운동에너지 곡선을 노려보는 것이 준정적 해석자의 국룰이다. 운동에너지가 튀는 순간, 그 결과는 물리가 아니라 질량 스케일링이 만든 환상이다. ↩
냉전의 산물이 오늘날 우리 자동차의 안전을 지키고 있는 셈이다. 무기 해석용으로 태어난 코드가 생명을 구하는 도구가 된 것은, 공학사에서 드물게 훈훈한 반전이다. ↩