41.10 강성 시스템의 수치 해석 기법

강성 시스템(stiff system)은 빠른 동역학과 느린 동역학이 공존하는 시스템이다. 로봇의 접촉, 강한 스프링 등이 강성을 야기하며, 특수한 수치 해석 기법이 필요하다. 본 절에서는 강성 시스템의 수치 해석 기법을 학술적으로 다룬다.

1. 강성의 정의

1.1 큰 고유값 비

고유값의 크기 비가 매우 큰 시스템이다.

1.2 시간 스케일

여러 시간 스케일이 공존한다.

1.3 수치적 어려움

수치적으로 어려운 시스템이다.

2. 명시적 방법의 한계

2.1 작은 간격

매우 작은 시간 간격이 요구된다.

2.2 비효율

매우 비효율적이다.

2.3 실무 한계

실무적으로 한계가 있다.

3. 음함수 방법

3.1 큰 간격

큰 시간 간격이 가능하다.

3.2 안정성

안정성이 우수하다.

3.3 강성에 적합

강성 시스템에 적합하다.

4. BDF 방법

4.1 Backward Differentiation Formula

후방 차분 공식이다.

4.2 다단계 음함수

다단계 음함수 방법이다.

4.3 강성 표준

강성 문제의 표준이다.

5. Radau 방법

5.1 음함수 RK

음함수 Runge-Kutta이다.

5.2 고차

고차 정확도이다.

5.3 A-안정

A-안정이다.

6. 로봇 예시

6.1 접촉

지면 접촉이 강성을 야기한다.

6.2 강한 스프링

강한 스프링도 원인이다.

6.3 감속기

관절 감속기도 영향이다.

7. 반암시적

7.1 부분 암시적

선형 부분만 암시적이다.

7.2 효율적

효율적 구현이다.

7.3 실무적 타협

실무적 타협이다.

8. 적응적 방법

8.1 자동 전환

명시적·음함수 자동 전환이다.

8.2 MATLAB ode23s

MATLAB이 제공한다.

8.3 범용 적용

범용적으로 적용된다.

9. 특수 기법

9.1 지수 적분기

지수 적분기(exponential integrator)도 있다.

9.2 특수 구조 활용

시스템 구조 활용이다.

9.3 고급 방법

현대 고급 방법이다.

10. 학술적 활용

본 절에서 다룬 강성 시스템의 수치 해석 기법은 로봇 시뮬레이션의 학술적·실무적 도전이다. 적절한 방법 선택이 신뢰할 수 있는 시뮬레이션의 핵심이다.

11. 출처

Hairer, E. and Wanner, G., Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, 2nd edition, Springer, 1996.
Shampine, L. F. and Reichelt, M. W., “The MATLAB ODE Suite”, SIAM Journal on Scientific Computing, Vol. 18, No. 1, pp. 1–22, 1997.
Butcher, J. C., Numerical Methods for Ordinary Differential Equations, 3rd edition, Wiley, 2016.
Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P., Numerical Recipes: The Art of Scientific Computing, 3rd edition, Cambridge University Press, 2007.
Todorov, E., Erez, T., and Tassa, Y., “MuJoCo: A physics engine for model-based control”, Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 5026–5033, 2012.

12. 버전

문서 버전: 1.0
작성일: 2026-04-18