潜水器的六自由度运动模型

潜水器的六自由度运动模型就是对潜水器列质心运动定理和相对于质心运动的动量矩定理,由于

  1. 传感器测量值与潜水器位姿相关
  2. 潜水器受到的水动力与其位姿关系紧密
  3. 执行机构与潜水器固联

为了方便分析,在潜水器运动模型建立时,引入了随体坐标系。

随体坐标系的定义

采用国际水池会议(ITTC)推荐的和造船与轮机工程学会(SNAME)术语公报的体系。坐标是下列两种右手系:一个是固定坐标系 E − ξ η ζ E-\xi\eta\zeta Eξηζ固定于地球;另一个是运动坐标系 G − x y z G-xyz Gxyz,固联于潜水器,随其一起运动。

潜水器运动与受力投影至随体坐标系下,用如下符号表示。

矢量 x 轴 x轴 x y 轴 y轴 y z 轴 z轴 z
速 度 V 速度V V u u u v v v w w w
角 速 度 Ω 角速度\Omega Ω p p p q q q r r r
力 F 力F F X X X Y Y Y Z Z Z
力 矩 M 力矩M M K K K M M M N N N

运动与外力关系

将质心运动定理与动量矩定理,由随体坐标系下的投影量表示出来
{ m ( u ˙ + q w − r v ) = X m ( v ˙ + r u − p w ) = Y m ( w ˙ + p v − q u ) = Z I x p ˙ + ( I z − I y ) q r = K I y q ˙ + ( I x − I z ) r p = M I z r ˙ + ( I y − I z ) p q = N \begin{cases} m(\dot{u}+qw-rv)=X\\ m(\dot{v}+ru-pw)=Y\\ m(\dot{w}+pv-qu)=Z\\ I_x\dot{p}+(I_z-I_y)qr = K\\ I_y\dot{q}+(I_x-I_z)rp = M\\ I_z\dot{r}+(I_y-I_z)pq = N \end{cases} m(u˙+qwrv)=Xm(v˙+rupw)=Ym(w˙+pvqu)=ZIxp˙+(IzIy)qr=KIyq˙+(IxIz)rp=MIzr˙+(IyIz)pq=N

外力的组成

潜水器受到的外力主要有两部分

  • 流体动力:由于潜水器在水中运动受到的流体对于水的反作用力
  • 非流体动力:重力和浮力
  • 执行机构附加力(舵力、推进器推力)

为了简化研究,将潜水器在实际流体中作非定常运动时所受的流体动力分为由于惯性引起的惯性类水动力和由于粘性引起的非惯性类水动力两类来考虑,并忽略其相互影响。

潜水器受到的惯性力
[ B 1 B 2 B 3 B 4 B 5 B 6 ] = [ λ 11   λ 12   λ 13   λ 14   λ 15   λ 16 λ 21   λ 22   λ 23   λ 24   λ 25   λ 26 λ 31   λ 32   λ 33   λ 34   λ 35   λ 36 λ 41   λ 42   λ 43   λ 44   λ 45   λ 46 λ 51   λ 52   λ 53   λ 54   λ 55   λ 56 λ 61   λ 62   λ 63   λ 64   λ 65   λ 66 ] ∗ [ u v w p q r ] \left[ \begin{matrix} B_1\\ B_2\\ B_3\\ B_4\\ B_5\\ B_6 \end{matrix}\right]= \left[ \begin{matrix} \lambda_{11}\ \lambda_{12}\ \lambda_{13}\ \lambda_{14}\ \lambda_{15}\ \lambda_{16}\\ \lambda_{21}\ \lambda_{22}\ \lambda_{23}\ \lambda_{24}\ \lambda_{25}\ \lambda_{26}\\ \lambda_{31}\ \lambda_{32}\ \lambda_{33}\ \lambda_{34}\ \lambda_{35}\ \lambda_{36}\\ \lambda_{41}\ \lambda_{42}\ \lambda_{43}\ \lambda_{44}\ \lambda_{45}\ \lambda_{46}\\ \lambda_{51}\ \lambda_{52}\ \lambda_{53}\ \lambda_{54}\ \lambda_{55}\ \lambda_{56}\\ \lambda_{61}\ \lambda_{62}\ \lambda_{63}\ \lambda_{64}\ \lambda_{65}\ \lambda_{66}\\ \end{matrix}\right]* \left[\begin{matrix} u \\ v \\ w \\ p \\ q \\ r \\ \end{matrix}\right] B1B2B3B4B5B6=λ11 λ12 λ13 λ14 λ15 λ16λ21 λ22 λ23 λ24 λ25 λ26λ31 λ32 λ33 λ34 λ35 λ36λ41 λ42 λ43 λ44 λ45 λ46λ51 λ52 λ53 λ54 λ55 λ56λ61 λ62 λ63 λ64 λ65 λ66uvwpqr
应用动量矩和动量定理
{ − X i f = B ˙ 1 + q B 3 − r B 2 − Y i f = B ˙ 2 + r B 3 − p B 1 − Z i f = B ˙ 3 + p B 2 − q B 1 − K i f = B ˙ 4 + ( q B 6 − r B 5 ) + ( v B 1 − w B 2 ) − M i f = B ˙ 5 + ( r B 4 − p B 6 ) + ( w B 1 − u B 3 ) − N i f = B ˙ 6 + ( p B 5 − q B 4 ) + ( u B 2 − v B 1 ) \begin{cases} -X_if=\dot B_1+qB_3-rB_2\\ -Y_if=\dot B_2+rB_3-pB_1\\ -Z_if=\dot B_3+pB_2-qB_1\\ -K_if=\dot B_4+(qB_6-rB_5)+(vB_1-wB_2)\\ -M_if=\dot B_5+(rB_4-pB_6)+(wB_1-uB_3)\\ -N_if=\dot B_6+(pB_5-qB_4)+(uB_2-vB_1) \end{cases} Xif=B˙1+qB3rB2Yif=B˙2+rB3pB1Zif=B˙3+pB2qB1Kif=B˙4+(qB6rB5)+(vB1wB2)Mif=B˙5+(rB4pB6)+(wB1uB3)Nif=B˙6+(pB5qB4)+(uB2vB1)
潜水器受到的粘性力
{ X v f = f X ( u , v , w , p , q , r ) Y v f = f Y ( u , v , w , p , q , r ) Z v f = f Z ( u , v , w , p , q , r ) K v f = f K ( u , v , w , p , q , r ) M v f = f M ( u , v , w , p , q , r ) N v f = f N ( u , v , w , p , q , r ) \begin{cases} X_{vf}=f_{X}(u,v,w,p,q,r)\\ Y_{vf}=f_{Y}(u,v,w,p,q,r)\\ Z_{vf}=f_{Z}(u,v,w,p,q,r)\\ K_{vf}=f_{K}(u,v,w,p,q,r)\\ M_{vf}=f_{M}(u,v,w,p,q,r)\\ N_{vf}=f_{N}(u,v,w,p,q,r)\\ \end{cases} Xvf=fX(u,v,w,p,q,r)Yvf=fY(u,v,w,p,q,r)Zvf=fZ(u,v,w,p,q,r)Kvf=fK(u,v,w,p,q,r)Mvf=fM(u,v,w,p,q,r)Nvf=fN(u,v,w,p,q,r)
将上式在等速直航基准运动为泰勒级数展开点展开,考虑潜水器对称性、运动特点、使用操纵性数学模型和水动力系数对操纵运动影响的大小,予以必要简化处理。
展开后的水动力系数可以通过模型试验或者CFD数值模拟得到。


将上述方程联立可以得到潜水器的控制量 u u u与潜水器运动的微分方程。该方程即为潜水器的六自由度运动方程。

  • 2
    点赞
  • 29
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值