MATALAB 作业
某三层钢筋混凝土结构,结构的各层特性参数为:第一层到第三层质量m 分别为2400kg ,1200kg ,1200kg ,第一层到第三层刚度k 分别为3.3*10^4N/m,1.1*10^4N/m,0.66^4N/m.。地震采用acc_ElCentro_0.34g ,采样周期为0.02。
M3=1200kg
K3=0.66*10^4N/m.
M2=1200kg
K2=1.1*10^4N/m
M1=2400kg
K1=3.3*10^4N/m
用振型分解法求解结构地震反应的MATLAB 层序如下,编制该程序的程序框图以下所示
%振型分解法求解结构地震反应;主程序 clear
开始 输入地震参数和结构参数 计算结构振型与自振型频率 计算振型参与系数 计算单自由度体系的地震反应
求解结构的地震反应
输出结果
结束
clc
%地震波数据
xs=2*0.287;
dzhbo=load('acc_ElCentro_0.34g_0.02s.txt');
ag=dzhbo*0.01*xs;
dt=0.02;
ndzh=400;
cn=3; %cn为结构的层数,即质点数
m0=[2.4 1.2 1.2]*1e+3; %结构各层质量
k0=[3.3 1.1 0.66]*1e+5; %结构各层刚度
l=diag(ones(cn));
m=diag(m0); %计算质量矩阵
[ik]=matrixju(k0,cn); %计算刚度矩阵
[x,d]=eig(ik,m); %结构动力特性求解
d=diag(sqrt(d)); %求解结构圆频率
for i=1:cn;
[d1(i),j]=min(d);
xgd(:,i)=x(:,j);
d(j)=max(d)+1;
end %以此循环对所求频率和振型进行排序w=d1; %所求自振频率
x=xgd; %所求结构主振型
a1=2*w(1)*w(2)*(0.05*w(2)-0.07*w(1))/(w(2)^2-w(1)^2); a2=2*(0.07*w(2)-0.05*w(1))/(w(2)^2-w(1)^2);
for j=1:cn
x(:,j)=x(:,j)/x(cn,j);
znb0(j)=(a1+a2*w(j)^2)/2/w(j);
zhcan(j)=(x(:,j))'*m*l/((x(:,j))'*m*x(:,j));
%求解振型参数
[dlt(j,:),dltacceler(j,:)]=zxzj(znb0(j),w(j),ag);
end
%求解结构各层的地震反应
for i=1:cn;
disp1=0;
accel1=0;
for j=1:cn
disp0=zhcan(j)*dlt(j,:)*x(i,j);
accel0=zhcan(j)*dltacceler(j,:)*x(i,j);
disp1=disp1+disp0;
accel1=accel1+accel0;
end
disp(i,:)=disp1;
accel(i,:)=accel1;
end
t=0:dt:ndzh*dt;
%subplot(2,2,1)
%plot(t,disp(3,:)*1e+3,'k-')
%subplot(2,2,2)
plot(t,accel(3,:),'k-')
%子程序
%用于求解单自由度结构体系的地震反应
function [bx,acceler]=zxzj(znb,w,dag)
dt=0.02;
n=400;
x(1)=0;
dx(1)=0;
ddx(1)=0;
s=1+znb*dt*w+w^2*dt^2/6; %中间参数s
for i=1:n
a(i)=x(i)+dx(i)*dt+ddx(i)*dt^2/3;
b(i)=dx(i)+ddx(i)*dt/2;
ddx(i+1)=-1*(dag(i+1)+1*znb*w*b(i)+w^2*a(i))/s;%加速度dx(i+1)=b(i)+ddx(i+1)*dt/2;%速度
x(i+1)=a(i)+ddx(i+1)*dt^2/6;%位移
end
bx=x;
acceler=ddx;
%子程序
%刚度和阻尼矩阵的聚合
function [kcju]=matrixju(korc,cn)
kcju=zeros(cn);
for i=1:cn-1
kcju(i,i)=korc(i)+korc(i+1);
kcju(i,i+1)=-korc(i+1);
kcju(i+1,i)=-korc(i+1);
end
kcju(cn,cn)=korc(cn);
经程序求解,该结构的自振频率为
W= 4.8683 11.0666 15.4541
结构的振型矩阵为
X=
0.1634 -0.9238 2.7604
0.5691 -1.2267 -3.3423
1.0000 1.0000 1.0000
移位/ mm
时间/s
(a) 顶层位移反应
加速度/(m/s^2)
时间/s
(b)顶层位移反应
第一主振型(低)
{}
(1)0.1630.5691????=??????Y 1 0.569 0.163
1
1.227
0.924 第二主振型 {}(2)0.9241.2271-????=-??????Y
1
3.3
42
2.
7
6
{}(3) 2.7603.3421????=-??????Y 第三主振型(高)
3.2.2(忽略剪重比验算) 解:(1)由∑∑===n i ji i n i ji i m m 121φφγ 得: 363.1024.33027.45118667.027334.027118667.027334.0272221==?+?+??+?+?= γ 428.0988.41991.17118)666.0(27)667.0(27118)666.0(27)667.0(272 222-=-=?+-?+-??+-?+-?=γ 063.0817.702568.44118)035.3(27)019.4(27118)035.3(27)019.4(272223==?+-?+??+-?+?= γ (2)查表3-2,3-3得35.0=g T ,16.0max =α,则 123.016.0467.035.09.0max 9.011=???? ??=???? ??=ααT T g 16.0max 2==αα
16.0max 3==αα (3)由ji j j i ji G F φγα=得 第一振型各质点水平地震作用为: kN F 16.148334.0363.1123.08.927011=????= kN F 88.295667.0363.1123.08.927012=????= kN F 73.2951363.1123.08.918013=????= 第二振型各质点水平地震作用为: kN F 86.120)667.0()428.0(16.08.927021=-?-???= kN F 68.120)666.0()428.0(16.08.927022=-?-???= kN F 80.1201)428.0(16.08.918023-=?-???= 第三振型各质点水平地震作用为: kN F 19.107019.4063.016.08.927031=????= kN F 95.80)035.3(063.016.08.927032-=-????= kN F 78.171063.016.08.918033=????= (4)由各振型水平地震作用产生的底部剪力为: kN F F F V 77.73973.29588.29516.14813121111=++=++= kN F F F V 74.12080.12068.12086.12023222121=-+=++= kN F F F V 02.4478.1795.8019.10733323131=+-=++= (5)振型组合求最大底部剪力: kN V V j j 85.75002.4474.12077.73922231211=++==∑= (6)由各振型水平地震作用产生的结构顶层位移为:
一、判断题(共223小题) 1。结构的类型若按几何特征可分为平面结构和空间结构。(A) 2、狭义结构力学的研究对象是板、壳结构(B)。 3 单铰相当于两个约束。(A) 4、单刚节点相当于三个约束。(A) 5、静定结构可由静力平衡方程确定全部约束力和内力。A 6、超静定结构可由静力平衡方程确定全部约束力和内力B。 7 无多余约束的几何不变体系是静定结构。A 8 三刚片规则中三铰共线为可变体系。B 9 两刚片用一个单铰和一个不通过该铰的链杆组成的体系为静定结构。A 10 两刚片用一个单铰和一个不通过该铰的链杆组成的体系为超静定结构B。 11链杆相当于两个约束。B 12 平面上的自由点的自由度为2 A 13 平面上的自由刚体的自由度为3 A 14 铰结点的特征是所联结各杆可以绕结点中心自由转动。A 15 有多余约束的几何不变体系是超静定结构。A 16 无多余约束的几何可变体系是超静定结构。B 17、无多余约束的几何可变体系是静定结构。B 18刚结点的特征是当结构发生变形时汇交于该点的各杆端间相对转角为零。A 19 三刚片规则中三铰共线为瞬变体系。A 20三个本身无多余约束的刚片用三个不共线的单铰两两相连,则组成的体系为静定结构。A 21 一个刚结点相当于3个约束。 22 一个连接3个刚片的复铰相当于2个单铰。A 23 一个铰结三角形可以作为一个刚片。A 24 一个铰结平行四边形可以作为一个刚片。B 25 一根曲杆可以作为一个刚片。A 26 一个连接4个刚片的复铰相当于2个单铰.B 27 任意体系加上或减去二元体,改变体系原有几何组成性质。B 28 平面几何不变体系的计算自由度一定等于零。B 29 平面几何可变体系的计算自由度一定等于零。B 30 三刚片体系中若有1对平行链杆,其他2铰的连线与该对链杆不平行,则该体系为几何不变体系。A 31 三刚片体系中,若有三对平行链杆,那么该体系仍有可能是几何不变的。B 32 三刚片体系中,若有2对平行链杆,那么该体系仍有可能是几何不变的。A 33 一个单铰相当于一个约束。B 34 进行体系的几何组成分析时,若体系通过三根支座链杆与基础相连,可以只分析体系内部。B 35 三刚片体系中,若有两个虚铰在无穷远处,则该体系一定为几何可变。B 36 有多余约束的体系为静定结构。B 37 静定结构一定几何不变。A 38 超静定结构一定几何不变.A 39 几何不变体系一定是静定结构。B 40几何不变体系一定是超静定结构。B 41力是物体间相互的机械作用。A 42 力的合成遵循平行四边形法则。A 43 力的合成遵循三角形法则。A 44 力偶没有合力。A 45 力偶只能用力偶来平衡。A 46 力偶可以和一个力平衡。B 47 力偶对物体既有转动效应,又有移动效应。B 48 固定铰支座使结构在支承处不能移动也不能转动。B 49 可动铰支座使结构在支承处能够转动,但不能沿链杆方向移动。A 50 结点法求解桁架内力应按照结构几何组成相反顺序来求解。A 51 将一个已知力分解为两个力可得到无数解答。A 52 作用力和反作用力是作用在同一物体上的两个力。B 53 作用力和反作用力是作用在不同物体上的两个力。A 54 两个力在同一轴上的投影相等,此两力必相等 B 55 力偶对平面内任一点的矩等于力偶矩A 56 力偶在坐标轴上的投影的代数和等于零A 57 一个固定铰支座相当于两个约束。A 58三个本身无多余约束的刚片用三个不共线的单铰两两相连,则组成的体系为超静定结构B 59 桁架是“只受结点荷载作用的直杆、铰结体系”。A 60桁架结构的内力有轴力。A 61 拱的合理拱轴线均为二次抛物线。B 62无铰拱属于超静定结构。A 63 三铰刚架和三铰拱都属于推力结构。A 64 简支刚架属于推力结构。B 65 三铰拱属于静定结构。A 66 相同竖向载荷作用下,同跨度拱的弯矩比代梁的弯矩大得多。B 67 桁架结构中,杆的内力有轴力和剪力。B 68 竖向载荷作用下,简支梁不会产生水平支反力.A 69 竖向载荷作用下,拱不会产生水平支反力。B 70 竖向载荷作用下,拱的水平推力与拱高成正比。B
结构动力学学习总结
通过对本课程的学习,感受颇深。我谈一下自己对这门课的理解: 一.结构动力学的基本概念和研究内容 随着经济的飞速发展,工程界对结构系统进行动力分析的要求日益提高。我国是个多地震的国家,保证多荷载作用下结构的安全、经济适用,是我们结构工程专业人员的基本任务。结构动力学研究结构系统在动力荷载作用下的位移和应力的分析原理和计算方法。它是振动力学的理论和方法在一些复杂工程问题中的综合应用和发展,是以改善结构系统在动力环境中的安全和可靠性为目的的。高老师讲课认真负责,结合实例,提高了教学效率,也便于我们学生寻找事物的内在联系。这门课的主要内容包括运动方程的建立、单自
由度体系、多自由度体系、无限自由度体系的动力学问题、随机振动、结构抗震计算及结构动力学的前沿研究课题。既有线性系统的计算,又有非线性系统的计算;既有确定性荷载作用下结构动力影响的计算,又有随机荷载作用下结构动力影响的随机振动问题;阻尼理论既有粘性阻尼计算,又有滞变阻尼、摩擦阻尼的计算,对结构工程最为突出的地震影响。 二.动力分析及荷载计算 1.动力计算的特点 动力荷载或动荷载是指荷载的大小、方向和作用位置随时间而变化的荷载。如果从荷载本身性质来看,绝大多数实际荷载都应属于动荷载。但是,如果荷载随时间变化得很慢,荷载对结构产生的影响与
静荷载相比相差甚微,这种荷载计算下的结构计算问题仍可以简化为静荷载作用下的结构计算问题。如果荷载不仅随时间变化,而且变化很快,荷载对结构产生的影响与静荷载相比相差较大,这种荷载作用下的结构计算问题就属于动力计算问题。 荷载变化的快与慢是相对与结构的固有周期而言的,确定一种随时间变化的荷载是否为动荷载,须将其本身的特征和结构的动力特性结合起来考虑才能决定。 在结构动力计算中,由于荷载时时间的函数,结构的影响也应是时间的函数。另外,结构中的内力不仅要平衡动力荷载,而且要平衡由于结构的变形加速度所引起的惯性力。结构的动力方程中除了动力荷载和弹簧力之外,还要引入因其质量产生的惯性力和耗散能量的阻尼力。而
、选择题(每小题3分,共18分) 1?图示体系的几何组成为:() A.几何不变,无多余联系; B.几何不变,有多余联系; C.瞬 变; 2?静定结构在支座移动时,会产生:() A.内力; B.应力; C.刚体位移; D.变形 3?在径向均布荷() A.圆弧线; 载作用下, B .抛物线 铰拱的合理轴线为: C .悬链线;D.正弦曲线。 4?图示桁架的零A. 6; B. 7杆数目为: ; C. 8 ; ( ) D. 9 。 D.常变。
5?图a结构的最后弯矩图为:() A.图b ; B .图c;C .图d; D .都不对。 6?力法方程是沿基本未知量方向的:() A.力的平衡方程; B.位移为零方程; C.位移协调方程;D ?力的平衡及位移为零方程。 :■、填空题(每题3分,共9分) 1.从几何组成上讲,静定和超静定结构都是_______________________________ 体系, 前者__________ 多余约束而后者______________________ 多余约束。 2.图b是图a结构_______________ 截面的 ____________ 影响线。 彳、亡A 卜 1 B K D —i |i li 11 行)f- 3._________________________________________________ 图示结构AB杆B端的转动刚度为_________________________________________________ ,分配系数为________ , 传递系数为 ___________ 。 三、简答题(每题5分,共10分) 1.静定结构内力分析情况与杆件截面的几何性质、材料物理性质是否相关?为什么? 2.影响线横坐标和纵坐标的物理意义是什么?
《结构动力学》读书报告 学院 专业 学号 指导老师 2013 年 5月 28日
摘要:本书在介绍基本概念和基础理论的同时,也介绍了结构动力学领域的若干前沿研究课题。既注重读者对基本知识的掌握,也注重读者对结构振动领域研究发展方向的掌握。主要容包括运动方程的建立、单自由度体系、多自由度体系、无限自由度体系的动力学问题、随机振动、结构动力学的前沿研究课题。侧重介绍单自由度体系和多自由度体系,重点突出,同时也着重介绍了在抗震中的应用。 1 概述 1.1结构动力学的发展及其研究容: 结构动力学,作为一门课程也可称作振动力学,广泛地应用于工程领域的各个学科,诸如航天工程,航空工程,机械工程,能源工程,动力工程,交通工程,土木工程,工程力学等等。作为固体力学的一门主要分支学科,结构动力学起源于经典牛顿力学,就是牛顿质点力学。质点力学的基本问题是用牛顿第二定律来建立公式的。牛顿质点力学,拉格朗日力学和哈密尔顿力学是结构动力学基本理论体系组成的三大支柱。 经典动力学的理论体系早在19世纪中叶就已建立,。但和弹性力学类似,理论体系虽早已建立,但由于数学求解上的异常困难,能够用来解析求解的实际问题实在是少之又少,能够通过手算完成的也不过仅仅限于几个自由度的结构动力体系。因此,在很长一段时间,动力学的求解思想在工程实际中并未得到很好的应用,人们依然习惯于在静力学的畴用静力学的方法来解决工程实际问题。 随着汽车,飞机等新时代交通工具的出现,后工业革命时代各种大型机械的创造发明,以及越来越多的摩天大楼的拔地而起,工程界日新月异的发展和变化对工程师们提出了越来越高的要求,传统的只考虑静力荷载的设计理念和设计方法显然已经跟不上时代的要求了。也正是从这个时候起,结构动力学作为一门学科,也开始受到工程界越来越高的重视,从而带动了结构动力学的快速发展。 结构动力学这门学科在过去几十年来所经历的深刻变革,其主要原因也正是由于电子计算机的问世使得大型结构动力体系数值解的得到成为可能。由于电子计算机的超快速度的计算能力,使得在过去凭借手工根本无法求解的问题得到了解决。目前,由于广泛地应用了快速傅立叶变换(FFT),促使结构动力学分析发生了更加深刻地变化,而且使得结构动力学分析与结构动力试验之间的相互关系也开始得以沟通。总之,计算机革命带来了结构动力学求解方法的本质改变。 作为一门课程,结构动力学的基本体系和容主要包括以下几个部分:单自由度系统结构动力学,;多自由度系统结构动力学,;连续系统结构动力学。此外,如果系统上所施加的动力荷载是确定性的,该系统就称为确定性结构动力系统;而如果系统上所施加的动力荷载是非确定性的,该系统就称为概率性结构动力系统。 1.2主要理论分析 结构的质量是一连续的空间函数,因此结构的运动方程是一个含有空间坐标和时间的偏微分方程,只是对某些简单结构,这些方程才有可能直接求解。对于绝大多数实际结构,在工程分析中主要采用数值方法。作法是先把结构离散化成为一个具有有限自由度的数学模
振型分解反应谱法
振型分解反应谱法 振型分解反应谱法是用来计算多自由度体系地震作用的一种方法。该法是利用单自由度体系的加速度设计反应谱和振型分解的原理,求解各阶振型对应的等效地震作用,然后按照一定的组合原则对各阶振型的地震作用效应进行组合,从而得到多自由度体系的地震作用效应。振型分解反应谱法一般可考虑为计算两种类型的地震作用:不考虑扭转影响的水平地震作用和考虑平扭藕联效应的地震作用。 适用条件 (1)高度不超过40米,以剪切变形为主且质量和刚度沿高度分布比较均匀的结构,以及近似于单质点体系的结构,可采用底部剪力法计算。(此为底部剪力法的适用范围) (2)除上述结构以外的建筑结构,宜采用“振型分解反应谱法”。 (3)特别不规则的建筑、甲类建筑和规范规定的高层建筑,应采用时程分析法进行补充计算。 刚重比 刚重比是指结构的侧向刚度和重力荷载设计值之比,是影响重力二阶效应的主要参数 刚重比=Di*Hi/Gi Di-第i楼层的弹性等效刚度,可取该层剪力与层间位移的比值Hi-第i楼层层高
Gi-第i楼层重力荷载设计值 刚重比与结构的侧移刚度成正比关系;周期比的调整将导致结构侧移刚度的变化,从而影响到刚重比。因此调整周期比时应注意,当某主轴方向的刚重比小于或接近规范限值时,应采用加强刚度的方法;当某主轴方向刚重比大于规范限值较多时,可采用削弱刚度的方法。同样,对刚重比的调整也可能影响周期比。特别是当结构的周期比接近规范限值时,应采用加强结构外围刚度的方法规范上限主要用于确定重力荷载在水平作用位移效应引起的二阶效应是否可以忽略不计。见高规5.4.1和5.4.2及相应的条文说明。刚重比不满足规范上限要求,说明重力二阶效应的影响较大,应该予以考虑。规范下限主要是控制重力荷载在水平作用位移效应引起的二阶效应不致过大,避免结构的失稳倒塌。见高规5.4.4及相应的条文说明。刚重比不满足规范下限要求,说明结构的刚度相对于重力荷载过小。但刚重比过分大,则说明结构的经济技术指标较差,宜适当减少墙、柱等竖向构件的截面面积。 长细比 长细比=计算长度/回转半径。 所以很显然,减小计算长度或者加大回转半径即可。 这里需要注意的是,计算长度并非实际长度,而是实际长度乘以长度系数,长度系数则与柱子两端的约束刚度有关。说白了就是
结构动力学历年试题(简答题) 1.根据荷载随时间的变化规律,动力荷载可以划分为哪几类?每一类荷载包括哪几种,请 简述每一种荷载的特点。P2 2.通过与静力问题的对比,试说明结构动力计算的特点。P3 3.动力自由度数目计算类 4.什么叫有势力?它有何种性质。P14 5.广义力是标量还是矢量?它与广义坐标的乘积是哪个物理量的量纲?P16 6.什么是振型的正交性?它的成立条件是什么?P105 7.在研究结构的动力反应时,重力的影响如何考虑?这样处理的前提条件是什么?P32 8.对于一种逐步积分计算方法,其优劣性应从哪些方面加以判断?P132 9.在对结构动力反应进行计算的思路上,数值积分方法与精确积分方法的差异主要表现在 哪里?第五章课件 10.利用Rayleigh法求解得到的振型体系的基本振型和频率及高阶振型和频率与各自的精确 解相比有何特点?造成这种现象的原因何在?P209 11.根据荷载是否预先确定,动荷载可以分为哪两类?它们各自具有怎样的特点?P1 12.坐标耦联的产生与什么有关,与什么无关?P96 13.动力反应的数值分析方法是一种近似的计算分析方法,这种近似性表现在哪些方面? P132及其课件 14.请给出度哈姆积分的物理意义?P81 15.结构地震反应分析的反应谱方法的基本原理是什么?P84总结 16.某人用逐步积分计算方法计算的结构位移,得到如下的位移时程的计算结果:。。。 17.按照是否需要联立求解耦联方程组,逐步积分法可以分为哪两类?这两类的优劣性应该 如何进行判断?P132 18.根据荷载随时间的变化规律,动力荷载可以划分为哪几类?每一类荷载又包括哪些类型, 每种类型请给出一种实例。P2 19.请分别给出自振频率与振型的物理意义?P103 20.振型叠加法的基本思想是什么?该方法的理论基础是什么?P111参考25题 21.在振型叠加法的求解过程中,只需要取有限项的低阶振型进行分析,即高阶振型的影响 可以不考虑,这样处理的物理基础是什么?P115 22.我们需要用数值积分方法求解一座大型的高坝结构的地震反应时程,动力自由度的总数 为25000个,我们如何缩短计算所耗费的机时?P103 23.什么是结构的动力自由度?动力自由度与静力自由度的区别何在?P11及卷子上答案 24.一台转动机械从启动到工作转速正好要经过系统的固有频率(又称为转子的临界转速), 为减小共振,便于转子顺利通过临界转速,通常采用什么措施比较直接有效?简要说明理由。详解见卷子上答案 25.简述用振型叠加法求解多自由度体系动力响应的基本原理及使用条件分别是什么?若 振型叠加法不适用,可采用何种普遍适用的方法计算体系响应?详解见卷子上答案 26.振型函数边界条件。。。 27.集中质量和一致质量有限元的差异和优缺点,采用这两种有限元模型给出的自振频率与 实际结构自振频率相比有何种关系?P242及卷子上答案 28.人站在桥上可以感觉到桥面的震动,简述当车辆行驶在桥上和驶离桥面的主要振型特征 有何不同? 29.简述用Duhamel积分法求体系动力响应的基本原理,以及积分表达式中的t和τ有何差
结构动力学分析 1静力分析与动力学分析的区别 静力分析是分析结构在承受稳定载荷作用下的受力特性。结构动力分析是分析结构在承受随时间变化的载荷作用下的动力学特性。 2动力学特性 动力学特性通常有下面几种类型: 2.1振动特性 即结构的振动形式和振动频率。 2.2随时间变化载荷的效应 例如,对结构位移和应力的效应。 2.3周期(振动)或随机载荷的效应 3四种动力学分析及举例 3.1模态分析 用于确定结构的振动特性,即固有频率和振型。在承受动态载荷的结构设计中,固有频率和振型是重要的参数。模态分析也是其他动力学分析前期必须完成的环节。 举例:如何避免汽车尾气排气管装配体的固有频率与发动机的固有频率相同? 3.2瞬态分析 用于确定结构在受到冲击载荷时的受力特性。 举例:怎样确保桥墩在受到撞击时的安全? 3.3谐响应分析 用于确定结构对稳态简谐载荷的响应。 举例:如何确定压缩机、电动机、泵、涡轮机械等旋转引起的轴承、支座、固定装置、部件应力? 3.4谱分析 用于确定结构在受到动载荷或随机载荷时的受力特性。 举例:如何确定房屋和桥梁承受地震载荷时的受力? 4四种动力学分析基本原理 4.1模态分析理论的基本假设 线性假设:结构的动态特性是线性的,即任何输入组合所引起的输出等于各自输出的组 合,其动力学特性可用一组线性二阶微分方程来描述。任何非线性特性,如塑性、接触单元
等,即使定义了也将被忽略。 时不变性假设:结构的动态特性不随时间而变化,微分方程的系数是与时间无关的常数。 可观测性假设:系统动态特性所需要的全部数据都是可测量的。 遵循Maxwell互易性定理:在结构的i点输入所引起的j响应,等于在j点的相同 输入所引起的i点响应。此假设使结构的质量矩阵、刚度矩阵、阻尼矩阵和频响矩阵都成了对称矩阵。 4.2谐响应分析基本原理 谐响应分析是一种线性分析,非线性特性被忽略。 输入:已知大小和频率的谐波载荷(力、压力和强迫位移);同一频率的多种载荷,可以是相同或不相同的。 输出:位移、应力、应变等。 已知动力学运动方程: [M]{u}+[C]{u}+[K]{u}={F(t)} 其中,[M] 为质量矩阵,[C]为阻尼矩阵,[K]为刚度矩阵,{u}为节点位移向量,{F(t)}载荷为时间的任意函数。对简谐运动而言,{u}和{F(t)}均为简谐形式。 4.3瞬态分析基本原理 瞬态分析也叫时间历程分析。载荷和时间的相关性使得惯性力和阻尼作用比较重要,如果惯性力和阻尼作用不重要,就可以用静力学分析代替瞬态分析。 输入:结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下随时间变化的载荷。 输出:随时间变化的位移、应力、应变等。 瞬态动力学的基本运动方程: [M]{u}+[C]{u}+[K]{u}={F(t)} 其中,[M] 为质量矩阵,[C]为阻尼矩阵,[K]为刚度矩阵,{u}为节点位移向量,{F(t)}载荷为时间的任意函数。 4.4谱分析基本原理 谱分析模态分析的扩展,是将模态分析的结果与一个已知的谱联系起来计算结构的位移和应力。 主要用于分析承受地震或其他随机载荷的建筑物及桥梁结构等。
5章 动力反应的数值计算 如果激励[作用力)(t p 或地面加速度)(t u g ]是随时间任意变化的,或者体系是非线性的,那么对单自由度体系的运动方程进行解析求解通常是不可能的。这类问题可以通过数值时间步进法对微分方程进行积分来处理。在应用力学广阔的学科领域中,有关各种类型微分方程数值求解方法的文献(包括几部著作中的主要章节)浩如烟海,这些文献包括这些方法的数学进展以及它们的精度、收敛性、稳定性和计算机实现等问题。 然而,本章仅对在单自由度体系动力反应分析中特别有用的很少几种方法进行简要介绍,这些介绍仅提供这些方法的基本概念和计算算法。尽管这些对许多实际问题和应用研究已经足够了,但是读者应该明白,有关这个主题存在大量的知识。 5.1 时间步进法 对于一个非弹性体系,欲采用数值求解的运动方程为 )(),(t p u u f u c u m s =++ 或者 )(t u m g - (5.1.1) 初始条件 )0(0u u = )0(0u u = 假定体系具有线性粘滞阻尼,不过,也可以考虑其他形式的阻尼(包括非线性阻尼),后面会明显看到这一点。然而由于缺乏阻尼信息.因此很少这样做,特别是在大振幅运动时。作用力)(t p 由一系列离散值给出: )(i i t p p = ,0=i 到N 。时间间隔 i i i t t t -=?+1 (5.1.2)
图5.1.1 时间步进法的记号 通常取为常数,尽管这不是必需的。在离散时刻i t (表示为i 时刻)确定反 应,单自由度体系的位移、速度和加速度分别为i u 、i u 和i u 。假定这些值是已知的,它们在i 时刻满足方程 i i s i i p f u c u m =++)( (5.1.3) 式中,i s f )(是i 时刻的抗力,对于线弹性体系,i i s ku f =)(,但是如果体系是非弹性的,那么它会依赖于i 时刻以前的位移时程和速度。将要介绍的数值方 法将使我们能够确定i +1时刻满足方程(5.1.1)的反应1+i u 、1+i u 和1+i u ,即在i +1时刻 1111)(++++=++i i s i i p f u c u m (5.1.4) 对于i =0,1,2,3,…,连续使用时间步进法,即可给出i =0,l ,2,3,… 所有瞬时所需的反应。已知的初始条件)0(0u u =)0(0u u =和提供了起动该方法的必要信息。 从i 时刻到i +1时刻的步进一般不是精确的方法,许多在数值上可以实现的近似方法是可能的。对于数值方法,有三个重要的要求:(1)收敛性一随着时间步长的减少,数值解应逼近精确解;(2)稳定性一在存在数值舍入误差的情况下,数值解应是稳定的;(3)精度一数值方法应提供与精确解足够接近的结果。这些重要的问题在本书中均作简要的讨论,全面的论述可在着重微分方程数值解法的书中找到。 本章介绍三种类型的时间步进法:(1)基于激励函数插值的方法;(2)基于速度和加速度有限差分表达的方法;(3)基于假设加速度变化的方法。前两类中各
振型分解反应谱法 振型分解反应谱法是用来计算多自由度体系地震作用的一种方法。该法是利用单自由度体系的加速度设计反应谱和振型分解的原理,求解各阶振型对应的等效地震作用,然后按照一定的组合原则对各阶振型的地震作用效应进行组合,从而得到多自由度体系的地震作用效应。振型分解反应谱法一般可考虑为计算两种类型的地震作用:不考虑扭转影响的水平地震作用和考虑平扭藕联效应的地震作用。 适用条件 (1)高度不超过40米,以剪切变形为主且质量和刚度沿高度分布比较均匀的结构,以及近似于单质点体系的结构,可采用底部剪力法计算。(此为底部剪力法的适用范围) (2)除上述结构以外的建筑结构,宜采用“振型分解反应谱法”。 (3)特别不规则的建筑、甲类建筑和规范规定的高层建筑,应采用时程分析法进行补充计算。 刚重比 刚重比是指结构的侧向刚度和重力荷载设计值之比,是影响重力二阶效应的主要参数 刚重比=Di*Hi/Gi Di-第i楼层的弹性等效刚度,可取该层剪力与层间位移的比值 Hi-第i楼层层高 Gi-第i楼层重力荷载设计值 刚重比与结构的侧移刚度成正比关系;周期比的调整将导致结构
侧移刚度的变化,从而影响到刚重比。因此调整周期比时应注意,当某主轴方向的刚重比小于或接近规范限值时,应采用加强刚度的方法;当某主轴方向刚重比大于规范限值较多时,可采用削弱刚度的方法。同样,对刚重比的调整也可能影响周期比。特别是当结构的周期比接近规范限值时,应采用加强结构外围刚度的方法 规范上限主要用于确定重力荷载在水平作用位移效应引起的二 阶效应是否可以忽略不计。见高规5.4.1和5.4.2及相应的条文说明。刚重比不满足规范上限要求,说明重力二阶效应的影响较大,应该予以考虑。规范下限主要是控制重力荷载在水平作用位移效应引起的二阶效应不致过大,避免结构的失稳倒塌。见高规5.4.4及相应的条文说明。刚重比不满足规范下限要求,说明结构的刚度相对于重力荷载过小。但刚重比过分大,则说明结构的经济技术指标较差,宜适当减少墙、柱等竖向构件的截面面积。 长细比 长细比=计算长度/回转半径。 所以很显然,减小计算长度或者加大回转半径即可。 这里需要注意的是,计算长度并非实际长度,而是实际长度乘以长度系数,长度系数则与柱子两端的约束刚度有关。说白了就是要看与柱相连的梁或者基础是否给力,如果这些构件的刚度越高,那么长度系数就越小,柱子的计算长度也就越短。 具体公式你可以去看钢结构规范,我记得长度系数的具体算法是附录
概念题 1.1 结构动力计算与静力计算的主要区别是什么? 答:主要区别表现在:(1) 在动力分析中要计入惯性力,静力分析中无惯性力;(2) 在动力分析中,结构的内力、位移等是时间的函数,静力分析中则是不随时间变化的量;(3) 动力分析方法常与荷载类型有关,而静力分析方法一般与荷载类型无关。 1.2 什么是动力自由度,确定体系动力自由度的目的是什么? 答:确定体系在振动过程中任一时刻体系全部质量位置或变形形态所需要的独立参数的个数,称为体系的动力自由度(质点处的基本位移未知量)。 确定动力自由度的目的是:(1) 根据自由度的数目确定所需建立的方程个数(运动方程数=自由度数),自由度不同所用的分析方法也不同;(2) 因为结构的动力响应(动力内力和动位移)与结构的动力特性有密切关系,而动力特性又与质量的可能位置有关。 1.3 结构动力自由度与体系几何分析中的自由度有何区别? 答:二者的区别是:几何组成分析中的自由度是确定刚体系位置所需独立参数的数目,分析的目的是要确定体系能否发生刚体运动。结构动力分析自由度是确定结构上各质量位置所需的独立参数数目,分析的目的是要确定结构振动形状。 1.4 结构的动力特性一般指什么? 答:结构的动力特性是指:频率(周期)、振型和阻尼。动力特性是结构固有的,这是因为它们是由体系的基本参数(质量、刚度)所确定的、表征结构动力响应特性的量。动力特性不同,在振动中的响应特点亦不同。 1.5 什么是阻尼、阻尼力,产生阻尼的原因一般有哪些?什么是等效粘滞阻尼? 答:振动过程的能量耗散称为阻尼。 产生阻尼的原因主要有:材料的内摩擦、构件间接触面的摩擦、介质的阻力等等。当然,也包括结构中安装的各种阻尼器、耗能器。阻尼力是根据所假设的阻尼理论作用于质量上用于代替能量耗散的一种假想力。粘滞阻尼理论假定阻尼力与质量的速度成比例。 粘滞阻尼理论的优点是便于求解,但其缺点是与往往实际不符,为扬长避短,按能量等效原则将实际的阻尼耗能换算成粘滞阻尼理论的相关参数,这种阻尼假设称为等效粘滞阻尼。 1.6 采用集中质量法、广义位移法(坐标法)和有限元法都可使无限自由度体系简化为有限自由度体系,它们采用的手法有何不同? 答:集中质量法:将结构的分布质量按一定规则集中到结构的某个或某些位置上,认为其他地方没有质量。质量集中后,结构杆件仍具有可变形性质,称为“无重杆”。 广义坐标法:在数学中常采用级数展开法求解微分方程,在结构动力分析中,也可采用相同的方法求解,这就是广义坐标法的理论依据。所假设的形状曲线数目代表在这个理想化形式中所考虑的自由度个数。考虑了质点间均匀分布质量的影响(形状函数),一般来说,对于一个给定自由度数目的动力分析,用理想化的形状函数法比用集中质量法更为精确。 有限元法:有限元法可以看成是广义坐标法的一种特殊的应用。一般的广义坐标中,广义坐标是形函数的幅值,有时没有明确的物理意义,并且在广义坐标中,形状函数是针对整个结构定义的。而有限元法则采用具有明确物理意义的参数作为广义坐标,且形函数是定义
研究生课程考核试卷 (适用于课程论文、提交报告) 科目:结构动力学大作业教师: 姓名:学号: 专业:岩土工程类别:专硕 上课时间:2015年9 月至2015 年11 月 考生成绩: 卷面成绩平时成绩课程综合成绩阅卷评语: 阅卷教师(签名)
重庆大学研究生院制 土木工程学院2015级硕士研究生考试试题 1 题目及要求 1、按规范要求设计一个3跨3层钢筋混凝土平面框架结构(部分要求如附件名单所示;未作规定部分自定)。根据所设计的结构参数,求该结构的一致质量矩阵、一致刚度矩阵; 2、至少采用两种方法求该框架结构的频率和振型; 3、输入地震波(地震波要求如附件名单所示),采用时程分析法,利用有限元软件或自编程序求出该框架结构各层的线性位移时程反应。
2 框架设计 2.1 初选截面尺寸 取所设计框架为3层3跨,跨度均为4.5m ,层高均为3.9m 。由于基础顶面离室内地面为1m ,故框架平面图中底层层高取 4.9m 。梁、柱混凝土均采用C30, 214.3/c f N mm =,423.010/E N mm =?,容重为325/kN m 。 估计梁、柱截面尺寸如下: (1)梁: 梁高b h 一般取跨度的 112 1 8 ,取梁高b h =500mm ; 取梁宽300b b mm =; 所以梁的截面尺寸为:300500mm mm ? (2)柱: 框架柱的截面尺寸根据柱的轴压比限值,按下列公式计算: ①柱组合的轴压力设计值...E N F g n β= 其中:β:考虑地震作用组合后柱轴压力增大系数; F :按简支状态计算柱的负荷面积; E g :折算在单位建筑面积上的重力荷载代表值,可近似取为 21214/KN m ; n :验算截面以上的楼层层数。 ②c N c N A u f ≥ 其中:N u :框架柱轴压比限值;8度(0.2g ),查抗震规范轴压比限值0.75N u =; c f :混凝土轴心抗压强度设计值,混凝土采用30C ,2 14.3/c f N mm =。
采用振型分解反应谱法时,对于不进行扭转耦联计算的结构,结构j 振型i 质点的水平地震作用标准值,按下列公式计算: G X i ji j j ji F γα= (i=1,2,…n , j=1,2,…n) (2-5) 式中: α j ——相应于j 振型自振周期 T j 的地震影响系数,按图2-1确定; X ji ——j 振型i 质点的水平相对位移; G i ——集中于质点i 的重力荷载代表值; γ j ——j 振型的参与系数, ∑∑=== n 1 i i 2 ji n 1i i ji j G X G X γ 对于进行扭转耦联计算的结构,各楼层可取两个正交的水平位移和一个转角共三个自由度。结构j 振型i 层的水平地震作用标准值,按下列公式计算: G X F i ji ij j xji γα= (2-6a ) G Y F i ji ij j yji γα= (2-6b ) tji j i tj ji 2i G F γ φγα= (2-6c ) 式中: F xji , F yji , tji F ——分别为j 振型i 层的x 方向、y 方向和转角方向的地震作用标 准值; X ji , Y ji ——分别为j 振型i 层质心在x 、y 方向的水平位移; ? ji ——j 振型i 层的相对扭转角; r i ——i 层转动半径,可取i 层绕质心的转动惯量除以该层质量的商的正二次方根: γ ij ——计入扭转的j 振型的参与系数,可按下式确定: 当仅取x 方向地震作用时: ∑∑==++= n 1 i i 2i 2 ji 2 ji 2 ji n 1 i i ji ij r *G Y X G X )(?γ (2-7)
《结构力学》作业参考答案 一、判断题(将判断结果填入括弧内,以 √表示正确 ,以 × 表示错误。) 1.图示桁架结构中有3个杆件轴力为0 。(×) 2.图示悬臂梁截面A 的弯矩值是ql 2。 (×) l l 3.静定多跨梁中基本部分、附属部分的划分与所承受的荷载无关。(√ ) 4.一般来说静定多跨梁的计算是先计算基本部分后计算附属部分。(× ) 5.用平衡条件能求出全部内力的结构是静定结构。( √ ) 6.求桁架内力时截面法所截取的隔离体包含两个或两个以上的结点。(√ ) 7.超静定结构的力法基本结构不是唯一的。(√) 8.在桁架结构中,杆件内力不是只有轴力。(×) 9.超静定结构由于支座位移可以产生内力。 (√ ) 10.超静定结构的内力与材料的性质无关。(× ) 11.力法典型方程的等号右端项不一定为0。 (√ ) 12.计算超静定结构的位移时,虚设力状态可以在力法的基本结构上设。(√) 13.用力矩分配法计算结构时,汇交于每一结点各杆端分配系数总和为1,则表明分配系 数的计算无错误。 (× ) 14.力矩分配法适用于所有超静定结构的计算。(×) 15.当AB 杆件刚度系数i S AB 3 时,杆件的B 端为定向支座。 (×)
二、单项选择题(在每小题的四个备选答案中选出一个正确答案,并将其代号填在题干后面的括号内。不选、错选或多选者,该题无分。) 1.图示简支梁中间截面的弯矩为( A ) q l A . 82ql B . 42ql C . 22 ql D . 2ql 2.超静定结构在荷载作用下产生的内力与刚度(B ) A . 无关 B . 相对值有关 C . 绝对值有关 D . 相对值绝对值都有关 3.超静定结构的超静定次数等于结构中(B ) A .约束的数目 B .多余约束的数目 C .结点数 D .杆件数 4.力法典型方程是根据以下哪个条件得到的(C )。 A .结构的平衡条件 B .结构的物理条件 C .多余约束处的位移协调条件 D .同时满足A 、B 两个条件 5. 图示对称结构作用反对称荷载,杆件EI 为常量,利用对称性简化后的一半结构为(A )。 6.超静定结构产生内力的原因有(D ) A .荷载作用与温度变化 B .支座位移 C .制造误差 D .以上四种原因
第五章 结构动力学中的常用数值方法 5.1.结构动力响应的数值算法 ... . 0()(0)(0)M x c x kx F t x a x v ? ++=??=??=?? 当c 为比例阻尼、线性问题→模态叠加最常用。但当C 无法解耦,有非线性存在,有 冲击作用(激起高阶模态,此时模态叠加法中的高阶模态不可以忽略)。此时就要借助数值积分方法,在结构动力学问题中,有一类方法称为直接积分方法最为常用。所识直接是为模态叠加法相对照来说,模态叠加法在求解之前,需要对原方程进行解耦处理,而本节的方法不用作解耦的处理,直接求解。(由以力学,工程中的力学问题为主要研究对象的学者发展出来的) 中心差分法的解题步骤 1. 初始值计算 (1) 形成刚度矩阵K ,质量矩阵M 和阻尼矩阵C 。 (2) 定初始值0x ,. 0x ,.. 0x 。 (3) 选择时间步长t ?,使它满足cr t t ?,并计算 02 1() a t = ?,112a t = ?,202a a = (4) 计算. .. 001 1122t x x x x a a -?=- + (5) 形成等效质量阵01M a M a C - =+ (6) 对M -阵进行三角分解T M LDL - = 2.对每一时间步长 (1) 计算时刻t 的等效载荷 2 1 ()( )t t t t t Q Q K a M x a M a C x - -?=-- -- (2) 求解t t +?时刻的位移 ()T t t t L D L x Q - +?= (3) 如需要计算时刻t 的速度和加速度值,则 . 1()t t t t t x a x x +?-?=- .. 0(2)t t t t t t x a x x x +?-?=-+
MIDAS/GEN六层框架结构的动力分析 工程概况 建筑地点:北京市 建筑类型:六层综合办公楼,框架填充墙结构。 地质条件:根据设计任务说明地震设防烈度为8度。 柱网与层高:本办公楼采用柱距为6.0m的内廊式小柱网,边跨为6.0m,中间跨为2.7m,层高取首层为4.5m,其余为3.3m,如下图所示: 框架结构的计算简图:
典型结构单元 梁、柱、板截面尺寸的初步确定: 1、梁截面高度一般取梁跨度的1/12至1/8。本方案取1/10×6000=600mm,截面宽度取600×1/2=250mm,可得梁的截面初步定为b×h=250*600。楼板取120mm,楼梯板及休息平台板为100mm,平台梁250×400。 2、框架柱的截面尺寸 梁截面尺寸(mm) 柱截面尺寸(mm)
结构动力学分析用来求解随时间变化的载荷对结构或部件的影响。与静力分析不同,动力分析要考虑随时间变化的力载荷以及它对阻尼和惯性的影响。MIDAS/GEN可进行的结构动力学分析类型包括:瞬态动力学分析、模态分析、屈曲分析、动力非线性分析等。本文将以一个六层框架结构为例对结构进行模态分析和谱分析。 一.模态分析 模态分析是用于确定设计中的结构或机器部件的振动特性。它也是其他更详细动力学分析的起点,例如瞬态动力学分析和谱分析等,可以通过模态分析确定结构部件的频率响应和模态。一般对于动力加载条件下的结构设计而言,频率响应和模态是非常重要的参数,即使在谱分析及瞬态分析中也是需要的。 1.1动力学求解方法 MIDAS目前提供了三种特征值分析方法,它们是子空间法、分块Lanczos 算法、多重Ritz向量法。本文采用子空间法进行计算求解。子空间法使用迭代技术,求出结构的前r阶振型,它内部使用广义Jacobi迭代算法。由于该方法采用了完整的质量和刚度矩阵,因此精度很高,但计算速度较慢,特别适用于大型对称特征值求解问题。分块Lanczos法特征值求解器采用Lanczos算法,Lanczos算法是用一组向量来实现递归计算。这种方法和子空间法一样精确,但速度较快。多重Ritz向量法以变分原理为基础,直接迭加法求出的是和激发荷载向量直接相关的振型,其收敛具有严格的理论基础,在物理和、力学的微分方程中占有很重要的位置,得到广泛的应用。 1.2本工程模态分析结果 1.2.1自振周期与振型: 使用MIDAS/GEN中的模态分析计算结构的自振周期和振型。模态分析所使用的方法是子空间迭代法。高层建筑结构振型多,分布规律很难掌握,扭转振动会对结构产生教大影响,因此不能简单的取前几阶进行计算。规范中规定对于高层结构一般取3}5阶振型。为使高层建筑的分析精度有所改进,其组合的 振型个数适当增加。考虑到MIDAS/GEN软件的强大快速的数据处理能力和精度的要求,本文取30阶振型。从国内高层建筑结构设计经验来看,建议基本自振周期按以下的几个公式估计,其中N为地面以上建筑物结构层数。经验公式表达简单,使用方便,但比较粗糙,而且只有基本周期,但经常用于对理论计算值的计算与评价。 框架:T1=(0.08~0.1) N 框架一剪力墙:T1=(0.06~0.09) N 钢结构:T1=0.1N 本工程得经过MIDAS/GEN的计算得到固有周期、固有频率、振型参与质 量等的数值结果;X方向振型参与达到总质量的95.57%, Y方向振型参与达到 总质量的94.83%,经过整理取前十阶列表可得到表1
华中科技大学土木工程与力学学院 《结构动力学》考试卷(B卷、闭卷) 2013~2014学年度第一学期成绩 学号专业班级姓名 一、简答题(每题5分、共25分) 1、刚度法和柔度法所建立的体系运动方程间有何联系?各在什么情况下使用方便? 答:从位移协调的角度建立振动方程的方法为柔度法。从力系平衡的角度建立的振动方程的方法为刚度法。这两种方法在本质上是一致的,有着相同的前提条件。在便于求出刚度系数的体系中用刚度法方便。同理,在便于求出柔度系数的体系中用柔度法方便。在超静定结构中,一般用刚度法方便,静定结构中用柔度法方便。 2、什么叫动力系数,动力系数大小与哪些因素有关?单自由度体系位移动力系数与内力动力系数是否一样? 答:动力系数是指最大动位移[y(t)]max与最大静位移yst的比值,其与体系的自振频率和荷载频率θ有关。当单自由度体系中的荷载作用在质量处才有位移动力系数与内力动力系数一样的结果。 3、什么叫临界阻尼?怎样量测体系振动过程中的阻尼比?若要避开共振应采取何种措施? 答:当阻尼增大到体系在自由反应中不再引起振动,这时的阻尼称为临界阻尼。根据公式即测出第k次振幅和第k+n次振幅即可测出阻尼比。 措施:○1可改变自振频率,如改变质量、刚度等。○2改变荷载的频率。○3可改变阻尼的大小,使之避开共振。 4、振型正交的物理意义是什么?振型正交有何应用?频率相等的两个主振型互相正交吗? 答:物理意义:第k主振型的惯性力与第i主振型的位移做的功和第i主振型的惯性力与第k主振型的静位移做的功相等,即功的互等定理。 作用:○1判断主振型的形状特点。○2利用正交关系来确定位移展开公式中的系数。 5、应用能量法求频率时,所设的位移函数应满足什么条件?其计算的第一频率与精确解相比是偏高还是偏低?什么情况下用能量法可得到精确解? 答:所设位移函数要满足位移边界条件,同时要尽可能与真实情况相符。第一频率与精确解相比偏高。如果所假设的位移形状系数与主振型的刚好一致,则可以得到精确解。