三轴压缩下粉砂质泥岩蠕变本构模型研究

依据粉砂质泥岩三轴压缩蠕变试验结果,分析了岩石的蠕变力学特性。选取线性黏弹性Burgers模型来描述岩石的衰减蠕变和稳定蠕变特性。引入非线性黏塑性元件,将其与Burgers模型串联起来,建立一个新的六元件非线性黏弹塑性Burgers蠕变本构模型。采用非线性Burgers模型来描述岩石的加速蠕变特性。应用Levenberg-Marquardt算法对蠕变试验曲线进行辨识,得出各级应力水平下岩石轴向与径向的蠕变模型参数。对比试验曲线和模型拟合曲线,二者吻合较好,表明线性Burgers模型与非线性Burgers模型可以准确描述岩石3个蠕变阶段的力学特性。模型参数辨识结果表明,在同一级应力水平下岩石轴向与径向的蠕变模型参数各自独立并不统一,模型参数取值时应考虑岩石蠕变的各向异性以及轴向与径向蠕变量、蠕变速率大小的差异。在前7级应力水平下,模型参数应取以轴向应变辨识得到的Burgers模型参数;在第8级应力水平下,模型参数应取以径向应变辨识得到的Burgers模型参数;在破裂应力水平下,模型参数应取以径向应变辨识得到的非线性Burgers模型参数。     

厦门海底隧道围岩流变特性及其特征曲线

厦门海底隧道是我国建设的第一条海底公路隧道,其穿越的F1,F2,F3,F4共4条断层破碎带,洞体围岩软弱、破碎,流变属性十分明显。对该类岩石进行了室内三轴压缩流变试验。试验结果表明,幂律型蠕变模型可以用来较好地模拟该类岩石除第三阶段蠕变以外的蠕变行为,并得到幂律型蠕变模型的相关流变参数。将幂律型蠕变方程和Drucker-Prager屈服准则相耦合来描述岩体的非线性黏弹塑性特性,应用数值方法研究了该隧道围岩的变形特性及其特征曲线。     

陕南泥岩三轴压缩蠕变试验及其数值模拟研究

在对泥岩进行室内三轴流变压缩蠕变试验的基础上,确定了泥岩Burgers流变模型本构方程中的各参数数值,描述了陕南勉宁地区泥岩的流变特征。结果表明:随着荷载的增加,泥岩出现明显的三个阶段的蠕变曲线,试块的总体蠕变量较大。利用FLAC3D软件对三轴压缩蠕变试验(围压1MPa)进行了模拟计算,确定泥岩计算机模拟时的流变计算参数。在此基础上,分别对泥岩不同围压和不同试件形态状态下的三轴压缩蠕变试验进行了数值模拟计算,结果表明:圆柱体试件的应变量数值整体小于立方体试验试件的应变量数值。     

软弱破碎围岩隧道炭质页岩蠕变特性试验研究

采集贵广铁路天平山隧道典型软弱破碎围岩岩样,在室内开展三轴常规压缩和三轴压缩蠕变试验。基于岩石试验成果,分析炭质页岩的蠕变力学特性,低应力水平条件下岩石仅出现减速蠕变和等速蠕变阶段,在更高的应力水平条件下,岩石并未出现加速蠕变的现象,但围压在一定程度上提高岩石的屈服强度,弹性模量亦有随着围压增加而增大的趋势。据此,建立依托工程炭质页岩蠕变全过程的黏弹塑性应变软化蠕变力学模型,推导该蠕变力学模型在三维应力状态下的本构方程,并确定模型的参数,运用最小二乘法对模型的参数进行辨识。通过与试验曲线对比显示,所建蠕变模型能够描述依托工程炭质页岩的蠕变力学特征,可以用来研究该隧道软弱破碎围岩的蠕变性质和稳定性。     

三轴压缩蠕变试验下石英云母片岩各向异性蠕变特性研究

对石英云母片岩进行三轴压缩蠕变试验,研究丹巴水电站石英云母片岩的三轴蠕变特性及其各向异性特性。按轴向荷载与层理面关系加工成平行组和垂直组2组试件,开展分级加载方式进行不同围压条件下的蠕变试验。试验研究表明,石英云母片岩具有蠕变特性,包括瞬时变形、衰减蠕变、稳定蠕变和加速蠕变阶段;采用带Kachanov蠕变损伤律的蠕变模型来描述石英云母片岩的蠕变特性,并进行蠕变参数辨识,拟合结果显示此模型能很好地描述石英云母片岩的蠕变特性;根据试验结果获得2组试件的长期屈服强度、破坏形态、瞬时变形参数和稳定蠕变速率,分析表明石英云母片岩的蠕变力学特性具有明显的各向异性特性。层理面与轴向荷载垂直的试件较层理面与轴向荷载平行的试件的强度、弹性剪切模量、体积模量和黏性系数相对较大,表现出较高的抗变形和抗破坏能力;平行组以剪切破坏为主,垂直组破坏时出现侧向鼓胀现象,显示出延性破坏的特点;2组石英云母片岩试件的瞬时应变量和稳定蠕变速率都随着应力水平的提高而增大。     

温度周期作用下大理岩三轴蠕变试验与理论模型研究

为了解温度周期作用对大理岩蠕变力学特性的影响,采用温度-应力-渗流耦合的岩石力学三轴流变试验系统对大理岩进行了不同温度及温度周期变化下三轴压缩蠕变试验。基于试验结果,研究了大理岩试样在不同温度及温度周期变化下轴向应变随时间的变化规律,分析了不同温度及温度周期变化下蠕变的变形机制,掌握了温度对大理岩三轴蠕变影响的基本规律,并且建立了考虑温度及温度周期变化影响的改进的Burgers模型。理论值与试验结果非常接近,表明考虑温度及温度周期变化的改进Burgers模型适合描述大理岩的黏弹性蠕变特征。为长期暴露于外部环境下的岩质边坡的长期稳定性分析提供了重要的参考依据。     

分散性土蠕变特性试验研究

采用三轴蠕变仪对宁木特水利枢纽工程黏土心墙的分散性土进行三轴固结排水试验,研究分散性土在不同初始条件下的蠕变变形规律,并获取蠕变模型及蠕变参数。试验结果表明,分散性土具有明显的非线性蠕变特征,蠕变过程可以分为瞬时弹塑性变形、衰减蠕变变形和稳定蠕变变形3个阶段;含水率、围压和偏应力水平对分散性土蠕变速率和蠕变量有显著的影响;应力-应变等时曲线具有明显的非线性特征。采用Singh-Mitchell蠕变模型和Mesri蠕变模型描述分散性土的蠕变特性,与试验数据对比表明,Mesri蠕变模型能够比较准确而方便地描述分散性土的蠕变特性,且模型具有参数少、适用性较强的特点。     

抛物线型屈服面人工冻土蠕变本构模型研究

利用自行研制的W3Z-200型冻土三轴压缩蠕变仪,对淮南地区深表土冻土进行了三轴蠕变试验。根据试验结果,提出以标准黏弹塑性模型为基础,并耦合温度自由度,应用抛物线型屈服准则对模型中的黏塑性项进行改进,建立了新的本构模型来描述高围压复杂应力状态下冻土的蠕变变形特征。基于ADINA有限元软件平台,开发了冻土本构子程序。对淮南某冻结井筒开挖过程进行了数值模拟,计算结果验证了冻土抛物线型黏弹塑本构模型的正确性和合理性。     

基于模型理论和遗传蠕变理论的改良膨胀土非线性流变研究EI北大核心CSCD

为研究膨胀土非线性流变,取水泥改良膨胀土样进行长期一维压缩试验,得到了改良膨胀土在不同应力水平下的应力-应变等时曲线簇。将曲线反映的非线性流变分解为线性黏弹性蠕变变形、线性黏塑性蠕变变形和非线性黏塑性蠕变变形三部分。基于模型理论和遗传蠕变理论建立了各部分流变模型,其中在线性黏弹性蠕变变形和线性黏塑性蠕变变形部分采用模型理论建立元件模型,在非线性黏塑性蠕变变形部分运用遗传蠕变理论建立积分蠕变方程,结合试验数据拟合出材料和方程参数,用该模型得到的蠕变曲线与实测结果能较好吻合。     

基于Weibull分布的红砂岩三轴蠕变试验及模型研究

当前,地下工程围岩蠕变问题仍然存在,蠕变理论需要进一步丰富。岩石蠕变实质上是损伤不断积累的过程,针对蠕变条件下岩石损伤演化情况,文章采用TAW-2000多功能三轴伺服试验系统对取自四川乐山依卜隧道的红砂岩进行三轴蠕变试验,分析不同围压下试样蠕变变形规律,同时以西原模型为基础,结合Weibull分布和Perzyna黏塑性理论,建立一种改进的可以描述岩石蠕变破坏全过程的黏弹塑性蠕变模型。通过划分蠕变阶段来定义临界点损伤变量,从而更为准确地确定加速蠕变启动时间。得出如下结论：（1）模型曲线与试验数据具有良好的一致性,验证了模型的准确性与合理性,说明基于Weibull分布建立的红砂岩黏弹塑性蠕变模型是可行的;（2）基于Perzyna黏塑性理论,建立了可以更加准确的描述加速蠕变的黏塑性应变表达式;（3）文章建立的基于Weibull分布和Perzyna黏塑性理论的三轴损伤蠕变模型能够较好的描述岩石蠕变全过程,克服了西原模型不能描述加速蠕变的缺点。本研究通过定义不同蠕变阶段的临界点损伤变量更好的反映了岩石蠕变变形与损伤之间关系,丰富了岩石类材料的蠕变本构理论。     

基于改进Bingham模型的软岩参数非定常三维非线性黏弹塑性蠕变本构研究

传统线性元件组合模型难以描述岩石非线性加速蠕变阶段特性,给工程应用带来了巨大困难。提出一种新的应变触发非线性黏壶元件,并将其与Bingham模型串联,建立了可以描述岩石等速和加速蠕变特性的一维本构模型,进而推导出相应的三维蠕变本构模型,再将三轴压缩过程中岩石弹性模量的衰减方程引入该三维蠕变本构方程,得到一个能反映岩石蠕变全过程特性的三维非线性黏弹塑蠕变本构模型,根据泥岩蠕变试验数据对该模型参数进行了辨识。模型理论计算值与泥岩蠕变试验数据的对比表明改进的Bingham模型不仅可以充分反映软岩的等速蠕变特性,还可以很好地描述软岩的衰减和加速蠕变规律,且模型元件少,组合形式简单,为软岩非线性蠕变本构模型研究提供了新的借鉴。     

基于应变软化指标的岩石非线性蠕变模型

研究并提出一种以应变软化指标为基础的岩石非线性蠕变模型。首先,基于岩石三轴压缩试验结果,分析提出一种峰后应变软化指标R1来描述岩石峰后力学参数软化与塑性变形间的关系。接着,通过分析加速蠕变阶段与峰后段的联系,建立加速蠕变应变软化指标R2,并以此为基础构建非线性黏塑性体。然后,将其与Hook体和Kelvin元件串联,组建非线性蠕变模型,并将模型嵌入ABAQUS有限元程序。最后,针对砂岩和泥岩三轴蠕变试验建立数值模型,模拟曲线与试验曲线吻合良好,说明所建模型适用于硬岩/软岩加速蠕变现象的模拟。其中,参数 对加速蠕变曲线形态起调节作用,显示模型模拟脆性?延性流变破坏的特点。另外,模型参数均可由常规压缩破坏试验以及蠕变试验确定,易于获取。     

冻土的蠕变及蠕变强度

基于试验，本本讨论分析了冻土三轴蠕变规律及温度和围压对冻土蠕变强度的影响，给出了土蠕变及蠕变强度随时间降低的方程式，进而提出了冻土蠕变强度的物型屈服准则。通过改变参数可将冻土的瞬时强度准则与冻土的瞬时强度准则与冻土的蠕变强度准则统一为同一个方程描述。     

基于分数阶导数的非定常蠕变本构模型研究

针对锦屏二级水电站深埋长大隧洞围岩在开挖过程中呈现明显的流变劣化特征,基于分数阶微积分理论,用Abel黏壶替换西原模型中的牛顿黏壶,推导出基于分数阶导数的蠕变本构方程,同时考虑岩石在蠕变过程中,尤其是所受应力超过其长期强度时,蠕变参数是非定常的特征,把蠕变参数的非定常性引入本构方程中,建立了基于分数阶导数的非定常蠕变本构模型。结合锦屏二级水电站引水隧洞主要围岩之一的大理岩的室内剪切蠕变试验的结果,对模型中的参数进行反演,结果表明,基于分数阶导数的非定常蠕变本构模型能够在蠕变初期和拐点处与试验结果很好地吻合,并且可以克服西原模型不能反映蠕变第3阶段的不足。通过对模型中参数的敏感性进行分析,得出分数阶阶数、非定常参数对岩石蠕变的影响规律。研究结果表明,提出的基于分数阶导数的非定常蠕变模型能够很好地反映大理岩的蠕变全过程。     

膨胀土卸荷蠕变特性及其非线性蠕变模型

鉴于膨胀土滑坡往往表现为长期性、渐进性等与时间相关的特性,利用GDS应力路径三轴仪对膨胀土进行了三轴卸荷蠕变试验。试验结果表明：当偏应力较小时,膨胀土的蠕变曲线仅出现瞬时变形和衰减蠕变;当偏应力达到一定值时,其蠕变曲线也呈现衰减蠕变、稳态蠕变和加速蠕变3个阶段,但其加速蠕变阶段的特征与一般岩土体不同,其蠕变速度近乎常数。膨胀土的应力?应变等时曲线显示,膨胀土卸荷蠕变具有非线性特征,且其非线性程度与蠕变时间和应力水平相关,蠕变时间越长、应力水平越高,其非线性程度越高。基于非线性流变力学理论,提出了一种非线性四元件蠕变模型,将标准线性体与一个非线性黏壶串联,该模型可描述等围压三轴压缩应力状态下膨胀土轴向应变随时间的演变规律。根据膨胀土卸荷蠕变试验结果,采用曲线拟合法对三维非线性模型的参数进行反演识别。拟合曲线和试验曲线对比显示,两者吻合良好,说明该模型可以很好地描述膨胀土的蠕变特性。此外,基于该蠕变模型获取了膨胀土的临界破坏应力,其与常规剪切破坏应力的比值随着固结压力的减小而减小,表明越接近坡面的土层越容易发生蠕变破坏。     

引红济石引水隧洞围岩力学特性研究

引红(红岩河)济石(石头河)引水隧洞围岩断层带破碎松软,岩石强度低,自稳能力差,易产生大变形。以上不良地质条件使其隧洞施工过程中经常会遇到TBM卡机等一系列特殊问题。针对引红济石引水隧洞施工中存在的软岩隧洞大变形问题,首先开展了X射线衍射、崩解试验,分析了试样的组成成分、黏土矿物含量对崩解性影响;再通过一系列单轴压缩试验、三轴压缩试验和蠕变试验研究了该类岩石不同应力条件下的变形破坏特征及蠕变特性。试验结果表明,该类岩石含有大量的黏土矿物(33.49%),对水比较敏感,遇水膨胀易崩解,导致岩体软化;试样具有较大塑性压缩变形,其应力-应变曲线为应变强化型,且没有明显的峰值。基于试验研究成果和现场监测数据,对该软岩隧洞大变形机制进行了分析,并提出了在围岩与管片之间安装聚氨酯缓冲层的新型支护方案,通过数值计算对支护方案的合理性进行了验证。分析结果表明:聚氨酯缓冲层可以很好地吸收围岩形变压力,避免应力集中带来的管片错台,从而大大减小管片上破坏区的产生。研究成果对该类岩体中隧洞的设计施工以及长期稳定性分析具有重要的参考作用。     

Anisotropic viscoplastic modeling of rate‐dependent behavior of clay

The objective of this study is to derive an effective stress‐based constitutive law capable of predicting rate‐dependent stress–strain, stress path and undrained shear strength and creep behavior. The flow rule used in the MIT‐E3 model and viscoplasticity theory is employed in the derivation. The model adopts the yield surface capable of representing the yield behavior of the Taipei silty clay and assumes that it is initially symmetric about the K₀‐line. A method is then developed to compute the gyration and expansion of the loading surface to simulate the anisotropic behavior due to the principal stress rotation after shear. There are 11 parameters required for the model to describe the soil behavior and six of them are exactly the same as those used in the Modified Cam‐clay model. The five additional parameters can be obtained by parametric studies or conventional soil tests, such as consolidation tests, triaxial compression and extension tests. Finally, verification of the model for the anisotropic behavior, creep behavior and the rate‐dependent undrained stress–strain and shear strength of the Taipei silty clay is conducted. Copyright © 2010 John Wiley & Sons, Ltd.