Characterization of hysteretic stress-strain behavior using the integrated Preisach density

In this paper, we introduce the concept of Integrated Preisach-Mayergoyz (IPM) density to analyze static uniaxial compression tests at values well below the critical strength, and to characterize the elasticity of materials with hysteresis in their stress-strain relationship. The IPM density can be deduced from a particular force protocol following basic data treatment. The advantage of the IPM density over prior approaches is that no second order differentiation of the data is required which reduces the errors and uncertainties typical for past practice in the specific context of rock elasticity using scanning curves and PM density analysis. The characterization of the elasticity of the material is established in terms of a non-hysteretic strain contribution in the form of a non-linear but reversible equation of state, and a hysteretic contribution represented by the IPM density. The IPM inversion procedure is tested for simulated stress-strain data subjected to additive noise, and the results are compared to the traditional methodology. In addition, we analyze the hysteretic and non-hysteretic characteristics of five natural building stones, and show evidence for a classification based on the inferred properties.     

Prediction of the softening and damage effects with permanent set in fibrous biological materials

Deformation induced softening is an inelastic phenomenon frequently accompanying mechanical response of soft biological tissues. Inelastic phenomena which occur in mechanical testing of biological tissues are very likely to be associated with alterations in the internal structure of these materials.In this study, a novel structural constitutive model is formulated to describe the inelastic effects in soft biological tissues such as Mullins type behavior, damage and permanent set as a result of residual strains after unloading. Anisotropic softening is considered by evolution of internal variables governing the anisotropic properties of the material. We consider two weight factors w_i (softening) and s_k (discontinuous damage) as internal variables characterizing the structural state of the material. Numerical simulations of several soft tissues are used to demonstrate the performance of the model in reproducing the inelastic behavior of soft biological tissues.     

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.     

Multi-axial behavior of shape-memory alloys undergoing martensitic reorientation and detwinning

Recently, a rate-independent, finite-deformation-based crystal mechanics constitutive model for martensitic reorientation and detwinning in shape-memory alloys has been developed by Thamburaja [Thamburaja, P., 2005. Constitutive equations for martensitic reorientation and detwinning in shape-memory alloys. Journal of the Mechanics and Physics of Solids 53, 825–856] and implemented in the ABAQUS/Explicit [Abaqus reference manuals. 2005. Providence, RI] finite-element program. In this work, we show that the aforementioned model is able to quantitatively predict the experimental response of an initially textured and martensitic polycrystalline Ti–Ni rod under a variety of uniaxial and multi-axial stress states. By fitting the material parameters in the model to the stress–strain response in simple tension, the constitutive model predicts the stress–strain curves for experiments conducted under simple compression, torsion, proportional-loading tension–torsion, and path-change tension–torsion loading conditions to good accord. Furthermore the constitutive model also reproduces the force–displacement response for an indentation experiment to reasonable accuracy.     

Creep and recovery of polypropylene/carbon nanotube composites

The creep and recovery of polypropylene/multi-walled carbon nanotube composites were studied. It was found for thermoplastics in general that the creep strain reduces with decreased temperature, and with enhanced content of carbon nanotubes. The incorporation of nanotubes improved the recovery property remarkably, especially at high temperature. The unrecovered creep strain of nanocomposites with content of 1 and 2.8 vol.% carbon nanotubes decreased by 53% and 73% compared to that of polymer matrix. To understand the mechanisms, the Burger’s model and Weibull distribution function were employed since the variations in the simulating parameters illustrated the influence of nano-fillers on the creep and recovery performance of the bulk matrix. To further study the recovery properties, the particular contribution of each Burger’s element to the total deformation was obtained and the recovery percentage was calculated. The time-temperature-superposition-principle was applied to predict the long-term creep behavior.     

Model identification and FE simulations: Effect of different yield loci and hardening laws in sheet forming

The bi-axial experimental equipment [Flores, P., Rondia, E., Habraken, A.M., 2005a. Development of an experimental equipment for the identification of constitutive laws (Special Issue). International Journal of Forming Processes] developed by Flores enables to perform Bauschinger shear tests and successive or simultaneous simple shear tests and plane strain tests. Flores investigates the material behavior with the help of classical tensile tests and the ones performed in his bi-axial machine in order to identify the yield locus and the hardening model. With tests performed on one steel grade, the methods applied to identify classical yield surfaces such as [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic materials. Proceedings of the Royal Society of London A 193, 281–297; Hosford, W.F., 1979. On yield loci of anisotropic cubic metals. In: Proceedings of the 7th North American Metalworking Conf. (NMRC), SME, Dearborn, MI, pp. 191–197] ones as well as isotropic Swift type hardening, kinematic Armstrong–Frederick or Teodosiu and Hu hardening models are explained. Comparison with the Taylor–Bishop–Hill yield locus is also provided. The effect of both yield locus and hardening model choices is presented for two applications: plane strain tensile test and Single Point Incremental Forming (SPIF).     

Thermally actuated shape-memory polymers: Experiments, theory, and numerical simulations

With the aim of developing a thermo-mechanically coupled large-deformation constitutive theory and a numerical-simulation capability for modeling the response of thermally actuated shape-memory polymers, we have (i) conducted large strain compression experiments on a representative shape-memory polymer to strains of approximately unity at strain rates of 10⁻³ and 10⁻¹ s⁻¹, and at temperatures ranging from room temperature to approximately 30 °C above the glass transition temperature of the polymer; (ii) formulated a thermo-mechanically coupled large-deformation constitutive theory; (iii) calibrated the material parameters appearing in the theory using the stress-strain data from the compression experiments; (iv) numerically implemented the theory by writing a user-material subroutine for a widely used finite element program; and (v) conducted representative experiments to validate the predictive capability of our theory and its numerical implementation in complex three-dimensional geometries. By comparing the numerically predicted response in these validation simulations against measurements from corresponding experiments, we show that our theory is capable of reasonably accurately reproducing the experimental results. As a demonstration of the robustness of the three-dimensional numerical capability, we also show results from a simulation of the shape-recovery response of a stent made from the polymer when it is inserted in an artery modeled as a compliant elastomeric tube.     

Visco-plastic constitutive modeling on Ohno-Wang kinematic hardening rule for uniaxial ratcheting behavior of Z2CND18.12N steel

Experimental results of monotonic uniaxial tensile tests at different strain rates and the reversed strain cycling test showed the characteristics of rate-dependence and cyclic hardening of Z2CND18.12N austenitic stainless steel at room temperature, respectively. Based on the Ohno-Wang kinematic hardening rule, a visco-plastic constitutive model incorporated with isotropic hardening was developed to describe the uniaxial ratcheting behavior of Z2CND18.12N steel under various stress-controlled loading conditions. Predicted results of the developed model agreed better with experimental results when the ratcheting strain level became higher, but the developed model overestimated the ratcheting deformation in other cases. A modified model was proposed to improve the prediction accuracy. In the modified model, the parameter m_i of the Ohno-Wang kinematic hardening rule was developed to evolve with the accumulated plastic strain. Simulation results of the modified model proved much better agreement with experiments.     

An updated version of the multimechanism model for cyclic plasticity

In this paper we consider the elastoplastic behavior of the 304L stainless steel under cyclic loading at room temperature. After the experimental investigations presented in Taleb and Hauet (2009), the present work deals with modeling in the light of the new observations. An improved version of the multimechanism model is proposed in which the isotropic variable is revisited in order to take into account the non-proportional effect of the loading as well as the strain memory phenomenon. A particular attention has been paid to the identification process in order to capture the main important phenomena: relative parts of isotropic and kinematic hardening, time dependent effects, non-proportionality effect, strain amplitude dependence. Only strain controlled tests have been used for the identification process. The capabilities of the model with “only” 17 parameters are evaluated considering a number of proportional and non-proportional stress and strain controlled tests.     

Constitutive equations for martensitic reorientation and detwinning in shape-memory alloys

A crystal-inelasticity-based constitutive model for martensitic reorientation and detwinning in shape-memory alloys (SMAs) has been developed from basic thermodynamics principles. The model has been implemented in a finite-element program by writing a user-material subroutine. We perform two sets of finite-element simulations to model the behavior of polycrystalline SMAs: (1) The full finite-element model where each finite element represents a collection of martensitic microstructures which originated from within an austenite single crystal, chosen from a set of crystal orientations that approximates the initial austentic crystallographic texture. The macroscopic stress-strain responses are calculated as volume averages over the entire aggregate: (2) The Taylor model (J. Inst. Metals 62 (1938) 32) where an integration point in a finite element represents a material point which consist of sets of martensitic microstructures which originated from within respective austenite single-crystals. Here the macroscopic stress-strain responses are calculated through a homogenization scheme.Experiments in tension and compression were conducted on textured polycrystalline Ti-Ni rod initially in the martensitic phase by Xie et al (Acta Mater. 46 (1998) 1989). The material parameters for the constitutive model were calibrated by fitting the tensile stress-strain response from a full finite-element calculation of a polycrystalline aggregate to the simple tension experiment. With the material parameters calibrated the predicted stress-strain curve for simple compression is in very good accord with the corresponding experiment. By comparing the simulated stress-strain response in simple tension and simple compression it is shown that the constitutive model is able to predict the observed tension-compression asymmetry exhibited by polycrystalline Ti-Ni to good accuracy. Furthermore, our calculations also show that the macroscopic stress-strain response depends strongly on the initial martensitic microstructure and crystallographic texture of the material.We also show that the Taylor model predicts the macroscopic stress-strain curves in simple tension and simple compression reasonably well. Therefore, it may be used as a relatively inexpensive computational tool for the design of components made from shape-memory materials.     

Tension–compression asymmetry of phosphor bronze for electronic parts and its effect on bending behavior

In-plane tension and compression experiments on copper alloy sheets (phosphor bronze) and 6000 series aluminum alloy sheets (AA6016-T4) were conducted using a specially designed testing apparatus. The apparatus is equipped with comb-type dies so that stress–strain curves of a sheet specimen subjected to tension followed by compression, and vice versa, can be measured without buckling of the specimen, as well as those for monotonic tension and compression. A difference was observed in the flow stresses between tension and compression for the as-received copper alloy, but not for the aluminum alloy. Moreover, stress reversal tests, such as tension followed by compression and compression followed by tension, were carried out in order to measure the Bauschinger effect. In the second part of the experiment, bending moment–curvature diagrams were measured for the as-received and pre-stretched specimens. The bending moment–curvature diagrams were compared with those calculated using the stress–strain curves obtained from the tension–compression tests, and were in good agreement with those calculated with the tension–compression asymmetry and the Bauschinger effect correctly reproduced.     

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.     

Micromechanical modeling of stress-induced phase transformations. Part 2. Computational algorithms and examples

Based on the theory developed in Part 1 of this paper [Levitas, V.I., Ozsoy, I.B., 2008. Micromechanical modeling of stress-induced phase transformations. Part 1. Thermodynamics and kinetics of coupled interface propagation and reorientation. Int. J. Plasticity. doi:10.1016/j.ijplas.2008.02.004], various non-trivial examples of microstructure evolution under complex multiaxial loading are presented. For the case without interface rotation, the effect of the athermal thresholds for austenite (A)–martensite (M) and martensitic variant M_I–variant M_II interfaces and loading paths on stress–strain curves and phase transformations was studied. For coupled interface propagation and rotation, two types of numerical simulations were carried out. The tetragonal–orthorhombic transformation has been studied under general three-dimensional interface orientation and zero athermal threshold. The cubic–tetragonal transformation was treated with allowing for an athermal threshold and interface reorientation within a plane. The effect of the athermal threshold, the number of martensitic variants and an interface orientation in the embryo was studied in detail. It was found that an instability in the interface normal leads to a jump-like interface reorientation that has the following features of the energetics of a first-order transformation: there are multiple energy minima versus interface orientation that are separated by an energy barrier; positions of minima do not change during loading but their depth varies; when the barrier disappears (i.e. one of the minima transforms to the local saddle or maximum points), the system rapidly evolves toward another stable orientation. Depending on the loading and material parameters, we observed a large continuous change in interface orientation, a jump in interface reorientation, a jump in volume fractions and stresses, an expected stress relaxation during the phase transition and unexpected stress growth during the transition because of large change in elastic moduli.     

Formulation of a damage internal state variable model for amorphous glassy polymers

The following article proposes a damage model that is implemented into a glassy, amorphous thermoplastic thermomechanical inelastic internal state variable framework. Internal state variable evolution equations are defined through thermodynamics, kinematics, and kinetics for isotropic damage arising from two different inclusion types: pores and particles. The damage arising from the particles and crazing is accounted for by three processes of damage: nucleation, growth, and coalescence. Nucleation is defined as the number density of voids/crazes with an associated internal state variable rate equation and is a function of stress state, molecular weight, fracture toughness, particle size, particle volume fraction, temperature, and strain rate. The damage growth is based upon a single void growing as an internal state variable rate equation that is a function of stress state, rate sensitivity, and strain rate. The coalescence internal state variable rate equation is an interactive term between voids and crazes and is a function of the nearest neighbor distance of voids/crazes and size of voids/crazes, temperature, and strain rate. The damage arising from the pre-existing voids employs the Cocks–Ashby void growth rule. The total damage progression is a summation of the damage volume fraction arising from particles and pores and subsequent crazing. The modeling results compare well to experimental findings garnered from the literature. Finally, this formulation can be readily implemented into a finite element analysis.     

A simple nonlinear model to simulate the localized necking and neck propagation

This paper deals with the equilibrium problem in nonlinear dissipative inelasticity of damaged bodies subject to uniaxial loading and its main purpose is to show the interesting potentialities offered by the damage theory in modeling the necking and neck propagation phenomena in polymeric materials. In detail, the proposed mechanical model is a two-phase system, with the same constitutive law but with different levels of damage for each phase. Despite its simplicity, it is shown that the model can straightforwardly reproduce the overall load–elongation curve provided by experimental tensile tests by involving only five parameters of clear physical meaning.     

Multi-scale defect interactions in high-rate brittle material failure. Part I: Model formulation and application to ALON

Within this two part series we develop a new material model for ceramic protection materials to provide an interface between microstructural parameters and bulk continuum behavior to provide guidance for materials design activities. Part I of this series focuses on the model formulation that captures the strength variability and strain rate sensitivity of brittle materials and presents a statistical approach to assigning the local flaw distribution within a specimen. The material model incorporates a Mie–Grüneisen equation of state, micromechanics based damage growth, granular flow and dilatation of the highly damaged material, and pore compaction for the porosity introduced by granular flow. To provide initial qualitative validation and illustrate the usefulness of the model, we use the model to investigate Edge on Impact experiments (Strassburger, 2004) on Aluminum Oxynitride (AlON), and discuss the interactions of multiple mechanisms during such an impact event. Part II of this series is focused on additional qualitative validation and using the model to suggest material design directions for boron carbide.