13.29 라그랑주 운동 방정식의 유도

1. 개요

라그랑주 운동 방정식(Lagrange’s equations of motion)은 라그랑지안에서 출발하여 시스템의 운동 방정식을 유도하는 표준 방법이다. 이는 1788년 조제프 루이 라그랑주(Joseph-Louis Lagrange)에 의해 Mécanique Analytique에서 제시되었으며, 해석 역학의 핵심이다. 본 절에서는 라그랑주 운동 방정식의 유도와 응용을 다룬다.

2. 라그랑주 방정식

2.1 정리

라그랑주 운동 방정식: 라그랑지안 $L = T - V$ 를 가진 시스템에 대해 다음의 방정식이 성립한다.

$\frac{d}{dt}\!\left(\frac{\partial L}{\partial\dot q_j}\right) - \frac{\partial L}{\partial q_j} = Q_j^{\text{non-cons}}, \quad j = 1, 2, \ldots, n$

여기서

$q_j$ : $j$ 번째 일반화 좌표
$\dot q_j$ : $j$ 번째 일반화 속도
$Q_j^{\text{non-cons}}$ : 비보존 일반화 힘 (보존력은 $V$ 에 포함되어 있음)

2.2 의미

이는 일반화 좌표와 같은 수의 운동 방정식을 제공한다. 자유도와 같은 수의 방정식이므로 시스템을 완전히 결정한다.

2.3 보존 시스템

비보존력이 없는 경우 ( $Q_j^{\text{non-cons}} = 0$ )

$\frac{d}{dt}\!\left(\frac{\partial L}{\partial\dot q_j}\right) - \frac{\partial L}{\partial q_j} = 0$

3. 유도의 출발점

3.1 달랑베르 원리

라그랑주 방정식의 유도는 달랑베르 원리에서 출발한다.

$\sum_i(\mathbf{F}_i^{\text{app}} - m_i\ddot{\mathbf{r}}_i) \cdot \delta\mathbf{r}_i = 0$

3.2 일반화 좌표로의 변환

가상 변위를 일반화 좌표로 표현한다.

$\delta\mathbf{r}_i = \sum_j\frac{\partial\mathbf{r}_i}{\partial q_j}\delta q_j$

3.3 관성력의 변환

관성력 항을 변환한다.

$\sum_i m_i\ddot{\mathbf{r}}_i \cdot \delta\mathbf{r}_i = \sum_j\!\left[\sum_i m_i\ddot{\mathbf{r}}_i \cdot \frac{\partial\mathbf{r}_i}{\partial q_j}\right]\delta q_j$

4. 운동 에너지의 처리

4.1 운동 에너지의 미분

운동 에너지를 일반화 좌표와 속도의 함수로 표현하면

$T = T(\mathbf{q}, \dot{\mathbf{q}}, t)$

다음의 항등식이 성립한다.

$\sum_i m_i\ddot{\mathbf{r}}_i \cdot \frac{\partial\mathbf{r}_i}{\partial q_j} = \frac{d}{dt}\!\left(\frac{\partial T}{\partial\dot q_j}\right) - \frac{\partial T}{\partial q_j}$

4.2 증명

이 항등식은 다음의 관계로부터 유도된다.

$\frac{\partial\mathbf{r}_i}{\partial q_j} = \frac{\partial\mathbf{v}_i}{\partial\dot q_j}$

(이는 운동학적 사실이다.)

이를 사용하여 식을 정리하면 위의 항등식이 얻어진다.

5. 라그랑주 방정식의 형태

5.1 기본 형태

위의 결과를 달랑베르 원리에 대입하면

$\sum_j\!\left[Q_j - \frac{d}{dt}\!\left(\frac{\partial T}{\partial\dot q_j}\right) + \frac{\partial T}{\partial q_j}\right]\delta q_j = 0$

여기서 $Q_j$ 는 모든 일반화 힘이다.

5.2 일반화 좌표의 독립성

일반화 좌표가 독립이면 각 $\delta q_j$ 의 계수가 0이어야 한다.

$\frac{d}{dt}\!\left(\frac{\partial T}{\partial\dot q_j}\right) - \frac{\partial T}{\partial q_j} = Q_j$

5.3 보존력의 분리

보존력은 위치 에너지로 표현되며 일반화 힘은 다음과 같다.

$Q_j^{\text{cons}} = -\frac{\partial V}{\partial q_j}$

따라서

$Q_j = Q_j^{\text{cons}} + Q_j^{\text{non-cons}} = -\frac{\partial V}{\partial q_j} + Q_j^{\text{non-cons}}$

5.4 라그랑지안의 도입

$L = T - V$ 를 사용하면 ( $V$ 가 $\dot q_j$ 에 의존하지 않는다고 가정)

$\frac{\partial L}{\partial\dot q_j} = \frac{\partial T}{\partial\dot q_j}$

$\frac{\partial L}{\partial q_j} = \frac{\partial T}{\partial q_j} - \frac{\partial V}{\partial q_j}$

5.5 결과

이를 대입하여 정리하면 라그랑주 운동 방정식이 얻어진다.

$\frac{d}{dt}\!\left(\frac{\partial L}{\partial\dot q_j}\right) - \frac{\partial L}{\partial q_j} = Q_j^{\text{non-cons}}$

6. 라그랑주 방정식의 사용

6.1 절차

시스템의 일반화 좌표를 선택한다.
운동 에너지 $T$ 를 일반화 좌표로 표현한다.
위치 에너지 $V$ 를 일반화 좌표로 표현한다.
라그랑지안 $L = T - V$ 를 구성한다.
각 일반화 좌표에 대해 라그랑주 방정식을 작성한다.
비보존 일반화 힘이 있으면 우변에 추가한다.

6.2 결과

각 일반화 좌표에 대한 미분 방정식이 얻어진다. 이것이 시스템의 운동 방정식이다.

7. 응용 예시: 단순 진자

7.1 라그랑지안

$L = \frac{1}{2}ML^2\dot\theta^2 + MgL\cos\theta$

7.2 편미분

$\frac{\partial L}{\partial\dot\theta} = ML^2\dot\theta$

$\frac{\partial L}{\partial\theta} = -MgL\sin\theta$

7.3 라그랑주 방정식

$\frac{d}{dt}(ML^2\dot\theta) - (-MgL\sin\theta) = 0$

$ML^2\ddot\theta + MgL\sin\theta = 0$

따라서 진자의 운동 방정식은

$\ddot\theta + \frac{g}{L}\sin\theta = 0$

8. 응용 예시: 매니퓰레이터

매니퓰레이터의 라그랑지안은 일반적으로

$L = \frac{1}{2}\dot{\mathbf{q}}^T\mathbf{M}(\mathbf{q})\dot{\mathbf{q}} - V(\mathbf{q})$

라그랑주 방정식을 적용하면 다음의 운동 방정식이 얻어진다.

$\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{C}$ : 코리올리/원심력 행렬
$\mathbf{g}$ : 중력 벡터
$\boldsymbol{\tau}$ : 관절 토크

이는 매니퓰레이터 동역학의 표준 형태이다.

9. 라그랑주 방법의 장점

9.1 좌표의 자유

라그랑주 방법은 일반화 좌표를 사용하므로 데카르트 좌표에 제한되지 않는다. 시스템에 적합한 좌표를 자유롭게 선택할 수 있다.

9.2 구속의 자동 만족

이상 구속의 경우 일반화 좌표가 구속을 자동으로 만족시키므로 구속력을 명시적으로 다루지 않아도 된다.

9.3 스칼라 양의 사용

라그랑주 방법은 스칼라 양(에너지)을 사용하므로 벡터 계산보다 단순할 수 있다.

9.4 대칭성과 보존 법칙

라그랑지안의 대칭성으로부터 보존 법칙이 직접 유도된다(뇌터 정리).

10. 본 절의 의의

본 절은 라그랑주 운동 방정식의 유도를 다루었다. 이는 해석 역학의 가장 강력한 도구 중 하나이며, 매니퓰레이터, 다체 시스템 등 복잡한 시스템의 동역학 분석에 표준이다.

11. 학습 권장사항

라그랑주 방정식의 형태를 정확히 이해한다.
유도 과정을 따라가 본다.
단순한 시스템에 적용해 본다.
매니퓰레이터의 운동 방정식 유도를 학습한다.
라그랑주 방법의 장점을 인식한다.

12. 참고 문헌

Lagrange, J. L. (1788). Mécanique Analytique.
Goldstein, H., Poole, C., & Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
Lanczos, C. (1986). The Variational Principles of Mechanics (4th ed.). Dover Publications.
Greenwood, D. T. (2003). Advanced Dynamics. Cambridge University Press.
Spong, M. W., Hutchinson, S., & Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.

version: 1.0