大型地下洞室三维弹塑性损伤动力有限元分析.pdfVIP

第 29 卷 第 5 期 岩石力学与工程学报 Vol.29 No.5 2010 年 5 月 Chinese Journal of Rock Mechanics and Engineering May ，2010 大型地下洞室三维弹塑性损伤动力有限元分析 张志国，肖 明，张雨霆，左双英 (武汉大学 水资源与水电工程科学国家重点实验室，湖北 武汉 430072) 摘要：采用动力时程分析法研究地下洞室地震响应特性。根据运动学动力平衡方程，以 Newmark 直接积分法为基 础，结合三维弹塑性损伤有限元理论，编写动力时程计算程序。采用塑性附加荷载、损伤附加荷载思想考虑单元 塑性及损伤情况，每时步计算对单元损伤状态进行继承，反映损伤不可逆性。采用增量变塑性刚度迭代法，求解 非线性动力方程。计算结果表明，该方法一般 3 步内收敛，计算快，精度高。阐述大型地下洞室动力响应时程分 析中，模型边界条件的设置方法及动荷载的输入方式。通过对渔子溪大型水电站地下厂房洞室群的模拟计算表明， 该程序计算速度较其他商业程序速度快，计算结果合理，准确，可以应用于地下工程动力时程分析。 关键词：岩石力学；地下洞室；动力时程分析；弹塑性；边界条件；荷载输入 中图分类号：TU 45 文献标识码：A 文章编号：1000 –6915(2010)05 –0982 –08 DYNAMIC FINITE ELEMENT ANALYSIS OF LARGE-SCALE COMPLEX UNDERGROUND CAVERNS WITH THREE-DIMENSIONAL ELASTOPLASTIC DAMAGE MODEL ZHANG Zhiguo ，XIAO Ming ，ZHANG Yuting ，ZUO Shuangying (State Key Laboratory of Water Resources and Hydropower Engineering Science ，Wuhan University，Wuhan，Hubei 430072 ，China) Abstract ：The seismic response characteristics of underground caverns are studied by employing dynamic time-history method. Based on the dynamic equilibrium equations of kinemato-dynamic time-history method and Newmark direct integration method ，the computational codes are programmed by adopting three-dimensional elastoplastic finite element theory. The concepts of additional plastic load and additional damage load are considered to simulate the plasticity and damage of elements. The damage state of rock mass within each time step is inherited by subsequent time steps. Therefore，the irreversible trend of rock damage can be reflected. The incremental variable plastic stiffness iteration method is emplo