17.9 준정적 해석의 원리와 가정

1. 개요

엄밀한 정적 해석은 시스템의 모든 질점 또는 강체의 속도와 가속도가 영이라는 전제 아래에서 수행된다. 그러나 실제 로봇 작업의 많은 경우 매니퓰레이터는 정지 상태가 아니라 매우 느리게 움직이며, 이때 속도와 가속도가 존재하지만 그 크기가 충분히 작아서 관성 효과가 상대적으로 무시할 수 있는 수준에 머무른다. 준정적 해석(quasi-static analysis)은 이러한 조건에서 관성 항을 생략하여 정적 평형 방정식과 실질적으로 동일한 형식의 해석을 수행하는 기법이며, 로봇 정역학과 동역학의 경계에 위치한다.

본 절에서는 준정적 해석이 적용 가능한 조건의 물리적 의미, 근사의 수학적 근거, 오차 한계, 그리고 준정적 해석이 상정하는 기본 가정을 학술적으로 기술한다. 준정적 해석의 구체적 운동 계획 응용은 별도의 절에서 다루며, 본 절은 해석 원리 자체에 초점을 맞춘다.

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{J}^T(\mathbf{q}) \mathbf{F}_e

여기서 \mathbf{M}(\mathbf{q})는 관성 행렬, \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}})는 코리올리와 원심력을 표현하는 행렬, \mathbf{g}(\mathbf{q})는 중력 토크, \boldsymbol{\tau}는 관절 토크, \mathbf{F}_e는 외부 렌치이다. 준정적 근사는 이 방정식에서 관성 항 \mathbf{M} \ddot{\mathbf{q}}와 속도 의존 항 \mathbf{C} \dot{\mathbf{q}}를 모두 영으로 간주하여 다음의 단순화된 방정식을 얻는 것에 해당한다.

\mathbf{g}(\mathbf{q}) = \boldsymbol{\tau} + \mathbf{J}^T(\mathbf{q}) \mathbf{F}_e

이 방정식은 외부 렌치와 관절 토크가 중력 토크를 균형 잡아야 한다는 조건으로 정적 평형 방정식과 동일한 형태이다. 차이점은 \mathbf{q}가 시간에 따라 천천히 변하는 변수이며, 방정식이 매 순간 성립한다고 가정한다는 점이다.

이 근사의 물리적 의미는 관성에 저장되는 운동 에너지가 위치 에너지와 외부 작업에 비하여 무시할 수 있을 정도로 작다는 것이다. 즉 시스템이 한 평형 자세에서 다른 평형 자세로 천천히 이동하는 동안, 각 순간의 상태는 사실상 평형 상태에 해당하는 것으로 간주된다.

3. 준정적 가정의 수학적 정식화

준정적 근사가 정당화되려면 관성 항과 속도 의존 항이 다른 항들에 비하여 충분히 작아야 한다. 이를 정량적으로 평가하기 위하여 대표적 시간 척도와 속도 척도를 도입한다. 매니퓰레이터가 특성 시간 T 동안 특성 변위 \Delta q를 이동한다고 할 때, 특성 속도는 \dot{q} \sim \Delta q / T이고 특성 가속도는 \ddot{q} \sim \Delta q / T^2이다.

관성 토크의 크기는 \mathbf{M}(\mathbf{q}) \ddot{\mathbf{q}}의 크기에 의하여 결정되며, 대표 스케일에서는 \mu_M \Delta q / T^2 수준이다. 여기서 \mu_M은 관성 행렬의 대표적 원소 크기이다. 반면 중력 토크의 대표 크기 \mu_g는 시간 척도에 무관하게 자세와 질량에 의하여 결정된다. 따라서 관성 기여와 중력 기여의 비율은 다음과 같이 나타난다.

\eta = \frac{\mu_M \Delta q / T^2}{\mu_g}

준정적 근사는 이 비율이 충분히 작을 때, 즉 \eta \ll 1일 때 정당화된다. 이는 운동이 특성 시간 T 안에서 \Delta q만큼 변화할 때 중력에 의한 토크가 관성에 의한 토크보다 훨씬 크다는 조건과 동등하다.

유사하게, 코리올리 및 원심력 항 \mathbf{C}(\mathbf{q}, \dot{\mathbf{q}}) \dot{\mathbf{q}}의 크기는 속도의 제곱에 비례하므로 \mu_C (\Delta q / T)^2의 수준이며, 이 또한 중력 토크보다 작아야 한다. 이 두 조건이 동시에 만족되면 준정적 근사는 충분한 정확성으로 성립한다.

시간 척도 분리와 특이 섭동

준정적 가정은 특이 섭동 이론(singular perturbation theory)의 관점에서 해석될 수 있다. 매니퓰레이터 동역학에서 작은 파라미터 \epsilon = T_{\text{slow}} / T_{\text{fast}}를 도입하여 느린 시간 척도와 빠른 시간 척도를 분리할 수 있다. 관성이 지배하는 빠른 진동은 특성 시간 T_{\text{fast}}에 속하며, 중력과 외력이 지배하는 느린 평형 이동은 T_{\text{slow}}에 속한다. 운동이 느린 시간 척도에서 분석될 때, 빠른 시간 척도의 동역학은 즉각적으로 평형에 도달하는 것으로 간주되며, 이는 준정적 근사의 수학적 기반을 형성한다.

이러한 시간 척도 분리는 준정적 해석이 원 동역학 해의 외부 해(outer solution)에 해당함을 의미한다. 초기 조건에 의하여 관성적 진동이 발생하는 짧은 경계층(boundary layer)을 제외하면, 느린 시간 척도에서의 움직임은 준정적 방정식의 해로 근사된다. 경계층 내부의 진동은 별도의 내부 해(inner solution)로 분석되며, 일반적 경우 마찰과 감쇠에 의하여 빠르게 감쇠한다.

준정적 해석의 표준 가정

준정적 해석을 엄밀하게 수행하기 위하여 다음의 기본 가정이 전제된다.

첫째, 모든 관절 가속도 \ddot{\mathbf{q}}가 충분히 작아서 관성 토크가 무시 가능하다. 수치적으로는 관성 항의 크기가 중력 항의 크기의 일정 비율(예: 5%) 아래에 머무른다는 조건이 설정된다. 둘째, 관절 속도 \dot{\mathbf{q}}가 충분히 작아서 코리올리 및 원심력 항 또한 무시 가능하다. 셋째, 외부 환경과의 상호작용력이 임펄스적 충격을 포함하지 않으며, 매끄럽게 변화한다. 넷째, 액추에이터의 대역폭이 준정적 변화를 추종하기에 충분하여, 토크 지령이 즉각적으로 실제 토크로 구현된다. 다섯째, 관절 유연성과 링크 탄성에 기인한 고주파 진동 모드는 감쇠되어 있거나 자극되지 않는다.

이러한 가정이 성립하는 작업 상황은 정밀 조립, 부드러운 교시, 저속 수동 안내, 느린 궤적의 그리기 및 용접 작업, 그리고 무중력 모사 실험 등이 포함된다. 반대로 급격한 가속, 충격 접촉, 고속 픽앤플레이스, 그리고 진동을 수반하는 작업은 준정적 해석의 적용 범위를 벗어난다.

관성과 중력 토크 비율의 정량 평가

주어진 매니퓰레이터와 궤적에 대하여 준정적 근사의 타당성을 정량적으로 평가하려면, 실제 궤적에서의 관성 토크와 중력 토크의 크기를 계산하여 비교한다. 한 가지 표준 지표는 다음과 같이 정의된다.

\eta_i(t) = \frac{|[\mathbf{M}(\mathbf{q}(t)) \ddot{\mathbf{q}}(t)]_i|}{|[\mathbf{g}(\mathbf{q}(t))]_i| + |[\mathbf{J}^T(\mathbf{q}(t)) \mathbf{F}_e(t)]_i| + \epsilon_0}

여기서 \epsilon_0는 영으로 나누는 것을 방지하기 위한 작은 양의 상수이다. 이 지표가 각 관절 i에서 사전 정의된 임계값(예: 0.05) 아래에 머무르면, 해당 관절에 대하여 준정적 근사가 유효하다고 판단한다. 유사한 지표가 코리올리 및 원심력 항에 대해서도 정의된다.

실무에서는 계획 궤적에 대하여 이러한 지표를 사전 시뮬레이션을 통해 평가하고, 지표가 큰 구간이 있으면 해당 구간에서는 완전 동역학 해석으로 전환하거나 궤적의 속도 프로파일을 조정하여 준정적 조건을 유지한다.

4. 준정적 해석의 해로서의 연속 자세

준정적 해석에서 자세 \mathbf{q}(t)는 시간의 매끄러운 함수로 주어지며, 매 순간 관절 토크와 외력이 중력과 균형을 이루는 조건 아래에서 결정된다. 시간에 따른 자세의 변화는 한 평형 자세에서 인접한 평형 자세로의 연속적 이동이며, 이러한 이동의 경로는 평형 자세들의 집합 안에서 정의된다.

외부 렌치 \mathbf{F}_e(t)가 시간의 함수로 주어지는 경우, 관절 토크는 다음과 같이 시간의 함수로 표현된다.

\boldsymbol{\tau}(t) = \mathbf{g}(\mathbf{q}(t)) - \mathbf{J}^T(\mathbf{q}(t)) \mathbf{F}_e(t)

이 관계는 준정적 궤적을 따라 각 순간의 관절 토크를 결정하는 시간 의존적 정적 방정식의 연속이며, 원 동역학 방정식의 느린 시간 근사에 해당한다.

준정적 해석이 예측하지 못하는 현상

준정적 해석은 단순성과 계산 효율에서 큰 이점을 가지는 반면, 관성과 운동 에너지가 본질적 역할을 하는 다음의 현상들은 예측할 수 없다.

첫째, 공진과 진동. 관성과 강성이 결합하여 고유 진동 모드를 형성하며, 준정적 해석에는 이러한 모드가 나타나지 않는다. 둘째, 속도에 의존하는 궤적 오차. 실제 로봇은 속도가 빨라질수록 궤적 추종 오차가 증가하지만, 준정적 해석은 이러한 속도 의존성을 포착하지 못한다. 셋째, 가속 기인 힘. 급격한 가속 단계에서는 관성력이 상호작용력의 지배 요인이 될 수 있으며, 준정적 모형은 이를 놓친다. 넷째, 에너지 저장과 회수. 빠른 운동에서 관성에 저장된 에너지가 이후 단계에서 회수되는 현상이 있으며, 이는 동역학적 분석에서만 나타난다.

이러한 한계는 준정적 해석이 상위 계획, 초기 설계, 그리고 저속 작업에는 적합하지만 정밀한 고속 동적 제어에는 충분하지 않음을 의미한다. 완전 동역학 해석과의 상호 보완이 필요하다.

마찰과 정상 상태에서의 고려

준정적 해석에서는 마찰 모델의 처리가 특별한 주의를 요한다. 정적 마찰(static friction)은 외력과 관절 토크가 작은 임계값 이하일 때 관절의 운동을 방지하여 매니퓰레이터가 한 자세에 머무르게 한다. 이 경우 관절 토크는 단일 값이 아니라 마찰 한계 내에서의 집합으로 결정된다. 반면 운동 중에는 쿨롱 마찰(Coulomb friction)과 점성 마찰(viscous friction)이 속도의 방향과 크기에 따라 토크에 기여한다.

준정적 해석의 맥락에서 속도가 매우 작으면 점성 마찰 기여는 무시 가능하지만, 쿨롱 마찰은 운동 방향에만 의존하고 크기는 일정한 값이므로 관절 토크 방정식에 비무시적으로 기여한다. 이러한 마찰 성분은 평형 방정식의 우변에 추가되어 다음과 같이 나타난다.

\boldsymbol{\tau} = \mathbf{g}(\mathbf{q}) - \mathbf{J}^T(\mathbf{q}) \mathbf{F}_e + \mathbf{f}_{\text{fric}}(\dot{\mathbf{q}})

마찰 항은 준정적 해석의 유일성에 영향을 미치며, 마찰 자체의 정밀한 모델링은 별도의 절에서 다루어진다.

5. 준정적 가정이 실패하는 경계 사례

준정적 해석이 실패하거나 오차가 큰 경우는 다음과 같다. 첫째, 접촉의 돌발적 발생 또는 해제. 접촉 상태가 갑자기 변할 때 접촉력이 순간적으로 변화하여 유한한 관성 반응을 유발한다. 둘째, 강직한 부하에 의한 충격. 유연성이 낮은 부하와 충돌할 때 고주파 진동과 큰 관성력이 발생한다. 셋째, 궤적의 급격한 방향 변화. 모서리 또는 코너를 포함하는 궤적에서는 순간적으로 큰 가속도가 요구된다. 넷째, 유연 관절과 링크의 저주파 모드가 자극되는 경우. 외력이 고유 진동수에 가까운 주파수로 가해지면 진동 응답이 증폭된다.

이러한 경우에는 준정적 해석이 부적합하며, 완전한 동역학 해석 또는 혼합 접근(경계층 수정)이 필요하다. 본 장의 동역학 관련 절들은 이러한 문제에 대한 정식화를 다룬다.

6. 본 절의 의의

본 절에서 다룬 준정적 해석의 원리와 가정은 로봇공학의 많은 분석과 설계 작업에서 암묵적으로 또는 명시적으로 전제된다. 정밀 조립, 힘 제어 실험, 협동 안내, 그리고 저속 작업에서의 경로 계획은 준정적 가정 아래에서 수행되며, 이 가정이 정확히 이해되지 않으면 설계된 제어기가 예상과 다르게 작동할 수 있다.

또한 준정적 해석은 완전 동역학 해석으로 이행하기 전의 첫 단계로 작용한다. 초기 설계 단계에서는 준정적 모형이 관절 토크 요구 사항과 구조 하중을 예비적으로 평가하는 데 사용되며, 결과가 만족스러울 때 비로소 동역학적 세부 해석이 수반된다. 본 절의 원리는 이러한 단계적 접근의 근거를 제공한다.

7. 학습 권장사항

본 절의 내용을 충분히 이해하기 위하여 상미분 방정식의 기본, 시간 척도 분석, 그리고 강체의 동역학 운동 방정식에 대한 선행 학습이 권장된다. 본 절의 내용을 심화 학습하기 위해서는 특이 섭동 이론, 느린-빠른 시스템의 분해, 그리고 마찰 모델링과 관련된 비매끄러운 동역학에 대한 추가 학습이 권장된다.

8. 참고 문헌

• Khalil, H. K. (2015). Nonlinear Systems (3rd ed.). Pearson.
• Kokotović, P. V., Khalil, H. K., and O’Reilly, J. (1986). Singular Perturbation Methods in Control: Analysis and Design. Academic Press.
• 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.
• 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.

version: 1.0