Orthotropic elastic–plastic material model for paper materials

Paper and paperboard generally exhibit anisotropic and non-linear mechanical material behaviour. In this work, the development of an orthotropic elastic–plastic constitutive model, suitable for modelling of the material behaviour of paper is presented. The anisotropic material behaviour is introduced into the model by orthotropic elasticity and an isotropic plasticity equivalent transformation tensor. A parabolic stress–strain relation is adopted to describe the hardening of the material. The experimental and numerical procedures for evaluation of the required material parameters for the model are described. Uniaxial tensile testing in three different inplane material directions provides the calibration of the material parameters under plane stress conditions. The numerical implementation of the material model is presented and the model is shown to perform well in agreement with experimentally observed mechanical behaviour of paper.     

Dual finite element methods in homogenization for elastic–plastic fibrous composite material

The problem of homogenization for a periodic, elastic–perfectly plastic, fiber reinforced, composite material is considered. The overall mechanical behavior of the material is described using the anisotropic model of elastic–plastic body with kinematic hardening. The appropriate initial–boundary value problem, set for one repeatable cell of the composite, is solved in order to find the effective constitutive relations. The cell problem is solved using the finite element method formulated in two dual forms: in displacements and in stresses. Stress functions are used in the latter formulation.     

Thermo-mechanical coupling strategies in elastic–plastic problems

Modeling strategies aimed at thermo-mechanical coupled problems has been developed for a wide range of engineering applications. Staggered-type coupling procedures have been largely used in materials processing operations, especially in commercial codes, owing to their simplicity and flexibility. The present work shows that, in thermo-plastic problems, the classical implementation of the most common coupling procedure may present accuracy issues and time-stepping dependency. Numerical experiments indicate that an iterative coupling scheme constitutes a viable and simple approach to this class of problems.     

A consistent equilibrium in a cross-section of an elastic–plastic beam

A phenomenon of inequality of equilibrium and constitutive internal forces in a cross-section of elastic–plastic beams is common to many finite element formulations. It is here discussed in a rate-independent, elastic–plastic beam context, and a possible treatment is presented. The starting point of our discussion is Reissners finite-strain beam theory, and its finite element implementation. The questions of the consistency of interpolations for displacements and rotations, and the related locking phenomena are fully avoided by considering the rotation function of the centroid axis of a beam as the only unknown function of the problem. Approximate equilibrium equations are derived by the use of the distribution theory in conjunction with the collocation method. The novelty of our formulation is an inclusion of a balance function that measures the error between the equilibrium and constitutive bending moments in a cross-section. An advantage of the present approach is that the locations, where the balance of equilibrium and constitutive moments should be satisfied, can be prescribed in advance. In order to minimize the error, explicit analytical expressions are used for the constitutive forces; for a rectangular cross-section and bilinear constitutive law, they are given in Appendix A. The comparison between the results of the two finite element formulations, the one using consistent, and the other inconsistent equilibrium in a cross-section, is shown for a cantilever beam subjected to a point load. The problem of high curvature gradients in a plastified region is also discussed and solved by using an adapted collocation method, in which the coordinate system is transformed such to follow high gradients of curvature.     

Modelling paper as a two-dimensional elastic–plastic stochastic network

Stochastic two-dimensional elastic–plastic network models are used to represent the inelastic deformation behavior of well-bonded paper. Linear kinematic hardening is employed with an initial non-zero back stress to represent anisotropic fiber yield. Network models are used to simulate simple monotonic tension and simple cyclic tension of paper materials. The performance of the models is compared to experimental results and found to perform reasonably well. The results suggest that interfiber bonding must be explicitly accounted for to adequately describe the material. Some discrepancy between the model and experimental cyclic tension results is believed to be due to time-dependent strain recovery in the material which is not represented in the network models. Experimental results are also presented which show that simple tension failure in these materials occurs along a line of localized deformation in a majority of the samples. This line is generally observed to form immediately prior to failure and is oriented at a well-defined angle with respect to the loading direction.     

Effect of strain hardening in elastic–plastic transition behavior in a hemisphere in contact with a rigid flat

An axisymmetrical hemispherical asperity in contact with a rigid flat is modeled for an elastic–plastic material on the lines of the Kogut–Etsion Model (KE Model) and the Jackson–Green Model (JG Model). The present work extends the previous KE and JG works, accounting for the effect of realistic material behavior in terms of the varying yield strengths and the isotropic strain hardening behavior. The predicted results show that the transition behavior of the materials from the elastic–plastic to the fully plastic case is influenced by the yield strength and the tangent modulus (E_t) and such transition do not take place at specific values of interference ratios as suggested by the KE model. New empirical relations are proposed to determine the contact load and the contact area based on the analysis. Numerical results from the finite element modeling are also validated with an experimental ball on flat configuration approach.     

Steady crack growth in elastic–plastic fluid-saturated porous media

An asymptotic solution is obtained for stress and pore pressure fields near the tip of a crack steadily propagating in an elastic–plastic fluid-saturated porous material displaying linear isotropic hardening. Quasi-static crack growth is considered under plane strain and Mode I loading conditions. In particular, the effective stress is assumed to obey the Drucker–Prager yield condition with associative or non-associative flow-rule and linear isotropic hardening is adopted. Both permeable and impermeable crack faces are considered. As for the problem of crack propagation in poroelastic media, the behavior is asymptotically drained at the crack-tip. Plastic dilatancy is observed to have a strong effect on the distribution and intensity of pore water pressure and to increase its flux towards the crack-tip.     

Energy relaxation of non-convex incremental stress potentials in a strain-softening elastic–plastic bar

We propose a fundamentally new concept to the treatment of material instabilities and localization phenomena based on energy minimization principles in a strain-softening elastic–plastic bar. The basis is a recently developed incremental variational formulation of the local constitutive response for generalized standard media. It provides a quasi-hyperelastic stress potential that is obtained from a local minimization of the incremental energy density with respect to the internal variables. The existence of this variational formulation induces the definition of the material stability of inelastic solids based on convexity properties in analogy to treatments in elasticity. Furthermore, localization phenomena are understood as micro-structure development associated with a non-convex incremental stress potential in analogy to phase decomposition problems in elasticity. For the one-dimensional bar considered the two-phase micro-structure can analytically be resolved by the construction of a sequentially weakly lower semicontinuous energy functional that envelops the not well-posed original problem. This relaxation procedure requires the solution of a local energy minimization problem with two variables which define the one-dimensional micro-structure developing: the volume fraction and the intensity of the micro-bifurcation. The relaxation analysis yields a well-posed boundary-value problem for an objective post-critical localization analysis. The performance of the proposed method is demonstrated for different discretizations of the elastic–plastic bar which document on the mesh-independence of the results.     

On the plastic wave speeds in rate-independent elastic–plastic materials with anisotropic elasticity

This paper presents a theoretical study of the speeds of plastic waves in rate-independent elastic–plastic materials with anisotropic elasticity. It is shown that for a given propagation direction the plastic wave speeds are equal to or lower than the corresponding elastic speeds, and a simple expression is provided for the bound on the difference between the elastic and the plastic wave speeds. The bound is given as a function of the plastic modulus and the magnitude of a vector defined by the current stress state and the propagation direction. For elastic–plastic materials with cubic symmetry and with tetragonal symmetry, the upper and lower bounds on the plastic wave speeds are obtained without numerically solving an eigenvalue problem. Numerical examples of materials with cubic symmetry (copper) and with tetragonal symmetry (tin) are presented as a validation of the proposed bounds. The lower bound proposed here on the minimum plastic wave speed may also be used as an efficient alternative to the bifurcation analysis at early stages of plastic deformation for the determination of the loss of ellipticity.     

Spherical impression of thin elastic films on elastic–plastic substrates

The mechanical behavior of thin elastic films deposited onto structural alloys plays a critical role in determining film durability. This paper presents analysis of an impression experiment designed to evaluate some of the relevant properties of these films. The modeling provides quantitative strain information which can be used to estimate the fracture toughness of the film, the static friction coefficient of the surface and the constitutive behavior of the substrate. Results are presented for radial and circumferential strain distributions in the film relevant to the interpretation of cracking patterns. Additionally, load-displacement curves are provided that may be used to evaluate the plastic properties of the substrate. To facilitate estimates of the film cracking strain through correlation with experiments, the radial strain distributions are presented as functions of impression depth, yield strain and hardening exponent.     

On establishing elastic–plastic properties of a sphere by indentation testing

Instrumented indentation is a popular technique for determining mechanical properties of materials. Currently, the evaluation techniques of instrumented indentation are mostly limited to a flat substrate being indented by various shaped indenters (e.g., conical or spherical). This work investigates the possibility of extending instrumented indentation to non-flat surfaces. To this end, conical indentation of a sphere is investigated where two methodologies for establishing mechanical properties are explored. In the first approach, a semi-analytical approach is employed to determine the elastic modulus of the sphere utilizing the elastic unloading response (the “unloading slope”). In the second method, reverse analysis based on finite element analysis is used, where non-dimensional characteristic functions derived from the force–displacement response are utilized to determine the elastic modulus and yield strength. To investigate the accuracies of the proposed methodologies, selected numerical experiments have been performed and excellent agreement was obtained.     

Nonlinear elastic wave propagation in a phononic material with periodic solid–solid contact interface

Phononic materials enable enhanced dynamic properties, and offer the ability to engineer the material response. In this work we study the wave propagation in such a structure when introduced with nonlinearity. Our system is comprised of pre-compressed material with periodic solid–solid contacts, which exhibit a quadratic nonlinearity for small displacements. We suggest a new approach to modeling this system, where we discretize the unit cell in order to derive an approximate analytical solution using a perturbation method, which we are then able to easily validate numerically. With these methods, we study the band structure in the system and the second harmonic generation originating from the nonlinearity. We qualitatively analyze the second harmonic response of the system in terms of the single-crack response with linear band structure considerations. Significant band structure manipulation by changing system parameters is demonstrated, including possible in-situ tuning. The system also exhibits effective frequency doubling, i.e. the transmitted wave is primarily comprised of the second harmonic wave, for a certain range of frequencies. We demonstrate very high robustness to disorder in the system, in terms of band structure and second harmonic generation. These results have possible applications as frequency-converting devices, tunable engineered materials, and in non-destructive evaluation.     

A note on divergence and flutter instabilities in elastic–plastic materials

Dynamic stability of uniform straining of a nonlinear elastic solid is known to require that all eigenvalues of the acoustic tensor associated with the tangent elastic moduli be real and nonnegative. The focus of this note is to what extent this conclusion applies to time-independent, elastoplastic materials. Nonlinearity of the elastic–plastic constitutive law imposes limits on validity of a solution to the linear problem for which the acoustic tensor is determined. The effect of those limits on the conclusions about instability is examined.     

Inhomogeneous ‘longitudinal’ plane waves in a deformed elastic material

By definition, a homogeneous isotropic compressible Hadamard material has the property that an infinitesimal longitudinal homogeneous plane wave may propagate in every direction when the material is maintained in a state of arbitrary finite static homogeneous deformation. Here, as regards the wave, homogeneous means that the direction of propagation of the wave is parallel to the direction of eventual attenuation; and longitudinal means that the wave is linearly polarized in a direction parallel to the direction of propagation. In other words, the displacement is of the form u = ncos k(n · x – ct), where n is a real vector. It is seen that the Hadamard material is the most general one for which a longitudinal inhomogeneous plane wave may also propagate in any direction of a predeformed body. Here, inhomogeneous means that the wave is attenuated, in a direction distinct from the direction of propagation; and longitudinal means that the wave is elliptically polarized in the plane containing these two directions, and that the ellipse of polarization is similar and similarly situated to the ellipse for which the real and imaginary parts of the complex wave vector are conjugate semi-diameters. In other words, the displacement is of the form u = {S exp i(S · x – ct)}, where S is a complex vector (or bivector). Then a Generalized Hadamard material is introduced. It is the most general homogeneous isotropic compressible material which allows the propagation of infinitesimal longitudinal inhomogeneous plane circularly polarized waves for all choices of the isotropic directional bivector. Finally, the most general forms of response functions are found for homogeneously deformed isotropic elastic materials in which longitudinal inhomogeneous plane waves may propagate with a circular polarization in each of the two planes of central circular section of the ⁿ -ellipsoid, where is the left Cauchy-Green strain tensor corresponding to the primary pure homogeneous deformation.     

Shakedown kinematic theorem for elastic–perfectly plastic bodies

The classical shakedown kinematic theorem due to Koiter for elastic–perfectly plastic bodies is re-examined and divided into separated shakedown and nonshakedown theorems. While the shakedown theorem is based on the set of Koiter's plastic strain rate cycles, the non-shakedown one involves a broader set of admissible plastic strain rate cycles, the end-cycle accumulated strains of which are deviatoric parts of compatible strain fields. For certain broad classes of practical problems the two statements are unified to yield the unique theorem in Koiter's sense.     

Shakedown theorems for elastic–plastic solids in the framework of gradient plasticity

Static and kinematic shakedown theorems are given for a class of generalized standard materials endowed with a hardening saturation surface in the framework of strain gradient plasticity. The so-called residual-based gradient plasticity theory is employed. The hardening law admits a hardening potential, which is a C¹-continuous function of a set of kinematic internal variables and of their spatial gradients, and is required to satisfy a global sign restriction (but not to be necessarily convex). The totally produced, the accumulated and the freely moving dislocations per unit volume, distinguished as statistically stored and geometrically necessary ones, are in this way accounted for. Like for a generalized standard material, the shakedown safety factor is found to depend on the (generally size dependent) yield and saturation limits, but not on the particular differential-type hardening law of the material.     

The elastic–plastic behaviour of foam under shock loading

An experimental investigation of the elastic–plastic nature of shock wave propagation in foams was undertaken. The study involved experimental blast wave and shock tube loading of three foams, two polyurethane open-cell foams and a low-density polyethylene closed-cell foam. Evidence of precursor waves was observed in all three foam samples under various compressive wave loadings. Experiments with an impermeable membrane are used to determine if the precursor wave in an open-cell foam is a result of gas filtration or an elastic response of the foam. The differences between quasi-static and shock compression of foams is discussed in terms of their compressive strain histories and the implications for the energy absorption capacity of foam in both loading scenarios. Through a comparison of shock tube and blast wave loading techniques, suggestions are made concerning the accurate measurements of the principal shock Hugoniot in foams.     

Effect of large elastic strains on cavitation instability predictions for elastic–plastic solids

For an infinite solid containing a void, the cavitation instability limit is defined as the remote stress–and strain state, at which the void grows without bound, driven by the elastic energy stored in the surrounding material. Such cavitation limits have been analysed by a number of authors for metal plasticity as well as for nonlinear elastic solids. The analyses for elastic–plastic solids are here extended to consider the effect of a large initial yield strain, and it is shown how the critical stress value decays for increasing value of the yield strain. Analyses are carried out for remote hydrostatic tension as well as for more general axisymmetric remote stress field, with an initially spherical void. Different levels of strain hardening are considered.     

Evolution of plastic zones in dynamically loaded plates using different elastic–viscoplastic laws

In the framework of viscoplastic theory many different laws were developed, accounting for material behaviors like creep, relaxation or evolution of overstresses. Though each model is able to predict in uni-axial material tests the values of stresses depending on plastic strains and plastic strain rates the question is if solutions of simulations are still realistic if the viscoplastic law is applied on structural deformations. In the present study strain rate sensitive metal plates are subjected to shock waves. The purpose is to compare simulation results obtained with different elastic–viscoplastic laws to experiments in order to determine the most appropriate material model. By subjecting circular metal plates experimentally to shock wave loadings realistic deformation histories are measured. The measurements are compared to simulation results obtained with different viscoplastic laws. The aim is to find out the accuracy of each model concerning the predictions of displacements, shape formings, spread of plastic zones and evolutions of inner bending moments.