浮体水动力分析的基本理论.docx

2 浮体水动力分析的基本理论 2.1 势流理论 流场中速度场是标量函数（即速度势）梯度的流称为势流（Potential Flow）。特点是无旋、无黏、不可压缩。 简谐传播的波浪中具有浮动刚体的流场速度势可以分为三个部分： ?x,y,z,t ?r为浮体运动产生的辐射势；波浪未经浮体扰动的入射势表示为?ω；?d为波浪绕射势 需要满足的边界条件有： 普拉斯方程（Laplace Equation）： ?2 底边界条件： ???z 由表面条件： ?2 没物体表面条件： ???n ⑤辐射条件： 辐射波无穷远处速度势趋近于0 limR→∞?=0 （ 2.1.1 波浪力的组成 浮体浸入水中受到的力和力矩分别为： （2-7） （2-8） S表示浮体湿表面，n的方向是由浮体内指向流场。 用线性化的伯努利方程以速度势表达压力： （2-9） 则 （2-10） （2-11） 辐射载荷表达为、，是由浮体强迫振动产生的；浮体固定时，入射波浪产生的载荷表示为、；浮体固定时，产生的绕射波载荷表示为、；静水力载荷表示为、。 2.1.2 附加质量与辐射阻尼 当浮体发生强迫振动时，其在j方向和k方向产生的耦合水动力包含附加质量和辐射阻尼两个部分： （2-12） （2-13） 如图2.1所示为波激力、附连质量力、阻尼力和回复力的叠加。 图2.1 波激力、附连质量力、阻尼力和回复力的叠加 2.1.3 格林第二公式 应用格林第二公式，两个单独的速度势关系可以表达为： （2-14） S’为封闭体积V’的封闭表面。体积V’由一个假定的、直径为R的圆形范围、深度为z=-h的海底平面以及浮体湿表面积包围而成。应用边界条件，最终可以得出如下结论： Mkj=M 这一结论说明对于六自由度运动的浮体，附加质量和辐射阻尼的6*6矩阵均为对称矩阵。 如图2.2所示为边界条件。 图2.2 边界条件 2.1.4 哈斯金德关系 浮体受到的波浪力/力矩可以表达为： （2-16） Φk为K方向的辐射势，Φd为绕射势是待求解的内容。在零航速条件下应用边界条件，根据格林第二定律可以给出辐射势与绕射势之间的关系： (2-17) (2-18) 则式（2-16）可变为由辐射势和入射势求解的方程： (2-19) 对于零航速的水动力求解问题，波浪激励可以由入射波和辐射波表达。在某些可以进行水动力计算的软件中（如WAMIT），求解辐射-绕射势得出的结果可以与应用哈斯金德关系求出的结果进行对比，这二者的结果应是一致的。 2.1.5 切片理论 切片理论是一种水动力问题求解的近似方法。对于长宽比较大（L/B≥3）、具有航速或零航速的船舶，在计算船体水动力时，可以假定船体由许多横剖面薄片组成，每片都认为是无限长柱体的一个横剖面，最终将三维水动力问题转换为二维水动力问题求解。通过计算船体每个典型剖面的水动力系数，沿着船长积分最终求出整体的附加质量、辐射阻尼和波浪力。如图2.3所示是船体的一个“切片”示意图。 切片理论在船舶水动力分析理论发展中占有非常重要的地位，主要包括三种基本方法： 厄塞尔法（Ursell method）：求解圆柱截面水动力问题。 保角变换法（Conformal mapping）：包括李维斯保角变换和渐进保角变换等，求解近似船体截面形状的水动力问题。 塔赛法（Tasai Theory）：结合厄塞尔法与保角变换的二维切片法。 弗兰克汇源法（Frank Theory）：用于船型截面水动力求解。 图2.3 船体的一个“切片”示意 2.1.6 面元法 当结构物比较大时，规则波下的载荷和运动响应可以用面元法来分析，这在研究过程中比较常见。一般来说，网格的单元质量会影响计算精度，且这种影响不可忽略，因此，一般遵循以下原则来保护计算精度： 面元不能过大，最大不能超过计算波长的1/7。 结构湿表面的几何尺度变化要能够被充分表征出来，例如：圆柱结构的单元布置应至少围绕圆周布置15个，这样才能充分捕捉几何尺度的变化；对于一些不规则的部分，应尽量降低网格的尺寸。 应尽量保证单元分布均匀，距离不易过近。 可以通过多次试算来保证精度。 对于面元法而言，单元之间可以不连续（如图2.4所示）。 应用三维势流源汇法的计算，有时候会出现明显的前后数据趋势的不一致，这属于正常现象，因为面元法模拟的