李亮; 李炜辰; 章佳荣; 王冠; Li Liang; Li Weichen; Zhang Jiarong; Wang Guan

引 en

浮标水下上浮过程建模与仿真 PDF

李亮

机构：

中国船舶工业系统工程研究院,北京 100094

×
，李炜辰

机构：

中国船舶工业系统工程研究院,北京 100094

×
，章佳荣

机构：

中国船舶工业系统工程研究院,北京 100094

×
，王冠

机构：

中国船舶工业系统工程研究院,北京 100094

×

中国船舶工业系统工程研究院,北京 100094；

Modeling and Simulation of Buoy Floating Process Underwater PDF

LI Liang

Affiliation：

Systems Engineering Research Institute,Beijing 100094 ,China

×
， LI Weichen

Affiliation：

Systems Engineering Research Institute,Beijing 100094 ,China

×
， ZHANG Jiarong

Affiliation：

Systems Engineering Research Institute,Beijing 100094 ,China

×
， WANG Guan

Affiliation：

Systems Engineering Research Institute,Beijing 100094 ,China

×

Systems Engineering Research Institute,Beijing 100094 ,China；

作者简介:

李亮(1983-),男,硕士,高级工程师,主要从事装备体系与电子系统研究。

中图分类号:TJ02

文献标识码:A

文章编号:2096-5753(2021)03-0227-06

DOI:10.19838/j.issn.2096-5753.2021.03.011

全文
评论
参考文献
出版信息

参考文献 1

纪金耀，夏天，肖汉华.“迪卡斯”声纳浮标探潜盲区与潜艇的规避分析[J].指挥控制与仿真，2009，31（2）：54-58.

查找原文

参考文献 2

孙明太.航空反潜战术[M].北京：军事科学出版社，2003.

查找原文

参考文献 3

高学强，杨日杰，杨春英.潜艇规避对主动声纳浮标作战效能影响研究[J].系统工程与电子技术，2008，30（2）：301-304.

查找原文

参考文献 4

李颖虹，王凡，任小波.海洋观测能力建设的现状、趋势与对策思考[J].地球科学进展，2010，25（7）：715-722.

查找原文

参考文献 5

吴运发，孙乃宏.反潜的关键设备——声纳浮标的现状和发展[J].声学技术，1999，18（2）：95-96.

查找原文

参考文献 6

梁志城.现代反潜武备[M].北京：海潮出版社，1995.

查找原文

目录contents

摘要Abstract
关键词Keywords
0 引言
1 浮标上浮过程动力学建模
1.1 坐标系与姿态角定义
1.2 动力学方程组
1.3 运动学方程组
2 上浮运动仿真分析
2.1 自由上浮运动
2.2 基准运动轨迹
2.3 速度影响
2.4 正浮力影响
2.5 浮重心距影响
2.6 海流影响
3 结束语
参考文献

摘要

浮标在水下上浮过程中的运动姿态会影响其接收阵在水面的展开效果，进而影响其对空中目标的探测能力。针对该问题，对浮标从开始上浮到上浮至水面的全运动过程进行了建模与仿真，详细分析了初速、重心、浮力等因素对浮标运动轨迹的影响，得到结论如下：在自由上浮的情况下，浮标作垂直上浮运动，在开始阶段速度逐渐增大，最终达平衡速度，并保持此平衡速度直至浮出水面；初速度和平台速度只影响浮标初始阶段的运动，最终其运动参数与基准运动轨迹趋于一致；增大正浮力可使浮标具有更高的出水速度，同时缩短上浮运动时间，而增大浮重心距对浮标运动姿态没有太大影响；海流对浮标上浮姿态影响不大，但会在水平方向产生较大的位移。

Abstract

The motion attitude of the buoy in the process of floating upwards underwater will affect the deployment effect of its receiving array on the water surface，and then affect its ability to detect air targets. In order to solve this problem，this paper models and simulates the whole motion process of the buoy from the beginning of floating upward to the status of surfacing. On analyzing the influences of initial velocity，gravity center，buoyancy and other factors on the trajectory in detail，we have drawn several conclusions as below. Firstly，under the condition of free floating，the buoy moves vertically，and its speed increases gradually in the beginning，while finally reaches the equilibrium speed and moves at this equilibrium speed until reaching the surface. Secondly，the initial velocity and platform velocity only affect the initial motion of the buoy，and finally the motion parameters tend to be consistent with the reference motion trajectory. Thirdly，increasing positive buoyancy can make the buoy more quick at surfacing speed and shorten the floating time，but increasing the distance to the floating center of gravity has little effect on the buoy's motion attitude. Fourthly，the current has little effect on the floating attitude of the buoy，but it will produce great displacement in the horizontal direction.

关键词

浮标；运动轨迹；仿真与建模

Keywords

buoy ； motion trajectory ； simulation and modeling

0 引言
潜艇探测空中目标的常用手段主要有雷达、光电、红外以及声传感器^[1-2]。雷达和光电探测手段只能在潜艇浮出水面或潜望深度时使用，极易暴露自身位置，而红外探测手段无法解决视场角小和水面海浪遮挡等问题^[3-4]。声探测手段是一种在保证自身隐蔽状态下实现对空探测的有效途径，主要有2种方式：一种是利用低频对空警戒声呐对入水的直升机辐射噪声进行探测预警，但由于所处介质不同，所能获取的目标信息极其有限；另一种是利用空气声探测浮标来获取空中目标的方位角、俯仰角和距离，这种方式可获取的信息要比前一种方式丰富，并且精度更高^[5-6]。为保证自身的隐蔽性，潜艇在使用空气声探测浮标时，需要在水下一定深度进行发射，浮标在上浮过程中的运动姿态会影响其接收阵在水面的展开效果，进而影响浮标对空中目标的探测能力。文章针对该问题，对浮标从开始上浮到上浮至水面的全运动过程进行了建模与仿真，详细分析了初速、重心、浮力等因素对浮标运动轨迹的影响，可为浮标的设计、加工提供理论支撑与数据支持。
1 浮标上浮过程动力学建模
1.1 坐标系与姿态角定义
按如图1所示建立地面坐标系 $O x_{o} y_{o} z_{o}$ 、体坐标系 oxyz、平移坐标系 $o x^{'} y^{'} z^{'}$ 。地面坐标系 $O x_{o} y_{o} z_{o}$ 原点 $O$ 置于发射时刻浮标浮心在水面上的投影点， $O x_{o} y_{o} z_{o}$ 坐标平面与静止水面重合， $y_{o}$ 轴垂直向上；体坐标系 oxyz原点位于浮标浮心所在横截面几何中心，x 轴指向浮标头部， y 轴在模型纵平面内， z 轴与 x 轴、 y 轴组成右手坐标系。
图1 模型坐标系
Fig.1 Model coordinate system
浮标在空间的位置由其浮心在地面坐标系内的坐标给出；为直观起见，浮标在空间的姿态仍由欧拉角系统中的3个欧拉角表示。常规的正欧拉角系统定义为：体坐标系从与平移坐标系重合的位置开始依次绕 oy、oz、ox 轴旋转 $ψ 、 θ 、 ϕ$ 3个角度得到的姿态用3个欧拉角组合来表示，分别称为偏航角、俯仰角和横滚角。在文章中，由于浮标垂直上浮，初始俯仰角度为90°，若采用正欧拉角进行描述，在俯仰角超过90°时，偏航角和滚动角均会发生畸变，故采用反欧拉角来描述浮标的姿态。反欧拉角系统定义为：体坐标系从与平移坐标系重合的位置开始依次绕 oz、 oy、 ox 轴旋转 $θ_{r} 、 ψ_{r} 、 ϕ_{r}$ 3个角度得到的姿态用3个欧拉角组合来表示，对应的坐标转换矩阵如式（1）。

$C_{b}^{0} = [\begin{matrix} c o s θ_{r} c o s ψ_{r} & s i n θ_{r} c o s ψ_{r} & - s i n ψ_{r} \\ c o s θ_{r} s i n ψ_{r} s i n ϕ_{r} - s i n θ_{r} c o s ϕ_{r} & c o s θ_{r} c o s ϕ_{r} + s i n θ_{r} s i n ψ_{r} s i n ϕ_{r} & c o s ψ_{r} s i n ϕ_{r} \\ s i n θ_{r} s i n ϕ_{r} + c o s θ_{r} s i n ψ_{r} c o s ϕ_{r} & s i n θ_{r} s i n ψ_{r} c o s ϕ_{r} - c o s θ_{r} s i n ϕ_{r} & c o s ψ_{r} c o s ϕ_{r} \end{matrix}]$

(1)

欧拉方程为

$\{\begin{matrix} θ_{r} = ω_{y} s i n ϕ_{r} s e c ψ_{r} + ω_{z} c o s ϕ_{r} s e c ψ_{r} \\ ψ_{r} = ω_{y} c o s ϕ_{y} - ω_{z} s i n ϕ_{r} \\ ϕ_{r} = ω_{x} + ω_{y} s i n φ_{r} t a n ψ_{r} + ω_{z} c o s ϕ_{r} t a n ψ_{r} \end{matrix}$

(2)

式中， $ω_{x} 、 ω_{y} 、 ω_{z}$ 为浮标角速度分量。
1.2 动力学方程组
根据动量定理和动量矩定理，得到浮标水下运动空间的动力学方程组一般形式为

(3)

式中：为浮标浮心处的线加速度分量；为浮标角加速度分量； $A_{m λ}$ 为惯性矩阵。 $A_{m λ}$ 可表示为

$A_{m λ} = [\begin{matrix} m + λ_{11} & λ_{12} & 0 & 0 & m z_{c} & - m y_{c} + λ_{16} \\ λ_{21} & m + λ_{22} & 0 & - m_{z c} & 0 & m x_{c} + λ_{26} \\ 0 & 0 & m + λ_{33} & m y_{c} + λ_{34} & - m x_{c} + λ_{35} & 0 \\ 0 & - m z_{c} & m y_{c} + λ_{43} & J_{x x} + λ_{44} & - J_{x y} + λ_{45} & - J_{x z} \\ m z_{c} & 0 & - m x_{c} + λ_{53} & - J_{y x} + λ_{54} & J_{y y} + λ_{55} & - J_{y z} \\ - m y_{c} + λ_{61} & m x_{c} + λ_{62} & 0 & - J_{z x} & - J_{z y} & J_{z z} + λ_{66} \end{matrix}]$

(4)

式中： m 为浮标质量； $x_{c}, y_{c} 、 z_{c}$ 为浮标质心位置； $λ_{i j} (i, j = 1,2, L, 6)$ 为附加质量； $J_{i j} (i, j = x, y, z)$ 为浮标转动惯量、惯性矩； $A_{v ω}$ 为速度矩阵； $A_{F M}$ 为力矩阵。 $A_{v ω}$ 可表示为

$A_{v ω} = [\begin{matrix} 0 & - ω_{z} & ω_{y} & 0 & 0 & 0 \\ ω_{z} & 0 & - ω_{x} & 0 & 0 & 0 \\ - ω_{y} & ω_{x} & 0 & 0 & 0 & 0 \\ 0 & - v_{o z} & v_{o y} & 0 & - ω_{z} & ω_{y} \\ v_{o z} & 0 & - v_{o x} & ω_{z} & 0 & - ω_{x} \\ - v_{o y} & v_{o x} & 0 & - ω_{y} & ω_{x} & 0 \end{matrix}]$

(5)

$A_{F M}$ 可表示为

$A_{F M} = [\begin{matrix} X_{α μ} + X_{B} + X_{G} \\ Y_{α μ} + Y_{ω μ} + Y_{B} + Y_{G T} \\ Z_{α μ} + Z_{ω μ} + Z_{B} + Z_{G} \\ M_{α μ x} + M_{ω μ x} + M_{G x} \\ M_{α μ y} + M_{ω μ y} + M_{G y} \\ M_{α μ z} + M_{ω μ z} + M_{G z} \end{matrix}]$

(6)

式中： $X_{B} 、 Y_{B} 、 Z_{B}$ 为浮力在体坐标系下的3个分量； $X_{G} 、 Y_{G} 、 Z_{G}$ 为重力在体坐标系下的3个分量； $M_{G x} 、 M_{G y} 、 M_{G z}$ 为重力矩在体坐标系下的3个分量； $X_{α μ} 、 Y_{α μ} 、 Z_{α μ}$ 为粘性位置力在体坐标系下的3个分量； $M_{α μ x} 、 M_{α μ y} 、 M_{α μ z}$ 为粘性位置力矩在体坐标系下的3个分量； $Y_{ω μ} 、 Z_{ω μ}$ 为粘性阻尼力在体坐标系下的分量； $M_{ω μ x} 、 M_{ω μ y} 、 M_{ω μ z}$ 为粘性阻尼力矩在体坐标系下的3个分量。
1.3 运动学方程组
描述浮标浮心平移运动的方程组可表示为

$[\begin{matrix} x_{o} \\ y_{o} \\ z_{o} \end{matrix}] = C_{b}^{o} [\begin{matrix} v_{o x} \\ v_{o y} \\ v_{o z} \end{matrix}]$

(7)

式中， $v_{o x} 、 V_{o y}, v_{o z}$ 为浮标浮心速度在地面坐标系下的分量。
描述模拟样机转动的运动学方程组为

$\{\begin{matrix} θ_{r} = ω_{y} s i n ϕ_{r} s e c ψ_{r} + ω_{z} c o s ϕ_{r} s e c ψ_{r} \\ ψ_{r} = ω_{y} c o s ϕ_{r} - ω_{z} s i n ϕ_{r} \\ ϕ_{r} = ω_{x} + ω_{y} s i n ϕ_{r} t a n ψ_{r} + ω_{z} c o s ϕ_{r} t a n ψ_{r} \end{matrix}$

(8)

水流攻角 $α$ 、侧滑角 $β$ 、速度 $v$ 可表示为

$α = - a r c t a n (\frac{v_{y}}{v_{x}})$

(9)

$β = a r c t a n (\frac{v_{z}}{\sqrt{v_{x}^{2} + v_{y}^{2}}})$

(10)

$v = \sqrt{v_{x}^{2} + v_{y}^{2} + v_{z}^{2}}$

(11)

2 上浮运动仿真分析
2.1 自由上浮运动
浮标在无海流、无初速影响下进行垂直自由上浮运动的速度变化如图2所示。
图2 垂直自由上浮过程速度变化
Fig.2 Velocity variation during vertical free floating process
初始时刻浮标在正浮力作用下开始加速运动，在t=5s时运动稳定在 v=2m/s。
2.2 基准运动轨迹
假设浮标开始上浮的初始速度为0.3m/s，发射平台的速度为3kn，在此条件下的速度变化如图3所示，运动参数变化如图4所示。浮标开始上浮约7s后各运动参数达到稳定，初始时刻较大攻角，在阻力作用下速度略有减小，正浮力作用使其具有垂直向上加速度，约7s后速度达到最大值，运动过程中模型最大偏转角度为11.5°。
图3 基准运动轨迹下的速度变化
Fig.3 Velocity change in reference motion trajectory
图4 基准运动轨迹下的运动参数变化
Fig.4 Changes of motion parameters in reference motion trajectory
2.3 速度影响
在基准运动轨迹基础上调整浮标初始时刻轴向速度，对 v_x=0.1m/s和 v_x=0.5m/s时的仿真结果进行对比，结果如图5和图6所示。结果表明初始速度只对浮标初始运动过程有一定影响，3s后浮标运动参数与基准运动轨迹趋于一致。
图5 不同初速下上浮过程速度变化
Fig.5 Surfacing velocity variation at different initial velocities
图6 不同初速下上浮过程运动参数变化
Fig.6 Motion parameters change during surfacing at different initial velocities
2.4 正浮力影响
在基准运动轨迹基础上进行衡重参数影响研究，增大模型浮力分别为实际模型的1.5倍和2倍，计算结果如图7和图8所示。随着浮力增大，浮标开始上浮后调整运动姿态时间缩短，可更快稳定于垂直上浮状态，且浮标稳定运动速度增大。
图7 不同正浮力下上浮过程速度变化
Fig.7 Surfacing velocity variation under different positive buoyancy
图8 不同正浮力下上浮过程运动参数变化
Fig.8 Changes of motion parameters during surfacing under different positive buoyancy
2.5 浮重心距影响
在基准运动轨迹基础上进行衡重参数影响研究，将浮重心距由0.03m分别增大为0.1m和0.2m，仿真结果如图9和图10所示。调整浮重心距对上浮姿态没有太大影响，在初始时刻存在简短振荡过程，最终稳定于基准运动轨迹。
图9 不同浮重心距下上浮过程速度变化
Fig.9 Surfacing velocity variation under different floating center-gravity distance
2.6 海流影响
在基准运动轨迹基础上开展海流对姿态影响的研究，分别设定海流为v=–2kn和v=2kn，仿真结果如图11和图12所示。浮标最终保持垂直向上的姿态上浮，但由于海流影响，浮标具有等同于海流大小的法向速度，上浮过程中会随着海流产生位移，如图13所示。
图10 不同浮重心距下上浮过程运动参数变化
Fig.10 Changes of motion parameters during surfacing under different positive buoyancy
图11 不同海流下上浮过程速度变化
Fig.11 Surfacing velocity variation under different currents
图12 不同海流下上浮过程运动参数变化
Fig.12 Changes of motion parameters under during surfacing different currents
图13 不同海流下浮标位移变化
Fig.13 Buoy Displacement variation under different currents
3 结束语
文章开展了空气声探测浮标水下上浮运动过程的建模，对初始条件、衡重参数、海流等因素对上浮过程的影响进行了仿真分析，得到结论如下：
1）在自由上浮的情况下，浮标作垂直上浮运动，在开始阶段速度逐渐增大，最终达平衡速度，并保持此平衡速度直至浮出水面；
2）初速度和平台速度只影响浮标初始阶段的运动，最终其运动参数与基准运动轨迹趋于一致；
3）增大正浮力可使浮标具有更高的出水速度，同时缩短上浮运动时间，而增大浮重心距对浮标运动姿态没有太大影响；
4）海流对浮标上浮姿态影响不大，但会在水平方向产生较大的位移。
参考文献
- [1] 纪金耀，夏天，肖汉华.“迪卡斯”声纳浮标探潜盲区与潜艇的规避分析[J].指挥控制与仿真，2009，31（2）：54-58.
- [2] 孙明太.航空反潜战术[M].北京：军事科学出版社，2003.
- [3] 高学强，杨日杰，杨春英.潜艇规避对主动声纳浮标作战效能影响研究[J].系统工程与电子技术，2008，30（2）：301-304.
- [4] 李颖虹，王凡，任小波.海洋观测能力建设的现状、趋势与对策思考[J].地球科学进展，2010，25（7）：715-722.
- [5] 吴运发，孙乃宏.反潜的关键设备——声纳浮标的现状和发展[J].声学技术，1999，18（2）：95-96.
- [6] 梁志城.现代反潜武备[M].北京：海潮出版社，1995.

参考文献

[1] 纪金耀，夏天，肖汉华.“迪卡斯”声纳浮标探潜盲区与潜艇的规避分析[J].指挥控制与仿真，2009，31（2）：54-58.
[2] 孙明太.航空反潜战术[M].北京：军事科学出版社，2003.
[3] 高学强，杨日杰，杨春英.潜艇规避对主动声纳浮标作战效能影响研究[J].系统工程与电子技术，2008，30（2）：301-304.
[4] 李颖虹，王凡，任小波.海洋观测能力建设的现状、趋势与对策思考[J].地球科学进展，2010，25（7）：715-722.
[5] 吴运发，孙乃宏.反潜的关键设备——声纳浮标的现状和发展[J].声学技术，1999，18（2）：95-96.
[6] 梁志城.现代反潜武备[M].北京：海潮出版社，1995.

首页

本刊简介

编委会

收录情况

期刊订阅

征稿简则

联系我们

浮标水下上浮过程建模与仿真 PDF

Modeling and Simulation of Buoy Floating Process Underwater PDF

摘要

Abstract

关键词

Keywords

0 引言

1 浮标上浮过程动力学建模

1.1 坐标系与姿态角定义

(1)

(2)

1.2 动力学方程组

(3)

(4)

(5)

(6)

1.3 运动学方程组

(7)

(8)

(9)

(10)

(11)

2 上浮运动仿真分析

2.1 自由上浮运动

2.2 基准运动轨迹

2.3 速度影响

2.4 正浮力影响

2.5 浮重心距影响

2.6 海流影响

3 结束语

参考文献

参考文献

首页

本刊简介

编委会

收录情况

期刊订阅

征稿简则

联系我们

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

浮标水下上浮过程建模与仿真 PDF

Modeling and Simulation of Buoy Floating Process Underwater PDF

摘要

Abstract

关键词

Keywords

0 引言

1 浮标上浮过程动力学建模

1.1 坐标系与姿态角定义

(1)

(2)

1.2 动力学方程组

(3)

(4)

(5)

(6)

1.3 运动学方程组

(7)

(8)

(9)

(10)

(11)

2 上浮运动仿真分析

2.1 自由上浮运动

2.2 基准运动轨迹

2.3 速度影响

2.4 正浮力影响

2.5 浮重心距影响

2.6 海流影响

3 结束语

参考文献

参考文献