3차원 접촉점 추적을 통한 마찰 하중 모델링

초록

본 논문은 두 물체가 접촉할 때 접촉점의 3차원 이동 경로를 추적하여 슬라이드·롤·스핀 마찰 하중을 동시에 계산하는 현상학적 모델을 제시한다. 작은 시간 간격에서 접촉점이 표면을 따라 이동하는 아크를 정의하고, 접촉 프레임을 재정의함으로써 정적·동적 마찰계수를 모두 반영한다. 탄성·점성 성분을 분리해 에너지 소산을 제어하고, 롤링·스핀 마찰까지 포함한 3D 마찰 해석을 가능하게 한다.

상세 분석

이 연구는 기존 2차원 마찰 모델이 갖는 한계를 극복하고, 복잡한 형상의 물체 간 3차원 접촉을 정밀히 기술하기 위해 ‘접촉점 궤적 추적’이라는 새로운 접근법을 도입하였다. 핵심 아이디어는 시간 단계 Δt(10⁻⁶~10⁻³ s) 동안 각 물체 표면 위에 그려진 접촉점(크로스·점) 사이의 이동을 지오데식(곡면상의 최단 경로)으로 가정하고, 이를 s_i, s_j 라는 아크 길이로 정량화한다. 이때 s_i와 s_j의 차이는 상대 슬립(s)으로 정의되며, 슬라이드 마찰의 기본 구동량이 된다.

접촉 프레임(n, u, w)은 접촉점의 법선 n을 기준으로 접선 평면에 두 개의 직교 단위벡터 u, w를 배치한다. 초기 프레임(u₀, w₀)은 임의로 정하고, 다음 단계에서는 전 단계의 프레임과 현재 프레임 사이의 코사인 유사도를 최대화하는 최적화(선형 비용 함수)로 u₁, w₁을 결정한다. 이 과정은 전역 좌표계에서의 회전 변환을 최소화함으로써 프레임 간 ‘스핀 각 ψ’를 명확히 정의한다. ψ는 물체 j가 물체 i에 대해 회전한 각도로, 스핀 마찰 토크 계산에 직접 사용된다.

모델은 세 가지 마찰 메커니즘—슬라이드, 롤링, 스핀—을 에너지 소산 경로로서 독립적으로 취급한다. 이는 실제 접촉에서 발생하는 복합적인 마찰 현상을 단순화하지만, 각 메커니즘에 대한 탄성(K_E)과 점성(K_D) 파라미터를 별도로 설정함으로써 실험적 보정이 가능하도록 설계되었다. 슬라이드 마찰력은 탄성 성분 E_f = K_E·S_ij와 점성 성분 D_f = K_D·ΔS_ij·Δt의 합으로 표현되며, S_ij는 누적 슬립 변위(기억 효과)이다. 정적·동적 마찰계수 μ_s, μ_k에 따라 슬립 변위의 임계값 S_s, S_k = μ·N / K_E 를 설정하고, 현재 슬립 모드(s 또는 k)를 판단한다. 슬립이 임계값을 초과하면 탄성 변위를 스케일링(α)하여 슬립 모드 전환을 강제한다. 이러한 ‘캡핑’ 절차는 쿠롱 마찰법칙을 유지하면서도 미소 변형에 의한 탄성 복원을 반영한다.

롤링 마찰은 접촉점이 이동한 거리 |p_i|와 표면 곡률 κ_i, κ_j를 이용해 반경 R_i=1/|κ_i|, R_j=1/|κ_j|를 정의하고, 롤링 저항 토크 τ_r = K_E·Δθ·R·n 등으로 모델링한다(구체식은 본문에 상세히 제시). 여기서 Δθ는 한 시간 단계 동안의 회전 각도이며, 롤링 저항은 접촉면의 곡률 부호에 따라 토크 방향이 달라진다. 스핀 마찰은 앞서 정의한 ψ와 스핀 토크 K_E·ψ·n 으로 표현되며, 점성 성분 K_D·ψ·Δt·n을 추가한다. 세 메커니즘이 모두 독립적으로 계산된 후, 전체 마찰 하중은 벡터 합산되어 DEM(Discrete Element Method) 시뮬레이션에 바로 적용될 수 있다.

모델의 가정—작은 변위, 메커니즘 간 비결합, 미시 탄성에 의한 정적 마찰 제한—은 수치적 안정성을 확보하면서도 실제 물리 현상을 충분히 포착한다는 장점을 가진다. 특히, 스핀·롤링 마찰을 명시적으로 포함함으로써 회전 구동체(스핀 탑, 볼링 공, 회전식 베어링 등)의 동역학을 보다 정밀히 예측할 수 있다. 파라미터 K_D의 선택은 실험적 보정이 필요하지만, 저자들은 가상의 질량 m_ij와 K_E를 이용한 임계 감쇠식 K_D = 2√(m_ij·K_E) 를 기본값으로 제시한다.

전체적으로 이 모델은 복잡한 3D 접촉을 다루는 DEM 구현에 필요한 모든 마찰 요소를 일관된 수학적 프레임워크 안에 통합했으며, 기존 2D 기반 마찰 모델 대비 높은 일반성과 적용 범위를 제공한다.