A weakened form of Ilyushin's postulate and the structure of self-consistent Eulerian finite elastoplasticity

From a general standpoint in terms of internal variables, we formulate a general theory of self-consistent Eulerian finite elastoplasticity based on the additive decomposition of the Eulerian strain rate, i.e., D=D^e+D^p, as well as two consistency criteria. In this theory, the elastic behaviour is characterized by an exactly integrable elastic rate equation for D^e with a general form of complementary elastic potential. It is assumed that the yield function depends in a general manner on the Kirchhoff stress and the internal variables. Moreover, the plastic rate equation for D^p and the evolution equation for each internal variable are allowed to assume general forms relying on the just-mentioned variables and the stress rate. It is indicated that two consistency criteria, i.e., the self-consistency for the elastic rate equation and Prager's yielding stationarity, lead to the unique choice of objective rates, i.e., the logarithmic rate.The structure of the above theory is further studied and examined by virtue of a weakened form of Ilyushin's postulate. In a spinning frame defining the logarithmic rate, we introduce the notion of standard elastoplastic strain cycle, which starts at a point not on but inside a yield surface and incorporates only one infinitesimal plastic subpath. We show that this type of strain cycle is always possible. Then, by ruling out strain cycles starting at points on yield surfaces we propose a weakened form of Ilyushin's postulate, which says that the changing rate of the stress work done along every standard strain cycle should be non-negative, whenever the incorporated plastic subpath tends to vanish. By virtue of simple, rigorous procedures, we demonstrate that this weakened form of Ilyushin's postulate is adequate to ensure direct results concerning the normality rule and the convexity of the yield surface in the context of the foregoing Eulerian finite elastoplasticity theory. Specifically, with an exactly integrable elastic rate equation defining D^e, we prove that, in the space of the Kirchhoff stresses, the difference (D−D^e) is just the gradient of the yield function multiplied by a plastic multiplier, and thus bears the very kinematical and physical feature of plastic strain rate. Furthermore, we prove that, in the space of the Kirchhoff stresses, the elastic domain bounded by each yield surface should be convex. The main results are derived in a self-contained manner within the context of an Eulerian theory of finite elastoplasticity, without involving issues concerning how to define intermediate stress-free states and plastic strains, etc.     

Détermination du comportement macroscopique d'un milieu à fissures frottantesDetermination of the macroscopic behavior of a medium with frictional cracks

A micromechanical modeling of closed cracks as flat ellipsoidal inhomogeneities has been recently proposed in [Deudé et al., C. R. Mecanique 330 (2002) 589–599]. The present Note extends this approach to the case of a frictional contact between the crack lips. For von Mises friction, a linear hardening is obtained at the macroscopic scale, the state equation in terms of stress and strain rates being identical to that derived for unfrictional cracks. For Drucker–Prager friction, the micromechanical approach predicts a macroscopic dilatant behavior. To cite this article: J.-F. Barthélémy et al., C. R. Mecanique 331 (2003).     

L'identification des propriétés mécaniques de matériaux avec la méthode des champs virtuels,une alternative au recalage par éléments finisIdentification of the mechanical properties of materials with the virtual fields method,an alternative to finite element model updating

A method allowing the identification of parameters that govern constitutive equations of materials is proposed in this Note. These parameters are identified from heterogeneous strain fields. The method is direct so that iterative calculations based on finite element model updating are avoided. To cite this article: M. Grédiac et al., C. R. Mecanique 330 (2002) 107–112.     

Mechanics of the rate-dependent elastic–plastic deformation of glassy polymers from low to high strain rates

A combined experimental and analytical investigation has been performed to understand the mechanical behavior of two amorphous polymers—polycarbonate and poly(methyl methacrylate)—at strain rates ranging from 10⁻⁴ to 10⁴ s⁻¹. This range in strain rates was achieved in uniaxial tension and compression tests using a dynamic mechanical analyzer (DMA), a servo-hydraulic testing machine, and an aluminum split-Hopkinson pressure bar. DMA tension tests were used to characterize the viscoelastic behavior of these materials, with focus on the rate-dependent shift of material transition temperatures. Uniaxial compression tests on the servo-hydraulic machine (10⁻⁴ to 1 s⁻¹) and the split-Hopkinson pressure bar (10³ to 10⁴ s⁻¹) were used to characterize the rate-dependent yield and post-yield behavior. Both materials were observed to exhibit increased rate sensitivity of yield under the same strain rate/temperature conditions as the β-transition of the viscoelastic behavior. A physically based constitutive model for large strain deformation of thermoplastics was then extended to encompass high-rate conditions. The model accounts for the contributions of different molecular motions which become operational and important in different frequency regimes. The new features enable the model to not only capture the transition in the yield behavior, but also accurately predict the post-yield, large strain behavior over a wide range of temperatures and strain rates.     

Plastic deformation modeling of AL-6XN stainless steel at low and high strain rates and temperatures using a combination of bcc and fcc mechanisms of metals

Combination of physically based constitutive models for body centered cubic (bcc) and face centered cubic (fcc) metals developed recently by the authors [Voyiadjis, G.Z., Abed, F.H., 2005. Microstructural based models for bcc and fcc metals with temperature and strain rate dependency. Mech. Mater. 37, 355–378] are used in modeling the plastic deformation of AL-6XN stainless steel over a wide range of strain rates between 0.001 and 8300 s⁻¹ at temperatures from 77 to 1000 K. The concept of thermal activation analysis as well as the dislocation interaction mechanism is used in developing the plastic flow model for both the isothermal and adiabatic plastic deformation. In addition, the experimental observations of AL-6XN conducted by Nemat-Nasser et al. [Nemat-Nasser, S., Guo, W., Kihl, D., 2001. Thermomechanical response of AL-6XN stainless steel over a wide range of strain rates and temperatures, J. Mech. Phys. Solids 49, 1823–1846] are utilized in understanding the underlying deformation mechanisms. The plastic flow is considered in the range of temperatures and strain rates where diffusion and creep are not dominant, i.e., the plastic deformation is attributed to the motion of dislocations only. The modeling of the true stress–true strain curves for AL-6XN stainless steel is achieved using the classical secant modulus for the case of unidirectional deformation. The model parameters are obtained using the experimental results of three strain rates (0.001, 0.1, and 3500 s⁻¹). Good agreement is obtained between the experimental results and the model predictions. Moreover, the independency of the present model to the experiments used in the modeling is verified by comparing the theoretical results to an independent set of experimental data at the strain rate of 8300 s⁻¹ and various initial temperatures. Good correlation is observed between the model predictions and the experimental observations.     

Thermodynamic laws and consistent Eulerian formulation of finite elastoplasticity with thermal effects

Recently it has been demonstrated that, on the basis of the separation D=D^e+D^p arising from the split of the stress power and two consistency criteria for objective Eulerian rate formulations, it is possible to establish a consistent Eulerian rate formulation of finite elastoplasticity in terms of the Kirchhoff stress and the stretching, without involving additional deformation-like variables labelled “elastic” or “plastic”. It has further been demonstrated that this consistent formulation leads to a simple essential structure implied by the work postulate, namely, both the normality rule for plastic flow D^p and the convexity of the yield surface in Kirchhoff stress space. Here, we attempt to place such an Eulerian formulation on the thermodynamic grounds by extending it to a general case with thermal effects, where the consistency requirements are treated in a twofold sense. First, we propose a general constitutive formulation based on the foregoing separation as well as the two consistency criteria. This is accomplished by employing the corotational logarithmic rate and by incorporating an exactly integrable Eulerian rate equation for D^e for thermo-elastic behaviour. Then, we study the consistency of the formulation with thermodynamic laws. Towards this goal, simple forms of restrictions are derived, and consequences are discussed. It is shown that the proposed Eulerian formulation is free in the sense of thermodynamic consistency. Namely, a Helmholtz free energy function in explicit form may be found such that the restrictions from the thermodynamic laws can be fulfilled with positive internal dissipation for arbitrary forms of constitutive functions included in the constitutive formulation. In particular, that is the case for the foregoing essential constitutive structure in the purely mechanical case. These results eventually lead to a complete, explicit constitutive theory for coupled fields of deformation, stress and temperature in thermo-elastoplastic solids at finite deformations.     

Phenomenological nonlinear modelling of glassy polymers

Compared with the numerous works into the constitutive equations for the mechanical behaviour of metals, very little attention has been devoted to those of polymers. However, a model is required to describe both the complex shape of the stress–strain curves and strain rate sensitivity of glassy polymers. In this Note, a unified viscoelastic-viscoplastic model is presented in which the nonlinear pre-yield behaviour, the strain softening and the strain hardening are described by internal variables, in analogy with the models developed for metals. In order to check the predictive capability of the model, the numerical results are compared with the experimental data (monotone, creep and relaxation tests) of a typical amorphous glassy polymer. To cite this article: F. Zaïri et al., C. R. Mecanique 333 (2005).     

Artificial neural network prediction to the hot compressive deformation behavior of Al–Cu–Mg–Ag heat-resistant aluminum alloy

The behavior of the flow stress of Al-Cu-Mg-Ag heat-resistant aluminum alloys during hot compression deformation was studied by thermal simulation test. The temperature and the strain rate during hot compression were 340-500 °C, 0.001 s⁻¹ to 10 s⁻¹, respectively. Constitutive equations and an artificial neural network (ANN) model were developed for the analysis and simulation of the flow behavior of the Al-Cu-Mg-Ag alloys. The inputs of the model are temperature, strain rate and strain. The output of the model is the flow stress. Comparison between constitutive equations and ANN results shows that ANN model has a better prediction power than the constitutive equations.     

Formulation DNLR intégrale des lois de comportement : généralisation du principe de superposition de BoltzmannDNLR integral formulation of behaviour laws: generalisation of Boltzmann's superposition principle.

A thermodynamic approach of relaxation phenomena called DNLR (Distribution of Non Linear Relaxations) has been employed for about fifteen years to describe the behaviour of materials. In this Note, we show that different DNLR constitutive equations, initially written in a rate formulation, can be put in a time-integrated (hereditary) form which generalises the superposition principle to non linear fields like viscoplasticity or damage behaviours. To cite this article: B. Martin et al., C. R. Mecanique 332 (2004).     

Mechanical and optical characterization of an anelastic polymer at large strain rates and large strains

Two constitutive relations have been determined from test results that characterize, respectively, the uniaxial and photomechanical behavior of a polyester-styrene copolymer for strain rates from 10^?5 to 3×10³ in./in./s and strains up to 40 percent. The high-strain-rate data were obtained by means of a split-Hopkinson-bar apparatus. Intermediate-strain-rate tests, performed with the aid of a drop tower, were reported in an earlier paper. Quasi-static experiments were conducted on a standard testing machine. A nonlinear, four-parameter, elastic-viscoplastic model was constructed which describes the mechanical behavior. The parameters were determined by a least-mean-squares curve-fitting procedure. The viscoplastic parameters were found to obey a power law in strain rate. The photomechanical model was found to be linear with strain well into the plastic-deformation region, while the slope of the strain-birefringence curve for each strain rate also varied by strain rate to a power.     

A Phenomenological Constitutive Model for Describing Thermo-Viscoplastic Behavior of Al-Zn-Mg-Cu Alloy Under Hot Working Condition

In order to predict the high-temperature deformation behavior of Al-Zn-Mg-Cu alloy, the hot compression tests were conducted in the strain rate range of (0.001–0.1)s⁻¹ and the forming temperature range of (573–723) K. Based on the experimental results, Johnson-Cook model was found inadequate to describe the high-temperature deformation behavior of Al-Zn-Mg-Cu alloy. Therefore, a new phenomenological constitutive model is proposed, considering the coupled effects of strain, strain rate and forming temperature on the material flow behavior of Al-Zn-Mg-Cu alloy. In the proposed model, the material constants are presented as functions of strain rate. The proposed constitutive model correlates well with the experimental results confirming that the proposed model can give an accurate and precise estimate of flow stress for the Al-Zn-Mg-Cu alloy investigated in this study.     

Duality in constitutive formulation of finite-strain elastoplasticity based on F=FeFp and F=FF decompositions

A constitutive theory for large elastic–plastic deformations is presented by employing F=F^pF^e decomposition of the total deformation gradient. A duality in constitutive formulation based on this and the well-known Lee's decomposition F=F_eF_p is established for isotropic polycrystalline and single crystal plasticity.     

Duality in constitutive formulation of finite-strain elastoplasticity based on F=FeFp and F=FpFe decompositions

A constitutive theory for large elastic–plastic deformations is presented by employing F=F^pF^e decomposition of the total deformation gradient. A duality in constitutive formulation based on this and the well-known Lee's decomposition F=F_eF_p is established for isotropic polycrystalline and single crystal plasticity.     

论塑性旋率本构方程的必要性

在有限变形弹塑性理论中,本构方程通常是以率型形式给出的。因此,应变率的分解将是一个十分基本的问题。在当前,较为流行的是基于中间构形的应变率分解,而这其中最有影响的有 Lee,E.H.等人的工作和 Dafalias 等人的工作。然而,本文的研究表明,至少在某些特殊情况下,我们可以得到与微观子结构定向旋率的有关表达式。这就使给出塑性旋率本构方程变得不必要了。显然,本文的结果既不同于 Dafalias 的工作,也不同于 E.H.Lee等人的工作。前者需要通过塑性旋率的本构方程来确定微观子结构定向的旋率,而后者则需要作出附加的隐含假设来避免给出塑性旋率的本构方程。可以相信,本文工作将可能为有限弹塑性变形的本构理论提供一种新的途径。     

Mechanically alloyed nanocrystalline iron and copper mixture: behavior and constitutive modeling over a wide range of strain rates

A comprehensive study on the response of a nanocrystalline iron and copper mixture (80% Fe and 20% Cu) to quasi-static and dynamic loading is performed. The constitutive model developed earlier by Khan, Huang & Liang (KHL) is extended to include the responses of nanocrystalline metallic materials. The strain rate and grain size dependent behaviors of porous nanocrystalline iron-copper mixture were determined experimentally for both static and dynamic loading. A viscoplastic model is formulated by associating the modified KHL model (representing the fully dense matrix behavior), and Gurson's plastic potential which provides the yield criteria for porous material. Simulations of uniaxial compressive deformations of iron-copper mixture with different initial porosity, grain size and at a wide range of strain rate (10⁻⁴ to 10³ s⁻¹) are made. The numerical results correlate well with the experimental observations.     

An implicit tensorial gradient plasticity model – Formulation and comparison with a scalar gradient model

Many rate-independent models for metals utilize the gradient of effective plastic strain to capture size-dependent behavior. This enhancement, sometimes termed as “explicit” gradient formulation, requires higher-order tractions to be imposed on the evolving elasto-plastic boundary and the resulting numerical framework is complicated. An “implicit” scalar gradient model was thus developed in Peerlings [Peerlings, R.H.J., 2007. On the role of moving elastic–plastic boundaries in strain gradient plasticity. Model. Simul. Mater. Sci. Eng. 15, 109–120] that has only C⁰ continuity requirements and its implementation is straightforward. However, both explicit and implicit scalar gradient models can be problematic when the effective plastic strains do not have smooth profiles. To address this limitation, an implicit tensorial gradient model is proposed in this paper based on the generalized micromorphic framework. It is also demonstrated that the scalar and tensorial implicit gradient models give similar results when the effective plastic strains fluctuate smoothly.     

A viscoplastic constitutive model for nickel-base superalloy,part 2: modeling under anisothermal conditions

In Part 2 of this study, extensive deformation tests were carried out on the nickel-base polycrystalline superalloy IN738LC under isothermal and anisothermal conditions between 450 and 950 °C. Under the isothermal conditions, the material showed almost no rate/time-dependency below 700 °C, while it showed distinct rate/time-dependency above 800 °C. Regarding the cyclic deformation, slight cyclic hardening behavior was observed when the temperature was below 700 °C and the imposed strain rate was fast, whereas in the case of the temperature above 800 °C or under slower strain rate conditions, the cyclic hardening behavior was scarcely observed. Unique inelastic behavior was observed under in-phase and out-of-phase anisothermal conditions: with an increase in the number of cycles, the stress at higher temperatures became smaller and the stress at lower temperatures became larger in absolute value although the stress range was approximately constant during the cyclic loading. In other words, the mean stress continues to evolve cycle-by-cycle in the direction of the stress at lower temperatures. Based on the experimental results, it was assumed that evolution of the variable Y that had been incorporated into a kinematic hardening rule in Part 1 of this study is active under higher temperatures and is negligible under lower temperatures. The material constants used in the constitutive equations were determined with the isothermal data, and were expressed as functions of temperature empirically. The extended viscoplastic constitutive equations were applied to the anisothermal cyclic loading as well as the monotonic tension, stress relaxation, creep and cyclic loading under the isothermal conditions. It was demonstrated that the present viscoplastic constitutive model was successful in describing the inelastic behavior of the material adequately, including the anomalous inelastic behavior observed under the anisothermal conditions, owing to the consideration of the variable Y.     

A material model for finite elasto-plastic deformations considering a substructure

Developing further the substructure models proposed by Mandel and Dafalias a thermodynamically consistent system of differential and algebraic equations is derived to describe anisotropic elasto-plastic material behavior at finite deformations. Based on the multiplicative split of the deformation gradient an appropriate material law is formulated applying the principle of the maximum of plastic dissipation. Generalized basic relations of this material model containing a relation of hyperelasticity, evolutional equations for the internal variables describing different kinds of hardening, and the yield condition are presented. The capacity of the proposed material model is demonstrated on the example of a sheet with a hole. Presenting the evolution of yield surfaces the capability of the model to describe anisotropic hardening behavior is shown.     

Experimental investigation and modeling of the tension behavior of polycarbonate with temperature effects from low to high strain rates

The effects of strain rate and temperature on the tension stress–strain responses of polycarbonate are experimentally investigated over a wide range of strain rates (0.001–1700 s⁻¹) and temperatures (0–120 °C). A modified split Hopkinson tension bar is used for high-rate uniaxial tension tests. Experimental results indicate that the stress–strain responses of polycarbonate at high strain rates exhibit the nonlinear characteristics including the obvious yielding and strain softening. The tension behavior is strongly dependent on the strain rate and temperature. The values of yield stress and strain at yield present a dramatic increase at higher strain rates and decrease with the increase in temperature. Moreover, there exists a significant rate-sensitivity transition in the polycarbonate tension yield behavior. Based on the experimental investigation, a physically based three-dimensional elastoplastic constitutive model for the finite deformation of glassy polymers is used to characterize the rate-temperature dependent yield and post-yield behavior of polycarbonate when subjected to tension loading. The model results are shown close to the experimental data within the investigated strain-rate and temperature ranges.