Numerical analysis and elastic–plastic deformation behavior of anisotropically damaged solids

The present paper is concerned with the numerical modelling of the large elastic–plastic deformation behavior and localization prediction of ductile metals which are sensitive to hydrostatic stress and anisotropically damaged. The model is based on a generalized macroscopic theory within the framework of nonlinear continuum damage mechanics. The formulation relies on a multiplicative decomposition of the metric transformation tensor into elastic and damaged-plastic parts. Furthermore, undamaged configurations are introduced which are related to the damaged configurations via associated metric transformations which allow for the interpretation as damage tensors. Strain rates are shown to be additively decomposed into elastic, plastic and damage strain rate tensors. Moreover, based on the standard dissipative material approach the constitutive framework is completed by different stress tensors, a yield criterion and a separate damage condition as well as corresponding potential functions. The evolution laws for plastic and damage strain rates are discussed in some detail. Estimates of the stress and strain histories are obtained via an explicit integration procedure which employs an inelastic (damage-plastic) predictor followed by an elastic corrector step. Numerical simulations of the elastic–plastic deformation behavior of damaged solids demonstrate the efficiency of the formulation. A variety of large strain elastic–plastic-damage problems including severe localization is presented, and the influence of different model parameters on the deformation and localization prediction of ductile metals is discussed.     

A crystal plasticity model that incorporates stresses and strains due to slip gradients

This work is concerned with incorporating the kinematic and stress effects of excess dislocations in a constitutive model for the elastoplastic behavior of crystalline materials. The foundation of the model is a three term multiplicative decomposition of the deformation gradient in which the two classical terms of plastic and elastic deformation are included along with an additional term for long range strain due to the collective effects of excess dislocations. The long range strain is obtained from an assumed density of Volterra edge dislocations and is directly related to gradients in slip. A new material parameter emerges which is the size the region about a continuum point that contributes to long range strains.Using Hookean elasticity, the stress at a point is linearly related to the sum of the elastic plus the long range strain fields. However, the driving force for slip is postulated to be due only to the elastic stress so that the long range stress is a back stress in the constitutive relationship for plastic deformation. A consistent balance of the total deformation rate with the three proposed mechanisms of deformation leads to a set of differential equations that can be solved for the elastic stress, rotation and pressure which then implicitly defines the material state and equilibrium stress. Results from the simulation of a tapered tensile specimen demonstrate that the constitutive model exhibits isotropic and kinematic type hardening effects as well as changes in the pattern of plastic deformation and necking when compared to a material without slip gradient effects.     

A theory of large-strain isotropic thermoplasticity based on metric transformation tensors

Summary A formulation of isotropic thermoplasticity for arbitrary large elastic and plastic strains is presented. The underlying concept is the introduction of a metric transformation tensor which maps a locally defined six-dimensional plastic metric onto the metric of the current configuration. This mixed-variant tensor field provides a basis for the definition of a local isotropic hyperelastic stress response in the thermoplastic solid. Following this fundamental assumption, we derive a consistent internal variable formulation of thermoplasticity in a Lagrangian as well as a Eulerian geometric setting. On the numerical side, we discuss in detail an objective integration algorithm for the mixed-variant plastic flow rule. The special feature here is a new representation of the stress return and the algorithmic elastoplastic moduli in the eigenvalue space of the Eulerian plastic metric for plane problems. Furthermore, an algorithm for the solution of the coupled problem is formulated based on an operator split of the global field equations of thermoplasticity. The paper concludes with two representative numerical simulations of thermoplastic deformation processes.     

Wedge indentation into elastic-plastic single crystals, 1: Asymptotic fields for nearly-flat wedge

Asymptotic stress and deformation fields under the contact point singularities of a nearly-flat wedge indenter and of a flat punch are derived for elastic ideally-plastic single crystals with three effective in-plane slip systems that admit a plane strain deformation state. Face-centered cubic (FCC), body-centered cubic (BCC), and hexagonal-close packed (HCP) crystals are considered. The asymptotic fields for the flat punch are analogous to those at the tip of a stationary crack, so a potential solution is that the deformation field consists entirely of angular constant stress plastic sectors separated by rays of plastic deformation across which stresses change discontinuously. The asymptotic fields for a nearly-flat wedge indenter are analogous to those of a quasistatically growing crack tip fields in that stress discontinuities can not exist across sector boundaries. Hence, the asymptotic fields under the contact point singularities of a nearly-flat wedge indenter are significantly different than those under a flat punch. A family of solutions is derived that consists entirely of elastically deforming angular sectors separated by rays of plastic deformation across which the stress state is continuous. Such a solution can be found for FCC and BCC crystals, but it is shown that the asymptotic fields for HCP crystals must include at least one angular constant stress plastic sector. The structure of such fields is important because they play a significant role in the establishment of the overall fields under a wedge indenter in a single crystal. Numerical simulations—discussed in detail in a companion paper—of the stress and deformation fields under the contact point singularity of a wedge indenter for a FCC crystal possess the salient features of the analytical solution.     

A two-dimensional dislocation dynamics model of the plastic deformation of polycrystalline metals

Two-dimensional dislocation dynamics (2D-DD) simulations under fully periodic boundary conditions are employed to study the relation between microstructure and strength of a material. The material is modeled as an elastic continuum that contains a defect microstructure consisting of a preexisting dislocation population, dislocation sources, and grain boundaries. The mechanical response of such a material is tested by uniaxially loading it up to a certain stress and allowing it to relax until the strain rate falls below a threshold. The total plastic strain obtained for a certain stress level yields the quasi-static stress-strain curve of the material. Besides assuming Frank-Read-like dislocation sources, we also investigate the influence of a pre-existing dislocation density on the flow stress of the model material. Our results show that - despite its inherent simplifications - the 2D-DD model yields material behavior that is consistent with the classical theories of Taylor and Hall-Petch. Consequently, if set up in a proper way, these models are suited to study plastic deformation of polycrystalline materials.     

Micromechanical determination of two-phase plasticity

A micromechanical theory is developed to predict the elastoplastic behavior of a two-phase alloy. Taking crystallographic slip to be the mechanism of plastic deformation, this theory also considers stress redistribution due to elastic and plastic heterogeneity in both phases. The corresponding self-consistent relation for two-phase plasticity was derived combining the spirit of Hill, Hutchinson, and Berveiller & Zaoui. It is found upon applications that both elastic and rigid particles may have a profound effect on the hardening behavior of two-phase systems. When applied to austenite-ferrite stainless steels, the theory also provides reasonable estimates as compared to experiments. The fictitious kink point commonly associated tith the continuum models is seen to be absent due to the gradual yielding of the constituent grains.     

考虑塑性和刚度不匹配的焊接残余应力估算

现有残余应力计算方法未能考虑材料塑性变形和焊接接头刚度不匹配的影响,使得焊接残余应力计算结果和实际残余应力存在较大偏差.在2219-T87铝合金钨极氩弧焊焊接头残余应力测试基础上,提出一种基于非线性有限元和材料弹性模量分区的残余应力—释放应变曲线的残余应力计算方法,研究了材料塑性变形和接头刚度不匹配对焊接残余应力计算的影响.结果表明,焊接接头中非均质材料塑性不匹配可以引起对于残余应力计算的较大误差;材料塑性变形对残余应力的影响大于接头刚度不匹配对残余应力的影响.所提出方法修正了传统方法在焊接接头的残余应力计算中由于未考虑接头非均质材料塑性不匹配而引起的误差.     

Dislocations and disclinations: continuously distributed defects in elasto-plastic crystalline materials

The paper deals with elasto-plastic models for crystalline materials with defects, dislocations coupled with disclinations. The behaviour of the material is described with respect to an anholonomic configuration, endowed with a non-Riemannian geometric structure. The geometry of the material structure is generated by the plastic distortion, which is an incompatible second-order tensor, and by the so-called plastic connection which is metric compatible, with respect to the metric tensor associated with the plastic distortion. The free energy function is dependent on the second-order elastic deformation and on the state of defects. The tensorial measure of the defects is considered to be the Cartan torsion of the plastic connection and the disclination tensor. When we restrict to small elastic and plastic distortions, the measures of the incompatibility as well as the dislocation densities reduced to the classical ones in the linear elasticity. The constitutive equations for macroforces and the evolution equations for the plastic distortion and disclination tensor are provided to be compatible with the free energy imbalance principle.     

Constitutive relations for isotropic or kinematic hardening at finite elastic–plastic deformations

The rate-type constitutive relations of rate-independent metals with isotropic or kinematic hardening at finite elastic–plastic deformations were presented through a phenomenological approach. This approach includes the decomposition of finite deformation into elastic and plastic parts, which is different from both the elastic–plastic additive decomposition of deformation rate and Lee’s elastic–plastic multiplicative decomposition of deformation gradient. The objectivity of the constitutive relations was dealt with in integrating the constitutive equations. A new objective derivative of back stress was proposed for kinematic hardening. In addition, the loading criteria were discussed. Finally, the stress for simple shear elastic–plastic deformation was worked out.     

A thermodynamic formulation of finite-deformation elastoplasticity with hardening based on the concept of material isomorphism

This paper presents a thermodynamic formulation of a model for finite deformation of materials exhibiting elastoplastic material behaviour with non-linear isotropic and kinematic hardening. Central to this formulation is the notion that the form of the elastic constitutive relation be unaffected by the plastic deformation or transformation in the material, as commonly assumed in particular in the context of crystal plasticity. When generalized to the phenomenological context, this implies that the internal variable representing plastic deformation is an elastic material isomorphism. Among other things, this requirement on the plastic deformation leads directly to the standard elastoplastic multiplicative decomposition of the deformation gradient. In addition, a dependence of the plastic part of the free energy on the plastic deformation itself yields a thermodynamic form for the centre of the elastic range of the material, i.e. the back stress. Finally, we show how this approach can be applied to formulate thermodynamic forms for linear, and non-linear Armstrong-Frederick, kinematic hardening models.     

激光衍射谱测量金属塑性变形的方法研究

本文讨论了弹性和塑性变形对金属表面衍射谱的不同影响,根据激光衍射谱的变化测量金属塑性变形的表面光反射率方法,对45~#钢材进行了测量,结果表明这种测试方法是可行的。同时测得了金属弹塑性变形与谱面衍射光强、逸出光强的关系曲线;对比了不同光路布置后提出了最佳方案;发现光强比与对数应变曲线的拐点可作为材料是否发生塑性变形的判据。     

模型复合材料弹塑性界面应力分析

由纤维增强弹塑性基体所产生的界面具有弹塑性力学行为。考虑到一般材料的塑性变形都遵循幂硬化规律,对模型复合材料的界面进行弹性和应变硬化状态下的变形规律及其应力分析。以纤维拔出试验为研究模型,将界面分成弹性区和塑性区。利用界面应力剪滞理论,分别建立弹性区和塑性区的界面力学基本方程。选择适当的位移函数满足基本方程及埋入纤维的边界条件,再按位移函数求出弹性区和塑性区的界面剪应力。推导出平均界面剪应力与纤维     

On stress rate and plasticity constitutive equations referred to a body-fixed coordinate system

Investigators have viewed the stress rate in two different ways: the material (body-fixed) point of view and the Eulerian point of view. We discuss the Zaremba–Jaumann rate and Oldroyd’s rate from the material viewpoint and apply them to the material formulation of a theory of plasticity for materials undergoing anisotropic plastic deformation. Significant advantages of the material formulation are that the derivation of equations is straight forward, the distortion of yield surface can be easily accounted for, and the issue of self-consistent elastic equation does not arise.     

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性-塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.      

The crack tip zone shielding effect for ductile particle reinforced brittle materials

The crack tip zone shielding effect for the ductile particle reinforced brittle materials is analyzed by using a micromechanics constitutive theory. The theory is developed here to determine the elastoplastic constitutive behavior of the composite. The elastoplastic particles, with isotropic or kinematical hardening, are uniformly dispersed in the brittle elastic matrix. The method proposed is based on the Mori-Tanaka's concept of average stress in the composite. The macroscopic yielding condition and the incremental stress strain relation of the composite during plastic deformation are explicity given in terms of the macroscopioc applied stress and the microstructural parameters of the composite such as the volume fraction and yield stress of ductile particles, elastic constants of the two phases, etc. Finally, the contribution of the plastic deformation in the particles near a crack tip to the toughening of the composite is evaluated. The project supported by National Natural Science Foundation of China     

Suppressed plastic deformation at blunt crack-tips due to strain gradient effects

Large deformation gradients occur near a crack-tip and strain gradient dependent crack-tip deformation and stress fields are expected. Nevertheless, for material length scales much smaller than the scale of the deformation gradients, a conventional elastic–plastic solution is obtained. On the other hand, for significant large material length scales, a conventional elastic solution is obtained. This transition in behaviour is investigated based on a finite strain version of the Fleck–Hutchinson strain gradient plasticity model from 2001. The predictions show that for a wide range of material parameters, the transition from the conventional elastic–plastic to the elastic solution occurs for length scales ranging from 0.001 times the size of the plastic zone to a length scale of the same order of magnitude as the plastic zone.     

A unified expression of elastic–plastic notch stress–strain calculation in bodies subjected to multiaxial cyclic loading

In this paper, a simplified thermodynamics analysis of cyclic plastic deformation is performed in order to establish an energy transition relation for describing the elastic–plastic stress and strain behavior of the notch-tip material element in bodies subjected to multiaxial cyclic loads. Based on the thermodynamics analysis, it is deduced that in the case of elastic–plastic deformation, Neuber’s rule inevitably overestimates the actual stress and strain at the notch tip, while the equivalent strain energy density (ESED) method tends to underestimate the actual notch-tip stress and strain. According to the actual energy conversion occurring in the notch-tip material element during cyclic plastic deformation, a unified expression for estimating the elastic–plastic notch stress–strain responses in bodies subjected to multiaxial cyclic loads is developed, of which Neuber’s rule and the ESED method become two particular cases, i.e. upper and lower bound limits of the notch stress and strain estimations. This expression is verified experimentally under both proportional and non-proportional multiaxial cyclic loads and a good agreement between the calculated and the measured notch strains has been achieved. It is also shown that in the case of multiaxial cyclic loading, the unified expression distinctly improves the accuracy of the notch-tip stress–strain estimations in comparison with Neuber’s rule and the ESED method. The unified expression of the notch stress–strain calculation developed in this paper can thus provide a more logical approximate approach for estimating the elastic–plastic notch-tip stress and strain responses of components subjected to lengthy multiaxial cyclic loading histories for local strain approach-based fatigue-crack-initiation life prediction.