不连续面对饱和–非饱和介质热–水–应力耦合影响的二维有限元分析.PDF

第25 卷 第 12 期 岩石力学与工程学报 Vol.25 No.12 2006 年 12 月 Chinese Journal of Rock Mechanics and Engineering Dec .，2006 不连续面对饱和–非饱和介质热–水–应力耦合 影响的二维有限元分析 张玉军 ( 中国科学院 武汉岩土力学研究所，湖北 武汉 430071) 摘要：从建立应力平衡方程、水连续性方程、能量守恒方程入手，利用 Galerkin 方法，对于岩体和不连续面分别 采用等参数单元和具有相同厚度的四节点节理单元，将各控制方程分别在空间域和时间域进行离散，建立了用于 分析存在不连续面的饱和–非饱和介质中的热–水–应力耦合弹塑性问题的有限元模型，并开发了相应的计算程 序。通过假定的有、无节理的核废料地下处置算例，表现了由于存在渗透系数很大的不连续面，推迟了缓冲材料 达到饱和的时间，并对岩体中的位移场、应力场和渗流场产生很大的影响。 关键词：岩石力学；饱和–非饱和介质；热–水–应力耦合；不连续面；二维有限元分析 中图分类号：TU 45 ；O 242 文献标识码：A 文章编号：1000 –6915(2006)12 –2579 –05 2D FE ANALYSIS FOR INFLUENCE OF DISCONTINUITIES ON COUPLED THERMO-HYDRO-MECHANICAL PROCESSES IN SATURATED-UNSATURATED MEDIA ZHANG Yujun (Institute of Rock and Soil Mechanics ，Chinese Academy of Sciences，Wuhan，Hubei 430071 ，China) Abstract ：In this paper ，firstly the momentum equilibrium equation，the continuity equation for water ，the energy balance equation and the elastoplastic matrix are established ；secondly by using the Galerkin method and isoparametric elements and 4-node elements with equal thickness for porous media and discontinuities respectively ，all the governing equations are discretized in the space domain and in the time domain respectively ； then a two dimensional FEM code for analysis of coupled thermo-hydro-mechanical elastoplastic processes in saturated-unsaturated porous media with discontinuities is developed tentatively. Through the numerical computation for an assumed coupled thermo-hydro-mechanical phenomena in nuclear waste repository considering the case without discontinuity and the case with discontinu