17.15 코리올리력 항과 원심력 항의 물리적 해석

1. 개요

매니퓰레이터의 동역학 운동 방정식에서 관성 가속 항 다음으로 등장하는 \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} 항은 속도 의존성 관성 효과를 포함한다. 이 항은 두 가지 서로 다른 물리적 기원을 가진 힘을 포함하며, 각각은 고전적으로 원심력(centrifugal force)과 코리올리력(Coriolis force)이라 불린다. 원심력은 단일 관절의 회전에 기인하는 속도 제곱 비례 기여이며, 코리올리력은 두 관절이 동시에 회전할 때 발생하는 교차 기여이다.

본 절에서는 이 두 유형의 힘이 어떻게 관성 행렬의 자세 의존성에서 비롯되는지, 각각의 물리적 의미는 무엇인지, 그리고 실제 매니퓰레이터에서 어떻게 관찰되고 어떤 영향을 미치는지를 학술적으로 기술한다. 코리올리-원심력 행렬의 반대칭 성질과 같은 대수적 특성은 별도의 절에서 다루어지며, 본 절은 물리적 해석과 실무적 함의에 초점을 맞춘다.

2. 속도 의존 관성력의 기원

동역학 운동 방정식의 유도 과정에서 관성 행렬이 자세 의존적이라는 사실은 불가피하게 속도 의존 힘을 발생시킨다. 구체적으로, 라그랑지안 형식의 운동 방정식

\frac{d}{dt} \left( \frac{\partial K}{\partial \dot{q}_i} \right) - \frac{\partial K}{\partial q_i} = \tau_i

에서 운동 에너지를 관성 행렬로 표현한 후 시간 미분과 편미분을 수행하면 다음의 형태가 얻어진다.

\sum_j M_{ij} \ddot{q}_j + \sum_{j,k} \left( \frac{\partial M_{ij}}{\partial q_k} - \frac{1}{2} \frac{\partial M_{jk}}{\partial q_i} \right) \dot{q}_k \dot{q}_j + \frac{\partial U}{\partial q_i} = \tau_i

여기서 둘째 항이 속도 의존 관성력이며, 이는 관성 행렬의 자세 편미분과 관절 속도의 곱에서 발생한다. 관성 행렬이 상수라면(즉 자세에 의존하지 않는다면) 이 항은 영이 되며, 속도 의존 관성력은 나타나지 않는다. 따라서 코리올리력과 원심력은 관성 행렬의 자세 의존성의 직접적 산물이다.

3. 크리스토펠 기호를 통한 표현

코리올리-원심력 항은 크리스토펠 기호(Christoffel symbol)를 이용하여 다음과 같이 표현된다.

c_{ijk}(\mathbf{q}) = \frac{1}{2} \left( \frac{\partial M_{ij}}{\partial q_k} + \frac{\partial M_{ik}}{\partial q_j} - \frac{\partial M_{jk}}{\partial q_i} \right)

그리고 이를 이용하여 코리올리 행렬이 다음과 같이 구성된다.

C_{ij}(\mathbf{q}, \dot{\mathbf{q}}) = \sum_{k=1}^{n} c_{ijk}(\mathbf{q}) \dot{q}_k

이 표현에서 크리스토펠 기호는 매니퓰레이터 관성 행렬의 1계 편미분으로 결정되며, 자세의 함수이다. 속도 의존 힘의 i번째 성분은 다음과 같이 이중합의 형태로 나타난다.

\sum_{j,k} c_{ijk}(\mathbf{q}) \dot{q}_j \dot{q}_k

이 이중합에서 j = k인 대각 항은 단일 관절 속도의 제곱에 비례하며 원심력에 해당하고, j \neq k인 비대각 항은 두 관절 속도의 곱에 비례하며 코리올리력에 해당한다.

원심력 항의 물리적 의미

원심력 항 c_{ijj}(\mathbf{q}) \dot{q}_j^2는 관절 j만 회전할 때 관절 i에 유발되는 속도 제곱 비례 힘을 나타낸다. 고전 역학에서 원심력은 회전 기준계에서 나타나는 관성력으로, 회전축으로부터 멀어지려는 방향의 힘으로 해석된다. 매니퓰레이터에서는 한 관절이 회전할 때 그 하류에 있는 링크가 회전축으로부터 멀어지는 경향을 갖게 되며, 이 경향이 해당 관절 또는 다른 관절에 반대 방향의 반응력을 유발한다.

구체적 예로, 어깨 관절만 회전하고 팔꿈치는 정지해 있는 자세를 생각하자. 팔꿈치가 정지 상태를 유지하려면 팔꿈치 관절이 하류 링크의 원심력에 대응하는 토크를 공급해야 한다. 이러한 토크는 어깨 관절 속도의 제곱에 비례하며, 관성 행렬의 편미분으로부터 유도된다. 어깨 속도가 2배로 증가하면 필요한 팔꿈치 토크는 4배로 증가하는 이차 관계가 성립한다.

원심력은 항상 단일 관절 속도의 제곱이므로 부호를 가지지 않으며, 속도의 방향과 무관하게 동일한 크기로 작용한다. 원심력의 방향과 크기는 관성 행렬의 편미분의 부호와 크기에 의하여 결정되며, 매니퓰레이터의 자세에 따라 크게 변한다.

코리올리력 항의 물리적 의미

코리올리력 항은 두 관절이 동시에 회전할 때 나타나는 결합 효과이다. 고전 역학에서 코리올리력은 회전 기준계 내부에서 움직이는 질점이 받는 관성력으로, 질점의 속도와 기준계의 각속도의 벡터 외적에 비례한다. 매니퓰레이터 동역학에서는 한 관절의 회전이 다른 관절의 링크 운동에 의하여 영향을 받는 경우 코리올리력이 발생한다.

구체적으로, 어깨와 팔꿈치가 동시에 회전하는 경우 팔꿈치 하류에 있는 링크는 어깨 회전에 의하여 가속되는 동시에 팔꿈치 회전에 의하여도 가속된다. 두 회전이 결합되면 단순한 중첩 이상의 추가 관성 효과가 발생하며, 이는 두 관절 속도의 곱에 비례하는 코리올리 토크로 나타난다.

코리올리력은 두 속도의 곱이므로 한 속도의 부호가 바뀌면 부호가 바뀐다. 따라서 원심력과 달리 방향성 있는 힘이며, 관절 속도 방향의 상대적 관계에 따라 부호가 결정된다.

3자유도 매니퓰레이터에서의 예시

3자유도 이상의 매니퓰레이터에서는 원심력과 코리올리력의 구조가 더 풍부하다. 세 관절 q_1, q_2, q_3이 있을 때, 관절 i에 작용하는 속도 의존 관성력은 다음의 항들을 포함한다.

c_{i11} \dot{q}_1^2 + c_{i22} \dot{q}_2^2 + c_{i33} \dot{q}_3^2 + 2 c_{i12} \dot{q}_1 \dot{q}_2 + 2 c_{i13} \dot{q}_1 \dot{q}_3 + 2 c_{i23} \dot{q}_2 \dot{q}_3

처음 세 항은 각 단일 관절 회전의 원심력 기여이고, 마지막 세 항은 쌍별 코리올리 결합이다. 각 기여의 크기는 자세에 따라 결정되며, 일부 자세에서는 특정 항이 지배적일 수 있다.

특히 회전축이 서로 평행하지 않은 관절들 사이에서 코리올리 결합이 크며, 회전축이 직교인 경우 극단적 기여가 나타날 수 있다. 반대로 회전축이 평행하고 서로 독립적으로 배치된 관절들 사이에서는 결합이 약하다.

4. 고속 운동에서의 지배력

속도 의존 관성력은 속도의 이차에 비례하므로, 고속 운동에서는 중력이나 외부 상호작용 힘보다 지배적이 될 수 있다. 대략 관절 속도가 \dot{q}_0 배로 증가하면 관성 토크는 관절 가속도에 비례하여 증가하는 반면, 속도 의존 관성력은 \dot{q}_0^2에 비례하여 증가한다. 따라서 빠른 가감속을 반복하는 작업에서 속도 의존 힘의 영향은 중요한 고려 대상이 된다.

전형적 산업용 매니퓰레이터의 최대 속도 영역에서 코리올리력과 원심력은 중력 토크와 유사한 수준이 될 수 있으며, 일부 자세에서는 중력을 초과하기도 한다. 이러한 조건에서는 속도 의존 힘을 보상하지 않는 제어기는 궤적 추종 오차가 커지며, 동역학 기반 보상이 필수적이 된다.

5. 저속과 준정적 조건에서의 근사적 무시

반대로 저속 운동이나 준정적 조건에서는 속도 의존 항이 다른 항에 비하여 충분히 작아 무시될 수 있다. 관절 속도가 매우 작을 때 코리올리력과 원심력의 크기는 속도 제곱에 비례하여 빠르게 감소하며, 중력 토크가 지배적이 된다. 이러한 조건에서 준정적 해석이 정당화되며, 그 구체적 근거는 별도의 절에서 다루어졌다.

6. 속도 의존 힘과 에너지 전달

속도 의존 관성력은 에너지를 시스템 내부에서 관절 사이에 재분배하는 역할을 한다. 구체적으로 다음의 성질이 성립한다. 코리올리-원심력 항이 관절 속도 벡터와 이루는 내적은 에너지 변화에 영향을 미치지 않으며, 이는 다음의 관계로 표현된다.

\dot{\mathbf{q}}^T \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} = \frac{1}{2} \dot{\mathbf{q}}^T \dot{\mathbf{M}}(\mathbf{q}) \dot{\mathbf{q}}

이 등식은 운동 에너지의 시간 미분이 \dot{K} = \dot{\mathbf{q}}^T \mathbf{M} \ddot{\mathbf{q}} + \frac{1}{2} \dot{\mathbf{q}}^T \dot{\mathbf{M}} \dot{\mathbf{q}}임을 이용하여 유도된다. 그 결과로 코리올리-원심력 항은 운동 에너지를 생성하거나 소산시키지 않고 단지 관절 사이에 재분배할 뿐이다. 이 성질은 동역학 방정식의 수동성(passivity)과 직결되며, 안정성 증명에서 중요하게 활용된다.

해석적 분해와 수치적 특성

코리올리 행렬은 수치적으로 명시적 편미분 계산을 필요로 한다. 관성 행렬의 각 원소에 대하여 모든 관절에 대한 편미분을 계산한 후, 크리스토펠 기호로 조립하여 행렬 원소를 얻는다. 이러한 직접 계산은 O(n^3) 또는 더 높은 복잡도를 가지며, 실시간 제어에서는 효율성이 중요하다.

효율적 대안으로 재귀적 뉴턴-오일러 알고리즘을 이용하는 방법이 있다. 이 방법은 \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}} 벡터를 직접 계산하며, 관성 행렬과 크리스토펠 기호를 명시적으로 구성하지 않고도 같은 결과를 얻는다. 이 재귀 계산은 O(n)의 복잡도를 가지며, 실시간 구현에 적합하다.

매니퓰레이터 설계와 관성 결합

코리올리력과 원심력은 매니퓰레이터의 관성 결합 구조에 의하여 결정되므로, 기계 설계 단계에서 이러한 결합을 최소화하는 설계 전략이 제안된다. 예를 들어 링크의 질량 중심을 회전축에 가깝게 배치하거나, 직렬형 대신 병렬형 구조를 채택하여 관성 결합을 분산하는 방법 등이 있다. 그러나 완전한 분리는 일반적으로 불가능하며, 어느 정도의 결합은 불가피하다.

설계 단계에서의 관성 결합 최소화는 이후 제어 설계의 부담을 줄이며, 특히 모델 기반 보상의 정확성에 대한 의존도를 낮춘다. 이는 특히 산업용 고속 매니퓰레이터와 의료 정밀 로봇에서 중요하게 고려된다.

제어기 설계에의 함의

속도 의존 관성력은 제어기 설계에서 다음의 접근으로 다루어진다. 첫째, 전향 보상은 코리올리-원심력 토크를 모델 기반으로 계산하여 지령 토크에 직접 더하는 방법이다. 이는 계산 토크 제어와 피드백 선형화의 일부를 이룬다. 둘째, 적응 제어는 관성 매개변수의 선형성을 이용하여 코리올리 토크의 기여를 동정한다. 셋째, 강건 제어는 코리올리 토크의 불확실성에 대하여 여유를 두어 설계 이득을 조정한다. 넷째, 수동성 기반 제어는 코리올리-원심력 항의 에너지 보존 성질을 이용하여 추가 보상 없이 안정성을 증명한다.

각 접근은 장단점을 가지며, 응용에 따라 선택된다. 모델의 정확도가 충분하면 전향 보상이 가장 효과적이지만, 불확실성이 큰 상황에서는 적응 또는 강건 기법이 선호된다.

본 절의 의의

본 절에서 다룬 코리올리력 항과 원심력 항의 물리적 해석은 매니퓰레이터 동역학의 속도 의존 비선형성을 이해하는 데 필수적이다. 이 항들은 관성 행렬의 자세 의존성에서 자연스럽게 유도되며, 단순한 수학적 산물이 아니라 명확한 물리적 해석을 가지는 힘이다. 고속 운동에서 이 힘들은 제어 성능의 결정적 요소가 되며, 저속 운동에서도 무시할 수 없는 경우가 있다.

본 절의 내용은 다음 절에서 다루어질 코리올리 행렬의 반대칭 성질과 안정성 증명의 기반이 된다. 또한 동역학 기반 제어, 매개변수 식별, 그리고 수치 시뮬레이션의 여러 단계에서 반복적으로 참조된다.

학습 권장사항

본 절의 내용을 충분히 이해하기 위하여 고전 역학의 회전 기준계, 라그랑주 방정식, 그리고 크리스토펠 기호의 정의에 대한 선행 학습이 권장된다. 본 절의 내용을 심화 학습하기 위해서는 리만 기하학, 다체 동역학의 재귀 알고리즘, 그리고 수동성 기반 제어 이론에 대한 추가 학습이 권장된다.

참고 문헌

• Spong, M. W., Hutchinson, S., and Vidyasagar, M. (2020). Robot Modeling and Control (2nd ed.). Wiley.
• Siciliano, B., Sciavicco, L., Villani, L., and Oriolo, G. (2009). Robotics: Modelling, Planning and Control. Springer.
• Craig, J. J. (2018). Introduction to Robotics: Mechanics and Control (4th ed.). Pearson.
• Murray, R. M., Li, Z., and Sastry, S. S. (1994). A Mathematical Introduction to Robotic Manipulation. CRC Press.
• Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
• Goldstein, H., Poole, C., and Safko, J. (2002). Classical Mechanics (3rd ed.). Addison-Wesley.
• Khalil, W., and Dombre, E. (2002). Modeling, Identification and Control of Robots. Hermes Penton.

version: 1.0