《地震工程作业.doc》由会员分享，可在线阅读，更多相关《地震工程作业.doc（10页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、作业1绘制1940 El Centro,NS分量地震动得绝对加速度、相对速度与相对位移反应谱。地震动:在PEER Ground Motion Database自行下载经典得1940 El Centro,NS分量。 要求:在此模板内完成,A4纸打印。自编程序与软件(Bispec或Seismosigna等)计算反应谱进行对比。提交自编写程序。Matlab程序:clearfid = fopen(E:EarthquakeEl centro、txt);Accelerate,count = fscanf(fid,%g); %count 读入得记录得量Accelerate=9、8*Accelerate; %
2、单位统一为 m与stime=0:0、02:(count1)*0、02; %单位 sDisplace=zeros(1,count); %相对位移Velocity=zeros(1,count); %相对速度AbsAcce=zeros(1,count); %绝对加速度DampA=0、00,0、02,0、05; %三个阻尼比TA=0、0:0、02:4; %TA=0、000001:0、02:4;结构周期Dt=0、02; %地震记录得步长%记录计算得到得反应,MDis为某阻尼时最大相对位移%MVel为某阻尼时最大相对速度,MAcc某阻尼时最大绝对加速度MDis=zeros(3,length(TA);MVe
3、l=zeros(3,length(TA);MAcc=zeros(3,length(TA);j=1; for Damp=0、00,0、05,0、1 t=1; for T=0、0:0、02:4 Frcy=2*pi/T ; DamFrcy=Frcy*sqrt(1Damp*Damp); e_t=exp(Damp*Frcy*Dt); s=sin(DamFrcy*Dt); c=cos(DamFrcy*Dt); A=zeros(2,2); A(1,1)=e_t*(s*Damp/sqrt(1Damp*Damp)+c); A(1,2)=e_t*s/DamFrcy; A(2,1)=Frcy*e_t*s/sqrt(
4、1Damp*Damp); A(2,2)=e_t*(s*Damp/sqrt(1Damp*Damp)+c); d_f=(2*Damp21)/(Frcy2*Dt); d_3t=Damp/(Frcy3*Dt); B=zeros(2,2); B(1,1)=e_t*(d_f+Damp/Frcy)*s/DamFrcy+(2*d_3t+1/Frcy2)*c)2*d_3t; B(1,2)=e_t*(d_f*s/DamFrcy+2*d_3t*c)1/Frcy2+2*d_3t; B(2,1)=e_t*(d_f+Damp/Frcy)*(cDamp/sqrt(1Damp2)*s)(2*d_3t+1/Frcy2)*(Da
5、mFrcy*s+Damp*Frcy*c)+1/(Frcy2*Dt); B(2,2)=e_t*(1/(Frcy2*Dt)*c+s*Damp/(Frcy*DamFrcy*Dt)1/(Frcy2*Dt); for i=1:(count1) Displace(i+1)=A(1,1)*Displace(i)+A(1,2)*Velocity(i)+B(1,1)*Accelerate(i)+B(1,2)*Accelerate(i+1); Velocity(i+1)=A(2,1)*Displace(i)+A(2,2)*Velocity(i)+B(2,1)*Accelerate(i)+B(2,2)*Accele
6、rate(i+1); AbsAcce(i+1)=2*Damp*Frcy*Velocity(i+1)Frcy2*Displace(i+1); end MDis(j,t)=max(abs(Displace); MVel(j,t)=max(abs(Velocity); if T=0、0 MAcc(j,t)=max(abs(Accelerate);else MAcc(j,t)=max(abs(AbsAcce); end Displace=zeros(1,count); Velocity=zeros(1,count); AbsAcce=zeros(1,count); t=t+1; end j=j+1;e
7、ndclose allfigure %绘制位移反应谱plot(TA,MDis(1,:),b,TA,MDis(2,:),r,TA,MDis(3,:),:k)title(Displacement)xlabel(Tn(s)ylabel(Displacement(m)legend(=0,=0、02,=0、05)gridfigure %绘制速度反应谱plot(TA,MVel(1,:),b,TA,MVel(2,:),r,TA,MVel(3,:),:k)title(Velocity)xlabel(Tn(s)ylabel(velocity(m/s)legend(=0,=0、02,=0、05)gridfigur
8、e %绘制绝对加速度反应谱plot(TA,MAcc(1,:),b,TA,MAcc(2,:),r,TA,MAcc(3,:),:k)title(Absolute Acceleration)xlabel(Tn(s)ylabel(absolute acceleration(m/s2)legend(=0,=0、02,=0、05)grid由此得到下图反应谱与由Seismosigna得到得反应谱作出比较相对位移反应谱(上)与Seismosigna计算相对位移反应谱(下)得比较相对速度反应谱(上)与Seismosigna计算相对速度反应谱(下)得比较绝对加速度反应谱(左图)与Seismosigna计算绝对加反
9、应谱(右图)得比较作业2、使用中心差分方法与Newmark方法分别计算下面单自由度结构得地震反应。题:考虑下面得具有2个自由度得一个简单系统,分析时间步长t=0、28sec得情况,求每时刻得加速度、速度、位移 求:要 求:在此模板内完成,A4纸打印。MATLAB程序:clear%求自振圆频率与振型m=2 0;0 1;k=6 2;2 4;v,d=eig(inv(m)*k);w=sqrt(d);fai1=v(:,1)、/v(1,1);fai2=v(:,2)、/v(1,2);fai=fai1 fai2;mm=fai*m*fai;kk=fai*k*fai; %等效质量矩阵及等效刚度矩阵p0=0;10;
10、%中心差分法s0=0;0;u0=0;0; %初速度及初位移均为0%以下过程均采用q、q1、q2得形式分别表示位移、速度、加速度q0=fai1*m*s0、/(fai1*m*fai1);fai2*m*s0、/(fai2*m*fai2);q01=fai1*m*u0、/(fai1*m*fai1);fai2*m*u0、/(fai2*m*fai2);P0=fai*p0;q02=mm(1)*(P0kk*q0);t=0、28;q(:,1)=q0t*q01+t2*q02/2; %q(:,1)表示q1,即初始时刻前一时刻得位移q(:,2)=0;0; %q(:,2)表示初始时刻得位移Kh=mm/(t2);a=mm/
11、(t2);b=kk2*mm/(t2);P=fai*p0;for i=2:7; Ph(:,i)=Pa*q(:,i1)b*q(:,i); q(:,i+1)=Kh(1)*Ph(:,i); q1(:,i)=(q(:,i+1)q(:,i1)/(2*t); q2(:,i)=(q(:,i+1)2*q(:,i)+q(:,i1)/(t2);endu=fai*q;u1=fai*q1;u2=fai*q2;for i=1:6 U(:,i)=u(:,i+1); U1(:,i)=u1(:,i+1); U2(:,i)=u2(:,i+1);endUU1U2%Newmark法(采用线加速度法)gamma=1/2;beta=1/
12、6;Kh1=kk+mm/(beta*t2);a1=mm/(beta*t);b1=mm/(2*beta);detaP=0 0 0 0 0;0 0 0 0 0;q1(:,1)=0;0;q11(:,1)=0;0;q12(:,1)=q02;for i=1:5 detaPh(:,i)=detaP(:,i)+a1*q11(:,i)+b1*q12(:,i); detaq(:,i)=Kh1(1)*detaPh(:,i); detaq1(:,i)=gamma*detaq(:,i)/(beta*t)gamma*q11(:,i)/beta+t*(1gamma/(2*beta)*q12(:,i); detaq2(:,
13、i)=detaq(:,i)/(beta*t2)q11(:,i)/(beta*t)q12(:,i)/(2*beta); q1(:,i+1)=q1(:,i)+detaq(:,i); q11(:,i+1)=q11(:,i)+detaq1(:,i); q12(:,i+1)=q12(:,i)+detaq2(:,i);endv=fai*q1;v1=fai*q11;v2=fai*q12;for i=1:6 V(:,i)=v(:,i); V1(:,i)=v1(:,i); V2(:,i)=v2(:,i);endVV1V2 表1中心差分法Time0、28s0、56s0、84s1、12s1、40s位移(m)00、0
14、3070、16750、48711、01700、39201、44512、83384、14415、0152速度(m/s)0、05490、29920、81491、51692、16752、58054、36034、81973、89541、9873加速度(m/s2)0、39201、35292、33122、68281、96428、43204028121、00015、60228、0268表2 Newmark法Time0、28s0、56s0、84s1、12s1、40s位移(m)0、00470、04440、18260、48500、97800、37261、38092、73174、04474、9744速度(m/s)0
15、、05020、27510、75551、42382、07202、59264、42444、98504、17712、3360加速度(m/s2)0、35861、24762、18392、58962、04038、51884、56540、56155、20887、9416作业3、如图1所示一3层钢筋混凝土框架结构FR3,各层层高均为4、2m。设计荷载:楼面恒荷5、0kN/m2,楼面活荷载2、0 kN/m2 ,屋面恒荷载5、0 kN/m2,屋面活荷载2、0 kN/m2,结构构件尺寸与配筋见表1,该结构遭受一地震动作用。使用非线性静力分析方法或时程分析方法求取该地震动作用下此结构(图1右图FR3A轴平面框架。不考
16、虑楼板贡献。)得顶点最大侧向位移值(参考使用得软件:OpenSees、Abaqus、Ansys、Sap2000、Idarc等)。(20分左右)地震动:1940 El Centro,NS分量。 要 求:在此模板内完成,A4纸打印。 图1 3层结构得平立面布置图表1a 框架结构柱得截面尺寸及配筋结构编号材料尺寸( mmmm)主筋面积( mm2)/箍筋混凝土钢筋(主筋)边柱中柱边柱中柱FR3 13层 C30 HRB400400400 400400 1608/81001608/8100表1b 框架结构梁得截面尺寸及配筋结构编号材料尺寸(宽高)(mmmm)主筋面积( mm2)/箍筋混凝土钢筋(主筋)边跨
17、中跨跨中支座跨中支座FR3 13层 C30 HRB400250500 942/ 82001140/ 8100942/ 82001140/ 8100利用SAP2000建立模型如下图所示,图中得grid point就是所要求得关键点。输入结构模型得物理数据,得到在El centro地震作用下框架侧移变形,如下图。单位:m表3 SAP2000输出得框架关键点位移表TABLE: Joint DisplacementsJointOutputCaseCaseTypeStepTypeU1U2U3TextTextTextTextmmm1POSITEbinationMax0001POSITEbinationMin00016POSITEbinationMax0、1931690、0001010、00060616POSITEbinationMin0、2026170、0001680、001481 表中U1、U2、U3分别表示所建模型得整体坐标系得x轴、y轴、z轴得正方向。由上述可得,所求FRA3框架得顶点最大侧移量就是0、2026m。