预览加载中,请您耐心等待几秒...
1/9
2/9
3/9
4/9
5/9
6/9
7/9
8/9
9/9

在线预览结束,喜欢就下载吧,查找使用更方便

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

《结构隔震与耗能减振》 课程作业1 姓名:XXX 学号:XXXXXX 作业要求 计算并绘制5%和30%阻尼比下的EICentro(1940,NS)地震记录的的绝对加速度反应谱、伪加速度反应谱及(此次作业中暂且称之为“近似加速度反应谱”),并比较三者的异同。 建立分析模型 对于图1所示的EICentro(1940,NS)地震波时程,采用图2所示的结构进行绝对加速度反应谱和伪加速度反应谱分析。 图1EICentro(1940,NS)波时程图2分析模型 设系统的自振频率为,阻尼比为,则此SDOF结构的标准运动方程为: 式(1) 转化成状态方程为: 式(2) 其中,,,, ,, 在Matlab中用函数y=lsim(A,B,C,D,u,t),即可求得系统的状态量。 3.地震反应谱计算方法 上面分析中的分别为结构相对地面的位移、速度和加速度。 绝对加速度反应谱: 由式(1)可得结构的绝对加速度为: 其绝对加速度反应谱值为: 绝对加速度反应放大系数为; 伪加速度反应谱: 按照抗震设计规范中的结构剪力公式: 得其伪加速度反应谱值为: 伪加速度反应放大系数为: 近似加速度反应谱(): 当激励为简谐荷载时: 求两次导得到加速度反应: 所以绝对加速度谱与伪加速度谱具有如下近似关系: 4.计算结果 根据前面分析,分别在阻尼比取为时计算得到的EICentro(1940,NS)地震波的绝对加速度反应谱、伪加速度反应谱与近似加速度反应谱的结果如图3所示。 a)b)图3绝对加速度反应谱、伪加速度反应谱与近似加速度反应谱的比较 为了更清楚的对比三者的差别,让三者两两相除(近似加速度反应谱与伪加速度反应谱比值恒定故未绘图),并将各比值绘于图4-7中。 图4伪加速度反应谱与绝对加速度反应谱之比() 图5近似加速度反应谱与绝对加速度反应谱之比() 图6伪加速度反应谱与绝对加速度反应谱之比() 图7近似加速度反应谱与绝对加速度反应谱之比()5.结果分析 5.1图像单独分析 由图3可以看出,在结构阻尼比较小(一般在0.01~0.05之间)时,三条曲线几近重合即伪加速度反应谱、绝对加速度反应谱及近似加速度反应谱差别很小(图3-a),但在结构阻尼比较大时(图3-b),三者有较大的差别。因此,对一般结构而言,绝对加速度反应谱和伪加速度反应谱差别不大。 由图4可以看出,当阻尼比为5%时,伪加速度反应谱与绝对加速度反应谱的比值绝大部分都小于1,但不是全部。说明绝对加速度反应谱在大部分情况下都比伪加速度反应谱值要大,当振动周期比较长时这个特点体现更明显。 由图5可以看出,当阻尼比为5%时,近似加速度谱总体峰值高于伪加速度反应谱,其与绝对加速度反应谱的比值在1附近跳动,当结构周期比较长时这个比值几乎都小于1。 由图6可以看出,当阻尼比为30%时,伪加速度反应谱的值都小于绝对加速度反应谱,且二者的差值比阻尼比为5%时大很多。这也证明了图3中的现象。 由图7可以看出,当阻尼比为30%时,在长周期段(T>0.08s)近似加速度反应谱的值几乎都小于绝对加速度反应谱,在短周期段(T<0.08s)近似加速度反应谱的值几乎都大于绝对加速度反应谱。 5.2结果总体分析 由所有的结果总体来看有如下特点: eq\o\ac(○,1)三种谱总体差别不大,其中伪加速度反应谱与绝对加速度反应谱更为接近。 eq\o\ac(○,2)无论阻尼比为多少,近似加速度反应谱和伪加速度反应谱二者与绝对加速度反应谱的比值总体上的趋势都是随着结构振动周期的增大而减小,阻尼比越大这个减小的幅度就越大。 eq\o\ac(○,3)结构振动周期越长近似加速度反应谱和伪加速度反应谱二者与绝对加速度反应谱偏离的越远。 附:MATLAB编程代码 clear;clc;closeall; loadelcen.dat; t=[0;elcen(:,1)]; xg_gal=[0;elcen(:,2)];%gal xg=xg_gal/100;%m/m^2 figure(4) plot(t,xg,'linewidth',2) xlabel('\fontsize{16}\fontname{TimesNewRoman}½á¹¹ÖÜÆÚ(sec)') ylabel('\fontsize{16}\fontname{TimesNewRoman}a_g(m/s^2)') set(gca,'fontsize',16,'fontname','TimesNewRoman') gridon kesi=0.3;%Dampingratio T=[0.01:0.01:5];%periods Saa=T;%