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

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

Micromorphic continuum. Part II: Finite deformation plasticity coupled with damage

It is demonstrated how a micromorphic plasticity theory may be formulated on the basis of multiplicative decompositions of the macro- and microdeformation gradient tensor, respectively. The theory exhibits non-linear isotropic and non-linear kinematic hardening. The yield function is expressed in terms of Mandel stress and double stress tensors, appropriately defined for micromorphic continua. Flow rules are derived from the postulate of Il’iushin and represent generalized normality conditions. Evolution equations for isotropic and kinematic hardening are introduced as sufficient conditions for the validity of the second law of thermodynamics in every admissible process. Finally, it is sketched how isotropic damage effects may be incorporated in the theory. This is done for the concept of effective stress combined with the hypothesis of strain equivalence.     

Nonlocal gradient-dependent modeling of plasticity with anisotropic hardening

This work addresses the formulation of the thermodynamics of nonlocal plasticity using the gradient theory. The formulation is based on the nonlocality energy residual introduced by Eringen and Edelen (1972). Gradients are introduced for those variables associated with isotropic and kinematic hardening. The formulation applies to small strain gradient plasticity and makes use of the evanescent memory model for kinematic hardening. This is accomplished using the kinematic flux evolution as developed by Zbib and Aifantis (1988). Therefore, the present theory is a four nonlocal parameter-based theory that accounts for the influence of large variations in the plastic strain, accumulated plastic strain, accumulated plastic strain gradients, and the micromechanical evolution of the kinematic flux. Using the principle of virtual power and the laws of thermodynamics, thermodynamically-consistent equations are derived for the nonlocal plasticity yield criterion and associated flow rule. The presence of higher-order gradients in the plastic strain is shown to enhance a corresponding history variable which arises from the accumulation of the plastic strain gradients. Furthermore, anisotropy is introduced by plastic strain gradients in the form of kinematic hardening. Plastic strain gradients can be attributed to the net Burgers vector, while gradients in the accumulation of plastic strain are responsible for the introduction of isotropic hardening. The equilibrium between internal Cauchy stress and the microstresses conjugate to the higher-order gradients frames the yield criterion, which is obtained from the principle of virtual power. Microscopic boundary conditions, associated with plastic flow, are introduced to supplement the macroscopic boundary conditions of classical plasticity. The nonlocal formulation developed here preserves the classical assumption of local plasticity, wherein plastic flow direction is governed by the deviatoric Cauchy stress. The theory is applied to the problems of thin films on both soft and hard substrates. Numerical solutions are presented for bi-axial tension and simple shear loading of thin films on substrates.     

A nonlocal strain gradient plasticity theory for finite deformations

Strain gradient plasticity for finite deformations is addressed within the framework of nonlocal continuum thermodynamics, featured by the concepts of (nonlocality) energy residual and globally simple material. The plastic strain gradient is assumed to be physically meaningful in the domain of particle isoclinic configurations (with the director vector triad constant both in space and time), whereas the objective notion of corotational gradient makes it possible to compute the plastic strain gradient in any domain of particle intermediate configurations. A phenomenological elastic–plastic constitutive model is presented, with mixed kinematic/isotropic hardening laws in the form of PDEs and related higher order boundary conditions (including those associated with the moving elastic/plastic boundary). Two fourth-order projection tensor operators, functions of the elastic and plastic strain states, are shown to relate the skew-symmetric parts of the Mandel stress and back stress to the related symmetric parts. Consistent with the thermodynamic restrictions therein derived, the flow laws for rate-independent associative plasticity are formulated in a six-dimensional tensor space in terms of symmetric parts of Mandel stresses and related work-conjugate generalized plastic strain rates. A simple shear problem application is presented for illustrative purposes.     

Cyclic stress-strain behaviour of alloy AA6060 T4, part II: Biaxial experiments and modelling

Laboratory tests have been conducted to investigate the inelastic behaviour of aluminium alloy AA6060 T4 subjected to non-proportional cyclic loading. The results of four tests with variable strain path shapes and strain amplitudes are reported in this paper. The tests were carried out by applying combined axial force and torque to thin-walled tubular specimens, using effective strain amplitudes in the range 0.4–0.8%. Major emphasis has been put on the two important material properties: plastic anisotropy and influence of strain range and strain path shapes on cyclic hardening. A constitutive model for cyclic plasticity is used to predict the stress response of the alloy for the non-proportional strain paths applied in the experiments. The model adopts a quadratic yield function and multi-component non-linear isotropic and kinematic hardening rules to describe plastic anisotropy, the shape of the hysteresis loops and the evolution of cyclic hardening. Good agreement is obtained between the physical and correlated stress response of the alloy.     

A mathematical basis for strain-gradient plasticity theory. Part II: Tensorial plastic multiplier

A phenomenological, flow theory version of gradient plasticity for isotropic and anisotropic solids is constructed along the lines of Gudmundson [Gudmundson, P., 2004. A unified treatment of strain-gradient plasticity. J. Mech. Phys. Solids 52, 1379-1406]. Both energetic and dissipative stresses are considered in order to develop a kinematic hardening theory, which in the absence of gradient terms reduces to conventional J₂ flow theory with kinematic hardening. The dissipative stress measures, work-conjugate to plastic strain and its gradient, satisfy a yield condition with associated plastic flow. The theory includes interfacial terms: elastic energy is stored and plastic work is dissipated at internal interfaces, and a yield surface is postulated for the work-conjugate stress quantities at the interface. Uniqueness and extremum principles are constructed for the solution of boundary value problems, for both the rate-dependent and the rate-independent cases. In the absence of strain gradient and interface effects, the minimum principles reduce to the classical extremum principles for a kinematically hardening elasto-plastic solid. A rigid-hardening version of the theory is also stated and the resulting theory gives rise to an extension to the classical limit load theorems. This has particular appeal as previous trial fields for limit load analysis can be used to generate immediately size-dependent bounds on limit loads.     

一般加载规律的弹塑性本构关系

将有关文献给出一般加载规律一维全量理论的简单模型推广到一般加载规律的一维增量理论,进而推广到一般加载规律的多维增量理论,在此基础上,建立了推导一般加载规律的多维增量理论的本构关系的一种途径。应用这种途径,从应力空间的加载函数和应变空间的加载函数出发,推导了等向强化材料和被加热的等向强化材料的一般加载规律的弹塑性本构关系的两种表示形式。理论和实例均表明,这种途径对等向强化材料、随动强化材料和理想弹塑性材料均适用。     

Plastic deformation modelling of tempered martensite steel block structure by a nonlocal crystal plasticity model

The plastic deformations of tempered martensite steel representative volume elements with different martensite block structures have been investigated by using a nonlocal crystal plasticity model which considers isotropic and kinematic hardening produced by plastic strain gradients. It was found that pronounced strain gradients occur in the grain boundary region even under homogeneous loading. The isotropic hardening of strain gradients strongly influences the global stress–strain diagram while the kinematic hardening of strain gradients influences the local deformation behaviour. It is found that the additional strain gradient hardening is not only dependent on the block width but also on the misorientations or the deformation incompatibilities in adjacent blocks.     

Anisotropic plasticity for sheet metals using the concept of combined isotropic-kinematic hardening

Hill's 1948 anisotropic theory of plasticity (Hill, R., 1948. A theory of yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. London A193, 281–297) is extended to include the concept of combined isotropic-kinematic hardening, and the objective of this paper is to validate the model so that it may be useful for analyses of sheet metal forming. Isotropic hardening and kinematic hardening may be experimentally observed in sheet metals, if yielding is defined by the proportional limit or by a small proof strain. In this paper, a single exponential term is used to describe isotropic hardening and Prager's linear kinematic hardening rule is applied for simplicity. It is shown that this model can satisfactorily describe both the yield stress and the plastic strain ratio, the R-ratio, observed in tension test of specimens cut at various angles measured from the rolling direction of the sheet. Kinematic hardening leads to a gradual change in the direction of the plastic strain increment, as the axial strain increases in the tension test; while in the traditional approach for sheet metal, this direction does not change due to the use of isotropic hardening.     

On the modelling of anisotropic elastic and inelastic material behaviour at large deformation

The purpose of this work is the formulation and discussion of an approach to the modelling of anisotropic elastic and inelastic material behaviour at large deformation. This is done in the framework of a thermodynamic, internal-variable-based formulation for such a behaviour. In particular, the formulation pursued here is based on a model for plastic or inelastic deformation as a transformation of local reference configuration for each material element. This represents a slight generalization of its modelling as an elastic material isomorphism pursued in earlier work, allowing one in particular to incorporate the effects of isotropic continuum damage directly into the formulation. As for the remaining deformation- and stress-like internal variables of the formulation, these are modelled in a fashion formally analogous to so-called structure tensors. On this basis, it is shown in particular that, while neither the Mandel nor back stress is generally so, the stress measure thermodynamically conjugate to the plastic “velocity gradient”, containing the difference of these two stress measures, is always symmetric with respect to the Euclidean metric, i.e., even in the case of classical or induced anisotropic elastic or inelastic material behaviour. Further, in the context of the assumption that the intermediate configuration is materially uniform, it is shown that the stress measure thermodynamically conjugate to the plastic velocity gradient is directly related to the Eshelby stress. Finally, the approach is applied to the formulation of metal plasticity with isotropic kinematic hardening.     

Plastic yielding in cyclically loaded porous materials

The present paper focuses on plastic yielding of cyclically loaded porous materials. Unit cell models are employed to observe the evolution of the yield surface of porous materials under cyclic loading conditions. Non-linear isotropic as well as non-linear kinematic hardening matrix materials are considered. The yield surfaces computed with the unit cell models are compared to predictions of a micro-mechanical porous plasticity model that incorporates hardening. It is found that, in the case of kinematic hardening, the porous plasticity model underestimates the yield strength for larger hydrostatic stresses. An improvement of the model is proposed, so that a reasonable micro-mechanical approach to model porous materials under cyclic loadings is found.     

A three-invariant cap model with isotropic–kinematic hardening rule and associated plasticity for granular materials

In this paper, a three-invariant cap model is developed for the isotropic–kinematic hardening and associated plasticity of granular materials. The model is based on the concepts of elasticity and plasticity theories together with an associated flow rule and a work hardening law for plastic deformations of granulars. The hardening rule is defined by its decomposition into the isotropic and kinematic material functions. The constitutive elasto-plastic matrix and its components are derived by using the definition of yield surface, material functions and non-linear elastic behavior, as function of hardening parameters. The model assessment and procedure for determination of material parameters are described. Finally, the applicability of proposed plasticity model is demonstrated in numerical simulation of several triaxial and confining pressure tests on different granular materials, including: wheat, rape, synthetic granulate and sand.     

Role of plasticity on indentation behavior: Relations between surface and subsurface responses

Surface and subsurface responses during frictionless indentation of elasto-plastic solids are investigated. Cases of monotonic and repeated loading are considered. It is shown that the role of plasticity parameters on indentation behavior cannot be well described in terms of surface response alone. It must be tied with subsurface response. In particular, the difference between the materials exhibiting isotropic hardening and that exhibiting kinematic hardening is less apparent when the force-displacement response and the deformed surface from the two are compared. Yet, it is quite apparent when subsurface response such as plastic strain, residual stress and plastic zone dimension are compared. An attempt has been made to characterize such surface and subsurface responses, for different plastic behavior, and to compare them with estimations obtained from analytical solutions.     

Constitutive modeling of cyclic plasticity deformation of a pure polycrystalline copper

Cyclic plasticity experiments were conducted on a pure polycrystalline copper and the material was found to display significant cyclic hardening and nonproportional hardening. An effort was made to describe the cyclic plasticity behavior of the material. The model is based on the framework using a yield surface together with the Armstrong–Frederick type kinematic hardening rule. No isotropic hardening is considered and the yield stress is assumed to be a constant. The backstress is decomposed into additive parts with each part following the Armstrong–Frederick type hardening rule. A memory surface in the plastic strain space is used to account for the strain range effect. The Tanaka fourth order tensor is used to characterize nonproportional loading. A set of material parameters in the hardening rules are related to the strain memory surface size and they are used to capture the strain range effect and the dependence of cyclic hardening and nonproportional hardening on the loading magnitude. The constitutive model can describe well the transient behavior during cyclic hardening and nonproportional hardening of the polycrystalline copper. Modeling of long-term ratcheting deformation is a difficult task and further investigations are required.     

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.     

A mathematical basis for strain-gradient plasticity theory—Part I: Scalar plastic multiplier

Strain-gradient plasticity theories are reviewed in which some measure of the plastic strain rate is treated as an independent kinematic variable. Dislocation arguments are invoked in order to provide a physical basis for the hardening at interfaces. A phenomenological, flow theory version of gradient plasticity is constructed in which stress measures, work-conjugate to plastic strain and its gradient, satisfy a yield condition. Plastic work is also done at internal interfaces and a yield surface is postulated for the work-conjugate stress quantities at the interface. Thereby, the theory has the potential to account for grain size effects in polycrystals. Both the bulk and interfacial stresses are taken to be dissipative in nature and due attention is paid to ensure that positive plastic work is done. It is shown that the mathematical structure of the elasto-plastic strain-gradient theory has similarities to conventional rigid-plasticity theory. Uniqueness and extremum principles are constructed for the solution of boundary value problems.     

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.     

Variational formulations for the plane strain elastic-plastic problem for materials governed by the von Mises criterion

Path-dependent materials, complying with Drucker's postulate requirements and governed by an internal variable rate plasticity model, are considered. A variational principle for the small strain, rate plasticity problem is established in this context and extended to cover finite loading steps. Results are subsequently specialized to plane strain solids made of elastically isotropic materials with a plastic behavior governed by the von Mises criterion, accounting for combined isotropic and kinematic hardening. By exploiting previous results, the formulation is fully reduced to the plane. Further generalizations of the statements are also provided, which can be regarded as extensions to the elastic-plastic, plane strain problem of the Hu-Washizu principle in elasticity.     

Cyclic multiaxial and shear finite deformation responses of OFHC Cu. Part II: An extension to the KHL model and simulations

In this part, the Khan–Huang–Liang (KHL) constitutive model was extended to account for kinematic hardening characteristic behavior of materials. The extended model is then generalized and used to simulate experimental response of oxygen free high conductivity (OFHC) copper under cyclic shear straining and biaxial tension–torsion (multiaxial ratchetting) experiments presented in Part I (Khan et al., 2007). In addition, a new modification for the non-linear kinematic hardening rule of Karim–Ohno (Abdel-Karim and Ohno, 2000) is proposed to simulate multiaxial ratchetting behaviors. Although, the kinematic hardening contributes the most to the response, it is shown that, the loading rate effect, and a coupled isotropic and kinematic hardening effect should also be considered while simulating the multiaxial ratchetting behavior of OFHC copper. Furthermore, the newly modified kinematic hardening rules is able to fairly well simulate the multiaxial ratchetting experiments under different loading conditions, irrespective of the value of applied axial tensile stress, shear strain amplitude, pre-cyclic hardening and/or loading sequence.