1、 式(1) 转化成状态方程为: 式(2) 其中,, 在Matlab中用函数y=lsim(A,B,C,D,u,t),即可求得系统的状态量。3. 地震反应谱计算方法上面分析中的分别为结构相对地面的位移、速度和加速度。绝对加速度反应谱:由式(1)可得结构的绝对加速度为:其绝对加速度反应谱值为:绝对加速度反应放大系数为;伪加速度反应谱:按照抗震设计规范中的结构剪力公式:得其伪加速度反应谱值为: 伪加速度反应放大系数为: 近似加速度反应谱():当激励为简谐荷载时:求两次导得到加速度反应:所以绝对加速度谱与伪加速度谱具有如下近似关系:4. 计算结果 根据前面分析,分别在阻尼比取为时计算得到的EI Cent
2、ro (1940,NS)地震波的绝对加速度反应谱、伪加速度反应谱与近似加速度反应谱的结果如图3所示。a) b) 图3 绝对加速度反应谱、伪加速度反应谱与近似加速度反应谱的比较为了更清楚的对比三者的差别,让三者两两相除(近似加速度反应谱与伪加速度反应谱比值恒定故未绘图),并将各比值绘于图4-7中。图4 伪加速度反应谱与绝对加速度反应谱之比()图5 近似加速度反应谱与绝对加速度反应谱之比(图6 伪加速度反应谱与绝对加速度反应谱之比(图7 近似加速度反应谱与绝对加速度反应谱之比(5. 结果分析 图像单独分析由图3可以看出,在结构阻尼比较小(一般在之间)时,三条曲线几近重合即伪加速度反应谱、绝对加速度
3、反应谱及近似加速度反应谱差别很小(图3-a),但在结构阻尼比较大时(图3-b),三者有较大的差别。因此,对一般结构而言,绝对加速度反应谱和伪加速度反应谱差别不大。由图4可以看出,当阻尼比为5%时,伪加速度反应谱与绝对加速度反应谱的比值绝大部分都小于1,但不是全部。说明绝对加速度反应谱在大部分情况下都比伪加速度反应谱值要大,当振动周期比较长时这个特点体现更明显。由图5可以看出,当阻尼比为5%时,近似加速度谱总体峰值高于伪加速度反应谱,其与绝对加速度反应谱的比值在1附近跳动,当结构周期比较长时这个比值几乎都小于1。由图6可以看出,当阻尼比为30%时,伪加速度反应谱的值都小于绝对加速度反应谱, 且二
4、者的差值比阻尼比为5%时大很多。这也证明了图3中的现象。由图7可以看出,当阻尼比为30%时,在长周期段(T)近似加速度反应谱的值几乎都小于绝对加速度反应谱,在短周期段(T)近似加速度反应谱的值几乎都大于绝对加速度反应谱。 结果总体分析由所有的结果总体来看有如下特点:三种谱总体差别不大,其中伪加速度反应谱与绝对加速度反应谱更为接近。无论阻尼比为多少,近似加速度反应谱和伪加速度反应谱二者与绝对加速度反应谱的比值总体上的趋势都是随着结构振动周期的增大而减小,阻尼比越大这个减小的幅度就越大。结构振动周期越长近似加速度反应谱和伪加速度反应谱二者与绝对加速度反应谱偏离的越远。附:MATLAB编程代码cle
5、ar; clc; close all;load ;t=0;elcen(:,1);xg_gal=0;,2);%galxg=xg_gal/100;%m/m2figure(4)plot(t,xg,linewidth,2)xlabel(fontsize16fontnameTimes New Roman (sec)ylabel(fontsize16fontnameTimes New Roman a_g (m/s2)set(gca,fontsize,16,fontname,Times New Romangrid onkesi=; % Damping ratioT=:5; %periodsSaa=T; %A
6、BS acceleration-response-spectrumSap=T; %pseudo acceleration-response-spectrumSae=T; %Sae=sqrt(1+4*kesi2)*Sapfor i=1:length(T) omiga=2*pi/T(i); A=0 1;-omiga2 -2*kesi*omiga; B=0;-1; C=eye(2); D=zeros(2,1); y=lsim(A,B,C,D,xg,t);%ABS acceleration-response-spectrum Saa(i)=max(abs(2*kesi*omiga*y(:,2)+omi
7、ga2*y(:,1)/max(abs(xg);%pseudo acceleration-response-spectrum Sap(i)=max(abs(y(:,1)*omiga2/max(abs(xg); %sqrt(1+4*kesi2)*Sap(i) Sae(i)=sqrt(1+4*kesi2)*Sap(i);endfigure(1)semilogx(T,Saa,k-,T,Sap,k:,T,Sae,k-fontsize16fontnameTimes New Roman (m/s2)legend(% set(gca,xtick,0,1,10,50)figure(2)fontsize16fontnameTimes New Roman xprimeprime_ao & S_APytick,0:3)% subplot(2,1,1),figure(3)semilogx(T,Sap./Saa,fontsize16fontnameTimes New Roman %set(gca,:,0.:1)% subplot(2,1,2),semilogx(T,Sae./Saa,figure(5)semilogx(T,Sae./Sap,
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1