0. 基础知识
0.1 线加速度
上一篇中我们讲到了在坐标系{A}和坐标系{B}原点重合时,BQ^BQBQ的速度矢量的表示方式:
AVQ^AV_QAVQ = BARBVQ^A_BR^BV_QBARBVQ + AΩB^A\Omega_BAΩB ×\times× BARBQ^A_BR^BQBARBQ (5.1)
我们可以将左侧改写为如下形式:
ddt\frac{\mathrm{d}} {\mathrm{d}t}dtd(BARBQ^A_BR^BQBARBQ) = BARBVQ^A_BR^BV_QBARBVQ + AΩB^A\Omega_BAΩB ×\times× BARBQ^A_BR^BQBARBQ (5.2)
记住这种形式的方程,会在求解加速度方程时很有用。
接下来推导线加速度的公式,大家一起动手推一下吧!
对式5.1求导,可得当坐标系{A}和坐标系{B}原点重合时,BQ^BQBQ的加速度在坐标系{A}中的表达式:
AV˙Q^A\dot{V}_QAV˙Q = ddt\frac{\mathrm{d}} {\mathrm{d}t}dtd(BARBVQ^A_BR^BV_QBARBVQ) + AΩ˙B^A\dot{\Omega}_BAΩ˙B ×\times× BARBQ^A_BR^BQBARBQ + AΩB^A\Omega_BAΩB ×\times× ddt\frac{\mathrm{d}} {\mathrm{d}t}dtd(BARBQ^A_BR^BQBARBQ) (5.3)
对式5.3右侧的第一项和最后一项代入式5.2,整理后可得,
AV˙Q^A\dot{V}_QAV˙Q = BARBV˙Q^A_BR^B\dot{V}_QBARBV˙Q + 2AΩB^A\Omega_BAΩB ×\times× BARBVQ^A_BR^BV_QBARBVQ + AΩ˙^A\dot{\Omega}AΩ˙ ×\times× BARBQ^A_BR^BQBARBQ + AΩB^A\Omega_BAΩB ×\times× (AΩB^A\Omega_BAΩB ×\times× BARBQ^A_BR^BQBARBQ) (5.4)
将上式扩展到原点不重合的情况,就得到了一般的表达式:
AV˙Q^A\dot{V}_QAV˙Q = AV˙BORG^A\dot{V}_{BORG}AV˙BORG + BARBV˙Q^A_BR^B\dot{V}_QBARBV˙Q + 2AΩB^A\Omega_BAΩB ×\times× BARBVQ^A_BR^BV_QBARBVQ + AΩ˙^A\dot{\Omega}AΩ˙ ×\times× BARBQ^A_BR^BQBARBQ + AΩB^A\Omega_BAΩB ×\times× (AΩB^A\Omega_BAΩB ×\times× BARBQ^A_BR^BQBARBQ) (5.5)
另外要指出,当计算旋转关节时,BQ^BQBQ为常量,那么,
BVQ^BV_QBVQ = BV˙Q^B\dot{V}_QBV˙Q = 0
因此相应地,式5.5就可以简化为:
AV˙Q^A\dot{V}_QAV˙Q = AV˙BORG^A\dot{V}_{BORG}AV˙BORG + AΩ˙^A\dot{\Omega}AΩ˙ ×\times× BARBQ^A_BR^BQBARBQ + AΩB^A\Omega_BAΩB ×\times× (AΩB^A\Omega_BAΩB ×\times× BARBQ^A_BR^BQBARBQ) (5.6)
0.2 角加速度
接下来我们推导角加速度公式。假设坐标系{B}以角速度AΩB^A\Omega_BAΩB相对于坐标系{A}转动,同时坐标系{C}以角速度BΩC^B\Omega_CBΩC相对于坐标系{B}转动。
AΩC^A\Omega_CAΩC = AΩB^A\Omega_BAΩB + BARBΩC^A_BR^B\Omega_CBARBΩC (5.7)
对式5.6求导可得,
AΩ˙C^A\dot{\Omega}_CAΩ˙C = AΩ˙B^A\dot{\Omega}_BAΩ˙B + ddt(BARBΩC)\frac{\mathrm{d}}{\mathrm{d}t}(^A_BR^B\Omega_C)dtd(BARBΩC) (5.8)
同样我们将式5.2带入上式最后一项,可以得到,
AΩ˙C^A\dot{\Omega}_CAΩ˙C = AΩ˙B^A\dot{\Omega}_BAΩ˙B + BARBΩ˙C^A_BR^B\dot{\Omega}_CBARBΩ˙C + AΩB^A\Omega_BAΩB ×\times× BARBΩC^A_BR^B\Omega_CBAR