波浪与结构物作用分析的一种高阶边界元方法_自由项和柯西主值积分的直接数值计算.doc

波浪与结构物作用分析的一种高阶边界元方法 ———自由项和柯西主值积分的直接数值计算 斌1 ,勾 莹1 ,宁德志1 滕 (11 大连理工大学 海岸和近海工程国家重点实验室 ,辽宁 大连 116024) 摘要 : 采用直接数值计算方法计算了势流问题高阶边界元方法中的自由项系数和柯西主值积分 , 建立了波浪与结构物作用的一种高阶边界元方法 . 通过算例研究了物体表面上固角系数的计算精 度和不同网格剖分 、不同阶高斯积分点对柯西主值积分的影响. 对截断圆柱上的波浪作用力与解析 解做了对比 ,发现本方法具有很高的计算精度 ,随网格的加密迅速收敛于解析解. 关键词 : 波浪力 ; 边界元方法 ; 奇异积分 文章编号 : 025324193 (2006) 0120132207 中图分类号 : U661 文献标识码 : A 由于源点在单元的节点上 ,格林函数的法向导数项 有时不再为 0 ,而且在源点周围某个单元中的积分 可能不再存在 ,因而对其计算需作特别处理. 另外 , 在结构物的非光滑部位 ,自由项系数也不再是 1/ 2 , 而需根据结构物的表面形状具体计算 . 实际计算中 一般尽可能地避免直接计算上述的 C PV 积分和自 由项系数 ,而是采用某种间接的方法加以计算. 对于 波浪与海洋工程结构物的作用问题 , Li u 等[ 1 ] 应用 了一个辅助势方法 ,在流域内间接地确定高阶边界 元方程中的主对角线元素 ; Taylo r 和 Cha u[ 2 ] 通过补 充一个物体内部积分的方法来消除固角值和积分函 数的奇异性 ,但需增加一个内水面积分 ,这带来了额 外的前期准 备 工 作 ; Te ng 和 Taylo r [ 3 ] 通 过 适 当 地 选取内部积分的格林函数来消除物体内部的水面积 分 ,但是对于带有阻尼层消波的数值水槽和应用多 极子展开加速的高阶边界元方法等 ,间接方法确定 主对角线元素往往不是十分便利 . 另外 ,对于波浪对 结构物的非线性作用问题以及波浪 、水流共同对结 构物作用问题 ,柯西主值积分还存在于物面强迫项 和自由水面强迫项的积分中 ,这也需要建立一个直 1 引言 对波浪与大尺度海洋结构物的相互作用可近似 地采用势流理论研究 ,在势流假定下可把规则波浪 对物体的绕射和辐射问题归结为拉普拉斯方程边值 问题的求解 . 对于无限水域中波浪与结构物的作用 , 边界元方法较之有限元方法或差分方法等有着显著 的优点. 若选用适当的格林函数 ,对边界元方法仅需 在物体表面上剖分网格和布置未知量 ,这将极大地 降低计算机的内存需求 、计算量和计算的前期准备 工作. 边界元方法的一个典型特点是要计算奇异积分 和自由项系数 ,在线性问题中这些奇异积分和自由 项系数构成了线性方程组的主对角线元素 . 当源点 与场点无限接近时格林函数和其空间导数将趋于无 穷 ,可以证明格林函数的奇异积分是存在的 ,但格林 函数空间导数的积分有时只在柯西主值 ( C PV) 意义 下存在. 在常数板元方法中物体表面被离散成平的 板块 ,点源布置在板块的中央 ,格林函数的法向导数 恰好为 0 ,自由项系数为 1/ 2 . 在高阶边界元方法中 收稿日期 : 2005203214 ; 修订日期 : 20052062101 基金项目 : 国家自然科学基金资助项目 ( . 作者简介 : 滕 斌 ( 1958 —) ,男 ,山东省日照市人 ,教授 ,从事波浪与结构物作用的研究. E2mail : bt eng @dl ut . edu . cn  接的求解方法. 本文实现了自由项和柯西主值积分的直接数值 计算 ,建立了波浪与结构物作用问题的一种高阶边 界元分析程序. 文中开展了算例的数值研究工作 ,研 究了方箱表面上不同点处的自由项系数 ; 考察了柯 西主值积分在不同单元剖分方 式下 的数 值积 分 结 果 ;对一截断圆柱上的波浪力与解析解做了对比研 究 ,结果表明本文方法的计算结果收敛速度很快 ,与 解析解吻合良好. 个刚体运动模态下的运动分量 , 为待定量 . 经过上述 分解 , 绕射势 7 ( x , y , z ) 和辐射势 j 在物面上满足 的边界条件为 nj , j = 1 , ?, 6 9 j 在 S b 上 , ( 7) = 9 0 , 9 n = 7 - j 9 n 式中 , n 为物面单位法向量 ,指出流体流动方向为正 向 ;在无穷远处散射势还满足 So mmerfeld 散射条件. 满足上 述 边 界 条 件 的 基 本 解 G ( x , x0 ) 已 由 J o h n 于 1950 年导出 . 当源点 x0 位于物体边界时 ( 图 2) ,对格林函数和速度势 j ( j