ANSYS在预应力钢筋混凝土结构非线性分析中的应用张耀庭文档格式.docx

1、体外预应力中图分类号:T U 378文献标识码:A 文章编号:1000-5730（200304-0020-04有限元法1,2是在连续体力学领域应用的一种数值分析方法.ANSYS 2是由美国ANSYS 公司开发的用于微机平台的大型有限元分析软件.其功能强大,已扩展到航空、航天、电子、汽车、土木工程等各领域35.近年来,AN SYS 逐步在土木工程6-8中得到应用.将AN SYS 应用于钢筋混凝土结构的分析时,在弹性及线弹性阶段一般可取得满意的结果,但考虑材料的非线性时,结果出入较大.因此,作者探讨用AN SYS 对钢筋混凝土结构进行非线性分析的方法.1材料单元及本构关系的选择结构的非线性包括材料

2、非线性、几何非线性和状态非线性.由于钢筋和混凝土都是弹塑性材料,因此混凝土结构非线性分析中最常见的是材料的非线性.AN SY S 中有很多种单元类型,但只有Solid 65能模拟混凝土的性质,因此采用8结点的3D-So lid65单元可以进行塑性、蠕变、应力刚化、大变形和大应变的分析.钢筋可以采用3D-Link 8单元来模拟.混凝土是弹塑性材料,因此需考虑其非线形.当定义了tb,concr 及matnum 后,仅仅是定义了混凝土的W -W 破坏准则,混凝土在开裂和压碎前均为线性的应力应变关系,开裂和压碎后采用程序给出的本构关系.非线性分析中也可以通过tb 和miso 输入混凝土的应力应变关系来

3、确定本构关系.单向应力-应变曲线常见的有双折线、二次曲线+水平直线段及二次曲线+下降直线段等,此外还有在试验中实测的应力-应变曲线.在ANSYS 计算中可以根据实际情况采用合适的曲线,为便于程序收敛可采用二次曲线+水平直线段的形式,其应力应变关系为0X X 0e =f c 2X /X 0-（X /X 02;（1X0X X cu e =f c ,（2式中,X 为混凝土的应变;X 0为相应于棱柱体抗压强度时的混凝土压应变,取0.002;X cu 为混凝土的极限压应变,取0.0033;e 为混凝土的应力;f c 为混凝土棱柱体的抗压强度设计值.其应力应变关系如图 1.图1混凝土本构关系在分析计算中,

4、通常把钢筋应力-应变曲线简化成双折线形式,因此在AN SYS 分析中使用经典的双线性随动强化（BKIN ,其包括弹性斜率和塑性斜率.由于使用了随动强化的Vo n Mises 屈服准则,所以有包辛格效应.第20卷第4期2003年12月华中科技大学学报（城市科学版J .o f HU ST.（U r ba n Science Editio nV o l.20No.4Dec.20032建模过程根据钢筋的处理方式,有限元模型主要有分离式、分布式和组合式.在建模过程中,钢筋的处理方法一般有两种.一种是体分割法,即用工作平面根据钢筋的位置不断地将体分割,分割后体上的一条线定义为力筋线,最终形成许多复杂的体和

5、多条力筋线.然后分别进行单元划分,施加荷载和边界条件,最后求解.此法基于几何模型的处理,优点是力筋位置比较准确,求解结果比较准确,特别是对体内普通钢筋和体内有粘结预应力的钢筋;用这种方法建模后,钢筋和混凝土能够直接共同工作,不需进行耦合等操作.缺点是当力筋形状比较复杂时工作量很大.另外一种方法是独立建模耦合法,把实体和力筋分别独立建模,划分单元,然后将力筋单元和实体单元耦合（CP或CPIN TF.此法基于有限元模型的处理,优点是建模比较简单,用来分析力筋形状复杂且数量较多的结构比较方便.缺点是当实体单元划分不够密时,力筋节点位置可能会发生一些偏移,但误差不是很大.施加荷载阶段,可以把荷载施加在

6、实体模型上或有限元模型上,但求解之前要保证所有的荷载施加在有限元模型上.荷载可以分为DOF约束、力、表面荷载、体积荷载、惯性力以及耦合场.非线性静态分析的结果主要有位移、应力、应变和反作用力等,用通用后处理器（PO ST1或时间历程后处理器（POS T26来考察结果.用POST1一次仅可以读取一个子步,且这个子步的结果应已被写入Jobnam e.rst.通过载荷步选项命令O U TRES可以控制哪一个子步的结果被存入Jo bname.rst.可以依据载荷步和子步号或者时间来读取需要的载荷步和子步结果,但不能依据时间识别出弧长结果.如果指定了一个没有结果可用的Time值,AN SY S程序将通过

7、线性内插的方法来计算出此Time处的结果,但在非线性分析中这种线性内插会产生一定误差.因此,应在一个精确地对应于所要求子步的Tim e处进行后处理.同样也可以使用POS T26考察非线性结构的载荷-历程响应.可用PO ST26比较一个AN-SYS变量对另一个变量的关系.3非线性分析中的一些问题利用AN SYS分析钢筋混凝土结构时,当混凝土没有开裂时程序比较容易收敛,但是伴随混凝土开裂和压碎,程序的收敛就越来越困难.3.1预应力混凝土结构分析方法a.等效荷载法.将预应力筋的作用以荷载的形式作用于混凝土结构,优点是建模简单,不必考虑预应力钢筋的具体位置,网格划分比较简单,程序收敛也比较容易,对结构

8、在预应力作用下的整体效应容易求得,在按杆系结构分析时应用较多.缺点在于对预应力混凝土构件特别是体内预应力构件,预应力钢筋对混凝土的作用各处不同,等效荷载法没有考虑这种分布和方向的不同;当预应力钢筋的布置比较复杂时,用等效荷载法模拟比较困难,而且可能产生较大误差;在外荷载作用下,难以考虑外荷载和预应力钢筋的共同作用,不能模拟预应力钢筋在外荷载作用下的应力增量;对体内有粘结预应力构件,无法模拟应力损失引起的钢筋各处应力不等.b.实体力筋法.实体力筋法是同时考虑预应力钢筋和混凝土的作用,用Solid单元模拟混凝土,用Link单元模拟预应力钢筋.优点为预应力的模拟可以采用初应变法或降温法,降温法就是给

9、体外筋单元施加温度值来模拟张拉预应力的大小,温度值通过张拉预应力的值反算得到.其中降温法比较简单,同时还可以模拟预应力筋的应力损失及其在外荷载作用下的应力增量;由于预应力钢筋的位置一定,可以比较真实的模拟预应力钢筋对结构的影响和作用;作用在梁上的任何荷载都是由力筋和混凝土共同承担,可以得到预应力钢筋在任何外荷载作用下的应力和应变;可以模拟预应力钢筋的应力损失的影响.3.2混凝土各参数的设定在AN SYS分析中,混凝土材料需要输入9个参数:裂缝开裂时剪力传递系数（Shrcfof-Op、裂缝闭合时剪力传递系数（Shrcfof-Cl、单轴抗拉强度（Un Tenss T、单轴抗压强度（UnCompsT

10、、双轴抗压强度（B iCo mpsT、周围静水压力（Hydro Prs、静水应力状态下双轴极限抗压强度（BiCompsT、静水应力状态下单轴压缩极限抗压强度（UnTenss T及断裂发生时刚度乘子（Tencr Fac.其中,两个剪力传递系数一般取为0.30.5,可根据结构的情况适当调整.前4个参数必须输入,第58个可以不输入,但如果输入其中一个就要全部输入.最后一个参数程序默认值为0.6,在计算中可根据程序的收敛情况适当调整以加速收敛.当混凝土的单轴抗压强度设为-1时,不考虑混凝土的压碎.当考虑混凝土21第4期张耀庭等:AN SYS在预应力钢筋混凝土结构非线性分析中的应用压碎时,即使没有达到压

11、碎应力,程序也不容易收敛,因此建议计算时关闭压碎开关.当进行混凝土结构极限承载能力的分析时,关闭压碎开关会给结果带来一定的误差,但可避免大量的试算.3.3关于程序收敛的问题如果考虑了混凝土和钢筋的材料非线性,开裂前程序比较容易收敛,当混凝土开裂后随着荷载的增大程序收敛就变得越来越困难.a.网格密度.选择适当的网格密度能够帮助程序的收敛.网格划分太稀会导致误差较大,划分太密会花费大量时间,计算费用很高,而且对程序的收敛所起得效果也不明显.网格密度可以根据结构的实际情况及精度的要求等进行调整.通过的分析认为单元尺寸在5070mm之间比较合适,能够满足精度要求且计算费用也不是太高.b.子步数.子步数

12、（N SU BS T的设置对程序的收敛性影响较大.目前对子步数设置没有很好的方法,如果收敛曲线的走形比较长,可以考虑增大NSUBST值,计算中一般通过试算的方法来确定,根据前一次计算情况来调整NSUBST的值.当荷载较大时,一般要增大荷载子步数来减小每级加载的值,防止程序过早发散.但是计算工作量很大,有时要经过多次反复调整才能达到收敛.c.收敛准则.AN SYS结构分析中的收敛准则主要有力的收敛（F、位移的收敛（U、弯矩的收敛（M和转动的收敛（RO T.在混凝土计算中,一般选择力的收敛而不同时使用位移的收敛,否则会给收敛带来困难.可以用CNV TO L命令来调整收敛精度,加速收敛减少计算时间.

13、收敛精度默认值是0.1%,根据计算精度需要一般可以放宽到5%10%,这样收敛速度将提高.d.结构形式.结构形式和受力情况越复杂,程序越不容易收敛.通过对体外预应力混凝土梁的分析发现,对没有设转向块直线布置的梁、一个转向块折线布置的梁和两个转向块折线布置的梁,收敛越来越困难,说明结构形式对收敛性也有一定的影响.3.4非线性分析的建议通过分析认为,对开裂后的混凝土构件,可以将弹性模量和惯性矩进行折减,然后当作弹性材料来计算.惯性矩的折减可以用有效惯性矩代替毛惯性矩来模拟由于裂缝的发展对梁刚度的影响.当混凝土应力较大时,混凝土已进入弹塑性阶段,初始的弹性模量已不能正确反映此时的应力应变性质,此时应该

14、用割线模量来代替.一般认为开裂截面的割线模量Ec是其弹性模量E c的0.5倍.对整根梁而言,有的截面已经开裂,有的还没有开裂,所以弹性模量折减系数v c应大于0.5小于 1.由于混凝土截面开裂的位置很难预先确定,各处混凝土的应变也不相同,所以很难精确的确定v c的值.根据对体外预应力混凝土梁的非线性分析结果,建议在混凝土开裂后v c取0.60.7间的一个数值,但这个系数对其它钢筋混凝土结构的实用性和精度还有待更多的分析来验证.4算例用ANSYS5.7对5根体外预应力混凝土梁进行了非线性分析,并与实验结果进行对比来验证模型的精度,试验情况参见文献9.分析中混凝土单元为So lid65,非预应力钢

15、筋和钢绞线采用Link8单元.混凝土应力-应变曲线采用二次曲线+水平直线段的形式,其关系见式（1和（2,各参数采用规范中给定的值（表1.表1体外预应力梁跨中短期挠度计算值与实测值梁号跨中弯矩/kNm-1跨中挠度实测值f0/mm计算值f/mmf的误差/%Pb-147.019.229.315 1.054.3613.1714.2578.358.0415.1315.1490.261.7118.1018.0120.5Pb-236.354.29 4.98316.140.175.87 6.098 3.943.927.357.478 1.747.608.948.8580.955.1012.4313.6129.

16、562.6016.1416.838 4.3Pb-344.104.82 4.20112.851.457.547.956 5.566.1514.9613.6418.880.8523.1423.2640.6Pb-443.342.69 2.76 2.654.365.72 6.3611.261.718.847.62613.769.0611.8211.114 6.076.4114.9114.276 4.383.7617.7217.528 1.191.1120.7419.961 3.8Pb-558.044.26 4.533 6.461.115.84 5.55 5.069.067.887.729 1.976.

17、4110.8811.422 5.083.7613.1113.907 6.191.1116.0416.407 2.398.4618.5418.942 2.2105.8121.4321.4120.122华中科技大学学报（城市科学版2003年图2转向块的体外预应力梁模型示意图根据体外预应力混凝土梁的受力特点,其模型采用整体式建模法.钢筋和梁体采用分别独立建模然后耦合的方法（图2,体外预应力采用降温法来模拟.为加速收敛打开线性搜索（LN SRCH ,1和预测（PRED ,ON 等功能,并将收敛精度放宽到8%.通过POST26处理结果,得出在各级荷载作用下梁的跨中挠度值（表1.由表1可见,AN-SYS

18、计算值与实测值间比较接近,误差一般在5%左右,精度较高,说明用ANSYS 建立的模型比较符合实际情况,在计算中的一些假定是正确的.参考文献1宋天霞.有限元法理论及应用基础教程M .武汉:华中工学院出版社,1987.2刘涛,杨凤鹏.精通AN SY S M .北京:清华大学出版社,2002.3陈精一,蔡国忠.电脑辅助工程分析AN SY S 使用指南M .北京:中国铁道出版社,2001.4夸克工作室.有限元分析基础篇AN SY S 与M a tlabM .北京:清华大学出版社,2002.5王国强.使用工程数值模拟技术及其在AN SY S 上的实践M .西安:西北工业大学出版社,1999.6S K P

19、admar ajaia ,Ana nth Ramasw amy .A finite el-ement assessment of flex ura l streng th o f pr e-str essed co nc rete beams w ith fibe r reinfo rcement J .Cement &Co ncr ete Co mpo sites ,2002,（24:229-241.7程进,江见鲸.用AN SY S 软件分析薄壁杆梁截面的几何特性的计算J .交通与计算机,2001,19（增刊:77-78.8毕桂平,等.AN SY S 在桥梁工程中的应用前景A .中国土木工程

20、学会桥梁界结构工程学会第十四届年会论文集C .上海:同济大学出版社,2000.9邱继生,等.体外预应力混凝土梁短期挠度的计算J .特种结构,2003,20（1:34-37.Application of A NSYS in Nonlinear Analysis of Prestressed Reinforced Concrete StructureZ H AN G Yao -ting 1QIU J i -sheng（1.School of Civ il Eng.&Mechanic ,HUST,W uhan 430074,China;2.Schoo l o f Architecture &Civi

21、l Eng.,X US T,Xi a n 710054,China Abstract :The metho d of no nlinea r a nalysis o f prestressed reinforced co ncrete structures in ANSYS is introduced .Some pro blems in the ana lysis are discussed and som e advice and ideas are giv en .At last the nonlinear analy sis of five ex ternal prestressed

22、co ncrete beams is made in ANSYS ,and the ca lcula-tion v alues in ANSYS a re compared with the test v alues.The result indicates that the precisio n of m odels is hig h,w hich provides a new metho d a nd approach to sim ula te the mechanic capability of re-info rced co ncrete structures and has practical v alue .Key words :reinfo rced co ncrete structure ;nonlinea r ana lysis;ex ter nal prestress23AN SYS 在预应力钢筋混凝土结构非线性分析中的应用