15.38 강성 시스템에 대한 수치 해석 고려사항
1. 개요
강성 시스템(stiff system)은 시간 척도의 차이가 큰 동역학 시스템이다. 강성 시스템의 수치 적분은 일반적인 적분 방법으로 어려우며 특별한 고려가 필요하다. 매니퓰레이터의 경우 유연 관절, 빠른 진동, 강한 감쇠 등이 강성을 유발한다. 본 절에서는 강성 시스템의 수치 해석 고려사항을 다룬다.
2. 강성 시스템의 정의
2.1 직관적 정의
강성 시스템은 시스템 안에 매우 빠른 동역학과 느린 동역학이 공존하는 시스템이다. 빠른 동역학의 특성 시간 척도가 느린 동역학에 비해 매우 짧다.
2.2 수학적 정의
선형 시스템 \dot{\mathbf{x}} = \mathbf{A}\mathbf{x}에서 \mathbf{A}의 고유값들 사이의 비율이 매우 크면 강성 시스템이다.
\frac{\max\vert\text{Re}(\lambda_i)\vert}{\min\vert\text{Re}(\lambda_i)\vert} \gg 1
2.3 비선형 시스템
비선형 시스템에서도 야코비안의 고유값으로 강성을 측정할 수 있다.
3. 강성의 원인
3.1 빠른 진동
매니퓰레이터의 일부에서 발생하는 매우 빠른 진동(예: 유연 관절의 진동)은 강성을 유발한다.
3.2 강한 감쇠
큰 감쇠 계수는 빠른 응답 시간을 만든다. 이는 강성 시스템이 될 수 있다.
3.3 큰 강성 계수
매니퓰레이터의 매우 큰 강성 계수(예: 거의 강체와 같은 행동)도 강성을 유발한다.
3.4 시간 척도의 분리
여러 시간 척도가 결합된 시스템(예: 빠른 모터와 느린 링크)은 강성이다.
4. 명시적 적분 방법의 한계
4.1 안정성 제약
명시적 적분 방법(오일러, RK4 등)은 강성 시스템에서 매우 작은 시간 단계가 필요하다. 안정성을 위해 시간 단계가 빠른 동역학의 시간 척도보다 작아야 한다.
h \leq \frac{C}{\vert\lambda_{\max}\vert}
여기서 C는 적분 방법에 의존하는 상수이다.
4.2 효율성 저하
시간 단계가 매우 작아지면 시뮬레이션이 매우 느려진다. 강성 시스템에서는 이는 실용적이지 않다.
4.3 발산
시간 단계가 너무 크면 수치 발산이 발생한다.
5. 암시적 적분 방법
5.1 안정성
암시적 적분 방법은 강성 시스템에 대해 무조건적 안정성을 가진다(또는 매우 큰 시간 단계에서도 안정).
5.2 후방 오일러
후방 오일러 방법은 가장 단순한 암시적 방법이다.
\mathbf{x}_{n+1} = \mathbf{x}_n + h\mathbf{f}(\mathbf{x}_{n+1}, t_{n+1})
이는 비선형 방정식의 풀이를 필요로 한다.
5.3 사다리꼴 방법
사다리꼴 방법은 2차 정확도이며 안정적이다.
\mathbf{x}_{n+1} = \mathbf{x}_n + \frac{h}{2}(\mathbf{f}_n + \mathbf{f}_{n+1})
5.4 BDF 방법
후방 미분 공식(BDF)은 강성 시스템에 효과적인 다단계 방법이다. BDF1, BDF2, …, BDF6 등이 있다.
5.5 Rosenbrock 방법
Rosenbrock 방법은 야코비안의 명시적 사용을 통해 효율적이다.
6. 비선형 방정식의 풀이
6.1 뉴턴 방법
암시적 방법에서 등장하는 비선형 방정식은 일반적으로 뉴턴 방법으로 풀린다.
\mathbf{x}^{k+1} = \mathbf{x}^k - [\mathbf{J}(\mathbf{x}^k)]^{-1}\mathbf{r}(\mathbf{x}^k)
여기서 \mathbf{r}은 잔차이고 \mathbf{J}는 야코비안이다.
6.2 야코비안의 계산
야코비안의 정확한 계산이 필요하다. 수치 미분 또는 자동 미분이 사용된다.
6.3 효율성
암시적 방법은 명시적 방법보다 단계당 비용이 높지만, 더 큰 시간 단계가 가능하므로 강성 시스템에서 더 효율적이다.
7. 매니퓰레이터에서의 응용
7.1 유연 관절 매니퓰레이터
유연 관절 매니퓰레이터는 일반적으로 강성 시스템이다. 모터의 빠른 동역학과 링크의 느린 동역학이 결합되어 있다.
7.2 시간 척도 분리
강성 시스템은 시간 척도 분리를 통해 처리할 수 있다. 빠른 동역학을 분리하여 별도로 처리하고 느린 동역학에 평균 효과로 통합한다.
7.3 특이 섭동
특이 섭동 이론(singular perturbation theory)은 강성 시스템의 분석에 사용된다.
8. 적응 시간 단계
8.1 동기
적응 시간 단계는 시스템의 변화율에 따라 시간 단계를 조정한다. 빠른 변화에서는 작은 단계, 느린 변화에서는 큰 단계를 사용한다.
8.2 오차 추정
오차 추정에 기반하여 시간 단계가 조정된다.
8.3 효율성
적응 시간 단계는 강성 시스템에서 효율성을 크게 향상시킬 수 있다.
9. 수치 발산의 처리
9.1 검출
수치 발산은 시뮬레이션의 결과가 비현실적인 값으로 발산하는 현상이다. 검출이 필요하다.
9.2 시간 단계의 감소
발산이 검출되면 시간 단계를 줄이거나 더 안정한 적분 방법을 사용한다.
9.3 모형의 재검토
발산이 지속되면 모형 자체에 문제가 있을 수 있다. 매개 변수와 식을 다시 확인한다.
10. 응용
10.1 유연 관절 매니퓰레이터 시뮬레이션
유연 관절 매니퓰레이터의 시뮬레이션은 강성 시스템을 다루는 적분기를 필요로 한다.
10.2 접촉 시뮬레이션
접촉 시뮬레이션에서는 매우 큰 강성이 등장한다. 강성 시스템 처리가 필요하다.
10.3 빠른 진동의 분석
빠른 진동을 정확히 시뮬레이션하기 위해 강성 시스템 처리가 필요하다.
11. 본 절의 의의
본 절은 강성 시스템의 수치 해석 고려사항을 다루었다. 강성 시스템의 정확한 시뮬레이션은 매니퓰레이터의 다양한 응용에서 중요하며, 적절한 적분 방법의 선택이 필수적이다.
12. 참고 문헌
- Hairer, E., Lubich, C., & Wanner, G. (1996). Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems (2nd ed.). Springer.
- Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. (2007). Numerical Recipes (3rd ed.). Cambridge University Press.
- Khalil, H. K. (2002). Nonlinear Systems (3rd ed.). Prentice Hall.
- Featherstone, R. (2008). Rigid Body Dynamics Algorithms. Springer.
version: 1.0