15.24 다자유도 시스템의 라그랑주 운동 방정식
1. 개요
다자유도 시스템의 라그랑주 운동 방정식은 단일 자유도 시스템의 일반화이다. 매니퓰레이터를 비롯한 대부분의 실제 시스템은 다자유도이므로 이러한 방정식의 유도와 분석이 핵심이다. 본 절에서는 다자유도 시스템의 라그랑주 운동 방정식을 다룬다.
2. 다자유도 시스템의 형식
2.1 일반화 좌표 벡터
자유도가 n인 시스템의 일반화 좌표는 벡터로 표현된다.
\mathbf{q} = (q_1, q_2, \ldots, q_n)^T
2.2 라그랑지언
라그랑지언은 일반화 좌표와 일반화 속도의 함수이다.
L(\mathbf{q}, \dot{\mathbf{q}}, t) = T(\mathbf{q}, \dot{\mathbf{q}}, t) - U(\mathbf{q}, t)
2.3 운동 에너지의 일반 형식
스클레로노믹 시스템의 운동 에너지는 일반화 속도의 이차 형식이다.
T(\mathbf{q}, \dot{\mathbf{q}}) = \frac{1}{2}\sum_{i,j=1}^{n}M_{ij}(\mathbf{q})\dot{q}_i\dot{q}_j = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{M}(\mathbf{q})\dot{\mathbf{q}}
3. 라그랑주 방정식
3.1 각 좌표에 대한 방정식
각 일반화 좌표 q_i에 대해 다음의 라그랑주 방정식이 성립한다.
\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_i} - \frac{\partial L}{\partial q_i} = Q_i
여기서 Q_i는 일반화 힘이다. 총 n개의 방정식이 있다.
3.2 시스템의 상태
시스템의 상태는 (\mathbf{q}, \dot{\mathbf{q}})의 2n차원 벡터로 표현된다. 라그랑주 방정식은 이를 결정하는 미분 방정식이다.
4. 운동 방정식의 명시적 형식
4.1 미분의 적용
라그랑주 방정식을 운동 에너지와 위치 에너지로 분해한 후 미분을 적용하면 다음의 형식이 얻어진다.
\sum_{j=1}^{n}M_{ij}(\mathbf{q})\ddot{q}_j + \sum_{j,k=1}^{n}c_{ijk}(\mathbf{q})\dot{q}_j\dot{q}_k + \frac{\partial U}{\partial q_i} = Q_i
여기서 c_{ijk}는 크리스토펠 기호의 1종이다.
4.2 행렬 형식
행렬 형식으로 표현하면
\mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q}) = \boldsymbol{\tau}
여기서
- \mathbf{M}(\mathbf{q}): 관성 행렬
- \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}): 코리올리/원심력 행렬
- \mathbf{g}(\mathbf{q}): 중력 벡터
- \boldsymbol{\tau}: 일반화 힘 벡터
이는 매니퓰레이터의 표준 동역학 방정식이다.
5. 운동 방정식의 성질
5.1 비선형성
다자유도 시스템의 운동 방정식은 일반적으로 강한 비선형이다. 관성 행렬과 코리올리/원심력 행렬이 자세에 의존하기 때문이다.
5.2 결합
각 자유도가 다른 자유도와 결합되어 있다. 한 관절의 운동이 다른 관절에 영향을 미친다.
5.3 비결합의 경우
특수한 경우에는 자유도 사이의 결합이 약하거나 없을 수 있다. 예를 들어 직선 관절만 가진 데카르트 매니퓰레이터의 경우 결합이 매우 약하다.
6. 운동 방정식의 풀이
6.1 순동역학
순동역학은 토크가 주어졌을 때 운동을 계산하는 문제이다.
\ddot{\mathbf{q}} = \mathbf{M}^{-1}(\mathbf{q})(\boldsymbol{\tau} - \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} - \mathbf{g}(\mathbf{q}))
이를 적분하여 위치와 속도를 계산한다.
6.2 역동역학
역동역학은 운동이 주어졌을 때 필요한 토크를 계산하는 문제이다.
\boldsymbol{\tau} = \mathbf{M}(\mathbf{q})\ddot{\mathbf{q}} + \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})\dot{\mathbf{q}} + \mathbf{g}(\mathbf{q})
이는 직접 계산이 가능하다.
6.3 수치 적분
순동역학의 수치 적분은 룽게-쿠타 등의 방법으로 수행된다.
7. 응용 예시: 2자유도 평면 매니퓰레이터
7.1 라그랑지언
L = \frac{1}{2}M_{11}\dot{q}_1^2 + M_{12}\dot{q}_1\dot{q}_2 + \frac{1}{2}M_{22}\dot{q}_2^2 - U(\mathbf{q})
7.2 운동 방정식
M_{11}\ddot{q}_1 + M_{12}\ddot{q}_2 + (\text{코리올리/원심력 항}) + g_1 = \tau_1
M_{12}\ddot{q}_1 + M_{22}\ddot{q}_2 + (\text{코리올리/원심력 항}) + g_2 = \tau_2
7.3 명시적 형식
이전 절에서 유도한 관성 행렬, 코리올리/원심력 행렬, 중력 벡터를 대입하면 명시적인 운동 방정식이 얻어진다.
8. 응용 예시: 두 입자의 결합 진동
8.1 모형
두 질량이 스프링으로 결합된 시스템에서 일반화 좌표는 두 변위 x_1, x_2이다.
8.2 라그랑지언
L = \frac{1}{2}m_1\dot{x}_1^2 + \frac{1}{2}m_2\dot{x}_2^2 - \frac{1}{2}k_1 x_1^2 - \frac{1}{2}k_{12}(x_2 - x_1)^2 - \frac{1}{2}k_2 x_2^2
8.3 운동 방정식
m_1\ddot{x}_1 + (k_1 + k_{12})x_1 - k_{12}x_2 = 0
m_2\ddot{x}_2 - k_{12}x_1 + (k_2 + k_{12})x_2 = 0
이는 결합 진동의 운동 방정식이다.
8.4 정상 모드
이 시스템의 정상 모드(normal mode)는 두 자유도가 같은 진동수로 진동하는 형태이다. 모드 분석을 통해 결정된다.
9. 본 절의 의의
본 절은 다자유도 시스템의 라그랑주 운동 방정식을 다루었다. 이는 매니퓰레이터를 비롯한 다양한 다자유도 시스템의 동역학 분석의 표준 방법이다.
10. 참고 문헌
- Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
- Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
- Greenwood, D. T. (2003). Advanced Dynamics. Cambridge University Press.
- Murray, R. M., Li, Z., & Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
version: 1.0