第53卷第2期2021年4月
Vol.53No.2
Apr.2021南京航空航天大学学报
Journal of Nanjing University of Aeronautics&Astronautics
舟船空投入水姿态与过载特征分析
李正达，郭海军，陈勇，洪淼
（航空工业航宇救生装备有限公司航空防护救生技术航空科技重点实验室，襄阳441003）
摘要：舟船空投入水时将承受较大的垂直和侧向过载，将对舟船的结构强度、稳定性产生不利影响，通过控制船体入水俯仰角可减小此过载。为研究船体入水角度对减载的影响规律，利用LS⁃DYNA对不同角度下船体着水过程进行了仿真，得到不同工况下着水过程加速度变化曲线，并通过水上跌落试验对仿真计算结果进行了验证，结果表明通过采取合理的入水角度，可有效降低船体入水过载。本文所做工作可为舟船空投系统着水减载提供依据。
关键词：流体力学；海上空投；入水减载；空投系统；跌落试验
中图分类号：V24文献标志码：A文章编号：1005⁃2615（2021）02⁃0194⁃08
Analysis of Entry Attitude and Overload Characteristics of Ships in Water
LI Zhengda，GUO Haijun，CHEN Yong，HONG Miao
（Aviation Key Laboratory of Science and Technology on Life⁃Support Technology，
AVIC Aerospace⁃Life⁃Support Industries Ltd.，Xiangyang441003，China）
Abstract:A ship will bear large vertical and lateral overloads when dropped into water，which impos adver effects on the structural strength and stability of the ship.The overloads can be reduced by controlling the pitching angle of the ship.In order to study the influence of the pitching angle on overload reducing，LS-DYNA is ud to simulate the landing process at different angles.The acceleration curves of landing process under different working conditions are obtained，and the simulation results are verified by the fall-down test. The results show that the reasonable entry angle can effectively reduce the ship’s water overloads.The work in this paper can provide the basis for the landing load reduction of the a airdrop system.
Key words:hydromechanics；airdrop on the a；load reduction；airdrop system；fall-down test
当前各国的水面力量除了各类大中型战舰，还有一批不起眼但很重要，也颇具技术含量的小船，即各类战斗快艇，用于执行巡逻、封锁、警戒护卫、搜索救援以及水上特种作战等任务。这类战斗快艇以体积小、吃水浅、机动灵活的特点，广泛活动于全球各地的江河内水和沿海港口。与运输机相结合，采用专用的舟船空投系统，以空投进行远距离、快速投送是其重要的使用方式，是应付远洋地区突发事件，执行远洋搜索救援行动的重要手段。
舟船空投系统与传统陆上空投系统最大差异体现在使用环境上，常规空投系统主要针对陆地着陆场设计，为保证装备和物资器材的空投安全，一般采用气囊、纸蜂窝等方式进行缓冲[1]；水上空投时，舟船经主伞减速稳降后直接着水，主伞脱离后即可进入工作状态。该方式在提升任务效率的同时对着水过载控制也提出了更高的要求。
目前国内外学者对陆上空投系统着陆过程，尤其是气囊缓冲着陆过程研究较多，如Cole等[2]建立了“火星探路者”气囊缓冲系统的解析模型，研究了气囊着陆反弹的动态响应问题。Esgar等[3]在不考
DOI：10.16356/j.1005⁃2615.2021.02.005
收稿日期：2021⁃01⁃20；修订日期：2021⁃03⁃28
通信作者：李正达，男，工程师，E-mail：158****************。
引用格式：李正达，郭海军，陈勇，等.舟船空投入水姿态与过载特征分析［J］.南京航空航天大学学报，2021，53（2）：194⁃201.LI Zhengda，GUO Haijun，CHEN Yong，et al.Analysis of entry attitude and overload characteristics of ships in water ［J］.Journal of Nanjing University of Aeronautics&Astronautics，2021，53（2）：194⁃201.
第2期李正达，等：舟船空投入水姿态与过载特征分析
虑气囊织物材料的弹性及囊内气体质量的前提下，建立了缓冲气囊的解析模型。戈嗣诚和施允涛
[4]
开展了固定排气孔缓冲气囊应用于无人机回收过程的可行性研究，并且探讨了不同气囊参数对缓冲性能的影响。温金鹏等[5⁃6]在考虑气囊织物材料发生弹性变形的情况下，对固定排气孔面积的气囊缓冲特点进行了探讨。徐保成等[7]建立了土壤双线性弹塑性材料模型并利用LS⁃DYNA 对火炮空投非线性着陆过程进行了仿真分析。洪煌杰等[8]对空降车⁃气囊系统着地缓冲过程进行了仿真分析。李建阳等[9]1米每秒等于多少千米每小时
对空降车着陆缓冲过程车体动态应力进行了仿真研究。唐晓慧[10]研究了车载武器系统空投着陆缓冲特性。
气囊数值模拟方法方面，目前常用的耦合计算方法主要有3种：（1）控制体积法；（2）任意拉格朗日欧拉法；（3）粒子法。Welch [11]采用控制体积法对乘员探测飞行器(CEV)的着陆缓冲过程进行了模拟分析。Dmitri 和Nitin [12]分别利用任意拉格朗日法和控制体积法对物体碰撞试验过程进行了模拟计算，对两种算法的计算结果进行了对比分析，发现任意拉格朗日法与试验数据更加吻合。代小芳等[13]利用有限体积法和任意拉格朗日法对折叠气囊折叠气囊展开过程进行了对比研究。
水上空投领域目前国内外研究相对较少，前苏联在其舟船空投系统中曾采用入水姿态控制技术，通过调整舟艇入水角度减少船体入水过载，保障船体着水安全。国内目前尚无成熟的舟船空投系统，但在水上应急救援方面做了一些工作。李名琦[14]建立了应急气囊的有限元模型，对应急气囊着水冲击过程进行了数值计算，并通过缩比模型试验和数值仿真数据进行了对比分析。
目前国内在陆地着陆缓冲和稳定技术已相当成熟，基本可保证空投着陆安全
[15]
，但水上空投着
水减载和稳定技术尚不具备，为有效指导舟船空投系统着水减载设计，有必要对舟船空投入水过程进行仿真分析，研究入水角度对舟船着水过载的影响规律，为水上空投减载研究提供帮助。
1
舟船入水过程仿真
1.1
舟船着水动力学方程介绍1.1.1基本假设
为了简化舟船运动方程的推导，作如下基本假设：
（1）舟船为刚体，不考虑船体弹性变形和旋转部件的影响；
（2）大气条件为国际标准大气；（3）地球为平面大地，忽略曲率和旋转。1.1.2
坐标系及转换矩阵
为了建立舟船运动方程，引入以下几种常用坐标系。各坐标系均为右手直角坐标系。
（1）地面坐标系Ox g y g z g
定义：固定于地球表面，原点O 位于回收初始时刻舟船质心在水面上的投影点；
Ox g 轴指向舟船初始运动方向；
Oy g 轴铅垂向下；Oz g 轴垂直于Ox g y g 平面，按右手定则确定。
（2）机体坐标系Ox b y b z b
定义：原点O 固联于舟船质心，Ox b 轴在舟船对称平面内，平行于船体轴线，指向前；Oz b 在对称平面内，垂直于Ox b 轴，指向下；
Oy b 轴垂直于对称平面指向右。
（3）气流坐标系Ox a y a z a
定义：原点O 位于舟船质心，Ox a 始终指向舟船空速方向；Oz a 轴位于对称平面内，垂直于Ox a 轴，指向下；
Oy a 轴垂直于Ox a z a 指向右。坐标系Ox a y a z a 通过按一定的顺序先绕y b 轴转过角-α，再绕当时的z 轴转过角β，就可以与重合。相应的转换矩阵为
L ab =éëê
ê
êùû
úúúcos αcos β
sin βsin αcos β-cos αsin β
cos β-sin αsin β-sin α
cos α
(1)
机体坐标系Ox b y b z b 相对于地面坐标系Ox g y g z g 的方位常用3个欧拉角表示，即：偏航角ψ，俯仰角θ，滚转角ϕ。
相应的Ox b y b z b 到Ox g y g z g 的转换矩阵为
L bg =é
ëê
ê
êùû
úú
úcos θcos ψcos θsin ψ
-sin θsin θsin ϕcos ψ-cos ϕsin ψsin θsin ϕsin ψ+cos ϕcos ψ
sin ϕcos θsin θcos ϕcos ψ+sin ϕsin ψ-sin ϕcos ψ
sin θcos ϕsin ψ
cos ϕcos θ(2)
1.1.3动力学方程
（1)舟船质心动力学方程
在机体坐标系Ox b y b z b 下，根据牛顿第二定律有
F b =m
d V b
d t
(3)式中：
F b 为船体所受合外力矢量F b =
G b +A b +T b
(4)
195
第53卷
南京航空航天大学学报式中：G b ,A b ,T b 分别为重力、气动力、连接绳拉力
矢量。由于重力、气动力分别定义在地面坐标系和气流坐标系下，所以其坐标转换关系如下
G b =L bg G g =L bg éëêêêêù
ûúúúú00mg (5)
A b =L ba A a =L T ab
éëêêêêùû
úúúú-D C -L (6)
式中：D ，C ，L 分别为舟船所受的阻力、侧力和升力。
质心的绝对加速度可表示为
d V b d t =δV b
δt
+ωb ×V b
(7)
联立即可得到机体坐标系下质心动力学方程为
L bg éëêêêêùû
úú
úú00mg +L T ab éëêêêêù
ûúúúú-D C -L +T b =m (δV b δt +ωb ×V b )(8)
（2）舟船绕质心转动动力学方程
在机体坐标系Ox b y b z b 下，根据动量矩定理
M b =δh b δt
+ωb ×h b
(9)
式中：
M b 为作用于船体质心的合外力矩，ωb 为绕质心的角速度矢量。机体坐标系下有
M b =éëêêêêù
û
úúúúL M N =M A +M T (10)
式中：
M A ,M T 分别为作用于船体的气动力矩和连接绳拉力产生的力矩。
连接绳拉力对质心产生力矩为
M T =r P ×T b
第三次长沙会战
(11)
式中：
r P 为舟船质心到连接绳拉力作用点的的矢径，
T b 为连接绳拉力。由于ωb =[p q r ]T
（其中p ，q ，r 分别表示3个方向上的角速度分量），故动量矩h b 可表示为
h b =éëêêêù
û
úúú
pI x -qI xy -rI zx qI y -pI xy -rI yz rI z -pI zx -qI yz (12)
式中I x ,I y ,I z 为分别舟船对x ,y ,z 轴的惯性矩，I xy ,I yz ,I zx 为惯性积。对于一般舟船，Ox b y b 平面通
常为对称面，此时I xy =I yz =0。
联立即可得船体绕质心转动的动力学方程如下
ìíîïïïï
ïïïïI x d p
d t +(I z -I y )qr -I zx (pq +d r d t )=L I y
d q
d t +(I x -I z )rp +I zx (p 2-r 2)=M I z d r d t
+(I y -I x )pq +I zx (qr -d p d t )=N
(13)（3）舟船质心运动学方程
在地面坐标系下，舟船速度矢量即为空间坐标的微分，即
V g =éëê
d x g妈妈不是我的佣人
d t d y g d t d z g d t ùû
úT
(14)
对V g 进行坐标变换可得
L bg éëêêêêù
û
úúúúu v w =éëêd x g
d t d y g d t
d z g d t ùû
úT
(15)
（4）舟船绕质心转动运动学方程带叶的诗句
根据机轴系形成过程可以写出旋转角速度在机体轴系上的投影为
éëêêêùû
úúúp q r =éëêêêêùûúúúúϕ00+éëêêêêùûúúúú0θcos ϕ-θsin ϕ+L bg éëêêêêùûúúúú00ψ(16)
解得
ìí
îïïïïïïïïd ϕ
d t =p +tan θ(q sin ϕ+r cos ϕ)d θ
d t =q cos ϕ-r sin ϕd ψd t
=1cos θ(q sin ϕ+r cos ϕ)
(17)
1.2
舟船入水过程建模
选取蜘蛛船作为研究对象，采用三维造型软件对蜘蛛船模型进行处理，保留部件整体外形及尺寸，去除较小的倒角、孔缝等特征，将几何模型导入Hypermesh 软件中进行网格划分，网格数量366万
个。蜘蛛船网格模型如图1所示。
流场域的大小直接决定了网格的数量，过大的流场域会导致计算资源的浪费，造成计算时间过长或求解器无法承受的大容量计算等问题[16]；过小的水域会导致计算结果受到边界的反射波的影响，因
此流场域在与蜘蛛船耦合的区域采用与船相似的
图1蜘蛛船网格模型Fig.1
Spider boat grid model
196
第2期李正达，等：舟船空投入水姿态与过载特征分析
网格尺寸，超出耦合区域的流场网格采用渐疏网格的方式处理，如图2所示。流场域尺寸27m×13m×8m ，建立的空气域和水域如图3所示。
水域周围，右侧采用全约束约束条件，底侧采取竖直方向位移约束条件，左侧采取无反射边界条件，前后两侧采取垂向位移约束条件，从而保证水域边界计算以及波浪稳定产生。1.3
求解参数设置
LS⁃DYNA 通过设置关键字进行求解参数设置[17]，其中主要涉及的关键字类型如下：
（1）部件参数设置。包括每个部件所用的材料、算法等参数。
（2）重力加速度设置。设置受到重力影响的物体。
（3）耦合控制关键字。本次计算主要考虑蜘蛛船与水之间的耦合作用。
（4）输出控制关键字。设置需要输出的主要参数。1.3.1
部件参数设置
部件主要包括蜘蛛船、空气域和水域，以蜘蛛船为例，主要参数设置如图4所示。
1.3.2
重力加速度设置
本次计算中主要考虑蜘蛛船的重力作用，重力加速度设为参数设置如图5所示。
1.3.3耦合控制关键字设置
蜘蛛船与水之间通过罚函数法进行耦合，相关参数如图6所示。
1.3.4输出控制关键字设置
为了在计算结果中提取蜘蛛船的加速度结果
和耦合力的结果，设置rbdout 和dbfsi 关键字进行输出，如图7所示。
1.4仿真计算结果
定义船体与水平方向的夹角为θ，计算不同夹
角情况下蜘蛛船着水过程中所受到的过载峰值，计算工况见表1，计算结果见表2
。
图2流场域网格图Fig.2
Field
grid
图3空气域和水域图Fig.3
Air and water
areas
图4蜘蛛船部件参数设置图Fig.4
Parameters of spider
ship
图5重力加速度设置图
Fig.5
Gravity acceleration
ttings
图6耦合控制关键字设置Fig.6
Coupling control key
t
图7
输出控制关键字设置Fig.7
Output control key t
197
第53卷
南京航空航天大学学报工况1蜘蛛船的加速度⁃时间曲线如图8所示，过载峰值48g ，着水过程如图9所示。
工况2蜘蛛船的加速度⁃时间曲线如图10所示，过载峰值6.2g ，着水过程如图11所示。
工况3蜘蛛船的加速度⁃时间曲线如图12所示，过载峰值2.8g ，着水过程如图13所示。
1.5
仿真结果分析从计算结果可以看出：
（1）水平着水过程历时约20ms ，远小于正常
着陆缓冲过程（约200ms ）。其着水过载很大，可达到48g （着陆过载一般不大于20g ），因此必须考虑舟船入水减载问题。
（2）随着船体与水平面夹角增大，其着水过程历时有大幅增加（30°时已达到400ms 以上），最大着水过载则相应减小。由此可见，通过调整舟船入水角度，可将其缓冲持续时间延长，使其过载变化趋于平缓，从而起到降低着水冲击载荷的目的。从
能量的角度来看，该措施使得能量经更长时间耗散，从而起到保护船体的作用。
表2
计算结果
Table 2
Calculation results
计算工况
123
与水平方向的夹角θ/（°）
03070
最大着水过载/g
486.22.8
表1
计算工况表
Table 1
Calculation conditions
计算工况
123
与水平方向的夹角θ/（°）
03070
着水垂直速度/（m·s -1）
888
着水水平速度/（m·s -1）
00
图8工况1加速度⁃时间曲线图
Fig.8
Acceleration -time curve of condition
1
图9工况1着水过程图信用等级查询
Fig.9
Entry process of condition冰火两重天是什么意思
1
图10工况2加速度⁃时间曲线图
Fig.10
Acceleration -time curve of condition
2
图11工况2着水过程图Fig.11
difficult的反义词Entry process of condition
2
图12工况3加速度⁃时间曲线图
Fig.12
Acceleration -time curve of condition
3蛋糕照片
图13
工况3着水过程图
Fig.13
Entry process of condition 3
198
