15.31 산일 함수의 정의와 감쇠 모델
1. 개요
산일 함수(dissipation function), 또는 레일리 산일 함수(Rayleigh dissipation function)는 라그랑주 역학에서 점성 마찰과 같은 비보존 감쇠력을 표현하는 데 사용되는 스칼라 함수이다. 본 절에서는 산일 함수의 정의와 다양한 감쇠 모델을 다룬다.
2. 산일 함수의 정의
2.1 기본 정의
레일리 산일 함수 \mathcal{F}는 일반화 속도의 이차 함수로 정의된다.
\mathcal{F}(\mathbf{q}, \dot{\mathbf{q}}) = \frac{1}{2}\sum_{i,j=1}^{n}d_{ij}(\mathbf{q})\dot{q}_i\dot{q}_j = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{D}(\mathbf{q})\dot{\mathbf{q}}
여기서 \mathbf{D}(\mathbf{q})는 감쇠 행렬이며 일반적으로 양의 반정행렬이다.
2.2 일반화 감쇠력
산일 함수로부터 일반화 감쇠력은 다음과 같이 유도된다.
Q_i^{\text{diss}} = -\frac{\partial \mathcal{F}}{\partial \dot{q}_i}
행렬 형식으로
\mathbf{Q}^{\text{diss}} = -\mathbf{D}(\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{F}}{\partial \dot{q}_i} = Q_i^{\text{nc}}
여기서 Q_i^{\text{nc}}는 산일 외의 비보존력이다.
3. 감쇠 모델
3.1 점성 마찰
가장 단순한 감쇠 모델은 점성 마찰이다. 산일 함수는 다음과 같다.
\mathcal{F}_v = \frac{1}{2}\sum_{i}\beta_i \dot{q}_i^2
여기서 \beta_i는 각 관절의 점성 마찰 계수이다. 일반화 감쇠력은
Q_i^{\text{diss}} = -\beta_i\dot{q}_i
3.2 결합 점성 감쇠
관절 사이의 결합이 있는 점성 감쇠는 일반화된 감쇠 행렬로 표현된다.
\mathcal{F} = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{B}\dot{\mathbf{q}}
여기서 \mathbf{B}는 대칭 양정 또는 양 반정 행렬이다.
3.3 비선형 감쇠
일부 응용에서는 비선형 감쇠 모델이 사용된다. 그러나 비선형 감쇠는 산일 함수로 항상 표현되지 않는다.
3.4 쿨롱 마찰
쿨롱 마찰은 산일 함수로 직접 표현되지 않는다. 그러나 다음과 같이 단순화된 산일 함수를 사용할 수 있다.
\mathcal{F}_C = \sum_{i}\tau_{C,i}\vert\dot{q}_i\vert
이는 미분 불가능하지만 부분 미분의 의미에서 다음을 만족한다.
Q_i^{\text{Coulomb}} = -\tau_{C,i}\,\text{sgn}(\dot{q}_i)
4. 감쇠의 응용
4.1 매니퓰레이터의 관절 마찰
매니퓰레이터의 관절 마찰은 산일 함수로 자주 표현된다. 점성 마찰과 쿨롱 마찰의 결합이 일반적이다.
4.2 공기 저항
매니퓰레이터의 빠른 운동에서 공기 저항도 감쇠 효과를 가진다. 일반적으로 속도의 제곱에 비례한다.
F_{\text{air}} = -c_d v\Vert v\Vert
이는 점성 산일 함수에 직접 들어가지 않는다.
4.3 구조 감쇠
매니퓰레이터의 링크 자체에서 발생하는 구조 감쇠는 일반적으로 작지만 정밀한 모형에서 고려된다.
5. 산일 함수의 성질
5.1 양수성
이상적인 산일 함수는 일반화 속도가 0이 아닐 때 항상 양수이다.
\mathcal{F}(\mathbf{q}, \dot{\mathbf{q}}) > 0, \quad \text{if } \dot{\mathbf{q}} \neq 0
이는 감쇠가 항상 운동을 감속시킨다는 물리적 사실을 표현한다.
5.2 에너지 손실율
산일 함수는 에너지 손실율과 관련된다.
\frac{dE}{dt} = -2\mathcal{F} = -\dot{\mathbf{q}}^T\mathbf{D}\dot{\mathbf{q}}
매니퓰레이터에서 에너지가 항상 감소함을 보여준다.
5.3 자세 의존성
감쇠 행렬 \mathbf{D}(\mathbf{q})는 일반적으로 자세에 의존할 수 있다. 그러나 많은 응용에서는 자세에 독립적으로 가정한다.
6. 매니퓰레이터의 동역학 방정식
6.1 감쇠를 포함한 형식
감쇠를 포함한 매니퓰레이터의 동역학 방정식은 다음과 같다.
\mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{D}\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau}
추가의 \mathbf{D}\dot{\mathbf{q}} 항이 점성 감쇠의 효과이다.
6.2 안정성에의 영향
감쇠는 시스템의 안정성을 향상시킨다. 진동을 줄이고 평형점으로의 수렴을 가속화한다.
7. 매개 변수의 식별
7.1 측정
감쇠 매개 변수는 일반적으로 실험적 측정으로 식별된다. 일정한 속도에서의 정상 상태 토크 또는 자유 응답의 감쇠율 등이 사용된다.
7.2 식별의 어려움
감쇠 매개 변수는 시간에 따라 변할 수 있으며(마모, 온도 등), 정확한 식별이 어려울 수 있다.
8. 응용
8.1 정확한 시뮬레이션
매니퓰레이터의 정확한 동적 시뮬레이션은 감쇠 모형을 포함해야 한다.
8.2 정밀 제어
정밀 제어에서 감쇠의 정확한 모형화와 보상이 중요하다.
8.3 안정성 분석
감쇠는 시스템의 안정성에 기여하며 안정성 분석에 포함된다.
9. 본 절의 의의
본 절은 산일 함수의 정의와 감쇠 모델을 다루었다. 산일 함수는 라그랑주 역학에서 비보존 감쇠력을 표현하는 우아한 도구이며, 매니퓰레이터의 정확한 동역학 모형에 필수적이다.
10. 참고 문헌
- 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.
- Greenwood, D. T. (2003). Advanced Dynamics. Cambridge University Press.
version: 1.0