基于精细时程积分方法的时域格林函数计算

文章提出了一种精确计算三维时域格林函数记忆项本身及其导数的方法。基于精细时程方法求解时域格林函数及其导数所满足的定常非齐次常微分方程组,得到时域格林函数数值解。文中通过与解析方法和其他方法的比较,验证了该方法的精确性、稳定性和高效性。     

海洋结构物时II边界积分方程的改进方法

文章基于三维时域格林函数理论,提出了一种时域边界积分方程的改进计算方法,用于预报海洋结构物在波浪中运动受到的时域波浪力.该方法通过利用物体内部构造的满足拉普拉斯方程的内部解,在时域边界积分方程中只需要计及时域格林函数记忆项本身,从而有效地避免了求解具有高震荡积分特性的时域格林函数记忆项的法向导数.对潜体和浮体作了计算,结果表明,这种方法十分有效,可以有效节省内存使用量和计算机空间.通过与相关已发表文献的结果进行比较,计算结果令人满意.     

船舶有航速时三维时域格林函数计算方法研究

船舶工业CAE软件势流求解器的开发对势流计算方法的可靠性提出要求,本文针对船舶有航速时深水三维时域格林函数计算方法开展研究。首先引入三维时域自由面格林函数,其瞬时项的求解采用Hess&Smith方法,自由面记忆项的求解采用美国密歇根大学Beck团队的方法。其次,分区域推导出三维时域自由面格林函数记忆项及其导数的编程所用的数学展开表达式。然后,采用间接法模型求解分布源强和速度势,并采用压力直接积分法求出时域辐射力。最后,针对Wigley I船型开展了辐射力计算,并与公开的试验及计算结果进行对比,验证本文方法的可靠性。     

基于时域格林函数法的船舶有航速问题数值求解

为了解决船舶以恒定速度在波浪中运动的水动力问题,文章开发了基于瞬时格林函数方法的三维时域计算程序。从数值精度和计算效率的角度出发,采用了精细时程法和一些数学上的推导计算时域格林函数的波动部分及其导数。并对Wigley型船进行了数值计算,通过水动力系数与试验数据及文献结果的比较,验证了三维时域计算程序的可靠性。     

船舶有航速辐射问题的新实用数值方法

为解决船舶在波浪中以恒定速度运动的辐射问题,基于势流理论,采用三维瞬时自由面格林函数方法对船舶水动力问题进行求解。提出三维瞬时自由面格林函数波动项的新数值计算方法,并在大时间区域内基于精细积分法对三维瞬时自由面格林函数的波动项进行插值调用,最终提高波动时域格林函数的数值计算精度。对Wigley船型的水动力系数进行数值计算,并将附加质量和阻尼系数的数值计算结果与试验结果和参考文献进行对比,验证新数值计算方法的有效性。     

基于三维时域格林函数高阶面元法的船体运动计算

针对波浪中航行船体的运动问题,基于三维时域自由面格林函数并采用高阶面元离散,编程计算了浮体的运动响应.采用自由面格林函数,只需在船体湿表面布源,减少系数矩阵维度,避开自由面网格划分;同时高阶面元法(HOBEM)在计算速度势空间分布时能够提高精度.本文通过计算Wigley-Ⅲ型船在波浪中有航速前行时的运动,与实验进行对比,发现结果符合较好,证明这种方法的有效性.     

时域有限水深格林函数及其导数的数值计算

精确快速地求取有限水深格林函数及它的导数是应用时域法求解有限水深水动力学问题的关键.通过对有限水深格林函数在不同区域渐近特性的分析,给出了使用多维Chebyshev多项式和渐近展开式快速近似计算有限水深格林函数及它的导数的方法.近似方法与直接积分方法的比较表明近似方法的计算工作量大幅减少,可以给出足够精度的计算结果.     

时域格林函数的计算（英文）

时域格林函数及其导数的准确计算对时域耐波性分析至关重要。本文基于Dawson积分和泰勒展开,提出了一种新的有效计算时域格林函数及其导数的方法。使用这种新方法计算了作用于一个进行大幅运动的潜没圆球上的垂荡力和一个以定速前进的潜没圆球的纵荡力,计算解和吴国雄给出的解析结果符合良好,说明本方法是准确的和有实用价值的。     

系泊Spar平台波浪中耦合运动的数值模拟及模型试验

为了研究恶劣环境载荷引起海洋结构物的非线性运动响应,文章基于三维时域格林函数理论,提出了采用时域物面非线性理论方法直接模拟系泊浮体时域耦合分析所需的水动力,建立了系泊浮体波浪中时域耦合运动的数学模型.该方法充分考虑了每时每刻浮体在瞬态波面下湿表面的变化,并计及平均水线面下瞬时空间位置的改变,系统研究了系泊Spar平台波浪中的耦合运动响应及系泊线张力.最后,将数值模拟结果与模型试验结果进行了比较,计算结果令人满意.     

斜壁结构的水动力收敛方法研究

在分析波浪中的船舶运动或者计算大型结构的水动力系数时,往往采用时域格林函数方法。时域格林函数的一个重要局限性在于它在计算具有倾斜壁面的水动力系数时,结果很容易发散。时域格林函数本身的奇异性以及高频振动特性显然是水动力发散的一个重要原因。但即使该奇异性通过增加粘性以及表面张力的方式使之消失,计算具有斜壁结构的水动力时,发散现象依然存在。因此,该文提出一种滤波方法,除去时域格林函数的高频部分,留下其低频部分,并定义滤波系数,从而使作用于斜壁结构上的水动力值收敛。通过文中结果与频域兴波格林函数法的比较来确定最优滤波系数。结果表明：最优滤波系数几乎与运动幅值无关,但会受到运动频率以及物面形状的影响。     

考虑航速及自由液面影响的声介质中三维结构水弹性力学研究

基于Wu建立的带航速三维水弹性理论,采用考虑自由液面效应的理想可压流体Green函数,提出了带航速和考虑自由液面的声介质中三维水弹性结构声辐射计算方法.对THAFTS三维水弹性计算程序进行修改,得到声介质中三维水弹性计算程序.针对无限大流体中的静态弹性球壳结构声辐射模型进行一系列数值计算,通过与解析结果的比较,验证了该计算方法和计算程序的正确性.进一步计算不同潜深、不同航速下弹性球壳的结构声辐射,以此为基础对该计算方法和计算程序进行进一步考核,同时分析了航速和自由液面对水下结构声辐射的影响.     

波浪中航行两船模池壁干扰效应的数值研究

对规则波中顶浪并行航行的两船模,基于频域移动脉动源格林函数建立了有效表达池壁效应的三维数学模型,该模型采用分布源积分方法求解近距航行两船模的波浪作用力和运动响应,并利用镜像原理将水池格林函数表示为开阔水域的格林函数和由镜像点源组成的无穷级数.分析移动脉动源格林函数的远场传播模态及其特性可知,当τ>0.272(τ为斯特劳哈尔数)时各传播波模态均不在点源的上游出现,故此时点源波系的影响范围可由环形—扇形波的半楔角确定,在此基础上计及船模辐射—绕射波的池壁反射作用对自身和另一船模的干扰,依据船型要素、水池宽度和两船相对位置,采用六个几何临界半楔角表征船体的辐射—绕射波,提出了并行航行两船模的池壁效应判别方法.数值研究表明:不同镜像次数下水池格林函数的数值差异主要体现在距船艏较远处,在兼顾计算精度与效率的原则下可选取有限项级数的水池格林函数作为积分内核,而池壁效应对近距航行两船运动响应的影响主要体现在长波区.     

三维锚泊系统时域计算分析

通过采用在频域内对浮体运动方程求解并将计算结果转换为时域结果,同时与锚链线运动时域方程耦合求解锚泊浮体在限制工作水深的运动位移与锚链线上的张力.采用一种考虑锚链线上各种受力的三维锚链线模型来计算锚链的位移与受力.通过对计算结果的分析表明,本方法简单、实用,可以提供具有工程精度的锚泊性能预报.     

近水面水翼影响的船舶兴波时域计算研究

近水面水翼的兴波由升力线的自由面兴波来表达,船体的兴波考虑升力线兴波的干扰影响,三维时域格林函数法用于船体的兴波分析计算.以Wigley船型为数值计算例,考虑近水面水翼的影响,分别对单体船和双体船进行了船舶兴波计算研究,提供了兴波波形和兴波阻力的计算结果.以系列60船型为数值计算例,把计算的结果与发表的有关结果进行对比,具有较好的吻合.本文还对船舶自由表面兴波波形进行了时域数值仿真.     

三体船波浪诱导载荷计算研究

基于三维势流理论和简单格林函数法,探索了三体船在规则波中运动的频响函数,并对航行于北大西洋海域的全船典型剖面波浪载荷进行了长期预报,分析了不同浪向角时典型剖面的剪力、弯矩和扭矩极值,研究了三体船剖面载荷的分布规律和趋势,为三体船结构设计提供了一定的依据.     

规则迎浪中船舶参数横摇的三维时域预报方法研究

参数横摇是船舶因复原特性改变而引起的典型非线性现象。文章采用三维时域方法预报规则迎浪中船舶的参数横摇运动。该方法引入匹配面将流域分为内域和外域,内域中采用Rankine源来满足物面条件和线性自由面条件,而外域中应用时域格林函数来满足线性自由面条件和远场辐射条件。数值方法中,Froude-Krylov力和恢复力是通过对船舶瞬时湿表面积积分获得,同时考虑了横摇、垂荡和纵摇三自由度之间的耦合作用,以及非线性横摇阻尼的影响。数值结果与试验结果吻合很好,说明该方法可以有效地预报参数横摇。     

三体船连接桥结构波浪载荷研究

根据三体船船型简要分析连接桥结构波浪诱导的剖面载荷特征。以三维势流理论为基础,应用简单格林函数法对三体船在波浪中运动与载荷的时历进行仿真计算,并将时域结果转换为频域,之后对航行于北大西洋的三体船的连接桥波浪载荷进行长期预报,最终探索三体船连接桥纵向剖面的各载荷分布规律和趋势。     

采用THOBEM方法研究波流与结构物相互作用的数值模拟(英文)

To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method (HOBEM). By assuming small flow velocities, the velocity potential could be expressed for linear and higher order components by perturbation expansion. A 4th-order Runge-Kutta method was applied for time marching. An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves. Validation of the numerical method was carried out on run-up, wave exciting forces, and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder. The results were in close agreement with the results of a frequency-domain method and a published time-domain method. The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.