随机波浪的模拟1

i1随机波数值模拟方法概述研究海浪及其对工程的作用有三种途径：一是现场观测研究；二是在实验室内进行模拟研究；三是理论分析研究。由于海浪的复杂多变性，加上现场环境恶劣，进行现场观测需花费大量的人力物力；理论研究目前也有较大的局限性，特别是对于不 规则 编码规则下载淘宝规则下载天猫规则下载麻将竞赛规则pdf麻将竞赛规则pdf 波浪，很多问题有赖于室内的模拟研究。模拟研究的方法可分为两大类。开始是在水槽或水池内利用风或造波机进行物理模拟，亦即进行波浪模型试验。在人们的精心设计下，可以把负责的现象分解为多个简单的模型，然后再把成果综合起来。过去已取得了大量的研究成果，目前仍是主要的研究方法之一。随着电子计算机的发展和普及，海浪的数值模拟得到迅速的发展，它具有经济方便等优点，日益受到人们的重视和广泛的应用。天然海浪是很复杂的，人们对它的认识和研究过程是由简到繁，由浅入深，及即由单向规则波—斜向规则波—单向不规则波—多向不规则波。不规则波浪的数值模拟—模拟频谱单向不规则波浪的数值模拟方法，大多建立在线性波浪理论的基础上。本文主要介绍利用线性叠加法和线性过滤法进行二维不规则波的模拟。线性波浪叠加法在工程中，如果已经得到了特征波的波参数如有效波高Hs、周期T等参数，如何得到一列不规则波面时间历程呢？一般通过模拟靶谱法来完成。将有效波高Hs、周期T等参数代入某波浪频谱形式中，得到的海浪谱即为靶谱。现在要模拟某波面不规则波面时间历程，使得模拟的波谱同靶谱一致。平稳海况下的海浪可视为平稳的具有各态历经性的随机过程，波动可以看作无限多个振幅不等、频率不等、初相位随机的简谐余弦波叠加而成，即Mtaicoskixiti(1式中，t为波动水面相对于静水面的瞬时高度；ai为第i个组成波的振幅；ki,i为第i个组成波的波数和圆频率；ki2Li,i2TiL，T分别为波长、周期；x，t分别 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示位置和时间，通常固定位置，可取x=0；i为第i个组成波的初位相，此处取在（0,2π）范围内均布的随机数。通过频谱来模拟海浪，设欲模拟的对象谱（靶谱）S的能量绝大部分分布在L~H范围内其余部分可忽略不计。把频率范围划分为M个区间，其间距ai为iii1，取?ii1i2，则第i个组成波的振幅为(2)则将代表M个区间内波能的M个余弦波动叠加起来，即得海浪的波面：M(3)t2S?iicos%itii1式中，~i为第i个组成波的代表频率。用波浪叠加法模拟海浪时应注意以下几点：频谱范围L~H的选取频谱范围L~H的选取，取决于所要求的精度。设在频谱高低侧各允许略去总能量的部分（例如取0.2%），对于可积分的谱，易于确定L和H若采用公式SL0.785exp3.114Hs2表示的P-M谱，可以得到3.11Hs2ln1/4，H1/43.11Hs2ln1(4)对不可积分的谱，可以采用数值计算的方法来确定H。首先采用数值积分的方法计算波浪频谱的总能量E，然后计算对应每个频率i的累积能量Ei，则Ei/E对应的频率即为下限L，Ei/E1对应的频率即为上限H应该看到，在M一定的情况下，不恰当地增大谱频范围，反而会使精度下降。一般取谱峰频率的3~4倍作为H已足够。图1划分波谱的频率区间示意图频率区间的划分划分频率区间的方法，有等分频率和等分能量法。2.1.2.1等分频率法下面简要介绍下等分频率法。取HLM（一般取M=50~100）。但若采用式?ii1i2中的?i作为i区间的代表频率，则由式（3）模拟所得的波浪将以周期2重复出现，除非值足够小；否则与实际的海浪情况不符。应在各区间内部随机选取频率作为该区间的代表频率~i。~i的选取方法对模拟结果有相当的影响。由于波能集中在谱峰部，如M值较小；只有少数位于谱峰处的组成波起主要作用，可能产生较大的误差。等分能量法定义累积谱为ES0d(5)如果按照等分能量法分成N份，则分界频率i可以用下式来确定。EiiEim0(6)NN对P-M等分可积分的谱，则B1/4(7)ilnN/i各组成波的振幅a相等2m0Nai2S?ii(8)此时式tancosi1ntn变为t2m0NcosNi1?iti(9)?i1i2(10)2.1.3随机相位的选取随机初位相i应在0~2区间内均布。如组成波数M不很大，则由计算机产生的随机数往往不够均布，影响模拟结果。我们采用人造的比较均匀的随机数，模拟结果较好。合田采用M=200，由计算机产生随机数（每次不同）进行多次重复计算，对结果进行统计分析，取其特征值。线性过滤法应用线性滤波法模拟海浪的基本思路是：以白噪声为一线性系统的输入，通过选择适当的系统函数使该系统输出的谱恰恰等于靶谱海浪等随机过程由多种不同频率的成分组成，他们可以通过不同的滤波器分离开来。如图2所示，只有高频信号能通过高通滤波器，通过低通滤波器的是低(15)频信号，允许一定频率范围内的信号通过的滤波器称做带通滤波器图2滤波原理示意图具有如图3中所示传递函数的滤波器称做成型滤波器。这些滤波器可以是数字式的，也可以由硬件组成。线性系统的输入谱S*f和输出谱S*yyf之间存在下列关系：2S*yyfTf2S*xxf(11)白噪声的谱密度为常数，且可等于1，如将它作为输入，通过按靶谱设计的成型滤波器后，即可得到谱形符合靶谱的随机波浪。因此线性过滤法的关键在于靶谱设计过滤器。过滤器的选择。输入白噪声的谱S*xxf1，要模拟的波浪靶谱为S*f(双xx侧谱)，由上式得过滤器的传递函数为TfS*fSf2(12)在时域，线性系统的输入。输出函数间有关系即txthd(13)h是脉冲响应函数，也是过滤器的权函数，其傅里叶变换即为传递函数T(f)，即(14)Tfei2fdf写成离散形式：Lajxtjtt0,t,2t,...,NtjL图3用过滤法模拟波浪示意图模拟不规则波浪。将上式代入可得到所要的波面。为便于计算，把它改写成LtA0xLtAjxLijxLij(1式中，xLt相当于xtLt，可取L=20~30白噪声x(t)可用一系列独立的正态分布的变量x1,x2,...来接近，这些变量的均值为零，方差为1。可按下式得到：nxk2RANi13n，k1,2,3,...N2Li1RANi为在(0,1)区间内均布的伪随机数，一般计算机可直接产生。可取n=30~50。程序实现3.1程序一：线性波浪叠加法模拟频谱%%不规则波浪的数值模拟—模拟频谱%%线性波浪叠加法t=0:0.01:1000;%时间间隔x=0;%初始尾椎a=abs(randn(1,3));%幅值T=abs(randn(1,3));%周期L=abs(randn(1,3));%波长c=abs(rand(1,3))*2*pi;%初相位A=0;fori=1:3A=A+a(i)*cos(2*pi*x/L(i)-2*pi*t/T(i)+c(i));endplot(A);set(gca,'xlim',[0100000]);set(gca,'ylim',[-55]);xlabel('\itTime');ylabel('\itAltitude');title('\bf线性叠加法模拟频谱');gridon%%end图4线性波浪叠加法模拟频谱图（N=3）图5线性波浪叠加法模拟频谱图（N=10）图6线性波浪叠加法模拟频谱图（N=50）图7线性波浪叠加法模拟频谱图（N=200）图8线性波浪叠加法模拟频谱图（N=1000）程序二：线性过滤法%%随机脉冲函响应数num=[0,0,25];den=[1,4,25];[y,~]=impulse(num,den);plot(randn(1)*y,'r:','linewidth',1.5);set(gca',xlim',[0127]);ylabel('\itAmplitude');title('\bfRandomImpulseResponcse)'gridon图9生成随机脉冲响应函数图%%生成传递函数T=fft(y,127);plot(real(T),':','linewidth',1.5);set(gca',xlim',[0127]);ylabel('\itAmplitude(Real/Imag)');holdonplot(imag(T),'r:','linewidth',1.5);title('\bfTransferFunction');gridon图10生成传递函数实部、虚部图%%A(j)M=127;dt=1/(127*2);t=0:dt:1000;F=1/(2*dt);df=F/M;L=20;A=0;forj=1:Lfori=1:MA=A+1/M*real(T(i*df))*cos(1*pi*i/M);B(j)=A;A=0endendstem(B,'filled');ylabel('\itAmplitude');xlabel('\itj');title('\bfA(j)');gridon图11生成A(j)数值图%%白噪声的模拟x=normrnd(0,1,1,167);plot(x,'linewidth',1.5);ylabel('\itAmplitude');title('\bfWhiteGaussianNoise)';set(gca,'xlim',[0167]);gridon图12模拟高斯白噪声图%%不规则波浪的模拟A0=randn(1);N1=0;N2=0;fori=1:Mforj=1:Lif(i-j<=0)N1=N1+A(1)*(x(i+j));elseN2=N2+A(1)*(x((i+j))+x((i-j)));endendendN=N1+N2;t1=t+20*0.039;z=round(t1);H=A0*x(z)+N;fori=1:25400;if(H(i+1)==H(i))H(i+1)=H(i);H(i)=0;endendH=H(H~=0);plot(H,'linewidth',1.5);ylabel('\itAmplitude');xlabel('\itTime');title('\bf线性过滤法模拟频谱');set(gca,'xlim',[0100]);gridon%%end图13线性滤波法模拟频谱图3线性海浪数值模拟的其他方法除上述的海浪模拟的主要方法外，还有一些其他方法。Hino等(1972)根据平稳过程线性预测的Wiener-Kolmogorov理论模拟海浪波面高度。在海面微波散射的研究中，有些作者提出将分形理论应用于海面特征的模拟，如Lo.T(1993)分析了分形海面散射特性，并将其应用于海面目标探测，姚纪欢等(1999)从海面的功率谱出发,建立了海面的Weierstrass分形函数模型。另外还有一些作者通过直接解算水动力学方程来对海浪进行数值模拟。如JimX.Chen等通过解算二维Nevier-Stokes求得水面网格节点的离散速度场，根据无旋理想流体的Bernoulli方程将速度场与水面高度场联系起来，改变速度场结点上的值就可达到模拟波面起伏的三维动画效果，并具有实时性。虽然这种方法是经人为处理用于视觉模拟的，但将流体动力学理论与随机海浪理论联系起来研究海浪的数值模拟是一个值得深入研究的方向。在海岸工程问题上，基于缓坡方程或Boussinesq方程对近岸海浪及海浪传播变形的模拟一般均采用有限元或边界元直接数值计算的方法，本文对此不再赘述结果分析利用线性波浪叠加法和线性过滤法都等实现对不规则波浪的数值模拟。线性叠加法较线性过滤法简洁些，但是产生的波浪的波动性比较大。线性过滤法相对复杂些，但是可以比较准确的模拟实际海浪。因此，我们可以利用线性叠加法来进行具体工程问题的简要估计，而用线性过滤法来进行比较准确的数值模拟，使更接近于实际情况。参考文献俞聿修.海浪的数值模拟[J].大连理工大学学报,1981(2):84-90.田子谦.波浪的数值模拟技术[J].建筑结构,2010(s1):414-416.王子健,肖盛燮.基于线性变幅叠加法的不规则波浪力数值模拟[J].人民黄河,2013,35(2):128-130.聂卫东,康凤举,褚彦军,等.基于线性海浪理论的海浪数值模拟[J].系统仿真学报,2005,17(5):1037-1039.刘思,张永良.多向不规则波群传播的数值模拟[J].清华大学学报:自然科学版,2015(12):1289-1295.洪广文.不规则波随机数值模拟及其工程应用[J].河海大学科技情报,1989(3):74-107.王修琼,贺国京.工程中随机过程的数值模拟方法及其评述[J].四川建筑科学研究,2005,31(1):26-29.环飞,冯卫兵,王爽.浅水非线性不规则波数值模拟[J].服装学报,2012,11(5):575-581.梁修锋,杨建民,李欣,等.不规则波浪的数值模拟[J].船舶力学,2010,14(5):481-486.竺艳蓉.海洋工程波浪力学[M].天津大学出版社,1991.