三浮式平台总体性能.pptVIP

基于三维分析的实用数值计算手段是存在的，他们能用来预估零Froude数时大体积结构物的线性波浪诱导运动和载荷。 面元法是分析大体积结构物在规则波中线性稳态响应的一种常用方法。 面元法有多种，其中一种方法是在平均湿表面上分布源（和汇），另一种方法是在物体的平均湿表面上混合分布源、汇和偶极子。面元法也称边界元法。 对于浮体的水动力载荷分析，主要使用奇点势函数（单独源或源和正规偶极子）分布船体边界形成诱导流场势。 这些奇性函数（格林函数）满足了拉普拉斯方程和所有边界条件（底部，自由面，无限远辐射），但船体上的不可穿透条件除外。不可穿透条件化为需数值求解的积分方程。 长久以来，此方程的求解遇到了困难：体现在格林函数计算和有效求解线性方程组的数值算法。现在，绕- 辐射水动力分析问题可在计算机上进行求解，变为工程工具。 大部分工业模型建立在积分方程上：将船体离散为许多小面元，在每个面元上源密度为常量。最广泛使用的软件是麻省理工学院最初开发的WAMIT。法国船级社开发的HYDROSTAR软件等。 　　二维极坐标下，拉普拉斯算子可写为： 二维点源（点汇） 　三维点源（点汇） 如果用连续的源系来覆盖一个表面，流体中的速度除了在物体的尖角处以外将会处处为有限。 首先以一个简单的问题来解释源的方法。考虑无限流体中的一个圆柱体，希望求出它的垂荡附连质量。 可以在物体表面分布源来求得速度势。于是有: (1) 将物体表面近似表达成N段直线。如图3.19所示，分成16段。 (2) 假设在每段上源密度为常数。这意味着将式(3.101)近似为所有分段的和，即 (3) 在每段的中点 上满足物面边界条件。首先将源密度标准化。标准化的目的是将时间相关的项与未知的垂荡运动分离开来。记作 积分符号中括号内的项为 。当对 微分的时候，应该认为 与 为变量。这表明 。 这里 和 分别为物体表面法向量 的 和 分量。在圆柱体的情况下 ， 。 式(3.104)表明我们必须求解一个方程组，其形式如下： 这表明如果 单元源强密度为1，它在单元i中点的诱导速度等于 。 此点的总法向速度为每个单元在此点的诱导速度的和，记为 。为了满足边界条件，此点的总法向速度应该等于前面所述的i 单元的速度 。 (5) 确定附连质量的压强部分可由下式计算： 式(3.102)和(3.104)的积分可以有解析解。通过研究在y轴上0与1之间的分段上源强分布的影响，可以展示怎样做到这一点。 本节将展示如何应用源汇法分析大体积结构物上的线性波浪诱导运动和载荷。以一条船作为例子，首先展示如何求得垂荡的附连质量和阻尼。假设船没有前进速度。速度势由下列方程决定： 辐射条件(式(3.117))没有写成数学条件，但在物理上必须保证波浪从船身向外传播。深水中三维辐射波系的一个例子是 可以证明边值问题的解能用满布船体平均湿表面的源来描述。然而，源速度势与无限流体原速度势即式(3.97)不同。速度势应该修正，使式 (3.114)、(3.116) 和(3.117)都能满足。修正显然也要满足Laplace方程。源密度的强度由满足船体边界条件 (3.115) 求得。Havelock(1942，1955)曾证明水深无限时源的速度势可记为下式的实部： 当r较大时，可以证明括号内第三项与最后两项相比可忽略不计。这就意味着已经证明源速度势在离点源远处描述的是辐射波，也可证明自由液面条件（见式(3.114)）能够满足。图3.21显示由源生成的波系图像。 强迫垂荡速度势的解可以用满布船体平均湿表面的源来表示: 对边界积分方程的离散化引出了与无限流场中圆柱的强迫垂荡问题相似的一组未知源密度的线性方程。要指出的不同点是： （1）在波浪问题中源的密度是复数。源密度与垂荡速度一般不同相也不是180°异相。这意味着如果用N 个单元来近似船体表面将产生N个未知的实部和N个未知的虚部，或者说N个复数未知量。 （2) 在波浪问题的数值计算中，源的表达式要比无限流场问题复杂得多。 (3) 在二维问题和三维问题中，在同一个单元上分布的源所引起的法向速度不同。 (4) 如果速度势由满布船体表面的兴波源来表达的话，解并不在所有的摇荡频率下都存在。对于浸入水表面的船体来说，