# American Institute of Mathematical Sciences

November  2021, 14(11): 3925-3952. doi: 10.3934/dcdss.2020459

## An experimentally-fitted thermodynamical constitutive model for polycrystalline shape memory alloys

 1 Faculty of Mathematics and Physics, Charles University in Prague, Sokolovská 83, CZ-186 75 Praha 8, Czech Republic 2 Institute of Thermomechanics, Czech Academy of Sciences, Dolejškova 5, CZ-182 00 Praha 8, Czech Republic 3 Institute of Physics, Czech Acad. Sci., Prague, Czech Republic, Na Slovance, CZ-18121 Praha 8, Czech Republic 4 Faculty of Nuclear Sciences and Physical Engineering, Czech Tech. Univ., Trojanova 13, CZ-120 00 Praha 2, Czech Republic

* Corresponding author: Barbora Benešová

Received  March 2020 Revised  June 2020 Published  November 2021 Early access  November 2020

Fund Project: This research has been partially supported from the grants 18-03834S (especially regarding the focus on shape memory alloys, experiments and modeling), 19-04956S (especially regarding the focus on the dynamic and nonlinear behavior) of Czech Science Foundation and LTAUSA18199 (especially regarding the focus on numerics) of MSMT CR, and also from the institutional support RVO: 61388998. Moreover, T.R. acknowledges a stay at Caltech in 2004 and discussions with Kaushik Bhattacharya about the isothermal variant of this sort of models

A phenomenological model for polycrystalline NiTi shape-memory alloys with a refined dissipation function is here enhanced by a thermomechanical coupling and rigorously analyzed as far as existence of weak solutions and numerical stability and convergence of the numerical approximation performed by a staggered time discretization. Moreover, the model is verified on one-dimensional computational simulations compared with real laboratory experiments on a NiTi wire.

Citation: Barbora Benešová, Miroslav Frost, Lukáš Kadeřávek, Tomáš Roubíček, Petr Sedlák. An experimentally-fitted thermodynamical constitutive model for polycrystalline shape memory alloys. Discrete & Continuous Dynamical Systems - S, 2021, 14 (11) : 3925-3952. doi: 10.3934/dcdss.2020459
##### References:

show all references

##### References:
Dependence of stress on strain for one isothermal simulation and adiabatic simulations at three different strain rates; see text for details
Evolution of temperature (left vertical axis) and volume fraction of martensite (right vertical axis) with relative progress of the stretching (horizontal axis: 0 = start of stretching, 0.5 = maximum tension, 1 = complete unstretching) at three different strain rates
Evolution of cumulative sum of three respective contributions to the total heat – latent heat, energy dissipated in rate-independent (RI) processes and in rate-dependent (RD) processes, respectively – integrated over the whole wire at three different strain rates. Relative progress of stretching on the horizontal axis (0 denotes start of stretching, 0.5 maximum tension, 1 complete unstretching). The RI contribution for $v_2$ basically coincides with the RI contribution of $v_1$
Comparison of in-situ experimental (left) and computational (right) evolution of spatial distribution of strain during loading at the total-strain rate $v_{\rm s} = 10^{-3}\, {\rm s}^{-1}$ of a thin NiTi wire
but under loading at the higher speed, $v_{\rm f} = 10^{-1}\, {\rm s}^{-1}$">Figure 5.  The same as in Fig. 4 but under loading at the higher speed, $v_{\rm f} = 10^{-1}\, {\rm s}^{-1}$
Main nomenclature for variables and for the data of the model
 $\varOmega$ a domain in $\mathbb{R}^d$ (a specimen) $K$ bulk modulus $\varGamma$ the boundary of $\varOmega$ $G=G(\xi)$ (a part of) shear modulus $u$ displacement $\mathcal{G}=\mathcal{G}(\xi, e)$ shear strain energy $\varepsilon(u)$ total (small) strain $=\frac12(\nabla u)^\top\!{+}\frac12\nabla u$ $c_{\rm v}$ heat capacity, $e$, $\pi$ elastic and inelastic strains $\mathbb{K}$ heat-transfer tensor $\xi$ volume fraction of martensite $\theta_{_{\rm D}}$ prescribed outer temperature $\theta$ temperature $u_{_{\rm D}}$ prescribed boundary displacement $s$ entropy $s_\text{AM}$ entropy of A/M-transformation $w$ heat part of the internal energy $\theta_\mathrm{tr}$ the transformation temperature $\eta$ a regularization parameter for approximation $\tau$ a time step for time discretisation
 $\varOmega$ a domain in $\mathbb{R}^d$ (a specimen) $K$ bulk modulus $\varGamma$ the boundary of $\varOmega$ $G=G(\xi)$ (a part of) shear modulus $u$ displacement $\mathcal{G}=\mathcal{G}(\xi, e)$ shear strain energy $\varepsilon(u)$ total (small) strain $=\frac12(\nabla u)^\top\!{+}\frac12\nabla u$ $c_{\rm v}$ heat capacity, $e$, $\pi$ elastic and inelastic strains $\mathbb{K}$ heat-transfer tensor $\xi$ volume fraction of martensite $\theta_{_{\rm D}}$ prescribed outer temperature $\theta$ temperature $u_{_{\rm D}}$ prescribed boundary displacement $s$ entropy $s_\text{AM}$ entropy of A/M-transformation $w$ heat part of the internal energy $\theta_\mathrm{tr}$ the transformation temperature $\eta$ a regularization parameter for approximation $\tau$ a time step for time discretisation
Material parameters used in simulations
 Parameter Value Unit Parameter Value Unit $l$ $25$ [mm] $\rho$ $0.1$ [mm] $f^{\rm tens}$ $0.08$ [1] $K$ $148$ [GPa] $G_{\rm A}$ $25$ [GPa] $G_{\rm M}$ $12$ [GPa] $c_{\rm v}$ $5$ [MPa/$^\circ$C] $E^{\rm nl}$ $80$ [MPa] $E^{\rm int}$ $10$ [MPa] $s_{\rm AM}^{}$ $360$ [kPa/$^\circ$C] $\beta$ $5$ [1] $\mu$ $100$ [kPa$\,$s] $\sigma^{\rm reo}_{\rm tr}$ $90$ [MPa] $\Sigma^{\rm reo}$ $-10$ [kPa/$^\circ$C] $A_{\rm s}$ $-15$ [$^\circ$C] $A_{\rm f}$ $-5$ [$^\circ$C] $M_{\rm s}$ $-20$ [$^\circ$C] $M_{\rm f}$ $-25$ [$^\circ$C] $T_0$ $-20$ [$^\circ$C] $\mathbb{K}_{11}$ $90$ [W/(m$^\circ$C)]
 Parameter Value Unit Parameter Value Unit $l$ $25$ [mm] $\rho$ $0.1$ [mm] $f^{\rm tens}$ $0.08$ [1] $K$ $148$ [GPa] $G_{\rm A}$ $25$ [GPa] $G_{\rm M}$ $12$ [GPa] $c_{\rm v}$ $5$ [MPa/$^\circ$C] $E^{\rm nl}$ $80$ [MPa] $E^{\rm int}$ $10$ [MPa] $s_{\rm AM}^{}$ $360$ [kPa/$^\circ$C] $\beta$ $5$ [1] $\mu$ $100$ [kPa$\,$s] $\sigma^{\rm reo}_{\rm tr}$ $90$ [MPa] $\Sigma^{\rm reo}$ $-10$ [kPa/$^\circ$C] $A_{\rm s}$ $-15$ [$^\circ$C] $A_{\rm f}$ $-5$ [$^\circ$C] $M_{\rm s}$ $-20$ [$^\circ$C] $M_{\rm f}$ $-25$ [$^\circ$C] $T_0$ $-20$ [$^\circ$C] $\mathbb{K}_{11}$ $90$ [W/(m$^\circ$C)]
 [1] Tomáš Roubíček. Modelling of thermodynamics of martensitic transformation in shape-memory alloys. Conference Publications, 2007, 2007 (Special) : 892-902. doi: 10.3934/proc.2007.2007.892 [2] Matthieu Hillairet, Alexei Lozinski, Marcela Szopos. On discretization in time in simulations of particulate flows. Discrete & Continuous Dynamical Systems - B, 2011, 15 (4) : 935-956. doi: 10.3934/dcdsb.2011.15.935 [3] Raghda A. M. Attia, Dumitru Baleanu, Dianchen Lu, Mostafa M. A. Khater, El-Sayed Ahmed. Computational and numerical simulations for the deoxyribonucleic acid (DNA) model. Discrete & Continuous Dynamical Systems - S, 2021, 14 (10) : 3459-3478. doi: 10.3934/dcdss.2021018 [4] Yongqiang Fu, Xiaoju Zhang. Global existence and asymptotic behavior of weak solutions for time-space fractional Kirchhoff-type diffusion equations. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021091 [5] Tong Li, Kun Zhao. Global existence and long-time behavior of entropy weak solutions to a quasilinear hyperbolic blood flow model. Networks & Heterogeneous Media, 2011, 6 (4) : 625-646. doi: 10.3934/nhm.2011.6.625 [6] Chungang Shi, Wei Wang, Dafeng Chen. Weak time discretization for slow-fast stochastic reaction-diffusion equations. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021019 [7] Carlos J. Garcia-Cervera, Xiao-Ping Wang. Spin-polarized transport: Existence of weak solutions. Discrete & Continuous Dynamical Systems - B, 2007, 7 (1) : 87-100. doi: 10.3934/dcdsb.2007.7.87 [8] Changchun Liu, Jingxue Yin, Juan Zhou. Existence of weak solutions for a generalized thin film equation. Communications on Pure & Applied Analysis, 2007, 6 (2) : 465-480. doi: 10.3934/cpaa.2007.6.465 [9] Marco Campo, José R. Fernández, Maria Grazia Naso. A dynamic problem involving a coupled suspension bridge system: Numerical analysis and computational experiments. Evolution Equations & Control Theory, 2019, 8 (3) : 489-502. doi: 10.3934/eect.2019024 [10] Peiying Chen. Existence and uniqueness of weak solutions for a class of nonlinear parabolic equations. Electronic Research Announcements, 2017, 24: 38-52. doi: 10.3934/era.2017.24.005 [11] Alain Bensoussan, Miroslav Bulíček, Jens Frehse. Existence and compactness for weak solutions to Bellman systems with critical growth. Discrete & Continuous Dynamical Systems - B, 2012, 17 (6) : 1729-1750. doi: 10.3934/dcdsb.2012.17.1729 [12] Md. Rabiul Haque, Takayoshi Ogawa, Ryuichi Sato. Existence of weak solutions to a convection–diffusion equation in a uniformly local lebesgue space. Communications on Pure & Applied Analysis, 2020, 19 (2) : 677-697. doi: 10.3934/cpaa.2020031 [13] Shijin Ding, Boling Guo, Junyu Lin, Ming Zeng. Global existence of weak solutions for Landau-Lifshitz-Maxwell equations. Discrete & Continuous Dynamical Systems, 2007, 17 (4) : 867-890. doi: 10.3934/dcds.2007.17.867 [14] Helmut Abels, Johannes Kampmann. Existence of weak solutions for a sharp interface model for phase separation on biological membranes. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 331-351. doi: 10.3934/dcdss.2020325 [15] M.J. Lopez-Herrero. The existence of weak solutions for a general class of mixed boundary value problems. Conference Publications, 2011, 2011 (Special) : 1015-1024. doi: 10.3934/proc.2011.2011.1015 [16] Shihui Zhu. Existence and uniqueness of global weak solutions of the Camassa-Holm equation with a forcing. Discrete & Continuous Dynamical Systems, 2016, 36 (9) : 5201-5221. doi: 10.3934/dcds.2016026 [17] Lukáš Poul. Existence of weak solutions to the Navier-Stokes-Fourier system on Lipschitz domains. Conference Publications, 2007, 2007 (Special) : 834-843. doi: 10.3934/proc.2007.2007.834 [18] Roman M. Taranets, Jeffrey T. Wong. Existence of weak solutions for particle-laden flow with surface tension. Discrete & Continuous Dynamical Systems, 2018, 38 (10) : 4979-4996. doi: 10.3934/dcds.2018217 [19] Shaoyong Lai, Qichang Xie, Yunxi Guo, YongHong Wu. The existence of weak solutions for a generalized Camassa-Holm equation. Communications on Pure & Applied Analysis, 2011, 10 (1) : 45-57. doi: 10.3934/cpaa.2011.10.45 [20] Carlos J. García-Cervera, Xiao-Ping Wang. A note on 'Spin-polarized transport: Existence of weak solutions'. Discrete & Continuous Dynamical Systems - B, 2015, 20 (8) : 2761-2763. doi: 10.3934/dcdsb.2015.20.2761

2020 Impact Factor: 2.425