人工地震波依据三角级数法Word下载.docx

人工地震波依据三角级数法Word下载.docx

《人工地震波依据三角级数法Word下载.docx》由会员分享，可在线阅读，更多相关《人工地震波依据三角级数法Word下载.docx（3页珍藏版）》请在冰点文库上搜索。

∆ω——频谱分度（rad/s）。

SaT（ω）——给定的目标加速度反应谱。

φk——均匀分布在0~2π之间的随机数。

具体matlab程序如下：

%形成人工波主程序

w=[0.04:

0.02:

0.1,0.15:

0.05:

3.0,3.2:

0.1:

5.0]'

;

%频谱范围

wn=length（w）;

TT=30;

%持时

dltw=2*pi/TT;

%Δw

ag=zeros（30/0.02+1,1）;

sw=0;

kist=0;

forn=0:

30/0.02

ckn=0;

fori=1:

wn

ck=sqrt（4*sw1（w（i））*dltw）;

ckn=ck*cos（w（i）*n*0.02+rand

（1）*2*pi）+ckn;

end

ag1（n+1）=ckn;

ag2（n+1）=ft（n*0.02）;

ag（n+1）=ft（n*0.02）*ckn;

t=0:

30;

subplot（221）

plot（t,ag）

title（'

地震波'

）

xlabel（'

30（s）'

ylabel（'

ag（m/s2）'

subplot（222）

plot（t,ag1）

功率谱密度函数（随机后）'

∑Ck*cos（）（m/s2）'

subplot（223）

plot（t,ag2）

强度包线'

ξ（m/s2）'

%强度包线子程序

functionksit=ft（t）

t0=2;

tn=10;

c=0.2;

ift>

=0&

&

t<

=t0

ksit=（t/t0）^2;

elseift>

t0&

tn

ksit=1;

else

ksit=exp（-c*（t-tn））;

%计算功率密度函数子程序

functionsw=sw1（w）

T=2*pi/w;

ksi=0.05;

r=0.9+（0.05-ksi）/（0.5+5*ksi）;

eit1=0.02+（0.05-ksi）/8;

eit2=1+（0.05-ksi）/（0.06+1.7*ksi）;

Tg=0.4;

%第二类场地第二组

amax=0.8;

%七度区多遇地震

if（T>

0&

T<

0.1）

st=0.45*amax/eit2+（amax-0.45*amax/eit2）/0.1;

elseifT>

0.1&

Tg

st=amax;

Tg&

5*Tg

st=（Tg/T）^r*amax;

st=（0.2^r-eit1/eit2*（T-5*Tg））*amax;

p=0.9;

%P为反应超越概率，一般取0.85<

P<

1

sw=ksi/（pi*w）*st^2/log（-pi/（w*T）*log（1-p））;