海上结构物周围波高的数值模拟研究.pdf

46卷增刊 中 国 造 船 V01．46 SpeciM 2005年11月 SHmBUⅡDING0FCHINA NOV2005 文章编号：1000-4882(2005)S一0227-05 海上结构物周围波高的数值模拟 滕斌，孙亮 (大连理工大学海岸和近海工程国家重点实验室，辽宁大连116024) 摘要 边界元法作为一种高效的数值方法，凭借其自身的优势，已经在土木工程和海洋工程等许多领域得到 了广泛的应用。在对海上结构物的设计过程中，可以采用边界元法对同一结构物的不同设计方案进行数值 模拟，对所得的计算结果进行比较、分析，根据工程的实际要求最终选择一个合理的结构型式。在波浪高 度的模拟中，由于积分方程的准奇异性，直接计算往往产生很大的计算误差。本文应用单元细分方法解决 了准奇异积分的问题，对物体周围的波高计算，可以得到高精度的计算结果。应用该方法本文计算了固定 在水中的均匀直立圆柱周围的波高分布，与解析解吻合得很好，验证了计算程序的正确性。接着计算了不 同截面型式圆柱周围的波高分布，比较了相同特征尺度、不同型式的柱状结构物对附近波高分布的影响。 关键词：波高；数值模拟；海上结构物设计 1前言 港口和近海工程中往往需要在岸边或海上修建一些结构物。这些结构物的优化设计是实际工程 中的关键问题。实验方法往往周期比较长，需要投入的经费也比较多，而且实验中也有许多不确定 因素；数值模拟的方法相对实验方法有很多优势，需要周期较短，只需投入少量经费，而且环境条 件可以人为给定。 为了满足实际工程的需要，众多研究者提出了许多数值模拟的方法，如有限元法、差分方法、 边界元法等。这些数值方法都各有优势。边界元法相对于有限元法和差分法这两种域内方法，具有 网格划分简单，未知数少等优点，已经在土木工程和海洋工程等许多领域得到了广泛的应用。 在海洋工程中，波高往往用来作为实际工程的一个重要参数。结构物的存在影响了波浪场的形 态，改变了海上的波高分布。本文针对工程中常见的柱状结构物，采用边界元方法计算了结构物附 近的波高分布，计算中应用单元细分的方法处理准奇异积分问题，得到了十分准确的结果。通过对 计算结果的比较，研究了几种不同型式的柱状结构物对其附近波高分布的影响。 2边界元法进行数值模拟的过程 在流体无粘、不可压缩、运动无旋的假定下，存在复速度势矽。对于规则波作用下的一阶绕射 问题n1(如图1所示)，缈=仍+仍，其中仍为入射势，仍为绕射势。采用直接边界元法把物面上 的绕射势仍作为基本未知量进行数值计算的过程如图2所示。 收稿日期：2005—07—28 中 国 造 船 学术论文 图I规则波对直立圆柱的作用 图2边界元法进行数值计算的过程 本文在计算物面上的绕射势仍的过程中，使用了如F积分方程H： 蠊，+舭磊，挈州元，挈卜2黔磊，鬻魏 求得了物面上的绕射势伤，就可以计算结构物附近的一阶波高分布。波浪场中的波高可以通过下面 公式求得： 日=2f詈缈+I=2I詈c仍++仍+，I 式中够+和仍+分别为自由水面上任一点处的入射势和绕射势，可根据该点的坐标使用以下公式求得： 仍+=一百igHi沙刚蛳肿) 仍瓴)=』』(嘶)挈+G(弧)警灿 在计算水面上任一点的绕射势仍+的过程中，当水面上的计算点距离物面较近时，积分 要嘶，誓出具有二阶准摒蝴分gG警出具有一阶准勰蛐脯接的高斯积分 对上述两个积分都很难获得准确的结果。本文应用单元细分的方法吲解决了准奇异积分的问题，可 3