船舶运动控制系统建模应用

1、第6章船舶运动控制系统建模应用6.1弓I言数学模型化(mathematical modelling)是用数学语言(微分方程式)描述实际过程动态 特性的方法。在船舶运动控制领域，建立船舶运动数学模型大体上有两个目的：一个目 的是建立船舶操纵模拟器(ship manoeuvring simulator)，为研究闭环系统性能提供一个基本的仿真平台；另一个目的是直接为设计船舶运动控制器服务。船舶运动数学模型主要 可分为非线性数学模型和线性数学模型，前者用于船舶操纵模拟器设计和神经网络控制器、模糊控制器等非线性控制器的训练和优化，后者则用于简化的闭环性能仿真研究和 线性控制器(PID, LQ, LQG,

2、 H :鲁棒控制器)的设计。船舶的实际运动异常复杂，在一般情况下具有6个自由度。在附体坐标系内考察，这种运动包括跟随3个附体坐标轴的移动及围绕3个附体坐标轴的转动，前者以前进速度(surge velocity) u、横漂速度(sway velocity) v、起伏速度(heave velocity) w 表述，后者以 艏摇角速度(yaw rate) r、横摇角速度(rolling rate) p及纵摇角速度(pitching rate) q表述；在 惯性坐标系内考察，船舶运动可以用它的3个空间位置 x0,y0,z0(或3个空间运动速度x0,y0,z0)和3个姿态角即方位角(heading an

3、gle) *、横倾角(rolling angle):、纵倾角(pitchi ng an gle)二(或3个角速度.', -)来描述，('-：U 称为欧拉角(见图 u,v,wT x0,y0,ZgT p,q,rT , ;:T，我们要把这6个自由度上的运动全部加以考虑。 数学模型是实际系统的简化，如何简化就有很大学问。太复杂和精细的模型可能包含难 于估计的参数，也不便于分析。过于简单的模型不能描述系统的重要性能。这就需要我 们建模时在复杂和简单之间做合理的折中。对于船舶运动控制来说，建立一个复杂程度 适宜、精度满足研究要求的数学模型是至关重要的。图0 X0Y0Z0， 为起始位置， O

4、X0指向正北，OYo指向正东，OZ0指向地心;o-xyz是附体坐标系，为船首尾之间连线的中点，ox沿船中线指向船首，oy指向右舷，oz指向地心；航向角'以正北为零度，沿顺时针方向取0360 ；舵角：以右舵为正。对于大多数船舶运动及其控制问题而言，可以忽略起伏运动、纵摇运动及横摇运动，而只需讨论前进运动、横漂运动和艏摇运动，这样就简化成一种只有3个自由度的平面运动问题。图船舶平面运动模型对于像航向保持、航迹跟踪、动力定位、自动避碰等问题，具有足够的精度；但在研究像舵阻摇、大舵角操纵等问题时，则必须考虑横摇运动。本章根据刚体动力学基本理论建立船舶平面运动基本方程，据此进一步导出状态空间型(

5、线性和非线性)及传递函数型船舶运动数学模型，并考虑了操舵伺服系统的动态特性和风、浪、流干扰的处理方法。这些结果将作为设计各种船舶运动控制器的基础。计及横摇的四自由度船舶运动数学模型参见文献5。1)惯性坐标系及与之相关的速度分量取O -XoY。为固定于地球的大地坐标系，原点0设为船舶运动始点或任取，地球的曲率在此可不考虑，不过在涉及大范围航行的航线设计问题时，需单独处理。设船舶运动速度向量V在OXo方向上的分量为Uo，V在OYo方向上的分量为Vo，船舶当前的位置是(Xo,y°)，时间变量以t表示，显有(6-2-1)Xo(t)-Xo(°)巳仏叭 yo(t) yo(O) = Ov

6、odt设船舶的艏摇角速度r顺时针方向为正，有(t)(0) = ；rdt(6-2-2)2) 附体坐标系及与之相关的速度分量取附体坐标系oxy位于满载水线面内。船舶运动速度V在ox方向上的分量为u,称为前进速度，V在oy方向上的分量为 V叫做横漂速度。 同一速度向量V在惯性坐标系的分量(Uo，v。)及附体坐标系的分量(u,v)有下列明显的关系uoco®' -sinuVo sin cos v(6-2-3)3) 两种坐标系内运动学变量之间的关系在惯性坐标系内船舶的位置和姿态由Xo (t), yo (t)；-' (t)T确定，在附体坐标系内船舶之运动速度和角速度由u(t), v

7、(t),r(t) T表示。由式(6-2-1)，式(6-2-2)和式(6-2-3)知*(t)(0厂 0r(t)dtxo(t) =xo(o)十 f)u(t)cosW(t) v(t) sin即(t)dt，yo(t) =yo(o)+ fu(t)s in 即(t) +v(t) cos即(t)dt-(6-2-4)可见，要确定船舶在任意时刻的位置和姿态，首先应该求出在附体坐标系内u,v,r的变化规律，为此需要建立船舶运动的动力学方程。平面运动中船舶各点上速度之间的关系1) 刚体运动分解为移动和转动从运动控制角度将船舶视为刚体是足够准确的，因此其运动是由移动(translation)和转动(rotation)

8、叠加而成；可以取船上任意一点为参考点，船舶一 方面整体地随该参考点平行移动，另一方面绕该参考点同时发生旋转运动；移动速度即参考点的速度，故与参考点选择有关，转动角速度则与参考点无关，即对任意的参考点均为同值，对于船舶平面运动，该转动角速度即为艏摇角速率r。2) 船舶任意点 P处的合速度 取o为参考点(图，船上任一点 P对o点向径为Po =xi - yj, i, j为ox及oy轴上的单位向量。以向量形式表示旋转角速度， 有3 =rk , k为沿 oz轴的单位向量，3即为艏摇角速度向量。由理论力学，因刚体转动而造成的速度为V r =3汉卩°,故P点的合速度是Vp =V +Vr =Vpo

9、=(u _ yr)i 十(v + xr) j(6-2-5)注意：单位向量X乘所得向量满足右手法则，如k i，右手从k的正方向逆时针握向i的正方向，大拇指所指方向即j的正方向，如果方向与 j的正方向相反，结果加负号。图移动与转动速度的合成考虑船舶质心C,其对o点之向径为Pc .pXc.i yc j，则C点之速度为VC =V 3 Pc. 他 一 yc)i (v 哭cJj = ui (v Xc r) j(6-2-6)I、 *上式最后一步是由于船舶配载对称于纵舯音剖面，yc =0。如果取质心 c为参考点，应该从oxy坐标系过渡到c坐标系，后者是前者沿 ox方向平行移动距离 ?c而得。P对c的向径为d

10、= i j，于是有Vp 二Vc 3 d =仏：.r)i (vc r)j(6-2-7)6.3船舶平面运动的动力学在推导船舶运动方程时，做下列假设：-船舶是一个刚体；二大地参照系是惯性参照系；'水动力与频率无关，水的自由表面做刚性壁处理。有了第一个假设就不用考虑每个质量元素之间的相互作用力的影响，而第二假设则可以消除由于地球相对于恒星参照系的运动 所产生的力。(1)平移运动方程的建立1) 刚体的动量刚体被看做无数质量微团的集合体，各微团保持其形状及彼此之间的距离不变。刚体动量 G为各微团动量Vpdm的积分，即上式最后一项按照质心的定义应为零，设m是刚体的总质量，则G = mV c(6-3-

11、1)2) 刚体动量定理牛顿运动定律指明，刚体动量的变化率等于其所受外力之和。以F =Xi Yj代表合外力，其中，X是作用于ox方向上的外力，Y是作用于oy方向上的外力， 有dG /dt = F(6-3-2)利用式(6-2-6)、式(6-3-1)和式(6-3-2)，且注意到di/dt二rj, dj/d-ri (因整个坐标系是 建立在附体坐标系基础上的，而附体坐标系是随着船舶的移动和转动而移动和转动的，故其导数存在。如果在惯性坐标系，则其导数为0)，参见图，经整理得m(u -vr -xCr2) = Xm(v ur xC r) =Y(6-3-3)图6.3.1单位向量微分关系式(6-3-3)即为船舶平

12、移的动力学基本方程，是由于建立船舶运动数学模型及mur是船舶宏观旋转中向丿心惯性力分量；附加项 原点o做旋转运动产生的口O注意其形状与熟知的牛顿方程有所差异，这的oxy是非惯性坐标系所致。 式(6-3-3)左端附加项-mvr-mxC r2及mxc r分别是由于质心 C对 刘生力及切向惯性力(离心惯性力)。(2)旋转运动方程的建立1) 刚体的动量矩刚体对质心C的动量矩H C为各微团对C动量矩d (Vpdm)的积分，即H C = (d VP)dm = (d VC)dm 亠 i(d 3 d)dm =(i j) rk ( i j)dm = r ( 2 2)dmk = rk(6-3-4)其中| = .2

13、 2)dm为船舶对过C点的垂直轴(o )的惯性矩。2) 对质心C的动量矩定理同样由牛顿运动定律，运动着的刚体对质心C的动量矩变化率等于其所受外力矩之和，以Mc =Nc k表示后者，Nc为外力矩之代数和，于是即I T = Nc(6-3-5)3) 对于坐标系oxy原点o的动量矩定理形式为式(6-3-5)的动量矩定理只适用于质心 C。现由该式出发对力矩和动量矩进行变换以导出适用于o点的动量矩定理表达式。以M °二Nk表示外力矩之和，其中 N是作用于船舶的绕 z轴的外力矩，以 表示船舶对oz 轴的惯性矩，由理论力学的力矩和惯性矩移轴公式，有M ° =M c - pc F及 Izz

14、=1m*，这样由式(6-3-4)和式(6-3-5)可推出I zzr mxC (v ur) = N(6-3-6)式(6-3-6)即为船舶转动的动力学基本方程，其形状与式(6-3-5)的区别在于，左端的附加项mxCv及mxCur分别代表由于质心 C对原点o做旋转运动所产生的离心惯性力矩和向心 惯性力矩。6.4船舶平面运动的线性化数学模型综合式(6-3-3)和式(6-3-6)，得下列形式的船舶平面运动基本方程 2m(u vr - xCr ) = Xm(v +ur +xCr) =Y>I zzr + mxC(v + ur) = N j(6-4-1)当附体坐标系原点取在质心 C时，XC =0，可得最

15、简形式的船舶平面运动基本方程m(u - vr) = X | m(V+ ur) = Y »I z=N(6-4-2)式(6-4-1)代表着3种力的平衡关系：左端是船体本身的惯性力和力矩，右端是流体对船体运动的反作用力，实际上包含了流体惯性力和力矩及黏性力和力矩。式(6-4-1)本质是非线性的，其左端显式地出现ur,vr等非线性项，尤其右端的X，丫，N将是运动变量和控制变量的多元非线性函数，结构异常复杂。(1) 船舶平面运动的非线性模型和线性模型船舶运动数学模型分线性化数学模型和非线性数学模型两大类。研究船舶数学模型通常有两种目的：一种目的是建立精密程度不同的船舶运动仿真器(又称船舶运动模

16、拟器)，用于通过仿真对船舶操纵特性进行研究，对船舶运动闭环控制系统进行研究，对船舶运动控制器性能进行评价。这种模型必须是非线性的，以包含尽可能多的机理细节；另一种模型目的是用于船舶运动控制器设计，这种模型主要是线性的，因为迄今为止，线性反馈控制理论仍是能够提供各种控制器设计系统性方法的惟一控制论分支。当引用神经网络控制或模糊控制 时，非线性船舶运动数学模型可以提供训练和学习的数据。1)船舶平面运动非线性数学模型为应用方程式(6-4-1)求解船舶平面运动的基本变量u,v,r，必须具体讨论流体动力X,Y和力矩N的结构形式。研究中把船体、螺旋桨和舵视为一个整体，此时X,Y,N将是移动速度(u,v)、

17、转动角速度(r)、它们的时间导数(u,v,r)、舵 角()以及螺旋桨转速(n)的非线性函数X = X(u,v,r,u,v,rQ, n厂Y =Y(u,v, r,u,v,r,® n)N = N (u,v,r ,u,v, r, 6 n)(6-4-3)完全从理论上确定式(6-4-3)的函数关系极为困难，迫使研究者不得不转向半理论半经验 的方法或多元数据回归方法。Abkowitz提出一种小扰动和 Taylor展开研究X,Y,N的表示式的方法，其主要思路是， 考虑船舶等速直线运动这一平衡状态：u二=V,v = r =0,、： =0,u=v = r = 0，这时在式(6-4-3)中的自变量n将不出

18、现；从该点出发，研究偏离平衡状态不远的运动： u 二 u0：u,v 二 v, r 二 r,u 二：u, v 二 v,r 二 二：，：u, v,，: 是小量； 将X,Y,N在平衡点附近展成 Taylor级数时，在展开式中将仅出现u,v和r的一次项，因为流体对船舶的惯性反作用力只取决于平移加速度u,v以及转动角加速度r本身，而与它们的各阶导数无关；至于和：u,v, r有关的黏性力各项及与 ：有关的舵力各项，则取至 3阶为止，更高阶的项全部略去。将式（6-4-3）的展开式代回式（6-4-1）并进行移项整理，可得到 Abkowitz非线性船舶运动方程。Norrbin发展一种非线性船舶运动数学模型，该模

19、型有两个特点，一是适用于运动变量（u,v,r）的整个变化范围；二是它不像Abkowitz模型那样，完全按数学方式处理流体动力，以至其 Taylor展开式的某些项缺乏物理意义，而是在更深的层次上依赖于流体动力学 的基本原理，构成一种半理论半经验的模型格局。以上所述的Abkowitz和Norrbin船舶运动非线性数学模型属于“整体式”模型，本节 将做较详细的介绍。与此相对应，日本船舶操纵数学模型小组（Manoeuvring Model Group,MMG）提出了一种分离式船舶运动数学模型，后者是在单独考虑船体、螺旋桨、舵的流体动力学特性的基础上再研究在它们构成一个推进和操纵系统时，各部分之间的相互

20、干扰。这种分离式模型的优势是具有完整的理论支持，易于进行实验研究从而获得较为通用的数据回归结果，对于希望建立自己的复杂程度不同的船舶操纵模拟器的各类研究人员均有裨益。有关MMG模型的结构和细节，有兴趣的读者可参考文献9。2）船舶平面运动线性数学模型 成Taylor级数时只保留一阶小量 而得到平面运动线性数学模型，有沿用Abkowitz的研究方案，在把流体动力X,Y,N展开S-Xu)00 画1-0(m -YQ(mxc -Yr)v=0(mxc NJ)(Izz Nr) 一一以矩阵形式表示之，有0YvNv0(Yr - mu。)11 _0+ Y汀(6-4-4(Nr -mxcUo)_ 屮6j，同时在船舶运

21、动基本方程左端也进行线性化处理，从）式（6-4-4）对研究平面运动稳定性有用。3）前进运动与横漂、转首运动的解耦式（6-4-4）表明，在线性化前提下，前进运动与其他两个自由度上的运动互相独立，从航速控制的角度，该自由度的运动可以单独考虑；横漂及转首运动之间存在着强耦合，这两个自由度上的运动与船舶航向、航迹控制密切相关， 是本章研究的重点，故而将式（6-4-4）重新写为（mX u ）u = X u - u(6-4-5)(m-Yv)(mxc -Yr) v = Yv (Yr -mu。)v Y；(mxc -Nv) (Izz-NJ.p. ILNv (Nr -mxcu°)N （6-4-6）4）流

22、体动力导数的无量纲化船舶线性化数学模型的进一步推演主要涉及10个流体动力导数Yv,Yr,Nv,Nr,Yv,Yr；Nv, Nr；Y, N .，前4个称为“速度导数”，第5第8个称为“加速度导数”，最后两个称为“舵力和舵力矩导数”。由于船舶（包括桨、舵）几何形状的复杂性，应用理论流体动力学方法计算这些流体动力导数是不可能的，因此它们的确定必须依赖于船模试验。为了数据处理的科学性以及使用的方便性，根据相似原理和量纲分析方法，应该采用无量纲的流体动力导数。为此选择一些基本的度量单位：长度L。一一 L（船长），速度V V（航速）,时间t。 L/V，质量m。（1/2）兀3 ,力F。（1/2）、2L2，力矩

23、M。质量:长度：Xc = Xc / L速度:v = v/V转首角速度：r =.rLZ力:惯性矩:' 1 2 2F =F/( V2L2)2l；z =lzz/(l 兀5)2力矩：N =N|mL2lzz =布'13m =m/( A.)2以此为基础，将得到无量纲速度导数"丨®黑V；=脅匕=£Yv，N； =_：N'/# =(2/C、VL4)Nr,，以此类推。以上介绍的无量纲化流体动力导数称为撇”系统(prime system)，由美国造船与轮机工程师协会(SNAME)于1950年提出；此后Norrbin又提出了“两撇”(bis system十7，其独

24、到之处是采用与上述不同的基本度量单位，如：长度L,速度 .gL，时间 .L/g，质量八为排水体积，力，力矩 心亠。由此得出的无量纲流体动力导数以Y'，N；,表示。5) 线性流体动力导数的估算公式Clarke整理大量船模试验数据，给出关于10个线性流体动力导数的回归公式10，汇集如下：Y； =-1 0.40CbB/T n(T/L)2(6-4-7)上式中B,Cb,A、分别指船宽、吃水、方形系数、舵叶面积。上式中的YV，Yr',N；,N；是船体尚需对这些流体动力本身的流体动力导数，在实际应用时应考虑舵对船体流体动力的干扰,导数做一定的修正，需修改的增量按下式确定 10(6-4-8)(

25、2) 状态空间型船舶平面运动数学模型状态空间型的船舶运动数学模型是船舶运动控制器设计的基础， 化方案，不同维数的模型用于不同的设计目的和精度要求，详见文献1)二自由度状态空间型船舶线性数学模型它可以有多层次的模型9。1 3在式(6-4-6)的第一行两端除以一L，第二21行除以2 L，并转化成无量纲流体动力导数，则有| m，L(mxc _Y,) IfYvm XC _Nv" L(lzz_Nr") _ rV(Y/-m) 1V ''''-Nv V(Nr mxc)_VV2'一工呵V2(6-4-9)I X =P(2)X(2)- Q(2)U上式可

26、简记为(6-4-10)其中分别是惯性力导数矩阵、黏性力导数矩阵及舵力导数矩阵，X二VJT是状态向量，U二：是控制输入。将式(6-4-10)化成标准的状态空间形式，得X = A(2) X (2) ' B(2)二(6-4-11)其中并且an =(l；z N；)Yv' (mxc Y；)nVV /S)印2 =(l；z _N；)(Y； _m')_(m'xc _Y、(N； _m Xc)LV / Sa21 =-(mxc nV)Y； +(m' Y?)nVV / L/S1a22 =(mxc N；)(Yr'm')+(m'Y)(N； mxc)V/S,

27、j''''''' nb11 =(lzz - Nr)Y&-(mxc - Yr)N§V / $bm =-(mxc -Nv')Y5+(m -Y;)N5V2 / L/S1S1 =(l； _N；)(m' _Y)一(mxc _N；)(m xc -Yr')L丿(6-4-12)2) 三自由度状态空间型船舶线性数学模型在式(6-4-11)的基础上，增加一个便于研究问题的状态变量：'(航向偏差)，且.V-'-'-；r；r为设定航向，使状态向量成为X二v r 八T。因附二r，可得X(3)= A(

28、3)X B©)(6-4-13)其中3阶模型是最基本的，由此可演化成其他更高阶的模型形式，直接利用3阶模型可进行线性二次型(Linear Quadratic, LQ)最优控制器设计。3) 四自由度状态空间型船舶线性数学模型在式(6-4-13)基础上再叠加以舵机伺服系统的模型，后者一般被视为一个1阶惯性环节，其时间常数为Tr，则有-1、r(6-4-14)TrTr其中：r为命令舵角，则状态变量成为X鬥vr 尹T，可得到X =A(4)X B(4) - r(6-4-15)其中4)考虑随机干扰时的线性船舶数学模型考虑海上环境干扰对船舶的影响，并把这种干扰简化为一种白噪声 w(2)= Wi W2

29、r,则船舶运动数学模型将从确定性系统变为随机系统,这样有=A(2) X(2) B(2)丄 W(2)(6-4-16)(6-4-17)X =A(3)X B(3)W(3)W=W2 W3 r.X =A(4) X B r W(4)W( 4)=W2W3 W4 f(6-4-18)其中白噪声w3代表航向角受到的高频噪声，w4代表海浪对舵叶驱动伺服系统的干扰作用。(3) 传递函数型的船舶运动数学模型传递函数型数学模型在经典控制论以至智能控制范畴内用于分析船舶运动的动态行为， 并且可作为设计航向、航迹控制器的基础。1)3阶传递函数模型对于船舶航向控制来说，采用3个自由度的状态空间数学模型式(6-4-12)，加上输

30、出方程-CX (6-4-19)其中m为量测航向，C= 00 1 1,将此状态空间模型转换成传递函数形式为G.(s)二 C sI -A4 B =K°(T3S 1)s(T1s1)(T2s 1)(6-4-20)这是一个3阶系统，具有两个非零极点和一个零点，且有由此不难解得3个时间常数T1,T2,T3以及一个系统增益系数Ko。2)2阶传递函数模型(Nomoto模型)野本(Nomoto)对3阶船舶模型式(6-4-20)做了一项出 色的简化工作，使之降为2阶11。论证的出发点在于，对于船舶这种大惯性的运载工具来说， 其动态特性只在低频段是重要的，故在式 (6-4-20)中令S二j= ； 0，且利

31、用一个熟知的近似关 系：当x > 0时有(x) 1/(1 - x)，并忽略2阶和3阶小量，由此导出著名的 Nomoto模型G.(s)=K。s(T°s 1)(6-4-21)其中增益K0与3阶模型相同，时间常数 T。TT3，或直接由下式求出式(6-4-21)广泛应用于船舶自动舵的控制器设计中。用Nomoto模型进行船舶运动控制器设计有两个好处：一是在低频范围，其频谱与高阶模型的频谱非常相近；二是设计出的控制器阶次低，易于实现。求解本节所述船舶运动数学模型需要已知8个船舶参数，即航速V ,两柱间长L ,船宽B，满载吃水T ,方形系数Cb，排水量I ,重心距中心距离xc，舵叶面积A 。

32、首先将这8个已知参数代入式(6-4-7)，求出10个流体动力导数，并用式(6-4-8)修正， 然后代入式(6-4-12)，即可求出各种自由度的数学模型。6.5船舶平面运动的一种简洁非线性数学模型(1)用于船舶运动闭环控制系统仿真的六自由度非线性模型各种线性船舶数学模型只用于在不同情况下进行控制器设计，当用于船舶闭环控制系统仿真研究时，必须以非线性模型表述被控过程的动态特性，并且还需考虑风、浪、流造成的环境干扰。从式(6-4-10)出发，在其右端加上非线性流体动力项F non、风力项F WIND、浪力项F wave，则无量纲的二自由度非线性船舶运动数学模型将呈下列形式IIIIII1 (2) X

33、=p(2) X(2) ' Q(2) u ' F NON，F WIND F WAVE(6-5-1)其中Ynon , Ywind , YWave 及N NON , NwIND , NwAVE 分别是非线性力、风力、浪力在y方向的合力及在绕z轴方向的合力矩。由式(6-5-1)和式(6-4-14)不难看出在X(4)=vr二才尹T的4个自由度上有非线性状态方 程(I '(2) ) F NON + F WIND + F WAVE X =A(4) X + B(4)6 +0(6-5-2)- 0 一考虑到船舶位置X0, y°T的两个自由度上的运动学关系(6-5-3)X0 二 U

34、COS- -vs inyo =usi n'vcos'-则式(6-5-2)与式(6-5-3)构成了六自由度的船舶运动非线性数学模型的基本框架，状态变量变为 各研究者关于式(6-5-2)中非线性流体动力F non的取法不同是区别到目前为止形形色色的非 线性船舶运动数学模型的主要标志。Norrbin关于非线性力的简化表示式F NON 的简洁表示Norrbin在研究船舶参数辨识问题时提出了一种关于非线性流体动力式712如下所示:'YNONCfY (v，r) I卜 NON'| N NON | Cf N (v，r)(6-5-4)其中fY(v,r)=< T .rr11v

35、2T rr一石 一”)212 L2 r1 1 v 2 1 v 33oC <2 L r 3 L311v2T rr石十p(-)212 L2 r1<21-<21vlr1vlr1vlr1< 21<-2<旳(6-5-5)1 1 v()6 L r11 1 v 2T一启一匚字()232 4 L2 r1 1 vT rr()6 L rT r rfN (v,r) “1 v 414(-)4-6 L4 r2121 v1<-< L r21 v 1< -<L r 21 v<- <=°L r(6-5-6)式(6-5-4)中的比例系数C为无量纲

36、横流系数，其值通常在0.30.8范围内。Norrbin关于Fnon的横流模型式(6-5-5)、式(6-5-6)的优越之处在于其表示式在各类非线性模型中最为简单，它 的导出具有比较明确的理论基础，并且公式中除了船舶吃水和船长之外，不需任何关于船体结构的数据，应用甚为方便。据笔者的经验，由式(6-5-2)式(6-5-6)组成简化的非线性船舶运动数学模型用于在自动舵控制下的闭环系统的仿真研究，结果是可信的1315。应指出，对式(6-5-5)和式(6-5-6)中同时出现v =0和r =0的情况应做专门处理。(3)风力干扰在式(6-5-2)中，风力 YwiND , N WIND 分成平均风力YwiND

37、, N WIND及脉动风力Nx6.5.1平均风力的表示式如下=C<C/R)7pAVRyAL/1PL3、>wind,NWind16。平均风力计算见图YwiNDONWind 二Cn( r)aV；AlL/今九4(6-5-7)式(6-5-7)中Cy( r),Cn( r)为无量纲的风力和风力矩系数，文献17给出了这两个系数的一密度；Vr为相对风速，r为相对风速与首向间的夹角，称为风舷角，由绝对风速Vt、绝对风向:-WIND以及航速V（u,v）按下式计算u _VT COS（C（WIND 即）uVR=-VrSin（o（wiND -屮）一V »Vr =Ur +VrJ（6-5-8）上式中:

38、WIND变动范围为0360 , R变动范围为0一 180，相对风从右舷来时R>0。脉动风力YwinD , N WIND是由大气的湍流所造成的，按文献 16，它们被认为是某种白噪声的实现，该白噪声的标准差与绝对风速Vt的平方成正比二Y = 0.2 ; AVT cy（ r） L oQ（6-5-9）二 N =02，AVt|CN（ R）L（4）浪力干扰浪力Ywave,Nwave分为两个组成部分：高频的一次力，它是与波浪宏观振荡运动同步的周期力，幅值可较船舶的推进力或因运动而产生的流体动力高一个数量级，但由于大惯性船舶本体的滤波作用，一次力产生的振荡运动（艏摇、横荡等）被限制在允许范围之内；低频的

39、二次力，数量级较小，数值变动缓慢，产生船位的漂移。1）一次力的计算采用文献16的结果，把波浪看成规则波，这种波浪只有一个频率 、一个周期T和一个波高h而把船舶看成一个简单的六面体；在小扰动假设下压力由波形抬高按Bernoulli公式求出，浪力是在船体水下表面上把压力积分而得，并表成封闭的解析形式。更准确地可采用不规则波概念，把不同风力下的波谱分解成一系列波谱段（例如10段）,每一段波谱对应着一定的频率和波高，这样不规则波就由一系列规则波叠加而成；船体也被分解成一系列六面体分段（例如20段）；分别计算各种波浪分量在每一分段上的力，最后按 频率和船长进行二维求和可得到总的浪力，但计算量大为增加18

40、，未予采用。规则波对于船的传播方向称为浪向角，以表之，参见图,有二冗一 G WIND -）(6-5-10)=0为顺浪，=二为顶浪,七为横浪（诗表示浪从右舷来）；船对波浪的遭遇频率是e 二 - kuc o skvs 'r n(6-5-11)八 wave propagationdirection 图 6.5.式(6-5-11)中为规则波本身的圆频率,wind direction/ 2 ,4n k =gT>(6-5-12)式中:L时为波长，T为波浪周期，与风速有大，其具体的依存大系视考察的海区而有所不大，同。Kallstrom根据海上观测数据进行最小二乘回归16，得到,(VT)和波高h

41、 ,(VT)公式如下：丁心)=.0014年 +0.042V2 +5.6、於T 2 T>(6-5-13)h (Vt) =0015Vt21.5注意式(6-5-13)只适用于Vt乞20m/s,对于 V20 m/s的情况应谨慎进行外推处理；并力表示式如下：sin b sin cYwave-2aLis(t)b2 .c cosc-sinc 门 . b cosb-sinb、 Nwave =ak(B sinb2- L sinc2) (t)cb且在Vt = 0时仍给出1.5 m的波高和5.6 s的波浪周期,这对大西洋上的情况可能是合适的,但用于我国近海海域，可能稍有误差。Kallstrom还对Zuidwe

42、g的工作19略加修改，给出浪(6-5-14)a - ：g(1 - e *T)/ k2b =kL/2 cos其中(6-5-15)c = kB/2 sinkhs(t)=( 2釣sin©et)式(6-5-14)中(t)代表在附体坐标系原点处波面的振荡，s(t)则表明沿波浪传播方向上的波面的斜率在原点处的值。在进行仿真时应对式(6-5-15) 的- 'e进行适当的滤波16，滤波的方法如下：ef(k) =(AS(k) ef(k-1) (B S(k) e(k)(6-5-16)S(k)A S(k -1)A S(k -1)(6-5-17)式(6-5-16)、式(6-5-17)表明所采用的是一

43、个时变滤波器。-'ef为经过滤波的遭遇频率；'e为未经滤波的遭遇频率； A,B是两个常数，取 A=0.999, B=0.001(即A+B=1); S(0)取为0.999,随着递推次数k的增加，由式(6-5-17)知S(k)下降，当k:时S(k) >0，此时式(6-5-16)趋于 一个定常滤波器，ef (k) = 0.999，ef (k -1)0.001，e(k)(6-5-18)滤波的结果是：在连续的采样周期实也连续变化。2) 二次力的计算目前尚无简捷而可靠的方法。海浪干扰的另一种简单模拟方法是用白噪声驱动一个典型的2阶振荡环节(相当于2阶低通滤波器)20。其中白噪声的带宽

44、为0.5 Hz, 2阶滤波器具有低阻尼比，参数为 0.05,自然频率 其中'0为海浪频谱的峰值频率，U为船速，为航向与海浪方向之间的夹角，g为重力加速度。例如，如果模拟的海况为5级风，中浪，参数可取 0为0.15 rad/s，船速为7 m/s(约14kn)，为 60。(5) 流干扰仿真时通常假定流是恒定并且均匀的，它只改变船舶运动的位置和速度，而不改变船舶的航向，有下列速度平衡方程：x0 =ucos'-vsiVc cos cy0 二usinF 】 vcos 亠Vc sin c其中Vc,Yc分别为流的绝对速度和绝对流向，见图N(6)船舶运动数学模型的总体结构x图6.5.3流(6-

45、5-19)图，还包括了模型和控)制器的连接：该图表述了仿真研究的信号流程。图船舶运动非线性数学模型总体结构.一V八，WIND (风速、风向 )操舵伺服訶服系统的数学模型模型6.6 0操舵伺'在式(6-4舵伺服系统是-14非线性力舵力，把舵机伺月服系统看成1阶惯性环节是一种较粗糙的风力浪力'以。实际上，操具有纯迟延、 船E区、滞运环动饱和等非线模特性的电动液压系统，这些因素在很大程度上影响到航向/航迹闭环控制系统的性能；换言之，制质量，除了要依赖各种“高级的”航附保持、航迹保持控制算伺服系统这一£(航向环和航迹环)之间近来已逐步为人们所时认识,一但是单从自动舵设计者的角

46、度进行努力不可能根要获得良好的航亍向和航迹控 法之外，还需十兮注意操舵的匹配船型这一点虽然 置显示本解决问题, 必须从自动舵与操舵电液伺服系统的结合上进行综合考虑在整个船舶运动控制的层次上， 在设备的选型、安装、管理以及控制方案的确定、 控制算法的设计等诸多方面进行细致的工 作，协调处理，方能收到良好效果。对于操舵伺服系统的分类以及性能比较、操舵伺服系统引起的船舶运动附加阻力等问题，Bianke曾进行过相当深入的研究21,22。按照Bianke的观点，操舵电液伺服系统可概分 为5类，其定义及模型化概述如下：第1类：单油路bang-bang控制伺服阀系统由命令舵角和实际舵角所形成的误差信号（：和

47、-'：）经功率放大，引起三位四通电液伺服阀一侧的电磁线圈（solenoid）通电，打开操舵主油缸（hydralic actuator）的通路，由定排量主油泵来的压力油驱动舵叶回转，直到实际舵角 与命令舵角一致为止。其原理和仿真模型参见图，其中DB为死区宽度，H为滞环宽度，N为最高转舵角速率。对于一艘250 000载重吨油船，典型数据为DB=2 , H=1 , N=2.3/s。solenoid第2类：双油路bang-bang控制液压操舵器系统它角误差信号 r 一 6）经功率放大，引起三位四通电液伺服阀一侧的电磁线圈通电，打开液压操舵器（telemotor, TM）的油路，由定排量辅油泵来

48、的压力油使TM的拉杆移动，这是一级放大；TM拉杆因而拉动变排量主油泵的油量控制杆，使主油泵排出与控制杆移动距离成比例的油量，这个压力油流被通至转舵主油缸，驱动舵叶回转，与此同时由舵柱带动的机械式三点追随机构产生位置反馈，把主油泵的油量控制杆拉回到零油量位置，此时：动态地停留在的位置上，这是二级放大。这类液压操舵系统在海船上应用相当广泛，其动态性能（操舵时间和舵角跟踪精度 ）明显优于第1类电液操舵伺服系统。图，图；K为一级放大系数，N为二级放大系数。典型的数据：DB=1 ，H=0.8，PB=7，K=4 /s，N=2.3 /s。第3类：单油路bang-bang控制主变量泵的主油路本系统实际上是把第

49、 2类系统的三点式追随机去除后形成的，因为没有了舵角位移的二级反馈，所以变量油泵的油量控制拉杆只能有3个位置：左满程、右满程、零位。这和（-、：）信号的符号是一致的；换句话说，此时的变量油泵只作为定量油泵使用。因而此类系统的功能类似于第1类系统，但其性能要优于后者，原因是在 3个位置之间的转换（即TM拉杆的移动）是逐步发生的，因而其动态过 程自然比油路的突然开、断控制的第1类要平滑，舵角超调也要小得多。图，其中比图，H=0.8，PB=7，K=4 /s，N=2.3 /s。辅泵I拉 油量控制杆杆信三位四通电磁阀液压（）v 操舵器a） '系位移主操舵回路柱塞'-HH功率 放大-DBa

50、uxiliaam耐仿真模型反馈位移信号严H1main loop机械式模反馈JPB电气式舵角反馈第3类电液操舵伺服系统仿真模型|DBJ附类系统的特点在于叩对用一种比例伺服阀而不是像类1系统中的那种位式伺服阀;2系统中的形式全同；这样本系统就存在着两个串联的连续运动环节，如图，其动态性能在第4类：双油路模拟控制操TM的油量控制采至于对主油路的控制则与类型各类中是最好的，但初置费明显增加。典型数据为PB1=1，K=4 /s， PB2=7，N=2.3 /s。N：图| 6.6.4+ TN伺服系统仿PB2此系统设计较为简捷,为只用了单油第5类：单油路；由于采用了比例控制的直线位移输出的伺服机构去拉动变量油

51、泵的油量控制杆；使舵叶 的运动快速，并且不会产生舵角的静态误差，其仿真模型见图，N=2.3 /s。u明显下降，此Abkowitz当船舶进行大舵角回旋操纵时，u, v, r之间的非线性耦合效应将使航速时上述的船舶运动数学模型已难于描述过程的动态，需要转向更为精密的模型。 非线性船舶运动数学模型是这方面的一个突出代表。(1) 模型结构形式的讨论综合船舶平面运动基本方程式(6-4-1)和作用于船舶上的流体动力式(6-4-3)，有m(u -vr xc2) =X(u,v,r,u,v,r,、Jm(v +ur +xc)=Y(u,v, r, u, v, r,6)£(6-7-1)Izzr +mxc &

52、amp; + ur) = N(u,v,r,u,v,r*)1) 船体惯性力式(6-7-1)左端的船体惯性力各项不做任何简化，其中与平移加速度u, v及回转角加速度F成正比的项保留于方程左端，其他各耦合相乘项和平方项全部移到方程右 端，与流体动力的有关项合并。2) 流体惯性力流体对船舶运动的惯性作用力只与u,v,r的一次幕有关，且由于船体形状的左右对称性，X展开式中不应存在横向加速度v及转首角加速度r的项；类似地，在Y,N展开式中也不应出现纵向加速度u的项。故有 X :Xuu; Y:Yvv Yrr; N : N Vv Nrr。这些项将被移至式(6-7-1)左端，与船体惯性力合并。3) 流体黏性力

53、这里既包含船体本身与周围介质的相对运动造成的升力和阻力，也包括各种控制面上产生的驱动力， 如舵叶和桨叶上的推力和阻力。 阻力直接与黏性有关； 升力虽 然数值上可用势流理论计算，但是物理上它的起因却是由于黏性造成初始的边界层分离从而产生环流，因而把这一部分流体动力含蓄地统称为黏性力还是适宜的。由于船形(包括置于正中位置的舵和桨)左右对称，X应是v,r,及其交叉乘积的偶函数，因为不管v,r,、：为正或为负，引起的 X分力应是同样的； X作为u的函数，一方面应包含:u2 ：u3类型的项，另一方面还应出现与v,rr构成的乘积的偶函数项。因此在 3阶精度范围内，X的黏性力部分应该包括下列各项，X : X

54、0, Xuu, Xu肿u2,Xuuu'u3; Xwv2, Xrr2, 乂加2公八，X-.V-, X"、：； XuVu, XrJu, X” 2UXvru vr ：u, Xv u. ：u, X r u r '。同样由于船形的左右对称性，Y, N应是yr,、.及其交叉乘积的奇函数，因为 v,r,、.的每一个变动方向时产生的横向流体动力Y及力矩N也要改变方向。同理还应考虑到纵向速度u的作用，一方面出现.：u：u2-：u3类型的项，另一方面应出现 Au与yr,、.之间的奇函数耦 合项。在3阶精度范围内，有233332N : No, Nu.：u, Nuu ：u , Nuuu ：u ,NvV, NjN .、，Nv ,心訂，N/，Nrv ,r22r2r2r2rrN vv "V ,N vrrvr , N “ ：；，N £：¥：；, N r:f , N vr：；.vr 'J ； N vu v- u , N ru u, N u,NvuuVg2, Nruur u2, N uu、U2。Xo,Yo, No是平衡状态下出现的不对称不平衡流体动力，部分是由于螺旋桨单向转动造成的。 注意：在上述关于X,Y,N黏性力的Taylor展开式中，为简便起见，已把各阶次对应的数字系_3数1/2,1/6并入到相应的流体动力导数中去，例