吊舱推进船舶运动数学模型及其操纵性仿真
2024-10-18
来源:威能网
第33卷第06期 文章编号:1006—9348(2016)06—0138—05 计算机仿真 2016年6月 吊舱推进船舶运动数学模型及其操纵性仿真 朴在吉,郭晨 (大连海事大学信息科学技术学院,辽宁大连1 16026) 摘要:研究吊舱推进(POD)船舶的操纵性问题,由于吊舱式推进器是一种新型的船舶电力推进系统,是利用推进电机驱动 可旋转螺旋桨推动船舶前进,所以建模方法和传统的船舶有所不同。针对上述问题,应用MMG建模机理,给出模型中多种 流体动力和力矩的合理计算方法,建立了该吊舱船舶操纵运动的综合仿真模型。采用吊舱推进的半潜船“泰安口”轮的实 船参数,在MATLAB平台下进行仿真,通过回转试验和z形试验仿真结果分析了吊舱船舶的操纵性。仿真结果表明所建模 型可较好地反映吊舱推进船舶操纵运动特性。 关键词:吊舱推进器;船舶操纵性;回转试验;操纵试验 中图分类号:TP391.9 文献标识码:B Mathematical Modeling and Maneuverability Simulation for POD——Driven Ship PIA0 Zai—ii.GU0 Chen (Dalian Maritime University,Information Science and Technology College,Dalian Liaoning 116026,China) ABSTRACT:The problem of POD~driven ship maneuverability is studied.The podded propulsor(POD)is a new type of ship electric propulsion system,which uses propulsion motor to drive rotatable propeller for moving the ship, SO the modeling method is different from the traditional ship.In view of the above problems,a MMG modeling mech— anism is applied to present reasonable calculation methods in a variety of fluid force and torque,and comprehensive simulation model of the POD—driven ship’s maneuvering motion is established.Based on the actual parameters of semi—submersible heavy lift vessel“TaiAnKou”,the method is simulated in MATLAB platform,and the POD— driven ship maneuverability is analyzed through the simulation results of the turning test and zig—zag test.The simu— lation results show that the model can well reflect POD—driven ship maneuvering motion characteristic. KEYWORDS:Podded propulsor;Ship maneuverability;Turning test;Maneuver test 电力推进技术中的亮点,它将使船舶电力推进的优越性得到 1 引言 随着船舶工业的不断发展,2O世纪90年代,出现了吊舱 更充分的体现。 国内外学者对于吊舱存在的问题进行了细致的研究。 国外,日本学者Tomihiro Haraguehi等采用船舶对航向的保持 和改变能力来评价吊舱推进船舶的操纵性。英美等国家也 (POD1推进方式,其作为一种新型的电力推进方式在经济 性、操纵性、节能性等诸多方面有着其特别的优越性 J。它 利用柴油机、燃气轮机或涡轮发电机组产生电能,再通过电 动机将电能转换成机械能,驱动螺旋桨推动船舶运动。POD 推进器提高了推进器的水动力特性,它提出了一种集推进和 操舵装置于一体的全新推进理念,极大地增加了船舶设计、 建造和使用的灵活性。POD推进器的出现和应用是舰船全 对吊舱有着较深入的研究 I4 J。国内方面,马骋 以POD 推进器为研究对象,运用数值计算和模型试验两种方法进行 了一系列研究,哈尔滨工程大学和大连海事大学等一些高校 也进行了对于吊舱推进器的研究工作,并取得了一些成 果 J。但是,吊舱推进器依然是很新的技术,在数学模型的 建立和其操纵性仿真等方面依然存在问题,其水动力参数和 基金项目:国家自然科学基金项目(61374114);中央高校基本科研业 务费专项资金资助项目(3132014321) 收稿日期:2015—07—20修回日期:2015—09—21 —吊舱的运行方式对于模型仿真结果的的准确性有很大的 影响。 为了更好的处理上述问题,本文主要应用MMG模型对 138一 装有一对SSP吊舱推进装置的”泰安口”号进行数学建模,给 船舶受到的纵向力、横向力和转艏力矩。下标为日的表示船 出模型中需要用到的多种流体动力和力矩的参数的计算方 法,建立了该吊舱船舶合理的运动数学模型。并基于MAT— LAB中的Simulink进行了较为精确的仿真,对其操纵性进行 研究。 体粘性类流体动力和力矩,下标为P的表示吊舱产生的力和 力矩。 。和 为船舶相对于大地的纵向和横向位移。吊舱 产生的力和力矩以及它和船体之间的互相影响问题一直是 建模的难点,本文给出了一种MMG模型的改进方法用以处 理这一问题。 2 吊舱推进船舶运动操纵系统原理 吊舱式推进器利用发电机把其它形式的能量转变成电 能,再通过电动机把电能装换成机械能,从而推动船舶的运 3 吊舱推进船舶的MMG建模方法 3.1坐标系 动,实现了能量的非机械模式传递。传统的船舶的推进方式 大多是由船舶尾部的螺旋桨产生向前或向后的推力,由舵产 生向左或向右的侧推力,从而实现船舶在水面的运动(如图 1(a)),而采用吊舱推进方式的船舶则改变了这一模式,在船 舶尾部装有滑环装置(由电动或液压马达来带动),能使吊舱 360度回转,从而起到了舵的作用(如图1(b))。采用吊舱 推进方式的船舶在同样回转角的情况下可以比采用舵来转 向的船舶产生更大侧向力和回转力矩,在转向的过程中纵向 的推进力有所减少,所以吊舱推进方式船舶的运动模型较传 统船舶的模型复杂度有所提高 j。本文以安装有2个sie— mens Schottel Propulsor,(SSP)推进器(如图1(C))的“泰安 口”号半潜船为例进行建模研究。 一 l . (a)水面运动 (b)舵作用 (c)半潜船 图1 传统船舶与吊舱船舶受到的推力方向 对于吊舱船舶,尾部没有舵,所以整体的受力分为裸船 体和吊舱的流体动力和力矩。按照MMG的建模机理得到吊 舱船舶的运动方程为 (m+m ) 一(m+m ) r=X +XP (m+m ) +(m+m )ur= +y 0I +J =NH+NP (1) =r 0=UCOS 一vsin ̄ =VCOS + 嘶 方程中m、mx、m 为船体质量和船舶在纵向和横向的附 加质量; 和 为船舶惯性矩和附加惯性矩; 、y、Ⅳ分别为 本文研究的是船舶在海洋上三自由度的运动,一般采用 两种坐标系统:惯性坐标系统与附体坐标系统(如图2)。图 中,惯性坐标系统为固定于地球表面的坐标系统,取作基准 参考系统,规定‰轴指向正北,y0轴指向正东;附体坐标系 为原点在船舶满载吃水平面的前后左右对称点上的一种随 船坐标系,规定 轴指向船首,Y轴指向右舷 。“为船舶纵 向速度; 为船舶横向速度;r为转首角速度; 为航向角;6 为舵角。 x |图2坐标系 3.2附加质量和附加惯性矩 吊舱推进船舶的运动数学模型表达式已在上文给出,其 中重要的一组参数为附加质量和附加惯性矩,这是由于船舶 作加速运动时,船体作用于周围的水,使之得到加速度,所以 船舶运动时除了要克服自身的惯性,还要克服水的惯性力。 流体惯性力的作用效果相当于物体的质量和惯性矩均增加 了一个数值,所以称为附加质量和附加惯性矩。附加质量和 附加惯性矩利用周昭明回归公式得到。周昭明回归公式为 周昭明对日本著名的元良图谱进行多元回归分析得出的公 式,见式(2)。其中, 为船长;B为船宽;d为吃水;C 为方 形系数。 一139~ m =1 ̄[0.398+11.97Cb( +s., 吾)-2.89Cb吉(1+1.13鲁)+0.175Cb(告) ( +。.s4 吾)-1.107百L m =m[0.882—0.54c (1—1.6吾)一o.156(1—0.673c )音+0.826百d百L(1—0.678 d)一0.638c 百d百L(1一o.619 ( /m)n = m[33—76.85c 784C6)+3.43吉(1_o_63c )] (2) 3.3船体粘性流体动力和力矩计算方法 为吊舱有效来流速度,W为伴流系数。 3.5伴流系数和推力减额系数的计算 船舶在航行过程中,附近的水受到船体的影响而产生运 目前能用于近似估算粘性流体动力的模型有两种,一种 是井上模型;另一种是贵岛模型,本文采用后者。贵岛模型 如下 fXH:X +X +X v "r+x rr { : + r+】,… l I +y…,I rI r+l,坤l r+ vr2 【Nn= +Ⅳ,r+Ⅳ… l } +Ⅳ… l rI r+Ⅳ r+Ⅳ r (3) 上述方程的各个水动力导数的计算方法请参考文献 [8],这里不再赘述。 3.4 吊舱力和力矩计算方法 本文研究的是安装有一对SSP吊舱推进器的船舶,即是 船尾有两个吊舱,每个吊舱有两个螺旋桨,其对于船舶运动 提供的力和力矩应当分为两部分,一部分为螺旋桨产生的推 力 ,另一部分为吊舱舱体本身在转向时产生的横向力Ⅳl9 J, 计算公式如下 rXP=2(1一tP)(Tp+ )COS 6~2(1一tR)NsinS {yP=2(1一t )( + )sin 6+2(1+ⅡH)Ncos ̄ (4) tNe=2(1一t )( 一 )C03 aO.5L 一 L0P 式中6为吊舱的转向角; 。为左右螺旋桨之间的横向距离; 。 吊舱距船舶中心的距离 。’” ; 和 分别为左右吊舱 螺旋桨产生的推力;tp为推力减额系数;t 为舵阻力减额系 数;a 为计入操舵诱导船体横向力后关于舵力的修正因子。 tR=0.2618+0.0539C^一0.1755C (5) aH=0.6784—1.3374C +1.8991C (6) 因为本文考虑转向时吊舱本身产生的横向力,所以相关 参数的计算方法参考舵的流体动力模型,t 和a 采用式(5) 和式(6)的近似公式计算。 71P=pn2P上,4 (.,P) (7) Kr=0.7+0.3589Jp+0.1875 (8) 式(7)中P为海水密度;n 为左螺旋桨转速;D为螺旋桨直 径;Jp=Up/(npD)为进速比; 为根据吊舱推进器的敞水性 能图谱,采用曲线拟合和插值方式求出的推力系数,具体公 式见式(8)。 的求法和 一致。 1 . Ⅳ= A sin(。 ) (9) 厶 up=(1一W)[U +(口+ Pr) ] (10) 式(9)中,A 为吊舱有效面积,吊舱的面积虽然不像舵那么 大,但是其产生的横向力在数值上并不能忽略。a = 一占 为吊舱有效攻角, :arctan[一( +Xpr)/u]。式(10)中 .--——140---—— 动,称其为伴流或迹流。伴流使螺旋桨产生的推力与敞水桨 也不同。当船舶纵向运动时,本文中的伴流系数采用对双螺 旋桨船的泰洛公式计算: W0=0.55C6+0.20 (11) 当船舶横向和旋转运动时,采用平野雅祥提出的一个实 用模型 W=Woexp(~4 ) (12) 为螺旋桨位置处横流速度,计算公式为 =12 +l rpr (13) 2 大概和螺旋桨位置的无量纲 坐标一致,取为0.5L。 螺旋桨在船后工作时,它的抽吸作用使桨盘前方的水流 速度增大。根据伯努利定理,该处压力必然下降,使螺旋桨 推力有所减少。推力额减系数采用对于双桨商船的汉克歇 尔公式 tp=0.5Cp一0.18 (14) C 为船舶的菱形系数¨ 。 4实验结果及分析 4.1船舶参数 本文以“泰安口”轮为例进行仿真,该船采用两套SSP电 力推动系统力推进系统。驱动较为轻巧的螺旋桨来实现推 动船舶行驶,从而代替了传统大型柴油机通过很长的尾轴直 接推动船舶前进的方式。这种推动系统装置在远洋船舶使 用,属于国际首例。运转中的螺旋桨可以360度回转,航行 时通过操纵螺旋桨回旋的角度,可以实现控制船舶前进的方 向,使船舶可以在极小的回旋半径和范围内灵活操纵。船舶 可以自行离靠码头和非常灵巧地自装自卸货物 。船舶主 尺寸如表1所示。 表1“泰安口”船舶主尺寸 船舶总长L (in) 两柱间长L什(m) 船宽 B(m) 设计吃水 d(m) 方形系数 cb 设计排水 m(Ton) 棱形系数 C “泰安口”轮SSP推进器采用的是同轴双桨式设计,主要 参数见表2。 表2螺旋桨参数 4.2基于Simulink的模型实现 利用MATLAB中的Smulink可以方便的实现船舶的建 模。用三个输入,六个输出的模式,三个输入分别是左螺旋 桨转速;吊舱转向角; 右螺旋桨转速,六个输出量为横向速度u,纵向速度v, 转首角速度r,航向角1lr,横向位移x,纵向位移Y,输出方程 采用前文提到的吊舱船舶运动方程。 本文搭建的仿真模型如图3所示。两个吊舱采用联动 方式,即给定两个吊舱一样的转向角。在S函数中设置状态 变量的初始值和船舶的基本参数,通过一系列的计算得到各 个力和力矩,从而得到上述各个输出量。 ≯ 卜 ,—一 ———————_- 图3基于Simu/iak的仿真框图 4.3旋回实验结果及操纵性分析 假设船舶在理想海面上航行,给定转向角为35。,即进行 全速右回旋实验,左右螺旋桨转速均为140r/min。初始速度 为15.3kn。仿真曲线如图4所示。 将转向角变为一35。,即进行左回旋实验,仿真曲线如图 5所示。船舶操纵性是船舶重要的航行性能之一,是船舶按 驾驶人员的意图并借助控制装置来改变或保持其运动速率、 姿态和方位的性能。船舶的回转性与船舶避让,离靠泊,灵 活调头等操纵密切相关。 图4右回旋实验船舶运动轨迹 / _一 — 、 。\ , 、 \ ; f 、 \ :;; :/+ y/m 图5左回旋实验船舶运动轨迹 从回转圈看,吊舱船舶的操纵性较传统船舶有较大的提 高,与各类船舶相对回转直径的比较见表3。 表3“泰安口”轮与各类船舶D0/L统计比较 船舶种类 最大舵角时定常回转D0/L 大型货船 中小型货船 油轮 中小型海上客船 “泰安口”轮 从表4中的仿真模型与实船的各项回转圈参数的对比 来看,该模型能基本反映“泰安口”轮的运动情况。模型可以 仿真出船舶在左右螺旋桨转速不一致时的运动姿态,以及航 向保持能力等等。但由于数据的不足等原因,还是存在一定 的误差。 一141— 表4“泰安口”轮右旋回圈数据 传统船舶的舵桨产生的力和力矩变为吊舱作用在船舶上的 力和力矩。基于吊舱推进的半潜船“泰安口”轮的实船参数, 仿真计算完成了船舶回转试验和z形试验,得到了船舶的仿 真运动轨迹,验证了所建吊舱推进船舶模型的可靠性。为进 一步研究POD推进船舶的控制问题奠定了初步的操纵数学 模型基础。 船舶的z形操纵试验是船舶操纵试验的一种,用来评价 船舶的艏摇抑制性能。同时也可以用z形试验的结果求取 船舶的操纵性指数K和T。通常可以根据强机动性和弱机 参考文献: [1] 咸屹,聂文天.吊舱式全回转电力推进器的现状及展望[J]. 江苏船舶,2008,24(6):28—29. 动性情况进行5。/5。、10。/10。、20。/20。的z形操纵试验,但 是通常是以10。/10。Z形操纵试验为准 。本文在高速满 载的情况下进行10。/10。(分子表示舵角,分母表示船舶进行 反向操舵时的航向角)的Z形试验,试验船舶轨迹结果如图 6所示。 甚 图6 Z形试验结果 表5 z形试验运动指标 分析上述图6和表5可得:船舶在做10。/10。Z形试验 时,船舶的仿真试验轨迹几乎 对称,符合理论上的研究推理,而船舶的第一超越角大 概为1.7。。第二超越角大概也为1.7。。相对于船舶的实际 尺度和船舶的结构来讲均在理论研究范围之内,进一步验证 了仿真模型的真实性 J。 5结论 POD推进器作为一种新型的船舶电力推进系统,提高了 推进器的水动力特性,实现了将推进和操舵装置集成于一体 的全新船舶推进理念。 作者应用MMG建模机理,以装备POD推进器的实船为 例,建立了该吊舱船舶操纵运动的综合仿真模型。建模时在 处理螺旋桨推力时做出了针对于吊舱船舶特性的修正,即把 一】42一 [2]M F Islam,et a1.Numerical Study of Hub Taper Angle on Podded Propeller Performance[J].Marine Technology,2006,43(1):1— 10(1O). [3]T Haraguchi,T Nimura.A Study on Maneuvrability Standards for a Ship with a Pod Propulsion System[C].International Conference on Marine Simulation and Ship Maneuverability,MARSIM.2003 —3. [4] R Shamsi,et a1.Numerical hydrodynamic evaluation of propeller (with hub taper)and podded drive in azimuthing conditions[J]. Ocean Engineering,2014,76(2):121—135. [5] 马骋.POD推进器的水动力性能研究[D].哈尔滨工程大 学。2006. [6] 董涛,王志新.吊舱推进器研究开发与应用新进展[J].船舶 工程,2007,29(1):65—68. [7] 马骋.吊舱推进技术[M].上海:上海交通大学出版社,2007. [8]贾欣乐,杨盐生.船舶运动数学模型:机理建模与辩识建模 [M].大连:大连海事大学出版社,1999. [9] 刘洁,李彦.吊舱推进船舶操纵性仿真研究[J].科学技术与 工程,201l,11(3):507—512. [10] 梁新莉.POD推进方式大型船舶的运动建模与仿真[D].大 连海事大学,2oo8. [11] 郭晨,雷正玲.海洋动力定位分层控制系统的建模与仿真研 究[J].系统仿真学报,2014,(5). [12]赵红,郭晨,吴志良.船舶电力推进系统的建模与仿真[J]. 中国造船,2006,47(4):51—56. [13] 廖再文.亚洲第一船——中远航运新型半潜船泰安口轮启航 [J].中国远洋航务公告,2003—1:018. [14]王珍旺.大型双桨双舵船舶操纵性及智能控制的仿真研究 [D].大连海事大学,2014. [15]奚庆潮.船舶航向自抗扰控制的研究[D].大连海事大 学,2014. [16]汪有良.半潜船运动建模与仿真[D].武汉理工大学,2011. 啊 【朴人郭士真在生,,硕智能控制理论与应用。导吉晨士师研(1,9究研50生6究一,方研),向究男作为方(汉者船朝向族舶 简为鲜)船介族自,江舶动]) 苏,电控内如力制蒙东推系古人进统扎,。 教兰系授屯统,市博仿