15.32 레일리 소산 함수와 점성 마찰 모델링

1. 개요

레일리 소산 함수(Rayleigh dissipation function)는 라그랑주 역학의 틀 안에서 점성 마찰을 다루는 표준 도구이다. 본 절에서는 레일리 소산 함수의 정의와 점성 마찰의 모델링 방법을 다룬다.

2. 레일리 소산 함수

2.1 정의

레일리 소산 함수는 일반화 속도의 양정 이차 형식이다.

$\mathcal{R}(\mathbf{q}, \dot{\mathbf{q}}) = \frac{1}{2}\sum_{i,j}b_{ij}(\mathbf{q})\dot{q}_i\dot{q}_j = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{B}(\mathbf{q})\dot{\mathbf{q}}$

여기서 $\mathbf{B}(\mathbf{q})$ 는 대칭 양정(또는 양 반정) 감쇠 행렬이다.

2.2 일반화 감쇠력의 유도

소산 함수로부터 일반화 감쇠력은 다음과 같이 유도된다.

$Q_i^{\text{damp}} = -\frac{\partial \mathcal{R}}{\partial \dot{q}_i} = -\sum_{j}b_{ij}\dot{q}_j$

행렬 형식으로

$\mathbf{Q}^{\text{damp}} = -\mathbf{B}(\mathbf{q})\dot{\mathbf{q}}$

2.3 라그랑주 방정식

레일리 소산 함수가 있는 경우 라그랑주 방정식은 다음과 같다.

$\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_i} - \frac{\partial L}{\partial q_i} + \frac{\partial \mathcal{R}}{\partial \dot{q}_i} = Q_i^{\text{ext}}$

여기서 $Q_i^{\text{ext}}$ 는 외부 비보존력이다.

3. 점성 마찰의 종류

3.1 단일 관절 점성 마찰

가장 단순한 점성 마찰은 각 관절에 독립적으로 작용한다.

$\tau_{i}^{\text{damp}} = -\beta_i\dot{q}_i$

여기서 $\beta_i > 0$ 는 점성 마찰 계수이다.

3.2 결합 점성 마찰

관절 사이의 결합이 있는 점성 마찰은 일반화된 행렬로 표현된다.

$\boldsymbol{\tau}^{\text{damp}} = -\mathbf{B}\dot{\mathbf{q}}$

여기서 $\mathbf{B}$ 는 일반적으로 대각 행렬이지만, 결합이 있는 경우 일반 양정행렬일 수 있다.

3.3 자세 의존 점성 마찰

일부 응용에서는 점성 마찰 계수가 자세에 의존할 수 있다.

$\mathbf{B}(\mathbf{q})$

이는 더 복잡한 모형이지만 일부 정밀한 응용에서 필요하다.

4. 점성 마찰의 물리적 배경

4.1 윤활 마찰

기계 부품이 윤활제(기름, 그리스 등)로 윤활되어 있으면 마찰력은 속도에 비례하는 점성 마찰이 된다. 이는 윤활제의 유체 흐름에서 기인한다.

4.2 페트로프의 법칙

페트로프(Petroff)의 법칙은 윤활된 베어링의 마찰 토크를 다음과 같이 표현한다.

$\tau_{\text{petroff}} = \frac{2\pi\mu r^2 L}{c}\omega$

여기서 $\mu$ 는 점도, $r$ 은 베어링 반지름, $L$ 은 길이, $c$ 는 간극, $\omega$ 는 각속도이다. 이는 각속도에 비례한다.

4.3 모터의 역기전력

DC 모터의 경우 역기전력에 의한 효과가 점성 마찰처럼 작용한다.

5. 매니퓰레이터의 점성 마찰 모형

5.1 표준 모형

매니퓰레이터의 표준 점성 마찰 모형은 각 관절에 점성 마찰 계수를 부여한다.

$\boldsymbol{\tau}^{\text{damp}} = -\text{diag}(\beta_1, \beta_2, \ldots, \beta_n)\dot{\mathbf{q}}$

5.2 동역학 방정식

점성 마찰을 포함한 매니퓰레이터의 동역학 방정식은

$\mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{B}\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau}$

6. 매개 변수의 식별

6.1 정상 상태 시험

일정 각속도에서의 정상 상태 토크를 측정하여 점성 마찰 계수를 식별한다. 모든 다른 항(관성, 코리올리/원심력, 중력)을 보정한 후 잔류 토크가 점성 마찰이다.

6.2 자유 응답 시험

자유 응답에서 진동의 감쇠율을 측정하여 점성 마찰 계수를 결정한다.

6.3 회귀 분석

매니퓰레이터의 동역학 방정식의 선형 매개 변수화를 활용하여 점성 마찰 계수를 다른 매개 변수와 함께 식별할 수 있다.

7. 에너지 손실율

7.1 일률

점성 마찰에 의한 에너지 손실율은 다음과 같다.

$P_{\text{loss}} = \dot{\mathbf{q}}^T\mathbf{B}\dot{\mathbf{q}} = 2\mathcal{R}$

이는 항상 양수이며, 시스템이 에너지를 손실함을 의미한다.

7.2 시스템의 안정성

점성 마찰은 시스템의 안정성을 향상시킨다. 진동을 감쇠시키고 평형점으로의 수렴을 가속화한다.

8. 한계와 확장

8.1 단순 모형의 한계

단순 점성 마찰 모형은 모든 마찰 효과를 정확히 표현하지 못한다. 정지 마찰, 스트라이벡 효과, 비선형 마찰 등은 별도의 모형이 필요하다.

8.2 더 복잡한 모형

스트라이벡 모형, 르그르 모형 등 더 복잡한 마찰 모형이 정밀한 응용에서 사용된다.

9. 응용

9.1 정확한 동역학 시뮬레이션

점성 마찰을 포함하면 매니퓰레이터의 동적 시뮬레이션이 더 정확해진다.

9.2 정밀 제어

정밀 제어에서 점성 마찰의 보상이 추종 성능을 향상시킨다.

9.3 안정성 분석

점성 마찰이 매니퓰레이터의 안정성에 어떻게 기여하는지 분석한다.

10. 본 절의 의의

본 절은 레일리 소산 함수와 점성 마찰의 모델링을 다루었다. 이는 매니퓰레이터의 정확한 동역학 모형의 핵심 요소이며, 시뮬레이션과 제어에 활용된다.

11. 참고 문헌

Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
Lanczos, C. (1986). The Variational Principles of Mechanics (4th ed.). Dover.
Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
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.

version: 1.0