《Part3-第14章-大跨度桥梁的稳定理论课件.ppt》由会员分享，可在线阅读，更多相关《Part3-第14章-大跨度桥梁的稳定理论课件.ppt（54页珍藏版）》请在三一办公上搜索。
1、第十二章,大跨度桥梁,稳定理论,同济大学桥梁工程系,大跨度桥梁研究室,第十二章,大跨度桥梁的稳定理论,本章主要内容,1,概,述,2,第一类弹性及弹塑性稳定分析,3,拱桥稳定分析和非保向力效应,4,材料非线性问题,5,第二类稳定问题和极限承载力全过程分析,6,小,结,1.,概述,1.1,稳定理论的发展,?,什么是结构失稳？,结构在外力增加到某一量值时，稳定性平衡状态开始丧失，稍有扰动，结构,变形迅速增大，使结构失去正常工作能力的现象,?,稳定问题的重要性,随着桥梁跨径的不断增大，桥塔高耸化、箱梁薄壁化以及高强材料的应用，,结构整体和局部的刚度下降，使得稳定问题显得比以往更为重要,?,桥梁结构的失
2、稳形态,桥梁结构的失稳现象表现为结构的整体失稳或局部失稳,局部失稳,是指部分子结构的失稳或个别构件的失稳，局部失稳常常导致整个,结构体系的失稳,1.1,稳定理论的发展,(,续,),?,桥梁失稳事故的发生促进了,桥梁稳定理论的发展,?,1744,年，欧拉,(L.Eular),就提出了,压杆稳定,的著名公式,?,彭加瑞,(A.Poincare,1885),明确了稳定概念，并推广到流体力学,的层流稳定问题中，即,稳定分支点,的概念,?,恩格塞,(Engesser),和卡门,(Karman),等根据大量中长压杆在压曲,前已超出弹性极限的事实，分别提出了,切线模量理论,和,折算模量,理论,?,普兰特尔和
3、米歇尔几乎同时发表了关于,梁侧倾问题,的研究成果,1.1,稳定理论的发展,(,续,),?,薄壁轻型结构的使用，提出了稳定新课题,?,瓦格纳,(H.Wagner,1929),及符拉索夫,(1940),等建立关于,薄壁杆件的,弯扭失稳理论,?,证明其临界荷载值大大低于欧拉理论值，且不能用分支点的概念来解释,?,引入了极值点失稳的观点以及跳跃现象的稳定理论,?,稳定理论与非线性理论的联系密不可分,?,只有通过对结构,几何非线性关系,以及,材料非线性本构关系,的研究，,才能深入揭示复杂稳定问题的实质,1.2,两类稳定问题,?,研究结构稳定问题的两种形式,1),第一类稳定：,分支点失稳,从小范围内观察，
4、以小位移理论为基础,2),第二类稳定：,极值点失稳,从大范围内研究，以大位移非线性理论为基础,?,由于第一类稳定问题是特征值问题，求解方便，在许多,情况下两类问题的临界值又相差不大，因此,研究第一类稳,定问题仍有着重要的工程意义,1.3,稳定问题的求解方法简介,?,静力平衡法,?,从平衡状态来研究压杆屈曲特征，即研究载荷达到多大时，弹,性系统可以发生不同的平衡状态,?,实质是求解弹性系统的平衡路径,(,曲线,),的分支点所对应的载荷,值,(,临界载荷,),?,能量法,?,求弹性系统的总势能不再是正定时的载荷值,1.3,稳定问题的求解方法简介,?,振动法,?,当压杆在给定的压力下，受到一定的初始
5、扰动之后，必将产,生自由振动,?,如果振动随时间的增加是收敛的，则压杆是稳定的,?,缺陷法,?,由于缺陷的影响，杆件开始受力时即产生弯曲变形,?,在一般条件下缺陷总是很小的，弯曲变形并不显著,?,当荷载接近临界值时，变形才迅速增大，由此确定失稳条件,1.3,稳定问题的求解方法简介,(,续,),?,对于欧拉压杆而言，所得到的临界荷载值是相同的,?,但它们的结论并不完全一样，表现在以下几个方面,(1),静力平衡法,?,当,P=P,1,、,P,2,.P,n,时压杆可能发生屈曲现象，无法判断何种情,况最可能失稳,?,在,P,?,P,1,、,P,2,.P,n,时，屈曲的变形形式不能平衡，无法回答直,线形
6、式的平衡是否稳定的问题,1.3,稳定问题的求解方法简介,(,续,),(2),缺陷法,?,当,P=P,1,、,P,2,.P,n,，杆件将发生无限变形,?,但对于,P,在,P,1,、,P,2,.P,n,各值之间时压杆是否稳定的问题也,不能解释,(3),能量法和振动法,?,PP,1,之后不论,P,值多大，压杆直线形式的平衡都是不稳定的,?,和事实完全一致,1.3,稳定问题的求解方法简介,(,续,),?,由于桥梁结构的复杂性，不可能单靠上述方法来解决,其稳定问题,?,大量使用的是,近似求解方法,：,?,从微分方程出发，通过数学上的各种近似方法求解,如逐次渐近法,?,基于能量变分原理的近似法,如,Rit
7、z,法，有限元方法可以看成是,Ritz,法的特殊形式,2.,第一类弹性及弹塑性稳定分析,2.1,第一类稳定问题的线弹性有限元分析,?,在发生第一类失稳前，结构在初始构形线性平衡，大位移,矩阵,0,K,L,为零,?,不论,T.L,还是,U.L,列式，表达形式是统一的,(,?,K,?,?,?,K,?,?,),?,?,u,?,?,?,?,R,?,?,按线性代数理论，必有：,(12-3),?,在结构处在临界状态下，即使,R0，,u,也有非零解,?,K,?,?,?,K,?,?,?,0,(12-4),2.1,第一类稳定问题的线弹性有限元分析,(,续,),?,发生第一类失稳前满足线性假设，应力与外荷载以及,
8、几何刚度为线性关系,?,若某种参考荷载,?,P,?,对应的结构几何刚度阵为,?,K,?,?,K,?,?,式,(12,4),可写成,?,?,K,?,?,?,?,(12-5),?,K,?,?,?,?,K,?,?,?,0,(12-6),?,稳定问题转化为求方程的,最小特征值问题,2.1,第一类稳定问题的线弹性有限元分析,(,续,),?,K,可以分成一期恒载的初内力刚度阵,K,1,?,和后期荷载（,二期、活载等）的初内力刚度阵,K,2,?,两部分,?,计算一期恒载稳定问题，,K,2,?,?,0,，,?,为,恒载稳定安全系数,?,计算后期荷载稳定问题，则恒载,K,1,?,可近似为一常量，式,(12,6)
9、,改写成：,K,?,K,1,?,?,K,?,0,?,2,?,(12-7),?,?,为,后期恒载稳定安全系数，相应的特征向量就是失稳模态,3.,拱桥稳定分析和非保向力效应,本节以,解析法,来阐述拱桥的第一类稳定计算,可分为以下两类问题：,?,面内稳定,?,侧向稳定,3.1,圆弧拱平面屈曲微分方程,圆弧拱的屈曲易于获得解析解，得到端支承情况和矢跨比对拱桥临界,荷载值的影响,图,12.2,均布径向荷载作用下的圆弧拱,?,在均布径向荷载,q,作用下，开始只有沿拱轴方向的弹性压缩变形,?,若忽略轴向变形的影响，拱轴线与压力线完全吻合，处于无弯矩,状态,3.1,圆弧拱平面屈曲微分方程,(,续,),当荷载达
10、到临界值时，拱发生微小的弯曲变形,v,，且在截,面上存在弯矩,M,，在这一变形状态下可以导出它的屈曲微分方,程为：,d,v,v,M,2,?,2,?,?,EI,x,ds,R,或：,2,2,2,(12,9),d,v,MR,?,v,?,?,(12,10),2,EI,x,d,?,3.2,等截面圆弧拱在均布径向荷载作用下的屈曲临界荷载,双铰圆弧拱在径向荷载,q,作用下,(,图,12.3),，其拱截面弯矩,M,?,N,?,v,?,qRv,(12,11),代入式,(12,10),，即得,d,2,v,2,?,k,v,?,0,(12,12),2,d,?,图,12.3,受均布径向荷载的等截面双铰圆弧拱,3.2,等
11、截面圆弧拱在均布径向荷载作用下的屈曲临界荷载,(,续,),边界条件,:,?,?,0,v,?,0,得,c,2,?,0,?,?,2,?,v,?,0,得,c,1,s,in,2,k,?,?,0,c,1,不能为零，则必须有,s,in2,k,?,?,0,(12,15),由此得到,2,k,?,?,n,?,(,n,?,1,2,3,.),由于拱的两端不能移动，圆弧拱轴也假定不发生伸缩。因此,n=1,相应的失稳模态是没有意义的。要求最小特征值时,n=2,，拱,的屈曲模态为,:,?,v,(,?,),?,c,1,sin,?,(12,16),?,3.2,等截面圆弧拱在均布径向荷载作用下的屈曲临界荷载,(,续,),临界荷
12、载值为：,EI,X,?,EI,x,(,2,?,1,),?,K,1,3,(12,17),q,cr,?,3,R,?,R,式中,2,?,K,1,?,2,?,1,(12,18),?,K,1,称为拱的临界荷载系数,(,或稳定系数,),，与夹角,有关。,式,(12,17),也可写成中心受压直杆的欧拉公式的标准形式,2,?,EI,X,?,EI,x,?,N,cr,?,q,cr,R,?,2,2,(,1,?,2,),?,(12,19),2,R,?,?,s,0,2,2,2,3.3,圆拱的面外稳定,平面拱轴侧倾后是一条空间曲线，其位移与几何关系用曲线坐标来描述。,图,12.4,侧倾变形后的拱,拱侧倾变形后,(,图,1
13、2.4),，任意截面,s,在垂直于拱平面,x,轴，指向拱轴法向的,y,轴和同拱轴切线重合的,z,轴三个方向分别发生了线位移,u,、,v,、,w,并绕这三个轴,发生了转角位移,、,、,。截面主轴,x,、,y,、,z,也随着拱的侧倾产生了变位。研,究相距,ds,截面的变形，可得拱绕,y,、,z,轴转动的曲率关系：,3.3,圆拱的面外稳定,(,续,),?,d,u,?,?,y,?,?,2,?,R,dS,?,?,(12,21),d,?,1,du,?,?,z,?,?,dS,R,dS,?,?,下面用能量法研究两端固结拱轴线长度为,L,的园拱的侧向稳定问题。,圆拱侧倾时，拱肋侧向弯曲变形能为：,V,B,?,2
14、,EI,y,2,?,L,/,2,?,L,/,2,?,ds,(12,22),2,y,拱肋扭转变形能为：,GJ,L,/,2,2,?,Z,ds,(12,23),V,T,?,?,2,?,L,/,2,拱肋轴力在侧倾时所作外力功为：,du,2,V,D,?,?,qR,?,?,L,/,2,(,),ds,(12,24),ds,L,/,2,3.3,圆拱的面外稳定,(,续,),结构势能为：,?,?,V,B,?,V,T,?,V,D,(12,25),设失稳模态为：,2,?,S,?,u,a,?,A,(,1,?,cos,),?,L,?,(12,26),2,?,S,?,a,?,B,(,1,?,cos,),?,L,?,将式,(
15、12,21),?,(12,24),、,(12,26),代入式,(12,25),，由,?,?,0,?,A,2,?,?,0,(12,27),?,B,2,2,EI,y,?,4,?,?,?,?,易得：,q,cr,?,3,2,2,2,(12,28),R,?,(,4,?,?,?,),3.3,圆拱的面外稳定,(,续,),EI,y,其中：,?,?,为弯、扭刚度比例系数。,GJ,当式,(12,28),确定的,q,cr,比相应面内失稳临界荷载,为小时，圆拱先出现侧倾失稳,对于宽跨比较小的拱桥，,侧向刚度相对较小和单承,重面拱桥，都有可能发生侧倾弯扭失稳，在设计时必须,对这类结构进行侧稳验算,3.4,拱桥稳定与非保
16、向力效应,1),拱桥的稳定问题,实际拱桥都带有拱上建筑，拱圈轴线形状各异，都受有竖向荷载，,研究拱桥的稳定问题可以抓住以下三点：,a),确定拱轴线时，力求拱肋受力以轴向受压为主，因此可将各类,拱肋的稳定问题通过一定的等代关系转换成圆拱的稳定问题来研究；,b),拱上建筑多以连续梁为主，梁的刚度增加了拱的稳定性，用能,量法计算这类结构的稳定性比较方便；,c),用解析法计算拱桥稳定问题比较复杂，一般都采用数值计算，,可以从大量数值计算结果的规律中总结出常用拱桥的稳定计算近似公,式。,3.4,拱桥稳定与非保向力效应,(,续,),例如拱桥的立柱刚度远比拱圈和梁的刚度小，为简化计算，可以,假定各立柱上下端
17、均系铰结,通过数值计算，可把这种简化结构的临界荷载近似地写成：,EI,a,?,q,cr,?,K,3,?,l,?,(12,29),2,?,?,f,?,f,?,?,I,b,?,?,?,?,?,K,?,K,?,1,?,?,0,.,95,?,0,.,7,?,?,?,?,?,?,l,l,I,?,?,?,?,?,?,a,?,?,?,?,?,式中：,K,为只有拱肋时的临界荷载系数；,I,b,加劲梁的抗弯刚度；,I,a,拱平面抗弯刚度。,3.4,拱桥稳定与非保向力效应,(,续,),2),非保向力效应,图,12.5a),是,上承式拱桥的面内失稳,，由于桥面抗弯刚度较小，当拱发生失,稳时，立柱受到梁施加的水平约束
18、而变成倾斜，产生的水平分力有,加速,其侧倾,的趋势,图,12.5b),是,系杆拱侧倾失稳,，吊杆受到梁施加的水平约束而变成倾斜，产,生的水平分力有,减缓,其发生失稳的趋势,图,12.5,非保向力系对拱稳定的影响,3.4,拱桥稳定与非保向力效应,(,续,),?,跟随结构变形而改变其方向，这种力系称为,非保向力系,?,非保向力系对结构稳定性有正面效应，也有负面效应,下面以单承重面系杆拱为例来讨论非保向力系对其侧向稳定的影响,作用均布荷载,q,，设吊杆的布置满足膜张力假定，则吊杆拉力为：,T,qa(12,31),式中：,a,?,吊杆间距。,3.4,拱桥稳定与非保向力效应,(,续,),拱肋侧倾后，吊杆
19、发生倾斜，其拉力,T,对桥面产生了一个向外的水平,分力，使之发生侧向弯曲变形,u,b,(x),，而对拱肋产生了一个向内的水平分力,H(x),，这个恢复力就是,非保向力效应,，相当于一个侧向水平弹簧支承效应。,u,a,?,u,b,H,(,x,),?,T,?,k,(,x,),u,a,(12,32),y,(,x,),qa,u,a,?,u,b,式中：,k,(,x,),?,(12,33),y,(,x,),u,a,考虑到桥面侧向刚度相对于拱肋要大得多，近似地取,EI,by,=,?,，则式,(12,33),简化成,qa,k,(,x,),?,(12,34),y,(,x,),3.4,拱桥稳定与非保向力效应,(,
20、续,),当系杆拱发生侧倾时，其总势能除了前面式,(12,22),?,(12,24),列出的三项,外，还增加了考虑非保向力效应的虚拟弹簧支承变形能,V,k,：,1,k,(,x,),2,u,a,ds,(12,35),V,k,?,?,2,a,将式,(12,35),增添到式,(12,25),中，由能量驻值原理可得系杆拱侧倾临界,荷载：,q,cra,L,2,L,?,2,?,?,非保向力效应系数,式中：,1,?,?,q,cr,?,q,cr,(12,36),1,?,C,q,cr,由式,(12-28),给出。,V,k,C,?,?,V,D,(12,37),对圆弧拱，偏安全地取,y(x)=f,，则,C,的下限为：
21、,3,?,2,R,C,?,(,),(12,38),4,?,f,3.4,拱桥稳定与非保向力效应,(,续,),根据不同的矢跨比,f/l,，可算得,c,和,?,值，列于表,12,1,。,非保向力效应系数一般在,2.5,?,3.5,之间。,不同矢跨比,f/l,时的,c,和,?,值,表,12,1,f/l,C,近似,?,值,精确,?,值,1/3,1/4,1/5,1/6,1/7,1/8,0.684,0.653,0.638,0.630,0.623,0.621,3.16,3.50,2.88,3.31,2.76,3.22,2.70,3.17,2.65,3.14,2.64,3.12,4.,材料非线性问题,4.1,概
22、,述,?,当构件应力超过弹性极限后，材料弹性模量,E,成为,应力的函数，导致基本控制方程的非线性，即材料,非线性问题,?,凡是在本构关系中放弃材料线性关系假定的理论，,均属材料非线性范畴,?,桥梁结构以钢和砼作为主要建材，因此涉及的材料,非线性主要是非线性弹塑性问题和砼徐变问题,4.2,弹塑性应力、应变关系与屈服准则,根据实验结果，单轴应力下材料的应力、应变关系如,图,12.7,所示，可归结为如下几点：,1),应力在达到比例,极限前，材料为线弹,性；应力在比例极限,和弹性极限之间，材,料为非线性弹性。,图,12-7,单轴应力下材料的应力、应变关系,4.2,弹塑性应力、应变关系与,屈,服准则,(
23、,续,),2),应力超过屈服点，材料应变中出现不可恢复的塑性,应变,:,?,?,?,e,?,?,p,(12-39),应力和应变间为非线性关系：,?,?,?,(,?,),(12-40),3),应力在某一应力下卸载，则应力增量与应变增量之,间存在线性关系，即：,d,?,?,Ed,?,为了判断是加载还是卸载，用如下加载准则：,(12-41),当,?,d,?,?,0,时为加载，满足,(12-40),当,?,d,?,?,0,时为卸载，满足,(12-41),4.2,弹塑性应力、应变关系与,屈,服准则,(,续,),4),在卸载后某应力,?,下重新加载，则：,?,?,?,0,时，,d,?,?,Ed,?,(12
24、-42),?,0,为卸载前材料曾经受到过的最大应力值，称后屈服应,力，若：,?,0=,?,s,材料称为理想塑性的；,?,0,?,s,称材料为硬化的。,5),从卸载转入反向力加载，应力、应变关系继续依式,(12-41),或,(12-42),，一直到反向屈服。在复杂应力状态,下，判断材料是否屈服，可以用应力的某种函数表示：,4.2,弹塑性应力、应变关系与,屈,服准则,(,续,),F,(,?,ij,),?,0,(12-43),若以,?,ij,为坐标轴建立一坐标空间，则式,(12-43),的几,何意义为空间超曲面。任一应力状态在此空间中代表一,个点，当此点落在屈服面之内时：,F,(,?,ij,),?,
25、0,，材料呈弹,性状态；,F,(,?,ij,),?,0,时，材料开始进入塑性。,各向同性材料的屈服条件与坐标轴选取无关，屈服函,数常以主应力函数形式表示：,F,(,?,1,?,2,?,3,),?,0,(12-44),4.2,弹塑性应力、应变关系与,屈,服准则,(,续,),常用的屈服条件有,：,?,屈雷斯卡,(Tresca),屈服条件：假定最大剪应力达,到某一极限值时，材料开始屈服，相当于材料力学,中的第三强度理论,?,密赛斯,(Von Mises),屈服条件：假定偏应力张量,的第二不变量达到某一极限时，材料开始屈服,相,当于材料力学中的第四强度理论,此外还有,Drucker-Prager,屈服
26、准则,Zienkiewicz-Pande,屈服准则等,4.3,弹塑性本构矩阵的增量表达式,设屈服函数用下式表示：,F,(,?,ij,K,),?,0,式中：,?,ij,应力状态；,K,硬化函数。,(12-45),在增量理论中，把材料达到屈服以后的应变增量,分为弹性增量和塑性增量两部分，即：,d,?,?,d,?,?,d,?,e,p,(12-46),4.3,弹塑性本构矩阵的增量表达式,(,续,),其中弹性应变增量部分与应力增量之间仍服从虎,克定律，即：,d,?,?,D,e,d,?,其中：,D,e,为弹性矩阵。,e,(12-47),塑性变形不是唯一确定的，对应于同一应力增量,，可以有不同的塑性变形增量
27、。若采用相关联的流,动法则，塑性变形大小虽然不能断定，但其流动方,向与屈服面正交。用数学公式表示这一假定，即可,得：,?,?,F,?,p,(12-48),d,?,?,?,?,?,?,?,?,4.3,弹塑性本构矩阵的增量表达式,(,续,),将,(12-47),、,(12-48),式代入,(12-46),式，则可得：,?,?,F,?,d,?,?,D,e,d,?,?,?,?,?,?,?,?,对式,(12-45),全微分得：,?,1,(12-49),?,F,?,F,?,F,dF,?,d,?,1,?,d,?,2,?,?,?,?,?,?,?,dK,?,0,?,1,?,2,?,K,(12-50),或,?,?
28、,F,?,?,?,d,?,?,A,?,?,0,?,?,?,T,(12-51),?,F,1,其中：,A,?,?,dK,?,K,?,(12-52),4.3,弹塑性本构矩阵的增量表达式,(,续,),?,?,F,?,将,?,d,?,?,D,e,前乘,(12-49),式，并利用,(12-51),式消去,?,?,?,T,可得：,T,?,?,?,F,?,F,?,?,F,?,?,?,?,?,A,?,?,?,D,e,?,?,?,?,?,D,e,d,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,?,T,(12-53),?,?,F,?,?,?,D,e,?,?,?,由此可得：,?,?,d,?,T,?
29、,?,F,?,?,?,F,?,A,?,?,?,D,e,?,?,?,?,?,?,?,?,T,(12-54),用,D,e,前乘,(12-49),式，移项后得,?,?,F,?,d,?,?,D,e,d,?,?,D,e,?,?,?,?,?,?,(12-55),4.3,弹塑性本构矩阵的增量表达式,(,续,),将,(12-54),式代入,(12-55),式，即可得：,T,?,?,?,?,F,?,?,?,F,?,D,e,?,?,?,?,?,D,e,?,?,?,?,?,?,?,?,?,d,?,?,(,D,?,D,),d,?,?,D,d,?,d,?,?,D,e,?,e,P,ep,T,?,?,?,?,F,?,?,?
30、,F,?,A,?,?,?,?,D,e,?,?,?,?,?,?,?,?,?,?,?,?,?,(12-56),T,?,?,?,?,F,?,?,?,F,?,D,e,?,?,?,?,D,e,?,?,?,?,?,?,?,?,?,?,D,?,D,?,其中,:,ep,?,e,T,?,F,?,?,F,?,?,?,?,?,A,?,?,?,D,e,?,?,?,?,?,?,?,?,?,?,?,?,?,(12-57),此即为,增量理论的弹塑性矩阵通式,。其具体的数,学表达式将由曲服函数确定。,4.4,弹塑性问题的有限元法,在弹塑性增量理论中，讨论仍限于小变形情况。,其应变位移几何运动方程和平衡方程相同于线性,问题，不
31、需要作任何变动。,需要改变的只是在塑性,区范围内用塑性材料的本构关系矩阵,D,ep,代替原来,的弹性系数矩阵,D,e,。,因此，可直接得到弹塑性分析,有限元平衡方程：,K,T,?,u,?,?,R,式中,:,t,K,T,?,?,?,B,T,D,ep,B,dv,v,t,t,t,t,t,t,t,t,(12-58),(12-59),?,R,?,?,F,?,?,T,?,?,F,c,?,?,F,I,(12-60),4.4,弹塑性问题的有限元法,(,续,),t,t,?,T,?,F,其中,和,分别表示与结构面荷载,t,及体荷载,f,对,?,t,F,c,应的等效节点力增量；,为节点集中外荷载增量；,t,?,F,
32、I,为初应力或初应变增量引起的外荷载增量，它,?,t,至,t,时间的增量为：,们在,t-,?,F,?,?,?,N,?,f,dv,t,T,t,v,(12-61),?,T,?,?,?,N,?,t,ds,t,T,t,v,(12-62),对于初应力问题：,t,T,?,F,I,?,?,?,B,?,?,I,dv,v,(12-63),对于初应变问题：,t,?,F,I,?,?,?,B,D,e,?,?,I,dv,(12-64),T,v,5.,第二类稳定和极限承载力全过程分析,?,传统的“强度设计”以构件最大工作应力乘以安全系数等于材,料的屈服应力为依据；,?,一般情况下，构件某截面开始屈服并不能代表结构完全破坏
33、，,结构所能承受的荷载通常较构件开始屈服时的荷载为大；,?,桥梁结构的极限承载力是指桥梁承受外荷载的最大能力；,?,可以准确地知道桥梁结构在给定荷载下的安全贮备或超载能力,，为其安全施工和营运管理提供依据和保障；,5.,第二类稳定和极限承载力全过程分析,(,续,),?,全过程分析是用于桥梁结构极限承载力分析的一种计算方,法，它通过逐级增加工作荷载集度来考察结构的变形和受力,特征，一直计算至结构发生破坏；,?,桥梁结构在不断增加的外载作用下，结构刚度发生不断变,化，当外载产生的应力使得结构切线刚度阵趋于奇异时，结,构承载能力就到达了极限，此时的外荷载即为极限荷载。,5.1,非线性方程的求解问题,
34、?,一般结构的结构刚度阵在,p-,?,曲线上升段是正定的，在下,降段为负定的；,?,进行“全过程”分析过程中，当荷载接近极限值时，很小,的荷载增量都会引起很大的位移，可能还未找到极限荷载,就出现了求解失效现象；,?,为了找到真实的极限荷载，克服下降段的不稳定现象，各,国学者提出了许多算法，下面就常用的两种方法作一介绍,5.1,非线性方程的求解问题,(,续,),1),逐步搜索法,对于只要求出极值荷载，而对,P-,?,下降段不感趣的情况，,可采用逐步搜索顶点的算法，其基本思想是：,加一荷载增量,?,P,，计算发散后，退回上级荷载状态并改,用荷载步长,?,P/2,；,若计算收敛，则再加一级荷载为,?
35、,P/4,；,若加,?,P/4,后计算发散，则再改用荷载步,长为,?,P/8,如此搜索，若原步长,?,P,预计为,5%,的破坏荷载，则,?,P/4,已,接近,1%,的极限荷载，对桥梁结构来说，已可满足精度要求。,当然还可向前再搜索一步到,?,P/8,。,5.1,非线性方程的求解问题,(,续,),2),位移控制法,如果在分析过程中不是控制荷载增量而是控制位移,增量，则,P-,?,曲线的下降段部分便不难求得。,对于一般结构，我们可将刚度矩阵重新排列，使得,要控制的位移,(,例如,=,?,u,2,),排到最后一项，同时将原刚,度矩阵分块，其有限元方程变为：,?,K,1,1,?,K,?,2,1,K,1
36、,2,?,?,?,u,1,?,?,P,?,R,1,?,1,?,?,?,?,?,?,?,?,?,?,?,?,K,2,2,?,?,?,u,2,?,2,?,?,P,?,R,2,?,(12-88),式中：,P,1,P,2,T,?,参考荷载向量；,?,?,控制荷载的步长系数；,R,1,R,2,T,?,求解迭代过程中的不平衡力向量,。,5.1,非线性方程的求解问题,(,续,),改写方程,(12-88),为：,?,K,11,?,P,1,?,?,?,u,1,?,?,R,1,?,?,K,12,?,?,?,?,?,?,u,2,?,K,?,?,?,?,2,?,?,?,?,?,?,21,?,P,?,R,2,?,?,K
37、,22,?,(12-89),这样，求解方程时可控制指定的值，求出相应的位移,?,u,1,及荷,载增量比例因子,?,。由于,K,ij,与位移有关，求解时需要迭代，使,得,R,1,R,2,T,值趋于零，以满足精度要求。,?,K,11,?,P,1,?,需要指出，方程,(12-89),中的系数矩阵,是不对称，,?,K,?,?,P,21,2,?,?,也不呈带状，求解时需要的存储单元较多，这是该方法的一大缺,点。,5.2,单元模式与破坏形态的选取,?,用于极限承载力分析的单元模式选取,桥梁结构分析以梁单元为主，用于极限承载力分析的,梁单元模式主要有三种：,?,带有塑性铰的一般梁单元,?,不分层的等参梁单元
38、,?,分层梁元,5.2,单元模式与破坏形态的选取,(,续,),?,破坏形态的模拟问题,下面以塑性铰法为例，说明桥梁结构的极限承载力计,算步骤：,1),确定成桥状态的内力与构形；,2),以成桥状态为初态，用单位计算荷载向量,p,进,行结构分析。根据计算结果和极限弯矩，估算第一,个塑性塑性铰出现时的荷载增量倍数；,3),以,?,P,1,p,作用于结构，按全非线性进行结,构分析，迭代形成第一个塑性铰和实际的荷载增量,倍数,?,1,；,5.2,单元模式与破坏形态的选取,(,续,),4),检验结构是否成为机构，若是，给出极限荷载，,计算结束。否则，估算出现下一个塑性塑性铰时的荷,载增量倍数；,5),以上
39、次计算结束时的结构状态为初态，以,?,P,i,p,作用于结构，按全非线性进行结构分析，迭代,形成第,i,个塑性铰和实际的荷载增量倍数,?,；,6),重复,4),5),的计算，直至第,n,个塑性铰出现时结,构成为机构。此时，结构的极限荷载为,:,P,p,?,?,i,?,1,n,?,i,p,(12-90),6.,小,结,本章首先介绍了桥梁结构稳定问题的分类及稳定理,论的的发展历程。,根据稳定与平衡的关系，建立了求解第一类稳定问,题的控制方程，并介绍了桥梁结构第一类稳定问题的,求解方法。,通过讨论圆弧拱稳定问题的解析解，给出了计算拱,桥稳定和考虑非保向力效应的实用方法。,通过讨论材料非线性问题，引出了桥梁极限承载力,及其全过程分析问题，指出只有通过对结构几何非线,性关系以及材料非线性本构关系的研究，才能深入揭,示复杂稳定问题的实质。,
