声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2409|回复: 0

[其他相关] 刚体姿态的数学表达之:对偶四元数

[复制链接]
发表于 2022-2-15 14:43 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
对偶数和对偶矢量
工程技术中涉及的刚体运动往往移动和转动同时发生,需要有表达刚体位置和姿态同时改变的更具普遍性的数学工具。1873年,英国数学家克利弗德 (Clifford,W K.) 创造出对偶数 (dual number) 的数学概念(图1)。对偶数与复数概念类似,是由实数单位1和对偶单位ɛ 组成的包含两个实元的一对有序组合[1。ɛ 称为Clifford 算符,被赋予特殊的运算规则:
1.png
图1 克利弗德 (Clifford,W K. 1845–1879)
用带顶标“˄”的字符表示对偶数,以 2.png 3.png 为例,令
4.png
基于Clifford 算符的基本规则,可推出对偶数的加法和乘法运算规则:
5.png
对偶数中的实数集合如改为矢量集合,则称为对偶矢量 (dual vector)。将式
6.png
中的标量换作矢量a、b,构成对偶矢量和:
7.png
对偶矢量的运算遵循对偶数的加法和乘法运算规则,以及矢量的加法、点积和叉积规则:
8.png

普吕克矢量
设a 为空间中的任意矢量,可沿作用线自由滑动。从空间中的固定点O 向矢量的作用线上任一点作矢量ra,将ra与a 的叉积定义为a 对O 点的矩,记作a ꞌ=ra×a。1869年德国数学家普吕克 (Plücker,J.) 提出(图2),矢量a 与矢量矩a ꞌ 二者的结合可完全确定矢量a 在三维空间内的位置。a 与a ꞌ 组成的矢量偶对称为普吕克矢量 (Plücker vectors),所对应的6个投影坐标称为普吕克坐标 (Plücker coordinates)。设矢量a 在空间中作任意位移到达矢量b 的位置,利用Clifford 算符ɛ,将矢量a 和矢量b 的位置以对偶形式的普吕克矢量表示为
9.png
图2 普吕克 (Plücker,J. 1801-1868)

在单位矢量a 和b 的公垂线方向作p 轴,以p 为p 轴的基矢量。设a 和b 的垂直距离为l,则a 沿p 轴移动l 距离后与b 的作用点重合于P 点,与b 的夹角为θ。可以认为矢量b 是矢量a 沿p 轴平移l 距离,且绕p 轴转动θ 角后到达的位置(图3)。1830年,法国物理学家沙勒 (Chasles,M.) 将这种运动称为螺旋运动 (screw motion)。他同时证明,刚体在空间中的任何位置移动均能通过绕某根轴的螺旋运动实现。换言之,螺旋运动是刚体运动的最普遍形式。

利用Clifford 算符ɛ 将转角θ 与平移l 组合成对偶数 10.png
11.png
称为a与b之间的对偶角,作为螺旋运动的数学描述。
12.png
图3 矢量a和b的平移和转动

将上式定义的对偶角的正弦和余弦函数展成泰勒级数,由于Clifford 算符ɛ 的特殊性质,ɛ 的二次以上的所有项均等于零。导出
13.png
计算式
14.png
的普吕克矢量 、之间的点积和叉积。利用rb=ra+lp,以及矢量的点积和叉积公式:
15.png
遵循对偶数的运算规则,利用适当的矢量变换,导出
16.png
其中, 17.png 为p的普吕克矢量
18.png
将式
19.png
与式
20.png
对比,可看出对偶矢量与一般矢量的点积和叉积规则在形式上完全相同。直接验算还可证明,对偶角也存在类似的三角恒等式,例如
21.png
由此可见,借助抽象的Clifford 算符,能使矢量代数和平面三角公式以同样的形式,发展为对偶矢量和对偶角的计算公式。

对偶四元数
将矢量b 沿正交的a 和p×a 方向分解(图3),得到
22.png
给定p 轴的基矢量p 和矢量a 绕p 轴的转角θ,即可利用上式算出矢量a 绕P 点作有限转动后的位置b。对于刚体作螺旋运动的更一般情形,矢量b 的空间位置应以式
23.png
定义的普吕克矢量 表示。将
24.png
代入
25.png
rb=ra+lp,遵循对偶数的运算规则,且利用适当的矢量变换公式,化作
26.png
依据式
27.png
将下式中括号内的变量用相应的对偶数和对偶矢量代替
28.png
得到螺旋运动的变换公式[2
29.png

30.png
对比可看出,将有限转动的变换公式 (上式) 中各字符增加顶标“˄”,即转换成螺旋运动变换公式。二者在形式上完全一致,仅上式中的矢量a、b 和转角θ 分别以普吕克矢量、和对偶角代替。

将对偶数 31.png 与对偶矢量 32.png 的组合定义为对偶四元数(dual quaternions),记作 36.png
33.png

34.png
可表示为
35.png
其中的空心圆点“ο”表示四元数的乘法运算。

对偶四元数也称为螺旋算子 (screw operator),是有限转动张量A 的扩展概念。利用螺旋算子可以很方便地确定矢量作任意螺旋运动后到达的新位置。多次螺旋运动对应的螺旋算子等于各次螺旋算子的四元数乘积:
37.png
刚体在空间中有6个自由度。仿照欧拉角的定义,刚体的任意空间位置变动也可分解为绕不同坐标轴的螺旋运动,所对应的3个对偶角可称为对偶欧拉角。所包含的6个标量可作为表达刚体位置和姿态的独立变量。

以上利用对偶数和四元数对刚体螺旋运动的推导结果,充分显示出数学表达的简洁和优美。利用对偶四元数处理刚体的螺旋运动在工程技术中,尤其在空间机构和机器人技术中有广泛的实际应用[3。用对偶形式描述的力学基本定律具有最简练的形式,在此基础上建立起来的运动学和动力学称为力学的旋量理论。

参考文献:
[1] Clifford M K. Preliminary sketch of biquaternions. Proc. London Math. Soc., 1873, 4: 381-395.
[2] Wittenburg J. Duale Quaternionen in der Kinematik rämlicher Getriebe, Eine anschauliche Darstellung. Ingenieur Archiv, 1981, 51: 17-20.
[3] Yang A T. Freidenstein F. Application of dual-number quaternion algebra to the analysis of spatial mechanisms. J. of Appl. Mech., ASME, Ser.E, 1964, N2.

原文注:改写自“刘延柱. 关于刚体姿态的数学表达. 力学与实践,2008,30(1):98-101

来源:刘延柱科学网博客,作者:刘延柱。

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-18 19:23 , Processed in 0.076023 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表