波浪理论以及工程应用08汇编.ppt

直接边界元法同其它算法的联合应用 边界元法 ?  间接边界元法同其它算法的联合应用 外域DII采用间接边界元法的边界积分方程： 对边界离散后，上式给出如下矩阵方程： 边界元法 ?  间接边界元法同其它算法的联合应用 利用交界面条件，计入分布源内外域的平衡条件，由上式： 及： 由此，联合内外域的矩阵方程，有 由此，可解得内域的?值和边界上的分布源P. 根据分布源P可以计算确定外域和边界上的?和q值。 边界元法 ?  * 波浪理论及其工程应用 船舶工程学院 孙 雷 * 1 边界积分方程 高斯公式：设空间闭区域?是由分片光滑的闭曲面?所围成，函数P、Q、R在?内具有一阶连续偏导数，则有 这里cosa 、cosb 、cosg是?上点(x, y, z)处的法向量的方向余弦． 这两个公式称为高斯公式． 矢量形式： * 1 边界积分方程 格林公式：令 ，有 代入高斯公式得 令 ，有 两式相减得 若 格林第二公式 * 1 边界积分方程 若令 ，其中P为域中任一场点，Q为域中变点。 由格林第二公式 表示在Q点放置的单位强度点源在P点诱导的速度势 表示在Q点放置的单位强度偶极在P点诱导的速度势 在 上： * 1 边界积分方程 当P在边界上，有 光滑边界 不光滑边界 ?为P点处曲面所包围的立体角 整理得积分方程 * 1 边界积分方程 * 1 边界积分方程 1/r称为Rankine源格林函数。通常使用的格林函数的形式都可以表示成Rankine源格林函数加上一个在计算域内处处调和的函数 。即 二维问题： 通过选择合适的 ，可以使G满足各种边界条件如自由水面边界条件、海底边界条件和无限远辐射边界条件，这样方程中的积分可以只在物体表面计算。 * 2. 间接边界元法 2.1 边界积分方程 边界上点 Q 作用单位源 (Q)，对域内一点P 的诱导场函数为： * 2.1 边界积分方程 * 2.1 边界积分方程 * 3. 间接边界元法 2.1 边界积分方程 若在全边界上分布源 ，则在求解域内产生的场函数为： * 2. 间接边界元法 2.1 边界积分方程 在边界上给出一类或二类边界条件，即给出速度势或速度势的法向导数，可以构建积分方程 * 边界元法 积分方程可以通过将边界面划分为有限个数的单元离散化，然后进行数值求解。 * 边界元法 Hess-Smith 法 Hess和Smith (1964) 首先使用平面四边形单元离散计算域求解势流问题。单元上源强分布为常数，积分方程在单元中心点上满足。 Hess-Smith 法易于实现，但是对复杂表面的物体，单元之间有缝隙，单元间源强分布不连续，影响计算精度。 * 对四边形常数单元影响系数矩阵C(格林函数)和D(导数) 奇异性处理 1 当源点和场点接近的时候， 格林函数会出现奇异性。 当场点接近单元时，可以 将计算单元划分为多个子 单元的方法进行计算 。 自适应加密网格，当 某个 子单元的计算精度满足要 求的时候就不再继续细分 标准：d/C a a=0.2 奇异性处理 2 当场点在单元内时，可以通过三角极坐标变换，将奇异性 降低一阶。 方向 积分方程离散 积分方程离散 右下 积分方程离散 积分方程离散 波浪与物体相互作用计算 均匀圆柱上的水平波浪力 圆柱半径为a 水深h=a 波浪与物体相互作用计算 漂浮半球上的波浪力 半球半径为a 水深h=3a 水平波浪力 垂直波浪力 水面上漂浮方箱的水动力计算 方箱的长和宽均为90m，吃水为40m 计算结果与Faltinsen的计算结果比较 纵荡附加质量 纵荡阻尼系数 水面上漂浮方箱的水动力计算 计算结果与Faltinsen的计算结果比较 纵荡波浪力 升沉波浪力 纵荡运动幅值 升沉运动幅值 实际工程问题是多种多样的和复杂的，通常可以联合几种方法来解决同一个工程问题。 直接耦合问题 两个子域应用的不同数值方法的联合，必须交界面条件来实现。对于势流问题，交界面条件为： 边界元法与其它算法的联合应用 直接边界元法同其它算法的联合应用 边界元法 讨论外边界为SII所围场势问题。 外域DII采用直接边界元法，内域DI采用有限元法。 两个子域的边界为SI, 为光滑的。 ? 外域DII采用直接边界元法的边界积分方程： ? ? 对边界SI和SII离散后，可以得到边界元方程： 直接边界元法同其它算法的联合应用 边界元法 对边界SI和SII离散后