基于弹塑性损伤本构模型的复合材料层合板破坏荷载预测

在连续损伤力学和塑性力学框架内,建立一个同时考虑塑性效应和损伤累积导致材料属性退化的复合材料弹塑性损伤本构模型。基于最近点投影回映算法,开发本构模型的应变驱动隐式积分算法以更新应力及与解答相关的状态变量,并推导与所开发算法相应的数值一致性切线刚度矩阵,保证有限元分析采用NewtonRaphson迭代法解答非线性问题的计算效率。采用断裂带模型对已开发的本构模型软化段进行规则化,以减轻有限元分析结果的网格相关性问题。对损伤变量进行粘滞规则化,并推导出相应的粘滞规则化数值一致性切线刚度张量,解决了在有限元隐式计算程序中采用含应变软化段本构关系的数值分析由于计算困难而提前终止的问题。开发包含数值积分算法的用户材料子程序UMAT,并嵌于有限元程序Abaqus v6.14中。通过对力学行为展现显著塑性效应的AS4/3501-6V型开口复合材料层合板的渐进失效分析,验证本文提出的材料本构模型的有效性。结果显示,预测结果与已报道的试验结果吻合良好,并且预测精度高于其他已有弹性损伤模型。表明已建立的弹塑性损伤本构模型能够准确预测力学行为,展现显著塑性效应的复合材料层合板的破坏荷载,为其构件和结构设计提供一种有效的分析方法。     

考虑体积塑性应变的岩石损伤本构模型研究

采用Eshelby等效夹杂方法建立岩石弹塑性损伤本构模型是一种有效方法,但相关文献目前还极少。在连续介质损伤力学框架内利用细观力学的Eshelby等效夹杂方法建立了考虑损伤相塑性体积变形的岩石的Helmholtz自由比能函数。利用连续介质损伤力学方法推导出了考虑损伤相塑性变形的岩石损伤本构关系,给出了损伤演化方程和塑性应变发展方程。并通过数值模拟证实该模型能够反映岩石体积塑性应变、损伤的变化规律和损伤部分不能承受拉应力等力学特性。     

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

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

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

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

三轴围压下砂浆弹塑性损伤变形过程的细观力学分析

基于对泛函势和Cauchy-Born准则,抽象出弹簧束构元和体积构元,组集两种构元的力学响应,给出了材料的弹性损伤本构关系;考虑滑移作为主要的塑性变形机制,提出了滑移构元,给出了材料的塑性本构关系;利用变形分解机制,得到了由三种构元共同描述的弹塑性损伤本构关系。阐述了在给定应变条件下弹塑性损伤本构关系的计算迭代流程。利用单轴拉伸算例详细阐述了模型参数的标定过程。对有围压作用下砂浆材料的压缩行为进行了模拟,从材料细观变形角度解释了随着围压增加,材料承载能力增加的现象。模型预测结果与试验结果符合较好,初步验证了模型具有处理非比例加载问题的能力。     

基于不可逆热力学的Ni-Ti合金动态本构模型及其有限元实现

为了准确描述Ni-Ti形状记忆合金在高应变率下的动态压缩力学行为,基于不可逆热力学理论框架假定了两个内变量表征Ni-Ti合金应力诱发马氏体相变与塑性屈服的不可逆变形过程,分别推导了马氏体相变与塑性屈服演化规律的主控方程,构建了Ni-Ti合金的三维动态本构模型。根据材料单轴动态压缩实验的应力-应变曲线并采用最小二乘法对本构参数进行了优化识别,然后采用应力补偿更新算法,通过隐式用户子程序接口UMAT将动态本构模型嵌入ABAQUS有限元软件,实现了Ni-Ti合金在高应变率下动态压缩力学行为的数值模拟。通过比对发现,模拟结果与实验数据吻合良好,验证了动态本构模型与UMAT子程序的准确性。本工作为Ni-Ti合金在高速冲击、切削等极端条件下的工程应用奠定了基础。     

梯度依赖的混凝土弹塑性非局部损伤的本构模型

以连续介质力学和不可逆热力学为基础,提出了新的混凝土弹塑性损伤本构模型.在该模型中采用了塑性应变εp、各向同性损伤标量D以及应变梯度(△)ε作为内变量,其中应变梯度反应了损伤的非局部性质,本模型是梯度依赖的非局部损伤模型,严格满足热力学的基本方程.应用梯度依赖弹塑性非局部损伤本构模型得出的结果比已有的混凝土损伤塑性模型更为合理.     

非热弹性马氏体相变的细观本构模型

为了更清楚地认识铁基合金经受非热弹性马氏体相变的本征特性,在细观尺度对非热弹性马氏体相变进行了研究.基于马氏体相变晶体学和内变量本构理论建立了非热弹性马氏体相变的细观本构模型.该模型采用微区相变应变、奥氏体及马氏体的塑性应变表征宏观的非弹性响应,把奥氏体和马氏体变体的等效塑性应变率和体积分数变化率作为内变量描述微观结构变化.模型采用J2流动理论描述微区塑性流动,与采用晶体塑性的描述方法相比模型更简单,且更适用于工程计算.单晶奥氏体单变体简单剪切的模拟结果表明:随着应变的增加,先发生奥氏体塑性变形,进而发生相变,马氏体体积分数与应变呈线性关系;温度较低时易发生马氏体相变并使得材料的强度提高.     

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

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

混凝土正交各向异性动态损伤本构模型研究

为建立混凝土的正交各向异性动态损伤本构,首先采用动力放大系数的形式考虑材料的应变率效应,然后基于Sidoroff能量等价原理,建立动态条件下单元体的损伤刚度矩阵。同时在损伤演化模型中,采用Mazars损伤模型描述主轴方向的损伤变量。另外采用适用于正交各向异性的Hoffman屈服破坏准则,并考虑损伤及动力放大系数对强度的修正。以建筑结构受爆破震动荷载作用的实验资料,对本模型进行了验证     

Failure analysis method of concrete arch dam based on elastic strain energy criterion

The arch dam is a type of massive water-retaining structure made of concrete. The overall failure mechanism and corresponding analysis criterion are key issues of concern in the dam engineering field. In this paper, the energy evolution of arch dams in the failure process is studied first, which can be decomposed as energy dissipation accompanied by concrete damage and elastic strain energy absorption and release during elastic deformation. An evaluation criterion for failure analysis of concrete arch dams is then established based on elastic strain energy. An orthotropic damage constitutive model for dam concrete is then proposed along with its numerical simulation method, which is established for structural failure analysis. Numerical simulations show that the elastic strain energy in elements increases with increasing overload safety coefficient and finally converges to the concrete material surface energy, at which time the locally plastic damage area develops rapidly and finally leads to cracking failure of structures. The proposed failure analysis criterion for concrete dams under integrated loads is suitable for analyzing dam instability failure, which has great operability and value in engineering applications in the future.     

一种新的混凝土各向异性弹塑性损伤本构模型及其数值实施

该文的目的是建立一种新的、相对简单的混凝土各向异性塑性损伤本构模型,以方便的模拟混凝土结构的破坏行为。为了更好地描述混凝土在拉、压荷载作用下的不同损伤机制,建立了拉、压不同的两种损伤演化方程,用于确定各向异性的拉、压损伤变量。另外,根据应变等效假设,假定有效构型和损伤构型的应变相等,该方法不仅大大简化了模型的推导过程,而且可方便的通过解耦算法进行有效应力和损伤及名义应力的计算,也即塑性部分计算可通过现有的隐式算法实现,损伤部分及名义应力的计算则可通过较为简便的显式算法实现,从而可大大提高计算效率。模型结果与试验结果的对比分析表明:该模型能较好地描述混凝土在三维应力状态下的非线性行为;对双边开口四点弯曲梁试件的模拟也表明:该模型能反应混凝土损伤各向异性的特点,计算结果相比ABAQUS软件自带的混凝土损伤塑性本构模型(CDP模型)更符合实际情况,计算效率也更高。     

Non‐linear seismic behaviour of concrete gravity dams using coupled finite element–boundary element technique

In this paper, the seismic response of concrete gravity dams is presented using the concept of Continuum Damage Mechanics (CDM) and adopting the hybrid Finite Element–Boundary Element technique (FE–BE). The finite element method is used for discretization of the near field and the boundary element method is employed to model the semi‐infinite far field. Because of the non‐linear nature of the discretizied equations of motion modified Newton–Raphson approach has been used at each time step to linearize them. Damage evolution based on tensile principal strain using mesh‐dependent hardening modulus technique is adopted to ensure the mesh objectivity and to calculate the accumulated damage. The methodology employed is shown to be computationally efficient and consistent in its treatment of both damage growth and damage propagation in gravity dams. Other important features considered in the analysis are: (1) realistic damage modelling for concrete that allows isotropic as well as anisotropic damage state and exhibits stiffness recovery upon load reversals. (2) softening initiation and strain softening criteria for concrete, and (3) proper modelling of semi‐infinite foundation using FE–BE method that allows to consider dam–foundation interaction analysis. As an application of the proposed formulation a gravity dam has been analysed and the results are compared with different foundation stiffnesses. The results of the analysis indicate the importance of including rock foundation in the seismic analysis of dams. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Koyna混凝土重力坝的塑性地震损伤响应分析

考虑了混凝土拉压异性损伤变量，根据相关流动法则推导了混凝土的损伤本构方程。通过有限元方法模拟坝体在地震激励下的非线性损伤演化过程，数值计算过程中考虑了库坝动力耦合及混凝土后继屈服的强化效应，采用四参数Ottosen屈服准则及相应的损伤型本构，算例验证了Koyna混凝土重力坝在竖向和横向实测地震激励下的非线性损伤响应特性。     

水工混凝土弥散型裂缝数值模型中开裂判据的研究

重大水利工程中,混凝土结构多处于高水压环境,坝体或多或少存在着裂缝,然而大部分开裂的混凝土坝仍能够安全运行,从而引申出是否可以允许混凝土结构出现裂缝,多大的裂缝才不至于引起破坏,以及用怎样一个合理的指标来对混凝土开裂程度进行描述评判的问题。该文以等效塑性应变表征混凝土弥散型裂缝模型数值模拟中的开裂破坏,通过混凝土单轴拉伸数值模拟,并根据试验资料及相关文献,建立了等效塑性应变与主拉应变、裂缝宽度之间的关系,得到对应于“有害裂缝”的等效塑性应变值。将方法和成果运用到混凝土重力坝的开裂研究中,具有较好的实际工程意义。     

考虑塑性徐变的高混凝土坝实测应变转换应力探讨

将应变计组实测应变转换为应力与混凝土的应力状态有关。针对高混凝土坝应力作用水平较高,当混凝土的应力超过一定的限度,混凝土将进入塑性徐变阶段,如果仍基于弹性徐变体的应力-应变关系进行实测应变的应力转换,获得的应力与实际情况不符。该文假设混凝土在高应力作用下将产生塑性流动,根据P.Perzyna假设计算黏塑性应变率,首先推导了最大拉应力屈服准则和Hsieh-Ting-Chen屈服准则的黏塑性应变率计算公式,接着推导了考虑塑性徐变的实测应变转换应力的计算公式,进而探讨了考虑塑性徐变的高混凝土坝实测应变转换为应力。实例分析表明:由于将实测应变转换为应力采用增量法进行计算,在转换过程中,某阶段的应力失真,必然导致后续转换应力的真实性,而考虑塑性徐变的实测应变转换的应力更符合实际情况。     

On analysis of the elasto-viscoplastic response of single crystals with anisotropic damage: constitutive modelling and computational aspects

Based on the concept of continuum damage mechanics, an anisotropic damage model for single crystals under the theory of crystal plasticity is presented. Damage and inelastic deformations are incorporated in the proposed model which is developed within the framework of thermodynamics with internal state variables. The dependence of the plastic anisotropy on the damage evolution has been considered. The anisotropic damage is characterized kinematically here through a second-order damage tensor which is physically based. The proposed model can successfully describe the interaction between the evolution of micro-structure of single crystals such as lattice orientation and the hardness development of each slip system and the process of material degradation. The Newton–Raphson iterative scheme is used to integrate the constitutive equations that work directly with the evolution equations for the elastic deformation gradient. The consistent algorithmic tangent stiffness for the present algorithm is formulated. The prescribed algorithm together with the consistent algorithmic tangent stiffness has been implemented into the ABAQUS finite element code by using user subroutine. Using the loading processes with homogeneous deformations and simulation of the classical tensile test of a notched bar illustrate the basic aspects of the model described. Numerical simulations show the validation and performance of the present model and algorithm. Copyright © 2004 John Wiley & Sons, Ltd.     

Seismic analysis of arch dams—a continuum damage mechanics approach

In this paper, the non‐linear seismic response of arch dams is presented using the concept of Continuum Damage Mechanics (CDM). The analysis is performed using the finite element technique and appropriate non‐linear material and damage models in conjunction with the α‐algorithm for time marching. Because of the non‐linear nature of the discretizied equations of motion, modified Newton–Raphson approach has been used at each time step. Damage evolution based on tensile principal strain using mesh‐dependent hardening modulus technique is adopted to ensure the mesh objectivity and to calculate the accumulated damage. The methodology employed is shown to be computationally efficient and consistent in its treatment of both damage growth and damage propagation. As an application of the proposed formulation, a double curvature arch dam has been analysed and the results are compared with the solutions from linear analysis and it is shown that the structural response of arch dams varies significantly in terms of damage evolution. Copyright © 1999 John Wiley & Sons, Ltd.