考虑微裂纹相互作用的的岩石微观力学弹塑性损伤模型研究

本文建立基于微裂纹扩展的岩石弹塑性损伤微观力学模型。用自洽方法考虑裂隙间相互影响,压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的微观损伤,基于应变能密度准则用Newton迭代法求复合型断裂的翼裂纹扩展长度,并采用微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映塑性内变量对硬化函数的影响;建立岩石的弹塑性损伤本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性损伤模型的程序。从围压和微裂隙长度等因素分析弹塑性损伤模型的岩石的损伤和宏观塑性特性。     

基于裂纹扩展脆性岩石动力压缩力学特性研究

动态压缩荷载作用下,脆性岩石内部动态细观裂纹扩展特性,对岩石宏观动态力学特性有着重要的影响。然而,对岩石内部动态细观裂纹扩展与宏观动态力学特性的关系研究较少。基于准静态裂纹扩展作用下的应力-应变本构模型、准静态与动态裂纹扩展断裂韧度关系、裂纹速率与应变率关系模型及应变率与动态断裂韧度关系,提出了一种基于细观力学的动态应力-应变本构模型。其中裂纹速率与应变率关系,是根据裂纹长度与应变关系的时间导数推出;应变率与动态断裂韧度关系,是根据推出的裂纹速率及应变率关系,与裂纹速率及断裂韧度关系相结合而得到。研究了应变率对应力-应变本构关系及动态压缩强度影响。并通过试验结果验证了模型的合理性。讨论了岩石初始损伤、围压、模型中参数m、ε0和R对应力-应变关系、动态压缩强度和动态弹性模量的影响。研究结果可为动态压缩荷载作用下深部地下工程脆性围岩稳定性分析提供了一定的理论支持。     

岩石化学腐蚀-围压细观损伤本构模型与离散元模拟研究

基于Lemaitre应变等价性假设理论,假定受水化学-力耦合损伤的岩石微元强度服从Weibull分布,考虑化学腐蚀与围压耦合作用对岩石力学参数的影响,通过核磁共振技术与损伤力学理论,引入细观化学损伤变量与力损伤变量,并认为微元破坏符合SMP准则,建立岩石化学腐蚀-力耦合损伤本构模型,并采用理论推导的方法得出所需的模型参数。同时基于颗粒离散元方法,引入参数半径乘子来改变颗粒间的黏结接触尺寸,从而模拟水化学损伤,采用平直节理模型对水化学作用后的岩石进行三轴压缩模拟,得到了水化学作用和不同围压下的岩石三轴应力-应变模拟曲线。通过对比所构建的岩石化学腐蚀-力耦合损伤本构模型理论曲线、离散元模拟曲线和试验曲线,结果表明三者吻合度较好,能够很好地反映岩石在化学腐蚀和围压耦合作用下的力学特性与破坏特征,并通过离散元方法得到了岩石在三轴压缩过程中裂纹的产生与分布情况。     

脆性岩石蠕变裂纹成核宏细观力学机理研究

脆性岩石内部细观裂纹扩展、贯通及成核影响下的脆性蠕变行为, 对深部地下工程围岩微震及岩爆事件评价有着重要意义. 然而, 目前能够解释裂纹成核损伤突变影响下, 脆性岩石完整蠕变宏细观力学机理模型的研究很少. 本文基于脆性岩石亚临界裂纹扩展模型、裂纹-应变-声发射事件相关的损伤模型及裂纹成核损伤时间演化路径函数, 提出了一种脆性岩石裂纹成核损伤突变影响下的蠕变宏细观力学模型. 裂纹成核损伤时间演化路径函数通过岩石内部裂纹成核损伤突变大小$\Delta D_{CN}$及相邻裂纹成核损伤时间差$\Delta t$进行定义, 该函数可以结合岩石声发射监测试验数据定义的岩石损伤数据确定. 通过与试验结果对比分析验证模型的合理性. 并讨论了裂纹成核损伤大小、相邻裂纹成核损伤时间间隔、及裂纹成核数量对脆性岩蠕变裂纹长度、裂纹速率、轴向应变及应变率的影响. 该结果对于更加合理、经济、高效的深部地下工程施工及设计提供了一定的理论支持.      

基于细观力学的混凝土弹塑脆性损伤研究

将混凝土假定为一种由硬化水泥砂浆、粗骨料、界面粘结带所组成的三相复合材料,在满足骨料级配曲线算法的基础上,采用细观单元的弹塑脆性损伤本构关系,考虑材料的非均质特性,建立了基于细观力学的混凝土弹塑脆性损伤数值模型;分别研究了单轴受拉预置裂纹试样和单轴受压混凝土试样的细观弹塑脆性损伤破坏行为,并揭示了混凝土的宏观表征强度存在明显的尺寸效应,通过将计算结果与 Bazant 尺寸效应公式、单轴受压物理实验曲线进行对比,验证了模型的正确性。数值试验表明：该模型可以清晰地模拟混凝土细观塑性屈服和失效裂纹的萌生和扩展。骨料与水泥砂浆间的界面粘结带相对薄弱,在混凝土试件形成宏观损伤局部化带前,试件的屈服和破坏首先发生在骨料边缘处的界面位置,并沿着界面粘结带扩展、贯通;同时,导致宏观裂纹形成和发展的因素仍以细观单元的拉伸破坏为主。     

综合考虑宏细观缺陷的岩体动态损伤本构模型

针对节理岩体同时含有节理、裂隙等宏观缺陷及微裂隙、微孔洞等细观缺陷的客观事实, 提出了在节理岩体动态损伤本构模型中应同时考虑宏细观缺陷的观点。为此, 首先对基于细观动态断裂机理的经典岩石动态损伤本构模型—TCK(Taylor-Chen-Kuszmaul)模型进行了阐述, 其次基于Lemaitre等效应变假设推导了综合考虑宏细观缺陷的复合损伤变量(张量), 进而在此基础上建立了相应的节理岩体动态损伤本构模型, 并利用该模型讨论了载荷应变率及节理条数对岩体动态力学特性的影响规律。结果表明, 在不同载荷应变率下试件在变形初始阶段是重合的, 而后随着应变的增加, 试件峰值强度、峰值应变及总应变均随载荷应变率的增加而增加; 随着节理条数的增加, 试件峰值强度逐渐降低, 但降低趋势逐渐变缓并趋于某一定值。上述研究结论与目前的理论及实验研究结果的基本规律是一致的, 说明了本模型的合理性。     

非贯通裂隙岩体损伤演化率相关性及变形特征

含非贯通裂隙岩体是自然界中岩体的主要赋存形式,其裂隙几何特征对岩体的强度及变形均产生显著影响。应变率对岩体的损伤演化及黏滞效应也具有显著的率相关性。首先,运用模型元件的方法,将非贯通裂隙岩体动态破坏过程视为具复合损伤、静态弹性特性、动态黏滞特性的非均质点组成,对黏弹性响应的Maxwell体进行改进,将细观损伤体与裂隙损伤演化的宏观损伤体根据等效应变假设并联组成宏细观复合损伤体,构建综合考虑岩体宏细观缺陷的动态损伤模型;其次,基于断裂力学及应变能理论,对岩体宏观裂隙动态扩展的能量机制进行分析,综合考虑初始裂隙应变能、裂隙动态损伤演化过程应变能、裂隙闭合应变能,得到裂隙岩体宏观动态损伤变量计算公式;最后,将模型计算结果与实验结果进行比较,模型计算结果与实验结果吻合较好,证明了模型的合理性,同时利用模型讨论了裂隙倾角、应变率、岩石性质对岩体变形特征的影响规律。     

基于细观力学的脆性岩石蠕变损伤特性研究

由于孔隙水长期应力侵蚀作用下的亚临界裂纹扩展对脆性岩石的蠕变特性有着重要的影响.论文基于翼型裂纹细观模型与裂纹扩展法则,得到了常压应力作用岩石的裂纹长度演化规律,然后结合宏观应变关系,分析了岩石的蠕变应变演化规律.基于上述关系,讨论了初始损伤一定情况下,裂纹角度对于岩石的蠕变失效特性及短期强度的影响,结果表明了裂纹角度对于脆性岩石的力学特性有着重要的影响.然后结合花岗岩细观力学试验参数,分析了花岗岩的裂纹尖端强度因子,裂纹长度,损伤,蠕变演化规律.并对比分析了蠕变应变的理论与试验结果,两者具有一定的相似性,验证了理论模型的合理性,进而对脆性岩石工程的设计与施工提供了一定的帮助.     

具有体膨胀和剪切效应的结构陶瓷相变塑性细观本构模型:Ⅰ.非比例加载历史

本文在对结构陶瓷的四方至单斜(t→m)马氏体相变进行细观力学、热力学和微观机制分析的基础上,导出了在非比例加载条件下考虑材料的体膨胀和剪切效应的相变塑性细观本构模型。作者首次采用 Mori-Tanaka 方法以自洽的方式导出了材料构元的 Helmho-ltz 自由能及余能函数的解析表达式,它是外加宏观应力(或应变)、温度、相变夹杂体积分数以及夹杂内平均相变应变的函数,其中夹杂体积分数和平均相变应变为描述材料构元微结构变化的内变量。最后按 Hill-Rice 本构理论框架导出相变塑性屈服面方程及增量本构关系。     

水泥基材料在宽围压范围的本构模型

本文首先进行了水泥石的三轴试验和静水试验,根据试验结果分析了宽围压范围下不同力学机理对水泥石宏观力学性能的影响。提出了基于热力学定律的剪切塑性、孔隙塑性与破坏准则相互作用的弹塑性损伤本构模型。该模型有效的考虑了围压对材料主要力学特征的影响,提出了与损伤准则相关联的塑性硬化函数。数值解答与试验结果比较表明,该本构模型可以很好地描述水泥石在低围压状态下的脆性、中围压下的脆-塑性转换以及高围压下的孔隙塑性为主的力学特征。     

单轴拉伸条件下脆性岩石微裂纹损伤模型研究

利用断裂力学、损伤力学和均匀化原理,对脆性岩石单轴拉伸条件下的力学特性进行分析,建立了脆性岩石的微裂纹损伤本构模型.首先对岩石内部微裂纹的统计分布规律进行分析,给出了理论分析过程中微裂纹分布的假设条件,在此基础上,参考已有研究成果,得到含细长微裂纹脆性岩石有效弹性参数的计算公式.然后,对岩石内部单一微裂纹进行断裂力学和损伤力学分析,得到了扩展裂纹尖端的应力强度因子计算公式,在一定微裂纹断裂扩展准则和断裂扩展速率的假设基础上,利用积分原理,得到了岩石整体的损伤变量和损伤演化方程,由此建立单轴拉伸条件下脆性岩石的微裂纹损伤本构模型.最后,通过一花岗岩的单轴拉伸试验结果对微裂纹损伤本构模型进行了验证.     

Asymptotic homogenization analysis for damage amplification due to singular interaction of micro-cracks

The paper investigates the overall damage amplification effect due to micro-crack interaction in a framework of two-scale modeling. A homogenization method based on asymptotic expansions is employed to deduce the macroscopic damage equations. The damage model completely results from energy-based micro-crack propagation laws. We consider a locally periodic microstructure with periods containing pairs of micro-cracks separated by small ligaments. The asymptotic solution in the ligament region allows the study of the effect of micro-crack interaction on the effective coefficients. The local macroscopic response expresses the collective coalescence of a periodic microstructure with interacting micro-cracks. We show that the slope of the homogenized coefficients is inversely proportional to the square root of the distance between the tips of the interacting micro-cracks, accounting for the singularity in the stress fields as the micro-cracks approach each other. This leads to damage amplification as the result of the interaction of micro-cracks.     

An interacting micro-crack damage model for failure of brittle materials under compression

A model is developed for brittle failure under compressive loading with an explicit accounting of micro-crack interactions. The model incorporates a pre-existing flaw distribution in the material. The macroscopic inelastic deformation is assumed to be due to the nucleation and growth of tensile “wing” micro-cracks associated with frictional sliding on these flaws. Interactions among the cracks are modeled by means of a crack-matrix-effective-medium approach in which each crack experiences a stress field different from that acting on isolated cracks. This yields an effective stress intensity factor at the crack tips which is utilized in the formulation of the crack growth dynamics. Load-induced damage in the material is defined in terms of a scalar crack density parameter, the evolution of which is a function of the existing flaw distribution and the crack growth dynamics. This methodology is applied for the case of uniaxial compression under constant strain rate loading. The model provides a natural prediction of a peak stress (defined as the compressive strength of the material) and also of a transition strain rate, beyond which the compressive strength increases dramatically with the imposed strain rate. The influences of the crack growth dynamics, the initial flaw distribution, and the imposed strain rate on the constitutive response and the damage evolution are studied. It is shown that different characteristics of the flaw distribution are dominant at different imposed strain rates: at low rates the spread of the distribution is critical, while at high strain rates the total flaw density is critical.     

A two-scale damage model with material length

The Note presents the formulation of a class of two-scale damage models involving a micro-structural length. A homogenization method based on asymptotic developments is employed to deduce the macroscopic damage equations. The damage model completely results from energy-based micro-crack propagation laws, without supplementary phenomenological assumptions.We show that the resulting two-scale model has the property of capturing micro-structural lengths. When damage evolves, the micro-structural length is given by the ratio of the surface density of energy dissipated during the micro-crack growth and the macroscopic damage energy release rate per unit volume of the material.The use of fracture criteria based on resistance curves or power laws for sub-critical growth of micro-cracks leads to quasi-brittle and, respectively, time-dependent damage models. To cite this article: C. Dascalu, C. R. Mecanique 337 (2009).     

A penny-shaped cohesive crack model for material damage

A novel micromechanics based damage model is proposed to address failure mechanism of defected solids with randomly distributed penny-shaped cohesive micro-cracks (Barenblatt–Dugdale type). Energy release contribution to the material damage process is estimated in a representative volume element (RVE) under macro hydrostatic stress state. Macro-constitutive relations of RVE are derived via self-consistent homogenization scheme, and they are characterized by effective nonlinear elastic properties and a class of pressure sensitive plasticity which depends on crack opening volume fraction and Poisson’s ratio. Several distinguished features of the present model are compared with Gurson model and Gurson–Tvergaard–Needleman (GTN) model, showing that the proposed model can better capture material degradation and catastrophic failure due to cohesive micro-crack growth and coalescence.     

Micro-crack tip fracture of commercial grade aluminum under mixed mode loading

In situ tensile tests were made in a scanning electron microscope (SEM) to investigate the deformation and micro-fracture in the immediate vicinity of a micro-crack tip in commercial pure aluminum with large-size crystal. Examined are the slip line field, stress intensity factor, strain energy density factor and crack tip opening displacement (CTOD) for mixed mode loading. Blunting and sharpening effects are observed. The latter is controlled by localized slip while the former by uniformed slip of the operating slip system with the highest crack tip Schmid factor. The operating slip system depends on the crystallographic orientation of crystal containing micro-cracks.The damage and fracture take place in the blunted region and depend on the coarsening and spacing of uniformed slip lines. The mixed mode micro-crack propagates along the direction where the voids grow and coalesce into the micro-crack. The direction also depends on the orientation of the applied loading. This suggests that the formation of macro-fracture mechanics could be applied. In particular, the minimum strain energy density criterion is suitable for determining the direction of micro-crack instability in the mixed mode. The in situ data were used to yield a nearly constant critical, minimum strain energy density factor for onset of micro-cracking.     

Effective elastic moduli of two dimensional solids with distributed cohesive microcracks

Effective elastic properties of a defected solid with distributed cohesive micro-cracks are estimated based on homogenization of the Dugdale–Bilby–Cottrell–Swinden (Dugdale–BCS) type micro-cracks in a two dimensional elastic representative volume element (RVE).Since the cohesive micro-crack model mimics various realistic bond forces at micro-scale, a statistical average of cohesive defects can effectively represent the overall properties of the material due to bond breaking or crack surface separation in small scale. The newly proposed model is distinctive in the fact that the resulting effective moduli are found to be pressure sensitive.     

基于数字散斑相关实验测量的材料构型力的计算方法

材料构型力学主要研究材料中的缺陷(夹杂、空穴、位错、裂纹、塑性区等)的构型(形状、尺寸和位置)改变时,所引起的系统自由能的变化。本研究将基于数字散斑相关技术,实验测量材料试件的位移场分布,随后通过材料构型力的定义式,计算求得弹塑性材料中缺陷构型力的分布。其方法概括如下:位移场通过数字图像相关技术测得;应变及位移梯度场利用三次样条拟合获得;线弹性材料应力通过简单线弹性本构方程获取,而塑性材料的表面应力场通过Ramberg-Osgood本构方程计算求得;弹塑性应变能密度分布则由应力-应变曲线数值积分获得。该方法对普通弹性材料或者弹塑性材料均适用,可以用于各种不同的缺陷及缺陷群的材料构型力测量。     

Stochastic damage model for concrete based on energy equivalent strain

Starting with the framework of conventional elastoplastic damage mechanics, a class of stochastic damage constitutive model is derived based on the concept of energy equivalent strain. The stochastic damage model derived from the parallel element model is adopted to develop the uniaxial damage evolution function. Based on the expressions of damage energy release rates (DERRs) conjugated to the damage variables thermodynamically, the concept and its tensor formulations of energy equivalent strain is proposed to bridge the gap between the uniaxial and the multiaxial constitutive models. Furthermore, a simplified coupling model is proposed to consider the evolution of plastic strain. And the analytical expressions of the constitutive model in 2-D are established from the abstract tensor expression. Several numerical simulations are presented against the biaxial loading test results of concrete, demonstrating that the proposed models can reflect the salient features for concrete under uniaxial and biaxial loading conditions.