An anisotropic damage–plasticity model for saturated quasi‐brittle materials

This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd.     

A macroscopic damage‐plastic constitutive law for modeling quasi‐brittle fracture and ductile behavior of concrete

A new phenomenological macroscopic constitutive model for the numerical simulation of quasi‐brittle fracture and ductile concrete behavior, under general triaxial stress conditions, is presented. The model is particularly addressed to simulate a wide range of confinement stress states, as also, to capture the strong influence of the mean stress value in the concrete failure mechanisms. The model is based on a two‐surface damage‐plastic formulation. The mechanical behavior in different domains of the stress space is separately described by means of a quasi‐brittle or ductile material response:

相似文献   

砂岩脆性临界破坏声发射信息应力比分析

根据岩石脆性破坏的本构方程和三维重整化理论,推导出了岩石破裂前声发射信号激增的临界点所对应的应力和岩石峰值应力比值的数学表达式。通过测试深部紫红色砂岩单轴加载变形破裂全过程中力学与声发射特征,得到了全过程力学特征曲线、声发射能量累计数等相关数,发现岩石破裂前声发射突增点所对应的应力与岩石临界点相对应。单轴刚性加载条件下岩石脆性破裂前声发射突增点所对应的应力与峰值应力比值大部分近似在74%左右,误差在±9%以内。     

Discrete element model for quasi‐brittle rupture under tensile and compressive loading

In the absence of initial cracks, the material behavior is limited by its strength, usually defined in homogeneous conditions (of stress and strain). Beyond this limit, in quasi‐brittle case, cracks may propagate and the material behavior tends to be well described by fracture mechanics. Discrete element approaches show consistent results dealing with this transition during rupture. However, the calibration of the parameters of the numerical models (i.e., stiffness, strength, and toughness) may be quite complex and sometimes only approximative. Based on a brittle rupture criterion, we analyze the biaxial response of uncracked samples. Thus, tensile and compressive strengths are analytically identified and become direct parameters of our discrete model. Furthermore, a physically reliable crack initiation (and subsequent propagation) is shown to be induced during rupture and verified by the simulation of three‐point bending and diametral compression tests. Copyright © 2016 John Wiley & Sons, Ltd.     

开挖与支护应力路径下硬岩破坏过程的真三轴与声发射试验研究

针对硬岩的脆性破坏问题,利用真三轴和声发射（AE）试验系统研究了灰岩试样在不同应力路径和应力水平下的力学特征与变形破坏机理、AE活动特征及能量释放规律。试验研究表明：（1）开挖卸荷应力路径下,试样强度受卸荷作用影响较大,且随中间主应力呈区间性变化,试样以脆性劈裂破坏为主;（2）开挖卸荷+支护的应力路径下,支护力显著改善了岩石的力学性质,破坏形式以张剪性破坏为主;（3）试样AE活动过程可以分为三个活跃期（裂隙压密阶段、卸荷阶段和失稳-破坏阶段）和两个相对平静期（弹性变形阶段和微破裂稳定发展阶段）,AE信号表征参数“absolute energy”作为能量度量指标能够较好地定量描述和研究硬岩的脆性破坏过程。试验结果和许多工程开挖和支护过程中岩体的力学行为和支护作用类似,可为工程现场开挖支护设计提供试验依据和有益指导。     

基于损伤过程的非均质岩体声发射及破坏规律数值模拟研究

声发射事件数大小及源分布特征是岩体损伤演化过程的直接表征。在前人研究的基础上,以数值模拟的方法,考虑细观非均匀特性,采用基于应变软化的弹-脆-塑性模型描述基元的损伤与破坏过程。借助FLAC3D软件中的FISH函数实现了以Weibull分布随机赋予基元材料属性及对声发射事件数的记录,从时间序列上探讨非均匀性对岩体声发射现象的影响,结果表明,均质度对于岩体的声发射规律及破坏模式影响显著。     

盐岩三轴压缩破坏与声波声发射特征试验研究

利用岩石声波-声发射一体化测试系统与应力加载装置,系统地研究了三轴压缩试验中盐岩的变形破坏特征与声波声发射活动规律。结果表明：（1）三轴加载试验中,盐岩的破坏过程分为压密与弹性压缩阶段、裂纹产生与稳定增长阶段、裂纹加速增长阶段和应变硬化阶段。阶段Ⅰ中,盐岩纵横波波速略微上升,并伴随少量的声发射事件;阶段Ⅱ中,盐岩的声发射活动开始活跃,超声波波速稳定;阶段Ⅲ和阶段Ⅳ中,岩石的声发射事件数量大幅度上升,纵波波速保持稳定,横波波速呈现明显地下降趋势。（2）盐岩的声波与声发射特征与试验围压应力密切相关。围压水平越低,受载盐岩破坏过程中的波速变化率越大;围压水平越高,盐岩超声波波速变化率越小。根据对试验数据的分析,当盐岩轴向应变达到10%,围压条件为5、10、15、20 MPa时盐岩的声发射事件数量分别为16 271、8 764、3 041、906个,围压对盐岩破坏过程的声发射活动有较大的影响,显示了围压致密效应。     

基于声发射的盐岩变形破坏过程的分形与损伤特征研究

为研究盐岩变形破坏过程中损伤变量和分形维数之间的关系,对取自某地的纯盐岩开展了单轴压缩和三轴压缩试验。在基于声发射振铃计数的损伤模型和基于声发射定位点空间演化的分形维数计算的基础上,对盐岩变化破坏过程中的损伤变量和分形维数进行了研究。研究表明：在盐岩变形破坏过程中,盐岩分形维数逐渐降低,损伤变量逐渐增加,且分形维数下降的各个阶段与损伤变量增加的各个阶段相对应;盐岩变形破坏前分形维数不再明显下降,且不同应力状态下,盐岩变形破坏前分形维数不同;随着围压增大,盐岩变形破坏前分形维数逐渐降低,其中单轴压缩和三轴压缩应力状态时,分形维数分别低于2.42和2.31、2.20时,预示着试件内部损伤开始汇集,逐渐形成宏观破裂面,导致试件变形破坏;围压对盐岩声发射活动有明显的抑制作用,随着围压的增大,盐岩变形破坏过程中声发射活动逐渐减少,且在此过程中可以发现：在应力加载初期,盐岩分形维数快速下降时所对应的应力百分比逐渐增大,在应力加载后期,当盐岩分形维数不再明显减小时对应的应力百分比逐渐增大;在盐岩变形破坏的加载初期,损伤变量较小,能量释放较少;盐岩变形破坏前损伤变量增加较快,能量快速释放。     

A rate‐dependent cohesive crack model based on anisotropic damage coupled to plasticity

In quasi‐brittle material the complex process of decohesion between particles in microcracks and localization of the displacement field into macrocracks is limited to a narrow fracture zone, and it is often modelled with cohesive crack models. Since the anisotropic nature of the decohesion process in separation and sliding is essential, it is particularly focused in this paper. Moreover, for cyclic and dynamic loading the unloading, load reversal (including crack closure) and rate dependency are essential features that are included in a new model. The modelling of degradation is based on a ‘localized’ version of anisotropic continuum damage coupled to inelasticity. The concept of strain energy equivalence between the states in the effective and nominal settings is adopted in order to define the free energy of the interface. The proposed fracture criterion is of the Mohr type, with a smooth transition of the failure and kinematics (slip and dilatation) characteristics between tension and shear. The chosen potential, of the Lemaitre‐type, for evolution of the dissipative processes is additively decomposed into plastic and damage parts, and non‐associative constitutive equations are obtained. The constitutive equations are integrated by applying the backward Euler rule and by using Newton iteration. The proposed model is assessed analytically and numerically and a typical calibration procedure for concrete is proposed. Copyright © 2006 John Wiley & Sons, Ltd.     

Homogenization‐based analysis of anisotropic damage in brittle materials with unilateral effect and interactions between microcracks

This paper is devoted to micromechanical modeling of induced anisotropic damage in brittle geomaterials. The formulation of the model is based on a proper homogenization procedure by taking into account unilateral effects and interactions between microcracks. The homogenization procedure is developed in the framework of Eshelby's inclusion solution and Ponte‐Castaneda and Willis (J. Mech. Phys. Solids 1995; 43 :1919–1951) estimate. The homogenization technique is combined with the thermodynamics framework at microscopic level for the determination of damage evolution law. A rigorous crack opening–closure transition condition is established and an energy‐release‐rate‐based damage criterion is proposed. Computational aspects on the implementation of micromechanical model are also discussed. The proposed model is evaluated by comparing numerical predictions with experimental data for various laboratory tests on concrete. Parametric studies on unilateral effects and influences of microcracks interactions are finally performed and analyzed. Copyright © 2008 John Wiley & Sons, Ltd.     

A bounding surface elasto‐viscoplastic constitutive model for non‐isothermal cyclic analysis of asphaltic materials

An elasto‐viscoplastic constitutive model for asphaltic materials is presented within the context of bounding surface plasticity theory, taking into account the effects of the stress state, void binder degree of saturation, temperature and strain rate on the material behaviour. A stress state dependent non‐linear elasticity model is introduced to represent time‐independent recoverable portion of the deformation. The consistent visco‐plasticity framework is utilised to capture the rate‐dependent, non‐recoverable strain components. The material parameters introduced in the model are identified, and their determination from conventional laboratory tests is discussed. The capability of the model to reproduce experimentally observed response of asphaltic materials is demonstrated through numerical simulations of several laboratory test data from the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Damage–viscoplastic consistency model with a parabolic cap for rocks with brittle and ductile behavior under low‐velocity impact loading

This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd.     

A coupled chemo‐thermo‐hygro‐mechanical model of concrete at high temperature and failure analysis

A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 John Wiley & Sons, Ltd.     

Coupling between progressive damage and permeability of concrete: analysis with a discrete model

A significant increase of the permeability of concrete upon micro‐cracking and a good correlation between the evolution of damage (material stiffness) and permeability are observed experimentally. The present contribution investigates this correlation theoretically, with the help of lattice analyses. Scaling analysis of lattices which contain elastic brittle bonds has shown that the material degradation should be described by the evolution of the material stiffness, or compliance, in a continuum setting (damage models). This result is reviewed and further documented in the first part of the paper. In the second part, hydro‐mechanical problems are considered with the construction of a hydraulic lattice, dual to the mechanical one. We observe that the average permeability upon micro‐cracking is the lattice scale‐independent controlling variable in the hydraulic problem. Additionally, results show that in a continuum poro‐mechanical approach, the evolution of the material permeability ought to be related to the elastic unloading stiffness, described e.g. with the help of continuum damage variables. Copyright © 2005 John Wiley & Sons, Ltd.     

Adaptive finite element analysis of mode I fracture in cement‐based materials

Microscopic studies using advanced experimental techniques have provided better insight into the fracture mechanisms in cement‐based materials. A clear understanding of fracture mechanisms is critical for the development of rigorous computational models for analysing fracture. Fracture analysis is usually carried out by finite element method. Accuracy of FE analysis depends upon the choice of mesh and for the predictions to be reliable, discretization errors are to be minimized. In cohesive crack approach, the non‐linearity is limited to the boundary conditions along the geometric discontinuity while the bulk of the material retains its elastic nature. The paper presents a mesh‐adaptive strategy based on ZZ error estimator to model discrete crack propagation in cement‐based materials. Examples of simulations have demonstrated the potential of the mesh‐adaptive technique in modelling the evolution of the localized strain profiles as well as failure of concrete test specimen. Copyright © 2001 John Wiley & Sons, Ltd.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Numerical modeling of dynamic rock fracture with a combined 3D continuum viscodamage‐embedded discontinuity model

This paper deals with numerical modeling of dynamic failure phenomena in rate‐sensitive quasi‐brittle materials, such as rocks, with initial microcrack populations. To this end, a continuum viscodamage‐embedded discontinuity model is developed and tested in full 3D setting. The model describes the pre‐peak nonlinear and rate‐sensitive hardening response of the material behavior, representing the fracture‐process zone creation, by a rate‐dependent continuum damage model. The post‐peak response, involving the macrocrack creation accompanied by exponential softening, is formulated by using an embedded displacement discontinuity model. The finite element implementation of this model relies upon the linear tetrahedral element, which seems appropriate for explicit dynamic analyses involving stress wave propagation. The problems of crack locking and spreading typical of embedded discontinuity models are addressed in this paper. A combination of two remedies, the inclusion of viscosity in the spirit of Wang's viscoplastic consistency approach and introduction of isotropic damaging into the embedded discontinuity model, is shown to be effective in the present explicit dynamics setting. The model performance is illustrated by several numerical simulations. In particular, the dynamic Brazilian disc test and the Kalthoff–Winkler experiment show that the present model provides realistic predictions with the correct failure modes and rate‐dependent tensile strengths of rock at different loading rates. The ability of initial embedded discontinuity populations to model the initial microcrack populations in rocks is also successfully tested. Copyright © 2016 John Wiley & Sons, Ltd.     

A constitutive model for the dynamic and high‐pressure behaviour of a propellant‐like material: Part I: Experimental background and general structure of the model

This paper is the first part of a work that aims at developing a mechanical model for the behaviour of propellant‐like materials under high confining pressure and strain rate. The behaviour of a typical material is investigated experimentally. Several microstructural deformation processes are identified and correlated with loading conditions. The resulting behaviour is complex, non‐linear, and characterized by multiple couplings. The general structure of a relevant model is sought using a thermodynamic framework. A viscoelastic‐viscoplastic‐compaction model structure is derived under suitable simplifying assumptions, in the framework of finite, though moderate, strains. Model development, identification and numerical applications are given in the companion paper. Copyright © 2001 John Wiley & Sons, Ltd.