喀斯特地区地层塌陷的离散元法数值模拟-绵阳师范学院学报.pdf

第３７卷　第５期 绵阳师范学院学报 Ｖｏｌ．３７　Ｎｏ．５ 　 ２０１８年５月 ＪｏｕｒｎａｌｏｆＭｉａｎｙａｎｇＴｅａｃｈｅｒｓ＇Ｃｏｌｌｅｇｅ Ｍａｙ．，２０１８ 　 ＤＯＩ：１０．１６２７６／ｊ．ｃｎｋｉ．ｃｎ５１－１６７０／ｇ．２０１８．０５．０２１ 喀斯特地区地层塌陷的离散元法数值模拟 梁冠军，尤　超 （滁州职业技术学院土木工程系，安徽滁州　２３９００００） 摘　要：运用ＬＳ－ＤＹＮＡ非线性动力学分析软件，采用离散单元法对喀斯特地区地层塌陷现象进行数值模 拟，并进一步分析了不同塌孔直径对塌陷效果的影响结果表明：在实际案例所涵盖的粒径及塌孔尺寸范围内，数 值模拟给出了地层塌陷漏斗母线和塌陷效果具有较好的一致性；岩层开洞尺寸的大小对土粒塌陷速度和塌陷形态 有一定的影响 关键词：喀斯特地区；离散元法；塌陷；数值模拟 中图分类号：ＴＵ４３３　　　　文献标志码：Ａ　　　　文章编号：１６７２６１２Ｘ（２０１８）０５０１１２０６ ０　引言 地层塌陷是指岩土体在自然或人为因素作用下，向下陷落，并在地面形成塌陷坑洞的一种地质现象而 根据地质的发育条件，地层塌陷可分为喀斯特塌陷和非喀斯特塌陷近几年，地面塌陷事件频多，对人们的 生命财产安全造成极大的影响 ［１－６］ 国内外诸多学者 通过理论分析、试验方法以及数值模拟对地层塌陷的机理和过程进行了研究欧 ［５］ 阳振华等人 分析了土洞形成的机理及其力学机制，研究了土洞扩展的条件，分析了地下水在地表塌陷形 ［６］ 成过程中的作用，从力学的角度揭示了地表塌陷的本质李涛等人 采用理论分析和数值模拟相结合的方 法研究了地下空洞引起土质地层地陷的形成机制，并提出了基于土层厚度和地陷漏斗最大沉降的地陷预测 与评估方法 上述部分文献虽然对地层塌陷的研究也才用数值模拟，但是运用ＬＳ－ＤＹＮＡ非线性动力学分析软件， 采用离散单元法进行数值模拟的研究却很少见，因此，本文的研究工作是对分析地层塌陷成因领域的补充， 也为认识喀斯特地区地层塌陷的防治提供参考 １　地层塌陷案例描述 如图１所示的两个案例案例一：２００１年１２月，新西兰北岛的怀赫镇出现宽５０ｍ、深１５ｍ的巨坑，导 致地面上一处房屋倒塌居住在此的５位居民，包括３名儿童，都受到重伤 案例二：２０１３年２月，四川省广元市毛坝村发生地陷，形成的“天坑”直径长２４．９ｍ“天坑”在一处岩溶 地貌上形成，该处地表在２０１２年９月曾连续６天下沉 上述两个案例均是由于喀斯特地区的岩溶作用形成地下巨大空洞，致使地面出现大面积的塌陷 收稿日期：２０１７－１２－２０ 基金项目：２０１５年安徽省高度学校自然科学重点项目（ＫＪ２０１５Ａ３６２）． 第一作者简介：梁冠军（１９８３－），男，安徽涡阳人，工程师，讲师，硕士，研究方向：土木工程数值仿真 ·１１２· 　 绵阳师范学院学报（自然科学版） 图１　喀斯特地区地层塌陷案例 Ｆｉｇ．１　ＣｏｌｌａｐｓｅｃａｓｅｉｎＫａｒｓｔａｒｅａ ２　计算模型 ２．１　数值模型 采用如图２所示的喀斯特地区地层模型，地层总厚度为１０ｍ，在目标土体外围的岩层用有限元建模，土 粒材料的本构关系用离散元粒子来描述，土体颗粒与岩层间接触采用合适的耦合算法，考虑到计算机的运算 能力及计算耗时，土体颗粒粒径取为０．１～１ｃｍ为了区分地层之间的位移变化，把地层分割为１０层，每层 厚度１ｍ，用不同的颜色标出在模型底部有限元岩层开洞尺寸为１ｍ×４ｍ 图２　喀斯特地区地层模型 Ｆｉｇ．２　ＭｏｄｅｌｏｆＫａｒｓｔａｒｅａ