基于malab的随机波浪模型仿真研究.docxVIP

基于malab的随机波浪模型仿真研究

如果在海域使用浮式结构，并不可避免地发生各种不良的海况，那么就不可避免地会发生各种不良的海况。浮式结构物受到海浪产生的绕流力、惯性力、冲击力、浮托力等作用,这些作用直接关系到结构物的生存安全。因此在浮式结构物设计时,必须要充分考虑海浪的影响,精确计算海浪载荷,校验和优化结构,研究浮式结构物的耐波性和稳性,保证其满足设计要求的稳定性和安全性。

目前评估海浪载荷对浮式结构物的影响主要采用物理模型试验的方法,然而这种方法费时费力,而数值模拟具有参数设置灵活、计算结果精确等优点,正逐步成为设计领域研究海浪载荷的重要手段。由于受到计算机硬件发展水平和波浪理论不成熟等的制约,早期的海浪数值仿真主要以二维为主,但是对于需要研究海浪和浮式结构物相互作用过程中产生的波浪折绕射、漩涡等现象时,显然不能满足工程的需要。因此研究三维海浪数值模型,实现对海浪现象更加真实准确的描述,是海上复式结构物设计领域中研究海浪和结构物相互作用的必然发展方向。

本文利用谱分析的方法,在MATLAB环境下对三维随机海浪模型进行了数值模拟仿真,并给出了三维随机海浪波面图,为浮式结构物设计中计算海浪载荷提供了参考。

1海洋频谱仿真

海浪现象十分复杂,具有随机性和非线性的特征,因此要建立精确的海浪模型是十分困难的。根据海浪谱和随机海浪理论,可将实际海浪看成为不同频率、不同传播方向、不同波高和不同初始相位的正弦波叠加的结果。Longuet-Higgins模型就是把长峰波海浪的波动用无数个随机余弦波的叠加来描述,其波幅的表达式为:

其中ξai、ωi和εi分别表示第i个余弦波的波幅、频率和相位。

在实际应用中,Longuet-Higgins模型常用海浪的频谱来表达,即:

ζ(t)=i=∑∞12SζΔωcos(ωi+εi)

式中Sζ(ωi)为海浪的频谱。目前应用较多的海洋频谱有Pierson-Moscowitz谱(PM谱)、Neumann谱(N谱)、ITTC双参数波谱等,其中PM谱应用最为广泛,能很好地代表实际充分发育的随机海浪,故本文采用PM谱。PM谱的表达式为:

其中v为离海面19.5m处的平均风速。

平均风速为10m/s、12m/s、15m/s的PM谱的仿真曲线如图1所示。从图1中可以看出,海浪的PM谱为窄带谱,能量集中于某一频段内,随着风速的增加,即高海情情况,PM谱峰值增大,能量更为集中,如当风速达到15m/s时(六级海情),谱峰值为3.6m2?s,能量集中于0.3-1.7Hz频带内。因此对海浪进行数值仿真时,可以选取能量集中的频段内的有限个谐波成分进行叠加,可减少计算量,从而加快仿真速度。

2不规则短峰波的三维模拟

2.1维不符短峰波海藻的波面方程

Longuet-Higgins模型假设海浪只沿一个固定方向传播,且假设其波峰和波谷线相互平行且垂直于前进方向。实际上,海面上的波浪不仅波高不同,频率不同,还会向各个方向传播。这些谐波除产生在主风向上的主浪向以外,在主浪向两侧的±π2角度范围内都会有谐波的扩散,这样的海浪称为三维不规则短峰波海浪。

与Longuet-Higgins模型不同,三维不规则短峰波海浪是由振幅不等、频率不等、初始相位不同、传播方向不同的谐波多重求和得到。双叠加模型为:

其中ζaij、ωi、ki、μj、εij分别为组成谐波的波幅、角频率、波数、方向角和随机初相位,εij为[0,π2]内的随机数,(ξ,η)为波面上某点的坐标。

根据谐波的波幅和频谱的关系,三维不规则短峰波海浪的波面方程可进一步写为:

其中Sζ(ωi,μj)为短峰波海浪的方向波谱。

通常认为波能的方向和频率是无关的,则短峰波海浪方向波谱可用两个独立函数的乘积来表示,即:

其中Sζ(ω)为海浪的PM谱,μ为海浪的方向角,?(μ)为波能的扩散函数。ITTC提出的建议性形式为:

波数ki表示谐波的传播方向上2π距离内波的个数,于是有ki=2πλi,其中λi为第i个谐波的波长。因为波长不易获得,由线性波动理论可知,在不考虑表面张力的情况下,数值上的波数和角频率有如下关系,即:

式中g为重力加速度,D为水深。

2.2海域模型材料仿真

某浮式结构物的作业和生存海域的海洋条件为:风向30°,风速13.8m/s,五级浪,为校验该结构物所能承受的波浪载荷,需要对该海情下的海浪进行模拟。

选择100m×100m大小的海域,利用上述理论进行三维随机波浪仿真。各参数设置为:仿真频段[0.3Hz,1.7Hz],海浪的方向角μ=30°,风速v=13.8m/s,时间t=10s。图2为该海域的模拟海浪图。