基于离散连续耦合方法的地下结构在地震中破坏过程模拟.pdfVIP

第34卷第1期 固体力学学报 V01．34No．1 CHINESE OFSOI。IDMECHANICS JOURNAL 2013 2013年2月 February 基于离散一连续耦合方法的地下结构 在地震中破坏过程模拟+ 金炜枫1’2“ 周 健2。 张姣4 (1浙江科技学院建筑工程学院，杭州，310023)(2同济大学地下建筑与工程系，上海，200092) (。同济大学岩土及地下工程教育部重点实验室，上海，200092)(4上海城市管理职业技术学院土木工程与交通学院，上海，200043) 摘 要基于改进后的土体离散一连续耦合算法，对坂神地震中神户大开地铁车站的坍塌破坏过程进行了数 值模拟．通过改进适用于Lagrange乘子法求解的优化模型，来强化离散和连续模型在耦合边界上力的相容性，并 将改进后的算法嵌入离散连续耦合动力模型．在耦合边界通过交换速度和力保证耦合模型的连续性，通过自振柱 模拟和循环双轴模拟使离散区域土体的宏观动力属性与连续土体模型一致．土体离散连续耦合方法可以从细观 尺度描述重点考察区域土体和结构的破坏过程，并有效降低离散元的模拟规模．模拟得到的车站破坏形式与文献 中现场情况非常接近，并从定量角度与文献中基于连续体的模拟结果进行了比较． 关键词 双尺度耦合模型，离散元，细观土力学 0 引言 实质上应该是个动力模型，一开始连续模型中的土 阪神地震中，神户大开地铁车站遭受破坏，中间 体是取静力属性的，如Cai等凹1用离散元模拟隧道 立柱顶部或底部的混凝土被压碎，这些部位的钢筋 附近土体，用连续模型模拟其它范围土体，Jin和 鼓出扭曲，顶板在中间立柱上方附近断裂下塌，造成 Zhou[71用离散元模拟靠近桩刺人的土体，其它区域 车站上方土体下陷[1。3]．这一现象为典型的地震中地 用连续模型模拟．周健等[83给出离散元的自振柱模 下结构破坏过程，震后学者们对大开地铁车站的地 拟方法，完善了联系颗粒集合体动力宏观属性和细 震响应进行了一系列的数值模拟．Huo等_·用弹塑观属性的方法，并基于耦合模型模拟离心机试验中 性模型模拟土体，曹炳政等H3模拟时考虑到土体剪 挡土墙的地震响应．金炜枫等凹]给出强化边界耦合 切模量随剪应变衰减的特性并用弹性模型模拟地下 力的优化模型，用嵌入此模型的耦合方法模拟列车 结构，庄海洋等口]采用土体非线性粘弹性动力本构 振动下隧道和土体的响应．耦合问题的核心在于[8]： 模型，并用黏塑性损伤模型模拟混凝土．以上的模拟 (1)离散颗粒集合体的宏观属性与连续模型一致， 工作都可以确定地震中地下车站的破坏位置或者地 这需要合适的模拟方法将颗粒集合体的宏细观属性 下结构的损伤演化过程，但都没有模拟实际地下结 联系起来，已有的模拟方法有获得静模量和塑性性 构破坏断裂下塌和土体下陷大变形的过程． 质的双轴压缩数值试验口…，获得颗粒集合体小应变 离散颗粒的动力模拟可以从细观尺度上模拟土 下动模量的自振柱模拟方法[8]，和获得大应变时动 体的大变形和破坏等力学行为．但由于计算机能力 模量的循环双轴试验[1lj；(2)离散一连续模型在耦合 的限制，模拟的颗粒数目是有限的，对于半无限土体 边界上保持连续性，目前的方法为保证速度和力的 中地下结构的地震响应，离散元模拟难以实现．而土 连续．对于提取边界耦合力有不同的算法，只在离散 体离散一连续耦合动力模型在重点考察的区域用离 模型中提取时有基于颗粒作用于耦合边界位置的提 散颗粒模拟，而其他区域用连