18.9 Dahl 마찰 모델

1. 개요

Dahl 마찰 모델은 1968년 Philip R. Dahl이 미사일 유도 시스템의 베어링 마찰을 기술하기 위해 제안한 상태 변수 기반 동적 마찰 모델이다. 이 모델은 사전 미끄럼 영역에서 관찰되는 접선력과 미세 변위 사이의 탄소성적 이력 관계를 단일 상태 변수의 상미분 방정식으로 체계화한 최초의 시도로 평가된다. Dahl은 구조 재료의 응력-변형 이력과 유사한 방식으로 마찰을 기술함으로써, 쿨롱 마찰 모델의 불연속성을 매끄러운 과정으로 대체하였다. 본 절은 모델의 물리적 발상, 수학적 구조, 정상 상태 거동, 매개변수의 의미, 한계와 후속 발전의 관계를 학술적으로 정리한다.

2. 물리적 발상

Dahl은 마찰을 구조적 관점에서 바라보았다. 두 표면 사이의 접촉점들은 외부 접선력이 가해질 때 미세한 탄성 변형을 겪으며, 접선력이 일정 한계를 초과하면 소성 변형과 미끄럼이 시작된다. 이는 금속의 응력-변형 관계가 초기에는 선형 탄성 구간을 이루다가 항복 이후 소성 구간으로 전이하는 구조와 유사하다. 이러한 유추는 마찰력이 단순한 속도의 함수가 아니라 접선 변위의 이력에 의존하는 상태 변수로 표현되어야 함을 시사하였으며, 이는 당시까지의 쿨롱 기반 정적 모델과는 근본적으로 다른 관점을 제공하였다.

3. 수학적 정식화

Dahl 모델의 핵심 방정식은 마찰력 F와 접선 변위 x 사이의 관계를 다음과 같이 기술한다.

\frac{dF}{dx} = \sigma\,\biggl(1 - \frac{F}{F_c}\,\mathrm{sgn}(v)\biggr)^{\alpha}

여기서 \sigma는 접선 강성, F_c는 쿨롱 마찰 계수와 법선력의 곱으로 주어지는 마찰 한계, \alpha는 응력-변형 곡선의 형태를 결정하는 지수, v = dx/dt는 상대 속도이다. 지수 \alpha는 일반적으로 \alpha \ge 1로 선택되며, \alpha = 1인 경우가 가장 단순한 표준 형태이다. 시간에 대한 미분으로 변환하면 다음의 상미분 방정식이 얻어진다.

\dot{F} = \sigma\,\biggl(1 - \frac{F}{F_c}\,\mathrm{sgn}(v)\biggr)^{\alpha}\,v

이 방정식은 마찰력이 정적 한계 F_c\,\mathrm{sgn}(v)를 향해 점진적으로 수렴하는 과정을 기술하며, 사전 미끄럼 영역의 탄성적 증가와 미끄럼 영역의 포화를 단일 구조로 통합한다.

4. 정상 상태와 과도 응답

정상 상태 미끄럼에서 \dot{F} = 0이 요구되므로, F_{ss} = F_c\,\mathrm{sgn}(v)가 성립한다. 이는 Dahl 모델이 고속 영역에서 쿨롱 모델로 수렴함을 의미하며, 이 모델은 스트리벡 효과를 직접 포함하지 않는다. 과도 응답에서는 마찰력이 접선 변위의 함수로 매끄럽게 증가하여 포화값에 접근하며, 이 과정이 사전 미끄럼과 돌파의 연속적 기술을 제공한다. 미세 진동이 가해지는 경우에는 부분 이력 루프가 자연스럽게 형성되며, 이는 마찰 소산과 이력 기억의 기본 기제를 포착한다.

5. 매개변수의 의미

Dahl 모델의 매개변수는 각각 뚜렷한 물리적 의미를 가진다. \sigma는 사전 미끄럼 영역의 접선 강성으로, 변위 대 마찰력의 초기 기울기를 결정한다. 이 값이 클수록 사전 미끄럼의 특성 변위가 작아지며, 전이가 날카롭게 나타난다. F_c는 미끄럼 상태의 최대 마찰력이며, 쿨롱 한계에 해당한다. 지수 \alpha는 전이의 형태를 결정하여, 큰 값은 곡선의 하부에서 거의 선형적인 증가와 상부에서 급격한 포화를 유발하며, \alpha = 1은 연속적이고 완만한 지수적 수렴을 제공한다. 이들 매개변수의 조합은 특정 접촉 시스템의 사전 미끄럼 거동을 재현하도록 선택된다.

6. 해석적 해와 단순 예제

\alpha = 1인 경우, 일정 방향의 미끄럼 (v > 0)에서 초기 조건 F(0) = 0에 대해 해석적 해가 존재한다.

F(x) = F_c\,\bigl(1 - e^{-\sigma x/F_c}\bigr)

이 해는 마찰력이 영에서 시작하여 F_c에 지수적으로 수렴함을 보이며, 특성 변위는 F_c/\sigma로 주어진다. 이 수식은 사전 미끄럼 영역의 매끄러운 증가를 간결하게 기술하며, Dahl 모델의 가장 명확한 수학적 예증이다. 역방향 미끄럼이 개시되면 초기 조건이 재설정되어 대칭적 응답이 얻어진다.

수치적 구현

Dahl 모델은 단일 상태 변수를 가지며 미분 방정식이 간단하므로 수치 구현이 용이하다. 명시적 오일러 기법이나 Runge-Kutta 기법이 모두 적용 가능하며, \sigma가 매우 큰 경우에 한하여 수치 강성이 나타날 수 있다. 이 경우에는 암시적 기법이나 작은 시간 단계를 사용하는 것이 권장된다. 초기 조건 F(0)의 설정은 시뮬레이션의 일관성을 위해 주의가 필요하며, 일반적으로 정적 상태의 외부 접선력과 평형을 이루는 값으로 초기화된다.

제어 설계와의 관계

Dahl 모델은 사전 미끄럼 영역의 정확한 기술이 요구되는 정밀 위치 제어 응용에서 활용되어 왔다. 모델 기반 피드포워드 보상은 궤적의 속도와 변위를 이용하여 마찰력을 예측하고, 이를 제어 입력에 더함으로써 저속 영역의 위치 오차를 감소시킨다. 또한 Dahl 모델은 관측기 설계의 기초가 되며, 내부 상태를 직접 측정할 수 없는 상황에서 마찰력을 추정하는 구조를 제공한다. Dahl 모델은 단순성 덕분에 실시간 응용에서 계산 부담이 작다는 장점을 가진다.

한계

Dahl 모델은 중요한 한계를 가진다. 첫째, 모델은 스트리벡 효과를 직접 포함하지 않으며, 저속 영역의 비단조적 마찰 감소를 기술하지 못한다. 둘째, 정상 상태에서 마찰력이 속도의 크기에 독립적으로 쿨롱 한계에 수렴하므로, 점성 마찰이나 속도 의존성이 기본 형태에서는 배제된다. 이를 보완하려면 모델의 출력에 점성 항을 추가로 결합해야 한다. 셋째, 모델은 마찰 기억과 이력의 일부 특징을 포착하지만, 고착-미끄럼 진동의 특정 거동을 완전히 재현하지 못하는 경우가 있다. 이러한 한계는 이후의 동적 마찰 모델이 Dahl 모델을 확장하는 방향으로 발전한 이유가 된다.

후속 발전과의 관계

Dahl 모델은 후속 동적 마찰 모델의 직접적 전구체로 평가된다. Canudas de Wit 등이 제안한 LuGre 모델은 Dahl의 상태 변수 기반 접근을 일반화하여 스트리벡 효과와 점성 성분을 통합하였으며, Dupont 등의 elasto-plastic 모델은 Dahl 모델의 소성 해석을 개선하여 사전 미끄럼의 순수 탄성 구간을 구분한다. Al-Bender 등의 generalized Maxwell-slip 모델은 여러 개의 Dahl 요소를 병렬로 결합하여 복잡한 이력 거동을 기술한다. 이러한 발전은 모두 Dahl의 기본 발상을 기반으로 구축된다.

본 절의 의의

본 절은 Dahl 마찰 모델을 체계적으로 정리하고, 상태 변수 기반 동적 마찰 모델의 역사적 출발점으로서 그 학술적 의의를 제시한다. Dahl 모델은 단순한 수학 구조에도 불구하고 사전 미끄럼과 이력을 연속적으로 기술하는 최초의 성공적 시도였으며, 이후 모든 동적 마찰 모델의 개념적 기반이 되었다. 이 모델의 이해는 마찰 모델링의 발전 맥락과 현대 모델의 구조적 선택을 이해하는 데 필수적이다.

학습 권장사항

독자는 \alpha = 1인 Dahl 모델의 해석적 해를 직접 구현하여, 다양한 \sigma 값에 대한 사전 미끄럼 곡선을 도시해 볼 것을 권장한다. 이어서 주기적 접선력 입력에 대해 마찰력의 이력 루프가 형성되는 과정을 시뮬레이션하면, 상태 변수 기반 모델이 단순한 정적 모델과 어떻게 다른지를 정량적으로 이해할 수 있다. 또한 Dahl 모델에 점성 항을 추가한 확장을 구현하여, 후속 모델로의 자연스러운 진화 과정을 체감해 볼 것을 권장한다.

참고 문헌

  • Dahl, P. R. (1968). A Solid Friction Model. Technical Report TOR-0158(3107-18)-1, The Aerospace Corporation, El Segundo, CA.
  • Dahl, P. R. (1976). Solid friction damping of mechanical vibrations. AIAA Journal, 14(12), 1675–1682.
  • Armstrong-Hélouvry, B., Dupont, P., & Canudas de Wit, C. (1994). A survey of models, analysis tools and compensation methods for the control of machines with friction. Automatica, 30(7), 1083–1138.
  • Olsson, H., Åström, K. J., Canudas de Wit, C., Gäfvert, M., & Lischinsky, P. (1998). Friction models and friction compensation. European Journal of Control, 4(3), 176–195.
  • Canudas de Wit, C., Olsson, H., Åström, K. J., & Lischinsky, P. (1995). A new model for control of systems with friction. IEEE Transactions on Automatic Control, 40(3), 419–425.
  • Bliman, P.-A. (1992). Mathematical study of the Dahl’s friction model. European Journal of Mechanics, A/Solids, 11(6), 835–848.
  • Dupont, P., Hayward, V., Armstrong, B., & Altpeter, F. (2002). Single state elastoplastic friction models. IEEE Transactions on Automatic Control, 47(5), 787–792.

version: 1.0