地震作用下砌体结构倒塌反应的数值模拟计算分析.docVIP

　第38卷第9期　2005年9月 土　木　工　程　学　报 CHINACIVILENGINEERINGJOURNAL Vol138　No19Sep1　2005 地震作用下砌体结构倒塌反应的数值模拟计算分析 苗吉军 (青岛理工大学,山东青岛266033) 3 顾祥林　张伟平　彭　斌 (同济大学,上海200092) 摘要:本文建立了能反映砌体结构出平面破坏的三维空间刚体离散单元模型,并编制了计算机仿真程序,模拟了地震作用下砌体结构的倒塌反应,模拟结果和振动台试验结果比较,效果良好。关键词:离散单元法;砌体结构;倒塌分析;数值模拟中图分类号:TU35211+2　　　　　文献标识码:A文章编号(2005)0920045208 NUMERICALSIMULATIONANALYSISJiJun GuXiangLin　ZhangWeiPing　PengBin (Qingdao,Qingdao266033,China)(TongjiUniversity,Shanghai200092,China) Abstract:Inthis,thediscreteelementmethod(DEM),isemployedforthefirsttimetoperformmasonrystructuralanalysisunderdynamicloading.A32dimentinalmodelisestablishedtosimulatethefailure.Thecomputerprogramdevelopedcansimulatethecollapseresponsesofthemasonrystructuresunderearthquakes.Thesimulationresultsareingoodagreementwiththeshakingtabletestdataforamasonrystructuralmodel. Keywords:discreteelementmethod;masonrystructure;collapseresponse;numericalsimulation 前　　言 砌体结构在地震作用下最容易遭受破坏甚至倒塌,但是关于砌体结构倒塌反应方面的研究却很少。目前对于砌体结构分析的研究大多数基于有限单元 ] 法,如Hillerbourg[1、Bocca[2]、Pande[3]都曾成功地模拟了砌体结构中裂缝的发展情况。朱伯龙[4]曾经利用接触面单元来模拟砂浆的作用,并对“混用体系”和粉煤灰砌块两个墙片的拟静力试验进行了有限元分析。发现,前者与试验的吻合程度远好于后者。究其原因在于前者由于混凝土芯柱的约束作用较强,使有限元所必须的位移协调条件得到较好的满足。但是在下降段,前者和试验依然有一定的误差。这说明当砌体结构进入下降段以后,随着裂缝的开展和变形的增大,应用有限单元法有着明显的不足。 离散单元法是P1A1Cundall[5]于1971年提出用于分析岩石节理系统或块体集合等不连续体系的大变形运动的数值方法。随着研究的深入和计算机技术的发展,离散单元法除了用于边坡、采矿和巷道的稳定性研究以及颗粒介质微观结构的分析外,已扩展到用于 收稿日期:2003207204 教育部留学回国人员科研启动基金资助项目 研究地震、爆炸等动力过程和地下水渗透、热传导等物理过程。尤其在结构低应力水平的平面应力和应变问题上逐步得到了完善,并推出了相应的计算机仿真软件。因此,本文尝试将离散单元法引入到地震作用下砌体结构的倒塌反应分析中来。 1　结构的离散化 当块体强度较高,砂浆强度相对较低时,砌体结构中的灰缝实际上是“薄弱层”。由于灰缝砂浆“薄弱层”的存在,砌体结构和Cundall最初提出的离散单元法的分析研究对象———有“节理”岩石———是非常类似的。因此,本文把砌块作为一个无限刚性的单元,灰缝砂浆离散成连接弹簧,见图1(在后续研究中将考虑构造柱和圈梁的影响,本文暂不考虑)。这和朱伯龙等用接触面单元来模拟砂浆作用的思路基本一致[4]。从材料层次上分析,连接单元的“组合弹簧”可以看成是图1e所示的抗拉弹簧和抗压弹簧(用Knz来表示),以及在垂直抗拉和抗压弹簧平面内的两个剪切弹簧(用Ksx和Ksy来表示)。为了保证计算的稳定性,并反映弹簧单元变形和砌块单元运动时的耗能特性,在图1e中的每个弹簧处,并联一个相应的阻尼器。刚体单元之间的连接关系示于图2。 　　　　　　　　　　　　　　　　　　　　　土　木　工　程　学　报?46?2005年 图1　砌体结构房屋的离散化过程 Fig11　Thediscretizationprocessofthemasonry structures 图2　典型块体单元之间