船舶操纵1要点.pdf
《船舶操纵1要点.pdf》由会员分享,可在线阅读,更多相关《船舶操纵1要点.pdf(13页珍藏版)》请在三一文库上搜索。
1、船舶操纵运动仿真报告 指导老师:王先洲 班级:船海 1303 小组成员:徐田丰U201312305 桂强 U201312311 张琳 U201312296 李坚宇U201312316 丁逸轩U201312323 王睿哲U201212135 2016年 5 月 7 日 一国内外研究现状 船舶操纵性是指船舶在驾驶者的操纵下,通过控制装置来保持或改变航速、 航向和位置等运动状态的能力, 他是船舶的运动性能之一。 它直接影响到船舶的 经济性和安全性。航运业的不断发展促进了船舶的大型化 多样化和高速化航 道航行密度的增加导致海上事故发生的概率上升使得人们对船舶操纵性的要求 提高。 船舶操纵运动仿真作为一
2、种在设计阶段常用的预报船舶操作性能的方法,在 现代计算机迅速发展,多媒体进步的今天,这种方法更是有了深厚的发展前景。 目前这种技术在航海人员的训练中得到了广泛地引用;并且在港口航道的设计和 管理上有一定作用。 同时在各工况下对船舶运动的仿真研究,各船型下运动的仿 真研究,各个辅助装置对船舶运动影响的仿真研究都多如繁星。 如大型船舶进坞的仿真: 大型油轮进出坞时依靠船坞引船系统或拖船拖曳运 动,由于船坞的坞口以及进出坞航道狭窄,使油轮的进出坞运动控制较困难,事故 风险较高。本文比较了油轮进出坞运动与深海域拖航系统运动的差异,从理论出 发,建立了油轮进坞运动方程 ,利用四阶 Runge-Kutta
3、 方法求解一阶微分方程组 ,得 到了油轮进坞的运动轨迹及瞬时速度。仿真结果表明,该油轮在较好的海况条件 下进坞是安全的。 如风浪流作用下的船舶操纵运动仿真:将 OpenGL 虚拟现实仿真技术运用到 船舶操纵与控制领域中, 利用 MMG 数学模型,建立了风浪流作用下船舶操纵运 动的三维动态仿真系统, 在系统开发过程中, 首先采用 MMG 分离式数学模型及 相关的系列化试验结果,建立单桨单舵海洋运输船舶在静水中的船舶操纵运动 方程,并编制计算程序,经与试验结果比较,证实了计算结果的正确性;为了解 MMG 数 学模型中模型参数变化对操纵性指数的影响程度在上述已有程序基础 上,对有关模型参数进行偏移修
4、正,探讨了相应参数变化后的操纵性指数,对 船舶操纵性指数对模型参数的灵敏度进行了详细的分析与探讨。 如半潜船的六维自由度仿真:采用分离建模的思想,通过对船舶运动以及受 力进行分析 ,建立了一个对半潜船在外力风、浪、流作用下的船舶运动数学模型 进行了 6 自由方向上的建模。并对船舶进行受力分析,建立了 6 自由度裸体船体 水动力及力矩、 6 自由度流体动力及力矩和6 自由度外力干扰力及力矩等计算模 型和吊舱螺旋桨推进器的推力以及推力矩的计算模型。最后通过计算机进行仿真, 利用数值分析方法四阶Runge-Kutta算法计算出船舶 6 自由度运动状态。 目前主流的数学模型是mmg,这种模型通过水动力
5、系数、桨和舵的相关系 数对水动力进行建模,因此数学模型与上面的系数有很大的关系。 而本次我们采用abkowitz 模型,这种模型需要船舶的所有水动力导数,进 而计算出该模型在各种操舵规律下, 各种运动参数随时间变化的曲线和模型的运 动轨迹。 二船舶运动方程 由于回转过程的速降、 各种激烈的操纵运动, 船舶操纵的线性模型因其局限 性而无法准确地预报各种情况下的操纵运动。因此,为了对船舶操纵运动进行更 为精确的数学模拟,一般采用非线性数学模型。 在此介绍 Abkowitz 非线性数学模型。该模型将方程中的流体动力函数X、 Y、N 在直航状态附近对各运动参数做泰勒级数展开,即把XYN 表示为各种状
6、态变量及耦合项的泰勒级数之和。 urXuvX uvrXrXvXvrmXuX uvXurXXrm XX vXuXuXuXXuvrf uruv vrurvvru vvurruGrr vvuuuuuu )( 2 1 2 1 2 1 2 1 ) 2 1 ( 2 1 6 1 2 1 )( 2 2222 232 01 uYvrY uYuYrYvYY YurYurYrYrvYrY rm uYruvYuvYvYvrY vYvYuYuYuYYuvrf uvr uuurrvv ruururrvvrrr vuuvuvvrr vvvvpouuupoupoup 3 2223 2223 222 332 02 6 1 2
7、1 2 1 2 1 6 1 2 1 2 1 2 1 6 1 )( 2 1 2 1 2 1 6 1 )( uNvrN uNuNrNvNN NurNurNrNrvNrN rum xNruvNuvNvNvrN vNvNuNuYuNNuvrf uvr uuurrvv ruururrvvrrr Gvuuvuvvrr vvvvpouuupoupoup 3 2223 2223 222 332 03 6 1 2 1 2 1 2 1 6 1 2 1 2 1 2 1 6 1 )( 2 1 2 1 2 1 6 1 )( 由上式可知,常系数非线性方程组的解如下: u Xm rvuf u dt du)(1、 )Y-)(
8、mxN-(mx-)N-)(IY-(m )()()()( rGvGrzv 32、rvufYmxrvufNI v dt dvrGrz )Y-)(mxN-(mx-)N-)(IY-(m )()()()( rGvGrzv 23、rvufNmxrvufYm r dt dvvGv 上式即为一阶微分方程组, 可利用计算机迅速求得其近似解。只要一直前一瞬时 t 的 u,v,r,值, 便可解出该瞬时的 rvu, 值。 进而推算出下一瞬间的速度状态变量。 因此微分后可以求得船舶相对于固定坐标系的线速度值,即为绝对速度: )(sin)()(cos)()( 0 ttvttutx )(cos)()(sin)()( 0 t
9、tvttuty 对上式进行积分: tt r trt 0 )()0()( tvuxtx tt r0 00 )(sin)()(cos)()0()( tvuyty tt r0 00)(cos)()(sin)()0()( 其中: )0()0()0( 00 yx、 为初始时刻状态值。 由各瞬时之 )()( 00 tytx、 值就可以绘 制其轨迹。如果船体重心G 与动坐标原点 O 不重合,则需要进行适当的修正从 而得到正确的重心轨迹。 以上即 Abkowitz 非线性模型的简要内容, 如果某船的全部水动力导数已知, 则运用该模型可以方便的计算出其在各种操舵情况下,运动参数随时间的变化及 其轨迹,由实际实验
10、证明, 该数学模型预报值和实船实验值符合较好,所以为了 计算某船的操纵运动,确定其水动力导数是非常重要的。Abkowitz 数学模型中 众多水动力导数可以通过各种约束船模实验而测得,而阻力实验和螺旋桨敞水试 验可以提供另一部分系数值。 三仿真结果及分析 本次测试的船型为mariner Matlab 代码如下 function xdot = mariner( x,ui ) if (length(x) = 7),error(x-vector must have dimension 7 !);end if (length(ui) = 1),error(u-vector must have dimen
11、sion 1 !);end U0 = 7.72; U = sqrt(U0+x(1)2+x(2)2); L = 160.93; delta_c = ui(1); u = x(1)/U; v = x(2)/U; delta = x(7); r = x(3)*L/U; psi = x(4); delta_max = 15; Ddelta_max = 5; m = 798e-5; Iz = 39.2e-5; xG = -0.023; Xudot=-840e-5;Yvdot=-1546e-5;Nvdot=23e-5; Xu=-184e-5;Yrdot=9e-5;Nrdot=-83e-5; Xuu=-11
12、0e-5;Yv=-1160e-5;Nv=-264e-5; Xuuu=-215e-5;Yr=-499e-5;Nr=-166e-5; Xvv=-899e-5;Yvvv=-8078e-5;Nvvv=1636e-5; Xrr=18e-5;Yvvr=15356e-5;Nvvr=-5483e-5; Xdd=-95e-5;Yvu=-1160e-5;Nvu=-264e-5; Xudd=-190e-5;Yru=-499e-5;Nru=-166e-5; Xrv=798e-5;Yd=278e-5;Nd=-139e-5; Xvd=93e-5;Yddd=-90e-5;Nddd=45e-5; Xuvd=93e-5;Yud
13、=556e-5;Nud=-278e-5; Yuud=278e-5;Nuud=-139e-5; Yvdd=-4e-5;Nvdd=13e-5; Yvvd=1190e-5;Nvvd=-489e-5; Y0=-4e-5;N0=3e-5; Y0u=-8e-5;N0u=6e-5; Y0uu=-4e-5;N0uu=3e-5; m11=m-Xudot;m23=m*xG-Yrdot;m33=Iz-Nrdot; m22=m-Yvdot;m32=m*xG-Nvdot; if abs(delta_c)=delta_max*pi/180, delta_c=sign(delta_c)*delta_max*pi/180;e
14、nd delta_dot=delta_c-delta; if abs(delta_dot)=Ddelta_max*pi/180, delta_dot=sign(delta_dot)*Ddelta_max*pi/180; end X=Xu*u+Xuu*u2+Xuuu*u3+Xvv*v2+Xrr*r2+Xrv*r*v+Xdd*delta2+Xudd*u*delt a2+Xvd*v*delta+Xuvd*u*v*delta; Y=Yv*v+Yr*r+Yvvv*v3+Yvvr*v2*r+Yvu*v*u+Yru*r*u+Yd*delta+Yddd*delta3+Yud *u*delta+Yuud*u2*
15、delta+Yvdd*v*delta2+Yvvd*v2*delta+(Y0+Y0u*u+Y0uu*u2) ; N=Nv*v+Nr*r+Nvvv*v3+Nvvr*v2*r+Nvu*v*u+Nru*r*u+Nd*delta+Nddd*delta3+ Nud*u*delta+Nuud*u2*delta+Nvdd*v*delta2+Nvvd*v2*delta+(N0+N0u*u+N0u u*u2); xdot=X*(U2/L)/m11 -(-m33*Y+m23*N)/(m22*m33-m23*m32)*U2/L (-m32*Y+m22*N)/(m22*m33-m23*m32)*U2/L2 r*U/L
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 船舶 操纵 要点
链接地址:https://www.31doc.com/p-5211841.html