硅藻质软岩的三维弹粘塑性模型分析

基于Perzyna的弹粘塑性理论的基本假设、粘塑性流动法则及殷建华和Graham的三维弹粘塑性模型,得出了在三轴应力状态下的本构关系表达式。殷建华和Graham的模型提出了“等效时间”,“瞬时时间线”,“参考时间线”,“极限时间线”等一系列新的概念。用三轴应力状态下的本构关系表达式对软岩的固结不排水三轴剪切试验进行了模拟计算,并与试验结果进行了比较。结果表明：数值模拟计算结果与试验结果几乎一致。该模型简单、实用、精度较高,可以模拟硅藻质软岩材料与时间有关的应力-应变关系、孔隙水压力与应变的关系,以及有效应力路径等。     

基于应变空间硅藻质软岩的软化本构模型

三轴试验结果表明,软岩具有显著的应变软化特征。正常固结软岩表现为峰值后较缓的软化特征,而超固结软岩峰值强度后的应力-应变曲线呈陡降软化特征。应变空间表述的弹塑性理论在解决大应变和软化问题时比应力空间表述的弹塑性理论更具有优越性。基于应变空间的基本的弹塑性本构方程式,采用Mises剪切屈服准则及相关联流动法则,导出固结不排水三轴条件下的应力-应变本构关系式,并采用不同的硬化函数表达式对软岩在不同围压下的应力-应变曲线进行数值模拟。结果表明,应变空间的弹塑性理论能较好地模拟软岩应变软化特征,其中硬化函数的确定是关键问题之一。通过研究,提出了用于固结不排水状态下正常固结软岩的硬化函数形式。     

超固结土热黏弹塑性本构关系

热黏弹塑性本构模型是描述土在温度（热）和时间（黏）耦合作用下的应力-应变关系的本构模型。在一些新型岩土工程诸如高放核废料地质处置、地热资源开发与贮存的建设中,需要同时考虑温度和时间对土的影响,所以建立一个热黏弹塑性本构模型具有理论和实际意义。将温度变化对黏土体积和强度参数的影响引入笔者之前提出的超固结土等向应力-应变-时间关系,建立了一个等向应力条件下的应力-应变-时间-温度关系。随后,基于该关系推导了屈服面硬化定律,并将其与超固结土统一硬化模型的屈服方程和流动法则结合,建立了超固结土的热黏弹塑性本构模型。最后,使用新模型预测室内试验,证明新模型能够反映时间和温度对土体积、一维压缩曲线和前期固结压力的耦合影响。     

考虑温度影响的软岩弹塑性本构模型

煤矿的深度开挖、核废料的地下储存和能量桩的广泛应用等诸多的岩土工程问题都需要考虑温度对软岩力学特性的影响。为了能较为全面地描述软岩的力学特性,基于上下负荷加载面,引入温度等价应力的概念,在tij应力空间下构建了一个可同时考虑温度效应、中间主应力影响、结构性、超固结性等软岩力学特性的弹塑性本构模型。新模型的所有参数都具有明确的物理含义。通过理论曲线与试验结果进行对比,验证了所提本构模型的正确性。最后通过改变模型参数,对新本构模型的性能进行了分析讨论。计算结果表明：（1）增大超固结比发展控制参数m或者减小结构状态发展控制参数m*,将会提升软岩的剪切强度。（2）随着温度的上升,软岩的剪切强度反而减弱。（3）初始的超固结比越大,软岩的剪胀特性更明显;而初始的结构性较大时,其体积应变在剪切的最后阶段表现为剪缩。     

海相沉积软黏土的弹塑性本构模型研究

在深入探讨海相沉积原状软黏土压缩、变形等力学特性和详细分析加载屈服面随荷载情况变化的基础上,确认了海相沉积原状软黏土的强度、变形特性与结构屈服应力密切相关。即当固结压力小于结构屈服应力时,其力学特性与超固结重塑土的力学特性类似;当固结压力大于结构屈服应力时,其力学特性与正常固结重塑土的力学特性类似。为描述海相沉积原状软黏土的上述力学特性,将姚仰平等提出的超固结重塑土本构模型引入到海相沉积软黏土弹塑性本构模型的构建中。在本构模型构建过程中,考虑了海相沉积原状软黏土具有的抗拉强度及其演化规律,软黏土强度包线的特点及其进一步修正的表达式,使模型更符合海相原状软黏土的强度、变形特性。最后,将3种不同海相沉积软黏土固结排水剪切试验得到的应力-应变-体变曲线与模型预测结果进行对比。比较结果显示,本文提出的弹塑性本构模型能很好地描述海相沉积原状软黏土的剪缩硬化、剪胀软化以及变形的应力水平依存性等力学特性。     

一个描述软黏土时效特性的一维弹黏塑性模型

为了描述软黏土一维应力-应变关系的时效特性,基于Bjerrum的等时间线体系,提出等黏塑性应变率线概念,推导了黏塑性应变率与黏塑性应变增量的关系,建立了软黏土的一维弹黏塑性模型;从理论上分析了新建模型与3种典型的一维弹黏塑性时效本构模型的内在联系,表明新建模型与其他3种模型在本质上是等效的,且形式更简洁,物理意义更明确;利用新建模型对软黏土的固结-蠕变耦合效应、应变率效应、应力松弛效应等时效特性进行了理论分析,并得到了相应的解析解;结合宁波软黏土的一维固结试验,阐述了模型参数的确定方法,并用新建模型对宁波软黏土的固结-蠕变试验、温州软黏土的一维多级等应变率试验、香港海相软黏土的一维应力松弛试验进行模拟,验证了新建模型的有效性。研究结果表明,新建模型能很好地模拟软黏土的一维时效特性。     

《岩土力学》第23卷总目次2002年1～6期(总75～80)

第1期（总75）2002年2月 并行遗传算法对硅藻软岩本构关系的求解……………………………………… 李邵军，冯夏庭，廖红建 (1) 土与结构物接触面损伤本构模型…………………………………………………………… 胡黎明，濮家骝 (6) 冻土单轴压缩动态试验研究……………………………………………………刘增利，李洪升，朱元林等 (12) 均布荷载作用下挡土墙后粘性填土的土压力计算…………………………………………年廷凯，栾茂田 (17) 水平荷载作用下单桩非线性ｍ法试验研究………………………………………王 梅，楼志刚，李建乡等 (23) 卸荷…     

《岩土力学》2002年1～3期EI收录论文情况

从2002年起，我国被美国“EI compendex”数据库检索的科技期刊原寄往美国或荷兰“EI”检索机构的期刊，全部转为由“EI”设立在我国北京“EI”信息部收录。 为了满足2002年在本刊发表论文的作者了解所发表的文章是否被“EI compendex”收录的要求，现将本刊编辑部查询到2002年1～3期被“EI”收录的论文公布如下。其它几期的收录情况一旦有了消息，本刊编辑部将继续在本刊刊出。请作者留意。 论文目次 作 者 页码 *** 第1期 *** 并行遗传算法对硅藻软岩本构关系…     

高地应力下硬岩的本构模型研究

传统的本构模型在模拟高地应力下硬岩破坏的范围和深度方面并不理想。针对高应力下的拉西瓦花岗岩,通过真三轴压缩试验模拟了开挖时应力路径的演化。在摩尔-库仑强度准则、格里菲斯强度准则、德鲁克-普拉格准则、双剪统一强度理论和霍克-布朗强度准则的基础上,考虑十二面单元体主剪面上主剪应力、正应力和静水压力的共同作用,提出了一个有关高应力下硬岩的三剪强度准则。根据试验结果,利用遗传算法全局寻优的功能,搜索出了三剪强度准则的参数。预测样本和与摩尔-库仑强度准则、德鲁克-普拉格准则、双剪统一强度理论的比较表明,该三剪强度准则与试验结果吻合,适用于高地应力下的硬岩。根据试验得到的应力-应变关系曲线,建立了基于三剪强度准则和应变软化的弹-脆-塑性本构模型,并用FLAC3D提供的基于C++的用户自定义模型工具UDM,创建了用户自定义模型的动态链接库,嵌入了FLAC3D软件。采用遗传算法-FLAC方法,搜索得到了基于岩石试样和岩体的本构模型参数。计算结果与实测情况吻合,表明建立的本构模型适用于高地应力下的硬岩,为高地应力下硬岩地下工程的安全性和稳定性分析奠定了基础。     

黏土岩时效变形特性试验与理论研究

高放废物地质处置岩体材料选取中,黏土岩因具有低渗透性、损伤自修复特性、对放射性核素具有良好的吸附作用等优点,被认为是一种合理的高放废物地质处置屏障。以黏土岩为研究对象,从黏土岩短期、长期力学特性等方面开展研究工作,主要内容如下：（1）进行黏土岩短期自然固结试验,确定了黏土岩的基本物理力学参数,主要包括前期固结压力、压缩指数等;（2）进行流-固耦合固结试验,试验结果表明盐水作用对黏土岩力学特性具有重要影响,且流-固耦合过程中的膨胀现象与黏土岩的黏土矿物含量和类型密切相关;（3）通过对黏土岩进行固结流变试验,研究黏土岩渗流-应力耦合作用下的长期力学特性,试验表明黏土岩具有明显的流变特性,且流变现象与载荷密切相关;（4）根据固结流变试验建立黏土岩一维流变本构模型,同时将模型计算结果与试验结果进行对比,分析表明,该模型能够很好地反映黏土岩一维固结流变特性。该研究对我国未来黏土岩高放废物处置库的规划、设计、选址和运营等具有重要的参考意义。     

基于改进遗传算法的岩体结构面幂函数模型参数辨识

沿结构面的剪切滑移是工程岩体主要的破坏模式之一,选择恰当的岩体结构面本构模型是分析岩体结构面破坏的重要一环。当前常用的结构面本构模型有指数型模型、曲线型模型和幂函数型模型,其模型参数都是根据工程实际经验或室内试验所得,不能很好的模拟实际变形情况。指数型模型得出的模拟值偏大,曲线型模型得出的模拟值偏小,幂函数型模型在较低的法向应力下能较好的模拟剪切变形过程,但在较高的法向应力下得出的模拟值偏小。为此,通过改进的遗传算法对岩体结构面幂函数模型参数进行辨识,得出由"等效参数"构成的幂函数模型。实例验证结果说明,通过改进的遗传算法辨识出的由"等效参数"构成的幂函数模型能够很好的模拟实际变形情况,能够很好的适应于不同的工程实际应用。      

岩石本构模型智能辨识方法研究

岩石本构模型是研究岩石力学特征和变形机制的基础,而本构模型或模型中相关参数的识别是本构模型研究中的热点和难点问题。本文基于红板岩室内力学实验数据,分别利用遗传算法、BP神经网络以及遗传规划对红板岩本构模型进行了模式识别,结果表明,遗传算法进行参数识别需要事先假定流变模型的形式,误差较大,而BP神经网络和遗传规划可以一次性同时确定流变模型的结构形式及参数,有效避免模型假定所带来的误差。而遗传规划与BP神经网络相比,具有精度高、收敛快,可视化程度高等特点,为岩石本构参数及模型的智能识别方法的选取提供参考。     

基于HPSO算法的岩石非定常蠕变本构模型辨识

复合微粒群优化(HPSO)是一类随机全局优化技术,具有搜索能力强、收敛速度快、搜索精度高的优点.针对岩石蠕变本构模型非定常参数的辨识问题,利用FLAC软件自带的fish语言实现了HPSO算法对非定常参数的辨识.该方法从非定常参数的随机值出发,以蠕变过程中试件变形的实验值与计算值的误差大小作为适应度函数来评价参数的品质,...     

基于Drucker-Prager准则的岩石弹塑性损伤本构模型研究

大多数岩石材料软化本构模型在硬化函数中引入塑性内变量来表示材料的硬化/软化性质,但并不能反映岩石微裂隙损伤对材料力学性能的影响及单轴拉伸和压缩所表现的初始屈服强度f0与屈服极限fu的差异。基于D-P准则同时考虑塑性软化及损伤软化,建立岩石类材料的弹塑性本构关系及其数值算法。塑性屈服函数采用Borja等的应力张量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;由于岩石损伤软化是微裂隙扩展所导致的体积膨胀引起的,因此,提出用体积应变表征岩石损伤变量的演化,并用回映隐式积分算法编制了岩石的弹塑性损伤本构程序。对单轴压缩及拉伸荷载作用下的岩石材料试验进行数值模拟,结果表明,所提出的岩石弹塑性损伤本构模型可以较好地符合岩石材料的力学特性。     

岩石破坏非连续变形分析弹脆性本构模型的二次开发

岩石非连续变形分析方法DDARF成功实现了模拟岩体裂纹萌生、扩展、贯通、破碎的全过程,但算法只考虑了岩体的线弹性本构模型,没有分析岩体的非线性应力-应变关系。为更符合真实岩体工程,同时扩展岩石非线性本构模型的应用范围,分别采用摩尔-库仑强度准则和最大拉应力强度准则对岩体进行剪切和拉伸破坏判断;对理想弹脆性本构模型进行算法分析,并在VC++平台下程序实现;对“自定义”的岩体非连续变形分析方法弹脆性本构模型与室内单轴压缩试验进行比较分析;将岩体非连续变形分析方法弹脆性本构模型应用于大型地下洞室开挖,并与线弹性本构模型进行对比。结果表明,非连续变形分析方法中自定义岩体弹脆性本构模型是可行的,它能够反映岩体变形的非线性特征,与室内试验吻合度较好;非连续变形分析方法弹脆性本构模型应用于大型地下工程,能更安全且真实地分析洞室围岩的稳定性,进而更好地指导地下洞室的防护措施。     

层状盐岩体的三维Cosserat介质扩展本构模

考虑到深部地下岩土工程一般是较复杂的三维问题,在层状盐岩体二维Cosserat介质扩展本构模型基础上,针对不同岩层交替而成的互层盐岩体,建立了宏观平均意义下考虑细观弯曲效应的三维Cosserat介质扩展本构模型,并给出了模型的初步试验验证。该模型既能分析硬（软）夹层对层状盐岩体刚度的强化（弱化）作用,又能考虑由于夹层的存在而导致的层状盐岩构造体宏观意义下力学特性的各向异性,为层状盐岩体及其他层状介质变形和破损分析提供了一种新的思路和方法。该模型可直接推广到非线性问题。结合数值分析软件可进行层状介质体内复杂洞室稳定性分析。     

Identification of non‐linear stress–strain–time relationship of soils using genetic algorithm

Recognition of non‐linear constitutive rock/soil model from experimental results is often multi‐modal in the large parameter space. A genetic evolution algorithm is thus proposed for its recognition, including that of structure of the model and coefficients in the model. The structure of the model can be firstly determined according to mechanical mechanism if the mechanism is clearly understood or searched by using evolutionary algorithm. The coefficients to be determined are then searched in global optional space. With the new evolutionary algorithm, the non‐linear stress–strain–time constitutive law to describe strain softening behaviours of diatomaceous soil under consolidated and undrained state was recognized by learning stress–strain–time behaviour of an intact sample under consolidated pressure of σ_c=0.1 MPa and strain velocity of_a=0.175%/min. This model gave reasonable prediction for diatomaceous soils under varying consolidated pressures (0.1–3.5 MPa) and strain velocities (0.0044–1.75%/min). It indicates that the methodology proposed in this paper is robust enough and strongly attractive for recognition of non‐linear constitutive model of soil and rock materials. Copyright © 2002 John Wiley & Sons, Ltd.     

岩石弹塑性损伤本构模型建立及在隧道工程中的应用

在实际隧道施工过程中,隧道开挖引起地下岩体应力重分布使得围岩的微裂纹扩展损伤,并伴随有塑性流动变形。在地下水环境中对于孔隙和微裂隙围岩介质受到应力作用时,在内部将产生高孔隙水压力影响岩石的力学性质,也改变了围岩的破坏模式。为了研究损伤引起的刚度退化和塑性导致的流动两种破坏机制的耦合作用,从弹塑性力学和损伤理论的角度出发,同时引入修正有效应力原理来考虑孔隙水压力的作用,建立基于Drucker-Prager屈服准则的弹塑性损伤本构模型;针对该本构模型推导了孔隙水压力作用下弹塑性损伤本构模型的数值积分算法-隐式返回映射算法,分别对预测应力返回到屈服面的光滑圆锥面或尖点奇异处两种可能的情况给出了详细的描述,隐式返回映射算法具有稳定性和准确性的特点;大多数弹塑性损伤模型中涉及参数多且不易确定的问题,采用反分析方法获得损伤参数,解决了损伤参数不易确定的难题;采用面向对象的编程方法,使用C++语言编制了弹塑性损伤本构求解程序,并对所建立的弹塑性损伤模型和所编程序进行了试验和数值两个方面的验证;最后将其在吉林抚松隧道工程中进行应用,模拟了塑性区和损伤区的发展变化。研究结果表明：所建立的弹塑性损伤本构模型能够较好地描述岩石的力学性能、塑性和损伤变化趋势,所编程序能够进行实际工程问题的模拟,对现场施工给予一定的指导。     

岩石时效损伤模型及其在工程爆破中应用

为了实现爆破智能化,降低爆破带来的负面影响,通过数值模拟预估出岩体中爆破诱发损伤的大小以及分布范围在工程上具有重要意义,其中科学合理的爆破损伤模型是关键。基于岩石介质是连续的、各向同性的、且存在初始裂纹的假设,把考虑时效的岩石拉裂损伤演化方程和双线形弹塑性本构结合起来,通过用户自定义简明地嵌入到大型LS-DYNA软件中,并引入了该软件中的侵蚀算法,对半无限岩体中临空面附近爆破漏斗问题进行了重点数值模拟,计算结果比较接近实际。本文方法在工程中将具有一定的参考价值。