🔢 mathematics

SoliDualSPHysics: An extension of DualSPHysics for solid mechanics with hyperelasticity, plasticity, and fracture

이 논문은 명시적인 균열 추적 없이 변형 가능한 고체의 초탄성, 유한 변형 소성 및 취성 파괴에 대한 통합된 대규모 시뮬레이션을 가능하게 하기 위해 전라그랑주(total Lagrangian) 정식화와 위상장(phase-field) 접근법을 활용하는 DualSPHysics의 새로운 오픈 소스 GPU 가속 확장판인 SoliDualSPHysics를 소개한다.

원저자: Mohammad Naqib Rahimi, George Moutsanidis

게시일 2026-06-23

📖 4 분 읽기🧠 심층 분석

CC BY 4.0

원저자: Mohammad Naqib Rahimi, George Moutsanidis

원본 논문은 CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

당신에게 거대하고 투명한 구슬 주머니가 있다고 상상해 보세요. 컴퓨터 시뮬레이션의 세계에서 이 구슬 주머니는 종-종 물의 흐름이나 기름 유출이 퍼지는 방식을 모델링하는 데 사용됩니다. 이는 구슬들이 서로 자유롭게 움직이고, 충돌하며, 엉키지 않고 서로를 압박할 수 있기 때문입니다(이는 전통적인 격자 기반 컴퓨터 모델에서 발생하는 '엉킨 실타래' 현상과는 다릅니다).

이 논문은 SoliDualSPHysics라고 불리는 강력하고 새로운 도구를 소개합니다. 이것은 그 "구슬 주머니"에 액체처럼 행동하는 법 대신 고체 물체(금속, 고무, 뼈와 같은)처럼 행동하는 법을 가르쳐주는 대대적인 소프트웨어 업그레이드라고 생각하면 됩니다.

다음은 이 소프트웨어가 하는 일을 쉬운 비유를 사용하여 설명한 것입니다.

1. 핵심 아이디어: 물에서 고체로

이전의 소프트웨어(DualSPHysics)는 유체를 시뮬레이션하는 데 탁월했습니다. 저자들은 이를 고체(늘어나고, 찌그러지고, 부서질 수 있는 물체)를 시뮬레이션할 수 있도록 만들고자 했습니다.

비유: 당신이 수영장을 시뮬레이션하고 있다고 상상해 보세요. 이제 당신은 고무줄이 끊어지거나 금속 빔이 휘어지는 것을 시뮬레이션하고 싶습니다. 이 새로운 소프트웨어는 "구슬"들이 유체처럼 자유롭게 움직이면서도, 동시에 고체 물체처럼 형태를 유지할 수 있게 해줍니다.

2. 늘어남과 굽힘을 처리하는 방법 (초탄성 및 소성)

고체는 당기거나 밀 때 서로 다르게 반응합니다.

초탄성 (고무줄): 고무와 같은 일부 재료는 늘어났다가 원래의 모양으로 완벽하게 되돌아옵니다. 이 소프트웨어는 이러한 "탱글탱글한" 동작을 시뮬레이션할 수 있습니다.
소성 (찰흙): 부드러운 금속이나 점토와 같은 다른 재료는 늘어난 상태 그대로 유지됩니다. 다시 돌아오지 않습니다. 이 소프트웨어는 재료가 영구적으로 변형되는 이 "영구적 변형" 현상을 이제 시뮬레이션할 수 있습니다.

3. 부서짐을 처리하는 방법 (파쇄)

이 부분이 가장 흥 есть 부분입니다. 과거에는 유리창에 금이 가는 것을 시뮬레이션하려면, 경주가 시작되기 전에 지도 위에 선을 긋는 것처럼 컴퓨터에게 정확히 어디에 금이 갈지를 미리 알려줘야 했습니다. 만약 금이 다른 곳으로 간다면, 시뮬레이션은 실패했습니다.

새로운 접근 방식 ("안개" 방법): 이 소프트웨어는 **위상장(Phase-Field)**이라 불리는 접근 방식을 사용합니다. 재료가 모든 곳에서 고체라고 생각하는 대신, 재료를 관통하는 얇은 "안개" 층이 있다고 상상해 보세요.
- 안개가 짙은 곳은 재료가 건강한 상태입니다.
- 안개가 옅은 곳은 재료가 손상된 상태입니다.
- 안개가 완전히 사라진 곳은 재료가 부서진 상태입니다.
멋진 이유: 당신은 컴퓨터에게 금이 어디로 갈지 알려줄 필요가 없습니다. 컴퓨터가 자연스럽게 스스로 찾아냅니다. 금은 실제 생활에서 그러하듯, 특별한 지시 없이도 시작되고, 자라나고, 두 갈래로 갈라지거나 다른 금과 합쳐질 수 있습니다.

4. 엔진: 속도와 힘

수백만 개의 구슬을 시뮬레이션하는 것은 컴퓨터에게 매우 힘든 일입니다.

비유: 백만 개의 당구공이 움직이는 경로를 계산한다고 상상해 보세요. 하나씩 계산한다면 시간이 엄청나게 오래 걸릴 것입니다.
해결책: 이 소프트웨어는 GPU(게이밍 컴퓨터에 들어있는 강력한 그래픽 카드)에서 실행되도록 구축되었습니다. 이는 수천 개의 작은 프로세서가 동시에 함께 작동하는 병렬 처리를 사용합니다. 이를 통해 시뮬레이션이 매우 빠르게 실행되어, 과학자들이 수백만 개의 입자를 가진 거대한 구조물을 합리적인 시간 내에 모델링할 수 있게 해줍니다.

5. 유연성: "리모컨"

이 소프트웨어에는 특별한 "리모컨"(XML이라는 언어로 작성됨)이 포함되어 있습니다.

기능: 당신은 수학 공식을 사용하여 특정 부분을 밀거나, 당기거나, 흔들도록 시뮬레이션에 명령할 수 있습니다. 심지어 빔의 왼쪽 부분에만 힘을 가하거나, 시뮬레이션의 첫 1초 동안에만 힘을 가하도록 설정할 수도 있습니다.
이점: 이를 통해 엔지니어들은 매번 프로그램 전체를 새로 작성하지 않고도 복잡한 테스트를 쉽게 설정할 수 있습니다.

아직 할 수 없는 것 (주의 사항)

저자들은 현재 이 도구의 한계에 대해 매우 솔직하게 밝히고 있습니다:

이 도구는 순수 금속이나 고무처럼 모든 방향에서 성질이 동일한(등방성) 재료에 가장 잘 작동합니다. 아직 탄소 섬유(한 방향으로 더 강한 재료)와 같은 복합 재료는 다루지 못합니다.
이 도구는 취성(유리나 세라믹처럼 깨지는 현상) 파쇄를 매우 잘 시뮬레이션합니다. 하지만 연성(금속 와이어가 늘어나다가 가늘어지며 끊어지는 현상) 파쇄는 아직 시뮬레이션하지 못합니다. 이는 다른 종류의 물리학이 필요하기 때문입니다.
현재 버전에서는 소성(영구적 굽힘)과 파쇄(부서짐)를 별개의 사건으로 취급합니다. 즉, 이 특정 버전에서는 재료가 구부러지면서 동시에 부서지는 현상을 시뮬레이션하지 않습니다.

요약

SoliDualSPHysics는 유체 시뮬레이션 프로그램을 강력한 고체 역학 시뮬레이터로 바꿔주는 무료 오픈 소스 도구입니다. 이 도구는 금이 자연스럽게 나타나도록 "안개" 방법을 사용하며, 현대적인 그래픽 카드에서 매우 빠르게 실행되고, 엔지니어가 극한 조건에서 재료가 어떻게 늘어나고, 휘고, 부서지는지 테스트할 수 있게 해줍니다. 이는 가상 실험을 위한 유연하고 고속인 실험실로 설계되었습니다.

기술 요약: SoliDualSPHysics

문제 정의

Smoored Particle Hydrodynamics (SPH)는 라그랑주적 특성 덕분에 유체 역학 및 대변형 고체 역학 분야에서 강력한 메시 프리(meshfree) 방법론으로 자리 잡았습니다. 그러나 현재의 고성능 오픈 소스 SPH 프레임워크 지형은 전산 유체 역학(CFD) 및 강체 또는 미소 변형 탄성을 포함하는 유체-구조 상호작용(FSI)에 크게 편중되어 있습니다. 고도화된 구성 방정식(constitutive behaviors), 구체적으로 유한 변형 소성(finite-strain plasticity)과 물리적으로 일관된 취성 파괴(brittle fracture)를 처리할 수 있는 통합 프레임워크 측면에서는 상당한 공백이 존재합니다.

기존의 SPH 파괴 정식화는 종종 임의적인 경험적 손상 모델, 입자 균열 접근법 또는 코헤시브 존 모델(cohesive zone models)에 의존하며, 이는 메시 의존성, 비수렴 결과, 또는 명시적인 균열 추적이나 국부적 보강(local enrichment)의 필요성 문제를 야기합니다. 또한, 이러한 고급 고체 역학 기능을 대규모 시뮬레이션에 적합한 확장 가능한 CPU/GPU 아키텍처로 통합하는 것 역시 과제로 남아 있습니다. 따라서 식별된 구체적인 필요성은 명시적인 균열 추적 없이도 초탄성(hyperelasticity), 유한 변형 소성 및 취성 파괴를 지원하면서, 병렬화를 통해 계산 효율성을 유지하는 통합된 SPH 정식화입니다.

방법론

저자들은 DualSPHysics 프레임워크의 오픈 소스 GPU 가속 확장판인 SoliDualSPHysics를 소개합니다. 이 방법론은 변형되지 않은 형상(undeformed configuration)에서 외부 하중과 경계 조건을 직접 적용할 수 있는 전체 라그랑주 SPH (Total Lagrangian SPH, TLSPH) 정식화에 기반합니다.

이론적 프레임워크

운동량 균형 (Momentum Balance): 솔버는 운동량 균형 방정식에 대한 강형(strong-form) TLSPH 이산화를 채택합니다. 일관성을 회복하기 위해 1차 커널 기울기 보정(kernel-gradient correction)을 사용하며, 수치적 불안정성과 가짜 제로 에너지 모드를 완화하기 위해 인공 점성 항을 포함합니다.
구성 모델링 (Constitutive Modeling):
- 초탄성 (Hyperelasticity): 본 프레임워크는 St. Venant–Kirchhoff 및 압축성 neo–Hookean 모델을 지원하며, 이는 취성 파괴를 시뮬레이션하기 위해 위상장(phase-field) 정식화와 결합됩니다.
- 유한 변형 소성 (Finite-Strain Plasticity): 등방성 경화가 포함된 rate-independent $J_2$ 탄소성 모델로의 확장이 구현되었습니다. 이는 변형 구배의 곱 분해( $F = F^e F^p$ )를 활용하며, 대변형 회전 상황에서도 프레임 객관성(frame objectivity)과 열역학적 일관성을 보장하기 위해 지수 맵 반환 알고리즘(exponential-map return algorithm)을 사용합니다.
- 결합 참고 (Coupling Note): 현재 구현에서 소성과 위상장 파괴는 서로 결합되어 있지 않습니다. 즉, 소성 모델은 파괴가 없는 탄소성 응답을 기술하며, 위상장 모델은 초탄성 재료에서의 취성 파괴를 처리합니다.
취성 파괴 (Brittle Fracture, Phase-Field): 파괴는 쌍곡형 위상장(hyperbolic phase-field) 정식화를 사용하여 모델링됩니다. 전역 선형 시스템을 풀어야 하는 타원형(elliptic) 정식화와 달리, 이 방식은 손상 변수 $s$ 에 대한 쌍곡형 진화 방정식을 사용합니다. 이를 통해 명시적인 추적이나 보강 없이도 균열의 생성, 전파, 분지(branching) 및 합체(coalescence)가 자연스럽게 나타나도록 합니다. 모델은 이력 함수(history functional)를 통해 비가역성을 강제하며, 탄성 변형 에너지를 분해함으로써 균열 전파가 인장 영역으로 제한되도록 합니다.
시간 적분 (Time Integration): 시스템은 Verlet 또는 Symplectic/Leapfrog와 같은 명시적 시간 적분 스킴을 사용하여 전진되며, 시간 단계는 파동 전파 및 가속도에 기초한 CFL 유형 조건에 의해 제한됩니다.

소프트웨어 아키텍처 및 구현

플랫폼: 코드는 C++ 및 CUDA로 작성되었으며, DualSPHysics의 네이티브 CPU (OpenMP) 및 GPU (CUDA) 병렬 아키텍처를 활용합니다.
입출력 (I/O): 설정은 XML 입력 파일을 통해 관리됩니다. 이 확장은 변형 구조물, 재료 표현식, 경계 조건 및 사용자 정의 표현식을 정의하기 위한 새로운 블록(<deformstrucs>, <mathexpressions>)을 도입합니다.
사용자 정의 표현식 (User-Defined Expressions): 핵심 기능 중 하나는 수학적 표현식 파서를 포함하여, 사용자가 내장된 변수(위치, 시간, 변위)와 표준 수학적/논리적 연산자를 사용하여 시간 및 공간 의존적 양(예: 속도장, 트랙션 하중)을 규정할 수 있도록 하는 것입니다.
경계 조건: 유연한 속도 및 트랙션 경계 조건을 전역적으로 또는 보조 기하 구조로 정의된 특정 표면 그룹에 적용할 수 있습니다.

주요 기여

통합 프레임워크: SoliDualSPHysics는 SPH 패러다임 내에서 초탄성, 유한 변형 소성 및 취성 파괴를 시뮬레이션할 수 있는 단일 오픈 소스 환경을 제공합니다.
SPH에서의 위상장 파괴: 명시적 균열 추적 없이도 효율적인 명시적 동역학을 가능하게 하여 균열 분지 및 합체를 포착할 수 있는 쌍곡형 위상장 정식화를 구현함으로써, ad hoc 기준 없이도 효율적인 계산을 실현했습니다.
유한 변형 소성: 지수 맵 반환 매핑을 갖춘 곱 분해 $J_2$ 소성 모델의 통합은 대변형 소성 문제에 대한 코드의 적용성을 확장합니다.
확장성: DualSPHysics의 CPU/GPU 인프라를 활용함으로써, 수백만 개의 입자가 포함된 대규모 시뮬레이션에서 상당한 가속을 달 수 있습니다.
유연성: 사용자 정의 표현식 파서와 유연한 경계 조건 정의를 포함함으로써 복잡한 하중 시나리오 및 연구 응용 분야에 대한 코드의 적응성을 높였습니다.
오픈 소스 가용성: 재현성을 보장하기 위해 전체 소스 코드, 문서 및 벤치마크 케이스를 공개 GitHub 저장소를 통해 무료로 제공합니다.

결과 및 검증

소프트웨어는 일련의 벤치마크 수치 문제 및 실험 데이터에 대해 검증되었습니다:

탄성 동역학 (Elasto-Dynamics): 캔틸레버 보(beam) 및 플레이트(2D 및 3D)의 자유 진동 실험은 해석적 해 및 참조 유한 요소법(FEM) 결과와 우수한 일치를 보였습니다. 시뮬레이션은 거의 보존적인 에너지 교환과 최소한의 수치적 소산을 보여주었습니다.
대변형 (Large Deformation): 대규모 분포 표면 트랙션을 받는 3D 캔틸레버 보는 기하학적 비선형성과 복잡한 경계 조건을 처리하는 능력을 확인시켜 주었습니다.
비틀림 동역학 (Torsional Dynamics): 3D 컬럼의 비틀림 사례는 프레임워크가 비틀림 파동 전파와 탄성 회복을 해결하는 능력을 검증하였으며, 결과는 참조 연구와 일치하였습니다.
파괴 (Fracture): 동적 균열 분지 벤치마크(인장 하의 노치 플레이트)는 확립된 위상장 문헌과 일치하는 방식으로 분지를 포함한 복잡한 파괴 패턴을 포착하는 모델의 능력을 입증했습니다.
성능: GPU 구현은 CPU 전용 실행에 비해 유의미한 가속을 제공하여, 대규모 시뮬레이션을 지원함을 보여주었습니다.

의의 및 범위

본 논문은 SoliDualSPHysics를 고성능 SPH 유체 솔버와 고급 고체 역학 능력 사이의 간극을 메우는 중요한 단계로 위치시킵니다. 이의 주요 의의는 취성 파괴 시뮬레이션에서 ad hoc 균열 추적의 필요성을 제거하고, 메시 프리 프레임워크 내에서 유한 변형 소성을 지원하는 통합된, 오픈 소스 기반의, GPU 가속 도구를 제공한다는 점에 있습니다.

저자들은 현재 릴리스의 한계를 명시적으로 언급하였습니다:

등방성(isotropic) 재료 응답(등방성 초탄성, 등방성 $J_2$ 소성)으로 제한됩니다.
소성과 파괴는 결합되어 있지 않으며, 연성 파괴(ductile fracture) 메커니즘은 포함되지 않았습니다.
이방성 구성 모델, 복합 재료 및 전용 아워글래스 제어(hourglass-control) 연산자는 아직 구현되지 않았습니다.

이러한 제한 사항에도 불구하고, 본 연구는 SPH 커뮤니티 내에서 고체 역학 및 파괴 모델링의 향후 발전을 위한 견고한 토대를 구축하며, 재현 가능하고 확장 가능한 플랫폼을 제공합니다.