[论文解读] 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 more|arXiv (Cornell University)|Jan 8, 2013

Robotic Mechanisms and Dynamics参考文献 13被引用 5

一句话总结

本文提出一种基于扩展双数的数值方法，可在单步计算中同时求解向量函数的一阶与二阶导数，从而实现球面四杆机构中速度与加速度的精确高效计算。该方法采用包含三个分量（函数值、一阶导数、二阶导数）的扩展双数，在 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.

研究动机与目标

• 开发一种无需近似误差的数值方法，用于计算向量函数的一阶与二阶导数。
• 克服球面四杆机构中解析导数方法的复杂性与低效性。
• 实现在机构综合与优化中高效引入速度与加速度的计算。
• 在 Fortran 中实现一个模块，提供基本函数与运动学函数的双数版本，以供实际应用。
• 通过在球面四杆机构上的应用，展示该方法的精度与效率。

提出的方法

• 该方法使用形式为 $\tilde{f} = f + \epsilon f' + \epsilon^2 f''$ 的扩展双数，其中 $\epsilon^2 = 0$ 且 $\epsilon^3 = 0$，从而实现函数值、一阶导数与二阶导数的同时计算。
• 通过操作符重载与链式法则传播，将基本函数（如 sin、cos、exp）以及向量/矩阵运算转换为双数形式。
• 该方法应用于球面四杆机构中连杆点的位置矢量，该矢量通过旋转矩阵与叉积表示。
• 连杆点轨迹通过绕由两根连杆矢量的叉积定义的轴旋转矢量来建模。
• 位置矢量的扩展双数形式可直接给出速度与加速度，分别为其一阶与二阶导数分量。
• 实现代码使用 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 中的实现提供了一个可重用、高效且精确的工具，适用于球面机构的运动学分析。