椭圆-抛物线双屈服面模型参数灵敏度分析

利用两种土石料三轴CD试验数据确定两组椭圆-抛物双屈服面模型参数,研究了模型参数? 、KG、n、h、m、M1、M2和pr分别单独变化（而其他参数保持不变）对模型所反映的土的应力-应变关系的影响。结果表明,土体的体积应变对参数h、m比较敏感;? 、KG、n对剪胀性土的体变影响较大,而对剪缩性土体变影响不明显;? 、KG、n、h和m对模型反映的强度没有影响;pr变化对体变和强度影响均较小。M1和M2对强度和体变均有显著影响,且M1和M2中的小值决定模型的强度。定义了灵敏度,即一种用于表示模型参数对其所反映的应力-应变的影响程度的参数。分析表明,不同的参数,其灵敏度差异很大;不同的土类,参数灵敏度差异也较大。实际应用时,对灵敏度大的参数取值应特别慎重。     

连云港海相软土流变特性试验及双屈服面流变模型

为研究连云港海相软土流变特性,进行了三轴流变试验。根据对试验结果的分析,可以将连云港海相软土作为弹-粘塑性体来研究;将广义Bingham模型和椭圆-抛物线双屈服面模型相结合,建立了新的流变模型来描述其流变特性,得出了各参数,并验证了该模型的适用性。     

基于SFG模型的统一屈服面本构模型与试验研究

将SFG模型中的加载湿陷屈服面(LC屈服面)与BBM模型中提出的吸力增加屈服面(SI屈服面)统一成为单独光滑的屈服面,基于弹塑性理论框架,推导出相应的弹塑性本构模型。通过非饱和固结直剪试验对江西红土的模型计算参数的应力-应变曲线与试验结果进行对比。结果表明:模型对多种应力路径下非饱和三轴试验的计算结果与相应试验的测定结果较为吻合,验证了模型能较好地模拟非饱和土在多种应力路径下的体应变行为。统一屈服面简化了加载湿陷屈服面与吸力增加屈服面之间的耦合分析,解决了因双屈服面的不连续而造成计算不便问题,拓宽了SFG模型的适用性。     

岩体结构面新型本构模型研究

对于涉及需要单独考虑岩体结构面的工程岩体结构分析,采用能反映岩体结构面主要力学特性的合理的本构模型是取得合理解答的关键问题之一。针对在经典连续介质力学理论框架内建立岩体结构面本构模型的缺点,基于岩体结构面的实际受力变形特性,采用直接法建立了一种新型岩体结构面本构模型。所建立的模型依据岩体结构面切向应力变形曲线及剪胀曲线的实际特征,将其分为峰前线性段、峰前非线性段以及峰后软化段,并分别给出了用于描述岩体结构面变形和强度等主要力学特性的数学模型,进而推导建立了结构面各变形阶段的增量型本构模型。最后,编写相关计算程序,采用所建立的新型本构模型以及被广泛采用的Plesha模型对经典的岩体结构面直剪试验成果进行拟合分析。结果表明,所建立的新型本构模型能更为合理的描述岩体结构面的主要力学特性,且模拟能力优于Plesha模型     

加卸载应力作用下煤岩渗透率演化模型研究

为研究开采过程中煤岩力学行为及渗透率演化规律,运用含瓦斯煤热-流-固耦合渗流伺服试验系统,进行了不同加卸载条件下原煤力学及渗流试验,分析了加卸载应力作用对煤岩变形及渗流规律的影响,得到了剪胀角随塑性剪切应变的变化关系,发现了塑性剪切应变在1.6%左右会出现剪胀角急剧变化的现象。根据试验现象和结果,考虑煤岩结构对渗透率的影响,对煤岩在弹性阶段和屈服损伤后的结构进行简化,基于火柴棍模型及渗流理论分析,从应变的角度出发,探讨了加卸载应力对煤岩渗透率的影响,建立了两个不同阶段的渗透率模型（即弹性阶段和损伤阶段）,基于不同阶段的渗透率模型和剪胀角规律,构建了煤岩全过程渗透率模型。所构建的渗透率模型与试验结果对比吻合效果较好,验证了该模型的适用性,可以为实现煤与瓦斯共采提供一定的理论依据。     

堆石料状态相关弹塑性本构模型

堆石料的强度变形特性与初始孔隙及应力状态等因素相关。建立了能够预测不同初始孔隙与初始围压影响的堆石料弹塑性本构模型。在剑桥类本构模型框架内,模型能够反映随着孔隙与围压的增大,变形特性由剪胀趋于剪缩的规律。模型采用了基于颗粒体材料细观结构变化的屈服函数和非关联流动准则,提出了能够反映堆石料正常固结线不唯一的硬化参数。为了反映状态相关性,假定堆石料存在唯一的临界状态面,探讨了考虑状态相关性需要满足的数学条件,从而对剪胀方程与硬化参数进行了修正。提出了基于粒子群优化算法的模型参数快速确定方法,将某筑坝堆石料不同初始孔隙比与围压条件下模型预测结果与三轴试验结果对比,验证了模型的合理性。     

土体弹塑性本构模型研究概述

土的本构模型是对土体进行研究和土工数值分析的关键。本文对土体弹塑性本构模型的发展过程加以介绍，并侧重介绍了近年来正在发展中的模型。     

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

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

基于椭圆-抛物双屈服面模型的广义泊松比研究

适用性较好的土体模型除了能较好地模拟土体的应力-应变关系外,还应该能较好地反映土体的侧向变形,即合理计算广义泊松比。在椭圆-抛物双屈服面模型的基础上,建立了反映侧向变形特性的广义泊松比公式,并利用广义泊松比来研究砂-聚苯乙烯颗粒轻质填料的变形规律。文中研究了包括材料配比在内的主要参数对广义泊松比的影响,发现与剪胀有关的参数对广义泊松比的影响较大。基于砂-聚苯乙烯颗粒轻质填料的三轴固结排水剪切试验数据,给出了相关模型参数,结合该材料的剪切强度、变形特性与建立的广义泊松比公式,对该泊松比能否恰当反映土体的变形特性（剪缩、剪胀和侧向变形）进行验证,结果表明广义泊松比可以反映土体的变形特性,尤其是侧向应变及体积应变规律。     

节理岩体等效流变损伤模型及其在卸载边坡中的应用

重大建设项目对施工过程中岩土体稳定性提出了更高的要求,常以一种动态设计施工模式来应对工程体的各类突发状况与病害。流变损伤模型既能反映岩土体在施工过程中变形的时效发展,又能反映其力学性质的时效劣化,继而能较准确地掌握工程体动态稳定性。基于以上考虑,为体现卸载边坡工程在卸载回弹阶段的瞬时塑性特征和时效演化阶段的黏塑性特征,因此,在流变模型中引入加载塑性元件和黏塑性元件,建立了复合黏弹塑（弹-黏-黏弹-黏塑-塑）模型,室内岩石压缩（卸载）蠕变试验证明了该流变模型的合理性,并对其参数进行辨识。在此基础上,从几何研究方法出发,引入反映节理分布的初始损伤张量及一种等效的依据黏塑性偏应变推导出的损伤演化方程,最终建立了一种新型的节理岩体等效流变损伤模型。将此模型用于川东红层某软硬岩互层型路堑边坡的卸载分析,结果表明：随不同的开挖阶段,易损部位（软岩集中段、软岩深埋段、软硬交接硬岩段）在瞬时卸载回弹阶段的塑性损伤和时效演化阶段的黏塑性损伤逐渐积累,边坡浅表部逐渐出现卸载损伤（松弛）带,在损伤累积中边坡各部位蠕变速率呈不同程度的增长。计算结果较好地反映了边坡变形、损伤发展与动态稳定性特征,其研究结果对于指导支护时机及相应的信息化施工具有一定的意义。     

锦屏大理岩卸荷本构模型与数值模拟研究

锦屏二级水电站深埋引水隧洞位于高地应力区,对取自隧洞内的大理岩进行常规三轴压缩试验、峰前卸围压试验。基于大理岩的峰前卸荷试验,研究峰前卸荷条件下岩石的变形、参数及破裂特征,建立峰前卸荷条件下岩石的幂函数型Mohr强度准则,给出数值仿真试验中大理岩的峰前卸载条件,利用有限差分程序FLAC3D建立了数值仿真模型,对模型进行了计算与分析,最后与室内试验结果对比,得到峰前卸荷状态下大理岩力学特性的相关规律。研究结果表明,幂函数型Mohr强度准则可较好地反映峰前卸荷条件下岩石的强度特性;卸围压初期,岩石试件呈现出明显的弹塑性特征;卸荷过程中岩石试件发生破坏时的峰值强度和围压降值随着围压增加而变大。研究结论揭示了锦屏大理岩卸载应力路径下的力学特性,具有一定的理论参考价值。     

软岩三轴卸荷流变力学特性及本构模型研究

以典型软岩－泥质粉砂岩为研究对象,进行了不同应力水平下的恒轴压、分级卸围压室内流变试验。试验结果表明,软岩在卸荷条件下的轴向及侧向流变变形较大,各向异性显著。在恒定围压较高时,轴向流变变形较大,而随着围压逐级降低,侧向流变变形发展较轴向更快,试样在破裂围压下的侧向变形远大于轴向。通过深入分析软岩卸荷流变试验成果可知,线性Burgers流变模型能较好地反映各级围压下流变曲线的衰减流变及稳态流变阶段,但对于破裂围压下的非线性加速流变阶段无法描述。基于上述分析,建立一个非线性损伤流变模型,采用Levenberg-Marquardt非线性优化最小二乘法对试验结果进行拟合,并获得软岩的非线性流变参数。经比较,拟合计算与试验曲线吻合程度很高,说明该流变模型能较好地反映软岩卸荷流变各阶段的曲线特征。     

岩体卸荷力学特性的数值仿真研究

岩体卸荷力学性质主要表现为非线性和各向异性等特征,与常规岩体试验结果和假定存在着差别,采用数值仿真的方法对卸荷岩体力学特性进行研究,阐述卸荷岩体力学特性研究的重要性;给出了数值仿真分析研究的岩体加、卸载条件;利用大型有限元软件ADINA建立数值仿真模型,对模型进行计算分析,并和已有试验成果进行对比,得出了卸荷岩体力学特性的有关规律。     

典型基坑开挖卸荷路径下软土三轴流变特性研究

模拟基坑开挖下的卸荷路径,进行了不同卸荷条件下包括围压、卸荷路径和排水的卸荷流变试验。试验结果表明：软土在竖向卸荷试验中的回弹变形可分为瞬时回弹变形和滞后回弹变形;在应力控制式三轴试验中,围压和轴压同时卸荷时的初始变形速率较大;仅卸轴压时回弹滞后效应明显,而围压、轴压同时卸荷时由于土体产生剪切变形引起负孔压,使土体回弹滞后效应不明显;在相同卸荷路径下,排水剪切时由于土体吸水膨胀使得流变特性更明显些,且产生的负孔压最终会消散。根据试验成果,建立Merchant流变模型参数与卸荷条件的线性函数关系,并将函数关系代入Merchant模型,最终得到能考虑围压及卸荷路径影响的软土卸荷流变经验模型。研究结果为软土流变变形模型研究及卸荷后基坑变形研究提供了参考。     

冻土模型试验试样冻结时间计算分析

快速简单地确定试样的冻结时间,对模型试验的设计具有重要意义。首先进行了不同干密度冻土模型试样整体冻结达到预设负温的试验,发现试样整体降到预设负温的时间随着干密度的增大而略有增加。然后利用COMSOL Multiphysics软件,在不考虑水分迁移的条件下,对不同干密度粉砂土试样和不同尺度条件下的正方体亚黏土试样的冻结时间进行数值模拟。将模拟结果与试验结果进行对比,发现简化条件下不同干密度的粉砂土试样能反应试样在冻结过程中冻结时间的主要特征,模拟结果中达到预设负温的时间与试验的实测结果基本吻合。正方体试样数值模拟的结果表明在同一干密度、同一含水量条件下,试样降到预设负温的时间随着体积的增大呈现幂函数增大的规律。所以,在进行模型试验设计时,可先根据试验条件进行简单的数值分析,以较合理、经济的安排试验时间,同时也可保证试样冻结质量。     

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

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

基坑土体侧向卸荷真三轴试验研究

对武汉地区具有代表性的粉质粘土,用基坑开挖过程中坑周土体的应力路径在真三轴上进行模拟试验,得到平面应变条件下固结不排水卸荷试验结果,将其整理成应力-应变关系曲线、孔隙水压力-竖向应变关系曲线、中主应力系数-竖向应变关系曲线和卸荷试验应力路径图。试验结果表明,竖向( )压力不变,侧向( )卸荷时,其应力-应变关系随固结压力的增加,由应变硬化型向应变软化型转化;固结压力较低时,表现为剪胀,固结压力较高时,表现为先剪缩,后剪胀;当中主应力( )始终大于但接近于侧向压力( ),按所定义的中主应力系数bd值较稳定;有效应力路径随固结压力增大产生偏转。这些结论为进一步研究粘性土中主应力的作用奠定了基础。     

基坑土体卸荷平面应变试验离散元数值分析

对基坑开挖影响范围内土体的应力路径进行平面应变试验离散元数值模拟,以研究结构性与卸荷形式对坑周土体宏微观力学特性的影响。首先,将一个描述土颗粒间胶结效应的简单三维胶结接触模型植入三维离散元软件PFC3D;其次,对初始K0固结状态的重塑土、结构性土试样分别进行常规三轴以及平面应变条件下4种不同卸荷应力路径的离散元模拟;最后,对经历不同卸荷形式的坑底土体单元进行再加荷模拟。模拟结果表明,在卸荷过程中,被动区土体峰值强度以及破坏时竖向应变随卸荷比增大而增大,且其强度小于主动区土体强度;在卸荷、再加荷过程中,被动区土体峰值强度随卸荷比增大而增大,但均小于不卸荷而直接加荷条件下的峰值强度;由于结构性的存在,土体由应变硬化向应变软化过渡,且强度增长;结构性与卸荷形式显著影响土体体积改变。在微观尺度,增大卸荷比或结构性均会增大垂直大主应力方向的平面上的法向接触力,进而提高其强度。