地震反应分析.docx
- 文档编号:12886149
- 上传时间:2023-04-22
- 格式:DOCX
- 页数:75
- 大小:1.25MB
地震反应分析.docx
《地震反应分析.docx》由会员分享,可在线阅读,更多相关《地震反应分析.docx(75页珍藏版)》请在冰豆网上搜索。
地震反应分析
结构地震反应分析
结构地震反应分析的主要工作是首先将结构简化成力学分析模型,然后输入地震作用,计算模拟结构的反应行为,包括内力和变形反应时程或最大值。
其目的是为结构抗震设计提供必要的数据资料;或为抗震安全鉴定和拟定抗震加固方案提供参考依据;或为研究结构破坏机理提供基本手段,从而改善设计,提高结构的抗震性能。
结构地震反应取决于地震动输入特性和结构特性。
随着人们对地震动特性和结构特性的了解越来越多,特别是技术手段越来越先进,结构地震反应分析方法也跟着有了飞跃的发展。
结构抗震分析方法的发展大体上可分为三个阶段,即静力法、拟静力法(通常指反应谱方法)和动力法阶段。
静力法是20世纪初首先在日本发展起来的。
该方法将结构物看成是刚体,并刚接于地面。
这样,结构在最大水平加速度绝对值为amax的地面运动激励下,受到的最大水平作用力P(即最大惯性力)为
其中,一是结构物的重量,k是地面最大水平加速度绝对值Amax与重力加速度g之比,称为地震系数。
在当时人们对地面运动的频谱和卓越周期的了解还不够多,以及房屋多为
低层建筑的情况下,应用上述地震荷载计算公式于抗震
设计还是可以的。
但是,随着地震资料的积累和城市与工业建设的发展,使
人们认识到作为静力法基础的刚性结构假定已明显地远离实际情况,于是考虑结构物的弹性性质、阻尼性质及相应动力特性的反应谱方法便发展起来
反应谱方法出现在20世纪40年代。
美国的一些学者在取得了一部分
强震地面运动记录之后,考虑地震动特性与结构动力特性共同对结构地震反应产生决定性影响的这一事实,提出了反应谱概念和相应的设计计算方法。
这一方法有动力法的内容,却具静力法的形式,故可称之为拟静力法。
该方法对结构地震反应分析产生巨大影响,至今仍是结构抗震设计的主要计算方法。
尽管反应谱方法取得的进步是实质性的,但它的应用还是受到一些限
制,如原则上只能用于线性结构体系;不能真实反映复杂结构体系的动力放大作用。
因此,随着重大工程的不断兴建和计算机技术的飞速发展,20世纪70年代,结构地震时程反应分析得到全面发展。
相对于反应谱方法而言,时程反应分析是一种动力分析方法,它求取的
不是结构的某种最大反应或其近似估计,而是结构在地震激励下的反应时间历程,即地震与结构相互作用的过程,其结果更为可靠。
另外,时程反应分析可以真正处理非线性问题,这是结构地震反应分析一个非常重要的方面。
随计算机和有限元技术的发展,结构分析模型也经历了一个由极其简化到相对较少简化的过程。
以前大家熟悉的一些简化分析模型,如剪切模型,考虑梁变形作用的D值法以及框架剪力墙协同工作体系模型等,在当前的研究与设计中已很少使用,取而代之的是三维空间有限元分析模型。
目前,各种大型有限元程序为结构地震反应分析提供了强有力的工具。
应用这些程序,结构弹性地震反应分析已不存在问题,无论多么复杂的结构体系,只要计算模型简化的合理都能得到满足一定工程精度要求的结果。
结构弹塑性地震反应性态极其复杂,尽管经科硏人员数十年的努力,发展了一些分析方法,但仅较规则的结构二维弹塑性分析可以取得基本令人满意结果,量大面广的复杂结构的分析方法至今未能很好解决,它是今后有尖科研人员需要重点解决的课题。
弹塑性地震反应可以分为静力弹塑性反应分析和动力时程反应分析。
静力弹塑性地震反应分析一般指近年为满足性态抗震设计而发展的pushover分析方法,该方法的主要步骤是首先将地震荷载等效成某种分布形式的静力荷载,用静力弹塑性分析方法求得结构的基底剪力与位移尖系曲线,即结构能力曲线,然后将结构等效成单自由度体系并将结构的能力曲线和地震输入谱曲线转换成相同坐标格式,根据两曲线的交点确定结构位移反应。
研究表明这一分析方法在分析中低层剪切型结构时,可以提供较满意的弹塑性位移反应估计结果,而分析高层结构时,则误差较大,基本不适用高层结构
地震反应分析
估计结构的弹塑性地震反应行为,较准确、可靠的方法无疑是弹塑性时程反应分析方法。
但目前仅二维分析方法发展的较为成熟,并在硏究中得到广泛应用,问题是大部分实际需要进行弹塑性分析的结构形式均较为复杂,难以简化成合理的二维分析模型,如勉强进行二维分析模型简化,必将导致分析结果的较大误差。
结构三维弹塑性地震反应分析一直是结构抗震分析中有待解决的难题。
主要困难是如果以构件作为单元,则构件三维受力状态下的恢复力模型难以确定;如果采用较精细的非线性有限元模型,则依然存在构件三维受力状态下的本构矢系难以确定问题,且计算量也难以接受。
尽管如此,有矢科研人员还是在努力探索,不断提出一些新的模型和计算方法,使结构三维弹塑性地震反应分析中存在的问题逐步得到解决。
结构动力反应分析
结构动力分析是计算分析结构在动荷载作用下的变形和内力,校核结构是否满足指定安全要求,它是结构工程领域的重要环节。
动荷载是指随时间而改变的荷载。
同样,动荷载作用下结构的反应(内力及变形)也是随时间而改变的。
结构动力反应分析主要内容涉及建立计算分析模型和体系运动方程,确定结构特性参数,选择合理的方法进行运动方程求解等。
分析方法可分为两大类:
确定性分析方法和不确定性即随机反应分析方法。
选取哪种分析方法取决于荷载和结构参数是否可以给定,如果它们完全是已知的,贝慄用确定性分析方法,如它们并不完全已知,但可从统计意义上定义,则一般采用随机反应分析方法。
从荷载大小和结构是否进入弹塑性状态角度考虑,分析方法可进一步分为弹性反应分析和弹塑性反应分析。
前者假定在动荷载作用下,结构始终处于初始弹性状态。
后者则要考虑随时间荷载的作用,结构参数不断改变情况下的反应。
在地震反应分析中,分析方法还可以分为静力法、反应谱方法和时程分析方法。
早期采用的静力法非常简单,即结构承受的侧向地震力等于结构质量乘以地面运动峰值加速度。
反应谱方法是首先将结构进行振型分解,然后根据指定的地震反应谱确定每个振型的反应,最后将这些叠加求出结构总体反应。
反应谱方法是一种伪动力方法,它没有考虑结构随时间的动力反应过程,但是考虑了结构动力特性对其反应的影响。
时程分析方法则完全是求解结构的振动反应过程。
运动方程
结构动力反应分析的目的是计算结构在给定随时间变化荷载作用下的变形和内力的反应过程。
我们知道,一个结构具有无限多个变形自由度,但在大多数情况下,应用包含有限个自由度的近似分析方法,就足够精确了。
这样,问题就变为求出这些所选定的有限个变形分量的时间过程。
描述动力变形的数学表达式称为结构的运动方程,而这些运动方程的解就提供了所求的变形过程。
理方法和哈密尔顿原理方法。
动力体系的运动方程可用上述三种不同方法中的任一种来建立.最简单明了的方法是采用直接平衡法建
立作用于体系上全部力(包括惯性力)的动力平衡方程,但对于更复杂的体系,特别是对那些质量和弹性只在有限区域是分布的体系,直接的矢量平衡可能是困难的,而应用仅包含功和能等标量来建立方程式的方法则更为方便,其中最直接的就是基于虚位移原理的方
法。
在这种方法中,首先计算作用于体系上的力,然后由它们在相应的虚位移上所作的功来导出运动方程。
哈密尔顿原理是利用能量
来建立运动方程,它不直接利用作用于体系内的惯性力或保守力,而是采用体系的动能和位能的变分来替代这些力的作用。
上述三种方法是完全相等的,采用哪种方法取决于是否方便和个人的喜爱以及动力体系的性质。
直接平衡法
任何动力体系的运动方程都可以用牛顿第二运动定律表示:
即
任何质量m的动量变化率等于作用在这个质量上的力。
这个尖系在数
学上可用微分方程来表达:
p(t)a(m空)
dtdt
其中p(t)为作用力矢量,V(t)为质量m的位置矢量•对于大多数的结
构动力学问题,可以假设质量是不随时间变化,这时方程(写作1)可改
d2v
P(t)m—2mv(t)(2
dt2)
它表示力为质量与加速度的乘积,
(2)式也可改写为:
p(t)mv(t)0
此时第二项mv(t)被称为抵抗质量加速度的惯性力
质量所产生的惯性力与它的加速度成正比,但方向相反。
这个概念称作为d'Alembert原理。
由于它可以把运动方程表示为动力平衡方程,因而是结构动力学问题中一个很方便的方法。
可以认为,力P(t)包括许多种作用于质量上的力:
抵抗位移的弹性约束力,抵抗速度的粘滞力,以及独立确定的外荷载。
因此,如果引人抵抗加速度的惯性力,则运动方程的表达式仅仅是作用于质量上所有力的平衡表达式。
在许多简单问题中,最直接而且方便的建立运动方程的方法就是直接平衡法。
虚位移原理
如果结构体系相当复杂,而且包含许多彼此联系的质量点或有限尺寸的质量块,则直接写出作用于体系上所有力的平衡方程可能是困难的。
通常所包含的各式各样的力可以容易地用位移自由度来
表示,而它们的平衡规律则可能是不清楚的。
此时,虚位移原理就可用来代替平衡规律建立运动方程。
虚位移原理可阐述如下:
如果一个平衡的体系在一组力的作用下承受一个虚位移,即体系约束所允许的任何微小位移,则这些力所作的总功将等于零•按这个原理,很明显,虚位移时所作的功为零是和平衡等价的。
因此,在建立动力体系的反应方程时,首先要搞清作用于体系质量上的所有的力,包括按照d'Alembert原理所定义的惯性力,然后引入相应于每个自由度的虚位移,并使所作的功等于零,这样就可以导出运动方程。
这个方法的主要优点是:
虚功为标量,可以按代数方式相加,而作用于结构上的力为矢量,只能按矢量叠加
利用虚位移原理建立运动方程的简单实例见单自由度体系运动方程。
哈密尔顿原理
避免建立平衡矢量方程的另一个方法是使用以变分形式表示的能量(标量)的方法。
通常最广泛应用的变分概念为哈密尔顿
(Handlton)原理。
此原理可表达为
tt
22
其中:
T体系的总动能;
V体系的位能,包括应变能及任何保守外力的势能;
Wnc作用于体系上的非保守力(包括阻尼力及任意外荷载)所作的功;在指定时间区间内所取的变分。
哈密尔顿原理说明在任何时间区间ti到t2内,动能和位能的变分加上所考虑的非保守力所作的功的变分必须等于零。
这个原理的应用直接导出任何给定体系的运动方程。
这个方法和虚功方法的不同在于:
在这个方法中,不明显使用惯性力和弹性力,而分别被动能和位能的变分项所代替。
因此,这种建立方程的方法的优点是,它只和纯粹的标量一能量有矢,而在虚功分析中,尽管功的本身是标量,但被用来计算功的力和位移却都是矢量。
哈密尔顿原理也可用于静力问题。
此时,动能项T消失,而方程(4)的积分中剩余的项是不随时间变化的,于是方程简化为
(VWnc)0
(2)这就是广泛应用在静力分析中著名的最小位能原理。
利用哈密尔顿原理建立运动方程的简单实例见单自由度体系运动方程。
单自由度体系运动方程
一个理想化单自由度运动体系如图1所示
(b)
IX
*PO)
k
图1理想化单自由度运动体系
(a)基本元件;(b)平衡力系
图中:
k■弹性刚度,m■质量,c・阻尼系数,P⑴■外力,x■位移坐标,fi■惯性力mx(x•加速度),fs■弹性恢复力kx,fD-
阻尼力CX(假设粘滞阻尼,x■速度)。
采用直接平衡法来建立运动方程。
如图
lb所示,根据力的平衡可以得到:
(1)
f|fofsP(t)
公式
(2)即为单自由度体系运动方程。
mxexkxp(t)
(2)
再采用虚位移原理来建立此单自由度体系的运动方程。
假设给
图lb所示体系一个虚位移X(仅仅是体系约束所允许的微小位移)则每个力都
将做功,体系所作的总功可写作
f|XfDxfsXp(t)X0
(3
)
式中的负号是由于力的方向和虚位移方向相反,将各种力的表达式代入方程
(3),并提取公因子x,可以得到
[mxexkxp(t)]x0(4)
因为X不等于零,所以可以得到同式
(2)—样的方程。
继续利用哈密尔顿原理来推导这个体系的运动方程。
根据定义,体系的动能为T1/2mx2,而仅由弹簧的应变能U表达的位能为VU1/2kx2,该体系的非保守力为阻尼力fD和外荷载p(t),这
些力所做功的变分为Wncp(t)XCXX5将以上各式代入哈密尔顿原
理的表达式,并经相应的变分和整理后可得
t2
(5)
(6)
[mxxexxkxxp(t)x]dt0
ti
上式中的第一项可以进行如下的分部积分:
"“t2mxxdt
mxxdtti
timxx
式中利用了尖系xd(x)/dt<>因为在哈密尔顿原理中假定变分x
在积分限ti和乜时为零,所以方程(6)右边第一项为零,将方程(6)代回方程(5),结果为
t2
[mxexkxp(t)]xdt0(7)
ti
因变分X的任意性,所以必须括号内的表达式等于零时才能使方程始终得以满
足。
这样就得到了同样的单自由度体系运动方程。
当受到地震作用时,不仅需要考虑体系相对地面的加速度,还
要考虑地面的加速度,因此,惯性力fim(xXg)0这时只需将mxg
移到方程的右端作为一种特殊的外荷载处理即可。
如果p(t)二。
,式
(2)贝y成为单自由度体系自由振动方程:
2x0
式中:
一八,称为阻尼比;.、k/m,称为圆频率,是表示体
2m
系自振特性的一个参数。
运动方程(8)的解为:
x(t)et_xo)_sinotx(0)cosDt(9)
D
式中:
x(0)和x(0)分别代表体系的初始速度和初始位移;
d.厂''称为阻尼振动圆频率。
式(4)表示在初始扰动下'
体系随时间变化逐渐衰减的振动过程。
体系在外力作用下的运动特性可以通过假设图1所示体系承受
幅值Po、圆频率为的简谐荷载P(t)作用下的运动微分方程求解来加
以说明,此时运动方程为:
x2x2x-p°sint(10)
m
该方程的解由通解和特解两部分组成:
x(t)通elASinDtBcosDt)(11)
示对所作用荷载的瞬态反应,常数A和B可由给定的初始条件算出,由于阻尼
使此项很快消失,因此一般不予考虑。
方程(10)的特解
为作用荷载同频率而不同相位的稳态反应。
多自由度体系运动方程
图1梁式结构的离散化
采用直接平衡法和图1所示的普通简支梁作为典型例子来说明多自由度体系的运动方程是如何建立的,这种方法对任何一种结构类型都同样适用。
假定这个结构的运动由梁上一系列离散点的位移
Ui(t),U2(t),LUn(t)所确定。
结构上的这些点可以任意设置,
所考虑的自由度(位移分量)数目取决于分析者,当然取较多时能更好地逼近结构真实的动力行为。
本例中梁上每一个节点只取一个位移分量。
然而,应该指出每一个节点也可以取几个位移分量,例如,可以取转角和水平向位移等。
针对体系的每一个自由度列出实际作用力的平衡方程就能写出图1所示体系的运动方程。
一般说来在任意一点i上,包含有四种力:
外
荷载Pi(t)、由于运动而产生的力(即惯性力fli)、阻尼力仏和弹性力
fsi。
这样,对于多自由度体系中的每一个自由度,动力平衡条件可写成
当力向量用矩阵形式表示时,亦可写成
fifofsp(t)
(2)
这就是多自由度体系的运动方程。
每一种抗力可以非常方便地用一组适当的影响系数来表示。
例
如,考虑在节点1上产生的弹性力分量,这个量一般依赖于结构所有节点产生
写成一般形式,
应于i坐标的力
用矩阵形式表达,弹性力与位移的矢系为
fs2
k21k22
k2i
u2
(4)
fsi
ki1ki2
kii
Ui
用符号表示:
fsku(5)
刚度系数矩阵k称为结构的刚度矩阵(针对指定的一组位移坐标)u是表示结构变形形状的位移向量。
若假定阻尼与速度有矢,即粘滞阻尼,则与所选择的自由度对应的阻尼力
就可以按同样的方式用阻尼系数表示。
类似于公式(4)和(5),阻尼力与
各点位移坐标的速度的矢系为:
fDCU(6)
这里U为速度向量,向量中的Ui表示i点位移坐标的时间变化率(速度),C为阻尼矩阵,矩阵中的Cij称作阻尼影响系数,定义为:
CiU由j坐标单位速度所引起的对应于i坐标的力
惯性力也可用一组影响系数表示,它们表示相应自由度的加速度与其产生
的惯性力之间的矢系。
类似于式(4)和(5),惯性力可表达成
fimu
(7
)质量矩阵,矩阵中的mj称为质量系数,定义为:
由坐标j的单位加速度所引起的对应于i坐标的力。
其中u为加速度向量,向量中的Ui表示i点位移坐标的加速度,m是
将式(5),(6)和(7)代入式
(2),可给出结构完整的动力平衡方程:
mucukup(t)⑻
这就是用以求解多自由度体系反应的运动方程。
当受到地震作用时,类似于单自由度体系,此时fimumlug(I为单位列
向量),将mlug移到式(8)的右端即得到体系在地震作用下的运动方程。
刚度矩阵
刚度矩阵k由刚度系数心组成,定义为:
kj=由j坐标单位位移所引起的对应于i坐标的力。
原则上,与任何一组指定的节点位移相矢的刚度系数都可直接应用它们的定义求得。
然而,实际应用中基本上都是采用有限元方法来计算刚度矩阵,因它非常简便直观。
其步骤是首先将结构分割成只在有限个接点处相互连接的离散单元体系,然后计算单个单元特性,最后适当地组合形成结构总刚度矩阵。
这样,确定结构刚度特性的问题基本上可简化为单元刚度计算问题。
单元刚度矩阵可以根据假定的单元内部变形插值函数和虚功原理求出。
面以图1所示二维弯曲梁单元为例简单说明如何建立单元刚
度矩阵°
当小发生单位位移而同时其它三个自由度为约束时,梁的挠曲
线可用多项式函数表示为:
!
3
2X
(la)
相应的梁的挠曲线
Xi(x)13・2r
与此类似,当出和e2分别发生单位变形时,
可表示为:
2(x)X(1
7)2
(lb)
2
X3(x)3・
3
I
(1C)
X2X
4(X)
11
1
(Id)
i(X)称为插值函数或形函数。
这样,
单元的挠曲形状就可以用它的
结点位移表示成:
V(x)1(x)lli2(x)13(X)U24(X)2
(2)
梁的应变能是
1HH
kjoEI(x)i(x)jdx
由式(4)可以推出二维等截面弯曲梁的单元刚度矩阵为:
12611261
以上只是介绍了最主要步骤,欲知具体推导过程可参见相矢有限元方法书籍。
建立单元刚度矩阵后,采用有限元规定的组装方法即可建立结构总体刚度矩阵。
质量矩阵
结构动力分析时需要考虑所有因振动而产生惯性力的质量。
这些质量一般以质量矩阵表示。
质量矩阵可分为集中质量矩阵和一致质量矩阵。
集中质量矩阵是假定全部质量积聚在某些需要计算平动位移的结点上。
采用静力学方法将单元质量等效到单元的各个节点上,结构任意节点积聚的总质量等于与该节点连接的各单元分配给此节点的质量之和,当然还要加上用同样方法等效到该节点的外部荷载质量。
这样形成的质量矩阵具有对角形式,即矩阵中只有对角线上有值,对角线以外值均为零。
一致质量矩阵是指根据有限元原理推导出的质量矩阵。
质量矩阵中的元素g定义是:
由坐标j的单位加速度所引起的对应于i坐标的力。
其形成过程同刚度矩阵一样也是先形成单元质量矩阵,然后再组装成总质量矩阵。
单元质量矩阵的建立亦类似于单元刚度矩阵,具体可参见刚度矩阵辞条中弯曲梁单元刚度矩阵的形成过程。
采用类似方法,可得到弯曲梁单元的质量系数
I
rriij°m(x)i(x)j(x)dx(l)
这里,(x)是插值函数。
事实上,式
(1)适用于任何一种单元,只不过不同的单元要选用不同的插值函数和积分区域而已。
集中质量矩阵非常简便,且大量计算结果表明具有较好的计算精度,因此,实际结构分析中较多采用这种质量矩阵。
阻尼矩阵
如果作用在结构上的各种阻尼力能够定量确定的话,那么有限元方法类似于在刚度矩阵辞条中介绍的刚度矩阵确定方法那样可以再一次用来确定体系的阻尼矩阵。
这样可以得到单元的阻尼系数为
cooc(X)i(X))(X)dX
(1)
采用不同的插值函数和相应的积分域即能够得到单元阻尼矩阵。
然而阻尼特性C(x)(或其它任何特殊的阻尼特性)实际上是算不出来的。
因此常常根据类似结构实验方法所确定的阻尼比来表示阻尼,而不是用一个显式的阻尼矩阵C。
当实际动力分析中需要显示阻尼矩阵时,常常采用瑞利阻尼矩阵或比
例阻尼矩阵。
瑞利阻尼矩阵
瑞利阻尼假设阻尼矩阵是结构质量矩阵和刚度矩阵的线性组合,即
(1)
式中和为比例常数。
根据振型正交条件,
与第1阶振型阻尼比之间应满足矢系
(2)
指定结构的两个自振频率
i和j以及相应的阻尼比:
和「可由式
(2)求出
瑞利阻尼在两个指定振型内,阻尼比小于或等于指定阻尼,其他
振型阻尼比随频率增高或减小不断增大,这样就削弱了这些振型对结构反应的
影响,见图1所示。
比例阻尼矩阵
假定ck,称为刚度比例阻尼矩阵,由瑞利阻尼矩阵中方程
(2)可以得到i72,当用第一振型的阻尼比确定后,贝卩其
它振型越高阻尼比越大,因此该模型限制了高阶振型对结构反应的影响。
假定cm,称为质量比例阻尼,由瑞利阻尼矩阵中方程
(2)可以得到i/
(2,),若用第一振型的阻尼比确定后,则与刚度比例阻尼相反,振型越高阻尼比越小,因此该模型突出了高阶振型
对结构反应的贡献。
自振特性
每个结构都有自己的振动特性,该特性以自振频率(或自振周期)和振型表示,它们可以通过振动测量方法或计算分析方法得到。
结构振动具有周期性,振动一个来回所需的时间称为振动周期,以T来表示,单位为秒。
每秒内振动的次数称为振动频率f,单位为赫兹。
另一个参数是圆频率,单位为弧度/秒,通常亦可简称为自
振频率。
三者之间的矢系是:
T=*/f=2/。
单自由度体系的自振频率为
\k/m
(1)
其中,k・弹性刚度,m•质量。
多自由度体系的自振频率和振型可以通过其运动方程mucukup(t)
(见多自由度体系运动方程)求得。
在上式中令P(t)0,另外因阻尼对结构自振特性影响不大,因此可进一步忽略阻尼力,这样可得到运动方程如下:
kumu0
(2
)
设结构作简谐运动uuxost,将其代入式
(2),可得到齐次方程
因自由振动振幅U。
不全为零,所以上式括号内矩阵的行列式之值必须等于零,由此得到结构自振频率方程:
k2m0
(4
)
结构的刚度矩阵k和质量矩阵m都是n阶方阵,其中n是节点自由度的数目,所以式(4)是尖于的n次代数方程,由此可以求出结构的n个自振频率。
对于结构的第i个自振频率,由式(3)可确定一组振幅值uo,u。
中各点幅值应保持固定比例,但绝对值可以任意变化,它们构成一个向量,称为特征向量,工程上通常称为结构的振型。
运动方程求解方法
运动方程求解方法可以分为解析方法和数值方法。
单自由度体系在简单荷载作用下可以采用解析方法
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地震 反应 分析