Strain gradient plasticity,strengthening effects and plastic limit analysis

Within the framework of isotropic strain gradient plasticity, a rate-independent constitutive model exhibiting size dependent hardening is formulated and discussed with particular concern to its strengthening behavior. The latter is modelled as a (fictitious) isotropic hardening featured by a potential which is a positively degree-one homogeneous function of the effective plastic strain and its gradient. This potential leads to a strengthening law in which the strengthening stress, i.e. the increase of the plastically undeformed material initial yield stress, is related to the effective plastic strain through a second order PDE and related higher order boundary conditions. The plasticity flow laws, with the role there played by the strengthening stress, are addressed and shown to admit a maximum dissipation principle. For an idealized elastic perfectly plastic material with strengthening effects, the plastic collapse load problem of a micro/nano scale structure is addressed and its basic features under the light of classical plastic limit analysis are pointed out. It is found that the conceptual framework of classical limit analysis, including the notion of rigid-plastic behavior, remains valid. The lower bound and upper bound theorems of classical limit analysis are extended to strengthening materials. A static-type maximum principle and a kinematic-type minimum principle, consequences of the lower and upper bound theorems, respectively, are each independently shown to solve the collapse load problem. These principles coincide with their respective classical counterparts in the case of simple material. Comparisons with existing theories are provided. An application of this nonclassical plastic limit analysis to a simple shear model is also presented, in which the plastic collapse load is shown to increase with the decreasing sample size (Hall–Petch size effects).     

Size effects in generalised continuum crystal plasticity for two-phase laminates

The solutions of a boundary value problem are explored for various classes of generalised crystal plasticity models including Cosserat, strain gradient and micromorphic crystal plasticity. The considered microstructure consists of a two-phase laminate containing a purely elastic and an elasto-plastic phase undergoing single or double slip. The local distributions of plastic slip, lattice rotation and stresses are derived when the microstructure is subjected to simple shear. The arising size effects are characterised by the overall extra back stress component resulting from the action of higher order stresses, a characteristic length l_c describing the size-dependent domain of material response, and by the corresponding scaling law lⁿ as a function of microstructural length scale, l. Explicit relations for these quantities are derived and compared for the different models. The conditions at the interface between the elastic and elasto-plastic phases are shown to play a major role in the solution. A range of material parameters is shown to exist for which the Cosserat and micromorphic approaches exhibit the same behaviour. The models display in general significantly different asymptotic regimes for small microstructural length scales. Scaling power laws with the exponent continuously ranging from 0 to −2 are obtained depending on the values of the material parameters. The unusual exponent value −2 is obtained for the strain gradient plasticity model, denoted “curl H^p” in this work. These results provide guidelines for the identification of higher order material parameters of crystal plasticity models from experimental data, such as precipitate size effects in precipitate strengthened alloys.     

Kinematic hardening in large strain plasticity

A finite strain hyper elasto-plastic constitutive model capable to describe non-linear kinematic hardening as well as non-linear isotropic hardening is presented. In addition to the intermediate configuration and in order to model kinematic hardening, an additional configuration is introduced – the center configuration; both configurations are chosen to be isoclinic. The yield condition is formulated in terms of the Mandel stress and a back-stress with a structure similar to the Mandel stress.It is shown that the non-dissipative part of the plastic velocity gradient not governed by the thermodynamical framework and the corresponding quantity associated with the kinematic hardening influence the material behaviour to a large extent when kinematic hardening is present. However, for isotropic elasticity and isotropic hardening plasticity it is shown that the non-dissipative quantities have no influence upon the stress–strain relation.As an example, kinematic hardening von Mises plasticity is considered, which fulfils the plastic incompressibility condition and is independent of the hydrostatic pressure. To evaluate the response and to examine the influence of the non-dissipative quantities, simple shear is considered; no stress oscillations occur.     

A unified residual-based thermodynamic framework for strain gradient theories of plasticity

A unified thermodynamic framework for gradient plasticity theories in small deformations is provided, which is able to accommodate (almost) all existing strain gradient plasticity theories. The concept of energy residual (the long range power density transferred to the generic particle from the surrounding material and locally spent to sustain some extra plastic power) plays a crucial role. An energy balance principle for the extra plastic power leads to a representation formula of the energy residual in terms of a long range stress, typically of the third order, a macroscopic counterpart of the micro-forces acting on the GNDs (Geometrically Necessary Dislocations). The insulation condition (implying that no long range energy interactions are allowed between the body and the exterior environment) is used to derive the higher order boundary conditions, as well as to ascertain a principle of the plastic power redistribution in which the energy residual plays the role of redistributor and guarantees that the actual plastic dissipation satisfies the second thermodynamics principle. The (nonlocal) Clausius-Duhem inequality, into which the long range stress enters aside the Cauchy stress, is used to derive the thermodynamic restrictions on the constitutive equations, which include the state equations and the dissipation inequality. Consistent with the latter inequality, the evolution laws are formulated for rate-independent models. These are shown to exhibit multiple size effects, namely (energetic) size effects on the hardening rate, as well as combined (dissipative) size effects on both the yield strength (intrinsic resistance to the onset of plastic strain) and the flow strength (resistance exhibited during plastic flow). A friction analogy is proposed as an aid for a better understanding of these two kinds of strengthening effects. The relevant boundary-value rate problem is addressed, for which a solution uniqueness theorem and a minimum variational principle are provided. Comparisons with other existing gradient theories are presented. The dissipation redistribution mechanism is illustrated by means of a simple shear model.     

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.     

Size effects on the plastic collapse limit load of thin foils in bending and thin wires in torsion

Following a previous paper by the author [Strain gradient plasticity, strengthening effects and plastic limit analysis, Int. J. Solids Struct. 47 (2010) 100–112], a nonconventional plastic limit analysis for a particular class of micron scale structures as, typically, thin foils in bending and thin wires in torsion, is here addressed. An idealized rigid-perfectly plastic material is considered, which is featured by a strengthening potential degree-one homogeneous function of the effective plastic strain and its spatial gradient. The nonlocal (gradient) nature of the material resides in the inherent strengthening law, whereby the yield strength is related to the effective plastic strain through a second order PDE with associated higher order boundary conditions. The peculiarity of the considered structures stems from their geometry and loading conditions, which dictate the shape of the collapse mechanism and make the higher order boundary conditions on the (microscopically) free boundary be accommodated by means of a boundary singularity mechanism. This consists in the formation of thin boundary layers with unbounded stresses, but bounded stress resultants which —together with the regular bulk stresses— contribute to the value of the collapse load. Closed-form solutions are provided for thin foils in pure bending and thin wires in pure torsion, and in particular the limit bending and torque moments are given as functions of an adimensionalized internal length parameter.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

A constitutive model for orthotropic elasto-plasticity at large strains

Summary The paper presents a thermodynamically consistent constitutive model for elasto-plastic analysis of orthotropic materials at large strain. The elastic and plastic anisotropies are assumed to be persistent in the material but the anisotropy axes can undergo a rigid rotation due to large plastic deformations. The orthotropic yield function is formulated in terms of the generally nonsymmetric Mandel stress tensor such that its skew-symmetric part is additionally taken into account. Special attention is focused on the convexity of the yield surface resulting in the nine-dimensional stress space. Of particular interest are new convexity conditions which do not appear in the classical theory of anisotropic plasticity. They impose additional constraints on the material constants governing the plastic spin. The role of the plastic spin is further studied in simple shear accompanied by large elastic and large plastic deformations. If the plastic spin is neglected, the shear stress response is characterized by oscillations with an amplitude strictly dependent on the degree of the plastic anisotropy.accepted for publication 2 March 2004     

Micromorphic approach for gradient-extended thermo-elastic–plastic solids in the logarithmic strain space

The coupled thermo-mechanical strain gradient plasticity theory that accounts for microstructure-based size effects is outlined within this work. It extends the recent work of Miehe et al. (Comput Methods Appl Mech Eng 268:704–734, 2014) to account for thermal effects at finite strains. From the computational viewpoint, the finite element design of the coupled problem is not straightforward and requires additional strategies due to the difficulties near the elastic–plastic boundaries. To simplify the finite element formulation, we extend it toward the micromorphic approach to gradient thermo-plasticity model in the logarithmic strain space. The key point is the introduction of dual local–global field variables via a penalty method, where only the global fields are restricted by boundary conditions. Hence, the problem of restricting the gradient variable to the plastic domain is relaxed, which makes the formulation very attractive for finite element implementation as discussed in Forest (J Eng Mech 135:117–131, 2009) and Miehe et al. (Philos Trans R Soc A Math Phys Eng Sci 374:20150170, 2016).     

Size effects and idealized dislocation microstructure at small scales: Predictions of a Phenomenological model of Mesoscopic Field Dislocation Mechanics: Part I

A Phenomenological Mesoscopic Field Dislocation Mechanics (PMFDM) model is developed, extending continuum plasticity theory for studying initial-boundary value problems of small-scale plasticity. PMFDM results from an elementary space-time averaging of the equations of Field Dislocation Mechanics (FDM), followed by a closure assumption from any strain-gradient plasticity model that attempts to account for effects of geometrically necessary dislocations (GNDs) only in work hardening. The specific lower-order gradient plasticity model chosen to substantiate this work requires one additional material parameter compared to its conventional continuum plasticity counterpart. The further addition of dislocation mechanics requires no additional material parameters. The model (a) retains the constitutive dependence of the free-energy only on elastic strain as in conventional continuum plasticity with no explicit dependence on dislocation density, (b) does not require higher-order stresses, and (c) does not require a constitutive specification of a ‘back-stress’ in the expression for average dislocation velocity/plastic strain rate. However, long-range stress effects of average dislocation distributions are predicted by the model in a mechanistically rigorous sense. Plausible boundary conditions (with obvious implication for corresponding interface conditions) are discussed in some detail from a physical point of view. Energetic and dissipative aspects of the model are also discussed. The developed framework is a continuous-time model of averaged dislocation plasticity, without having to rely on the notion of incremental work functions, their convexity properties, or their minimization. The tangent modulus relating stress rate and total strain rate in the model is the positive-definite tensor of linear elasticity, and this is not an impediment to the development of idealized microstructure in the theory and computations, even when such a convexity property is preserved in a computational scheme. A model of finite deformation, mesoscopic single crystal plasticity is also presented, motivated by the above considerations.Lower-order gradient plasticity appears as a constitutive limit of PMFDM, and the development suggests a plausible boundary condition on the plastic strain rate for this limit that is appropriate for the modeling of constrained plastic flow in three-dimensional situations.     

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

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

A finite element investigation of the effect of crack tip constraint on hole growth under mode i and mixed mode loading

In this work, the effect of constraint on hole growth near a notch tip in a ductile material under mode I and mixed mode loading (involving modes I and II) is investigated. To this end, a 2-D plane strain, modified boundary layer formulation is employed in which the mixed mode elastic K–T field is prescribed as remote boundary conditions. A finite element procedure that accounts for finite deformations and rotations is used along with an appropriate version of J₂ flow theory of plasticity with small elastic strains. Several analyses are carried out corresponding to different values of T-stress and remote elastic mode-mixity. The interaction between the notch and hole is studied by examining the distribution of hydrostatic stress and equivalent plastic strain in the ligament between the notch tip and the hole, as well as the growth of the hole. The implications of the above results on ductile fracture initiation due to micro-void coalescence are discussed.     

A conventional theory of mechanism-based strain gradient plasticity

There exist two frameworks of strain gradient plasticity theories to model size effects observed at the micron and sub-micron scales in experiments. The first framework involves the higher-order stress and therefore requires extra boundary conditions, such as the theory of mechanism-based strain gradient (MSG) plasticity [J Mech Phys Solids 47 (1999) 1239; J Mech Phys Solids 48 (2000) 99; J Mater Res 15 (2000) 1786] established from the Taylor dislocation model. The other framework does not involve the higher-order stress, and the strain gradient effect come into play via the incremental plastic moduli. A conventional theory of mechanism-based strain gradient plasticity is established in this paper. It is also based on the Taylor dislocation model, but it does not involve the higher-order stress and therefore falls into the second strain gradient plasticity framework that preserves the structure of conventional plasticity theories. The plastic strain gradient appears only in the constitutive model, and the equilibrium equations and boundary conditions are the same as the conventional continuum theories. It is shown that the difference between this theory and the higher-order MSG plasticity theory based on the same dislocation model is only significant within a thin boundary layer of the solid.     

Interfacial energy effects within the framework of strain gradient plasticity

In the framework of strain gradient plasticity, a solid body with boundary surface playing the role of a dissipative boundary layer endowed with surface tension and surface energy, is addressed. Using the so-called residual-based gradient plasticity theory, the state equations and the higher order boundary conditions are derived quite naturally for both the bulk material and the boundary layer. A phenomenological constitutive model is envisioned, in which the bulk material and the boundary layer obey (rate independent associative) coupled plasticity evolution laws, with kinematic hardening laws of differential nature for the bulk material, but of nondifferential nature for the layer. A combined global maximum dissipation principle is shown to hold. The higher order boundary conditions are discussed in details and categorized in relation to some peculiar features of the boundary surface, and their basic role in the coupling of the bulk/layer plasticity evolution laws is pointed out. The case of an internal interface is also studied. An illustrative example relating to a shear model exhibiting energetic size effects is presented. The theory provides a unified view on gradient plasticity with interfacial energy effects.     

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.     

Extremum principles for certain hyperelastic materials

A hyperelastic material is here said to be of class H_m if the elastic potential is a homogeneous function of order m + 1 in the components of the Lagrangian displacement gradient. It is shown that a single solution to a boundary value problem generates an infinite family of solutions to a family of related boundary value problems. Assuming that a solution to a boundary value problem exists, it is shown that it is unique provided that the material is stable in the sense of Hill in a deleted neighbourhood of the stress-free state. A minimum theorem concerning the strain energy and the virtual work of the prescribed forces is established for the equilibrium configurations, and a maximum theorem concerning the virtual work of the prescribed surface displacements and the complementary stress energy is established for compatible stress fields. As an application, upper and lower bounds are found for the torsional stiffness of a cylindrical bar of square cross section under infinitesimal deformation.     

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

A finite strain kinematic hardening constitutive model based on Hencky strain: General framework, solution algorithm and application to shape memory alloys

The logarithmic or Hencky strain measure is a favored measure of strain due to its remarkable properties in large deformation problems. Compared with other strain measures, e.g., the commonly used Green-Lagrange measure, logarithmic strain is a more physical measure of strain. In this paper, we present a Hencky-based phenomenological finite strain kinematic hardening, non-associated constitutive model, developed within the framework of irreversible thermodynamics with internal variables. The derivation is based on the multiplicative decomposition of the deformation gradient into elastic and inelastic parts, and on the use of the isotropic property of the Helmholtz strain energy function. We also use the fact that the corotational rate of the Eulerian Hencky strain associated with the so-called logarithmic spin is equal to the strain rate tensor (symmetric part of the velocity gradient tensor). Satisfying the second law of thermodynamics in the Clausius-Duhem inequality form, we derive a thermodynamically-consistent constitutive model in a Lagrangian form. In comparison with the available finite strain models in which the unsymmetric Mandel stress appears in the equations, the proposed constitutive model includes only symmetric variables. Introducing a logarithmic mapping, we also present an appropriate form of the proposed constitutive equations in the time-discrete frame. We then apply the developed constitutive model to shape memory alloys and propose a well-defined, non-singular definition for model variables. In addition, we present a nucleation-completion condition in constructing the solution algorithm. We finally solve several boundary value problems to demonstrate the proposed model features as well as the numerical counterpart capabilities.