脆性固体非弹性和软化性能的损伤本构模拟

本文概要地介绍了作者近年来在脆性固体粘塑性损伤本构模型方面的工作。首先，通过假设损伤率依赖于粘塑性应变率和应力变化率，提出了一个考虑损伤演化及其对材料力学性能的影响的粘塑性损伤本构模型。     

基于OTTOSEN准则的混凝土粘塑性力学模型

研究了混凝土的弹．粘塑性动态本构关系,提出一个应变率相关的非线性混凝土模型．模型改进Bicanic的混凝土塑性间断面变化规律和Ottosen的四参数屈服准则,建立冲击荷载下的混凝土本构方程,可以应用于地震和爆炸作用下混凝土材料响应的研究．根据混凝土动态实验结果对应变率和材料强度的关系提出合理假设,考虑以往被忽略的材料峰值应变和泊松比随应变率的变化．模型包含静水压力参数和动态塑性损伤因子,可以有效地反映混凝土的动态力学行为,为其动态问题的研究提供有益的思路和有效的工具．     

含损伤材料的热粘塑性本构关系和柱壳破裂研究

以含内变量的本构关系理论为基础 ,结合材料损伤演化方程 ,并考虑了温度和损伤对材料参数的影响 ,得到了增量形式的热粘塑性本构关系的普适显式表达式。然后使用Bodner幂函数型粘塑性模型 ,具体推导了其增量形式的热粘塑性本构方程。接着结合在实践中有重要意义的内部爆炸载荷作用下的柱壳破裂问题 ,建立了含损伤热粘塑性柱壳破裂问题的完备方程组 ,使用有限差分方法 ,完成了对问题的数值模拟 ,并对结果进行了分析。计算结果与实验结果符合良好。     

耦合蠕变损伤的Chaboche粘塑性模型的研究

Chaboche粘塑性统一本构模型由于没有包含损伤变量,在应力恒定时只能得到不变的塑性应变速率,因此不能描述蠕变第三阶段的加速过程.该文按照Lemaitre有效应力的概念,采用Kachanov损伤演化方程,推导了耦合损伤的Chaboche粘塑性流动方程和硬化方程,并链接到有限元软件ANSYS中,将之用于高温合金钢P91的蠕变寿命计算,结果表明耦合损伤模型的模拟值与实验数据基本吻合;为该模型描述蠕变损伤和疲劳交互作用奠定了基础.     

混凝土黏塑性动力损伤本构关系

从静力弹塑性损伤本构关系的基本框架出发, 综合考虑塑性应变与损伤 演化的率敏感性, 建立了能够较为全面地描述混凝土在动力加载条件下非线性性能的混凝土 黏塑性动力损伤本构模型. 为了考虑塑性应变的率敏感性, 基于Perzyna理论推导了有效应 力空间黏塑性力学基本公式, 采用改进的Perzyna型动力演化方程, 将损伤静力演化方程推 广到动力加载情形. 基于并联弹簧模型, 从概率论的角度推导给出了一维损伤静力演化方程, 并基于能量等效应变的基本概念将其推广到多维损伤演化. 利用数值模拟, 计算得到了 混凝土在不同应变率下的应力应变全曲线, 同时得到了一维动力提高因子和二维动力强度包 络图, 数值结果与试验结果的对比表明了该模型的有效性.     

描述大应变率范围下材料响应的粘塑性本构模型

以位错动力学理论中的Ｏｒｗａｎ和Ｇｉｌｍａｎ关系为基础建立描述率相关材料非弹性响应的基本方程，选择材料准静态实验的单轴响应作为强化演化的规律，并考虑应变率敏感程度随变形产生变化的特性，建立了适用于大应变率范围内率相关材料的统一型粘塑性本构模型。对铝１１００－０在应变率范围１０－５～１０４ｓ－１内产生的有限塑性应变的单轴响应进行了理论预测，与Ｋｈａｎ和Ｈｕａｎｇ［１］的实验数据及模型预测结果进行了比较，结果表明本文模型具有较高的预测精度，在高应变率和较大应变下不容忽视率敏感参数随变形的变化。     

粘塑性薄壁管中复合应力波的传播特性研究

以本构关系一般理论为基础,导出了计及材料功硬化效应和应变率硬化效应的粘塑性薄壁管的本构关系及管中复合应力波的控制方程,应用有限差分方法研究了在压扭复合冲击载荷作用下粘塑性薄壁管中复合应力波的传播特性与演化规律,分析了复合应力波的耦合效应以及薄壁管中粘塑性参数和功硬化效应对复合应力波传播与演化规律的影响,并对有关现象进行了分析和解释。     

非比例循环塑性和循环粘塑性本构描述的某些新进展

金属材料循环塑性本构方程和循环粘塑性本构方程是固体力学中近１０多年来一个十分重要的领域。本文评述了金属材料非比例循环塑性界限面本构理论、内时理论和循环粘塑性本构理论及其某些进展，对某些模型中非比例度定义，材料在复杂应变幅值历史、非比例循环加载历史以及其它历史下的强化规则和流动规则进行了分析与评价，在此基础上对循环本构理论的发展趋势提出自己的看法。     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

幂硬化粘塑性材料动态扩展裂纹尖端场的渐近解

采用弹性－粘塑性本构模型，对幂硬化粘塑性介质中反平面剪切动态扩展裂纹尖端的应力，应变场进行了渐近分析，给出了反平面剪切动态扩展纹尖端场的渐进方程。分析结果表明，在裂纹法端应力具有（ｌｎＲ／ｒ）１／ｎ－１的奇异性，应变具有（ｌｎＲ／Ｒｎ／ｎ－１的奇异性。从而提示了幂硬化粘塑材料反平面剪动态扩展裂纹尖端场的渐近行为。     

Application of corotational rates of the logarithmic strain in constitutive modeling of hardening materials at finite deformations

In this paper a finite deformation constitutive model for rigid plastic hardening materials based on the logarithmic strain tensor is introduced. The flow rule of this constitutive model relates the corotational rate of the logarithmic strain tensor to the difference of the deviatoric Cauchy stress and the back stress tensors. The evolution equation for the kinematic hardening of this model relates the corotational rate of the back stress tensor to the corotational rate of the logarithmic strain tensor. Using Jaumann, Green–Naghdi, Eulerian and logarithmic corotational rates in the proposed constitutive model, stress–strain responses and subsequent yield surfaces are determined for rigid plastic kinematic and isotropic hardening materials in the simple shear problem at finite deformations.     

高周疲劳的损伤-硬化模型

两级循环加载条件下,材料的剩余寿命强烈地依赖于加载历史。究其原因,不同加载历史将引起材料的微结构发生不同的变化,使得材料的硬化效果和变形行为表现出明显的差异,从而影响了损伤的演化过程。本文引入硬化状态变量来表征加载历史对疲劳损伤演化过程的影响。通过对两级循环加载下损伤演化规律和剩余寿命的研究,认为在两级（或多级）加载条件下,材料的损伤演化和剩余寿命强烈地依赖于加载历史造成的损伤和硬 化状态。     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

Recovery experimental techniques of tensile impact

Based on loading-unloading test, tensile impact recovery experimental techniques have been developed to obtain the isothermal stress-strain curves of materials under high strain rates. The thermal softening effect can be decoupled by comparing the isothermal stress-strain curves with the adiabatic stress-strain curves at the same strain rate. In the present paper, recovery experiments of brass have been carried out on a self-designed rotating disk tensile impact apparatus. According to the parabolic strain hardening power-law thermo-viscoplastic constitutive model, strain hardening parameter, strain rates strengthening parameter and thermal softening synthetical parameter have been decoupled from experimental results. Furthermore, from these parameters, one can determine the theoretical isothermal curves and adiabatic curves at high strain rates well-coinciding the experimental results respectively. It indicates that the recovery experimental techniques of tensile impact are effective and reliable and are important means for the study of thermo-mechanical coupling. The experimental results also reveals that brass is a typical thermo-viscoplastic material. The project supported by the National Natural Science Foundation of China     

Constitutive modeling of temperature and strain rate dependent elastoplastic hardening materials using a corotational rate associated with the plastic deformation

In this paper, a constitutive model with a temperature and strain rate dependent flow stress (Bergstrom hardening rule) and modified Armstrong-Frederick kinematic evolution equation for elastoplastic hardening materials is introduced. Based on the multiplicative decomposition of the deformation gradient,new kinematic relations for the elastic and plastic left stretch tensors as well as the plastic deformation-dependent spin tensor are proposed. Also, a closed-form solution has been obtained for the elastic and plastic left stretch tensors for the simple shear problem.To evaluate model validity, results are compared with known experimental data for SUS 304 stainless steel, which shows a good agreement with the results of the proposed theoretical model.Finally, the stress-deformation curve, as predicted by the model, is plotted for the simple shear problem at room and elevated temperatures using the same material properties for AA5754-O aluminium alloy.     

Predictions of bifurcation and instabilities during dynamic extension

Dynamic bifurcation and flow instabilities of cylindrical bars, made of an incompressible strain hardening plastic material, are investigated. A Lagrangian linear perturbation analysis is performed to obtain a fourth order partial differential equation which governs the evolution of the perturbation. The analysis shows that inertia slows down the growth of long wavelengths while bidimensional effects conjugated to strain hardening extinct short wavelengths. The present approach is applied successfully to the analysis of bifurcation and instabilities in (i) a rectangular block during plane strain extension, (ii) a circular bar during uniaxial extension. New results are obtained in the case of rate independent materials and a synthetical point of view is obtained for rate dependent behaviors.     

A constitutive theory for creep behavior of initially isotropic materials sustaining unilateral damage

The goal of the present work is to modify structure of the creep constitutive equations existing in the literature, and simultaneously to incorporate both damage induced anisotropy and unilateral damage into the constitutive model. The proposed nonlinear-tensor constitutive equation for creep together with the damage evolution equation take into account the secondary and tertiary creep of the initially isotropic materials. The material parameters of the model are determined using basic experiments. It is shown that the creep model is capable of describing available experimental data for the lateral creep responses under uniaxial compression.     

Thermodynamic based model for the evolution equation of the backstress in cyclic plasticity

A nonlinear kinematic hardening rule is developed here within the framework of thermodynamic principles. The derived kinematic hardening evolution equation has three distinct terms: two strain hardening terms and a dynamic recovery term that operates at all times. The proposed hardening rule, which is referred in this paper as the FAPC (Fredrick and Armstrong–Phillips–Chaboche) kinematic hardening rule, shows a combined form of the Frederick and Armstrong backstress evolution equation, Phillips evolution equation, and Chaboche series rule. A new term is incorporated into the Frederick and Armstrong evolution equation that appears to have agreement with the experimental observations that show the motion of the center of the yield surface in the stress space is directed between the gradient to the surface at the stress point and the stress rate direction at that point. The model is further modified in order to simulate nonproportional cyclic hardening by proposing a measure representing the degree of nonproportionality of loading. This measure represents the topology of the incremental stress path. Numerically, it represents the angle between the current stress increment and the previous stress increment, which is interpreted through the material constants of the kinematic hardening evolution equation. This new kinematic hardening rule is incorporated in a material constitutive model based on the von Mises plasticity type and the Chaboche isotropic hardening type. Numerical integration of the incremental elasto-plastic constitutive equations is based on a simple semi-implicit return-mapping algorithm and the full Newton–Raphson iterative method is used to solve the resulting nonlinear equations. Experimental simulations are conducted for proportional and non-proportional cyclic loadings. The model shows good correlation with the experimental results.     

考虑微缺陷形状的细片层状珠光体钢的损伤本构模型

基于非经典塑性理论和连续介质损伤力学,利用在一个特殊坐标系下基于椭球形孔洞模型得到的可考虑孔洞形状变化混合强化材料的损伤演化率得到了铁素体相的损伤本构方程,通过混合物理论利用铁素体和渗碳体相各自本构关系并考虑其几何特征得到了珠光体团的损伤本构模型。进而采用Hill自洽方法,得到了珠光体材料的宏观损伤本构描述,发展了相应的数值方法与程序。讨论了孔洞形状对材料损伤的影响,并对典型珠光体双相材料BS11在非对称循环加载史下的弹塑性响应特性进行了分析,得到了与实验较为一致的结果。