水下滑翔机稳态螺旋运动的时空尺度分析

于鹏垚，沈 聪，甄春博，王天霖

(大连海事大学 船舶与海洋工程学院，辽宁 大连 116026)

0 引 言

水下滑翔机是一种新型的水下航行器，其通过改变自身净浮力和质心位置来实现水下滑翔运动。与传统AUV相比，水下滑翔机具有成本低、能耗低、续航力强等优点，更适合大范围、长航程的工作任务，不足之处在于水下滑翔机的航速低、机动性较差[1,2]。水下滑翔机可通过调节尾舵舵角或改变横滚姿态实现转向，相比于传统的AUV，其在转向调节过程通常需要较大的空间和较长的时间[3,4]。尤其对于利用横滚姿态进行转向调节的水下滑翔机，所需的时间和空间尺度更大。将滑翔机的转向速度作为摄动参数，Mahmoudian等[5]推导了滑翔机稳态转向时转向速度的近似解析表达式，Yang等[6]推广此方法分析了装有尾舵的滑翔机转向性能。顾建农等[7]基于柯西霍夫方程实现了水下滑翔机的动力学建模和转向运动仿真。朱时雨等[8]分析了水柱测量中的水下滑翔机转向性能。为降低内部机构调节的能耗，在大多数工作时间内，水下滑翔机处于稳态运动。水下滑翔机的稳态螺旋运动，作为一种典型运动形式，常用于水下滑翔机航向的调节。Zhang等[9]建立了一种稳态螺旋运动的快速求解方法，并与Seawing滑翔机的试验结果进行对比验证。赵宝强等[10]采用李雅普诺夫方法进行了螺旋运动稳定性分析。Liu[11]采用理论与试验方法分析了水下滑翔机的反螺旋运动。

当水下滑翔机在近海底或近壁面等限制水域运动时，水域的空间范围会限制其转向运动。水下滑翔机的转弯运动通常可视为稳态螺旋运动的一部分，因此，稳态螺旋运动在时间、空间上的运动尺度特征可以为水下滑翔机在限制水域内的转向运动提供指导。目前，在现有研究中尚未看到螺旋运动时空尺度特征的系统分析，文献[9]也仅是给出了下潜过程中回转半径与滑块转角的变化规律，并未分析其他控制参数和尺度特征。本文从水下滑翔机的运动学和动力学微分方程入手，建立了稳态螺旋运动的时间、空间尺度特征变化规律的研究思路，并以沈阳自动化研究所研制的Seawing水下滑翔机[9,12]（该水下滑翔机通过内部滑块旋转来实现转向）为例开展研究。

1 运动学建模

引入惯性坐标系 E-ξηζ和体坐标系o-xyz来描述水下滑翔机的空间运动。如图1所示，坐标轴 E-ζ与重力方向一致，坐标系o -xyz的 原点位于滑翔机的浮心，o-x轴沿着滑翔机的机身轴线指向首部方向， o-y轴指向右侧机翼的方向。定义水下滑翔机在惯性坐标系下的位置和姿态分别为和定义水下滑翔机在体坐标系下的线速度和角速度分别为则水下滑翔机的运动学方程可表达为[13]：

图1 坐标系与质量分布Fig. 1 Coordinate systems and mass distribution

其中，REB和RBΩ可用水下滑翔机的姿态角来表达。采用简 化形式 c·=cos(·)， s·=sin(·)， t·=tan(·)表达三角函数，则有

2 动力学建模

水下滑翔机可以视为不同种结构组成的多体系统，包括机身及与机身固结在一起的所有结构mrb，可在机身内部平移和旋转的滑块结构mp，位置不变质量可变的浮力调节机构mb。mrb，mp，mb的质心在体坐标系中的位置分别为向量rrb， rp，rb。rrx为mp的质心沿着o-x 轴的位置，γ为mp绕 o-x轴的转角。其中，γ=0对应mp的质心处在体坐标的o-xz平面内，此时滑翔机将在垂直平面内做直线运动；γ≠0将导致滑翔机所受的升力偏离垂直平面，升力在水平面中的分量将作为向心力使其产生螺旋转向运动，如图2所示。

图2 螺旋转向运动的力学机理Fig. 2 Mechanisms of turning spiraling motion

水下滑翔机在水下运动时，会带动其周围的附连水共同运动，水下滑翔机及其附连水系统的总体动能可以表达为：

其中：M为系统的广义质量矩阵；MA， JA和CA分别为体坐标系下附连水的质量矩阵，转动惯量矩阵和耦合项矩阵；为体坐标系固定质量的转动惯量矩阵；为单位对角矩阵。

进一步表达有

其中：G 和B分别为水下滑翔机的重力和浮力，rCG为体坐标系下的重心位置向量，FB和MB分别为体坐标系下水下滑翔机的粘性水动力和水动力矩。对于FB和MB，可结合滑翔机运动速度和水动力系数来计算[9,13]。

将式（6）、式（7）代入式（8）、式（9），即得到水下滑翔机稳态运动的动力学方程：

3 稳态螺旋运动的时空描述

基于水下滑翔机的运动学方程和动力学方程，当给定其控制输入量和 b ，Ω，υ，ω的初值后，即可实现水下滑翔机的三维空间运动仿真。对于依靠质心移动实现转向的水下滑翔机，其运动姿态的控制输入量包括移动滑块沿轴的位置、移动滑块绕轴的转角和净浮力质量mb。本文以Seawing水下滑翔机为例进行分析，水下滑翔机自身的质量参数和水动力参数见文献[9]。将运动状态量 b，Ω，υ，ω的初值取为零向量，并将控制参数取值为 γ=20°，rrx=0.42m，mb=0.3 kg，得到水下滑翔机的运动轨迹、运动速度变化规律分别如图3和图4所示。从图3可以看出，可采用回转半径，螺距和回转周期来表达水下滑翔机的空间和时间尺度特征。其中，回转半径 R的方向与滑翔机绕 o-ζ轴的角速度方向一致，螺距 H的方向与滑翔机沿 o-ζ轴的线速度方向一致，回转周期 T为滑翔机一个周期回转运动的时间。若水下滑翔机需要转向某一角度，可以利用转向角数值和螺旋运动的时空尺度特征，换算得到对应转向运动所需的时间、空间量。

图3 水下滑翔机的螺旋运动轨迹Fig. 3 The spiral motion trajectory of an underwater glider

图4 运动速度随时间的变化规律Fig. 4 Variation of the velocity with the time

从图4可以看出，随着时间的变化，水下滑翔机的运动速度 υ， ω将逐渐稳定，相应的水下滑翔机运动也将逐渐达到稳态。为进一步定量地验证本文稳态运动仿真方法的正确性，将本文回转半径的分析结果与文献[9]中的结果进行对比，如图5所示，其中, mb=0.3 kg。可以看出，2种方法结果的变化趋势一致，但量值上存在较小的差别，原因为文献[9]中为采用递归方法实现稳态螺旋运动的快速分析，在推导过程中忽略了质量矩阵中影响较小的耦合项，而本文方法从滑翔机的运动微分方程入手，更为详细地计入了质量矩阵。总体看来，本文的运动仿真方法可用于水下滑翔机稳态螺旋运动时空尺度的分析。

图5 本文方法的验证Fig. 5 Validation of the method presented in this paper

4 时空尺度的分析

根据水下滑翔机沿垂向的运动方向进行分类，水下滑翔机的稳态螺旋运动可分为上浮螺旋运动和下潜螺旋运动。当重力大于浮力，即mb＞0时，滑翔机将进行下潜螺旋运动，此时需调节 rrx使水下滑翔机的重心在浮心的前方，从而保证滑翔机沿着机头的方向滑行。而当重力小于浮力，即mb＜0时，滑翔机将进行上浮螺旋运动，此时需调节 rrx使滑翔机的重心在浮心的后方。以Seawing水下滑翔机为例进行研究，通过改变滑翔机下潜和上浮运动中的控制参数，系统地设计一系列稳态运动计算工况，进而分析水下滑翔机的时间、空间尺度特征随控制输入参数间的变化规律。

在下潜和上浮运动中保持 rrx和mb的取值不变，通过改变进行稳态运动的计算，可得螺旋运动的时空尺度特征随转角的变化规律，如图6所示。其中，, mb=0.3 kg，对应下潜螺旋运动； rrx=0.42m,mb=0.3 kg，对应上浮螺旋运动。可以看出，回转半径R 关于坐标原点中心对称，其方向与转角的方向一致；螺距 H 和回转周期关于位置左右对称。随着转角 γ绝对值的增大，上浮和下潜螺旋运动中回转半径 R 、螺距 H 和回转周期的绝对值均减小。这表明在有限尺度的空间中转向时，可考虑增大来降低转向过程中所需的时间、空间尺度。

在下潜和上浮运动中保持 γ和mb的取值不变，通过改变 rrx进行稳态运动的计算，可得到时空尺度特征随 rrx的变化规律，如图7所示。其中， mb=0.3 kg,γ=90°和mb=-0.3 kg, γ =90°，分别对应下潜螺旋运动和上浮螺旋运动。可以看出，在下潜螺旋运动中，回转半径、螺距和回转周期 T均存在最大值；在上浮螺旋运动中，回转半径 R存在最大值，而且上浮过程中的回转周期要小于下潜过程。 rrx在0.4～0.405 m附近时，水下滑翔机的重心和浮心的纵向位置比较接近。尽管在最值位置两侧的单调范围中调节 rrx，均可降低所需的时空尺度，但总体来看，当调节 rrx来减小水下滑翔机的重心和浮心的纵向距离（使 rrx接近0.4～0.405 m范围），更有助于快速地降低螺旋转向过程中的时间、空间尺度。

图6 时空尺度随转角的变化规律Fig. 6 Relationship between the temporal and spatial scales and the rotation angle

图7 时空尺度随纵向位置的变化规律Fig. 7 Relationship between the temporal and spatial scales and the longitudinal position of the movable block

图8 时空尺度随净浮力的变化规律Fig. 8 Relationship between the temporal and spatial scales and the net buoyancy

5 结 语

本文通过对水下滑翔机稳态螺旋运动的时间、空间尺度进行分析，得到如下结论：1）建立水下滑翔机的稳态螺旋运动仿真方法，并验证了该方法的正确。2）系统地分析了稳态螺旋运动的时空尺度特征随控制输入参数间的变化规律，发现了上浮和下潜运动中变化规律的不同之处。3）当通过调节移动滑块的转角或纵向位置来降低螺旋运动中的时间、空间尺度时，可考虑增大滑块转角、减小水下滑翔机的重心和浮心的纵向距离来实现。4）当通过调节净浮力来降低螺旋运动中的时间、空间尺度时，在下潜过程中可考虑降低净浮力，在上浮过程中应根据转向时刻mb的具体量值来调节。