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

人工地震波依据三角级数法武汉@桥梁隧道 799084759采用shinozuka的方法来模拟平稳化后的随机地面运动加速度其中：强度包线：P——反应超越概率。S——功率谱密度函数。?ω——频谱分度（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;%Δwag=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 (s))ylabel(ag (m/s2))subplot(222)plot(t,ag1)title(功率谱密度函数（随机后）)xlabel(t=0:30 (s))ylabel(∑Ck*cos() (m/s2))subplot(223)plot(t,ag2)title(强度包线)xlabel(t=0:30 (s))ylabel(ξ (m/s2))%强度包线子程序function ksit=ft(t)t0=2;tn=10;c=0.2;if t=0t=t0 ksit=(t/t0)^2; elseif tt0ttnksit=1;elseksit=exp(-c*(t-tn));endend%计算功率密度函数子程序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;%七度区多遇地震if(T0T0.1)st=0.45*amax/eit2+(amax-0.45*amax/eit2)/0.1;elseif T0.1TTgst=amax;elseif TTgT5*Tgst=(Tg/T)^r*amax;elsest=(0.2^r-eit1/eit2*(T-5*Tg))*amax; endp=0.9;%P为反应超越概率，一般取0.85P1 sw=ksi/(pi*w)*st^2/log(-pi/(w*T)*log(1-p));end