近水面单体双体升沉运动有限元研究与分析.pdfVIP

近水面单体双体升沉运动有限元研究与分析 李宗民 邢殿录 16023) (集美大学轮机工程学院福建厦门361021)(大连理工大学船舶工程系辽宁大连1 摘要本文采用有限元方法对近水单体和双体在做升沉运动时的附加质量系数和阻 尼系数进行计算。用Fortran语言编制程序，求解速度势，得到此种运动状态下的附加质量 系数和阻尼系数的数值解。文中分别对二维圆柱和方柱进行计算，并对计算结果进行分析。 关键词有限元附加质量阻尼 中图分类号：U661．1 1 引言 人们研究船舶水动力性能、计算流体动力系数目的在于了解和掌握船舶运动规律和影响因 素，以便在设计初期为耐波性的优化设计提供信息、寻求有效的方法，改善船舶在波浪中的航 行|生能，避免或减轻摇荡运动引起的灾害，使船舶更安全地航行。 计算流体动力系数的方法通常有三种，即源汇分布法、多极展开法和有限单元法。1965年 力学中的有限元方法起点。有限单元法在流场计算中比上述前两种方法具有更加广泛的适用性， 在海洋工程领域有着越来越广泛地应用。在使用有限元法求解流场的偏微分方程时，通常使用变 分法和伽辽金法。变分法中，要用到泛函。但是带有各种各样边界条件的泛函是不易求取的，甚 至不存在泛函。因此在实际应用时，伽辽金法是很方便的。对于泛函存在的问题，用伽辽金法和 用变分法可获得同样的结果，但用伽辽金法可免去求取泛函的困难。本文计算采用伽辽金法，取 插值函数作为权函数，对拉普拉斯方程在全域内进行积分，形成矩阵方程式，具体积分时采用 Gauss数值积分，选取三个高斯点。 2流体动力系数的计算 由波浪的振荡运动产生的扰动势力可写作 6 办(x，Y，z，f)=≯，+≯。+∑白≯／ (2．1) 』=l 其中矽，为人射势，九为绕射势，办(j=1，2，…，6)为物体作第J相单位振幅振荡运动 的辐射势。确定船体上的流体动力，首先要在下列条件下确定矽J(／=1,2，．．…．，6) 119 d矽． 在船体表面平均位置上 (2．2) 1’考2‘绷J+‰／ 其中门J为广义法线，h=(，2l，，22，以3)为外法线方向，(，24，胛5，聆6)=尹×历 f 0 当j=1,2，3，4 朋／={n3 J=5 --n 2 J=6 在z=0上 (2．3) 2)(i圳杀)2办+g知=。 将班分为两部分，吮=秽+兰∥，则 000． a矽iU ． ——=109n：—二=i国in： (2．4) ’ ’ On 加 硝是与船速无关的部分，由此，只要求出硝便可求出死。 咖1=虻 J=1,2，3，4 矽，=心+ 鳄 (2．5) u五u 矽。=簖一 艿 i 在计算出二维流场的速度势后，在物体表面上用线性化伯努力方程，可得到剖面的附加质量 系数和阻尼系数。 咿翱嚣弛抛sNjdt 刨6、 I