斜向波与长棱柱型浮式结构物相互作用边界元模拟.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 斜向波与长棱柱型浮式结构物相互作用的边界元模拟 郑永红 沈永明 游亚戈 中国科学院广州能源研究所 广东广州 大连理工大学海岸与近海工程国家重点实验室 辽宁大连 摘要 采用第二类变形 函数对斜向线性入射波与任意截面形状的长棱柱型浮式海洋结构物的相互作用问题 进行了边界元法模拟 给出了常数元和线性元的奇异积分处理方法以及开边界的处理方式 对矩形和圆形结构物 的水动力系数和波浪力进行了计算 并与解析解及有限元计算结果进行了比较 结果表明 本文方法能得到满意 的结果 且未出现不规则频率现象 关键词 斜向波浪 浮式结构物 边界元法 水动力学系数 波浪力 中图分类号 文献标识码 随着海洋开发事业的发展 斜向线性波与任意截面的长棱柱型浮式结构物的相互作用问题已受到 海洋工程研究人员的关注 文献 应用有限元方法 文献 应用边界元法以及文献 采用试 验方法对此类问题进行了研究 本文将采用边界元法来进一步研究该类问题 与过去常用的边界元方 法不同的是格林函数的选取及对辐射开边界条件的处理 过去常用的格林函数是满足自由面条件 水 平海底边界条件和辐射侧开边界条件的复杂格林函数 其优点是只需在物面边界上布置计算单元 在其 他边界上不需要布置计算单元 如果海底为非平底地形 则在该海底边界上仍然需要布置计算单元 从 而可以节省大量的存储空间 但其缺点是格林函数及奇异积分的处理比较复杂 而且还可能出现不规则 频率现象 因此本文仅采用简单形式的格林函数 即第二类变形 函数 其优点是奇异积分处理 和编程实现十分简单 并且通过数值计算表明 未出现不规则频率现象 采用该方法对算例的水动力学 系数和波浪力进行了计算 与解析解及其他数值结果的比较表明 本文研制的边界元法能够得到满意的 结果 问题描述与数学模型 考虑海域中 漂浮在自由面或淹没在流体中的任意 形状棱柱型结构物 视作刚体 与斜向线性入射波的相互 作用问题 包括散射和辐射问题 其中辐射运动又分垂 荡 横荡 和横摇 种运动模式 结构 物的几何特征和坐标布置如图 所示 坐标系 为 右手坐标系 其原点定义在静水面上 轴垂直向上为 正 轴向右为正 浮体在 方向足够长 相对于入射波 图 斜向波与浮式棱柱型结构物 长 入射波的传播方向与 轴的夹角为 收稿日期 基金项目 国家重点基础研究发展计划 计划 项目 国家自然科学基金项目 辽宁省自然科学基金项 目 作者简介 郑永红 男 安徽潜山人 博士 研究员 主要从事水动力学方面的研究 假定流体不可压缩 流动无旋且浮体在波浪作用下作周期性的微幅运动 则流体状态可用速度势 进行描述 其中 表示取复数表达式的实部 为圆频率 为时 间 为与时间无关的空间复速度势 满足三维拉普拉斯方程 由于假定浮体在 方向足够长 因此可以认为流体在 方向也呈周期性变化 即有 这种假定被许多学者用于斜向波与浮体的相互作用研究 将式 代入式 得 式 即为本文待求解的基本控制方程 对于线性波 可分解为 式中 和 分别为与时间和 无关的入射波速度势 散射势和浮体做 模式运动产生的辐射 势 它们均满足式 这里 分别表示垂荡 横荡和横摇 由于 是预先给定的 因此求解 也就是求解 和 对在水深为 的水域中传播的振幅为 圆频率为 传播方向与 轴的夹角为 的斜向入射波 其复速度势 为 式中 为重力加速度 为波数 由色散关系