1、一、物理背景无论是船舶还是海洋平台在海洋开发中都起着关键的作用，而开发海洋首先需要对海洋结构物进行深入地研究。这其中，水动力学中的附加质量是研究的重要方面，掌握物体附加质量的计算无疑具有重要的意义。附加惯性力的存在使物体在理想流体中的变速运动相当于物体自身质量上增加了一个附加质量而在真空中运动，换句话说，理想流体增大了物体的惯性，使物体很难加速也难减速。计算机是求解附加质量的重要工具，本课程设计主要依据分布源模型的面元法等知识来对圆球、椭球、圆柱、双椭球的附加质量进行数值模拟计算，并进行相关讨论。二、理论依据用 表示无界流中的物体表面，来流为均匀流，其未扰动速度或无穷远处s的速度为(2.1.1
2、),1xyzVijVk用 表示定常速度势，它在物体外部空间域中适合拉普拉斯方程，x,yz在物面上适合不可进入条件，在无穷远处，应该与均匀来流的速度势吻合，即（物体外） (2.1.2)20（物面 上） (2.1.3)ns（无穷远处）xyzV其中,单位法线向量 指向物体内部。在速度势 中分出已知的均匀来流项，记(2.1.4)xyzV这里的 是扰动速度势， 应适合以下定解条件：（物体外） (2.1.5)20（物面 上） (2.1.6)Vns（无穷远处） (2.1.7)易知过物面 的通量为零，即ssdn0所以远方条件(2.1.7)可进一步具体化为（ ） (2.1.8)21rO22xyz用 表示点 和
3、之间的距离，对函数 和 在物面 外部和远方控pqrqq1/pqrs制面 的内部之空间域内用格林公式，当点 在上述空间域内时c(2.1.9)14 qpqqpsc dsrnnr从 的远方条件(2.1.8)可知， 上积分趋于零，式(2.1.9)成为(2.1.10)114 qpqqps dsrnnr其中, 是物面 外的任意一点。p在物体的内部域中构造一个合适的内部解 ，它在 内部适合拉普拉斯方is程，在物面 上适合某种物面条件，其具体形式将在下面给出。对于上述物体s外部的点 函数 在物体内部域中没有奇点，在内部域中对函数 和p1/pqr ()iq用格林公式，得到1/qr(2.1.11)10ii qpq
4、 qps dsrnnr式(2.1.10)和(2.1.11)中的 是物体外部同一个点，把两式相减，得到(2.1.12)qs pqiqipq dsrnnr 114在物面 上取 适合下述两种物面条件，得到两种 的定解条件，一种是：si i(2.1.13)20(is内 部 )上 ）定解问题(2.1.13)是拉普拉斯方程的第一类边值问题，它的解是存在且唯一的。取式(2.1.12)中的内部解 为式(2.1.13) 所决定的函数，则式(2.1.12)成为i(2.1.14)qspdr其中(2.1.15)nqi41式(2.1.14)表示扰动速度势 可以用物面 上的分布源表示，其中分布源密s度 是未知函数，将由扰
5、动势的物面条件（2.1.6）来决定。物体的附加质量 ，表示物体沿 方向运动引起的 方向的附加质量，公ijmij式如下:(2.2.1)(,12,6)jjiijisbsbnddsijn利用式(2.1.14)，再结合物面条件，得到(2.3.1)12qsppspqdVnnr这就是分布源密度 所适合的线性积分方程。把积分方程(2.3.1)转换成线性代数方程组，即用离散量代替连续变量。把物面 分成 小块，记sN(2.3.2)1Njjs用平面四边形或三角形来近似代替小曲面 。具体做法如下，取第 小块jsj的四个顶点坐标之算术平均值，得到中心点 的坐标。计算对角线向量的向量jp积（指向与曲面法线指向相符合）
6、，用 表示该方向上的单位向量，形成以jn为法线且通过中心点 的平面，再把四个顶点向该平面作投影，以四个投影jnjp点为顶点组成平面四边形 ，用 代替原来的小曲面 ，称 为单元。jQj jsjQ通常把小范围内的分布源密度 作为常数，因此只要分割不太粗，可以认为在单元 上为常数，记作 ，从而jj(2.3.3)qpQjqps dsrndsrnqj 11因此物面 上的积分可以用 个平面四边形（三角形）上积分之和来近似，N即(2.3.4) qpQjqps dsrndsrnq111上式左端的未知量 是连续型变量，而上式右端的未知量是 个离散量() N。为了求解这 个未知数，须要 个方程。取积分方程(2.3
7、.1)中(1)jN:NN的动点 为 个单元 的中心点 ，称之为控制点，即控制物面pjQ(1)jp:条件使之成立的点。用近似式(2.3.4)代替积分方程(2.3.1)的左端，便可以写出的 阶线性代数方程组：j(2.3.5)1(1,2)NijiabN其中()iijij qpQdsjinr2,iijipabV称 为影响系数，即第 个单元上的分布源在第 个控制点上的影响。ija求解线性代数方程组(2.3.5)得到 的值以后，便可以得到速度势 在控jp制点 处的值，即ip(2.3.6)jiNjicp1(2.3.7)ijijqpQdsr另外，物面 上的诱导速度为s(2.3.8)12Nijiipnvs其中
8、表示求和是不计 这一项。 ，这里的曲面法线 指向物体内部。j n三、数值模型将物体表面划分成四边形面面元，物面为 ，每一个四边形面面元为 。SjS为了简化计算，将面网格投影到各自对应的平面上，使曲面网格 变为平面j网格 。投影的方法为：jQ取 四个顶点坐标之平均值，作为中心点 的坐标。计算 对角线连S jpj线向量的向量积并使得积的方向与流域法向相同。用 表示该方向上的单位向jn量。设： (,)jjxyjzn的坐标为：jp(,)jxyjzp则取投影面为过 并以 为法向量的平面：jpjn()()()0jxjxjyjyjzjzpnp设在该平面上的投影点为： ,z而曲面四边形某个顶点为： 0()则有
9、： 0(,)/,jxyjzxyzn因此得到由 顶点坐标求解投影点（ 顶点）坐标的线性方程组：jSjQ00jxjyjzjxjyjzyxjzjxjzjpnpn由此线性方程组可解出投影点坐标。假设速度势和分布源在 上是不变的，其值为该单元中点（控制点）处jQ的速度势或分布源。由于所求速度势和速度等物理量均为物面上的物理量，因此要令 点落在p物面之上。式右端分布源的法向导数极限由两部分组成，一部分是 P 点附近小曲面 的贡献，另一部分是屋面其余部分 贡献。当 所趋近于的物面上的 sp点 作为控制点的单元，积分时需要考虑奇异性；其余部分为 。设其中一p s单元为单元 ，其余模型为单元 。对于每一个控制点
10、 ，令 循环一次ij i1jN求得前述方程的积分项（包括奇异积分） 。再由 可以得到 组方程，进1而形成求解各个控制点处物理量的矩阵。点表示控制点（编号 ） ，对于每一个控制点的物理量，通过在 和pi 上积分得到。s即：对于任意一 , ，有i qpQjNqps dsrndsrnq1111(,2)Nijiab其中， iijij qpQdsjinr ,iiipabVn将模型控制点数据导入到程序中计 ij和 ib，可以得到方程组：11NNNba 求解上面线性代数方程组得到 的值以后，便可以得到速度势 在控j p制点 处的值，即ip其中，jiNjicp11ijijqpQcdsr每一个 i 点（控制点）
11、处诱导速度为（是向量）： 11()()Nijiiipvps其中， (d)()1lim)jjpiqQjipiqsijrvij在常分布单元假设条件下： li(d)2jpiqiQsnr可得， 11()2()Nijijiiijnvps于是便可求解出速度势和物面速度。四、几何模型4.1 椭球、圆球、有限长圆柱、平行椭球利用计算机编程来完成以上四类几何模型的建立并划分网格如图 4.1、图4.2、图 4.3、图 4.4 所示。其中，椭球长短轴之比为 5：1，有限长圆柱柱体长和截面直径之比为 5，平行椭球的两个椭球相同且长轴平行，间距为短半轴的3、5、7 倍。图 4.1 椭球面网格划分图 4.2 圆球面网格划
12、分图 4.3 有限长圆柱面网格划分图 4.4 平行椭球面网格划分五、计算参数及结果讨论5.1 椭球5.1.1 椭球附加质量系数 M11经查阅资料可知，椭球附加质量系数 M11 的理论值为 0.059。不同面元数时，M11 的计算结果见表 5.1。表 5.1 椭球的附加质量系数 M11 随面元数的变化面元数 椭球的附加质量系数 M11300 0.06204400 0.06138500 0.06117600 0.06101700 0.06090800 0.06074900 0.06064M11 随面元数变化曲线见图 5.1。图 5.1 M11 随面元数变化曲线5.1.2 椭球附加质量系数 M33椭
13、球附加质量系数 M33 的理论值为 0.894。不同面元数时，M33 的计算结果见表 5.2。表 5.2 椭球的附加质量系数 M33 随面元数的变化面元数 椭球的附加质量系数 M33300 0.9706541400 0.9676541500 0.9582145600 0.9518741700 0.9453254800 0.9422548900 0.9406142M33 随面元数变化曲线见图 5.2。 图 5.2 M33 随面元数变化曲线5.1.3 结果讨论由计算结果可以看出，椭球在附加质量系数 M11 较 M33 小很多，这点由椭球几何形状可以容易看出，与实际有较好的符合。随着面元数量的增加附
14、加质量系数 M11 和 M33 均更接近于各自理论值，误差都有减小。在计算过程中发现，面元的划分形式存在质量好坏的区别，对于同样多数量的面元，面元形式的不同会造成不一样的结果，比如同是数量为 600 的单元，长轴方向划分为 20 份，周向 30 份得出的 M33 为 0.964419，而长轴方向划分为30 份，周向为 20 份，得出的附加质量 M33 为 0.9514185，更接近与理论值，误差更小，这是由于对于本细长的椭球模型而言，周方向划分 20 份已经较为紧密，而细长的长轴方向需要划分更多的数量以使面元更接近于真实物面，所以在划分面元时，需要考虑物体的实际情况，以得到更高质量的面元。5.
15、2 圆球5.2.1 圆球附加质量系数 M11经查资料可知，圆球附加质量系数 M11 的理论值为 0.5。不同面元数时，M11 的计算结果见表 5.3。表 5.3 球的附加质量系数 M11 随面元数的变化面元数 球的附加质量系数 M11324 0.5637298360 0.5467839432 0.5368716540 0.5266547648 0.5238728784 0.5142735M11 随面元数变化曲线见图 5.3。图 5.3 M11 随面元数变化曲线5.2.2 圆球附加质量系数 M33由于对称性，圆球附加质量系数 M33 的理论值也为 0.5。不同面元数时，M33 的计算结果见表 5
16、.4。表 5.4 球的附加质量系数 M33 随面元数的变化面元数 球的附加质量系数 M33324 0.5385478360 0.5363544432 0.5337452540 0.5314524648 0.5268541784 0.5234581M33 随面元数变化曲线见图 5.4。图 5.4 M33 随面元数变化曲线5.2.3 结果讨论由计算结果可以看出，随着圆球面元数的增加，计算而得的附加质量系数逐渐接近解析解，即理论值，误差逐渐减小，由面元法知识可知，这是由于面元数量增加，面元密度增加，单元也就更接近真实物面，从而结果更准确。M11 与 M33 计算结果基本一致，差别较小。5.3 有限长
17、圆柱5.3.1 圆柱附加质量系数 M11不同面元数时，M11 的计算结果见下表。表 5.5 有限长圆柱的附加质量系数 M11 随面元数的变化面元数量 m11200 0.137835240 0.139154300 0.140278360 0.140923400 0.141215480 0.141617600 0.141974720 0.142188800 0.142287900 0.142383有限长圆柱 M11 随面元数变化曲线图 5.5图 5.5 M11 随面元数变化曲线5.3.2 圆柱附加质量系数 M33不同面元数时，圆柱 M33 的计算结果见下表。表 5.6 有限长圆柱的附加质量系数 M
18、11 随面元数的变化面元数量 m33200 1.021685240 0.994427300 0.968839360 0.952627400 0.944758480 0.933241600 0.92205720 0.914769800 0.911179900 0.907623有限长圆柱 M33 随面元数变化曲线图 5.6图 5.6 M33 随面元数变化曲线5.3.3 结果讨论由计算结果可以看出，对于柱体长与直径长之比为 5：1 的有限长圆柱体而言，其沿柱体方向附加质量系数 M11 要小于垂直柱体方向附加质量系数M33，并且随着面元数量增加而增大并且有收敛趋势，而垂直于柱体长轴方向的附加质量系数
19、M33 随着面元数量增加而减小并且也呈现出收敛趋势。5.4 平行 椭球5.4.1 两椭球平行时附加质量系数 M11每个椭球面元划分数量为 400，并且划分形式与单个椭球时一样，以便比较二者的区别。单个椭球在以上计算中得到的附加质量系数 M11 的值为0.06149，以该值作为计算两椭球并行时相互影响的参考值。在不同轴间距时时，M11 的计算结果与对比如下。表 5.7 两椭球并行时 M11 随长轴间距的变化间距(倍数) 并行时附加质量系数 M11 影响3 0.07562 24.08%5 0.06754 10.11%7 0.06459 5.10%M11 随长轴间距变化曲线见下图。图 5.7 两椭球
20、并行时 M11 随长轴间距的变化曲线5.4.2 两椭球平行时附加质量系数 M33单个椭球在以上计算中得到的附加质量系数 M33 的值为 0.9675557，以该值作为计算两椭球并行时相互影响的参考值。不同轴间距下的算结果见下表。表 5.8 两椭球平行时 M33 随长轴间距变化间距 并行时附加质量系数 M33 影响3 1.115148 14.90%5 1.004621 3.87%7 0.983651 1.59%双椭球并行时 M33 随长轴间距变化时的曲线如下，图 5.8 双椭球并行时 M33 随长轴间距变化曲线5.4.3 结果讨论由计算结果可以看出，两个椭球并行时，二者之间会产生干扰，使单个椭球
21、附加质量系数增加，然而其影响会随着椭球之间轴间距增大而减小，以使每个椭球附加质量系数逐渐接近于单个椭球运动时的值，可以推测当其间距足够大时，每个椭球可以看成单独运动的效果。参考文献1. 戴遗山,段文洋. 船舶在波浪中运动的势流理论M. 国防工业出版社，2008.2. 陈科. FORTRAN完全自学手册M. 机械工业出版社 . 2009.3. 孙寒冰,邹劲, 庄家园,吴恭兴.深 V 型无人艇附加质量计算及其影响D.2011 年中国智能自动化学术会议论文集（第一分册）.2011.4. 李刚,段文洋, 郭志彬.复杂构型潜器附加质量的研究D.哈尔滨工业大学学报.2010.5. 林超友.潜艇近海底航行附加质量数值计算 D.SHIP ENGINEERING,2003理想流体力学课程设计 应用 HessSmith 方法求附加质量
