Computational modeling of size-dependent superelasticity of shape memory alloys

We propose a nonlocal continuum model to describe the size-dependent superelastic responses observed in recent experiments of shape memory alloys. The modeling approach extends a superelasticity formulation based on the martensitic volume fraction, and combines it with gradient plasticity theories. Size effects are incorporated through two internal length scales, an energetic length scale and a dissipative length scale, which correspond to the gradient terms in the free energy and the dissipation, respectively. We also propose a computational framework based on a variational formulation to solve the coupled governing equations resulting from the nonlocal superelastic model. Within this framework, a robust and scalable algorithm is implemented for large scale three-dimensional problems. A numerical study of the grain boundary constraint effect shows that the model is able to capture the size-dependent stress hysteresis and strain hardening during the loading and unloading cycles in polycrystalline SMAs.     

Finite deformation thermo-mechanical behavior of thermally induced shape memory polymers

Shape memory polymers (SMPs) are polymers that can demonstrate programmable shape memory effects. Typically, an SMP is pre-deformed from an initial shape to a deformed shape by applying a mechanical load at the temperature T_H>T_g. It will maintain this deformed shape after subsequently lowering the temperature to T_L<T_g and removing the externally mechanical load. The shape memory effect is activated by increasing the temperature to T_D>T_g, where the initial shape is recovered. In this paper, the finite deformation thermo-mechanical behaviors of amorphous SMPs are experimentally investigated. Based on the experimental observations and an understanding of the underlying physical mechanism of the shape memory behavior, a three-dimensional (3D) constitutive model is developed to describe the finite deformation thermo-mechanical response of SMPs. The model in this paper has been implemented into an ABAQUS user material subroutine (UMAT) for finite element analysis, and numerical simulations of the thermo-mechanical experiments verify the efficiency of the model. This model will serve as a modeling tool for the design of more complicated SMP-based structures and devices.     

Simulations of localized thermo-mechanical behavior in a NiTi shape memory alloy

Previous experiments have shown that stress-induced martensitic transformation in certain polycrystalline NiTi shape memory alloys can lead to strain localization and propagation phenomena when loaded in uniaxial tension. The number of nucleation events and kinetics of transformation fronts were found to be sensitive to the nature of the ambient media and imposed loading rate due to the release/absorption of latent heat and the material's inherent temperature sensitivity of the transformation stress. A special plasticity-based constitutive model used within a 3-D finite element framework has previously been shown to capture the isothermal, purely mechanical front features seen in experiments of thin uniaxial NiTi strips. This paper extends the approach to include the thermo-mechanical coupling of the material with its environment. The simulations successfully capture the nucleation and evolution of fronts and the corresponding temperature fields seen during the experiments.     

Coupled thermo-mechanical analysis of shape memory alloy circular bars in pure torsion

Pure torsion of shape memory alloy (SMA) bars with circular cross section is studied by considering the effect of temperature gradient in the cross sections as a result of latent heat generation and absorption during forward and reverse phase transformations. The local form of energy balance for SMAs by taking into account the heat flux effect is coupled to a closed-form solution of SMA bars subjected to pure torsion. The resulting coupled thermo-mechanical equations are solved for SMA bars with circular cross sections. Several numerical case studies are presented and the necessity of considering the coupled thermo-mechanical formulation is demonstrated by comparing the results of the proposed model with those obtained by assuming an isothermal process during loading–unloading. Pure torsion of SMA bars in various ambient conditions (free and forced convection of air, and forced convection of water flow) subjected to different loading–unloading rates are studied and it is shown that the isothermal solution is valid only for specific combinations of ambient conditions and loading rates.     

A macroscopic multi-mechanism based constitutive model for the thermo-mechanical cyclic degeneration of shape memory effect of NiTi shape memory alloy

A macroscopic based multi-mechanism constitutive model is constructed in the framework of irreversible thermodynamics to describe the degeneration of shape memory effect occurring in the thermo-mechanical cyclic deformation of NiTi shape memory alloys(SMAs). Three phases,austenite A, twinned martensite Mtand detwinned martensite M~d, as well as the phase transitions occurring between each pair of phases( A → M~t, M~t→ A, A → M~d,M~d→ A, and M~t→ M~d) are considered in the proposed model. Meanwhile, two kinds of inelastic deformation mechanisms, martensite transformation-induced plasticity and reorientation-induced plasticity, are used to explain the degeneration of shape memory effects of NiTi SMAs. The evolution equations of internal variables are proposed by attributing the degeneration of shape memory effect to the interaction between the three phases(A, M~t, and M~d) and plastic deformation. Finally, the capability of the proposed model is verified by comparing the predictions with the experimental results of NiTi SMAs. It is shown that the degeneration of shape memory effect and its dependence on the loading level can be reasonably described by the proposed model.     

Dynamic multi-axial behavior of shape memory alloy nanowires with coupled thermo-mechanical phase-field models

The objective of this paper is to provide new insight into the dynamic thermo-mechanical properties of shape memory alloy (SMA) nanowires subjected to multi-axial loadings. The phase-field model with Ginzburg–Landau energy, having appropriate strain based order parameter and strain gradient energy contributions, is used to study the martensitic transformations in the representative 2D square-to-rectangular phase transformations for FePd SMA nanowires. The microstructure and mechanical behavior of martensitic transformations in SMA nanostructures have been studied extensively in the literature for uniaxial loading, usually under isothermal assumptions. The developed model describes the martensitic transformations in SMAs based on the equations for momentum and energy with bi-directional coupling via strain, strain rate and temperature. These governing equations of the thermo-mechanical model are numerically solved simultaneously for different external loadings starting with the evolved twinned and austenitic phases. We observed a strong influence of multi-axial loading on dynamic thermo-mechanical properties of SMA nanowires. Notably, the multi-axial loadings are quite distinct as compared to the uniaxial loading case, and the particular axial stress level is reached at a lower strain. The SMA behaviors predicted by the model are in qualitative agreements with experimental and numerical results published in the literature. The new results reported here on the nanowire response to multi-axial loadings provide new physical insight into underlying phenomena and are important, for example, in developing better SMA-based MEMS and NEMS devices     

Modeling the dynamic behavior of shape memory alloys

The paper studies the single degree of freedom vibration of a rigid mass suspended by a thin-walled shape memory alloy tube under torsional loading. The behavior is analyzed for the cases of quasiplasticity (low temperatures) and pseudoelasticity (high temperatures) on the basis of an improved version of the Müller–Achenbach model. To illustrate the strong hysteresis-induced damping capacity and the non-linear vibration characteristics, both, free and forced vibrations are considered in the first part of the paper. This is done on the basis of an isothermal version of the model, while the second part of the paper focuses on the effect of non-constant temperature caused by the rate-dependent release and absorption of latent heats.     

Study on behaviors of functionally graded shape memory alloy cylinder

For better controllability in actuations,it is desirable to create Functionally Graded Shape Memory Alloys(FG-SMAs)in the actuation direction.It can be achieved by applying different heat treatment processes to create the gradient along the radius of a SMA cylinder.Analytical solutions are derived to predict the macroscopic behaviors of such a functionally graded SMA cylinder.The Tresca yield criterion and linear hardening are used to describe the different phase transformations with different gradient parameters.The numerical results for an example of the model exhibit different pseudo-elastic behaviors from the non-gradient case,as well as a variational hysteresis loop for the transformation,providing a mechanism for easy actuation control.When the gradient disappears,the model can degenerate to the non-gradient case.     

Constitutive modeling of shape memory polymer based self-healing syntactic foam

In a previous study, it was found that the shape memory functionality of a shape memory polymer based syntactic foam can be utilized to self-seal impact damage repeatedly, efficiently, and almost autonomously [Li G., John M., 2008. A self-healing smart syntactic foam under multiple impacts. Comp. Sci. Technol. 68(15–16), 3337–3343]. The purpose of this study is to develop a thermodynamics based constitutive model to predict the thermomechanical behavior of the smart foam. First, based on DMA tests and FTIR tests, the foam is perceived as a three-phase composite with interfacial transition zone (interphase) coated microballoons dispersed in the shape memory polymer (SMP) matrix; for simplicity, it is assumed to be an equivalent two-phase composite by dispersing elastic microballoons into an equivalent SMP matrix. Second, the equivalent SMP matrix is phenomenologically assumed to consist of an active (rubbery) phase and a frozen (glassy) phase following Liu et al. [Liu, Y., Gall, K., Dunn, M.L., Greenberg, A.R., Diani J., 2006. Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modeling. Int. J. Plasticity 22, 279–313]. The phase transition between these two phases is through the change of the volume fraction of each phase and it captures the thermomechanical behavior of the foam. The time rate effect is also considered by using rheological models. With some parameters determined by additional experimental testing, the prediction by this model is in good agreement with the 1D test result found in the literature. Parametric studies are also conducted using the constitutive model, which provide guidance for future design of this novel self-healing syntactic foam and a class of light-weight composite sandwich structures.     

Analysis of the rate-dependent coupled thermo-mechanical response of shape memory alloy bars and wires in tension

In this paper, the coupled thermo-mechanical response of shape memory alloy (SMA) bars and wires in tension is studied. By using the Gibbs free energy as the thermodynamic potential and choosing appropriate internal state variables, a three-dimensional phenomenological macroscopic constitutive model for polycrystalline SMAs is derived. Taking into account the effect of generated (absorbed) latent heat during the forward (inverse) martensitic phase transformation, the local form of the first law of thermodynamics is used to obtain the energy balance relation. The three-dimensional coupled relations for the energy balance in the presence of the internal heat flux and the constitutive equations are reduced to a one-dimensional problem. An explicit finite difference scheme is used to discretize the governing initial-boundary-value problem of bars and wires with circular cross-sections in tension. Considering several case studies for SMA wires and bars with different diameters, the effect of loading–unloading rate and different boundary conditions imposed by free and forced convections at the surface are studied. It is shown that the accuracy of assuming adiabatic or isothermal conditions in the tensile response of SMA bars strongly depends on the size and the ambient condition in addition to the rate dependency that has been known in the literature. The data of three experimental tests are used for validating the numerical results of the present formulation in predicting the stress–strain and temperature distribution for SMA bars and wires subjected to axial loading–unloading.     

面向火星探测的环氧形状记忆聚合物复合材料的性能研究

国家火星探测任务是建设航天强国进程中的重大标志性工程,是中国航天走向更远深空的里程碑工程.智能材料这种集材料、结构和功能于一体的先进材料将会对火星探测任务有所助力.形状记忆聚合物及其复合材料作为一种典型的智能材料,可在有效减轻载荷的同时实现自主变形,已经在地球同步轨道航天器的应用中崭露头角.因此有必要研究这种新型环氧基形状记忆聚合物复合材料应用于火星探测工程的可能性.首先,针对“天问一号”火星探测器的任务需求,设计了一个具有自释放功能的着陆平台国旗装置.其中的锁紧释放装置由碳纤维增强的形状记忆聚合物复合材料制成,分别从静态拉伸力学性能、动态热机械性能和形状记忆性能三个角度评估了空间辐照和长期存储对形状记忆聚合物复合材料的影响.其中,空间辐照包括γ射线和紫外射线,辐照剂量分别为5×10⁵ rad和23.6 kCal.长期存储分为低温-196℃、室温25℃和高温85℃存储30天,和低温-196℃存储457天两组实验.最后,从“祝融号”火星车所携带相机拍摄的照片可以看到五星红旗被成功释放,旗面平整、图案清晰.这说明所研究的环氧基形状记忆聚合物复合材料可成功应用于火星探测...     

显式方法精确模拟形状记忆聚合物热力学行为

通过构建一个热耦合的多轴可压缩应变能函数,得到应力-应变、应力-温度和应变-温度之间的函数关系,建立形状记忆聚合物的本构方程.本文引入三个基于对数应变的不变量使得模型（i）可以模拟可压缩情况;（ii）适用于单轴拉伸和等双轴拉伸至少两个基准实验;（iii）多轴有效.通过显式方法（i）给出自由能和熵的具体表达,证明模型热力学定律;（ii）给出应变-应力,温度-应力以及,温度-应变的形函数具体表达.多轴模型在特定的情况下可以自动退化到各自的单轴情况. 通过调节形函数的参数,最终得到的模型结果和实验结果能够精确匹配.新方法建立的本构模型得到的结果能更加准确地指导形状记忆聚合物的工程设计。     

显式方法精确模拟形状记忆聚合物热力学行为

通过构建一个热耦合的多轴可压缩应变能函数,得到应力-应变、应力-温度和应变-温度之间的函数关系,建立形状记忆聚合物的本构方程.本文引入三个基于对数应变的不变量使得模型（i）可以模拟可压缩情况;（ii）适用于单轴拉伸和等双轴拉伸至少两个基准实验;（iii）多轴有效.通过显式方法（i）给出自由能和熵的具体表达,证明模型热力学定律;（ii）给出应变-应力,温度-应力以及,温度-应变的形函数具体表达.多轴模型在特定的情况下可以自动退化到各自的单轴情况. 通过调节形函数的参数,最终得到的模型结果和实验结果能够精确匹配.新方法建立的本构模型得到的结果能更加准确地指导形状记忆聚合物的工程设计。     

Predicting thermal shape memory of crosslinked polymer networks from linear viscoelasticity

The viscoelastic behavior of an amorphous shape memory polymer network and its dependence on time and temperature were measured by dynamic mechanical analysis. The resulting thermo-mechanical behavior was modeled and implemented in a commercial finite element code. The ability of the resulting thermomechanical model to simulate and, eventually, predict the shape storage and shape recovery of the material was evaluated against experimental shape memory thermomechanical torsion data in a large deformation regimen. The simulations showed excellent agreement with experimental shape memory thermomechanical cycle data. This demonstrates the dependence of the shape recovery on time and temperature. The results suggest that accurate predictions of the shape recovery of any amorphous polymer networks under any thermomechanical conditions combination solely depends on considering the material viscoelasticity and its time–temperature dependence.     

Modeling thermomechanical processes in shape memory polymers under finite deformations

A model taking into account finite deformations is constructed for the behavior of a shape memory polymer which undergoes a transition from the highly elastic to the vitreous state and back during deformation and temperature change. The obtained relations are tested on problems which have experimental support.     

Hysteresis in magnetic shape memory composites: Modeling and simulation

Magnetic shape memory alloys are characterized by the coupling between the reorientation of structural variants and the rearrangement of magnetic domains. This permits to control the shape change via an external magnetic field, at least in single crystals. Composite materials with single-crystalline particles embedded in a softer matrix have been proposed as a way to overcome the blocking of the reorientation at grain boundaries.We investigate hysteresis phenomena for small NiMnGa single crystals embedded in a polymer matrix for slowly varying magnetic fields. The evolution of the microstructure is studied within the rate-independent variational framework proposed by Mielke and Theil (1999). The underlying variational model incorporates linearized elasticity, micromagnetism, stray field and a dissipation term proportional to the volume swept by the twin boundary. The time discretization is based on an incremental minimization of the sum of energy and dissipation. A backtracking approach is employed to approximately ensure the global minimality condition.We illustrate and discuss the influence of the particle geometry (volume fraction, shape, arrangement) and the polymer elastic parameters on the observed hysteresis and compare with recent experimental results.     

Modeling elastic-inelastic processes in shape memory alloys at finite deformations

An equation of state for a shape memory alloy is derived using a formalized approach to the construction of finite-deformation constitutive equations for complex media. The obtained equations were tested for coupled elastic-inelastic boundary-value problems of deformation of a sample of a shape-memory during forward and reverse martensitic transformations.     

Thermomechanical behavior of shape memory polymer beams reinforced by corrugated polymeric sections

With the ever-increasing applications of smart actuators, great attention is drawn to the SMP-based composites. In this paper the flexural behavior of corrugated SMP composite beams is studied. Employing a widely-accepted one-dimensional SMP model, based on the assumptions of Euler–Bernoulli beam theory, the governing equations are obtained. Further using a finite difference scheme the equations are solved. In this regard, different types of corrugated sections (rectangular, sinusoidal, …) with equal SMP content are studied and the mechanical properties of interest (load capacity, shape fixity, …) are compared. It is observed that reinforced single-cell patterns have completely different mechanical behaviors. For the sake of generality, the single-cell reinforced composite sections are studied in detail. Numerical results show an increase in load capacity of the structure. However, like any other reinforcing method, an inevitable small decrease in shape fixity is observed. In addition to the properties studied here, other desirable characteristics can be achieved by introducing a reinforcing cover. The results can be utilized in designing SMP-based actuators since the present work proposes a simple and efficient method for enhancing the load capacity of the SMP-based composites.