热一水一应力耦合弹塑性二维有限元程序的研制.pdfVIP

2005年 1月 焦作大学学报 No．1 第 1期 JOURNAL 0F JIAOZUO UNIVERSITY Jan．2005 热一水一应力耦合弹塑性二维有限元程序的研制 张玉军 (中国科学院武汉岩土力学研究所 ，湖北 武汉 430071) 摘要 ：从建立应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手，使用Galerkin方法， 将各控制方程分别在空间域和时间域进行离散，初步开发 出了一个用于分析饱和岩土介质 中 热一水一应力耦 台弹塑性问题的二维有限元程序。通过一维固结问题的有限元解与太沙基解析 解的对比，以及对一个简单的热一水一应力耦合问题的数值计算，验证 了本程序的正确性。 关键词 ：热一水一应力耦合 ；弹塑性分析 ；二维有限元程序 中图分类号 ：TU45 文献标识码：A 文章编号：1008—7257(2005)01—0069—05 十多年来 ，由于高放射性核废料地下处置需要 的 do=do一m (1) 推动 ，欧美及 日本等发达 国家对岩土体 中的流体一温 如 =DT(dF—d8 一dsp—d8r—d8．) (2) 度一力学的耦合现 象作 了大量 的理论 和试 验工作。 对于一处置场地而言 ，岩体 中的裂隙是地下水 和可能 而d8一 妾，d8r-m争 ，如-c 泄漏的放射性核素运动的主要通道 。在一定的热力、 式 中，D，——切线弹塑性模量矩阵； 水力梯度及地应力场 的作用下 ，地下水及部分泄漏 的 如—— 总应变增量 ； 核素会在一定范 围的岩体迁移 ，同时产生着物理、化 dsp——孔 隙水压力 引起 的应变增量 ： 学反应 (如相变 ，溶蚀 、析 出和再结 晶)。因此 ，为 了 如 ，——温度引起 的应变增量 ； 保证核废料处置 的绝对安全 ，预测岩体 中伴随着地下 如 ——蠕变引起 的应变增量 ； 水流动 的可能的核 素扩散是特别重要 的。在这一点 如。——其他原因引起 的应变增量 ； 上 ，有 限元方法等数值模拟是非常有效的分析工具 。 m——法 向应力单位列阵； 为此 ，欧美及 日本等 国纷纷开发 出了不少著名 的相关 卢—— 固相 的热膨胀 系数 ； 计 算 程 序… ，如 MOTIF 、THAMES 、FRACON 、 ， — — 固相的体积模量 ； ABAQUS 、UDEC 等 。鉴于 目前 国内高放射性核废 — — 孔 隙水压力增量 ； 料地下处置 已经开始进 行但还缺少前述类计算软件 d——温度增量 ； 的状况 。作者有 意在这方面做一些开发性 的工作 。故 c——蠕变参数 。 利用在 日本京都大学工学部从 事访 问学 习的有利机 根据虚功原理 ，增量 型的应力平衡方程可写为 会 ，初步开发 出了一个用于分析饱和孔 隙介质 中热一 f88dodO—f8udbOd 一』8udidF=0 (3) 水一应力耦合弹塑性 问题 的二维有限元程序 ，并对其 力 力 r ’ 式 中，d——总应力增量 ； 正确性进行 了试算验证 。现将 主要 的程序开发工作 、 — — 分别为体积力和面力增量 ； 介绍如下 。