岩石流变损伤模型及春有限元程序编制

基于试验资料，考虑了由裂纹扩展导致的岩石扩容，给出了损伤变量及其演化方程，建立了岩石流变损本构关系和有限元迭代格式。用Ｆｏｒｔｒａｎ语言编制了程序，并进行了算例分析。     

岩石疲劳损伤模型的参数估计方法研究

岩石是一种非均质的复杂地质材料,其疲劳损伤演化具有明显的3阶段发展规律。倒S型非线性疲劳累积损伤模型可以很好地描述这一规律。应用倒S型模型开展损伤计算和疲劳寿命分析的先决条件是确定出不同应力水平下的模型参数,而运用Levenberg-Marquardt法进行倒S型曲线拟合的结果表明,不用任何优化方法而随意选取参数初始值时迭代收敛到最优解的概率很低。因此,提出运用有限组合法估算初始值,即在参数可行域中均匀布点,以这些点为起始点分别实施一次Levenberg-Marquardt优化计算,取残差平方和最小的一组参数作为第2次Levenberg-Marquardt拟合计算的初始值。试验结果表明,以有限组合法估算参数初始值,Levenberg-Marquardt迭代计算总能收敛到最优解。     

岩石弹脆性分维损伤本构模型

本文定义岩石构元中破裂面的分维值为各向同性损伤变量,而各个方向上裂纹面的累加量定义为各向异性损伤变量,并根据裂纹发育特征提出了损伤变量演化方程,从而建立起岩石脆性变形破坏过程的分维损伤本构模型。最后,利用该模型对大理岩单轴压缩应力应变曲线进行了模拟,结果说明本文提出的模型是较为合理的。     

拉伸荷载作用下岩石的细观损伤力学模型

本文建立了拉伸荷载作用下岩石的二维、轴对称和三维细观损伤力学的有效场模型，分析了岩石中的裂纹扩展过程，求得了损伤柔度的表达式，并给出了相应的算例，数值模拟结果与Ｔａｙｌｏｒ模型和自洽模型的计算结果以及实验结果进行了对比，结果表明Ｔａｙｌｏｒ模型和自洽模型的数值模拟结果与实验结果有较大的误差，而本文提出的数值模拟结果则与实验结果具有很好的一致性。     

参数型岩石几何损伤模型的构建及其应用研究

岩石几何损伤模型是建立统计损伤本构模型的重要基础。在岩石变形力学特性基础上对现有岩石几何损伤模型进行了回顾性分析,针对它们难以较好地反映初始损伤特征和峰后变形破坏特征的缺陷与不足,首先,将岩石材料视为由未损伤部分、初始损伤部分和后继损伤部分组成,提出了考虑初始损伤的岩石几何损伤模型;然后,通过探讨Weibull分布参数m和F₀对损伤变量变化产生的影响,构建了参数型岩石几何损伤模型,在此基础上建立了应变软化类岩石统计损伤本构模型并进行了修正,给出了模型参数的确定方法;最后,通过模型验证和参数分析表明,修正后模型能够较好地模拟岩石变形破裂全过程,参数λ和η对损伤变量变化产生的影响与m和F₀等效,解决了现有岩石几何损伤模型存在的共性问题,该模型和方法具有一定的合理性和可行性。     

岩石统计渗流模型和统计损伤本构模型研究

采用Eshelby等效夹杂方法建立考虑渗流的岩石损伤本构模型是一种有效方法,但相关文献目前还极少。利用细观力学的Eshelby等效夹杂方法, 探索了考虑渗流和损伤的Helmholtz自由比能函数的确定,用连续介质损伤力学方法建立了相应的考虑渗流的岩石损伤统计本构模型,提出了考虑渗流和损伤过程的岩石破坏准则,建议了考虑损伤和应变的岩石渗透系数演化方程。与试验结果比较表明,所建立模型是合理的。     

基于影响因子修正的岩石损伤力学模型探讨

针对现有岩石损伤模型的局限性与不足,在深入探讨岩石变形损伤机制的基础上考虑岩石损伤的过程,提出损伤部分可分为损伤材料和因材料损伤而变化的微缺陷。在该基础上提出新的损伤定义,还提出了损伤材料与因材料损伤而变化的微缺陷之间的影响因子?,同时提出其参数确定方法,研究其基本性质。通过岩石微观受力分析,建立了岩石损伤模型。通过深入研究岩石有效应力与名义应力的关系,并引进统计损伤理论,建立岩石损伤统计本构模型。该模型能够充分地反映岩石的应变软化特性和应变硬化特性,且参数物理意义明确,应用将较为方便。最后,与实测结果及前人成果的比较分析,验证了该模型的合理性、可行性。     

岩石爆破损伤模型的比选与改进

归纳了目前应用较为广泛的爆破损伤模型关于损伤变量的定义方法,基于FORTRAN与LS-DYNA自定义接口,将5种典型的爆破损伤模型成功导入LS-DYNA。根据具体的工程实例,对5种爆破损伤模型计算的精确性进行了对比计算和验证。结果表明,KUS模型与RFPA模型对应的结果与实测值更接近。选取计算结果较为准确的KUS模型进行改进,考虑压缩损伤,并修正了宏观弹性常数的确定方法,建立了拉压损伤模型,并对该模型的计算精确性进行验证,结果表明,该拉压损伤模型可以更好地定量描述岩体爆破损伤范围。     

基于岩石统一能量屈服准则的硬岩损伤模型

为了更准确、便捷地确定深埋地下工程中硬岩的力学行为,开展了如下研究：根据损伤与塑性演化存在相似机制,假设损伤阈值与初始屈服采用相同的确定方法,将岩石统一能量屈服准则作为损伤阈值的判定准则;基于Rabotnov对损伤变量的定义和Lemaitre应变等效性假设,并考虑未损伤部分的应力-应变关系符合广义虎克定律,提出了简化形式的硬岩损伤本构方程;为了更好地表达硬岩从低围压到高围压条件下的脆-延转化特性,提出了改进的Mazars损伤演化方程;基于锦屏T2b大理岩常规三轴压缩试验结果,对损伤模型参数进行求解,并分别与基于Mazars损伤演化方程和Mohr-Coulomb准则为损伤阈值判定准则的损伤模型进行对比,结果表明：提出的损伤模型可更好地表达硬岩从低围压的脆性到高围压的延性转变过程,特别是对损伤阈值后岩石的损伤演化过程的模拟更为精确,对硬脆性岩石的工程计算有一定的借鉴意义。     

基于新型损伤定义的岩石损伤统计本构模型探讨

针对几何损伤理论中损伤定义的局限性与不足,对岩石损伤定义进行了改进;在此基础上,将在应力作用下的岩石材料抽象为破坏与未破坏两部分,并根据这两部分不同的受力情况,建立出新型损伤模型;其次,利用统计损伤理论,建立出岩石统计损伤演化方程,进而建立起基于特定围压下的岩石损伤统计本构模型;然后,通过探讨模型参数与围压的经验关系,对模型进行合理修正,从而建立了反映不同围压条件下的岩石损伤统计本构模型,并通过探讨模型参数的物理意义,初步建立了反映岩石软化与硬化特性相互转化的临界物理力学参数的确定方法。该模型能够反映岩石软化与硬化特性相互转化的特征,与试验结果进行比较分析表明,该模型与实测结果吻合良好。     

Investigation and modeling on fatigue damage evolution of rock as a function of logarithmic cycle

The fatigue damage behavior of granite under constant and variable amplitude loadings is studied. The experimental analysis reveals that there is a three‐stage law for the fatigue damage evolution as a function of absolute or relative cycle and the inverted‐S damage model proposed by the author, in this case, is capable of representing the damage behavior of rock. However, when the logarithmic cycle is considered, there are only two stages, i.e. steady and accelerated stages and the fatigue damage evolution greatly depends on the properties of rock and stress level. Accordingly, the fatigue damage evolution curves have been categorized into three types. Then, the effect of maximum stress, amplitude and fatigue initial damage on the damage evolution of rock is investigated. The analysis reveals that the damage evolution greatly depends on these influencing factors. The fatigue life decreases with the increase in the maximum stress, amplitude and fatigue initial damage due to the decrease in the proportion of the first stage to the whole fatigue process and the increase in the damage rate in the first stage. Meanwhile, a linear‐exponential formula is used to model the fatigue damage behavior of rock subjected to cyclic loading. This damage model is superior to the inverted‐S damage model in the convenience of establishment of critical instability point. The physical meanings of its constants have been illuminated and the applicability of this model to constant and variable amplitude cyclic loading explored. The fitting results for the test data show that this damage model can properly represent the fatigue damage behavior of rock. Copyright © 2010 John Wiley & Sons, Ltd.     

周期荷载下盐岩的疲劳变形及损伤特性研究

利用RMT－150C型岩石力学多功能试验机,进行了盐岩单轴循环荷载作用下的疲劳试验,研究了盐岩的疲劳强度、变形及损伤特性。试验结果表明,当上限应力大于“门槛值”时,盐岩疲劳破坏时的轴向应变可以分为初始变形、等速变形和加速变形3个阶段,呈疏－密－疏的发展过程。改变上限应力和平均应力会显著影响疲劳的进程,提高上限应力值和平均应力值,初始轴向变形和循环轴向变形的比率都会提高,疲劳破坏时的总循环次数显著减小。盐岩疲劳破坏终点的变形量同样受静态轴向应力-应变全过程曲线的控制,控制误差范围在10%左右。循环荷载作用下,盐岩的变形模量经历了一个先逐渐增加,后缓慢降低,最后加速减小的过程。对循环荷载作用下盐岩疲劳损伤演化规律进行了初步探讨     

巷道软岩蠕变损伤模型

介绍了西原及修正西原模型的蠕变特性、物理机制和优点。根据软岩加速蠕变的力学特性,指出西原模型不能较好的描述蠕变加速阶段,并引入损伤变量,认为蠕变引起的应变是由无损材料及已损材料蠕变应变之和,推导出材料非稳定蠕变阶段的蠕变方程。同时,利用有限元方法,建立了带有损伤变量蠕变单元方程,可以此编制有限元程序,进行蠕变分析。     

考虑初始损伤和蠕变损伤的岩石蠕变全过程本构模型

考虑到天然岩石存在不同程度初始损伤以及蠕变过程中岩石受载后裂隙扩展而导致的新损伤,对具有初始损伤的岩石蠕变特性进行全面描述。根据不闭合结构面应力与法向变形之间的关系,提出裂隙岩石塑性变形体元件,描述岩石蠕变过程中的瞬时塑性变形。引入初始损伤影响因子,建立具有初始损伤的岩石损伤变量演化方程,构建模拟岩石加速蠕变的蠕变损伤体元件。将裂隙岩石塑性变形体和蠕变损伤体与描述瞬时弹性变形和黏弹性变形的广义开尔文模型进行串联组合,形成能够反映具有初始损伤的岩石瞬时弹-塑性变形、稳定蠕变和加速蠕变的蠕变全过程本构模型,提出了进行少量蠕变试验既能解析模型参数的方法,在不同应力水平下模型理论曲线与蠕变试验曲线吻合。     

综合考虑宏微观复合损伤的节理岩体本构模型

针对目前节理岩体损伤本构模型中仅考虑节理等宏观缺陷造成的损伤,而没有考虑岩块中微裂纹等微观损伤的不足,提出了综合考虑宏微观损伤的节理岩体本构模型。其中微观损伤模型采用基于应变强度理论和岩石微元强度服从Weibull分布假定的统计损伤模型,把其应用于被节理切割而成的岩块。而宏观损伤模型主要基于连续介质力学原理,把节理对岩体性质的劣化作为损伤来考虑。在考虑节理传压及传剪系数的基础上,提出了综合考虑宏微观复合损伤的节理岩体本构模型。算例表明宏观节理的存在大大削弱了岩体的强度,降低了其刚度,增大了其柔性。所提出的模型能够较好地反映宏微观两类损伤对岩体应力-应变关系及强度的影响,较为合理。     

岩石损伤过程中的热-流-力耦合模型及其应用初探

岩石损伤过程热-流-力（THM）耦合问题的研究对于深部采矿等许多工程领域都具有重要的理论意义。以岩石的损伤为主线,在多场耦合分析方程中引入损伤变量,基于质量守恒和能量守恒原理,提出岩体损伤过程中的THM耦合模型。通过把均匀弹性介质THM耦合响应的模拟结果与理论分析结果进行对比,验证了程序及有限元实施的正确性。然后,用该耦合模型进行了不同地应力条件下流固耦合过程的数值模拟,探讨了水压力对于岩石损伤过程的作用机制。数值模拟表明,水压力导致了拉伸损伤范围的扩大和损伤程度的加剧,同时亦对剪切损伤具有抑制作用。     

Microplane damage model for jointed rock masses

The paper presents a new microplane constitutive model for the inelastic behavior of jointed rock masses that takes into account the mechanical behavior and geometric characteristics of cracks and joints. The basic idea is that the microplane modeling of rock masses under general triaxial loading, including compression, requires the isotropic rock matrix and the joints to be considered as two distinct phases coupled in parallel. A joint continuity factor is defined as a microplane damage variable to represent the stress‐carrying area fraction of the joint phase. Based on the assumption of parallel coupling between the rock joint and the rock matrix, the overall mechanical behavior of the rock is characterized by microplane constitutive laws for the rock matrix and for the rock joints, along with an evolution law for the microplane joint continuity factor. The inelastic response of the rock matrix and the rock joints is controlled on the microplane level by the stress–strain boundaries. Based on the arguments enunciated in developing the new microplane model M7 for concrete, the previously used volumetric–deviatoric splits of the elastic strains and of the tensile boundary are avoided. The boundaries are tensile normal, compressive normal, and shear. The numerical simulations demonstrate satisfactory fits of published triaxial test data on sandstone and on jointed plaster mortar, including quintessential features such as the strain softening and dilatancy under low confining pressure, as well as the brittle–ductile transition under higher confining pressure, and the decrease of jointed rock strength and Young's modulus with an increasing dip angle of the joint. Copyright © 2014 John Wiley & Sons, Ltd.