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

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

轴对称荷载下饱和地基热-水-力耦合动力响应

基于Biot波动理论,通过建立饱和地基的热-水-力耦合动力响应的控制方程,对轴对称荷载作用下半无限地基的热-水-力耦合动力响应问题进行探讨。利用Hankel变换技术,得到外荷载作用下地基中温度增量、应力、位移和孔隙水压力积分形式的解答。利用Hankel数值逆变换得到计算结果,对热-水-力耦合条件下地基土体中温度增量、应力、位移和孔隙水压力响应的分布进行分析,并与水-力耦合动力响应情况下的结果进行比较。     

尤卡山坑道规模试验T-H-M耦合过程数值模拟

简述了笔者建立的求解热-水-应力耦合问题的数学模型,针对美国尤卡山高放射性核废料地质处置坑道规模试验使用笔者研制的有限元程序进行了热-气-应力耦合过程的数值模拟,分析了坑道周边的位移、温度及气压力的变化;并将该计算的部分结果和实测数据及国外有限元/有限差分方法的成果作了比较,看到三者温度的分布与变化较为接近,位移分布与LBNL的结果在定性和定量上有某种一致性,位移随时间的变化与实测结果的规律性大体相当.从而可以认为使用该模型和有限元程序对核废料地质处置中的热-水/气-应力耦合过程进行模拟是可行的.     

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

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

考虑变形及滞回效应影响的三维土-水特征曲面模型

为研究变形和滞回效应对非饱和土水-力耦合特性的影响,以边界面理论为基础,建立一个同时考虑变形及滞回效应影响的三维非饱和土土-水特征曲面模型,所建立模型可完整地描述非饱和土的水-力耦合特性。提出的本构方程以吸力和孔隙比为自变量,以饱和度为因变量,建立三维土-水特征曲面模型,通过编程实现了本模型的预测功能。通过程序的预测结果与一系列不同应力和水力路径下的试验结果进行对比,在饱和度-吸力/孔隙比二维平面及饱和度-吸力-孔隙比三维空间中,均验证了所建立模型的适用性及预测精度。     

高放废物工程屏障的热-水-力耦合作用分析

采用多相多类方法，提出热-水-力耦合作用分析的综合模型，分析考虑多孔介质的三相（固，液，气）和三类（矿物，水分，空气）模型，该模型是根据控制方程，本构关系和平衡条件来建立的，并提出用于模拟高放废物地下处置体系和附近场地性状的热-水-力耦合分析通过程序。热-水-力耦合分析可将工程屏障中不同过程间发生的主要相互作用，以及附近岩石的散热和水化现象都考虑进去，通过试验结果的验证，得出一些关于热-水-力方面的重要结论。     

高放废物地质处置库中非饱和缓冲层的热-水-力耦合数值模拟

作为高放废物地质处置库中的缓冲/回填材料,膨润土在核废料放出的热量和周围天然岩体中地下水共同作用下,其热-水-力耦合作用效果对处置库的稳定性、安全性有着极为重要的影响。以非饱和土多场耦合理论为基础,基于ABAQUS有限元分析平台,建立二维轴对称模型,对高放废物地质处置库中膨润土缓冲层的热-水-力耦合过程开展了有限元模拟,给出了处置库近场缓冲层中不同位置的温度、饱和度、竖向应力及位移、孔隙水压力及吸力的分布及变化规律,并重点分析了温度场的影响,通过与已有研究成果的对比,证实了所计算结果的合理性。     

核废料地质贮存介质黏土岩的三维各向异性热-水-力耦合数值模拟

为了解黏土岩在放射性废料长期贮存中的热-水-力耦合过程,结合MontTerri核废料贮存地下岩石试验工程中黏土岩各种物理量的各向异性特点,应用多孔介质力学耦合理论研究了该黏土岩在加热和冷却全过程中由于热荷载引起的耦合效应场。研究过程考虑温度升高引起的孔隙水黏滞性改变对渗透系数的影响。研究结果表明,岩体力学参数、水力学参数和热传导参数的各向异性特性是影响岩体的温度场、孔隙压力场和应力场分布的最主要因素。各向异性耦合模型与各向同性耦合模型的数值模拟对比研究结果表明：各向异性模型数值结果能更加客观地反映该地下岩石试验工程中黏土岩在受热状态下的热-水-力耦合效应;同时,也表明岩体在加热过程中一直处于受压状态,而在冷却过程中局部会出现拉应力,从而有可能导致拉裂缝的产生。     

CO_2地层注入的热流力耦合模型及有限元分析

为考察CO2地层注入过程中岩体中发生的热—流—固耦合现象,笔者借用R.Leiws等建立的可变形孔隙介质中非等温空气流和水流模型,但在其水连续性方程中加入了温度梯度引起的水分扩散项,研制出相应的热—水—气—应力耦合二维弹塑性有限元程序。针对一个假定的CO2深部地质埋置库模型,在一定的CO2注入速率和注入时间条件下进行了热—水—气—应力耦合数值模拟,考察了岩体中的位移、应力、气体流速、孔隙压力和温度的分布与变化。结果显示,CO2的注入对围岩中的位移和应力影响强烈。     

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

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