Constitutive framework optimized for myocardium and other high-strain, laminar materials with one fiber family

Central to this analysis is the identification of six rotation invariant scalars α_1-6 that succinctly define the strain in materials that have one family of parallel fibers arranged in laminae. These scalars were chosen so as to minimize covariance amongst the response terms in the hyperelastic limit, and they are termed strain attributes because it is necessary to distinguish them from strain invariants. The Cauchy stress t is expressed as the sum of six response terms, almost all of which are mutually orthogonal for finite strain (i.e. 14 of the 15 inner products vanish). For small deformations, the response terms are entirely orthogonal (i.e. all 15 inner products vanish). A response term is the product of a response function with its associated kinematic tensor. Each response function is a scalar partial derivative of the strain energy W with respect to a strain attribute. Applications for this theory presently include myocardium (heart muscle) which is often modeled as having muscle fibers arranged in sheets. Utility for experimental identification of strain energy functions is demonstrated by showing that common tests on incompressible materials can directly determine terms in W. Since the described set of strain attributes reduces the covariance amongst response terms, this approach may enhance the speed and precision of inverse finite element methods.     

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 simple framework for full-network hyperelasticity and anisotropic damage

A formulation of a constitutive behaviour law is proposed for hyperelastic materials, such that damage induced anisotropy can be accounted for continuously. The full-network approach with directional damage is adopted as a starting point. The full-network law with elementary strain energy density based on the inverse Langevin is chosen as a reference law which is cast into the proposed framework. This continuum formalism is then rewritten using spherical harmonics to capture damage directionality. The proposed formalism allows for an efficient (and systematic) expansion of complex non-linear anisotropic constitutive laws. A low order truncated expression of the resulting behaviour is shown to reproduce accurately the stress-strain curves of the exact behaviour laws.     

Homogenization-based constitutive models for porous elastomers and implications for macroscopic instabilities: II—Results

In Part I of this paper, we developed a homogenization-based constitutive model for the effective behavior of isotropic porous elastomers subjected to finite deformations. In this part, we make use of the proposed model to predict the overall response of porous elastomers with (compressible and incompressible) Gent matrix phases under a wide variety of loading conditions and initial values of porosity. The results indicate that the evolution of the underlying microstructure—which results from the finite changes in geometry that are induced by the applied loading—has a significant effect on the overall behavior of porous elastomers. Further, the model is in very good agreement with the exact and numerical results available from the literature for special loading conditions and generally improves on existing models for more general conditions. More specifically, we find that, in spite of the fact that Gent elastomers are strongly elliptic materials, the constitutive models for the porous elastomers are found to lose strong ellipticity at sufficiently large compressive deformations, corresponding to the possible onset of “macroscopic” (shear band-type) instabilities. This capability of the proposed model appears to be unique among theoretical models to date and is in agreement with numerical simulations and physical experience. The resulting elliptic and non-elliptic domains, which serve to define the macroscopic “failure surfaces” of these materials, are presented and discussed in both strain and stress space.     

Homogenization-based constitutive models for porous elastomers and implications for macroscopic instabilities: I—Analysis

The purpose of this paper is to provide homogenization-based constitutive models for the overall, finite-deformation response of isotropic porous rubbers with random microstructures. The proposed model is generated by means of the “second-order” homogenization method, which makes use of suitably designed variational principles utilizing the idea of a “linear comparison composite.” The constitutive model takes into account the evolution of the size, shape, orientation, and distribution of the underlying pores in the material, resulting from the finite changes in geometry that are induced by the applied loading. This point is key, as the evolution of the microstructure provides geometric softening/stiffening mechanisms that can have a very significant effect on the overall behavior and stability of porous rubbers. In this work, explicit results are generated for porous elastomers with isotropic, (in)compressible, strongly elliptic matrix phases. In spite of the strong ellipticity of the matrix phases, the derived constitutive model may lose strong ellipticity, indicating the possible development of shear/compaction band-type instabilities. The general model developed in this paper will be applied in Part II of this work to a special, but representative, class of isotropic porous elastomers with the objective of exploring the complex interplay between geometric and constitutive softening/stiffening in these materials.     

On the r-value of textured sheet metals

The r-value of a sheet metal is a measure of plastic anisotropy frequently used for prediction of performance in deep-drawing. It has also figured prominently in the literature for validation of theories where the predicted angular dependence of r is compared with the measured dependence. As plastic anisotropy in sheet metals is caused mainly by the preferred orientations of grains within the polycrystalline metal, it is natural to ask how r would depend on the orientation distribution function (ODF) w which defines the crystallographic texture of the polycrystal. In this paper a general formula relating r to w is derived for textured sheet metals whose plastic flow behavior is governed by a plastic potential f(σ, w), the anisotropic part of which depends linearly on the texture coefficients; here σ denotes the deviator of the Cauchy stress. Specific forms of this formula for orthorhombic sheets of cubic and of hexagonal metals are explicitly given.     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

A general hyperelastic model for incompressible fiber-reinforced elastomers

This work presents a new constitutive model for the effective response of fiber-reinforced elastomers at finite strains. The matrix and fiber phases are assumed to be incompressible, isotropic, hyperelastic solids. Furthermore, the fibers are taken to be perfectly aligned and distributed randomly and isotropically in the transverse plane, leading to overall transversely isotropic behavior for the composite. The model is derived by means of the “second-order” homogenization theory, which makes use of suitably designed variational principles utilizing the idea of a “linear comparison composite.” Compared to other constitutive models that have been proposed thus far for this class of materials, the present model has the distinguishing feature that it allows consideration of behaviors for the constituent phases that are more general than Neo-Hookean, while still being able to account directly for the shape, orientation, and distribution of the fibers. In addition, the proposed model has the merit that it recovers a known exact solution for the special case of incompressible Neo-Hookean phases, as well as some other known exact solutions for more general constituents under special loading conditions.     

Three-dimensional micromechanical modeling of voided polymeric materials

A three-dimensional micromechanical unit cell model for particle-filled materials is presented. The cell model is based on a Voronoi tessellation of particles arranged on a body-centered cubic (BCC) array. The three-dimensionality of the present cell model enables the study of several deformation modes, including uniaxial, plane strain and simple shear deformations, as well as arbitrary principal stress states.The unit cell model is applied to studies on the micromechanical and macromechanical behavior of rubber-toughened polycarbonate. Different load cases are examined, including plane strain deformation, simple shear deformation and principal stress states. For a constant macroscopic strain rate, the different load cases show that the macroscopic flow strength of the blend decreases with an increase in void volume fraction, as expected. The main mechanism for plastic deformation is broad shear banding across inter-particle ligaments. The distributed nature of plastic straining acts to reduce the amount of macroscopic strain softening in the blend as the initial void volume fraction is increased. In the case of plane strain deformation, the plastic flow is observed to initiate across inter-particle ligaments in the direction of constraint. This particular mode of deformation could not have been captured using a two-dimensional, plane strain idealization of cylindrical voids in a matrix.The potential for localized crazing and/or cavitation in the matrix is addressed. It is observed that the introduction of voids acts to relieve hydrostatic stress in the matrix material, compared to the homopolymer. It is also seen that the predicted peak hydrostatic stress in the matrix is higher under plane strain deformation than under triaxial tension (with equal lateral stresses), for the same macroscopic stress triaxiality.The effect of void volume fraction on the macroscopic uniaxial tension behavior of the different blends is examined using a Considère construction for dilatant materials. The natural draw ratio was predicted to decrease with an increase in void volume fraction.     

On a geometrically nonlinear damage model based on a multiplicative decomposition of the deformation gradient and the propagation of microcracks

We aim to derive a damage model for materials damaged by microcracks. The evolution of the cracks shall be governed by the maximum energy release rate, which was recently shown to be a direct consequence of the variational principle of a body with a crack (Arch. Appl. Mech. 69 (5) (1999) 337). From this, we get the path of the growing crack by introducing a series of thermodynamically equivalent straight cracks. The equivalence of the energy dissipated by microcrack growth and the damage dissipation leads to our damage evolution law. This evolution law will be embedded in a finite deformation framework based on a multiplicative decomposition into elastic and damage parts. As a consequence of this, we can present the anisotropic damaged elasticity tensor with the help of push and pull operations. The connection of this approach to other well known damage theories will be shown and the advantages of a finite element framework will be worked out. Numerical examples show the possibilities of the proposed model.     

Modeling transients in the mechanical response of copper due to strain path changes

When copper is deformed to large strains its texture and microstructure change drastically, leading to plastic anisotropy and extended transients when it is reloaded along a different strain path. For predicting these transients, we develop a constitutive model for polycrystalline metals that incorporates texture and grain microstructure. The directional anisotropy in the single crystals is considered to be induced by variable latent hardening associated with cross-slip, cut-through of planar dislocation walls, and dislocation-based reversal mechanisms. These effects are introduced in a crystallographic hardening model which is, in turn, implemented into a polycrystal model. This approach successfully explains the flow response of OFHC Cu pre-loaded in tension (compression) and reloaded in tension (compression), and the response of OFHC Cu severely strained in shear by equal channel angular extrusion and subsequently compressed in each of the three orthogonal directions. This new theoretical framework applies to arbitrary strain path changes, and is fully anisotropic.     

A micromechanical model of hardening, rate sensitivity and thermal softening in BCC single crystals

The present paper is concerned with the development of a micromechanical model of the hardening, rate-sensitivity and thermal softening of bcc crystals. In formulating the model, we specifically consider the following unit processes: double-kink formation and thermally activated motion of kinks; the close-range interactions between primary and forest dislocations, leading to the formation of jogs; the percolation motion of dislocations through a random array of forest dislocations introducing short-range obstacles of different strengths; dislocation multiplication due to breeding by double cross-slip; and dislocation pair annihilation. The model is found to capture salient features of the behavior of Ta crystals such as: the dependence of the initial yield point on temperature and strain rate; the presence of a marked stage I of easy glide, specially at low temperatures and high strain rates; the sharp onset of stage II hardening and its tendency to shift towards lower strains, and eventually disappear, as the temperature increases or the strain rate decreases; the parabolic stage II hardening at low strain rates or high temperatures; the stage II softening at high strain rates or low temperatures; the trend towards saturation at high strains; the temperature and strain-rate dependence of the saturation stress; and the orientation dependence of the hardening rate.     

Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheets: Application to sheet springback

The constitutive model for the unusual asymmetric hardening behavior of magnesium alloy sheet presented in a companion paper (Lee, M.G., Wagoner, R.H., Lee, J.K., Chung, K., Kim, H.Y., 2008. Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheet, Int. J. Plasticity 24(4), 545–582) was applied to the springback prediction in sheet metal forming. The implicit finite element program ABAQUS was utilized to implement the developed constitutive equations via user material subroutine. For the verification purpose, the springback of AZ31B magnesium alloy sheet was measured using the unconstrained cylindrical bending test of Numisheet (Numisheet ’2002 Benchmark Problem, 2002. In: Yang, D.Y., Oh, S.I., Huh, H., Kim, Y.H. (Eds.), Proceedings of 5th International Conference and Workshop on Numerical Simulation of 3D Sheet Forming Processes, Jeju, Korea) and 2D draw bend test. With the specially designed draw bend test the direct restraining force and long drawn distance were attainable, thus the measurement of the springback could be made with improved accuracy comparable with conventional U channel draw bend test. Besides the developed constitutive models, other models based on isotropic constitutive equations and the Chaboche type kinematic hardening model were also considered. Comparisons were made between simulated results by the finite element analysis and corresponding experiments and the newly proposed model showed enhanced prediction capability, which was also supported by the simple bending analysis adopting asymmetric stress–strain response.     

Experimental investigation and modeling of non-monotonic creep behavior in polymers

The deformation behavior of two unfilled engineering thermoplastics, ultra high molecular weight polyethylene (UHMWPE) and polycarbonate (PC), has been investigated in creep test conditions. It has been found that a loading history (prior to the creep test) comprising of loading to a maximum stress or strain value followed by partial unloading to arrive at the target stress value can greatly modify the strain-time behavior. Under such a test protocol, while the expected increase in strain during creep (constant tensile load) is observed, at relatively low creep stresses specimens have also demonstrated a monotonic decrease in strain. In an intermediate stress range, specimens have demonstrated time dependent behavior comprising of a transition from decreasing to increasing strain during creep in tension. This paper presents experimental results to delineate these findings and explore the effect of prior strain rate on the qualitative and quantitative changes in the output (strain-time) behavior. Furthermore, modification of the viscoplasticity theory based on overstress (VBO) model into a double element configuration is introduced. These changes confer upon the model the ability to yield non-monotonic behavior in creep, and supporting simulation results have been included. These changes, therefore, allow the model to simulate strain rate sensitivity, creep, relaxation, and recovery behavior, but more importantly address the issue of non-monotonic changes in creep and relaxation when a loading history involves some degree of unloading.     

Modeling aspects of low plastic strain amplitude multiaxial cyclic plasticity in conventional and ultrafine grain nickel

Aspects of the cyclically saturated responses of initially annealed, conventional grain size (average grain diameter of approximately 50 μm) and electrodeposited, ultrafine grain (grains from 20 to 500 nm) nickel to reversed proportional and 90° out-of-phase axial-torsional, strain-controlled cycling at a nominally constant equivalent inelastic strain amplitude of approximately 100 μ strain are reported. An anisotropic, axial-torsional subspace version of Abdel-Karim and Ohno’s kinematic hardening model is presented. Within the framework of conventional small strain, rate-independent plasticity, this approach is used to model the responses. An anhysteretic, phenomenolically based, magnetomechanical model is coupled to the rate-independent plasticity model to include the cyclic magnetostriction response. The kinematic hardening parameter determination scheme, using the proportional path responses, is described. The model correlations achieved are presented and the ability of the resulting models to capture the 90° out-of-phase responses is examined. The model parameter sets, as determined from the proportional responses, require small changes to result in more accurate correlation of the 90° out-of-phase responses and the implications of this are discussed. The relative values of the model parameters between the two materials reflect the initial microstructures. Persistent mean stresses associated with mean total strains imposed are successfully modeled for the proportional strain path responses but not for the 90° out-of-phase responses.     

On the shear and bending of a degrading polymer beam

In this paper we develop a model, within a general framework that has been developed to describe the response of dissipative systems, for the strain induced degradation of polymeric solids, due to scission. The theory can be generalized to include degradation due to ultraviolet radiation, oxygen diffusion etc., by incorporating an appropriate form for the rate of dissipation associated with these processes. We study the simple shear and pure bending of such degrading polymer beams.