접촉 해석

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

접촉 해석 Contact Analysis
분야	고체역학 × 비선형 구조해석
핵심 조건	비침투 조건 (KKT 조건)
대표 기법	페널티법, 라그랑주 승수법, 확대 라그랑주법
마찰 모델	쿨롱 마찰

1. 개요[편집]

두 물체가 서로를 통과하지 못한다. 이 당연한 사실을 수식으로 우기는 것이 접촉 해석이다.

접촉 해석(Contact Analysis)은 서로 맞닿는 두 개 이상의 물체가 하중을 받으며 힘을 주고받을 때, 그 접촉면에서의 응력·변형·미끄러짐을 계산하는 비선형 구조해석의 한 분야이다. 핵심은 “두 물체는 서로 파고들 수 없다” 는 지극히 상식적인 물리 조건을, 미리 알 수 없는 접촉 영역과 접촉력이라는 미지수와 함께 유한요소법 방정식에 강제로 집어넣는 데 있다.

문제는 이게 말처럼 쉽지 않다는 거다. 해석을 시작하기 전에는 두 물체가 어디에서, 얼마만큼의 면적으로 닿을지 아무도 모른다. 접촉 영역 자체가 미지수인 데다, 접촉하느냐 떨어지느냐에 따라 강성이 확 바뀌는 경계 비선형성(boundary nonlinearity)이 생긴다. 그래서 접촉 해석은 비선형 구조해석 중에서도 유독 수렴이 안 되기로 악명 높다.¹

2. 비침투 조건[편집]

접촉 해석의 지배 조건은 수학적으로 아주 깔끔하게 쓸 수 있다. 두 면 사이의 간극(gap)을 $g$ , 접촉 방향의 수직 접촉력(압력)을 $p_N$ 이라 하면 다음 세 조건이 동시에 성립해야 한다.

g \geq 0, \quad p_N \leq 0, \quad p_N \, g = 0

풀어 말하면 이렇다. (1) 간극은 음수가 될 수 없다(파고들면 안 된다). (2) 접촉력은 잡아당길 수 없고 밀어내기만 한다(접착제가 아니다). (3) 둘 중 하나는 반드시 0이다 — 떨어져 있으면( $g>0$ ) 접촉력이 0이고, 접촉력이 있으면( $p_N<0$ ) 간극이 0이다. 이 조합을 최적화 이론에서는 카루쉬-쿤-터커(KKT) 조건이라 부른다.² 즉 접촉 문제는 본질적으로 부등식 제약이 걸린 최소화 문제이며, 이 상보성(complementarity) 구조가 접촉 해석을 비선형·비평활(non-smooth)하게 만드는 근원이다.

3. 두 진영: 페널티법 vs 라그랑주 승수법[편집]

이 부등식 제약을 유한요소법 정식화에 어떻게 녹여 넣느냐를 두고 두 가지 대표 접근이 경쟁한다.

3.1. 페널티법 (Penalty Method)[편집]

물체가 살짝 파고드는 것을 허용하되, 파고든 침투량 $g$ 에 비례해 강력한 반발력을 주는 방식이다.

p_N = \varepsilon_N \, g \quad (g < 0)

여기서 $\varepsilon_N$ 이 접촉 강성(penalty stiffness)이다. 침투부를 아주 뻣뻣한 스프링으로 채운다고 생각하면 된다. 장점은 미지수가 늘지 않아 구현이 간단하고 방정식 구조가 깨지지 않는다는 것. 단점은 두 가지다. 접촉 강성이 너무 작으면 물체가 서로 뭉텅이로 파고들어 비물리적이고, 너무 크면 강성행렬의 조건수가 나빠져 수렴이 박살난다. 결국 “적당한 값”을 찾는 것이 사용자의 몫으로 남는데, 이 튜닝이 접촉 해석 초심자를 좌절시키는 주범이다.

3.2. 라그랑주 승수법 (Lagrange Multiplier Method)[편집]

접촉력 자체를 새로운 미지수(라그랑주 승수)로 도입해, 비침투 조건 $g=0$ 을 정확히 만족시키는 방식이다. 침투가 원천적으로 0이라 물리적으로 깔끔하다. 대신 미지수가 늘어나고, 강성행렬의 대각 성분에 0이 들어가 안장점(saddle point) 구조가 되면서 방정식 풀기가 까다로워진다.³ 두 방법의 장점을 절충한 것이 확대 라그랑주법(Augmented Lagrangian)으로, 페널티 항으로 대략 잡고 승수를 반복 갱신해 침투를 허용 오차 이내로 밀어 넣는다. 상용 코드(Abaqus, ANSYS)가 기본값으로 선호하는 방식이다.

4. 마찰: 쿨롱 모델[편집]

수직 방향 접촉이 정리되면 다음 골칫거리는 접선 방향, 즉 마찰이다. 가장 널리 쓰이는 것은 고전적인 쿨롱 마찰(Coulomb friction) 모델이다.

\|\mathbf{t}_T\| \leq \mu \, |p_N|

접선력 $\mathbf{t}_T$ 의 크기가 마찰계수 $\mu$ 와 수직 접촉력의 곱 이하이면 두 면은 붙어서 안 움직이고(stick), 이 한계에 도달하면 미끄러진다(slip). 문제는 이 stick-slip 전환이 또 하나의 상보성 조건이라는 점이다. 수직 접촉이 “붙었나 떨어졌나”를 판정하는 것도 모자라, 접선 방향으로 “붙었나 미끄러지나”까지 판정해야 하니 비선형성이 이중으로 겹친다. 게다가 쿨롱 마찰은 접선력-변위 관계가 비평활해서 뉴턴-랩슨법의 이차 수렴성을 갉아먹는다.⁴

5. 수치적 이산화: Node-to-Surface vs Mortar[편집]

접촉면을 이산화하는 방식도 정확도에 크게 영향을 준다. 전통적인 점-대-면(Node-to-Surface) 방식은 한쪽 면의 절점이 상대 면의 요소를 뚫는지 검사하는데, 구현이 쉽지만 접촉 압력이 절점에 집중되어 진동하고, 두 면의 절점이 어긋나면 압력 전달이 불연속해진다. 이를 개선한 것이 모르타르(Mortar) 방법으로, 접촉 조건을 절점이 아니라 면에 걸쳐 약형식(weak form)으로 적분해 부과한다. 매끈한 압력 분포와 격자 비정합에 대한 강건성 덕분에 현대 상용 코드가 점점 이쪽으로 이동하고 있다.

6. 왜 이렇게 수렴이 안 되는가[편집]

접촉 해석 실무자라면 누구나 잔차 그래프가 튀는 광경을 지겹도록 본다. 근본 원인을 정리하면 다음과 같다.

상태 급변: 접촉/분리, stick/slip이 사이클마다 켜졌다 꺼지면 강성이 계단처럼 점프해 뉴턴-랩슨법이 진동한다.
채터링(chattering): 한 절점이 접촉과 분리를 매 반복마다 오가며 해가 정착하지 못하는 현상. 접촉 안정화(stabilization)나 감쇠를 넣어 달래야 한다.
초기 접촉 부정확: 시작 시 두 물체가 붕 떠 있으면 강체 운동(rigid body motion)이 발생해 방정식이 특이(singular)해진다. 초기 간섭이나 접촉 안정화로 잡아준다.
하중 스텝 과대: 한 번에 너무 큰 하중을 걸면 접촉 영역이 폭증한다. 하중을 잘게 나눠(incrementation) 조금씩 접촉을 늘려가는 것이 정석.

결국 접촉 해석은 물리 이해도 중요하지만, 접촉 강성·안정화·하중 증분 같은 수치 파라미터를 얼마나 잘 어르고 달래느냐가 승패를 가르는 경우가 많다. 그래서 현업에서 접촉 해석은 “예술의 경지”라는 자조 섞인 농담이 붙어 다닌다.

7. 관련 문서[편집]

접촉 해석 사용자들 사이엔 “접촉이 수렴하면 그날은 로또를 사라”는 농담이 있다. 물론 실제로는 격자·강성·증분을 제대로 잡으면 잘 수렴한다. 다만 그 “제대로”에 도달하기까지의 삽질이 길 뿐. ↩
그래서 접촉역학을 최적화 이론과 변분부등식(variational inequality)의 언어로 재정식화하는 연구 분야가 따로 존재한다. “부딪히지 마”라는 유치원생도 아는 규칙이 수학적으로는 꽤 깊은 물이다. ↩
대각에 0이 들어간 시스템은 일반 촐레스키 분해가 통하지 않아 특수한 풀이 전략이 필요하다. 이 안장점 구조 때문에 라그랑주 승수법은 이론적으로 우아해도 대규모 문제에서 구현이 성가시다. ↩
쿨롱 마찰의 “정지↔미끄러짐” 전환점은 미분 불가능한 꺾인 점이라, 매끈한 함수를 가정하는 뉴턴법이 이 근처에서 헛발질을 한다. 그래서 실무에선 이 꺾임을 살짝 둥글린 정규화(regularized) 마찰 모델을 몰래 쓰기도 한다. ↩