人工地震波依據(jù)三角級數(shù)法

1、人工地震波依據(jù)三角級數(shù)法武漢橋梁隧道799084759采用shinozuka的方法來模擬平穩(wěn)化后的隨機(jī)地面運(yùn)動加速度Xg t f t t其中:nf (t) Ckcos ktk 1Ck .4 S kSaTln ln 1 PT2t t0t0t tntn ttt/t0強(qiáng)度包線：t 1t exp c t tnP 反響超越概率.S功率譜密度函數(shù).?3頻譜分度(rad/s).SaT (3)給定的目標(biāo)加速度反響譜.(j)k均勻分布在02兀之間的隨機(jī)數(shù).具體matlab程序如下:%形成人工波主程序w=0.04:0.02:0.1,0.15:0.05:3.0,320.1:5.0'% 頻譜范圍 wn=len

2、gth(w);TT=30;% 持時dltw=2*pi/TT;% Awag=zeros(30/0.02+1,1);sw=0;kist=0;for n=0:30/0.02ckn=0;for i=1:wnck=sqrt(4*sw1(w(i)*dltw);ckn=ck*cos(w(i)*n*0.02+rand(1)*2*pi)+ckn; endag1(n+1)=ckn;ag2(n+1)=ft(n*0.02);ag(n+1)=ft(n*0.02)*ckn;endt=0:0.02:30;subplot(221)plot(t,ag)title('地震波')xlabel('t=0:30

3、 (s)')ylabel('ag (m/s2)')subplot(222)plot(t,ag1)title('功率譜密度函數(shù)(隨機(jī)后)')xlabel('t=0:30 (s)')ylabel(' E Ck*cos() (m/s2)')subplot(223)plot(t,ag2)title('強(qiáng)度包線')xlabel('t=0:30 (s)')ylabel(' E (m/s2)')%強(qiáng)度包線子程序function ksit=ft(t)t0=2;tn=10;c=0.2;if t

4、>=0&&t<=t0ksit=(t/t0)A2;elseif t>t0&&t<tnksit=1;elseksit=exp(-c*(t-tn);endend%計(jì)算功率密度函數(shù)子程序function sw=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;%七度區(qū)多遇地震if(T>0&&T<0.1)st=0.45*amax/eit2+(amax-0.45*amax/eit2)/0.1;elseif T>0.1&&T<Tg st=amax;elseif T>Tg&&T<5*Tg st=(Tg/T)Ar*amax;elsest=(0.2Ar-eit1/eit2*(T-5*Tg)*amax;endp=0.9;%P為反響超越概率,一般取 0.85<P<1 sw=ksi/(pi*w)*st