18.8 LuGre 마찰 모델 (LuGre Friction Model)

1. 개요

LuGre 마찰 모델은 Lund 대학과 Grenoble 대학의 공동 연구에서 1995년 Canudas de Wit, Olsson, Åström, Lischinsky에 의해 제안된 상태 변수 기반 동적 마찰 모델이다. 이 모델은 접촉면 사이의 미시적 강모(bristle)의 평균 변형을 내부 상태 변수로 도입하여, 사전 미끄럼, 고착-미끄럼, 스트리벡 효과, 마찰 이력, 마찰 지연과 같은 다양한 현상을 단일한 상미분 방정식 체계로 통합적으로 기술한다. LuGre 모델은 정적 마찰 모델의 한계를 극복하고 제어 설계에 적합한 미분 가능 구조를 제공한다는 점에서 로봇 공학과 정밀 제어에서 광범위하게 사용되어 왔다. 본 절은 모델의 물리적 동기, 수학적 구조, 수동성과 안정성, 매개변수의 의미, 식별 절차, 수치적 구현, 한계와 확장을 학술적으로 정리한다.

2. 물리적 동기

LuGre 모델의 출발점은 두 접촉 표면 사이의 미시적 접촉을 유연한 강모의 집합으로 모델링하는 것이다. 각 강모는 상대 변위에 따라 탄성적으로 변형되며, 충분한 변형에 도달하면 뒤틀림과 미끄럼이 발생한다. 거시적 마찰력은 이러한 강모들의 평균적 변형 상태에 의해 결정된다고 가정된다. 이 관점은 Dahl의 초기 연구에서 제시된 사전 미끄럼의 개념을 일반화한 것으로, 미시적 변형과 거시적 미끄럼을 단일한 연속 과정으로 통합한다. 강모 집단의 평균 변형을 하나의 상태 변수 $z$ 로 축약함으로써, 복잡한 미시 거동을 단순한 거시 방정식으로 치환한다.

3. 수학적 정식화

LuGre 모델의 내부 상태 변수 $z$ 는 다음의 상미분 방정식에 의해 지배된다.

$\dot{z} = v - \frac{|v|}{g(v)}z$

여기서 $v$ 는 접촉면의 상대 속도이며, $g(v)$ 는 정상 상태 마찰 곡선을 포착하는 양의 함수이다. 거시적 마찰력은 상태 변수와 그 도함수, 상대 속도에 의존하며 다음과 같이 표현된다.

$F_f = \sigma_0 z + \sigma_1 \dot{z} + \sigma_2 v$

여기서 $\sigma_0$ 은 강모의 강성, $\sigma_1$ 은 강모의 감쇠, $\sigma_2$ 는 점성 마찰 계수이다. 정상 상태 함수 $g(v)$ 는 전형적으로 다음과 같이 스트리벡 형태로 주어진다.

$g(v) = F_c + (F_s - F_c)\exp\!\bigl(-(v/v_s)^2\bigr)$

이 함수는 속도의 크기에 따라 쿨롱 마찰력과 정지 마찰력 사이에서 매끄럽게 전이되며, 스트리벡 효과를 자연스럽게 포함한다. 정상 상태에서 $\dot{z} = 0$ 이므로 $z_{ss} = g(v)\,\mathrm{sgn}(v)/\sigma_0$ 이 성립하며, 이에 대응하는 마찰력은 $g(v)\,\mathrm{sgn}(v) + \sigma_2 v$ 이다.

포착되는 현상

LuGre 모델은 여러 마찰 현상을 단일 구조로 포착한다. 사전 미끄럼 영역에서는 $|v|$ 가 매우 작고 $z$ 가 거의 선형으로 변하므로, 마찰력은 변위에 근사적으로 비례하는 탄성적 응답을 보인다. 정상 상태 미끄럼에서는 $z$ 가 $g(v)$ 에 의해 결정되는 값으로 수렴하며, 속도 의존적 마찰 곡선이 재현된다. 외부 힘이 주기적으로 진동하는 경우, $z$ 의 진폭 변화에 의해 자연스러운 마찰 이력 루프가 형성되며, 이는 마찰 기억(friction memory)과 마찰 지연(frictional lag)을 기술한다. 또한 고착-미끄럼 진동은 $g(v)$ 의 음의 기울기와 시스템 동역학의 상호 작용으로 자연스럽게 발생한다.

수동성과 에너지 해석

LuGre 모델은 특정 매개변수 조건에서 수동성(passivity)을 만족한다. 즉, 마찰이 상대 운동으로부터 에너지를 유입받지만 공급하지는 않는다는 물리적 조건이 $\sigma_0$ , $\sigma_1$ , $g(v)$ 사이의 관계에 의해 보장된다. 이러한 수동성은 LuGre 모델과 수동 제어기의 결합에서 폐회로의 에너지 관점의 안정성을 보장하는 데 유용하다. Barahanov와 Ortega는 LuGre 모델의 수동성 조건을 해석적으로 연구하였으며, Lischinsky 등은 수동성이 실험적 매개변수 범위에서 성립함을 보였다. 이러한 이론적 뒷받침은 LuGre 모델이 제어 설계에서 표준적 선택이 된 이유 중 하나이다.

매개변수의 의미

LuGre 모델의 매개변수는 각각 뚜렷한 물리적 의미를 가진다. $\sigma_0$ 은 사전 미끄럼 영역의 접선 강성으로, 작은 변위에 대한 마찰력의 기울기를 결정한다. $\sigma_1$ 은 사전 미끄럼의 점성 감쇠를 나타내며, 진동 거동에서 마찰 이력의 감쇠 특성을 결정한다. $\sigma_2$ 는 고속 영역의 점성 마찰 계수이다. $F_c$ 와 $F_s$ 는 각각 쿨롱 마찰과 정지 마찰에 해당하며, $v_s$ 는 스트리벡 특성 속도이다. 이러한 매개변수들의 합리적 선택은 특정 접촉 시스템의 거동을 정확히 재현하는 데 필수적이다.

매개변수 식별

LuGre 모델의 매개변수는 정적 식별과 동적 식별의 두 단계로 진행된다. 정적 식별에서는 일정 속도 실험을 통해 $F_c$ , $F_s$ , $v_s$ , $\sigma_2$ 를 추정한다. 이는 정상 상태 관계 $F_{ss}(v) = g(v)\,\mathrm{sgn}(v) + \sigma_2 v$ 에 대한 비선형 적합으로 수행된다. 동적 식별에서는 미소 진폭의 주기 운동 실험이나 사전 미끄럼 영역의 접선력-변위 이력 측정을 통해 $\sigma_0$ 과 $\sigma_1$ 을 추정한다. 사전 미끄럼의 선형 기울기가 $\sigma_0$ 을, 진동 이력의 감쇠 특성이 $\sigma_1$ 을 제공한다. Canudas de Wit과 Lischinsky는 이러한 식별 절차를 단계적으로 제시하였다.

수치적 구현

LuGre 모델의 수치적 구현은 상태 변수 $z$ 의 상미분 방정식을 통상의 적분기로 풀어내는 과정을 포함한다. 명시적 오일러 기법은 단순하지만, $g(v)$ 가 매우 작은 영역에서 $z$ 의 시간 상수가 작아져 수치 강성(stiffness)이 증가할 수 있다. 이러한 경우에는 암시적 방법이나 고차 적분기가 선호된다. 또한 $\sigma_0$ 이 매우 큰 경우에는 마찰력의 관측 방정식이 거의 대수 관계에 가까워져 수치 민감성이 증가하므로, 적분기의 선택과 단계 크기에 주의가 필요하다. 실시간 응용에서는 간소화된 적분이나 명시적 방법이 사용된다.

제어 설계와의 관계

LuGre 모델은 관측기 기반 마찰 보상 기법의 표준 기반으로 사용된다. 내부 상태 $z$ 는 직접 측정되지 않으므로, 속도 측정을 이용하여 $z$ 의 추정치 $\hat{z}$ 를 구동하는 관측기가 설계되며, 이 추정치를 이용하여 마찰력을 예측하고 제어 입력에 더하는 보상 구조가 일반적이다. Canudas de Wit과 Lischinsky는 이러한 관측기 기반 보상의 수렴성과 안정성을 분석하였으며, 이후의 많은 연구에서 로봇 관절의 저속 정밀 제어 성능 향상을 확인하였다. 적응 마찰 보상은 매개변수의 온라인 추정을 결합하여 작동 조건 변화에 대응한다.

한계와 확장

LuGre 모델은 강력한 포괄성에도 불구하고 몇 가지 한계를 가진다. 첫째, 극저속 영역에서 사전 미끄럼의 거동이 이력 곡선의 일부 세부 특성을 정확히 재현하지 못할 수 있으며, 이는 Dupont 등이 분석한 바 있다. 둘째, 수동성 조건이 매개변수 범위에 따라 엄격하게 성립하지 않는 경우가 있어, 이를 보장하는 수정된 버전이 제안되었다. 셋째, 온도와 윤활 상태의 변화에 따른 매개변수 변동을 포착하려면 확장된 정식화가 필요하다. 이러한 한계를 극복하기 위해 elasto-plastic 모델, generalized Maxwell-slip 모델, Leuven 모델과 같은 확장 모델들이 제안되었다.

본 절의 의의

본 절은 LuGre 마찰 모델을 체계적으로 정리하고, 상태 변수 기반 동적 마찰 모델의 대표 예로서 그 수학적 구조와 실무적 의의를 제시한다. LuGre 모델은 정적 마찰 곡선, 사전 미끄럼, 마찰 이력, 고착-미끄럼을 단일 체계로 통합하는 능력 덕분에 로봇 공학과 정밀 제어의 실무에서 표준적 선택이 되었다. 이 모델의 이해는 이후의 동적 마찰 모델과 이력 기반 마찰 기술의 학습에 필수적이다.

학습 권장사항

독자는 1자유도 질량 시스템에 LuGre 모델을 구현하고, 정상 상태 속도 스윕, 사전 미끄럼 변위 실험, 저속 고착-미끄럼 시뮬레이션을 차례로 수행해 볼 것을 권장한다. 각 실험에서 정상 상태 마찰 곡선, 사전 미끄럼 이력, 주기 운동의 마찰 지연 루프가 재현됨을 확인하면, 모델의 포괄성을 직접 체감할 수 있다. 또한 관측기 기반 마찰 보상을 단순 PI 제어와 결합하여 저속 궤적 추종 성능을 비교해 보는 것도 유익하다.

참고 문헌

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.
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., & Lischinsky, P. (1997). Adaptive friction compensation with partially known dynamic friction model. International Journal of Adaptive Control and Signal Processing, 11(1), 65–80.
Dupont, P., Hayward, V., Armstrong, B., & Altpeter, F. (2002). Single state elastoplastic friction models. IEEE Transactions on Automatic Control, 47(5), 787–792.
Barahanov, N., & Ortega, R. (2000). Necessary and sufficient conditions for passivity of the LuGre friction model. IEEE Transactions on Automatic Control, 45(4), 830–832.
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.
Al-Bender, F., Lampaert, V., & Swevers, J. (2005). The generalized Maxwell-slip model: a novel model for friction simulation and compensation. IEEE Transactions on Automatic Control, 50(11), 1883–1887.

version: 1.0