模拟热–水–应力耦合作用的三维节理单元及其数值分析

对于具有相同厚度的空间8节点等参数节理单元,考虑温度和孔隙水压力的影响,将其的应力–位移关系加以拓展,建立了相应的应力平衡方程、水连续性方程、能量守恒方程,使之可用于分析饱和–非饱和岩体中存在不连续面时的热–水–应力耦合问题,并结合实体单元,研制出相应的三维有限元程序。针对一个假定的核废物地质处置模型,通过有、无节理的热–水–应力耦合过程数值分析,显示与无节理的情况相比,由于节理的存在使得近场岩体中的温度、应力及水流速度的分布与数值都发生了明显的变化,特别是节理附近的岩体中有某种水流"冷却"作用和一定程度的应力集中现象,其它部位岩体的应力也有所增加。     

热-水-应力耦合弹塑性二维有限元程序研制

通过从已有的应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,使用Galerkin方法。将各控制方程分别在空间域和时间域进行离散,初步开发出了一个用于分析饱和岩土介质中热-水-应力耦合弹塑性问题的二维有限元程序。通过一维固结问题的有限元解与太沙基解析解的对比,以及对一个简化的核废料地下处置的热-水-应力耦合问题的数值计算,验证了本程序的正确性。     

裂隙岩体的热-水-应力耦合模型及二维有限元分析

尝试将Barton—Bandis模型和Oda的裂隙张量理论结合起来应用于饱和-非饱和裂隙岩体中热-水-应力耦合过程的计算分析,研制了相应的二维有限元程序,并列举算例比较了无、有多组裂隙岩体中的应力场、渗流场和温度场的分布及变化情况。由此看到:由于不连续面的存在减弱了母岩的刚度和增大了其透水能力,使得裂隙岩体中的应力集中程度降低和渗流速度提高,从而热源产生的热量可较快地被水流传输到周围区域中去,岩体中的温度和负孔隙水压力量值也相对较低。     

遍有节理岩体的双重孔隙-裂隙介质热-水-应力耦合模型及有限元分析

建立一种饱和-非饱和遍有节理岩体的双重孔隙-裂隙介质热-水-应力耦合模型,其特点是应力场和温度场是单一的,但具有不同的孔隙渗流场和裂隙渗流场,以及可考虑裂隙的组数、间距、方向、连通率和刚度对本构关系的影响,并研制出相应的二维有限元程序.针对一个假定的高放废物地质处置库,就岩体为双重介质和单重介质2种工况进行数值分析,考察缓冲层和岩体中的温度、孔隙水压力、饱和度、地下水流速和主应力的变化、分布情况.结果显示,地下水由双重介质进入缓冲层中要快得多,2种工况的计算域中温度差别不大,但缓冲层及附近部位的主应力大小及分布有显著不同,单重介质的应力集中程度要大.     

核废料地质处置近场热-水-应力-迁移耦合二维有限元分析

引入渗透迁移方程,将所开发的饱和-非饱和孔隙介质中热-水-应力耦合弹塑性模型及其二维有限元程序进行了扩展和改进,使之可以同步地对温度场、渗流场、应力场和放射性核素浓度场的变化进行解析.然后针对一个假想的核废料处置库算例,考察了在热-水-应力-迁移耦合作用条件下近场温度、孔隙水压力及放射性核素浓度的分布及变化.     

考虑水-应力耦合作用的地下洞室的锚杆支护效果

从建立应力平衡方程、水连续性方程入手,使用Galerkin方法,将各控制方程分别在空间域和时间域进行离散,开发出了一个用于分析饱和岩土介质中水–应力耦合弹性问题的二维有限元程序。通过处于地下水位以下的洞室施工时锚杆支护效果问题的模拟计算表明:在有地下水赋存的条件下,锚杆的作用仍主要是约束围岩位移,但比无水时有所减弱,其对岩体中渗流的抑制效果不大,而渗透系数的取值方式将明显地影响孔隙水压力的分布。     

核废料处置概念库近场热-水-应力耦合二维有限元模拟

建立了分析饱和–非饱和介质中热–水–应力耦合现象的应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵,并使用所开发的有限元程序对一个核废料处置概念库近场的热–水–应力耦合过程进行了数值模拟,考察了缓冲层及岩体中若干点的温度、饱和度、孔隙水压力及水平方向正应力随时间的变化情况,并将部分结果与国外类似软件的计算数据作了对比,看到二者的规律在定性上有某种一致性,但定量上仍存在一定差别,笔者试着分析了原因并得出了一定的认识。     

饱和-非饱和介质水-应力耦合弹塑性二维有限元分析

从建立应力平衡方程、水连续性方程和弹塑性矩阵入手，使用Galerkin方法，将各控制方程分别在空间域和时间域进行离散，初步开发出了一个用于分析饱和-非饱和孔隙介质中水-应力耦合弹塑性问题的二维有限元程序。通过引入一组特定的与非饱和状态有关的计算公式，对一个假定的非饱和土体中水-应力耦合问题进行了数值计算，考察了不同时间土体中的位移、孔隙水压力、有效主应力、流速和塑性区的分布，定性验证了该程序的正确性。     

不同场热-水-应力耦合过程二维有限元分析

从建立应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,提出了一个饱和-非饱和孔隙介质中的热-水-应力耦合模型和开发了相应的有限元程序.为验证本耦合模型及程序的可靠性和重点探讨不同场的热-水-应力耦合过程,以高放射性核废料地质处置的Prototype Repository Project(PRP)原位试验为模拟对象,对一个处置试验坑道近场进行了水单场、热-水耦合、水-应力耦合、热-应力耦合和热-水-应力耦合条件下的数值分析,考察了工程屏障及围岩中的温度、饱和度及应力的变化、分布情况,并得出了一定的认识.     

遍有节理岩体的双重孔隙–裂隙介质热–水–应力耦合模型及有限元分析

 建立一种饱和–非饱和遍有节理岩体的双重孔隙–裂隙介质热–水–应力耦合模型,其特点是应力场和温度场是单一的,但具有不同的孔隙渗流场和裂隙渗流场,以及可考虑裂隙的组数、间距、方向、连通率和刚度对本构关系的影响,并研制出相应的二维有限元程序。针对一个假定的高放废物地质处置库,就岩体为双重介质和单重介质2种工况进行数值分析,考察缓冲层和岩体中的温度、孔隙水压力、饱和度、地下水流速和主应力的变化、分布情况。结果显示,地下水由双重介质进入缓冲层中要快得多,2种工况的计算域中温度差别不大,但缓冲层及附近部位的主应力大小及分布有显著不同,单重介质的应力集中程度要大。     

孔隙率对岩体中THM耦合作用的弹塑性有限元分析

笔者使用Taron等提出的颗粒聚集体的压力溶解模型和摩尔-库伦准则,假设在饱和的石英颗粒聚集岩体中有一实验室尺度的高放废物地质处置库模型,针对其拟定两种计算工况:(1)内摩擦角和黏聚力均为常数;(2)内摩擦角为常数,黏聚力是孔隙率的负指数函数,进行4a处置时段的热-水-应力(THM)耦合有限元数值模拟,就岩体中的温度场、渗流场、应力场和浓度场的变化及分布情况进行了考察。结果显示:与黏聚力是常数的情况相比,黏聚力随孔隙率而变化时,岩体相同部位进入屈服阶段的时机滞后,塑性区减小,并推迟了塑性部位的溶质浓度、迁移/沉淀质量、反应体积和颗粒贯穿深度的突变时间;弹塑性分析中由于应力调整和增大了分子扩散系数,使得塑性区的颗粒介质的溶解、迁移和沉淀相比于弹性区有明显的变化,并对渗流场(孔隙水的压力及流速)和应力场产生显著的影响。     

饱和–非饱和介质水–应力耦合弹塑性二维有限元分析

从建立应力平衡方程、水连续性方程和弹塑性矩阵入手,使用Galerkin方法,将各控制方程分别在空间域和时间域进行离散,初步开发出了一个用于分析饱和–非饱和孔隙介质中水–应力耦合弹塑性问题的二维有限元程序。通过引入一组特定的与非饱和状态有关的计算公式,对一个假定的非饱和土体中水–应力耦合问题进行了数值计算,考察了不同时间土体中的位移、孔隙水压力、有效主应力、流速和塑性区的分布,定性验证了该程序的正确性。     

热–水–应力耦合模型及FEBEX原位试验二维有限元分析

从建立应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,提出了一个饱和–非饱和孔隙介质中且可考虑膨胀应力的热–水–应力耦合模型,并编制了相应的二维有限元程序。以高放射性核废料地质处置的FEBEX原位试验为模拟对象进行了数值分析,比较了实测值和计算值,从而验证了程序的正确性并从分析结果中得到了一些有益的认识。