地震荷载作用下地下厂房围岩稳定的离散元计算方法-长江科学院院报.pdfVIP

　 第２６卷第１２期 长　江　科　学　院　院　报 Ｖｏｌ．２６　Ｎｏ．１２　　 　　２００９年１２月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ｄｅｃ．２００９　 　 　　文章编号：１００１－５４８５（２００９）１２－００５８－０５ 地震荷载作用下地下厂房围岩稳定的 离散元计算方法研究 １，２ ３ １ １ １ １ １ 王　涛 ，熊　将 ，郭武祥 ，李宏明 ，胡少华 ，李　杨 ，胡　楠 （１．武汉大学 水资源与水电工程国家重点实验室，武汉　４３００７２；２．长江科学院水利部岩土力学与工程重点实验室， 武汉　４３００１０；３．北京国电水利电力工程有限公司，北京　１０００２４） 摘要：对离散元动力计算方法的理论基础进行了探索，主要是离散元动力计算的实现过程，边界条件和动力荷载的 输入问题。在此基础上以某大型水电工程的地下厂房为例，计算在静力和地震荷载作用下，围岩的变形破坏情况， 并监测了岩体不同时间速度和位移的变化情况。计算方法与成果可以为地下厂房的抗震工程设计提供指导。 关　键　词：地震；离散单元法；地下厂房；围岩 中图分类号：Ｏ２４２．１；ＴＵ４５　　　文献标识码：Ａ １　概　述 ２　动力离散元计算原理 岩体是一种非线性非连续的复杂介质，其主要 离散元中块体的运动采用牛顿第二定律来表 ［４］ 特点之一是分布着大量随机的各种结构面。如果工 示 ，即  程岩体稳定问题与岩体结构面密切相关，采用基于 ｄ　ｕ Ｆ ＝ 　， （１） ｄｔ ｍ 连续介质理论的计算方法不能得到合理的解答，采 ［１］ 进行中心差分，则 用非连续节理岩体计算方法是必需的 。    ｄ　ｕ 　ｕ（ｔ＋ ｔ／２）－　ｕ（ｔ－ ｔ／２） Δ Δ 离散单元法（ＤＥＭ）是专门用来解决不连续介 ＝ 　， （２） ｄｔ Δｔ 质问题的数值模拟方法，该法把节理岩体看成是由 速度方程变为 离散的岩块和岩块间的节理面所组成。岩块能移 Ｆ（ｔ）   Σｉ 　ｕ（ｔ＋ ｔ／２）＝　ｕ（ｔ－ ｔ／２）＋（ ＋ｇ）ｔ， 动、转动和变形，而节理面可被压缩、分离或滑动。 ｉ Δ ｉ Δ ｉ Δ