船舶螺旋桨水动力性能的快速预报方法.docxVIP

船舶螺旋桨水动力性能的快速预报方法 正确预测海洋螺旋桨的水动力性能一直是研究的主题。到目前为止，已经研究了许多理论和方法来计算海洋螺旋桨的水动力性能。例如，螺旋提升力线理论、升力面理论和面元法。升力线理论包含许多假设，仅适用于轻度载荷螺钉。升力面理论可以更精确地预测螺旋桨的动力系数和旋转系数，但由于它只满足于螺丝刀的参考表面，因此预测值的桨叶面的压力分布不正确。面元法没有假设杠杆的几何形状，而是在实际的杆表面上满足物面的边界条件。面元法不仅可以准确地预测螺旋桨的外力和旋转，还可以正确地预测螺旋桨表面的压力分布。 面元法有基于速度势的基本公式和基于扰动速度的基本公式.采用不同的基本公式,不同形状的面元,都可得到不同的数值预报方法.而对Kutta条件的处理和面元的影响函数的计算方法,都将影响数值预报方法的准确性.本文采用了计算较为简捷的基于扰动速度势的基本公式及双曲面形状的面元,在桨叶随边满足更趋合理的压力Kutta条件,并用Morino导出的解析公式计算面元的影响系数,形成了一套快速有效的数值预报方法. 1 nqnq的边界面 考虑在速度为V0的无旋、非粘性、不可压缩来流中的任意升力体,取一外部控制面将其封闭在内.如图1所示,流域的边界面由物面SB,尾涡面SW和外边界面S∞组成.在该流域中可由扰动速度势?来表示升力体的扰动.?满足Laplace方程 ?2?=0(1)?2?=0(1) 根据Green定理,任意一场点P(x,y,z)的扰动势可由以下积分表示 式中 E={01/21Ρ在S之内Ρ在S上Ρ在S之外E=?????01/21P在S之内P在S上P在S之外 R(P,Q)是场点P和边界上点Q(x0,y0,z0)之间的距离,??nQ??nQ是在Q点边界面S的法向导数. 在边界面S的每一部分上,还应满足如下边界条件: 1)当外控制面距升力体极远时,其上的扰动速度趋于零,即 ??→0?当S∞→∞(3)??→0?当S∞→∞(3) 2)在物面上满足法向速度是零的运动边界条件 ???nQ=-V0?nQ在SB上(4)???nQ=?V0?nQ在SB上(4) 式中,nQ是边界面上的单位法向量. 3)假设尾涡面的厚度为零并通过它没有法向速度跳跃和压力跳跃,即 p+-p-=0(???nQ1)+-(???nQ1)-=0在SW上(5)p+?p?=0(???nQ1)+?(???nQ1)?=0在SW上(5) 式中,Q1是尾涡面上的点;上标“+”和“-”分别表示在尾涡面上、下表面的值. 引用以上3个条件,在物面上积分方程(2)可展开为 式中,Δ?为通过尾涡面的速度势跳跃,可记为 Δ?=?+-?-在SW上(7)Δ?=?+???在SW上(7) 对于非定常问题,由于物面的速度势是时间的函数,速度势跳跃也是连续变化,而对定常问题它沿流线是常量. 对于升力问题,必须在升力体的尾缘满足Kutta条件以确定环量强度.本文应用压力Kutta条件求取速度势跳跃Δ?,该条件要求在升力体的尾缘处上下表面的压力差为零,即 (Δp)ΤE=p+ΤE-p-ΤE=0(8)(Δp)TE=p+TE?p?TE=0(8) 积分方程(6)中的3项积分可以分别被看成是在物体表面上强度为-?(Q)的法向偶极子分布,在尾涡面上强度为-Δ?(Q1)的法向偶极子分布和生物体表面上强度为-V0·nQ)的源分布.这是一个第二类Fredholm积分方程,其未知量是物体表面上的偶极子强度,也就是物体表面上强度势的值.结合压力Kutta条件(8),可以求得该积分方程的数值解. 2 柱坐标系和几何螺距角 对于在速度为Va的来流中以角速度ω转动螺旋浆,采用如图2所示的固定于浆叶的坐标系o-xyz.原点o固定于桨盘面的中心点,x轴沿螺旋桨的旋转轴指向下游,y轴沿螺旋浆某一叶片的母线,z轴与x轴和y轴组成右手坐标系.同时也采用一柱坐标系o-xrθ作为参考坐标系.其θ角度量由y轴起始逆时针旋转的角度.这样螺旋桨的相对进流速度为 在柱坐标系下,如图3所示螺旋桨叶面上点的坐标可表为 o-xyz坐标系下的相应坐标为 {x=xr+(-c1+s)?sinβ-(ybyf)?cosβy=r?cosβz=r?sinβ(11)???????????x=xr+(?c1+s)?sinβ?(ybyf)?cosβy=r?cosβz=r?sinβ(11) 式中,s为叶剖面上点至导边的弦向距离;c1为叶剖面上导边至母线的距离;sr为剖面的纵倾;θs为剖面的侧斜角;β为几何螺距角;yb,yf分别为叶背、叶面上的点到弦线的距离. 3 叶片面元数的积分方程转化 为获得积分方程(6)的数值解,将螺旋桨表面和尾涡面分成一系列小单元,并用双曲面元替代每一单元.这里在弦向和展向分别采用如下余弦分割 si=12(1-cosβci)i=1?2???Νc+1(12)rj=12(R+R