On phase transformations in shape memory alloy materials and large deformation generalized plasticity

A new version of rate-independent generalized plasticity, suitable for the derivation of general thermomechanical constitutive laws for materials undergoing phase transformations, is proposed within a finite deformation framework. More specifically, by assuming an additive decomposition of the finite strain tensor into elastic and inelastic (transformation induced) parts and by considering the fractions of the various material phases as internal variables, a multi-phase formulation of the theory is developed. The concepts presented are applied for the derivation of a three-dimensional thermomechanical model for shape memory alloy materials. The ability of the model in simulating several patterns of the extremely complex behavior of these materials, under both monotonic and cyclic loadings, is assessed by representative numerical examples.     

A three-dimensional constitutive model for shape memory alloys

Shape memory alloys (SMAs) are materials that, among other characteristics, have the ability to present high deformation levels when subjected to mechanical loading, returning to their original form after a temperature change. Literature presents numerous constitutive models that describe the phenomenological features of the thermomechanical behavior of SMAs. The present paper introduces a novel three-dimensional constitutive model that describes the martensitic phase transformations within the scope of standard generalized materials. The model is capable of describing the main features of the thermomechanical behavior of SMAs by considering four macroscopic phases associated with austenitic phase and three variants of martensite. A numerical procedure is proposed to deal with the nonlinearities of the model. Numerical simulations are carried out dealing with uniaxial and multiaxial single-point tests showing the capability of the introduced model to describe the general behavior of SMAs. Specifically, uniaxial tests show pseudoelasticity, shape memory effect, phase transformation due to temperature change and internal subloops due to incomplete phase transformations. Concerning multiaxial tests, the pure shear stress and hydrostatic tests are discussed showing qualitatively coherent results. Moreover, other tensile–shear tests are conducted modeling the general three-dimensional behavior of SMAs. It is shown that the multiaxial results are qualitative coherent with the related data presented in the literature.     

A constitutive theory for shape memory polymers. Part II: A linearized model for small deformations

A constitutive theory is developed for shape memory polymers. It is to describe the thermomechanical properties of such materials under large deformations. The theory is based on the idea, which is developed in the work of Liu et al. [2006. Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modelling. Int. J. Plasticity 22, 279-313], that the coexisting active and frozen phases of the polymer and the transitions between them provide the underlying mechanisms for strain storage and recovery during a shape memory cycle. General constitutive functions for nonlinear thermoelastic materials are used for the active and frozen phases. Also used is an internal state variable which describes the volume fraction of the frozen phase. The material behavior of history dependence in the frozen phase is captured by using the concept of frozen reference configuration. The relation between the overall deformation and the stress is derived by integration of the constitutive equations of the coexisting phases. As a special case of the nonlinear constitutive model, a neo-Hookean type constitutive function for each phase is considered. The material behaviors in a shape memory cycle under uniaxial loading are examined. A linear constitutive model is derived from the nonlinear theory by considering small deformations. The predictions of this model are compared with experimental measurements.     

A constitutive theory for shape memory polymers. Part I: Large deformations

A constitutive theory is developed for shape memory polymers. It is to describe the thermomechanical properties of such materials under large deformations. The theory is based on the idea, which is developed in the work of Liu et al. [2006. Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modeling. Int. J. Plasticity 22, 279-313], that the coexisting active and frozen phases of the polymer and the transitions between them provide the underlying mechanisms for strain storage and recovery during a shape memory cycle. General constitutive functions for nonlinear thermoelastic materials are used for the active and frozen phases. Also used is an internal state variable which describes the volume fraction of the frozen phase. The material behavior of history dependence in the frozen phase is captured by using the concept of frozen reference configuration. The relation between the overall deformation and the stress is derived by integration of the constitutive equations of the coexisting phases. As a special case of the nonlinear constitutive model, a neo-Hookean type constitutive function for each phase is considered. The material behaviors in a shape memory cycle under uniaxial loading are examined. A linear constitutive model is derived from the nonlinear theory by considering small deformations. The predictions of this model are compared with experimental measurements.     

Finite Element Analysis of Shape Memory Alloy Adaptive Trusses with Geometrical Nonlinearities

This contribution deals with the nonlinear analysis of shape memory alloy (SMA) adaptive trusses employing the finite element method. Geometrical nonlinearities are incorporated into the formulation together with a constitutive model that describes different thermomechanical behaviors of SMA. It has four macroscopic phases (three variants of martensite and an austenitic phase), and considers different material properties for austenitic and martensitic phases together with thermal expansion. An iterative numerical procedure based on the operator split technique is proposed in order to deal with the nonlinearities in the constitutive formulation. This procedure is introduced into ABAQUS as a user material routine. Numerical simulations are carried out illustrating the ability of the developed model to capture the general behavior of shape memory bars. After that, it is analyzed the behavior of some adaptive trusses built with SMA actuators subjected to different thermomechanical loadings.     

Three-dimensional constitutive model for shape memory alloys based on microplane model

A new model for the behavior of polycrystalline shape memory alloys (SMA), based on a statically constrained microplane theory, is proposed. The new model can predict three-dimensional response by superposing the effects of inelastic deformations computed on several planes of different orientation, thus reproducing closely the actual physical behavior of the material. Due to the structure of the microplane algorithm, only a one-dimensional constitutive law is necessary on each plane. In this paper, a simple constitutive law and a robust kinetic expression are used as the local constitutive law on the microplane level. The results for SMA response on the macroscale are promising: simple one-dimensional response is easily reproduced, as are more complex features such as stress-strain subloops and tension-compression asymmetry. A key feature of the new model is its ability to accurately represent the deviation from normality exhibited by SMAs under nonproportional loading paths.     

The decomposition F = FeFp,material symmetry,and plastic irrotationality for solids that are isotropic-viscoplastic or amorphous

This study develops a general framework for discussing both isotropic-viscoplastic materials and amorphous materials. The framework, which allows for large deformations, is based on the Kröner–Lee decomposition of the deformation gradient into elastic and inelastic parts, a system of microforces consistent with its own balance, and a mechanical version of the second law that includes, via the microforces, work performed during inelastic flow. The constitutive theory allows for dependences on the elastic and inelastic parts of the deformation gradient and on the inelastic stretch-rate, but dependences on the inelastic spin are not included. The constitutive equation for the microstress T^p conjugate to inelastic flow – suitably restricted by the second law – and the microforce balance are shown to be together equivalent to a flow rule that includes a back stress due to the variation in the free energy with inelastic deformation. The introduction of a concept of material microstability reduces this flow rule to one of classical Mises-type.In a theory based on the Kröner–Lee decomposition, there are two classes of symmetry transformations available: transformations of the reference configuration and transformations of the relaxed spaces. We discuss the notion of material symmetry for a general class of materials that includes, as special cases, isotropic-viscoplastic solids, and amorphous solids. Essential to this discussion of symmetry is a general constitutive relation for the microstress T^p.The symmetry-based framework allows us to show that for typical boundary-value problems involving isotropic, viscoplastic solids or amorphous solids, if a problem has a solution, then every time- and space-dependent rotation of the relaxed spaces also yields a solution, and it is possible to choose this rotation such that the transformed solution is inelastically spin-free: W^p ≡ 0. Thus, when discussing such materials, we may, without loss in generality, restrict attention to flow rules that are inelastically irrotational.     

Differential and integrated form consistency in 1-D phenomenological models for shape memory alloy constitutive behavior

Shape memory alloys are being explored increasingly for developing smart structures and devices in aerospace, automotive and other application areas. The material behavior is highly nonlinear with coupled thermomechanical response involving temperature and/or stress induced phase transformations. Modeling the constitutive behavior of these materials poses several challenges and a few phenomenological models exist that provide a quick and reasonable approach to assess their behavior. Due to phenomenological approach, several assumptions are made in order to simplify the model and some of them introduce inconsistencies or anomalies into the model. In this paper, a frequently used approach, namely, Brinson [Brinson, L.C., 1993. One dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable. J. Intell. Mater. Syst. Struct. 4(2), 229–242.] model, is investigated. The constitutive equation is usually expressed at the outset in the differential form and the integrated form of the same is obtained. It is shown that the two forms of equations are not consistent in the Brinson [Brinson, L.C., 1993. One dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable. J. Intell. Mater. Syst. Struct. 4(2), 229–242.] model, given the assumed form of material functions. In the present work, the nature and implications of the inconsistency are highlighted. The cause of incompatibility is the inconsistent material definitions. A modified consistent constitutive model is proposed by redefining the material function which satisfies the compatibility condition. The advantages in using the proposed modified model are highlighted with numerical case studies involving hysteretic stress–strain behavior.     

A finite strain kinematic hardening constitutive model based on Hencky strain: General framework, solution algorithm and application to shape memory alloys

The logarithmic or Hencky strain measure is a favored measure of strain due to its remarkable properties in large deformation problems. Compared with other strain measures, e.g., the commonly used Green-Lagrange measure, logarithmic strain is a more physical measure of strain. In this paper, we present a Hencky-based phenomenological finite strain kinematic hardening, non-associated constitutive model, developed within the framework of irreversible thermodynamics with internal variables. The derivation is based on the multiplicative decomposition of the deformation gradient into elastic and inelastic parts, and on the use of the isotropic property of the Helmholtz strain energy function. We also use the fact that the corotational rate of the Eulerian Hencky strain associated with the so-called logarithmic spin is equal to the strain rate tensor (symmetric part of the velocity gradient tensor). Satisfying the second law of thermodynamics in the Clausius-Duhem inequality form, we derive a thermodynamically-consistent constitutive model in a Lagrangian form. In comparison with the available finite strain models in which the unsymmetric Mandel stress appears in the equations, the proposed constitutive model includes only symmetric variables. Introducing a logarithmic mapping, we also present an appropriate form of the proposed constitutive equations in the time-discrete frame. We then apply the developed constitutive model to shape memory alloys and propose a well-defined, non-singular definition for model variables. In addition, we present a nucleation-completion condition in constructing the solution algorithm. We finally solve several boundary value problems to demonstrate the proposed model features as well as the numerical counterpart capabilities.     

COARSE-GRAINED ATOMISTIC MODELING AND SIMULATION OF INELASTIC MATERIAL BEHAVIOR

This paper presents a new methodology for coarse-grained atomistic simulation of inelastic material behavior including phase transformations in ceramics and dislocation mediated plasticity in metals. The methodology combines an atomistic formulation of balance equations and a modified finite element method. With significantly fewer degrees of freedom than those of a fully atomistic model and without additional constitutive rules but the interatomic force field, the new coarse-grained (CG) method is shown to be feasible in predicting the nonlinear constitutive responses of materials and also reproducing atomic-scale phenomena such as phase transformations (diamond →β-Sn) in silicon and dislocation nucleation and migration, formation of dislocation loops and stacking faults ribbons in single crystal nickel. Direct comparisons between CG and the corresponding full molecular dynamics (MD) simulations show that the present methodology is efficient and promising in modeling and simulation of inelastic material behavior without losing the essential atomistic features. The potential applications and the limitations of the CG method are also discussed.     

Coupled thermomechanical dynamics of phase transitions in shape memory alloys and related hysteresis phenomena

In this paper the nonlinear dynamics of shape memory alloy phase transformations is studied with thermomechanical models based on coupled systems of partial differential equations by using computer algebra tools. The reduction procedures of the original model to a system of differential-algebraic equations and its solution are based on the general methodology developed by the authors for the analysis of phase transformations in shape memory materials with low dimensional approximations derived from center manifold theory. Results of computational experiments revealing the martensitic-austenitic phase transition mechanism in a shape-memory-alloy rod are presented. Several groups of computational experiments are reported. They include results on stress- and temperature-induced phase transformations as well as the analysis of the hysteresis phenomenon. All computational experiments are presented for Cu-based structures.     

Thermodynamic and relaxation-based modeling of the interaction between martensitic phase transformations and plasticity

This paper focuses on the issue plasticity within the framework of a micromechanical model for single-crystal shape-memory alloys. As a first step towards a complete micromechanical formulation of such models, we work with classical J₂-von Mises-type plasticity for simplicity. The modeling of martensitic phase transitions is based on the concept of energy relaxation (quasiconvexification) in connection with evolution equations derived from inelastic potentials. Crystallographic considerations lead to the derivation of Bain strains characterizing the transformation kinematics. The model is derived for arbitrary numbers of martensite variants and thus can be applied to any shape-memory material such as CuAlNi or NiTi. The phase transition model captures effects like tension/compression asymmetry and transformation induced anisotropy. Additionally, attention is focused on the interaction between phase transformations and plasticity in terms of the inheritance of plastic strain. The effect of such interaction is demonstrated by elementary numerical studies.     

Stress-resultant plasticity theories for composite laminated plates

Constitutive laws are presented for the inelastic analysis of laminated composite plates. The implications of using an elastoplastic theory, applied in a stress-resultant formulation, are discussed and investigated. Two different stress-resultant plasticity theories are proposed, both of which overlook the matrix and fiber inelastic behavior and describe the inelastic response of the laminate as a function of overall laminate properties. Results from numerical experiments with the proposed models are compared with results obtained using a micromechanical elastoplastic composite constitutive model.     

A 3-D constitutive model for shape memory alloys incorporating pseudoelasticity and detwinning of self-accommodated martensite

A 3-D constitutive model for polycrystalline shape memory alloys (SMAs), based on a modified phase transformation diagram, is presented. The model takes into account both direct conversion of austenite into detwinned martensite as well as the detwinning of self-accommodated martensite. This model is suitable for performing numerical simulations on SMA materials undergoing complex thermomechanical loading paths in stress–temperature space. The model is based on thermodynamic potentials and utilizes three internal variables to predict the phase transformation and detwinning of martensite in polycrystalline SMAs. Complementing the theoretical developments, experimental data are presented showing that the phase transformation temperatures for the self-accommodated martensite to austenite and detwinned martensite to austenite transformations are different. Determination of some of the SMA material parameters from such experimental data is also discussed. The paper concludes with several numerical examples of boundary value problems with complex thermomechanical loading paths which demonstrate the capabilities of the model.     

Mechanics of soft active materials with phase evolution

This paper studies the mechanics of soft active materials where the actuation is generated due to the formation of phases that are stress-free at the moment of their creation and therefore experience no deformation in the associated configuration. Phase formation is a continuous time-dependent process, which results in individual phases forming at different times and in different configurations of the material body, and thus it is coupled with mechanical deformation. Subsequent deformation of the material body results in individual phases experiencing different states of deformation and the overall material response results from the combined responses of the individual phases weighted by their respective volume fractions. Therefore, a great challenge in modeling the mechanics of soft active materials with evolving phases is to track the deformation and evolution of individual phases formed at different times and in different configurations. In this paper, a generalized one-dimensional model framework is presented to address this challenge. However, this model proves to be computationally inefficient. In response, an effective phase model is developed that tracks the combined deformation histories of new phases through a single, effective deformation. Both the general and effective phase models are evaluated with two fundamentally distinct phase evolution rules for three common mechanical problems: extension, stress relaxation, and creep. The first evolution rule represents a discrete transition from one phase to another while the second rule corresponds to a general transition from several phases into one phase. The effective phase model demonstrates excellent agreement with the generalized theory for all three mechanical problems considered under both types of evolution rules.