Relationship between strain localization and catastrophic rupture

In order to explore a prior warning to catastrophic rupture of heterogeneous media, like rocks, the present study investigates the relationship between surface strain localization and catastrophic rupture. Instrumented observations on the evolution of surface strain field and the catastrophic rupture of a rock under uniaxial compression were carried out. It is found that the evolution of surface strain field displays two phases: at the early stage, the strain field keeps nearly uniform with weak fluctuations increasing slowly; but at the stage prior to catastrophic rupture, a certain accelerating localization develops and a localized zone emerges. Based on the measurements, an analysis was performed with local mean-field approximation. More importantly, it is found that the scale of localized zone is closely related to the catastrophic rupture strain and the rupture strain can be calculated in accord with the local-mean-field model satisfactorily. This provides a possible clue to the forecast of catastrophic rupture.     

非均质脆性材料灾变性破坏的一种敏感前兆

非均质脆性材料的失效可能有两种模式:渐进式破坏与灾变性破坏. 灾变性破坏具有更大的危险性,灾变预测是一个有重要意义的难题. 基于统计细观损伤力学分析了两种破坏模式的差异,发现在每一加载步中系统卸载所释放的弹性能与损伤所耗散的能量之比值可为灾变破坏提供一种敏感的前兆. 研究发现在灾变性破坏和非灾变性的渐进式破坏两种失效模式中,虽然系统中卸载所贡献的弹性能增量与损伤消耗和释放的能量增量的比值都有趋向于1的趋势,但是两者的趋近方式不同:灾变性破坏时表征能量比随控制变量变化特征的响应函数R有一个上升段,在灾变破坏点$R$趋向于无穷大,而非灾变性破坏中$R$始终是减小的.      

静应力作用下的层合板接头损伤失效分析

基于单层板理论,结合有限元分析技术即应力分析、Hashin三维失效判定准则、包含4种基本损伤类型相互关联作用的材料性能退化方法及结构最终破坏判据等,建立了含辅助铺层层合板接头静载损伤失效分析方法. 同时,对层合板接头损伤扩展进行了模拟分析,损伤计算结果与试验分级加载试样X光进行了对比. 通过多种类型层合板接头静强度预测结果与试验结果对比及静载累积损伤规律分析表明,建立的静载三维累积损伤分析的强度预测方法可方便模拟不同结构尺寸层合板接头内部各铺层损伤起始、发展及结构最终破坏整个累积过程,同时获得其最终破坏强度及破坏模式. 该方法的预测结果与试验结果吻合较好.      

岩石变形演化诱致灾变破坏过程的同步实验观测

将试样变形场演化特征与试样的宏观载荷位移曲线的演化特征结合起来研究是揭示非均匀脆性介质变形演化诱致灾变破坏的一个重要途径.本文发展了一套实验系统,通过对试样表面变形场的演化、宏观载荷和位移信号的同步观测,对单轴加载下岩石试样变形演化和灾变破坏的过程进行了实验研究.揭示了试样变形场由加载初期的随机涨落到灾变破坏前出现明显的变形局部化的演化特征现象,试样最终在变形局部化区内形成宏观破裂面.     

Particle fracture in high-volume-fraction ceramic-reinforced metals: Governing parameters and implications for composite failure

Weibull parameters of angular alumina particles are determined from experimental tensile test data on high-ceramic-content metal matrix composites using a micromechanical model that accounts for internal damage in the form of particle cracking, the dominant damage mode in these composites. The fraction of broken particles is assessed from the drop of Young's modulus and particle fracture is assumed to be stress controlled. Two extreme load-sharing modes, namely a purely local and a global load-sharing mode, are considered to account for the load redistribution due to particle fracture. Consistent powder strength parameters can be thus “back-calculated” for particles that are embedded in different Al-Cu matrices. On the other hand, this calculation fails for pure Al matrix composites, which exhibit a much larger strain to failure than Al-Cu matrix composites. It is shown that for Al matrix composites, the role of plastic (composite) strain on particle fracture constitutes a second parameter governing particle damage. This finding is rationalized by particle-particle interactions in these tightly packed ceramic particle-reinforced composites, and by the increase of matrix stress heterogeneity that is brought with increasing plastic strain. Failure of the alloyed matrix composites is well described by the (lower bound) local load-sharing micromechanical model, which predicts a catastrophic failure due to an avalanche of damage. The same model predicts failure of pure aluminium matrix composites to occur at the onset of tensile instability, also in agreement with experimental results once the role of plastic strain on damage accumulation is accounted for.     

The growth of unstable thermoplastic shear with application to steady-wave shock compression in solids*

The catastrophic growth of unstable thermoplastic shear following the transition from homogeneous deformation to heterogeneous localized deformation through distributed shear banding is studied through approximate analytic and computational methods. The calculations provide expressions for shear band widths, spacing, catastrophic growth times and the rate of stress communication between shear bands. The optimum shear band width and spacing are found to be consistent with a minimum work principle. The model predicts that the product of the energy dissipated and the localization time in the shear localization process is invariant with respect to changes in the driving strain rate. Such behavior has been noted in the steady-wave shock compression of a number of solids. The calculations are applied to heterogeneous shear localization observed in the shock compression of aluminum.     

Two scale damage model and related numerical issues for thermo-mechanical High Cycle Fatigue

On the idea that fatigue damage is localized at the microscopic scale, a scale smaller than the mesoscopic one of the Representative Volume Element (RVE), a three-dimensional two scale damage model has been proposed for High Cycle Fatigue applications. It is extended here to anisothermal cases and then to thermo-mechanical fatigue. The modeling consists in the micromechanics analysis of a weak micro-inclusion subjected to plasticity and damage embedded in an elastic meso-element (the RVE of continuum mechanics). The consideration of plasticity coupled with damage equations at microscale, altogether with Eshelby–Kröner localization law, allows to compute the value of microscopic damage up to failure for any kind of loading, 1D or 3D, cyclic or random, isothermal or anisothermal, mechanical, thermal or thermo-mechanical. A robust numerical scheme is proposed in order to make the computations fast. A post-processor for damage and fatigue (DAMAGE_2005) has been developed. It applies to complex thermo-mechanical loadings. Examples of the representation by the two scale damage model of physical phenomena related to High Cycle Fatigue are given such as the mean stress effect, the non-linear accumulation of damage. Examples of thermal and thermo-mechanical fatigue as well as complex applications on real size testing structure subjected to thermo-mechanical fatigue are detailed.     

连续损伤力学基临界奇异指数与破坏时间预测

响应量在临近破坏时呈现出临界幂律奇异性加速特征,是一种被广泛证实的灾变破坏前兆,并被火山、滑坡和岩石破坏实验等后验预测结果证实为一种对破坏时间进行短临期预测的可行方法.但是,奇异性指数测量值的较大分散性导致了对其具体取值的争议和预测效果的不确定性.因此,理解奇异性指数取值特征及其内在物理控制因素,成为了一个核心问题.本文基于连续介质损伤力学和材料时间相关失效特征,构建了刻画损伤加速发展通向破坏过程的力学模型.导出了恒名义应力蠕变加载和控制名义应力随时间线性增大两种典型加载方式下,损伤和应变率加速发展通向破坏的临界幂律奇异性前兆特征.阐明了临界幂律奇异性指数取值依赖于材料损伤与承受真应力之间的非线性关系这一内在物理根源,表明了实际测量中奇异性指数的分散性不完全归结于测量数据误差,而是有着内在物理控制因素.针对破坏前奇异性指数的不确定性,建议了在未知奇异性指数条件下预测破坏时间的方法,并基于花岗岩脆性蠕变破坏实验进行了验证和说明.      

Length scale effects on the shear localization process in metallic glasses: A theoretical and computational study

Some recent experiments on sub-micron and nano-sized metallic glass (amorphous alloy) specimens have shown that the shear localization process becomes more stable and less catastrophic when compared to the response exhibited by large sample sizes. This leads to the discovery that the shear localization process and fracture can be delayed by decreasing sample volume. In this work we develop a non-local and finite-deformation-based constitutive model using thermodynamic principles and the theory of micro-force balance to study the causes for the aforementioned observations. The constitutive model has also been implemented into a commercially available finite-element program by writing a user-material subroutine. With the aid of finite-element simulations, our constitutive model predicts that metallic glass samples have the intrinsic ability to exhibit: (a) the delaying of (catastrophic) shear localization with decreasing sample size, and (b) homogeneous deformation behavior for sample volumes smaller than the shear band nucleus.The cause for the observations listed above is the increasing influence of a non-local interaction stress with decreasing sample volume. This interaction stress has energetic origins and it affects plastic deformation due to the strong coupling between plastic shearing and free-volume generation. Akin to strain-gradient plasticity theory, the role of the interaction stress is to strengthen the material at locations where the defect density/free volume is higher compared to the rest of metallic glass sample.     

On the equivalence between traction- and stress-based approaches for the modeling of localized failure in solids

This work investigates systematically traction- and stress-based approaches for the modeling of strong and regularized discontinuities induced by localized failure in solids. Two complementary methodologies, i.e., discontinuities localized in an elastic solid and strain localization of an inelastic softening solid, are addressed. In the former it is assumed a priori that the discontinuity forms with a continuous stress field and along the known orientation. A traction-based failure criterion is introduced to characterize the discontinuity and the orientation is determined from Mohr's maximization postulate. If the displacement jumps are retained as independent variables, the strong/regularized discontinuity approaches follow, requiring constitutive models for both the bulk and discontinuity. Elimination of the displacement jumps at the material point level results in the embedded/smeared discontinuity approaches in which an overall inelastic constitutive model fulfilling the static constraint suffices. The second methodology is then adopted to check whether the assumed strain localization can occur and identify its consequences on the resulting approaches. The kinematic constraint guaranteeing stress boundedness and continuity upon strain localization is established for general inelastic softening solids. Application to a unified stress-based elastoplastic damage model naturally yields all the ingredients of a localized model for the discontinuity (band), justifying the first methodology. Two dual but not necessarily equivalent approaches, i.e., the traction-based elastoplastic damage model and the stress-based projected discontinuity model, are identified. The former is equivalent to the embedded and smeared discontinuity approaches, whereas in the later the discontinuity orientation and associated failure criterion are determined consistently from the kinematic constraint rather than given a priori. The bi-directional connections and equivalence conditions between the traction- and stress-based approaches are classified. Closed-form results under plane stress condition are also given. A generic failure criterion of either elliptic, parabolic or hyperbolic type is analyzed in a unified manner, with the classical von Mises (J₂), Drucker–Prager, Mohr–Coulomb and many other frequently employed criteria recovered as its particular cases.     

三维编织复合材料渐进损伤的非线性数值分析

基于考虑纤维束相互挤压的八边形纤维束截面单胞模型,引入周期性位移边界条件,采用 细观非线性有限元方法,建立了三维四向编织复合材料的渐进损伤拉伸强度模型. 该模型考 虑了增强体纤维束纵向非线性剪切应力-应变关系,采用Hashin型损伤失效准则定义了纤维 束的典型损伤类型,并根据纤维束和纯基体相应损伤类型所造成的材料性能退化,模拟了不 同编织角试件各类损伤产生、扩展及材料最终破坏的整个过程. 模型数值结果与实验数据吻 合较好,证明了该模型的合理有效性. 探讨了组分材料剪切非线性、损伤对材料宏观非线性 本构行为的影响,结果表明：随着编织角增大,纤维束剪切非线性效应和累积损伤对材料非 线性力学行为的影响明显增强.     

Multiresolution modeling of ductile reinforced brittle composites

Tungsten carbide-cobalt (WC-Co) is an important ductile reinforced brittle composite used in a range of important applications. The relationship between microstructure and mechanical properties of WC-Co is truly multiscale; micromechanical processes interact at different scales, resulting in permanent plastic deformation, damage accumulation and final failure of the composite. The goal of the current paper is to develop a continuum-based model, which captures the progressively finer scales of strain localization observed in WC-Co composites during plastic deformation and failure. This is achieved via a set of multiresolution governing equations; a microstress is introduced at each scale of strain localization, which represents the resistance to inhomogeneous strain localization at that scale. The extra constitutive models associated with these microstresses can be elucidated from the average response of separate computational cell models of a representative microstructure. The final multiresolution continuum model is capable of capturing the important length scales of deformation during the plastic stage of deformation without resorting to modeling microstructural scale features directly. The result is a more realistic continuum model; in particular the fracture toughness prediction is more physical when these length scales are incorporated compared to a conventional continuum approach.     

固体结构损伤破坏统一相场理论、算法和应用

固体开裂引起的损伤和断裂是工程材料和结构最为普遍的破坏形式.为了防止这种破坏,结构设计首先必须了解裂缝在固体内如何萌生、扩展、分叉、汇聚甚至破碎;更重要的是,还需要准确量化这些裂缝演化过程对于结构完整性和安全性降低的不利影响.针对上述固体结构损伤破坏问题,本工作系统地介绍了笔者提出的统一相场理论、算法及其应用.作为一种...     

Influence of damage on properties of strain localization in geomaterials at plane stress and plane strain

Summary By regarding geomaterials under loading as a mixture of intact and damaged parts, we investigate the influence of damage on the properties of strain localization in elastoplastic geomaterials at plane stress and plane strain. Conditions for the onset of strain localization including the effects of damage are derived for the cases of plane strain and plane stress. Discussed are the inclination of the localized band and the hardening modulus corresponding to the onset of strain localization. It is shown that the properties of the strain localization are dependent on the damage and the capacity of bearing hydrostatic pressure by the damaged part, and that damage may induce an earlier onset of strain localization and lead to instability of a geomaterial.accepted for publication 11 March 2004     

岩石单轴压缩作用下变形局部化的梯度塑性解

采用梯度塑性理论研究单轴压缩作用下岩石变形局部化,得到了单轴压缩作用下岩石变形局部化带宽度的一维、二维解析解,为实验测定内部材料长度参数提供了理论依据.     

Modeling adiabatic shear failure from energy considerations

The numerical simulation of dynamic structural failure by localized shear is quite complex in terms of constitutive models and choice of adequate failure criteria, along with a pronounced mesh-sensitivity. As a result, the existing numerical procedures are usually quite sophisticated, so that their application for design purposes is still limited. This study is based on the implementation of a simple energy-based criterion, which was developed on experimental considerations (Rittel et al., 2006), and uses a minimal number of adjustable parameters. According to this criterion, a material point starts to fail when the total strain energy density reaches a critical value. Thereafter, the strength of the element decreases gradually to zero to mimic the actual structural behavior. The criterion was embedded into commercial finite element software and tested by simulating numerically four typical high-rate experiments. The first is the dynamic torsion test of a tubular specimen. The second concerns the failure mode transition in mode II fracture of an edge crack in plain strain. The last two involve dynamic shear localization under high rate compression of a cylindrical and a shear compression specimen. A very good adequation was found both qualitatively and quantitatively. Qualitatively, in terms of failure path selection, and quantitatively, in terms of local strains, temperatures and critical impact velocity. The proposed approach is enticing from an engineering perspective aimed at predicting the onset and propagation of dynamic shear localization in actual structures.     

Adaptive coupling between damage mechanics and peridynamics: A route for objective simulation of material degradation up to complete failure

The objective (mesh-independent) simulation of evolving discontinuities, such as cracks, remains a challenge. Current techniques are highly complex or involve intractable computational costs, making simulations up to complete failure difficult. We propose a framework as a new route toward solving this problem that adaptively couples local-continuum damage mechanics with peridynamics to objectively simulate all the steps that lead to material failure: damage nucleation, crack formation and propagation. Local-continuum damage mechanics successfully describes the degradation related to dispersed microdefects before the formation of a macrocrack. However, when damage localizes, it suffers spurious mesh dependency, making the simulation of macrocracks challenging. On the other hand, the peridynamic theory is promising for the simulation of fractures, as it naturally allows discontinuities in the displacement field. Here, we present a hybrid local-continuum damage/peridynamic model. Local-continuum damage mechanics is used to describe “volume” damage before localization. Once localization is detected at a point, the remaining part of the energy is dissipated through an adaptive peridynamic model capable of the transition to a “surface” degradation, typically a crack. We believe that this framework, which actually mimics the real physical process of crack formation, is the first bridge between continuum damage theories and peridynamics. Two-dimensional numerical examples are used to illustrate that an objective simulation of material failure can be achieved by this method.     

水岩耦合变形破坏过程及机理研究进展

坝址水力劈裂、煤矿和隧道突水以及水力压裂工程都是岩体损伤演化积累诱发破坏 (灾变) 的过程,然而能否正确认识和评价这些过程取决于人们对水岩耦合变形破坏过程和机理的研究程度. 阐述了岩石渗流 --变形 --破坏特性及机理研究的进程, 侧重于物理实验成果. 综述了水岩耦合变形破坏过程分析模型和数值计算, 特别是有限单元法处理裂纹的思路, 并介绍了作者课题组近年来所取得的一些代表性成果. 给出了当前理论 (数值) 研究中所面临的几个具体问题和挑战, 并对今后研究工作进行了力所能及的展望.