Modelling large deformation behaviour under loading–unloading of semicrystalline polymers: Application to a high density polyethylene

In this work, the large deformation behaviour under monotonic loading and unloading of a high density polyethylene (HDPE) is studied. To analyze the nonlinear time-dependent response of the material, mechanical tests were conducted at room temperature under constant true strain rates and stress relaxation conditions. A physically-based inelastic model written under finite strain formulation is proposed to describe the mechanical behaviour of HDPE. In the model, the inelastic mechanisms involve two parallel elements: a visco-hyperelastic network resistance acting in parallel with a viscoelastic–viscoplastic intermolecular resistance where the amorphous and crystalline phases are explicitly taken into consideration. The semicrystalline polymer is considered as a two-phase composite. The influence of the crystallinity on the loading and unloading behaviour is investigated. Numerical results are compared to experimental data. It is shown that the model is able to accurately reproduce the experimental observations corresponding to monotonic loading, unloading and stress relaxation behaviours at different strain levels. Finally, the model capabilities to capture cyclic loading–unloading behaviour up to large strains are discussed. To demonstrate the improved modelling capabilities, simulations are also performed using the original model of Boyce et al. [Boyce, M.C., Socrate, S., Llana, P.G., 2000. Constitutive model for the finite deformation stress–strain behavior of poly(ethylene terephthalate) above the glass transition. Polymer 41, 2183–2201] modified by Ahzi et al. [Ahzi, S., Makradi, A., Gregory, R.V., Edie, D.D., 2003. Modeling of deformation behavior and strain-induced crystallization in poly(ethylene terephthalate) above the glass transition temperature. Mechanics of Materials 35, 1139–1148].     

A physically-based constitutive model for anisotropic damage in rubber-toughened glassy polymers during finite deformation

The present work focuses on the development of a physically-based model for large deformation stress-strain response and anisotropic damage in rubber-toughened glassy polymers. The main features leading to a microstructural evolution (regarding cavitation, void aspect ratio, matrix plastic anisotropy and rubbery phase deformation) in rubber-toughened glassy polymers are introduced in the proposed constitutive model. The constitutive response of the glassy polymer matrix is modelled using the hyperelastic-viscoplastic model of [Boyce et al., 1988] and [Boyce et al., 2000]. The deformation mechanisms of the matrix material are accounted for by two resistances: an elastic-viscoplastic isotropic intermolecular resistance acting in parallel with a visco-hyperelastic anisotropic network resistance, each resistance being modified to account for damage effects by void growth with a variation of the void aspect ratio. The effective contribution of the hyperelastic particles to the overall composite behaviour is taken into account by treating the overall system in a composite scheme framework. The capabilities of the proposed constitutive model are checked by comparing experimental data with numerical simulations. The deformation behaviour of rubber-toughened poly(methyl methacrylate) was investigated experimentally in tension at a temperature of 80 °C and for different constant true strain rates monitored by a video-controlled technique. The reinforcing phase is of the soft core-hard shell type and its diameter is of the order of one hundred nanometers. The particle volume fraction was adjusted from 15% to 45% by increments of 5%. The stress-strain response and the inelastic volumetric strain are found to depend markedly on particle volume fraction. For a wide range of rubber volume fractions, the model simulations are in good agreement with the experimental results. Finally, a parametric analysis demonstrates the importance of accounting for void shape, matrix plastic anisotropy and rubber content.     

A crystal plasticity model for strain-path changes in metals

A model is proposed that deals with the transient mechanical anisotropy during strain-path changes in metals. The basic mechanism is assumed to be latent hardening or softening of the slip systems, dependent on if they are active or passive during deformation, reflecting microstructural mechanisms that depend on the deformation mode rather than on the crystallography. The new model captures the experimentally observed behaviour of cross hardening in agreement with experiments for an AA3103 aluminium alloy. Generic results for strain reversals qualitatively agree with two types of behaviour reported in the literature – with or without a plateau on the stress–strain curve. The influence of the model parameters is studied through detailed calculations of the response of three selected parameter combinations, including the evolution of yield surface sections subsequent to 10% pre-strain. The mathematical complexity is kept to a minimum by avoiding explicit predictions related directly to underpinning microstructural changes. The starting point of the model is a combination of conventional texture and work hardening approaches, where an adapted full-constraints Taylor theory and a simple single-crystal work-hardening model for monotonic strain are used. However, the framework of the model is not restricted to these particular models.     

Modelling the behaviour of polycrystalline austenitic steel with twinning-induced plasticity effect

A micromechanical model using the scale transition method in elastoviscoplasticity has been developed to describe the behaviour of those austenitic steels that display a TWIP effect. A physically based constitutive equation at the grain scale is proposed considering two inelastic strain modes: crystallographic slip and twinning. The typical organizations of microtwins observed in electron microscopy are considered, and the twin–slip as well as the twin–twin interactions are accounted for. The parameters for slip are first fitted on the uniaxial tensile response obtained at intermediate temperatures (when twinning is inhibited). Then, the parameters associated with twinning are identified using the stress–strain curve at room temperature. The simulated results in both macro and micro scales are in good agreement with experimentally obtained results.     

Constitutive modeling of biodegradable polymers: Hydrolytic degradation and time-dependent behavior

A large range of biodegradable polymers has been used to produce implantable medical devices. Apart from biological compatibility, these devices shall be also functional compatible and perform adequate mechanical temporary support during the healing process. However, the mechanical behavior of biodegradable materials during its degradation, which is an important aspect of the design of these biodegradable devices, is still an unexplored subject. Based on the literature, the mechanical behavior of biodegradable polymers is strain rate dependent and exhibits hysteresis upon cyclic loading. On the other hand, ductility, toughness and strength of the material decay during hydrolytic degradation. In this work, it is considered a three-dimensional time-dependent model adapted from the one developed by Bergström and Boyce to simulate the performance of biodegradable structures undergoing large deformations incorporating the hydrolysis degradation. Since this model assumes that the mechanical behavior is divided into a time independent network and a non-linear time-dependent network, it enables to simulate the monotonic tests of a biodegradable structure loaded under different strain rates. The hysteresis effects during unloading–reloading cycles at different strain levels can be predicted by the model. A parametric study of the material model parameters evolution during the hydrolytic degradation was conducted to identify which parameters are more sensible to this degradation process. The investigated model could predict very well the experimental results of a blend of polylactic acid and polycaprolactone (PLA–PCL) in the full range of strains until rupture during hydrolytic degradation. From these results and analyses, a method is proposed to simulate the three-dimensional mechanical behavior during hydrolytic degradation.     

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

Determination of the effective elastic–plastic response of metal–ceramic composites

The mechanical response of metal–ceramic composites is analysed through a homogenization model accounting for the mechanical behaviour of the constituent materials. In order to achieve this purpose a nonlinear homogenization method based on the phase field approach has been suitably implemented into a numerical code. A prescribed homogenized strain state is applied to a unit volume element of a metal–ceramic composite with proportional loading in which all components of the strain tensor are proportional to one scalar parameter. The mechanical response of the material has been modeled by considering a von Mises plasticity model for the metal phase and a Drucker–Prager associative elastic–plastic material model for the ceramic phase. A two stages plasticity has been obtained in which inelastic strain develops in the metal phase followed by a fully plastic response. A comparison with a finite element model of the stress–strain response of an axisymmetric unit cell has been carried out with the purpose to validate the homogenization based modeling presented in the paper. Plastic parameters of a Drucker–Prager yield surface for the homogenized composite have been calculated at different materials compositions. Associative Drucker–Prager plasticity has been found to be accurate for high ceramic content.     

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.     

New flow rules in elasto-viscoplastic constitutive models for spheroidal graphite cast-iron

A specific flow rules and the corresponding constitutive elasto-viscoplastic model combined with new experimental strategy are introduced in order to represent a spheroidal graphite cast-iron behaviour on a wide range of strain, strain rate and temperature. A “full model” is first proposed to correctly reproduce the alloy behaviour even for very small strain levels. A “light model” with a bit poorer experimental agreement but a simpler formulation is also proposed. These macroscopic models, whose equations are based on physical phenomena observed at the dislocation scale, are able to cope with the various load conditions tested – progressive straining and cyclic hardening tests – and to correctly describe anisothermal evolution. The accuracy of these two models and the experimental databases to which they are linked is estimated on different types of experimental tests and compared with the accuracy of more standard Chaboche-type constitutive models. Each test leads to the superiority of the “full model”, particularly for slow strain rates regimes. After developing a material user subroutine, FEM simulations are performed on Abaqus for a car engine exhaust manifold and confirm the good results obtained from the experimental basis. We obtain more accurate results than those given by more traditional laws. A very good correlation is observed between the simulations and the engine bench tests.     

Constitutive model for uniaxial transformation ratchetting of super-elastic NiTi shape memory alloy at room temperature

The transformation ratchetting of super-elastic NiTi shape memory alloy was observed by the uniaxial stress-controlled cyclic tests [Kang, G.Z., Kan, Q.H., Qian, L.M., Liu, Y.J, 2009a. Ratchetting deformation of super-elastic and shape memory NiTi Alloys. Mech. Mater. 41, 139–153]. It is concluded that the NiTi alloy presents apparent ratchetting behaviour, and the ratchetting is collectively caused by the cyclic accumulation of residual induced-martensite and the transformation-induced plastic deformation (i.e., namely transformation ratchetting). Based on the experimental results, a cyclic constitutive model was constructed in the framework of generalized plasticity [Lubliner, J., Auricchio, F., 1996. Generalized plasticity and shape memory alloys. Int. J. Solids Struct. 33, 991–1003] to describe the transformation ratchetting of super-elastic NiTi alloy. The proposed model simultaneously accounts for the evolutions of residual induced-martensite and transformation-induced plastic strain during the stress-controlled cyclic loading by introducing an internal variable z_c, i.e., cumulated induced-martensite volume fraction. The dependence of transformation ratchetting on the applied stress levels and the phase transformation hardening behaviour of the NiTi alloy are also considered in the developed model. The anisotropic phase transformation behaviours of the alloy presented in the tension and compression cases are described by employing a Drucker–Prager-typed transformation surface. It is shown that the simulated results of transformation ratchetting obtained by the proposed model are in good agreement with the corresponding experiments, since the typical features of transformation ratchetting are reasonably captured by the proposed model.     

Saturated and unsaturated behaviour modelling of Meuse–Haute/Marne argillite

Poromechanical behaviour modelling of the Callovo–Oxfordian argillite under saturated and partially saturated conditions is proposed using the equivalent stress concept. In comparison with the previous works on this rock, the particular form of the yield criterion and the plastic flow potential proposed here help to better describe the rock behaviour in tension–stress paths. The evolution of the poroelastic parameters due to the induced cracks is also considered in a simple way. Due to its physical nature, different from classical soils, the partially saturated behaviour of this rock could not be correctly described by any of the pre-existing theories usually used for the partially saturated porous media (soils). Based upon experimental results on this rock and developing some ideas proposed by other authors, an extension of the saturated elastic–plastic model in unsaturated field is proposed. The key hypothesis of this extension is the evolution of the Biot’s coefficient as a function of the suction, justified by laboratory results and micromechanical analyses. The predictions of the model in saturated and partially saturated conditions are compared with laboratory results and a good general agreement is found.     

Micromechanical modelling of the effect of plastic deformation on the mechanical behaviour in pseudoelastic shape memory alloys

Except for the recoverable strain induced by phase transformation, NiTi alloys are very ductile even in the martensite phase. The purpose of the present paper is to study the influence of permanent deformation, which results from plastic deformation of martensite, on the mechanical behaviour of pseudoelastic NiTi alloys. Based on phenomenological theory of martensitic transformation and crystal plasticity, a new three dimensional micromechanical model is proposed by coupling both the slip and twinning deformation mechanisms. The present model is implemented as User MATerial subroutine (UMAT) into ABAQUS/Standard to study the influences of plastic deformation on the stress and strain fields, and on the evolution of martensite transformation. Results show that with the increasing of plastic deformation the residual strain increases and the phase transformation stress–strain curves from the martensite to austenite become steeper and less obvious. Both characteristics, stabilisation of martensite and impedance of the reverse transformation, due to plastic deformation are captured.     

Modelling and experimental characterisation of the Lüders strain in complex loaded ferritic steel compact tension specimens

The use of 3D digital image correlation (DIC) has been used to capture the Lüders strains in a low carbon ferritic steel. Results were used to calibrate and compare with finite element (FE) results based on a constitutive plasticity model, capable of yield drop behaviour and therefore Lüders strains, by Zhang et al. (2001). Tensile tests were carried out at several strain rates to characterise the material behaviour. The results of these tests were used to fit parameters in the constitutive plasticity model. The FE model was then tested on a complex loading situation of in-plane compression of a compact tension (CT) specimen. The FE model predicts the shape and formation of the Lüders bands well. This FE model, using Zhang’s constitutive plasticity model, was used to predict the residual stress profile to compare with standard elastic–plastic isotropic hardening models with no yield point. The yield point reduced both the predicted peak tensile stress, at the notch root, and the amount of plastic strain. In regions where the plastic strain was of a similar size to the Lüders strain the stress profiles were perturbed from flat profiles predicted by the standard elastic–plastic hardening models.     

Déformation cumulée tensorielle dans le référentiel en rotation logarithmiqueCumulated tensorial strain measure in the logarithmic rotating frame.

The time integration of strain rate tensor $\stackrel{?}{\mathit{D}}$ is a central problem in large transformations even if it is often an underlying one. The cumulated tensorial strains, obtained by the time integration of strain rate tensor $\stackrel{?}{\mathit{D}}$, allow the tackling of this problem from a geometrical point of view, and independently of material behaviour considerations. The time integration here takes place in the local objective frame defined by the logarithmic spin proposed by Lehmann et al. and Xiao et al. The numerical results obtained in a closed deformation path are presented here. The advantages and drawbacks of this novel integration for the development of behaviour laws are described. To cite this article: V. Mora et al., C. R. Mecanique 332 (2004).     

An evaluation of the Larson-Doi model for liquid crystalline polymers using recoil

The mesoscopic models for the rheological properties of liquid crystalline polymers proposed by Larson and Doi in 1991 and Kawaguchi and Denn in 1999 are based on phenomenological expressions that describe the evolution of the defect density and the contribution of the “texture” to the stress. In the present work, we attempt to assess some of these assumptions by monitoring how the energy stored in the texture of liquid crystalline materials evolves during shear flows. For that purpose, strain recovery is measured as a function of the applied strain for flow reversal and intermittent flow. Solutions of poly-benzylglutamate in m-cresol, hydroxypropylcellulose in water and a nematic surfactant solution are used as model systems. Although the behaviour is described qualitatively by the model, discrepancies between the predictions and the experiments are observed, especially when the shear history includes rest periods. Received: 14 July 1999 /Accepted: 30 August 1999     

A variational formulation for the incremental homogenization of elasto-plastic composites

This work addresses the micro–macro modeling of composites having elasto-plastic constituents. A new model is proposed to compute the effective stress–strain relation along arbitrary loading paths. The proposed model is based on an incremental variational principle (Ortiz, M., Stainier, L., 1999. The variational formulation of viscoplastic constitutive updates. Comput. Methods Appl. Mech. Eng. 171, 419–444) according to which the local stress–strain relation derives from a single incremental potential at each time step. The effective incremental potential of the composite is then estimated based on a linear comparison composite (LCC) with an effective behavior computed using available schemes in linear elasticity. Algorithmic elegance of the time-integration of J₂ elasto-plasticity is exploited in order to define the LCC. In particular, the elastic predictor strain is used explicitly. The method yields a homogenized yield criterion and radial return equation for each phase, as well as a homogenized plastic flow rule. The predictive capabilities of the proposed method are assessed against reference full-field finite element results for several particle-reinforced composites.     

Necking in glassy polymers: Effects of intrinsic anisotropy and structural evolution kinetics in their viscoplastic flow

Two recently proposed developments of the Glass–Rubber constitutive model for glassy polymers treat the viscoplastic deformation as intrinsically anisotropic, and incorporate the kinetics of structural evolution. These features enable the model to capture better the distinctive features of glassy polymers’ constitutive response: post-yield strain-softening and strain-hardening and effects of pre-existing molecular orientation. They have been combined to form a new variant of the model, and the consequences for necking have been explored. Uniaxial extension of prismatic bars was simulated using the finite element method, employing a numerical implementation of the new model, with material parameters of polystyrene. Strain localization predicted with the new model was found to be systematically retarded as compared to predictions with the original (intrinsically isotropic) version of the model, for the same conditions. In particular, the effect of frozen-in molecular orientation was examined. This was found to retard strain localization for stretching parallel to the orientation direction, for both models. But the localization predicted with the new model was always significantly less pronounced than with the original model. Indeed, for sufficiently high pre-orientation (e.g. a uniaxial stretch of 2.2), localization could be effectively prevented with the new model, under conditions when otherwise failure by necking is predicted. Such results can all be explained in terms of a linear stability analysis. They suggest that all previous simulations of necking in glassy polymers made using intrinsically isotropic representations of polymer viscoplasticity may have over-predicted the rate of strain localization.     

A coupled elastoplastic-damage constitutive model with Lode angle dependent failure criterion

A coupled elastoplastic-damage constitutive model with Lode angle dependent failure criterion for high strain and ballistic applications is presented. A Lode angle dependent function is added to the equivalent plastic strain to failure definition of the Johnson–Cook failure criterion. The weakening in the elastic law and in the Johnson–Cook-like constitutive relation implicitly introduces the Lode angle dependency in the elastoplastic behaviour. The material model is calibrated for precipitation hardened Inconel 718 nickel-base superalloy. The combination of a Lode angle dependent failure criterion with weakened constitutive equations is proven to predict fracture patterns of the mechanical tests performed and provide reliable results. Additionally, the mesh size dependency on the prediction of the fracture patterns was studied, showing that was crucial to predict such patterns.     

Micromechanical quantification of elastic,twinning, and slip strain partitioning exhibited by polycrystalline,monoclinic nickel–titanium during large uniaxial deformations measured via in-situ neutron diffraction

We draw upon existing knowledge of twinning and slip mechanics to develop a diffraction analysis model that allows for empirical quantification of individual deformation mechanisms to the macroscopic behaviors of low symmetry and phase transforming crystalline solids. These methods are applied in studying elasticity, accommodation twinning, deformation twinning, and slip through neutron diffraction data of tensile and compressive deformations of monoclinic NiTi to ~18% true strain. A deeper understanding of tension–compression asymmetry in NiTi is gained by connecting crystallographic calculations of polycrystalline twinning strains with in situ diffraction measurements. Our analyses culminate in empirical, micromechanical quantification of individual elastic, accommodation twinning, deformation twinning, and slip contributions to the total macroscopic stress–strain response of a monoclinic material subjected to large deformations. From these results, we find that 20–40% of the total plastic response at high strains is due to deformation twinning and 60–80% due to slip.