A crack-mechanics based model for damage and plasticity of brittle materials under dynamic loading

A rate-dependent model for damage and plastic deformation of brittle materials under dynamic loading is presented. The model improves upon a recently developed micromechanical damage model (Zuo et al., 2006) by incorporating plastic deformation of the material. The distribution of the microcracks in the material is assumed to remain isotropic, and the damage evolution is through the growth of the average crack size. Plasticity is considered through an additive decomposition of the total strain rate, and a rate-independent, von Mises model is used. The model was applied to simulate the response of a model material (SiC) under uniaxial strain loading. To further examine the behavior of the model, cyclic loading and large-strain compressive loading were considered. Numerical results of the model predictions are presented, and comparisons with those from a previous model are provided.     

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.     

脆性材料在双向应力下的断裂实验与理论分析

研究了脆性材料在双向应力下的断裂特性和失效机理，特别是在平行于裂纹的应力对临界断裂参数的影响方面进行了实验上和理论上的研究．采用玻璃、陶瓷等脆性材料进行了平面双向拉伸和单向拉伸试验，并对实验结果进行比较．观测直通裂纹的启裂和扩展过程，证明了双向应力对裂纹驱动力有明显影响，讨论了裂纹扩展的应变准则．     

砂岩脆性变形破坏过程中声发射信息试验研究

针对砂岩试样进行了单轴刚性加载试验,并测试整个过程的声发射信息。根据试验情况,可将砂岩在整个变形阶段的声发射信号特征划分为压密、弹性、稳定破裂、非稳定破裂和峰后情况五个阶段。通过对稳定破裂和非稳定破裂分界点对应的临界点处的应力和应变值与峰值应力及其对应的应变比值分析发现,大部分应变比在74%与78%之间;应力比则相对具有较大的离散性,均值在73%左右。上述研究,通过对砂岩脆性破坏过程中声发射信息与应力比和应变比之间的关系分析,增进了对岩石破裂过程中生发现象的认识,为相应岩石脆性破坏导致的地质灾害分析提供了新的研究思路。     

Modelling of crack growth by fracture damage zone

The strain energy density (SED) criterion is applied for analyzing the full range of mixed mode fracture from tensile to shear loading. A fracture damage zone (FDZ) local to the crack tip is defined and discussed in connection with the influence of crack geometry, loading and local material property. The size of FDZ tends to change continuously from statically to cyclically applied load conditions. It can be estimated from the uniaxial mechanical properties of the material. Both experimental and analytical results are examined for subcritical crack growth under static loading that depends on the type steel structures the fracture behavior of which could be represented by a single curve for the given specimen geometry.     

Discrete element modeling on the crack evolution behavior of brittle sandstone containing three fissures under uniaxial compression

Based on experimental results of brittle, intact sandstone under uniaxial compression, the micro-parameters were firstly confirmed by adopting particle flow code(PFC2D). Then, the validation of the simulated models were cross checked with the experimental results of brittle sandstone containing three parallel fissures under uniaxial compression. The simulated results agreed very well with the experimental results, including the peak strength, peak axial strain, and ultimate failure mode. Using the same microparameters, the numerical models containing a new geometry of three fissures are constructed to investigate the fissure angle on the fracture mechanical behavior of brittle sandstone under uniaxial compression. The strength and deformation parameters of brittle sandstone containing new three fissures are dependent to the fissure angle. With the increase of the fissure angle, the elastic modulus, the crack damage threshold,and the peak strength of brittle sandstone containing three fissures firstly increase and secondly decrease. But the peak axial strain is nonlinearly related to the fissure angle. In the entire process of deformation, the crack initiation and propagation behavior of brittle sandstone containing three fissures under uniaxial compression are investigated with respect tothe fissure angle. Six different crack coalescence modes are identified for brittle sandstone containing three fissures under uniaxial compression. The influence of the fissure angle on the length of crack propagation and crack coalescence stress is evaluated. These investigated conclusions are very important for ensuring the stability and safety of rock engineering with intermittent structures.     

Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones, notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode II component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.     

Plastic strip model of circular crack under cyclic loading with damage accumulation

A cylindrical solid containing a penny-shaped crack in its mid plane is subjected to a remote tensile stress field. In the plastic region near the crack border, damage due to continuous deterioration of the material accumulates, and may lead to crack instability and crack growth. A damage model relating the crack opening displacement in the plastic zone to the fraction of the cross-sectional area occupied by voids is used to describe the conditions necessary for the onset of crack instability, fatigue crack propagation by cyclic loading, and rates of crack growth for time-dependent environmental effects.     

An improved atomistic simulation based mixed-mode cohesive zone law considering non-planar crack growth

A novel and improved atomistic simulation based cohesive zone law characterizing interfacial debonding is developed which explicitly accounts for the non-planarity of the crack propagation. Group of atoms in the simulation constituting cohesive zones which are used to obtain local stress and crack opening displacement data are determined dynamically during the non-planar crack growth as they cannot be determined apriori. The methodology is used to study the debonding of Σ5 (2 1 0)/[0 0 1] symmetric tilt grain boundary interface in a Cu bicrystal under several mixed mode loading conditions. Simulations show that such bicrystalline specimen exhibits three types of energy dissipative mechanisms – shear coupled GB migration (SCM) away from the crack-tips, change in spacial orientation of GB structural units rendering highly disordered grain boundary near the crack tips and brittle intergranular fracture. Which combination of these three deformation mechanism will be active influencing the degree of non-planarity of the crack propagation at various stages of loading depends on the loading mode-mixity. As the ratio of shear component of the loading parallel to the GB plane and normal to the tilt axis with respect to the normal loading increases (thereby increasing the mode-mixity), overall strain-to-failure also increases and SCM tends to become the dominant deformation mechanism. Through this framework, analytical functional forms and parameters describing cohesive laws for both normal and shear traction as a function of the mode-mixity of the loading and crack opening displacement are predicted.     

Damage accumulation for growth of anti-plane crack in work-hardening material

Elastic–plastic solutions of an anti-plane crack in an infinite body are used in conjunction with a continuum damage model to describe the conditions necessary for the onset of crack instability, fatigue crack propagation due to cyclic loading, and rates of crack growth due to time dependent events. A power law relates the stress to the strain of the material. The damage, which invokes nucleation, growth and coalescence of microvoids due to elevated strain, is confined to the plastic zone surrounding the crack tip. For applied loading below the yield stress, the small-scale and large-scale yielding solutions are used to determine the influence of strain hardening on crack instability and failure. Crack growth due to cyclic loading and time-dependent deformations are studied using the small-scale yielding solution of the deformation theory of plasticity.     

A micro–macro method for predicting the shear strength of brittle rock under compressive loading

Microcracks have great significance for shear strength of brittle rock in compression. A major challenge of this area is to establish the correlation of microcracks and macroscopic shear strength. A new micro–macro method is presented to predict the shear strength of brittle rock in compression. This method incorporates the microcrack model suggested by Ashby, Mohr–Coulomb failure criterion and a crack-strain relation. This crack–strain relation is presented to link the crack growth and axial strain by combining the micro and macro definitions from rock damage. The shear strength and stress–strain relationship of Jinping marble are theoretically investigated in detail. The rationality of this suggested method is verified by using the experimental results founded on Jinping marble. Effects of the initial microcrack size, friction coefficient and confining pressure on internal friction angle, cohesion, and shear strength are also discussed.     

单轴压缩下页岩类脆性材料细观损伤演化的离散元分析

了解广泛存在的类似页岩的脆性材料各向异性对工程安全具有重要意义。本研究将页岩视为粘结颗粒材料,基于离散单元方法研究了横观各向同性脆性页岩的损伤演化。再现了不同层理角的页岩试样的破坏模式,并对比了实验和数值模拟的抗压强度和弹性模量。引入微裂纹的概念,通过定义裂纹密度函数,系统地研究了单轴压缩条件下,页岩层理角对细观结构的影响。此外基于平均配位数建立了配位数变化与细观损伤的联系,并根据配位数的变化与裂纹数量将加载过程分为三个阶段,分析了不同阶段配位数与裂纹数量的对应变化关系。研究表明,页岩的裂纹密度随着层理角的增加而增加,而试样的平均配位数在加载过程中先上升后剧烈下降,颗粒集合体在单轴压缩条件下的应力应变及裂纹数量曲线与平均配位数曲线有良好的一致性。该研究揭示了横观各向同性脆性岩石的破坏过程和内在机理,将为页岩类脆性材料的工程应用提供理论指导。     

Constitutive modeling of ice in the high strain rate regime

The objective of the present work is to propose a constitutive model for ice by considering the influence of important parameters such as strain rate dependence and pressure sensitivity on the response of the material. In this regard, the constitutive model proposed by Carney et al. (2006) is considered as a starting basis and subsequently modified to incorporate the effect of brittle cracking within a continuum damage mechanics framework. The damage is taken to occur in the form of distributed cracking within the material during impact which is consistent with experimental observations. At the point of failure, the material is assumed to be fluid-like with deviatoric stress almost dropping down to zero. The constitutive model is implemented in a general purpose finite element code using an explicit formulation. Several single element tests under uniaxial tension and compression, as well as biaxial loading are conducted in order to understand the performance of the model. Few large size simulations are also performed to understand the capability of the model to predict brittle damage evolution in un-notched and notched three point bend specimens. The proposed model predicts lower strength under tensile loading as compared to compressive loading which is in tune with experimental observations. Further the model also asserts the strain rate dependency of the strength behavior under both compressive as well as tensile loading, which also corroborates well with experimental results.     

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

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

An augmented multicrack elastoplastic damage model for tensile cracking

In this paper we present the general formulation and numerical aspects of an augmented multicrack elastoplastic damage model aiming to reflect the crack induced anisotropy in concrete like quasi-brittle materials. Consistent evolution laws for the involved internal variables are derived based on the augmented Lagrangian method. The (time) discrete formulation and the corresponding variational structure are investigated, with the Euler–Lagrangian equations defining the closest-point projection approximation of the proposed model. The numerical aspects, such as the stress updating algorithm and the algorithmic consistent tangent moduli, are also discussed in details. It is found that in the developed numerical algorithm the active loading surfaces are determined in such a posterior manner that potential numerical problems due to the iteratively updating procedure in classical algorithms can be avoided. The proposed model is applied to the modeling of tensile cracking in concrete. The behavior of a single crack is characterized by an elliptical cracking surface and a hyperbolic softening function, with the orientations of potential cracks determined by Mohr’s postulate. The model is verified by calculating the single point stress vs. strain relations of concrete under several typical proportional and non-proportional loading cases. Finally, two benchmark tests of concrete structures, i.e. four-point bending beam under cyclic loading (Hordijk, 1992) and double edge notched specimens under mixed tension/shear forces (Nooru-Mohamed, 1992), are numerically simulated. Both predicted load vs. displacement curves and crack patterns agree well with the experimental data.     

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

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

Parameters controlling fracture resistance in functionally graded materials under mode I loading

In this investigation the fracture behavior of functionally graded materials (FGMs) was studied by means of experiments carried out on model polymer-based FGMs. Model graded materials were manufactured by selective ultraviolet irradiation of ECO [poly(ethylene carbon monoxide)], a photo-sensitive ductile copolymer that becomes more brittle and stiffer under exposure to ultraviolet light. The mechanical response of the graded material was characterized using uniaxial tensile tests. Single edge notched tension graded ECO specimens possessing different spatial variations of Young’s modulus, failure stress and failure strain were tested under remote opening loading. A full-field digital image correlation technique was used to measure in real-time the displacement field around the crack tip while it propagated through the graded material. The measured displacement field was then used to extract fracture parameters such as stress intensity factor and T-stress, and thus construct resistance curves for crack growth in the FGMs. For this loading configuration it was found that the nonsingular T-stress term in the asymptotic expansion for stresses needs to be accounted for in order to accurately measure the fracture resistance in FGMs. In addition, the influence of local failure properties (i.e., failure stress and failure strain) on crack growth resistance was investigated in detail. It was found that depending on the combined effects of the spatial variation of these two failure parameters, regardless of the spatial variation of the Young’s modulus, the FGM fracture resistance can either increase, decrease or remain constant with continued crack growth.     

Asymptotic character of mixed mode in plane deformation of crack in rubber-like material

The asymptotic stress and strain distribution near a crack tip in rubber-like materials is determined by finite element for in-plane mixed mode loading. For large strain, the crack tip field is always in a state of uniaxial tension. The shear load affects only the orientation of the deformed near tip field in the space. A good agreement is obtained between the theoretical and numerical results.     

PMMA材料的动态压缩力学特性及应变率相关本构模型研究

利用INSTRON万能试验机和分离式Hopkinson压杆(SHPB)对PMMA试件在较宽应变率范围内进行了单轴压缩实验,研究加载应变率对PMMA材料力学性能的影响.利用扫描电子显微镜对回收的试样进行了显微观察,重点分析不同加载应变率下PMMA的微观损伤破坏模式.结果表明:随着应变率的增大,PMMA的流动应力显著地增加,且冲击加载条件下,峰值应力的应变率敏感性明显高于准静态;在准静态加载条件下,PMMA试样呈现明显的延性破坏特征,在动态加载条件下则表现为脆性破坏.最后,对PMMA材料的ZWT粘弹性本构模型参数进行了拟合,拟合结果与实验结果吻合较好,表明该本构模型能够较好地描述较宽应变率范围内PMMA材料的应力应变关系.     

不同加载状态下TA2钛合金绝热剪切破坏响应特性

一般认为绝热剪切现象在宏观上表现为材料动态本构失稳,即热软化大于应变硬化.本文采用帽型受迫剪切试样研究TA2钛合金的动态力学特性和本构失稳过程.首先对剪切区加载应力状态进行理论和数值分析,通过合理设计帽型试样,剪切区变形可近似按剪切状态处理;结合二维数字图像相关法(two-dimensional digital image correlation,DIC-2D)直接测试试样剪切区应变演化,给出帽型受迫剪切实验的等效应力-应变响应曲线.进一步,利用Hopkinson压杆对TA2钛合金开展动态压缩及帽型剪切对比试验研究,比较压缩、剪切试验得到的等效应力-应变曲线,采用"冻结"试样方法分析试样中绝热剪切局域化演化过程,探讨不同加载状态下TA2钛合金的绝热剪切破坏现象及其动态力学响应特性.实验结果表明,在塑性变形初始阶段,动态压缩及剪切加载下的等效应力-应变曲线符合较好,但随塑性损伤发展及绝热剪切带形成,两者出现分离,表明损伤及绝热剪切演化过程与应力状态相关.剪切试样实验得到的本构"软化"特性能够反映绝热剪切带起始、破坏演化过程的力学响应特性,而在动态压缩实验中,即使试样中已出现双锥形的绝热剪切带及局部裂纹分布,其表观等效应力-应变曲线并不出现软化特征,动态压缩实验无法得到关于绝热剪切起始、发展以及破坏的本构软化响应特性.