盐岩非线性Burgers模型及其参数识别

基于非线性流变力学理论,对常规黏滞体元件进行改进,提出一种新的变系数非线性黏滞体,其黏滞系数随时间增加而逐渐增大。将该非线性黏滞体替换常规Burgers模型中的线性黏滞体,建立非线性Burgers模型即NBurgers模型。该模型在一定的条件下可退化为常规Burgers模型。对盐岩试件开展三轴压缩蠕变试验,并利用试验结果,根据非线性最小二乘法基本原理对NBurgers模型的参数进行反演识别。拟合曲线和试验曲线对比显示:Nburgers拟合效果良好,误差较小,拟合曲线非常逼近于试验曲线,且其与试验曲线的吻合程度要明显高于常规Burgers模型,说明该模型能够更好地描述盐岩的黏弹性蠕变特性。     

基于改进分数阶黏滞体的岩石非线性蠕变模型

通过将表征应力水平对岩石非线性蠕变特性影响的函数引入到常规分数阶黏滞体的本构关系中,提出一种改进的分数阶黏滞体;基于岩石非线性流变力学理论,提出一个能够描述岩石加速蠕变的非线性黏滞体,其黏滞系数随应力水平提高和蠕变时间延长而逐渐减小;将改进的分数阶黏滞体、非线性黏滞体与基本弹性体和塑性体进行组合,建立一个新的4元件非线性黏弹塑性蠕变模型,给出模型的蠕变方程,并利用岩石蠕变试验结果对模型合理性进行验证。研究结果表明:该模型能够较好地描述岩石蠕变全过程的3个阶段,且拟合曲线和试验曲线吻合良好,误差较小。     

冻土蠕变的非线性模型研究

为了研究冻土的蠕变特性,在传统的伯格斯模型中引入一种考虑时间和应力水平影响的非线性黏滞体,代替与弹簧并联的理想牛顿流体,并在此基础上串联一个由非定常、非牛顿黏壶和塑性体并联的新的非线性黏塑性体,从而建立了新的非线性模型,并推到了三维应力状态。随后采用该模型借助Origin软件对冻土蠕变的试验数据进行拟合。结果表明:1该模型所拟合的曲线与试验数据吻合性较好,误差很小,而且适用于不同的温度。2该模型不仅可以很好地描述冻土蠕变的前两个阶段,而且对于加速阶段也适用,这对于高应力下的冻土蠕变分析很有意义。     

武汉地区软土固结蠕变耦合性状模型研究

运用数理统计方法对武汉地区软黏土固结蠕变试验资料进行处理,得出了各试验曲线的经验公式;通过对经验公式的分析,建立了非线性弹性黏性固结蠕变模型,其中,以Duncan-Chang模型描述土体的瞬时弹性,以Burgers模型的一部分描述土体的黏弹性;对Burgers模型中的参数进行了改进,根据固结蠕变模型特征及试验资料,用数理统计方法求出改进的Burgers模型中参数的经验值.     

基于广义Kelvin模型的非定常盐岩蠕变模型

针对川东北油气田盐岩地层普遍出现的套损现象,对该地层段盐岩的蠕变特性进行研究。用非定常黏壶元件替换广义Kelvin模型中的线性黏壶元件,建立非定常广义Kelvin模型;将其与能够描述盐岩稳态蠕变特征的Heard模型进行串联,构建新的四元件NGKH盐岩非线性蠕变本构模型,可以描述盐岩的衰减蠕变阶段和稳态蠕变阶段特征。基于室内三轴蠕变试验结果,利用Levenberg-Marquardt算法对NGKH模型蠕变参数进行识别。基于FLAC~(3D)的UDM接口程序,采用C++编程语言对其进行二次开发;并通过三轴蠕变试验的数值模拟对开发NGKH模型进行验证。研究结果表明:NGKH模型的拟合结果与室内试验结果吻合良好,该蠕变模型能够准确描述盐岩的非线性蠕变特性;通过黏弹性与非线性计算结果,验证了开发NGKH模型程序的正确性和合理性。     

排土场工后沉降及蠕变规律

采用Burgers模型对排土场散粒体蠕变的减速蠕变和等速蠕变两阶段进行描述,利用分层总和法的思想,将垂直填筑的排土场进行分层处理,底层在上覆“自重”变荷载作用下发生沉降变形,分别采用定常和非定常Burgers蠕变模型从理论解析角度推导排土场填筑动态过程中沉降、工后沉降及累计沉降计算公式.以齐大山铁矿排土场监测数据进行实例验算,应用结果表明：非定常Burgers模型和定常Burgers模型拟合的相关系数均较高,非定常Burgers模型沉降最终收敛于5.07 m,定常Burgers模型沉降曲线具有发散性,可见非定常Burgers模型能更好地表述排土场沉降真实工况;将排土场分为十个分层,结合FLAC3D软件的蠕变数值分析计算,得出各单层沉降率的变化规律,即各单层工后沉降量上层沉降值小于下层和中间层,且上层沉降量呈现单调递减变化,越接近排土场顶部单层沉降量越小;中间层沉降量相对下层沉降量要大,其中第5单层的沉降量最大.     

考虑系数非定常性的成都黏土非线性蠕变本构模型

为了研究成都黏土蠕变规律,展开固结不排水三轴蠕变试验,分析了成都黏土变形的非线性特性及参数的非定常特性。结果表明：当应力小于黏土屈服强度时,蠕变以线性变形为主,包括瞬时弹性变形和黏弹性变形;当应力大于屈服强度时,蠕变以非线性变形为主,包括黏弹塑性变形,具有显著的非线性特性;蠕变过程中,黏土的弹性模量和黏滞系数均为应力和时间的函数,具有显著的非定常特性。结合元件模型理论和分数阶导数模型理论的优点,构建了考虑弹性模量和黏滞系数非定常特性的成都黏土非线性蠕变本构模型,并对模型进行了拟合验证分析,发现蠕变模型对蠕变各阶段的拟合度较高,充分发挥了元件模型和分数阶导数模型的优点,可以很好地反映成都黏土的蠕变全过程。     

基于颗粒流程序的非定常西原体模型

为实现岩石试样蠕变全过程的准确模拟,并从细观角度探究蠕变过程中微裂隙的发生和发展规律,在二维颗粒流程序(PFC2D)中开发出具有黏弹塑性特征的西原体流变接触本构模型,进一步提出包含两种非定常元件的非定常西原体模型,推导了模型本构关系和蠕变方程.在PFC2D中调用自定义西原体流变模型,通过参数调试,获得与真实试样具有相同强度特性的数值试样.以室内单轴压缩蠕变试验数据为基础,在Matlab中对模型非定常参数进行拟合反演分析.在此基础上,进行单轴压缩蠕变试验的模拟,计算过程中分别采用定常和非定常两种模型,并对微裂隙进行监测.对比分析结果表明:定常模型仅适用于衰减和稳定蠕变阶段;非定常模型也可用于描述加速蠕变阶段,从而准确模拟蠕变全过程;加速蠕变阶段主要是由微裂隙的加速发展而产生,加速蠕变将导致试样剪切破坏.     

考虑塑性变形的分数阶岩石蠕变损伤模型研究

目前较为经典的岩石蠕变模型元件都是线性的,虽然可以较好地模拟岩石蠕变的前两个阶段,但是却无法对岩石加速蠕变阶段进行模拟,且因为忽略了蠕变初期岩石压密产生的瞬时塑性变形,模型所反映出的瞬时变形小于真实值.本文引用了裂隙压密体元件描述岩石蠕变过程中的瞬时塑性变形,通过使Riemman-Liouville分数阶微积分算子理论构建分数阶黏滞体元件,同时考虑到岩石加速蠕变阶段的岩石损伤对黏滞体的黏滞系数的影响,使得构建的分数阶黏滞体元件可以反映岩石加速蠕变阶段.将裂隙压密体和考虑损伤的分数阶黏滞体与经典的Poyting-Thomson模型进行串联组合组成新的模型,推导了新模型的蠕变方程,最后,将模型的蠕变方程与相关文献中的蠕变试验曲线进行拟合确定模型参数,同时将理论蠕变曲线与蠕变试验曲线进行对比分析,结果表明了本文所建模型的正确性.     

基于互层岩体的非定常黏弹塑性蠕变模型及有限元分析

为研究互层岩体在考虑时效变形作用时的蠕变特性，提出了一种考虑应力、应变阈值的非定常黏弹塑性五元件模型。该模型能够同时描述岩体的瞬时弹性应变、黏弹性蠕变、线性黏塑性蠕变和非线性黏塑性蠕变（加速蠕变），推导得出了其在三维应力状态下的蠕变方程；基于ABAQUS有限元软件完成了UMAT子程序的研发，并通过对比岩体的室内蠕变试验与该模型数值模拟结果来验证模型的适用性。结果表明：蠕变试验结果与数值计算结果相吻合，非定常黏弹塑性模型不仅可以精确描述减速蠕变和稳态蠕变过程，还可以较好地描述岩体的加速蠕变过程，从而验证了该模型的适用性与有效性。