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

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

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

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

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

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

气液二相非饱和岩体热–水–应力耦合模型及二维有限元分析

针对高放射性核废料地质处置中水及水蒸汽二相共存的非饱和孔隙岩体,根据相关的理论,推导了应力平衡方程、流体连续性方程及能量守恒方程,使用Galerkin方法将各控制方程分别在空间域和时间域进行离散,初步开发出了相应的分析热–水–应力耦合弹塑性问题的二维有限元程序。通过一个核废料处置概念库的计算认识到:高放废物地质处置后的几十年内,缓冲层中的温度及液相饱和度均上升,同时该区域内要出现明显的应力、孔隙水压力、温度、液相饱和度和气相饱和度的集中现象。     

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

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

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

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

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

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

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

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

核废料贮存裂隙岩体水热耦合迁移及其与应力的耦合分析

 博士学位论文摘要　核废料深层地质贮存的安全性分析, 一个重要指标就是要充分体现在贮库裂隙围岩介质中地下水流动、附加应力和热载耦合作用下的岩体性能稳定。为了达到这一耦合过程对岩体行为的全面理解, 试从各自耦合过程特征出发, 把裂隙岩体视为等效连续介质, 对裂隙岩体介质THM 耦合参数特性进行分析, 从而建立描述裂隙岩体介质THM 耦合的数学模型。概括起来, 从以下几方面来完成这项研究。在裂隙岩体力学模型方面, 利用O ′Connell 建立的干(或饱和水) 裂隙岩体等效弹性模量与岩石模量的关系式, 并由式中的裂隙密度概念意义, 建立了温度作用下的裂隙密度与O da 提出的裂隙张量之间的关系式。在裂隙岩体渗透模型方面, 利用O da提出的描述裂隙岩体渗透特性的附加裂隙张量, 并以裂隙结构面的开度、岩体裂隙数(包括受温度影响开通裂隙数)、裂隙连通率、附加应力、剪切膨胀和化学成份为研究对象, 建立了具有THM 耦合特性的渗流系数张量。在理论分析方面, 建立了THM三方面满足的本构方程式和描述核废料贮库裂隙岩体介质热2液2力耗散过程的定解方程。在实验基础上, 给出了温度、饱和水下的单裂隙岩体应力2应变、抗压强度回归拟合关系表达式以及岩体裂隙结构面的温度2应力2水力耦合本构关系式。在数值分析方面, 利用加权残数法理论, 导出了求解所建立的THM 耦合数学模型的有限元计算公式, 并编制了二维有限元计算程序。用BM T1 问题的算例, 获得了较满意的计算结果, 从而显示了其数值模拟的成功性。     

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

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

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

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

考虑层状岩体强度异向性的地下洞室平面有限元分析

由于层状岩体在平行和垂直层面的方向上具有不同的强度特性,故进行非线性分析时必须考虑岩体强度的异向性。据此,笔者研制了一个可以处理层状岩体强度异向性的弹塑性及无拉分析的平面有限元程序,并举例进行了计算和结果分析,得到了一些有益的认识。     

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

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

一种模拟热-水-应力耦合作用的节理单元及数值分析

对于具有相同厚度的4 节点等参数节理单元,从建立应力平衡方程、水连续性方程、能量守恒方程入手,利用Galerkin方法,将各控制方程分别在空间域和时间域进行离散,使之可用于分析岩体饱和不连续面中热-水-应力耦合弹性问题,并将这一单元模型加入到相应的二维有限元程序中。通过假定的有、无节理的稳态渗流算例和核废料地下处置算例,表现了由于不连续面的存在,使得岩体中不连续面两侧的位移、水头不再连续,而温度场虽然连续但与无不连续面时相比仍有一定差别。     

Cosserat理论在互层岩体洞室有限元分析中的应用

Cosserat理论考虑了偶应力对介质变形特性的影响,能反映出层状材料的弯曲变形。首先对基于平面应变问题的偶应力理论及扩展模型简略介绍,然后导出了模型的Mohr-Coulomb塑性屈服条件。再利用Matlab平台编制这种模型的有限元程序,并对地下洞室工程进行数值模拟验证程序的有效性。结果表明:Cosserat介质扩展模型与参数混合法互补是解决互层岩体问题的有效方法之一。