非饱和土的非正交弹塑性本构模型

从应变增量大小和方向的角度出发,以有效饱和度和Bishop有效应力作为基本本构变量,呈现了一个非饱和土的非正交弹塑性本构模型建模框架.建模框架采用非正交的塑性流动法则,直接通过屈服函数确定塑性应变增量的方向,避免了构造塑性势函数的困难.基于所提建模框架,利用屈服函数、硬化参数和土水特征曲线的具体表达式,发展得出一个非饱和土的非正交弹塑性本构模型.所提模型通过材料参数μ可以描述具有不同刚度的非饱和土的力学与水力行为.模型所需材料参数均具有明确的物理意义,且可以通过室内试验确定.利用两种经典应力路径(即恒基质吸力的常规三轴剪切路径和恒净应力的吸湿路径)分析了具有不同刚度的非饱和土的力学与水力行为.模型预测结果与非饱和土试验数据的对比验证了所提模型的优越性.     

三维双重孔隙-裂隙介质热-水-应力-迁移耦合模型及其有限元分析

建立了一种饱和-非饱和的双重孔隙-裂隙介质热-水-应力-迁移耦合三维模型,其特点是应力场和温度场是单一的,但具有不同的孔隙渗流场、裂隙渗流场和孔隙浓度场、裂隙浓度场,以及可考虑裂隙的组数、间距、方向、连通率和刚度对本构关系的影响,并研制出相应的三维有限元程序.通过与已有算例的对比,验证了该模型和程序的可靠性.针对一个假定的高放废物地质处置库,就岩体和缓冲层均为非饱和介质和放射性核素泄漏的情况进行了数值分析,考察了岩体中的温度、负孔隙水压力、地下水流速、核素浓度和正应力的状态.结果显示,缓冲层中的温度、负孔隙水压力及核素浓度呈现非线性的变化、分布;尽管裂隙水饱和度平均仅为孔隙水饱和度的1/9,但因裂隙的渗透系数比孔隙的渗透系数大4个数量级,故裂隙中地下水的流速约是孔隙中相应值的6倍;应力分布集中的区域位于缓冲层和处置孔壁交界两侧附近.     

应力腐蚀和压力溶解对双重孔隙岩体中THM耦合影响的有限元分析

引入裂隙开度的应力腐蚀和压力溶解模型,针对一个假设的位于非饱和双重孔隙-裂隙岩体中的高放废物地质处置库,拟定2种计算工况：（i）裂隙开度随应力腐蚀和压力溶解而变化（基岩的孔隙率亦是应力的函数）;（ii）裂隙开度和基岩的孔隙率均为常数,进行热-水-应力耦合的二维有限元分析,考察了岩体中的温度、裂隙开度的闭合速率、闭合量、孔（裂）隙水压力、地下水流速和应力的变化、分布等情况.结果主要显示：应力腐蚀引起的闭合速率要高于压力溶解引起的闭合速率6个数量级,且2种因素产生的闭合速率随时间先增加后减少,并趋于稳定;随着时间的推移裂隙开度由初始值下降并趋于残余值,而粗糙面接触率亦由初始值上升并趋于其名义值;当考虑应力腐蚀和压力溶解时,近场的负裂隙水压力上升很高;2种工况的岩体中的应力量值及分布差别很小.