물질점법

시뮬레이션 고체역학 컴퓨터 그래픽스 마지막 수정: 2026-07-19 04:21:44

1. 개요[편집]

물질점법 Material Point Method
약칭	MPM
제안	Sulsky, Chen & Schreyer (1994)
계보	PIC → FLIP → MPM
표현	라그랑주 입자 + 오일러 배경격자
대표 응용	눈·모래·젤리, 산사태, 고속 충돌

유한요소법은 격자가 찌그러지면 죽고, SPH는 이웃 탐색과 경계에서 운다. 그래서 둘의 약점을 서로의 강점으로 메운 잡종이 태어났다.

물질점법(Material Point Method, MPM)은 연속체를 질량·속도·변형구배·응력을 들고 다니는 라그랑주 물질점(material point, 입자)의 집합으로 표현하고, 운동량 방정식은 매 스텝 새로 초기화되는 오일러 배경격자(background grid) 위에서 푸는 하이브리드 연속체 해석 기법이다. 1994년 설스키(Sulsky) 등이 FLIP-PIC 유체 계열의 입자-격자 전송 구조를 고체역학의 약형식(weak form)에 접목하면서 정식화되었다.¹

핵심 아이디어는 역할 분담이다. 물질 정보(질량, 변형 이력, 소성 상태)는 입자가 영구히 들고 다니므로 이류(advection) 오차와 격자 왜곡이 원천적으로 없다. 반면 공간 미분이 필요한 계산 — 응력의 발산, 즉 내력 계산과 운동량 갱신 — 은 격자에서 한다. 격자는 매 스텝 계산이 끝나면 버려진다. 다음 스텝에는 깨끗한 새 격자가 나온다. 그래서 아무리 변형이 커도 메시가 꼬여 자코비안이 음수가 되는 유한요소법의 고질병이 발생하지 않는다.

2. 한 타임스텝의 알고리즘[편집]

MPM 한 스텝은 다음 다섯 단계로 정리된다. 앞의 두 글자가 P2G(particle-to-grid), 뒤가 G2P다.

P2G 전송: 형상함수 $N_i(\mathbf{x}_p)$ 를 가중치로 입자의 질량과 운동량을 격자 노드로 뿌린다.

m_i = \sum_p N_i(\mathbf{x}_p)\, m_p, \qquad m_i \mathbf{v}_i = \sum_p N_i(\mathbf{x}_p)\, m_p \mathbf{v}_p

격자 내력 계산: 입자가 들고 있는 코시 응력 $\boldsymbol{\sigma}_p$ 로부터 노드 내력을 조립한다. 이 항이 사실상 유한요소의 $\mathbf{B}^{T}\boldsymbol{\sigma}$ 적분이며, 입자가 곧 적분점(quadrature point) 역할을 한다.

\mathbf{f}_i^{\text{int}} = -\sum_p V_p\, \boldsymbol{\sigma}_p \nabla N_i(\mathbf{x}_p)

격자 운동량 갱신: $\mathbf{f}^{\text{int}} + \mathbf{f}^{\text{ext}}$ 로 노드 속도를 전진시키고, 여기서 경계 조건(마찰, 충돌면)을 부과한다.
G2P 전송: 갱신된 노드 속도를 입자로 되돌린다. 속도구배 $\nabla \mathbf{v}_p = \sum_i \mathbf{v}_i \otimes \nabla N_i$ 도 함께 얻는다.
입자 갱신 및 격자 폐기: 변형구배와 위치를 갱신하고 격자를 초기화한다.

\mathbf{F}_p^{n+1} = \left(\mathbf{I} + \Delta t\, \nabla \mathbf{v}_p\right)\mathbf{F}_p^{n}

이 변형구배 $\mathbf{F}_p$ 가 MPM을 유체 전용 FLIP과 갈라놓는 결정적 장치다. 입자가 자기 변형 이력을 기억하므로, 탄성 복원력이든 소성 항복이든 응력과 변형률 관계를 그대로 얹을 수 있다.

3. 형상함수와 셀 교차 노이즈[편집]

초기 MPM은 유한요소법처럼 선형(삼선형) 형상함수를 썼는데, 여기서 유명한 병이 나온다. 선형 함수의 기울기 $\nabla N_i$ 는 셀 경계에서 불연속이다. 입자가 셀을 가로지르는 순간 내력이 튀어, 아무 일도 없어야 할 정지 상태에서 진동이 생긴다. 이것이 셀 교차 노이즈(cell-crossing instability)다.

처방은 형상함수를 매끄럽게 만드는 쪽으로 수렴했다.

GIMP(Generalized Interpolation MPM): 입자를 점이 아니라 유한한 부피를 가진 영역으로 보고 형상함수를 그 영역에서 적분해 매끄럽게 만든다.
B-스플라인 MPM: 2차·3차 B-스플라인 커널을 쓴다. 기울기가 연속이라 노이즈가 거의 사라진다. 대가는 지지 영역(support)이 넓어져 P2G/G2P 비용이 커지는 것. 오늘날 그래픽스 구현의 사실상 표준.
CPDI(Convected Particle Domain Interpolation): 입자 영역이 변형구배를 따라 함께 변형되게 해, 큰 전단 변형에서도 정확도를 유지한다.

전송 방식도 그대로 물려받았다. 순수 PIC 전송은 각운동량과 에너지를 잃어 물체가 축 처지고, 순수 FLIP 전송은 노이즈가 쌓인다. 그래서 요즘은 입자에 국소 아핀 속도장까지 실어 보내는 APIC나, 최소자승 이동(moving least squares) 구조로 전송과 응력 계산을 통합해 성능을 끌어올린 MLS-MPM이 널리 쓰인다.²

4. 재료 모델: 진짜 강점[편집]

MPM이 진짜 빛나는 지점은 구성방정식이다. 격자는 계산용 스크래치패드일 뿐이므로, 입자에 어떤 재료 모델을 얹느냐에 따라 같은 코드가 전혀 다른 물질을 시뮬레이션한다. 탄소성 재료는 변형구배를 탄성부와 소성부로 곱셈 분해해 다룬다.

\mathbf{F} = \mathbf{F}_E \mathbf{F}_P

탄성부 $\mathbf{F}_E$ 로 초탄성 재료 모델(네오-후크, 픽스드 코로테이티드 등)에서 응력을 계산하고, 항복조건을 넘으면 초과분을 $\mathbf{F}_P$ 로 넘긴다. 이 한 줄짜리 틀 위에서:

눈: 탄소성 + 경화/연화 파라미터. 2013년 스토마킨(Stomakhin) 등의 눈 MPM 논문이 디즈니 애니메이션에 쓰이면서 이 기법이 그래픽스 업계에 폭발적으로 퍼졌다.³
모래·지반: 마찰 재료의 드러커-프라거 모델 항복면. 산사태·토석류 해석의 주력.
젤리·고무: 소성 없이 초탄성만.
점성 유체·꿀: 응력을 압력 + 점성 항으로 두면 그냥 유체가 된다.
파단·절단: 항복 후 응집력을 떨어뜨려 재료를 분리시킨다.

서로 다른 재료의 입자를 같은 격자에 함께 뿌리면 연성(coupling)이 공짜로 해결된다. 눈 위에 젤리 덩어리를 떨어뜨릴 때 접촉 알고리즘을 따로 짤 필요가 없다. 노드에서 운동량이 합쳐지는 순간 이미 상호작용이 처리된다. 접촉 해석이 얼마나 지긋지긋한 작업인지 아는 사람에게는 이것만으로도 MPM을 쓸 이유가 된다.

5. 약점[편집]

물론 세상에 공짜는 없다.

수치적 응착(numerical cohesion): 격자가 자동 접촉을 해주는 그 성질이 반대로, 떨어져야 할 두 표면이 한 셀 안에 들어오면 억지로 붙어버리는 부작용을 낳는다. 미세한 균열이나 얇은 막 분리가 잘 표현되지 않는다.
격자 의존성: 계산이 배경격자에서 이뤄지므로 셀 크기가 해상도를 결정한다. 입자만 늘려도 물리 해상도는 오르지 않는다.
체적 잠김: 비압축에 가까운 재료에서 저차 요소와 같은 체적 잠김(volumetric locking)이 나타나, F-bar류의 보정이 필요하다.
비용: 셀당 입자 4~8개(2D) 또는 8개 이상(3D)을 깔고 매 스텝 두 번 전송하므로 무겁다. 대신 P2G/G2P가 국소 연산이라 GPU 컴퓨팅과 궁합이 좋아, 최근 구현은 대부분 GPU 위에 올라간다.
명시적 시간적분의 제약: 대개 명시적으로 풀기 때문에 CFL 조건에 준하는 시간 간격 제한을 받는다. 강성이 큰 재료를 다루려면 암시적 MPM으로 가야 하고, 그러면 매 스텝 비선형 방정식을 풀어야 한다.

6. 어디에 쓰이나[편집]

영화·게임 쪽에서는 눈, 모래, 진흙, 젤리, 파괴되는 물체를 한 솔버로 처리할 수 있다는 점 때문에 물리 엔진과 별개의 오프라인 파이프라인으로 자리 잡았다. 공학 쪽에서는 산사태·토석류·눈사태 같은 대변형 지반 문제, 발파와 관입 같은 고속 충돌 문제(명시적 동해석의 영역), 그리고 이산요소법으로 개별 입자를 다 세기엔 규모가 큰 분립체 유동에 적용된다. “격자가 찢어지지 않는다”는 단 하나의 장점이 이렇게 넓은 응용을 만들어낸 셈이다.⁴

7. 관련 문서[편집]

Sulsky, D., Chen, Z. & Schreyer, H. L. (1994). “A particle method for history-dependent materials”. Comput. Methods Appl. Mech. Engrg. 제목의 “history-dependent”가 핵심이다. 소성처럼 이력을 기억해야 하는 재료를 격자 왜곡 없이 다루겠다는 것이 처음부터의 목표였다. ↩
MLS-MPM은 전송 커널과 응력 적분을 하나의 최소자승 구조로 합쳐 상수 계수를 줄인 것으로, 100줄 남짓한 코드로 동작하는 데모가 공개되면서 “MPM은 어렵다”는 인식을 상당히 깨뜨렸다. 물론 그 100줄을 이해하는 건 별개의 문제다. ↩
Stomakhin, A. et al. (2013). “A material point method for snow simulation”. ACM TOG (SIGGRAPH). 논문 한 편이 특정 기법을 업계 표준으로 밀어 올린 대표 사례로 꼽힌다. 이후 모래·거품·머리카락·점탄성으로 응용이 줄줄이 이어졌다. ↩
정리하자면 MPM은 “라그랑주로 물질을 기억하고, 오일러로 미분을 계산한다”는 한 문장으로 요약된다. 이 문장을 이해하면 나머지는 형상함수 고르기와 재료 모델 갈아 끼우기다. 물론 그 나머지가 논문 수천 편이다. ↩