基于Davenport风速谱两点时程模拟
来源:互联网 发布:java算法题 编辑:程序博客网 时间:2024/04/30 19:34
转自:http://blog.sina.com.cn/s/blog_a1d5b9ba0102wfrm.html
clear all;
clc;
N=500;
omegaup=4*2*pi;
dm=omegaup/N;
dt=3.2*pi/(2*omegaup);
k=0.00464;
d=0.001;
f=d:d:10;
v10=26.8;
x=1200*f/v10;
s=4*k*v10*v10.*x.^2./f./(1+x.^2).^(4/3);
%x=1200*f/v10;
%Sv_Davenport=4*x.^2./(6*(1+x.^2).^(4/3));
%x1=z*f/(v10*(z/10)^2);
%Sv_Simiu=200*x1./(6*(1+50*x1).^(5/3));
%Sv_Kaimal=200*fn./(6*(1+50*fn).^(5/3));
%x2=1800*f/v10;
%Sv_Harris=4*x2./(6.677*(2+x2.^2).^(5/6));
%Lu=100*(z/30)^(1/2);
%x3=Lu*f/v0;
%Sv_Karman=4*x3./((1+70.8*x3.^2).^(5/6));
z1=10;
z2=50;
r=0.22;
v5=v10*(z2/z1)^(r);
Cz=10;
v1=zeros(2*N,1);
v2=zeros(2*N,1);
thta1=rand(N,1);
thta2=rand(N,1);