船舶运动水动力线性解二维、二维半与三维方法的比较.pdfVIP

船舶运动水动力线性解二维、二维半 和三维方法的比较 段文洋 马山 (哈尔滨工程大学船舶工程学院，啥尔滨，150001) 摘要本文应用二维时域格林函数方法求解基于高速细长体理论的排水式船舶在波浪 上运动的水动力定解问题。该理论的控制方程和物面条件是二维的，自由面条件是三维的，因 此又被称为二维半理论。应用该理论计算了WilgeyI船模的垂向辐射水动力系数和波浪力， 同三维线性时域有航速方法和基于二维切片理论的STF法的预报结果和模型试验结果进行 了比较。通过比较，可以发现二维半方法给出的结果同三维时域有航速方法给出的结果预报 精度相当。由于二维半方法的计算效率较高，故特别适合于高航速下船舶的酎波性和波浪载 荷预报。 关键词船舶水动力。二维半理论，三维时域有航速理论，STF切片理论，数值结果比较 1 前言 时至今日，有航速特另4是高航速下，船舶在波浪中运动水动力问题，即使是线性解问题，仍 然存在诸多问题。一方面，传统的基于二维水动力解的切片方法虽简单，但由于计及航速效应 甚少，故精度差，特别是对水动力运动模态耦合项和纵摇项预报与试验结果差别较大；另一方 面，较完整的三维问题求解，虽然频域解和时域解理论和方法繁多，但应用尚不普及，这是因为 这类方法普遍费时费工，对某些实际船型还存在算法稳定性和收敛性问题。为了改善有航速 船舶水动力的预报精度，近十年来，提出了称为二维半理论的耐波性实用计算方法。二维半方 法在保留细长体假设前提下，抓住了船舶有航速这一主要特征，近两年来，我们通过发展二维 半流场快速解法和水动力求解稳定算法，逐步使其实用化。 大学完成的标准Wigley系列船模试验数据为标准，采用标准STF切片方法、我们提出的二维 半方法和美国麻省理工学院发展的三维线性时域解方法(TIMIT)进行水动力系数的比较发 现，在中高航速下，二维半理论预报结果优于切片理论结果，和三维方法结果比较，除纵摇阻尼 在高频时有一定差别外，水动力系数均相当甚至更优。而二维半理论和切片法计算时间相当， 因此．在缺乏高效率和可靠三维有航速理论的现状下，我们提出的二维半理论不失为改善耐波 性理论预报精度的实用选择． 2船舶运动水动力线性解理论简介 船舶在波浪中沿固定方向恒速前进，由于受到波浪的作用，船舶除了以恒速U前进外，还 有六个自由度的摇摆运动。假设六个自由度的运动都是微幅的，作为一阶无穷小，我们将在线 性假设下讨论船舶受力与摇荡问题。 建立随物体平动的直角坐标系。一a》z，D一掣平面与静水面重合，∞轴与u同向，∞轴 垂直向上。船舶在波浪中的摇荡运动属于刚体6自由度运动，可由船舶重心G点的3个线位移 (纵荡、横荡、垂荡)和绕G点的3个角位移(横摇、纵摇、首摇)表示。在船舶的摇荡作用下，船 体表面上任一点的位移可以表示为 孑(￡)=季(￡)+百(￡)×≠ {(￡)表示为坐标原点处船舶运动的线位移矢量，即 {(￡)；(rh，珑，rid (2) 5(￡)表示重心处船舶运动的转动位移矢量，即 口(￡)=(仉，r／s，琉) (3) 对于上述在波浪中运动的船舶，按照线性势流理论，其周围的流场总速度势为垂(z，，，z r)。由定常势与非定常势两部分组成如下： (4) 垂(z，Y，z，￡)=一[一Ux+熟(z，Y，z)]+4叶(z，y，z，￡) 定常部分是相应于船舶在静水中前进时的流场的稳态定常势。它的量级是O(U)或O(1)。 秭＆，y，z，f)是非定常部分流场势。假设波浪是微幅的，其波幅为一阶无穷小，船的响应也是 一阶无穷小，对应的秭o，y，=，t)也是一阶无穷小。在线性势流理论假设下，听o，y，z，f)可 写成人射势办、绕射势如和辐射势“三部分。绕射势和辐射势统称为流场扰动势如。 6