Due to the effects of pseudoelasticity and pseudoplasticity, shape memory alloys (SMAs) are very promising materials for the industrial usage. However, applications of SMA are still challenging due to the functional degradation during cyclic loading. The related effect of functional fatigue which occurs during pseudoelastic loading is modeled by subdividing the diffusionless solid/solid phase transformation into a reversible and an irreversible process which is an experimentally motivated ansatz for the complex material behavior, see e.g. [1]. Thus, we take into account a reversible and an irreversible volume fraction for the austenitic and several martensitic phases. To consider the material's polycrystalline structure and hence, differently oriented grains, we use an orientation distribution function which depends on three Euler angles and affects a high numerical efficiency as presented in [2].
Micromechanical modelThe transformation from austenite to martensite in SMAs is accompanied by a formation of dislocations which causes a stabilization of martensite and thus, a functional degradation during cyclic loading, see e.g. [1]. The direct modeling of dislocations comes along with a high numerical effort. For a better numerical stability and computational performance, we skip the consideration of dislocations and directly take into account a stabilized martensitic volume fraction, see also [3]. By subdividing the phase transformation into a reversible and an irreversible process, we consider a reversible and an irreversible part of the volume fractions λ and ρ. Therefore, the quantity ρ accounts for the stabilized martensite during cyclic loading until a maximum value ρ max . In addition, we introduce a set of Euler angles α = {ϕ, ϑ, ω} (see [2]) which describes the averaged orientation of the transforming grains and hence, takes into account the polycrystalline structure. The presented model is based on the principle of the minimum of the dissipation potential, see e.g. [4,5]. By minimization of a Lagrange function L that consists of the rate of the Helmholtz free energy Ψ and a dissipation function D,(1) the variational method directly results in evolution equations for the used internal variables. The Helmholtz free energy is formulated with use of a Reuß energy bound. It depends on the total strain ε and the three internal variables λ, ρ and α, and is given byThe formulation(·) denotes the effective quantities which are a result of the chosen energy bound. The effective reversible and irreversible transformation strainη andῡ, the effective material's stiffness matrixĒ and the effective caloric energyc are calculated viāThe rotation tensor Q in (3) depends on the set of Euler angles and is used to rotate the effective transformation strains in the orientation of the transforming grains. The formulation of the dissipation function D is subdivided into a coupled part for the reversible and irreversible transformation and a part related to the orientation distribution -more precisely the evolution of α. In th...