




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、作业1绘制1940 El Centro,N-S分量地震动的绝对加速度、相对速度和相对位移反应谱。地震动:在PEER Ground Motion Database自行下载经典的1940 El Centro,N-S分量。 要求:在此模板内完成,A4纸打印。自编程序与软件(Bispec或Seismosigna等)计算反应谱进行对比。提交自编写程序。Matlab程序:clearfid = fopen('E:EarthquakeEl centro.txt');Accelerate,count = fscanf(fid,'%g'); %count 读入的记录的量Accele
2、rate=9.8*Accelerate' %单位统一为 m和stime=0:0.02:(count-1)*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某阻尼时最大绝对加速度
3、MDis=zeros(3,length(TA);MVel=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(1-Damp*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(1-Damp*Damp)+c); A(1,2)=e_t*s/Da
4、mFrcy; A(2,1)=-Frcy*e_t*s/sqrt(1-Damp*Damp); A(2,2)=e_t*(-s*Damp/sqrt(1-Damp*Damp)+c); d_f=(2*Damp2-1)/(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)*(c-D
5、amp/sqrt(1-Damp2)*s)-(2*d_3t+1/Frcy2)*(DamFrcy*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:(count-1) 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)*Vel
6、ocity(i)+B(2,1)*Accelerate(i)+B(2,2)*Accelerate(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
7、); AbsAcce=zeros(1,count); t=t+1; end j=j+1;endclose 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')gridfigur
8、e %绘制速度反应谱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')gridfigure %绘制绝对加速度反应谱plot(TA,MAcc(1,:),'-b',TA,MAcc(2,:),'-r
9、9;,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计算相对速度反应谱(下)的比较绝对加
10、速度反应谱(左图)与Seismosigna计算绝对加反应谱(右图)的比较作业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
11、=fai'*k*fai; %等效质量矩阵及等效刚度矩阵p0=0;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)*(P0-kk*q0);t=0.28;q(:,1)=
12、q0-t*q01+t2*q02/2; %q(:,1)表示q-1,即初始时刻前一时刻的位移q(:,2)=0;0; %q(:,2)表示初始时刻的位移Kh=mm/(t2);a=mm/(t2);b=kk-2*mm/(t2);P=fai'*p0;for i=2:7; Ph(:,i)=P-a*q(:,i-1)-b*q(:,i); q(:,i+1)=Kh(-1)*Ph(:,i); q1(:,i)=(q(:,i+1)-q(:,i-1)/(2*t); q2(:,i)=(q(:,i+1)-2*q(:,i)+q(:,i-1)/(t2);endu=fai*q;u1=fai*q1;u2=fai*q2;for i
13、=1:6 U(:,i)=u(:,i+1); U1(:,i)=u1(:,i+1); U2(:,i)=u2(:,i+1);endUU1U2%Newmark法(采用线加速度法)gamma=1/2;beta=1/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)*detaP
14、h(:,i); detaq1(:,i)=gamma*detaq(:,i)/(beta*t)-gamma*q11(:,i)/beta+t*(1-gamma/(2*beta)*q12(:,i); detaq2(:,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
15、=1:6 V(:,i)=v(:,i); V1(:,i)=v1(:,i); V2(:,i)=v2(:,i);endVV1V2 表1中心差分法Time0.28s0.56s0.84s1.12s1.40s位移(m)00.03070.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.4320402812-1.0001-5.6022-8.0268表2
16、Newmark法Time0.28s0.56s0.84s1.12s1.40s位移(m)0.00470.04440.18260.48500.97800.37261.38092.73174.04474.9744速度(m/s)0.05020.27510.75551.42382.07202.59264.42444.98504.17712.3360加速度(m/s2)0.35861.24762.18392.58962.04038.51884.5654-0.5615-5.2088-7.9416作业3、如图1所示一3层钢筋混凝土框架结构FR3,各层层高均为4.2m。设计荷载:楼面恒荷5.0kN/m2,楼面活荷载
17、2.0 kN/m2 ,屋面恒荷载5.0 kN/m2,屋面活荷载2.0 kN/m2,结构构件尺寸和配筋见表1,该结构遭受一地震动作用。使用非线性静力分析方法或时程分析方法求取该地震动作用下此结构(图1右图FR3-A轴平面框架。不考虑楼板贡献。)的顶点最大侧向位移值(参考使用的软件:OpenSees、Abaqus、Ansys、Sap2000、Idarc等)。(20分左右)地震动:1940 El Centro,N-S分量。 要 求:在此模板内完成,A4纸打印。 图1 3层结构的平立面布置图表1a 框架结构柱的截面尺寸及配筋结构编号材料尺寸( mm×mm)主筋面积( mm2)/箍筋混凝土钢筋
18、(主筋)边柱中柱边柱中柱FR3 13层 C30 HRB400400×400 400×400 1608/81001608/8100表1b 框架结构梁的截面尺寸及配筋结构编号材料尺寸(宽×高)(mm×mm)主筋面积( mm2)/箍筋混凝土钢筋(主筋)边跨中跨跨中支座跨中支座FR3 13层 C30 HRB400250×500 942/ 82001140/ 8100942/ 82001140/ 8100利用SAP2000建立模型如下图所示,图中的grid point是所要求的关键点。输入结构模型的物理数据,得到在El centro地震作用下框架侧移变形,如下图。单位:m表3 SAP2000输出的框架关键点位移表TABLE: Joint Displacem
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 【正版授权】 ISO 13140:2025 EN Electronic fee collection - Evaluation of on-board and roadside equipment for conformity to ISO 13141
- 【白城】2025年吉林白城市洮北区面向下半年应征入伍高校毕业生招聘事业单位工作人员17人笔试历年典型考题及考点剖析附带答案详解
- 【乐山】2025年上半年四川峨眉山市赴四川师范大学招聘事业单位工作人员15人笔试历年典型考题及考点剖析附带答案详解
- 【湘潭】2025年湖南湘潭市岳塘区事业单位引进急需紧缺专业人才55人笔试历年典型考题及考点剖析附带答案详解
- 古迹作文教学课件
- 2025年绿色规划视角下低碳城市建设实践案例分析报告
- 2025年共享出行平台信用积分兑换机制与市场响应策略研究报告
- 跨境电商保税仓仓储资源整合与共享模式研究报告
- 绿色金融产品创新与2025年绿色金融科技应用报告
- 企业合同履约管理办法
- 企业人力资源管理创新与发展趋势分析
- 2022城市轨道交通列车驾驶员技能及素质要求第1部分:地铁、轻轨和单轨
- 《人工智能基础与应用(第2版)》全套教学课件
- 2025年宿迁市公需考试试题
- 建设工程项目的组织协调保障措施
- 2025年春季学期国家开放大学行政管理专业毕业论文选题指南
- CO2还原合成甲醇反应器研究进展
- DB3713T-291-2023工程建设项目“多测合一”技术规程
- 苏教版四年级下册数学计算题每日一练带答案(共20天)
- 江苏连云港某公司“12.9”爆炸事故报告
- 人教版三年级上下数学试卷合集-综合素质训练
评论
0/150
提交评论