Midas静力弹塑性分析.docx
《Midas静力弹塑性分析.docx》由会员分享,可在线阅读,更多相关《Midas静力弹塑性分析.docx(22页珍藏版)》请在冰豆网上搜索。
Midas静力弹塑性分析
■简介
Pushover分析是考虑构件的材料非线性特点,分析构件进入弹塑性状态直至到达极限状态时结构响应的方法oPushover分析是最近崔地震研究及耐震设计中经常釆用的基于性能的耐震设计(Performance-BasedSeismicDesign,PBSD)方法中最具代表性的分析方法。
所谓基于性能的耐震设计就是由用户及设计人员设定结构的目标性能(targetperformance),并使结构设计能满足该目标性能的方法。
Push""分析前要经过一般设计方法先进行耐爲设计使结构满足小爲不坏、中震可修的规范要求,然后再通过pushover分析评价结构在大震作用下是否能满足预先设定的目标性能。
计算等效地震靜力荷栽一般釆用如图所示的方法。
该方法定遇过反应修正系数(R)莎设计荷载降低并使结构能承受该荷栽的方法。
在这里使用反应修正系数的原因是为了考虑结构进入弹塑性阶段时吸收地震能量的能力,即考虑结构具有的延性使结构超过弹性极限后还可以承受较大的塑性变形,所以设计时的地馬作用就可以比对应的弹性结构折减很多,设计将会更经济。
目前我国的抗震规范中的反应谱分析方法中的小震彩响系数曲线就是反应了这种设计思想。
这样的设计方法可以说是基于荷载的设计(force-baseddesign)方法。
一般来说结构刚度建大釆用的修正系数R越大,一般在广10之问。
但是这种基于荷裁与抗力的比较进行的设计无法预测结构实际的地震咱应,也无法从各构件的抗力推测出整体结构的耐震能力,设计人员在设计完成后对结构的耐震性能的把握也是模糊的。
基于性能的耐馬设计中可由开发商或设计人员预先设定目标性能,即在预想的地震作用下事先设定结构的破坏程度或者耗能能力,并使结构设计满足该性能目标。
结构的耗能能力与结构的变形能力相关,所以要预测到结构的变形发展情况。
所以基于性能的耐震设计经常通过评价结构的变形来实现,所以也可称为基于位移的设计(disp1acement-baseddesign)<>
图基于稜載的设计方法中地震作用的计算
Pushover分析是评价结构的变形性能的方法之一,分析后会得到如图所示的荷
载一位移能力谱曲线。
另外,根据结构耗能情况会得到弹塑性需求谱曲线。
两个曲线的交点就是针对该地震作用结构所能发挥的最大内力以及最大位移点。
当该交点在目标性能范围内,则表示该结构设计满足了目标性能要求。
图使用基于位移的设计方法评价结构的用震性能
■分析方法
结构设计必须满足规范的一系列规定和要求,在完成满足规范要求的设计之后,结构的目标性能具体控制在哪个水准上,则由建筑物的使用者和设计老决定。
为了评价结构性能需要进行结构分析,基于性能的耐震设计方法一般有下列四种。
线性静力分析方法(LinearStaticProcedure.LSP)线性动力分析方法(LinearDynamicProcedure.LDP)非线性棘力分析方法(NonlinearStaticProcedure・NSP)非线性动力分析方法(NonlinearDynamicProcedure.NDP)
MIIWS/Gen中提供了上述四种分析方法,其中Pushover分析属于非线性前力分析方法。
Pushover分析又称为挣力弹塑性分析,是评价结构进入弹塑性状态后的结构极限状态和稳定状态的有效而简捷的方法。
该方法主要适用于低频结构彩响较大的结构中。
Pushover分析中可以考虑材料和几何非线性,材料非线性特性是通过定义构件裁面的荷载一位移的非线性待性实现的。
Pushover分析是通过逐渐加大预先设定的荷敦直到最大性能控制点位置,获得荷载一位移能力曲线(capacitycurve).多自由度的荷載一位移关系转换为使用单自由度体系的加逮度一位移方式表现的能力谱(capacityspectrum)t地爲作用的响应谱裟换为用ADRS(Acceleration-DisplacementResponseSpectrum)方式表现的需求谱(demandspectrum)<>通过比较两个谱曲线,评价结构在弹塑性状态下的最大需求内力和变形能力,通过与目标性能的比较,决定结构的性能水平(performancelevel)0
在HIDAS/Gen中使用ATC-40(1996)和FEMA-273(1997)中提供的能力谱法(CapacilySpectrumMethod.CSM)评价结构的耐震性能。
能力谱法(CSM)的原理如图所示0
(a)计算结构物的能力曲埃(capacitycurve)和能力谱(capacityspectruo)
S能力谱法(CapacitySpectrumMethod,CSH)的原理
Pushover分析是为了评价结构所拥有的耐爲性能,其前提条件是巳经完成了初步的分析和设计,即对于混凝土结构必须巳经完成了配筋设计。
Push""分析的优点如下:
可以评价结构进入弹塑性阶段的响应以及所拥有的抵抗能力
可以掌握结构的耗能能力和位移需求
可以掌握各构件屈服的顺序
对确定需要维修和加固的构件提供计算依据
■分析中适用的单元类型
MIDAS/Gen中Pushover分析中适用的单元类型有二维梁单元(2-dimensionalbeamelement)三维梁一柱单元(3-dimensionalbeam-columnelement).三维埴单元(3-dimensionalwallelement)N桁架单元(trusselement)<>各单元的特性如下。
二维梁单元和三维梁一柱单元
梁单元和梁一柱单元釆用的模型如图所示•其位移和荷载如下,其中适用于梁单元时无轴力项。
{〃}『={",M,„凡,Mvl,Ftl,Md,岛,M’2,坨2,M^,Fz2,M』
{"}'={"“,Qi,片I,G1,0】,Ql,“r2,%,Vv2,Q"徨2,
s二堆秦草元和三维梁一柱单元模型
三维墙单元模型
如图所示墻单元模型由中问的线单元,上下两端的刚性杆构成。
中问的线单元与三维梁一柱单元相同,刚性杆在XZ平面内做刚体运动。
S堆単元的节点力和节点位移
桁架单元模型
如图所示,桁架单元釆用轴向(x方向)的弹簧模型。
图析架卑元的节点力
■非线性弹黄的特性
在各单元模型中表现的弹簧并非表示弹簧的存在,而是表现分析的方法,即在弹簧位置将发生塑性变形。
弹簧具有的特性如下。
梁单元模型的弹黃特性用荷载一位移、轴力一单向弯矩一位移角.剪力一剪切变形、扭矩一扭转角等关系来表现
柱以及墙体单元模型的弹簧特性用荷载一位移轴力一双向弯矩一位移角、剪力一剪切变形.扭矩一扭转角等关系来表现
桁架单元模型的弹簧待性用荷载一位移关系来表现
单元的变形可用下面的各式来表现。
弯曲变形
节点上构件的变形角可用下列三项之和来表现。
在此,涉、歹、分别为弹性弯曲变形角.塑性弯曲变形角、因剪切产生的弯曲变形角。
另外,如图所示弯矩引起的塑性变形将假设集中在Q厶区段内。
图形中阴彩部分表示发生塑性变形的区段。
因此包含塑性变形和剪切变形的柔度矩阵(flexibilitymatrix)如下。
阳弯曲励度的分布假定
构件的荷载一位移关系可用柔度矩阵表现如下。
(4)
在此,
【/]=[/「+I/F+I/T⑸
塑性变形角、因剪切引起的
如图所示,式(5)中各项分别表示弹性弯曲变形角.弯曲变形角。
图芳矩一变形角关廉
轴向变形、扭转变形.剪切变形弹簧
在MIDAS/Gen的Pushover分析中假定轴力、扭矩、剪力在构件内不变,塑性较发生在构件中央。
其荷裁一位移关系可参照弯曲变形中的各式。
双向弯曲弹簧
双向受弯且受抽力作用时,先计算各向的屈服弯矩后使用下面关系式建立双向受弯相关公式。
上式适用与钢筋混礙土和钢材等所有构件。
■塑性较(plastichinge)特性
随着荷载的増加,结构构件将产生塑性钱,结构的刚度会发生变化,横向位務也将逐渐加大。
MIDAS/Gen中采用的塑性钱特性如下。
校待性:
多折线类型(Multi-LinearType)
-采用切向刚度矩阵(tangentstiffnessmatrix)
-荷载控制(loadcontrol)和位移控制(displacementcontrol)均可
-可考虑P-Delta效果
较特性:
FEMA类型(FEMAType)时
-割线刚度矩阵(secantstiffnessmatrix)
-釆用位移控制(displacementcontrol)
-可考虑P-Delta效果和大位移(largedeformation)效果
因为结构承受的荷载大小为巳知条件,所以一般采用荷裁控制方法。
荷载控制方法就是将荷载从零开始逐浙加载到极限荷載的方法。
位移控制是在基于性能的耐震设计中采用比较多的方法。
虽然不知道加载的荷载大小,但是可以通过预先设定满足目标性能的位移进行分析。
分析过程中可以获得荷载传递能力(load-carryingcapacity)和失稳(unstable)状态。
采用位移控制和割线刚度矩阵(secantstiffnessmatrix)时,在最大荷栽附件可以获得稳定的解。
多折线较类型
多折线钱类型可以用于荷载控制和位移控制方法中。
-荷栽一位移关系釆用双折线(Bilinear)和三折线(Trilinear)两中形式-屈服后刚度和抗裂刚度用与初始刚度的比值(stiffnessratio)来表现-能表现构件的刚度降低,但不能表现材料的强度降低
StifFncssafteryield
Bi-lincaryieldpoint
Tri-linear2ndyieldpoint
Stiffnessaftervield
Stiffnessaftercrack
Tri-linear
1styieldpoint(crack)
Initialstiffness
LateralDeformation
fS多折圾钱类型的集性钱转性
FEMA较类型
图FEMA较类型的建性较转性
FE5IA较类型是将钢筋混濮土构件和钢构件的循环加载试菽(reversedcyclicload)获得的资料理想化的结果,其特性如下。
MIDAS/Gen的FEMA钱待性只能使用位移控制方法。
点A位置:
未加载状态
AB区段:
具有构件的初始刚度(initialstiffness).由材料.构件尺寸、配筋率、边界条件、应力和变形水准决定。
点B位置:
公称屈服强度(nominalyieldstrength)状态
BC区段:
强度硬化(strainhardening)区段,刚度一般为初始刚度的5-10%,对相邻构件问的内力重分史有较大彩响。
点C位置:
由公称强度(nominalstrength)开始构件抵抗能力开始下降
CD区段:
构件的初始破坏(initialfailure)状态,報筋遏凝土构件的主筋断裂(fracture)^混凝土压碎(spalling)状态,钢构件的抗剪能力急剧下降区段。
DE区段:
残余抵抗(residualresilance)状态,公称强度的20%左右
-点E位置:
最大变形能力位置,无法继续承受重力荷载的状态。
■Pushover分析方法
MIDAS/Gen中提供两种Pushover分析方法,即基于荷载瓚分的荷载控制法和基于目标位移的位移控制法。
基于荷载增分的荷載控制法
MIDAS/Gen的荷栽控制法采用全牛顿一拉普森(FullNewton-Raphson)方法。
牛顿一拉普森方法是釆用徽分原理求解的方法,其优点是速度快。
釆用荷栽増的Pushover分析方法的图形接介绍如下。
荷載釆用具有一定分布模式的横向荷载。
荷载分布模式既可以采用地震荷载(QJ也可以釆用任意的荷载分布模式。
另外,也可以釆用包含节点荷截在内的用户定义的任何荷载工况。
(1)第1阶段:
计算弹性极限
首先使用用户定义的水平荷载计算构件的应力,然后计算各构件的应力与屈服应力的比值几。
将各构件的比值中的最小值乘以加载的荷载工况重新定义荷载。
入=(U
在此兄:
各较计算的屈服荷載系数(最大值
较的屈般应力
L:
初始应力
荷栽工况计算的钱的应力
(2)笫2阶段:
基于等差级数的增分分析
由弹性极限到预估的坍塌荷载(Q卅X)之问的荷载增量由下面的等差级数计算。
P=[(n+l)-/](PxXy^/等差增分步骤数⑻
在此乂:
第i步的荷載增量
d总荷栽
X:
预估的坍塌荷载与总荷载的比值(基本值为
(3)笫3阶段:
狈估坍塌荷載之后的荷載增分
使用最终计算的S+1)步菠的荷载増分。
(4)终止分析的条件
-到达最大増分步时
-层问位移角到达极限层问位移角时
-分析中计算的水平内力到达指定的大小时
-刚度矩阵为负(negative)时
基于目标位移的位移控制法
MIDAS/Gen的位移控制法是由用户定义目标位移,然后逐渐增加荷载直到达到目标位移的方法。
目标位移分为整体控制和主节点控制两种,整体控制是所有节点的位移都要满足用户输入最大位移,位移也定整体位移,不设置某一方向的位移控制。
主节点控制是用户指定待定节点的特定方向上的最大位移的方法。
基于性能的耐震设计大部分是先确定可能发生最大位移的节点和位移方向后给该节点设定目标位移的方法。
初始的目标位移一般可假定为结构总高度的1叙2%.4%o这些数值一般相当于最大层问位移值,于结构的破坏情况相关。
ATC-40或FEMA-273中将层问位移为1%时定义为直接居住水准(ImmediateOccupantLevel).2%时定义为生命安全水准(LifeSafetyLevel).4%时定义为坍塌预防水准(CollapsePreventionLevel)o这些值在构件级别上的意义可能会稍有不同。
作用荷栽
作用荷载应该采用能反映各层惯性力的横向荷栽。
所以推荐至少使用两种以上的横向荷载分布模式。
在MIDAS/Gen中提供了三种类型的荷载分布模式,即的力荷载分布模式、振型形状分布模式、与各层质呈成比例关系的等加速度分布模式。
釆用静力荷载分布模式时,用户可以定义任意形状的待力荷载分布。
采用振型形状荷载分布模式时必须先做特征值分析。
■基于性能点的耐農性能评价
在MIDAS/Gen中使用能力谱(CSM)原理评价结构的保有内力和耐震性能。
结构的保有
内力可通过Pushover分析获得的能力曲线和能力谱进行评价。
对地農作用的需求谱可以适用有效阻尼的弹性设计谱来评价。
将这两个谱表现在相同的坐标系上将获得意味着结构非线性最大雷求内力的交点,即性能点(performancepoint)o利用性能点位置的变形程度和保有内力来评价结构的耐蔑性能和性能水准。
能力谱和需求谱
评价结构的耐震性能和性能水准时会使用能力谱和需求谱的槪念。
通过Pushover分析将获得荷载一位移关系(V-U).响应谱也可荻得加速度一周期(A-T)的相关关系。
为了比较两个谱.需要将其該换为加遶度一位移(acceleration-displaccmentresponsespectrum,ADRS)。
图将荷栽-位移关廉转换为加速度-位移谱
图将加速度-局期谱转换为加速度-位移蔚
如图所示,荷载一位移关系转换为加速度一位移关系的方法如下:
(10)
在此,厂‘和AT为各自方向的斤阶振型的振型参与系数和有效质量系数,计算方法如下:
.V
振型参与系数厂'=——(11)
£巴尤
/-I
.V\2
乞叱J
振型参与质量(12)
式(9)和(10)为动力学理论的多自由度(MD0F)体系和单自由度(SD0F)体系之问
的关系。
即A和D为单自由度体系响应谱上的响应加速度和响应位移,V和U为多自由度体系的基底剪力和位移。
如图,弹性响应谱可以利用单自由度体系的位移和加速度关系式(13)进行转换
T1
D=(13)
4”
性能点(performancepoint)的评价
能力谱和笛求谱的交点称为性能点。
在MIDAS/Gen中提供的计算性能点的方法为ATC-40的能力(CSM)中提供的Procedure-A和Proccdu】p-B两种方法。
两种方法的基本原理相同,通过计算有效阻尼反复计算获得性能点的方法为Procedure-A方法,利用延性比和有效周期原理计算性能点的方法为Procedure-Bo
⑴计算等效阻尼(equivalentdamping)
在能力谱法(CS\D中,通过pushover分析获得能力谱后如下图所示使用具有相同面枳的双折线(biloinear)线来表现。
在CSM中使用具有5%阻尼的弹性响应谱和能力谱计算结构的等效阻尼。
因为结构的阻尼而耗散的能量等于双折线滞回曲线的面积,可按式(14)计算。
图裁用滞回曲践计算等效阻尼
(M)
]爲_63・7(%心,_〃卫严)r瓦讽—
在此.&=结构阻尼引起的耗散能
Ew=结构的最大变形能
将式(M)使用百分率的形式表现如下。
在此,0“表示阻尼比俚),在ATC-40中阻尼比超过25%时,需要谨慎的判断,且不许超过50尖。
⑵计算有效阻尼(effectivedamping)
地震作用作用下的钢筋混凝土结构的滞回曲线中没有表现刚度退化(stiffnessdegradation)和强度退化(strengthdeterioration).滑移或握憂(sliporpinching)的待性的理想化的滞回模型。
所以在ATC-40中为了反映钢筋混凝土的这些滞回待性,使用阻尼调整系数(dampingmodificationfactor)来调整等效阻尼。
调整后的等效阻尼称为有效阻尼系数,按下式计算。
式左侧的阻尼比5%为弹性体系的阻尼,不同结构特性的阻尼调整系数如下表。
结构待性
等效阻尼ft.(%)
阻尼调整系数(K)
类型A
(完全滞回转性)
<
>
)n0.51仏厶-心)
类型B
(一段滞回特性}
<25
>25
0.446仏心一
0.845,,,;
类型c
(弱滞回待性)
所有值
(3)非线性需求谱
使用前面计算的有效阻尼系数决定非线性咱应谱。
即利用有效阻尼系数计算响应谱的谱折减系数(spectrumreductionfactor,SR)。
如图所示加速度和速度的谱折减系数不同。
谱折减系数釆用了Newmark和Hall(1982)的地基运动扩大系数,加速度的谱折减系数(SRJ和速度速度的谱折减系数(SRJ的计算式如下。
根据结构的滞回特性,ATC-40中给出了谱折减系数的下限值。
項目
K
SRi
SR、
类型A(完全滞回特性)
类型B(一般滞回转性)
类型C(朝滞回转性)
表结构的滞回伶性对应的谦折截系数下限值
根据上述的计算过程可以获得设计地爲作用或线弹性反应谱对应的非线性需求谱。
将获得的非线性地爲需求谱和通过Pushover分析获得的结构的能力谱进行比较,可以获得结构的性能点。
(4)计算性能点
利用Pushover分析得到的结构的能力谱和非线性设计响应谱的比较,可以获得表现结构的非线性最大位移和保有内力的性能点,并且利用其来评价结构的性能水准。
■确定性能点的方法
MIDAS/Gen中根据能力谱(CSM)确定性能点的方法采用ATC-40中提供的两种方法。
其基本原理为使用有效阻尼系数评价需求谱并求其与能力谱的交点作为性能点。
Procedure-A
是ATC-40中提供的基本方法,首先将能力谱中斜率为初始刚度的切线和阻尼比为5%的弹性设计响应谱的交点作为初始的性能点。
然后确定初始性能点位置的等效阻尼,然后求使用有效阻尼系数的非线性设计响应谱,然后重新计算交叉点作为性能点。
重复上述过程,克到在使用有效阻尼系数的非线性设计响应谱和能力谱的的交点位置上位移响应和加速度响应的变化量在误差范因内,将此时的交点视为性能点。
采用Procedure-A方法确定性能点的方法参见图。
图使用Procedure-A方法计算性能点(ATC-40)
Procedure-B
ATC70中计算性能点的第二种方法是首先假设位移延性比,然后计算对应延性
比的结构的结构的有效周期,将有效周期直践和5%弹性设计响应谱的交点作为初始的性能点。
对弈于假定的位移延性比的放射线状的有效周期和非线性设计响应谱的交点将形成一个轨迹线,该轨迹线与结构的能力谱的交点为最终的性能点。
利用Procedure-B方法计算性能点的原理如图所示。
DemandCurvesforp^=5%,10%,15%,20%,25%and30%
SpectralDisplacement,inches
图刑用Procedure-B方法计算性能点(ATC-40)
该方法長首先假定位移薙性比,然后逐步计算有效阻尼系数,所以在交叉点计算的响应误差发敬的槪率较低。
前面介绍的Procedure-A方法在寻找性能点的过程当中收敛性不是很好,而Procedure-B方法不仅收敛性能好,而且不必建立针对多个阻尼比的弹性反应谱,而是根据变化的阻尼比和振动周期获得响应谱的轨迹即可获得性能点,所以Procedure-B方法是相对比较简单的方法。
在MIDAS/Gen中提供的两种方法的操作界面如下图所示。
图科用Procedure-A的方法决定性館点(MIDAS/Gen)
X利用ProcedurQ-B的方法决戾性能点(MIDAS/Gen)
■性能评价
如果结构的位務在目标性能范国内,則继续评价各构件的性能。
在MIDAS/Gen中采用与FEMA-273或ATC-40中推荐的方法类似的方法评价构件的性能。
如图所示性能状态分为三个阶段。
1()=直接居住极限状态(ImmediateOccupancy)
LS=安全极限状态(LifeSafety)
CP=坍塌防止极限状态(CollapsePrevention)
图构件的性能评价
■Pushover分析过程
1.棘力分析及完成设计
在做Pushover分析之前首先要对结构进行一般的静力分析及设计。
2.输入Pushover分析控制数据
在设计〉屜汾""夕祈控村灯话框中输入pushover分析的最大迭代计算步骤数、各步骤内迭代计算次数和收敛条件。
3.定义Pushover荷载.工况
杜设计>Pushover荷栽工况试话框中,输入Pushover分析前的初始荷载和Pushover荷载工况。
首先选择荷我控制或位移控制。
作为初始荷栽输入自重‘Pushover荷载工况可选靜力荷载工况、等加速度、振型,各荷载类型之问也可进行组合。
4・定义校数据
隹设计>定义校待性值对话框中输入反映材料非线性特性的较数据。
在MIDAS/Gen中提供多折线和FEMA两种类型的较特性数据。
除了默认的特性,用户也可以自定义较特性。
5.分配钱给构件
在设计>分配较轄性值戏话框中选择要分配的较特性值并分配给适当的单元。
一般来说给梁单元分配弯矩钱,给柱单元和剪力埴单元分配PH或PMM絞,给桁架单元分配轴力较。
6.运行Push