地震工程学-反应谱和地震时程波的相互转化matlab编程doc资料

精品文档

图2 位移反应谱

图3 速度反应谱

图4 加速度反应谱

精品文档

精品文档

一、 反应谱生成地震波

1 所取的反应谱为上海市设计反应谱

图5 上海市设计反应谱

2 反应谱取值程序为：

%%规范反应谱取值程序 参照01年抗震规范

function rs_z=r_s_1(pl,zn,ld,cd,fz) %%%pl 圆频率,zn阻尼比,ld烈度,cd场地类型,场地分组fz %%%%烈度选择 if ld==6 arfmax=0.11; end if ld==7 arfmax=0.23; end if ld==8 arfmax=0.45; end if ld==9 arfmax=0.90; end

%%%%场地类别，设计地震分组选择 if cd==1 if fz==1 Tg=0.25; end

if fz==2 Tg=0.30; end if fz==3 Tg=0.35; end end if cd==2 if fz==1 Tg=0.35;

精品文档

精品文档

end

if fz==2 Tg=0.40; end if fz==3 Tg=0.45; end end if cd==3 if fz==1 Tg=0.45; end

if fz==2 Tg=0.55; end if fz==3 Tg=0.65; end end if cd==4 if fz==1 Tg=0.65; end

if fz==2 Tg=0.75; end if fz==3 Tg=0.90; end end %%%%%%%%%

ceita=zn; %%%%%阻尼比 lmt1=0.02+(0.05-ceita)/8; if lmt1<0 lmt1=0; end

lmt2=1+(0.05-ceita)/(0.06+1.7*ceita); if lmt2<0.55 lmt2=0.55; end

sjzs=0.9+(0.05-ceita)/(0.5+5*ceita); %%%%%分段位置 T1 T2 T3 T1=0.1; T2=Tg;

精品文档