[논문 리뷰] A Dual Number Approach for Numerical Calculation of Velocity and Acceleration in the Spherical 4R Mechanism

F. Peñuñuri, R. Peón-Escalante외 3명|arXiv (Cornell University)|2013. 01. 08.

Robotic Mechanisms and Dynamics참고 문헌 13인용 수 5

한 줄 요약

이 논문은 한 번의 계산 단계로 벡터 함수의 일阶 및 이阶 도함수를 계산하기 위한 이중 수 기반의 수치적 방법을 제안한다. 이는 구면 4R 기구에서 속도 및 가속도를 정확하고 효율적으로 계산할 수 있도록 한다. 이 방법은 함수, 일阶 도함수, 이阶 도함수를 포함하는 세 구성 요소를 가진 확장된 이중 수를 사용하며, Fortran 모듈에 구현되어 있으며, 높은 정밀도로 커퍼 포인트 운동학을 계산하는 데 성공적으로 적용되었다.

ABSTRACT

This paper proposes a methodology to calculate both the first and second derivatives of a vector function of one variable in a single computation step. The method is based on the nested application of the dual number approach for first order derivatives. It has been implemented in Fortran language, a module which contains the dual version of elementary functions as well as more complex functions, which are common in the field of rotational kinematics. Since we have three quantities of interest, namely the function itself and its first and second derivative, our basic numerical entity has three elements. Then, for a given vector function $f:\mathbb{R} o \mathbb{R}^m$, its dual version will have the form $ ilde{f}:\mathbb{R}^3 o \mathbb{R}^{3m}$. As a study case, the proposed methodology is used to calculate the velocity and acceleration of a point moving on the coupler-point curve generated by a spherical four-bar mechanism.

연구 동기 및 목표

• 일계 및 이계 도함수의 근사 오차 없이 수치적 방법으로 계산하는 것을 목적으로 한다.
• 구면 4R 기구에서 분석적 도함수의 복잡성과 비효율성을 해결하고자 한다.
• 기구 설계 및 최적화 과정에서 속도 및 가속도를 효율적으로 통합하고자 한다.
• 실용적 사용을 위해 일반적인 운동학 함수의 이중 형태를 포함하는 Fortran 기반 모듈을 구현하고자 한다.
• 구면 4바디 기구에 대한 적용을 통해 이 방법의 정확성과 효율성을 입증하고자 한다.

제안 방법

• 이 방법은 $\tilde{f} = f + \epsilon f' + \epsilon^2 f''$ 형태의 확장된 이중 수를 사용하며, 여기서 $\epsilon^2 = 0$ 이고 $\epsilon^3 = 0$ 이다. 이를 통해 함수, 일계 도함수, 이계 도함수를 동시에 계산할 수 있다.
• 삼각함수(sin, cos, exp 등)와 벡터/행렬 연산은 연산자 오버로딩과 체인 룰 전파를 통해 이중 형태로 구현된다.
• 이 방법은 구면 4R 기구의 커퍼 포인트 위치 벡터에 적용되며, 이는 회전 행렬과 외적을 통해 표현된다.
• 커퍼 포인트의 궤적은 두 링크 벡터의 외적을 통해 정의된 축을 중심으로 벡터를 회전시켜 모델링된다.
• 확장된 이중 형태의 위치 벡터는 직접적으로 속도와 가속도를 그 일계 및 이계 도함수 성분으로 제공한다.
• 이 구현은 Fortran으로 작성되었으며, 일반적인 운동학 함수의 이중 형태를 제공하는 다운로드 가능한 모듈을 포함한다.

실험 결과

연구 질문

• RQ1 일계 및 이계 도함수를 동시에 한 번의 계산 단계로 계산할 수 있는가?
• RQ2 이중 수 접근법은 복잡한 운동학 함수를 포함한 구면 기구에 어떻게 적용될 수 있는가?
• RQ3 이 방법은 분석적 도함수 유도 없이도 커퍼 포인트 운동의 정확한 속도 및 가속도 값을 산출할 수 있는가?
• RQ4 유한차분 방법과 비교할 때 이중 기반 방법의 계산 효율성과 수치 정밀도는 어떠한가?
• RQ5 이 방법은 기구 설계 응용 분야에서 고성능 언어인 Fortran에 실제로 구현될 수 있는가?

주요 결과

• 이중 수 방법은 유한차분 방법에서 흔히 발생하는 절삭 오차나 상쇄 오차 없이 정확한 일계 및 이계 도함수 계산을 가능하게 한다.
• 여러 입력 각도 $\theta$ 에서의 속도 성분은 정확하게 계산되었으며, 결과는 표 2에 제공되어 있다.
• 입력 각도 $\theta = 0.00000$ 에서 속도 성분은 $\dot{x}_{\text{gen}} = -0.14255$, $\dot{y}_{\text{gen}} = -0.06884$, $\dot{z}_{\text{gen}} = 0.59467$ (SI 단위)였다.
• 입력 각도 $\theta = 3.14159$ 에서 가속도 성분은 $\ddot{x}_{\text{gen}} = -0.00803$, $\ddot{y}_{\text{gen}} = 0.19437$, $\ddot{z}_{\text{gen}} = -0.09747$였다.
• 확장된 이중 표현을 통해 도함수를 직접 추출할 수 있어 별도의 도함수 계산이 필요 없어진다.
• Fortran에 구현된 이 방법은 재사용 가능하고 효율적이며 정확한 구면 기구 운동학 분석 도구를 제공한다.