Generalized elastic damage theory and its application to composite plate

Based on continuum mechanics, a generalized damage theory for elastic material which can be used for anisotropic composite is presented in this paper. This theory for anisotropic elastic material has been proposed here from the stress-strain relation of the actual damaged material. Introducing a fourth order damage operator that may be formed by a symmetrical second order damage factor tensor, the constitutive equation of the damaged material has been set up. The expressions of components of both the stress tensor and the strain tensor of the damaged material and their first order invariants have been also derived. The application of this theory to the 2-dimensional composite laminate, including the technique estimating the components of the damage factor tensor and the damage variable tensor and also the practical measure technique of the damage in the whole process, have been explained in detail. Finally, the changes of the anisotropic elastic properties and the actual stress state of damaged material have been discussed and some interesting results have been obtained in this paper.     

The effect of crack face contact on the anisotropic effective moduli of microcrack damaged media

The generalized self-consistent method (GSCM) in conjunction with a computational finite element method is used to calculate the anisotropic effective moduli of a medium containing damage consisting of microcracks with an arbitrary degree of alignment. Since cracks respond differently under different external loads, the moduli of the medium subjected to tension, compression and an initially stress-free state are evaluated and shown to be significantly different, which will further affect the wave speed inside the damaged media. There are four independent material moduli for a 2-D plane stress orthotropic medium in tension or compression, and seven independent material moduli for a 2-D plane stress orthotropic cracked medium, which is initially stress free. When friction exists, it further changes the effective moduli. Numerical methods are used to take into account crack face contact and friction. The wave slowness profiles for microcrack damaged media are plotted using the predicted effective material moduli.     

An anisotropic model of the Mullins effect

The Mullins effect in rubber-like materials is inherently anisotropic. However, most constitutive models developed in the past are isotropic. These models cannot describe the anisotropic stress-softening effect, often called the Mullins effect. In this paper a phenomenological three-dimensional anisotropic model for the Mullins effect in incompressible rubber-like materials is developed. The terms, damage function and damage point, are introduced to facilitate the analysis of anisotropic stress softening in rubber-like materials. A material parametric energy function which depends on the right stretch tensor and written explicitly in terms of principal stretches and directions is postulated. The material parameters in the energy function are symmetric second-order damage and shear-history tensors. A class of energy functions and a specific form for the constitutive equation are proposed which appear to simplify both the analysis of the three-dimensional model and the calculation of material constants from experimental data. The behaviour of tensional and compressive ground-state Young’s moduli in uniaxial deformations is discussed. To further justify our model we show that the proposed model produces a transversely anisotropic non-virgin material in a stress-free state after a simple tension deformation. The proposed anisotropic theory is applied to several types of homogenous deformations and the theoretical results obtained are consistent with expected behaviour and compare well with several experimental data.     

复合固体推进剂黏弹性微裂纹损伤本构模型

为建立复合固体推进剂的损伤本构模型,在介观尺度上视其为微裂纹损伤,选取微裂纹密度为损伤内变量。在Abdel-Tawab本构方程的基础上,基于微裂纹均匀化理论,推导了损伤映射张量的一般形式。该张量通常具有非完全对称性,其物理意义是将真实应力空间中各向异性材料的多轴加载映射为等效应力空间中各向同性材料的更为复杂的多轴加载。其次,基于黏弹性动态裂纹扩展模型和裂纹扩展阻力曲线的概念,建立了损伤内变量的演化方程。该演化方程仅含4个物理意义明确的细观参数,并且参数的取值规律与宏观应力曲线的变化规律相一致。数值结果表明,建立的模型能够有效反映材料损伤的应变率、温度依赖性及各向异性特征,并且具有一定的蠕变损伤预测能力。     

A comparative study of continuum damage models for crack propagation under gross yielding

In this paper, the problem of fracture initiation in an aluminum alloy thin plate containing a central crack is examined by employing several phenomenological continuum damage mechanics models. These models differ mainly in the selection of the kind of tensorial property the damage variable assumes, the nature of the equivalence postulate between damaged and pseudo undamaged material states, and the way damage evolution laws are formulated. Two formulations of damage effect tensor based on the engineering notation and the normative notation of stress and strain, respectively, are compared. In addition, the hypothesis of strain equivalence is compared to that of stress working equivalence. The error in the assumption of isotropic damage development in the crack tip process zone is also checked against that of anisotropic damage. In the numerical algorithm, both updated Lagrangian formulation and small displacement formulation of material non-linearity only are adopted and compared. The influence of non-proportionality in stress histories present in the crack tip region is accounted for by introducing a dynamic coordinate system of principal damage such that the principal direction of damage rotates in accordance with that of the loading. The calculated fracture initiation loads are finally compared with those determined experimentally.     

Euclid空间中损伤变量表示的Riemann空间中的连续性方程

损伤作为一种缺陷,宏观上通常是在Euclid空间中通过虚拟构形的方式,以连续分布的损伤变量加以描述。但如果描述更复杂的缺陷,处理变形非协调性问题,仅停留在Euclid空间中是不够的。同时,针对工程中不同材料以及同种材料的不同损伤机制,目前尚未建立一个统一的损伤模型。根据Euclid空间中的四阶损伤变量张量,定义了处于自然状态中的损伤变形体在Riemann空间中的三阶拟塑性张量、四阶异物张量,并用其描述损伤缺陷。并给出Riemann空间中异物张量所满足的连续性方程。从而建立了损伤缺陷与Riemann空间的对应关系,以Riemann空间中Bianchi恒等式刻划损伤变形体的非协调性。使得可以在Riemann这样一个弯曲空间中讨论损伤所引起的材料力学性能的劣化。最后给出一个各向异性损伤的算例。     

CONDITIONS FOR ADAPTATION OF DAMAGED ELASTIC–PLASTIC BODIES TO CYCLIC LOADING

Conditions for adaptation of isotropically damaged elastic–plastic bodies with isotropic strain hardening are investigated in the framework of the energy-based coupled elastic–plastic damage model by Ju. The yield function is assumed to be a homogeneous function of the first order in the stress tensor components. Due to this assumption, the notion of effective yield stress can be introduced. The loading program is supposed to be prescribed. Features of the stress path at the post-adaptation stage are considered, which lead to new necessary shakedown conditions expressed by a set of inequalities, and, in turn, to a problem of mathematical programming whose solution yields lower estimates for the damage and strain-hardening parameters. In the event, if the calculated value of the damage parameter is greater than its critical value, an adaptation to a given loading program is impossible. This condition is also necessary for adaptation in the case if only bounds for applied loads are prescribed. A correction of the constitutive material model is proposed which possibly could be good for ductile damage. The derived shakedown condition is not only necessary, but also sufficient for the plastic adaptation. The developed method is expounded in an example.     

修正偶应力理论层合薄板稳定性模型及尺度效应

该文将各向同性修正偶应力理论推广到各向异性,提出各向异性的细观尺度的复合材料层合板的本构方程,基于虚功原理建立了各向异性修正偶应力理论并用于建立复合材料层合薄板偶应力理论稳定性模型。该理论的偶应力部分的转角不是独立变量(称为C1理论),对于各单层引入纤维和基体材料的不同的两个材料细观参数,建立了适用于层合板/夹层板的偶应力理论模型。该理论的应变不对称,但是,用于各向同性材料与修正偶应力理论等价。为了便于工程应用,忽略基体材料的细观长度参数,建立了各单层只含一个材料细观参数的偶应力层合薄板理论稳定性模型。算例表明建立的偶应力层合板模型能用于分析层合板稳定性的尺度效应。     

基于二阶组构张量的各向异性屈服准则

将四个屈服准则:Tresca准则、Mises准则、Mohr-Coulomb准则以及Drucker-Prager准则归类为剪切屈服准则。Tresca准则和Mohr-Coulomb准则是关于最不利截面的剪切屈服准则,而Mises准则和Drucker-Prager准则是关于各方向截面的剪应力和正应力的某种综合度量的八面体剪应力和八面体正应力的剪切屈服准则。从方向函数(ODF)的概念入手,将各方向截面的剪应力和正应力综合度量直接取为所有方向截面上的剪应力和正应力的平均。对各向同性材料,提出了平均剪切屈服度准则:当平均剪应力和平均正应力的组合达到某一极限值时,材料开始屈服。研究表明,平均剪切屈服准则与Drucker-Prager准则具有相同的形式,当不考虑平均正应力对屈服的影响时,它与Mises准则具有相同的形式。针对由各向异性损伤导致的材料各向异性强度问题,定义截面上的有效正应力和有效剪应力则分别为截面上的法向力和切向力与有效承载面积之比,基于截面上的有效应力提出了各向异性材料的平均剪切屈服准则。各向异性损伤引起的截面上有效应力放大系数为方向函数,可以采用二阶组构张量来近似表示,在任意坐标系中,各向异性屈服准则为应力分量的二次齐次式,导出了其中的系数与二阶组构张量之间的显式关系式。在二阶组构张量的主轴坐标系内,各向异性屈服准则与殷有泉的拓展Hill准则形式完全相同,当不考虑正应力对屈服的影响时,它与Hill准则具有相同的形式。     

一种新的疲劳损伤演化模型

考虑到金属构件的疲劳损伤主要属于机械损伤,则造成损伤的外因主要是应力幅,而导致损伤的内因是材料本身的性质。基于损伤力学基础理论推导的一般损伤演化方程,并结合断裂力学中经典的帕里斯公式,以有效当量应力幅和材料损伤的特性参数为控制变量得到了一种新的疲劳损伤演化模型,并以12Cr1MoV钢为例进行了实验分析。结果表明:新模型形式简单、参数少,且比Lemaitre模型与实验结果符合得更好。     

Anisotropic Gurson‐Tvergaard‐Needleman plasticity and damage model for finite element analysis of elastic‐plastic problems

Implementation and analysis of the anisotropic version of the Gurson‐Tvergaard‐Needleman (GTN) isotropic damage criterion are performed on the basis of Hill's quadratic anisotropic yield theory with the definition of an effective anisotropic coefficient to represent the elastic‐plastic behavior of ductile metals. This study aims to analyze the extension of the GTN model suitable for anisotropic porous metals and to investigate the GTN model extension. An anisotropic damage model is implemented using the user material subroutine in ABAQUS/standard finite element code. The implementation is verified and applied to simulate a uniaxial tensile test on a commercially produced aluminum sheet material for three‐dimensional and plane stress test cases. Spherical and ellipsoidal micro voids are considered in the matrix material, and their effects on the uniaxial stress‐strain response of the material are analyzed. Hill's quadratic anisotropic yield theory predicts substantially large damage evolution and a low stress‐strain curve compared with those predicted by the isotropic model. An approximate model for anisotropic materials is proposed to avoid increased damage evolution. In this approximate model, Hill's anisotropic constants are replaced with an effective anisotropy coefficient. All model‐generated stress‐strain predictions are compared with the experimental stress‐strain curve of AA6016‐T4 alloy.     

A finite element approach to anisotropic damage of ductile materials in large deformations Part I: an anisotropic ductile elastic-plastic-damage model

During past decades, many material models using the continuum damage mechanics (CDM) approach have been proposed successfully in the small deformation regime to describe inelastic behaviors and fracturing phenomena of a material. For ductile materials, large deformation takes place at the level of damage appearance. Damage is anisotropic in nature. In this paper, the ductile damage at finite deformations is modeled as an anisotropic tensor quantity. Then, a fourth-order symmetric stress correction tensor is proposed for computationally efficient and easy implementation in the finite element formulations. Consequently, an explicit form of the fourth-order constitutive equations of the proposed elastic-plastic-damage model is derived. Both isotropic and kinematic hardening effects are included in the formulation. The new constitutive model can predict not only the elastic-plastic behaviors, but also the sequential variations of ductile materials. An evaluation of the constitutive and damage evolution equations is presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

Numerical implementation of anisotropic damage mechanics

This paper describes implementation of anisotropic damage mechanics in the material point method. The approach was based on previously proposed, fourth‐rank anisotropic damage tenors. For implementation, it was convenient to recast the stress update using a new damage strain partitioning tensor. This new tensor simplifies numerical implementation (a detailed algorithm is provided) and clarifies the connection between cracking strain and an implied physical crack with crack opening displacements. By using 2 softening laws and 3 damage parameters corresponding to 1 normal and 2 shear cracking strains, damage evolution can be directly connected to mixed tensile and shear fracture mechanics. Several examples illustrate interesting properties of robust anisotropic damage mechanics such as modeling of necking, multiple cracking in coatings, and compression failure. Direct comparisons between explicit crack modeling and damage mechanics in the same material point method code show that damage mechanics can quantitatively reproduce many features of explicit crack modeling. A caveat is that strengths and energies assigned to damage mechanics materials must be changed from measured material properties to apparent properties before damage mechanics can agree with fracture mechanics.     

A localizing gradient damage enhancement with micromorphic stress-based anisotropic nonlocal interactions

This article presents a localizing gradient damage model with evolving micromorphic stress-based anisotropic nonlocal interactions. The objective is to model mesh independent fracture behavior of quasi-brittle materials, and to avoid the issues associated with the existing gradient-enhanced damage models. In the proposed model, an evolving anisotropic nonlocal interaction domain governs the spatial diffusive behavior, which helps to maintain a localized damage bandwidth during the final stages of loading. The anisotropy in nonlocal interactions is captured through an anisotropic gradient tensor, which defines the orientation of the diffusive interaction domain based on the principal stresses at a given material point. In this article, a smooth micromorphic stress tensor is utilized for the determination of principal stress states, to enforce a properly oriented interaction across the bandwidth of the damage process zone throughout the loading process. The proposed approach also enables the usage of low order finite elements without any oscillatory micromorphic or nonlocal equivalent strain response in the later stages of deformation. The accuracy and performance of the proposed model are demonstrated numerically in plane strain/stress for mode-I, mode-II, and mixed-mode loading conditions.     

Anisotropic dynamic damage and fragmentation of rock materials under explosive loading

This paper describes the development of a constitutive model for predicting dynamic anisotropic damage and fragmentation of rock materials under blast loading. In order to take account of the anisotropy of damage, a second rank symmetric damage tensor is introduced in the present model. Based on the mechanics of microcrack nucleation, growth and coalescence, the evolution of damage is formulated. The model provides a quantitative method to estimate the fragment distribution and fragment size generated by crack coalescence in the dynamic fragmentation process. It takes account of the experimental facts that a brittle rock material does not fail if the applied stress is lower than its static strength and certain time duration is needed for fracture to take place when it is subjected to a stress higher than its static strength. Numerical results are compared with those from independent field tests.     

The Effect of Initial Defects on Overall Mechanical Properties of Concrete Material

Considering the fact that the initial defects, like the imperfect interfacial transition zones (ITZ) and the micro voids in mortar matrix, weaken the mechanical properties of concrete, this study develops corresponding constitutive models for ITZ and matrix, and simulates the concrete failure with finite element methods. Specifically, an elastic-damage traction-separation model for ITZ is constructed, and an anisotropic plastic-damage model distinguishing the strength-difference under tension and compression for mortar matrix is proposed as well. In this anisotropic plastic-damage model, the weakening effect of micro voids is reflected by introducing initial isotropic damage, the distinct characteristic of tension and compression which described by decomposing damage tensor into tensile and compressive components, and the plastic yield surface which established on the effective stress space. Furthermore, by tracking the damage evolution of concrete specimens suffering uniaxial tension and compression, the effects of imperfect status of ITZ and volume fraction of initial voids on the concrete mechanical properties are investigated.     

An anisotropic damage model with dilatation for concrete

An anisotropic damage model for concrete is developed within the general framework of the internal variable theory of thermodynamics. The rate of change of the compliance tensor is described in terms of kinetic relations involving a damage parameter whose increment is governed by the consistency equation associated with a pressure-dependent damage surface in stress space. The use of the compliance tensor implies that damage is reflected through a fourth-order tensor. Dilatation is obtained as a consequence of damage, and permanent deformation due to damage is addressed via a simple evolution equation. The theory is capable of accommodating the anisotropy induced by microcracking and is very suitable for computer implementation.     

On constitutive equations for thermoelastic dielectric continuum in terms of invariants

Electro-thermomechanical behavior of a thermoelastic dielectric body subject to external loading has been investigated theoretically in the present analysis. The theory is formulated in the context of continuum electrodynamics. The reaction of the body subject to external loads is expressed in symmetric stress, electrical polarization and heat flux. The solid medium is assumed to be linear, dielectric, isotropic, incompressible and dependent on temperature gradient. It has been observed that, as a result of thermodynamic constraints, the stress potential function is dependent on the deformation tensor, the electric field vector and the temperature, while the heat flux vector function is dependent on the deformation tensor, the electric field vector, the temperature and temperature gradient. To determine arguments of the stress potential and the heat flux vector functionals, findings of the theory of invariants have been used as a method because of that isotropy constraint is imposed on the material. As a result, constitutive equations of symmetric stress, polarization field and heat flux vector have been obtained in both material and spatial coordinates and asymmetric stress has been found using the expressions of symmetric stress and polarization field.     

On fibre bundle approach to a damage analysis

In this paper the geometric structure of a damaged state in the medium is investigated based on the fibre bundle technique. The damage (a failure) in some region of the medium under loads is regarded as breakdown of the holonomicity in this region. A measure of reduction of load carrying area elements caused by the development of microcracks or microvoids in the damaged medium is identified with the damage tensor. Multiplicative representations of the damage tensor formulated both on the tangent bundle with an affine connection and a fibre bundle with a non-linear connection over the damaged medium are obtained. In the former case the generalization of the damage tensor to the tangent bundle rests on the technique of prolonging the deformation vector to include not only the independent variables and dependent variables appearing in the damage tensor, but also the derivatives of the dependent variables. The so-called lifting technique is used to express the damage tensor by initial, additional and direct (deformation-induced) damages or by direct and transferred ones. In the latter one the higher-order contact geometry approach identified with the Finslerian geometry is employed to analyze the influences of elastic, inelastic, direct, transferred and initial damages on the total damage tensor. Time-dependent relationships of the damage tensor are then presented.