하이드로코드

구조해석 고체역학 시뮬레이션 마지막 수정: 2026-07-18 04:43:27

1. 개요[편집]

하이드로코드 Hydrocode
분야	고속 충격 · 폭발 · 관통 해석
핵심 발상	고압에서 고체를 유체처럼 취급
필수 요소	상태방정식(EOS), 인공점성
이산화	라그랑주 / 오일러 / SPH
대표 코드	LS-DYNA, AUTODYN, CTH, ALE3D

총알이 강판을 뚫는 순간, 강철은 더 이상 고체가 아니라 걸쭉한 유체다.

하이드로코드(hydrocode)는 고속 충격·폭발처럼 압력이 재료의 강도를 압도하는 상황에서, 고체를 사실상 유체처럼 다루어 충격파(shock wave) 전파를 푸는 수치해석 코드다. 이름의 “하이드로”가 바로 그 발상을 담고 있다. 총알의 관통, 성형 작약탄(shaped charge)의 제트, 운석 충돌, 폭약의 폭굉(detonation) 같은 문제가 주 무대다.¹

일상적인 구조해석이 재료의 전단강성과 탄성계수를 중심에 둔다면, 하이드로코드의 세계에서는 압력이 수 기가파스칼(GPa)에 이르러 재료의 항복강도(수백 MPa)를 압도한다. 이 극한에서 물질의 거동을 지배하는 것은 강성이 아니라 압축률과 충격 특성이며, 그래서 압축성 유동의 방법론이 고체역학으로 넘어온다.

2. 왜 고체를 유체로 다루나[편집]

응력 텐서는 항상 두 부분으로 쪼갤 수 있다. 부피 변화를 담당하는 정수압(hydrostatic pressure) 성분과, 형상 변화를 담당하는 편차응력(deviatoric stress) 성분이다.

\boldsymbol{\sigma} = -p\,\mathbf{I} + \mathbf{s}

여기서 편차응력 $\mathbf{s}$ 의 크기는 재료의 항복강도로 제한된다. 문제는 고속 충격에서 정수압 $p$ 가 이 항복강도보다 수십~수백 배 커진다는 것이다. 상대적으로 편차응력을 무시할 만해지면, 응력 상태는 사실상 $\boldsymbol{\sigma} \approx -p\,\mathbf{I}$ , 즉 압력만 있는 유체와 같은 상태로 근사된다. “하이드로”라는 이름은 이 근사에서 나왔다. 실제 하이드로코드는 압력은 상태방정식으로, 편차응력은 강도 모델로 따로 다루어 둘을 합친다 — 강도를 끄면 순수 유체, 켜면 탄소성 고체가 되는 구조다.

편차응력을 담당하는 강도 모델(strength model)은 유체 근사를 고체로 되돌리는 보정항이다. Johnson-Cook 모델처럼 변형률·변형률 속도·온도에 따라 항복강도가 변하는 관계식을 써서, 충격 초반의 과압 국면과 충격이 지나간 뒤의 잔류 변형·파단을 모두 담는다. 특히 인장 방향의 압력이 임계값을 넘으면 재료가 내부에서 찢어지는 스폴(spall) 파괴가 일어나는데, 관통·폭발에서 뒷면이 떨어져 나가는 배면 박리(scabbing)가 그 결과다. 순수 유체 근사만으로는 이런 인장 파괴를 놓치므로, 강도·손상 모델이 EOS와 짝을 이루어야 실제 관통 현상이 재현된다.

3. 상태방정식(EOS)[편집]

하이드로코드의 심장은 상태방정식(Equation of State, EOS)이다. 압력을 밀도(또는 비체적)와 내부에너지의 함수로 주는 관계식으로, 충격 상태의 압력-압축 관계를 결정한다. 일반적 형태는

p = p(\rho, e)

이다. 충격 압축에서는 질량·운동량·에너지 보존이 충격면을 가로질러 성립해야 하며, 이를 정리한 것이 랭킨-위고니오(Rankine-Hugoniot) 관계다. 많은 금속에서 충격파 속도 $U_s$ 와 입자 속도 $u_p$ 가 선형 관계

U_s = c_0 + s\,u_p

를 따른다는 실험 사실을 바탕으로, Mie-Grüneisen EOS가 널리 쓰인다. 폭약에는 폭굉 생성물의 팽창을 기술하는 JWL(Jones-Wilkins-Lee) EOS가 표준이다. EOS를 잘못 고르면 압력이, 나아가 충격파 전체가 틀어지므로 하이드로코드에서 EOS 선정은 재료 모델 선정만큼 중대하다.

4. 충격파와 인공점성[편집]

압축성 유동과 마찬가지로, 하이드로코드가 마주하는 근본 난제는 충격파의 불연속이다. 이상적 충격면은 두께가 0인 수학적 불연속이라, 격자로 이산화하면 그 지점에서 수치 진동(oscillation)이 폭발한다.

두 갈래의 처방이 있다.

충격 포착(shock capturing) — 오일러 서술에서 리만 솔버와 고해상도 기법(TVD, WENO 등)으로 불연속을 안정적으로 잡아낸다. 압축성 유동 CFD의 유산을 그대로 가져온다.
인공점성(artificial viscosity) — 폰 노이만과 릭트마이어가 1950년 도입한 고전적 기법으로, 충격면에 인위적인 점성 압력을 더해 불연속을 몇 개 격자에 걸쳐 부드럽게 퍼뜨린다.² 라그랑주·SPH 코드에서 지금도 표준으로 쓰인다.

인공점성 계수는 너무 작으면 진동이 남고 너무 크면 충격이 뭉개지는, 전형적인 튜닝 지옥의 대상이다.

5. 라그랑주 · 오일러 · SPH[편집]

하이드로코드는 격자와 물질의 관계를 어떻게 잡느냐에 따라 세 갈래로 나뉜다.

라그랑주 — 격자가 물질을 따라 움직인다. 물질 경계와 이력(history) 추적이 정확하지만, 극단적 변형에서 격자가 뒤틀려 죽는다.
오일러 — 격자를 공간에 고정하고 물질이 그 사이로 흐른다. 아무리 찢겨도 격자가 안 죽지만, 물질 경계가 격자에 번지고 이력 추적이 어렵다.
SPH — 격자 자체를 버리고 입자로 물질을 표현한다. 파편 비산, 극단적 파쇄, 자유표면 문제에 강하다.

이 셋의 장단이 상보적이라, 현대 하이드로코드는 ALE 기법(격자와 물질을 부분적으로 분리)이나 라그랑주-오일러 커플링으로 절충한다. 관통체는 라그랑주로, 튀는 파편은 SPH로, 폭발 기체는 오일러로 — 한 해석에서 서술을 섞는 것이 요즘의 국룰이다.

6. 대표 코드와 활용[편집]

하이드로코드는 별도의 소프트웨어 장르라기보다, 위 기능(EOS·인공점성·다중 서술)을 갖춘 솔버를 부르는 이름에 가깝다.

LS-DYNA — 명시적 크래시 코드지만 EOS와 SPH·ALE를 갖춰 폭발·고속 충격에서 하이드로코드로 동작한다.
AUTODYN / ALE3D / CTH / Sandia의 코드들 — 국방·충격물리 연구에서 전통적으로 하이드로코드로 불려온 계보.

활용 분야는 태생부터 격렬하다. 방탄·관통 해석, 성형 작약탄과 관통자 설계, 폭발 하중과 폭풍파, 운석·우주 파편 충돌, 그리고 파괴역학이 감당 못 하는 극단적 파쇄까지. 공통점은 하나다 — 재료가 강도를 잊고 유체처럼 흐르는 찰나를, 컴퓨터로 붙잡아 두는 일이다.³

7. 관련 문서[편집]

“hydrocode”라는 표현은 냉전기 미국 무기연구소(로스앨러모스·리버모어·샌디아)에서 굳어진 관용어다. 초기 핵무기 설계 코드들이 고압 유체역학 계산에서 출발했기 때문에, 강도 항이 들어간 뒤에도 이름은 “하이드로”로 남았다. ↩
von Neumann, J. & Richtmyer, R. D. (1950). “A Method for the Numerical Calculation of Hydrodynamic Shocks.” 컴퓨터 여명기에 나온 이 논문이 오늘날의 모든 명시적 충격 해석의 조상이다. 물리를 정확히 푸는 대신 “진동이 안 나게 문지른다”는 실용주의가 70년을 살아남았다. ↩
재미있게도, 이 극한을 다루는 사람들 사이에서는 “고압에서는 물질의 종류가 별로 안 중요하다”는 말이 돈다. 충분히 세게 때리면 강철이든 알루미늄이든 결국 EOS로 기술되는 압축성 유체로 수렴하기 때문. 강도 모델은 충격 초반과 후반의 잔류 거동에서나 빛을 본다. ↩