冲击压缩下准脆性材料含微裂纹损伤的本构模型

基于准脆性材料中翼型拉伸裂纹的成核准则,运用细观损伤理论推导了翼型裂纹损伤对材料弹性模量的弱化作用.考虑裂纹扩展对材料动态断裂的滞后效应,建立了动态裂纹扩展准则,并给出损伤演化方程,在此基础上建立了准脆性材料单轴冲击压缩下的动态损伤本构模型.结合氧化铝陶瓷材料独特的力学响应和破坏特性,讨论了模型中微裂纹成核参数、微裂纹尺寸对动态断裂强度的影响,并用该模型计算了单轴压缩下氧化铝陶瓷的应力应变曲线,数值结果与实验结果吻合良好.     

基于应变空间的碾压混凝土各向异性损伤本构模型

根据碾压混凝土材料的力学特性和损伤拉压显著不同的特点,分别在拉应变和压应变空间建立了碾压混凝土的本构关系和损伤演化方程.在拉应变空间,碾压混凝土的变形特性表现为脆弹性,考虑弹性与损伤耦合,应用正交各向异性损伤理论描述碾压混凝土的刚度退化和应变软化;在压应变空间,考虑弹塑性与损伤耦合,应用内时理论来描述碾压混凝土的弹塑性特性,正交各向异性损伤理论来描述微裂缝扩展引起的刚度退化和应变软化,内时理论没有屈服面,使模型的参数和方程大大减少,从而简化了非线性计算过程.计算结果表明,该模型能够较好地描述碾压混凝土在单轴和多轴加载下的性质.     

2.5D C/SiC复合材料连续损伤本构模型

基于连续损伤力学建立了一种包含拉伸与剪切损伤变量的2.5D C/SiC复合材料连续损伤本构模型。分别开展了拉伸和剪切试验,获得应力-应变曲线,并通过拟合试验曲线获得各损伤变量的演化参数。采用子程序技术将本构模型嵌入商用有限元软件ANSYS,应用有限元法计算了材料的应力-应变曲线。考虑了拉剪损伤耦合效应,计算了偏轴拉伸情况下的应力-应变曲线。结果表明:沿经纱拉伸、沿纬纱拉伸以及面内剪切的应力-应变曲线与试验结果吻合,最大偏差依次为4.30%、3.09%及3.73%;偏轴拉伸计算与试验应力-应变曲线也吻合较好。      

准脆性材料中椭圆形微裂纹的生长与演化

微裂纹的生长与演化是导致准脆性材料损伤及破坏的本质,本工作对承载过程中准脆性材料内部的椭圆形微裂纹的生长与演化规律进行研究。采用复势函数求解受远场载荷作用下代表性单元中椭圆形微裂纹的变形,讨论了椭圆形微裂纹初始取向的变化对微裂纹尺寸增长的影响,并结合微裂纹扩展准则推导出损伤开始时的临界应力条件。基于翼型裂纹扩展过程的能量守恒关系,建立了损伤阶段的本构关系。     

循环荷载作用下超高性能混凝土的轴拉力学性能及本构关系模型

对具有不同拉伸应变特性(应变强化和应变软化)的超高性能混凝土(Ultra high performance concrete, UHPC)进行了单调和循环荷载作用下的直接拉伸试验。试验结果表明:应变强化UHPC基体开裂后进入多点微裂纹分布的应变强化段,达到极限抗拉强度后进入单缝开裂的应变软化段;应变软化UHPC基体开裂后直接进入单缝开裂的应变软化段;循环荷载下两种类型UHPC的轴拉应力-应变曲线包络线与单调荷载下的应力-应变曲线基本一致;基于刚度退化过程建立了两种类型UHPC的轴拉损伤演化方程,根据实测应力-应变曲线和试件的裂缝分布形态建立了两种类型UHPC的轴拉本构关系模型,与试验结果基本吻合;采用能量法研究了应变强化UHPC两阶段轴拉本构关系在数值计算时的等效方法。最后,通过无筋应变强化UHPC抗弯试验梁的数值模拟对本文建立的应变强化UHPC轴拉本构关系模型和损伤演化方程及相关假定进行了验证,结果表明本文建立的应变强化UHPC轴拉本构模型能较好地预测UHPC弯拉构件的极限承载力,轴拉损伤变量能在宏观层面上较好地反应试件的裂缝分布状态。      

基于能量的弹塑性损伤实用本构模型

建立一个实用的弹塑性损伤本构模型。在有效应力空间采用经验公式计算塑性变形,能将模型塑性部分与损伤部分解耦,降低模型的数值处理复杂性,同时大大简化模型塑性应变的计算。结合不可逆热力学理论,基于损伤能量释放率建立损伤准则,损伤能量释放率由修正后的弹性Helmholtz自由能导出,模型中将弹性Helmholtz自由能分解为应力球量部分和应力偏量部分,将其应力球量部分产生的损伤取为零,同时根据应力状态引入折减系数对其应力偏量部分进行修正,使得模型能较为准确的模拟混凝土材料在双轴加载下的本构行为。将应力张量谱分解为正、负两部分以分别定义材料受拉、受压损伤演化,并采用受拉损伤变量、受压损伤变量分别模拟混凝土材料在拉、压加载下的本构特性。引入一个加权损伤变量使得模型能较准确的反映混凝土材料的“拉-压软化效应”。最后该文给出初步试验验证,证明了该文模型的有效性。     

三维正交机织玻璃纤维/环氧树脂复合材料动态力学性能的实验和理论研究

通过实验系统研究了三维正交机织玻璃纤维/环氧树脂复合材料厚度方向和面内方向的动、 静态压缩力学性能。结果表明, 动、 静态压缩载荷作用下该材料响应表现出明显的各向异性、 非线性和应变率敏感性。针对三维正交机织玻璃纤维/环氧树脂复合材料高速变形过程中不同形式的内部缺陷和微损伤的演化, 提出了一个依赖应变、 应变率的宏观损伤量, 建立了一种含损伤的非线性黏弹性本构模型。采用数据处理方法拟合了其本构方程材料参数, 在加载过程中, 模型计算值与实验结果吻合较好。     

冲击荷载作用下混凝土动态本构模型的研究

基于混凝土冲击荷载作用下的实验研究,以修正Ottosen四参数破坏准则为屈服法则,引入损伤,构造了一个动态本构模型用于描述混凝土材料的冲击特性。宏观上,假设混凝土材料是一个均匀连续体;而从细观角度来看,混凝土材料内部存在大量随机分布的微裂纹损伤。假设微裂纹均匀分布,且符合理想微裂纹体系统条件,定义含裂纹材料中单位体积内微裂纹所占的比例来表征微裂纹损伤所引起的混凝土材料宏观力学性能的劣化,并给出了损伤的演化方程。通过模型计算模拟结果与实验结果比较发现,模拟曲线与实验曲线拟合良好,因而可以用该模型模拟混凝土材料在冲击荷载下的动态特性。     

考虑剪切非线性影响的复合材料连续损伤模型及损伤参数识别

在复合材料正交各向异性连续损伤力学的基础上,引入由应变描述的连续综合变量,损伤演化方案采用综合变量的指数函数。考虑复合材料剪切非线性,通过非线性损伤因子在全局损伤出现之前将其引入本构方程。编写了VUMAT实施损伤模型,介绍连续损伤模型中材料参数的取值方法。对无法通过试验获得的损伤参数,采用遗传算法原理,定义个体适应度为试验及模拟载荷-位移曲线误差,通过编写MATLAB程序调用ABAQUS有限元模型实施参数反演方法。选择了常用铺层结构的复合材料带孔层合板进行仿真参数识别分析,验证了基于遗传算法参数识别方法的可行性。     

二维编织C/SiC复合材料非线性损伤本构模型与应用

基于二维编织C/SiC复合材料的基本力学性能试验, 建立了该材料的宏观正交各向异性非线性损伤本构模型。模型中以可检测的应变作为变量, 采用形式简单的函数分别描述了单轴拉伸和剪切加载下的材料损伤演变下的应力-应变关系, 以及卸载状态的刚度变化规律。同时, 考虑了材料的单边效应以及拉压应力状态转换时的损伤钝化行为。将此本构模型编写成UMAT子程序并引入ABAQUS有限元软件, 可以完整描述该材料的加载非线性和卸载线性的应力-应变关系特征, 及其加卸载历史。通过对带孔板的拉伸模拟, 孔边应变分布与试验结果吻合较好, 验证了本构模型的有效性。     

Low-cycle fatigue model of damage accumulation - The strain approach

The paper presents a model of damage accumulation designed to analyse fatigue life of structural elements exploited in multiaxial, non-proportional, low-cycle loading conditions. The discussed approach consists of two calculation blocs. In the first bloc the components of stress and strain tensor are determined. This module, in which Mroz’s multisurface model was used, contains constitutive relations and the law of kinematic hardening. The second bloc contains the dependencies which determine the growth of anisotropic measure of damage accumulation (associated with the physical plane) and crack initiation criterion. The growth of the damage accumulation measure was associated with the loading damage accumulation function and the increment of non-dilatational plastic strain on the physical plane. It was assumed that crack initiation occurs when stress or a measure of damage accumulation on any physical plane reaches critical values.     

A Damage-Mode Based Three Dimensional Constitutive Model for Fibre-Reinforced Composites

This article presents a three dimensional constitutive model for anisotropic damage to describe the elastic-brittle behavior of unidirectional fibrereinforced laminated composites. The primary objective of the article focuses on the three dimensional relationship between damage of the material and the effective elastic properties for the purpose of stress analysis of composite structures, in extension to the two dimensional model in Matzenmiller, Lubliner and Taylor (1995). A homogenized continuum is adopted for the constitutive theory of anisotropic damage and elasticity. Damage initiation criteria are based on Puck failure criterion for first ply failure and progressive micro crack propagation is based on the idea of continuum damage evolution. Internal variables are introduced to describe the evolution of the damage state under loading and as a subsequence the degradation of the material stiffness. Emphasis is placed on a suitable coupling among the equations for the rates of the damage variables with respect to the different damage modes.     

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 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.     

Multi-axial thermo-mechanical analysis of power plant components from 9–12% Cr steels at high temperature

The aim of this paper is to analyze local changes of stress and strain states in a power plant component under a transient thermal environment. A robust constitutive model is developed to describe inelastic behavior of advanced 9–12% Cr heat-resistant steels at high temperature and in a multi-axial stress state. The model includes the constitutive equation for the inelastic strain rate tensor, the evolution equation for a tensor-valued state variable to reflect hardening/recovery processes and two evolution equations for two scalar-valued variables that characterize softening and damage states. The model is calibrated against experimental creep curves and verified for inelastic responses under different isothermal and non-isothermal loading paths. Steam temperature and loading profiles that correspond to an idealized start-up, holding and shut-down sequence of a power plant component are assumed. To estimate the thermal fields, transient heat transfer analysis is performed. The results are applied in the subsequent structural analysis using the developed inelastic constitutive model. The outcome is a multi-axial thermo-mechanical fatigue loop which can be used for damage assessment.     

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.     

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.     

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.     

Modeling of continuum damage for application in elastic-viscoplastic constitutive equations

A procedure is described for including isotropic and directional damage as load-history dependent softening variables in a set of elastic-viscoplastic constitutive equations. The evolution equation proposed for isotropic damage integrates to an exponential form for the case of constant stress. Directional damage is represented as a second-order symmetric tensor with a scalar effective value used in the constitutive equations. A method is proposed for treating directional damage in the case of non-proportional loading histories. Comparisons are given of uniaxial creep test results for an alloy at high temperatures with calculations based on the constitutive equations with the inclusion of isotropic damage.     

Simplified analysis of steady-state crack growth of piezoelectric ceramics based on the continuum damage mechanics

Summary A simplified analysis of crack growth of piezoelectric ceramics was performed by using the double cantilever beam (DCB) model. The analysis was intended for the crack growth in a steady state by using some assumptions to simplify the analysis. In the present paper, damage within piezoelectric ceramics was represented by a damage variable based on the continuum damage mechanics, and the effect of the damage on material properties was taken into account in a constitutive equation of the piezoelectric ceramics. The constitutive equation and an evolution equation of the damage variable were applied to the DCB model. Then, the crack growth in the DCB model was dealt with within the framework of the continuum damage mechanics. As a result of the analysis, a constant crack growth rate and distribution of stress, strain and damage variable in front of the crack tip were elucidated, and the effects of the applied electric field on the crack growth were discussed in particular.