skip to main content


This content will become publicly available on August 1, 2024

Title: A gradient regularized model for shape memory alloys
In this work a transformation strain gradient enhancement is introduced into a phenomenological constitutive model for the pseudoelastic behavior of shape memory alloys. The constitutive model is able to capture several unique features of the constitutive response of these materials during the transformation between austenite and martensite during the pseudoelastic response. These features include the asymmetry in the initial transformation stresses in tension versus compression, the asymmetry in the transformation strains in tension and compression, and finally the asymmetry in the hardening behavior in tension and compression. In fact, experiments have shown that untrained NiTi exhibits hardening during its transformation in compression, but softening for tensile loading. It is this softening behavior that motivates the need for the introduction of the transformation strain gradient into the constitutive modeling. Transformation strain gradient effects are introduced via a phase variable that describes the extent of transformation. The free energy of the material then depends on gradients of the phase variable, which introduces a material length scale into the theory. The governing equation for the phase variable is developed from a microforce balance and continuum thermodynamics analysis. The model is implemented in the commercial finite element software Abaqus through user defined subroutines and several numerical simulations are performed to illustrate the model response and lack of numerical mesh-dependency of the results.  more » « less
Award ID(s):
1762389
NSF-PAR ID:
10479637
Author(s) / Creator(s):
;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Mechanics of Materials
Volume:
183
Issue:
C
ISSN:
0167-6636
Page Range / eLocation ID:
104689
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This paper describes the main findings from an experimental investigation into local and overall strength and fracture behavior of a microstructurally flexible, quadruplex, high entropy alloy (HEA), Fe42Mn28Co10Cr15Si5 (in at%). The alloy consists of metastable face-centered cubic austenite (g), stable hexagonal epsilon martensite (ε), stable body-centered cubic ferrite (a), and stable tetragonal sigma (σ) phases. The overall behavior of the alloy in compression features a great deal of plasticity and strain hardening before fracture. While the contents of diffusion created a and σ phases remain constant during deformation, the fraction of ε increases at the expanse of g due to the diffusionless strain induced γ→ε phase transformation. High-throughput nanoindentation mapping is used to assess the mechanical hardness of individual phases contributing to the plasticity and hardening of the alloy. Increasing the fraction of the dislocated ε phase during deformation due to the transformation is found to act as a secondary source of hardening because g and ε exhibit similar hardness at a given strain level. While these two phases exhibit moderate hardening during plasticity, significant softening is observed in σ owing to the phase fragmentation. While the phase transformation mechanism facilitates accommodation of the plasticity, the primary source of strain hardening in the alloy is the refinement of the structure during the transformation inducing a dynamic Hall-Petch-type barrier effect. Results pertaining to the evolution of microstructure and local behavior of the alloy under compression are presented and discussed clarifying the origins of strain hardening. While good under compression, the alloy poorly behaves under tension. Fracture surfaces after tension feature brittle micromechanisms of fracture. Such behavior is attributed to the presence of the brittle σ phase. 
    more » « less
  2. A novel experimental setup is presented that allows for precise control of thermal and mechanical loads, and simultaneous monitoring of the temperature and the full-field deformation of small SMA structures that undergo phase transformations. The facility is used to conduct two experi- ments in which NiTi tubes are taken through a temperature cycle under constant load that leads to phase transformations in the form of helical localization bands that propagate along the specimen. The latent heat of transformation causes a complex interaction with the prescribed load and thermal environment. By changing the rate of the airflow through the environmental chamber it is revealed that the velocities of the transformation fronts depend on the rate at which heat is removed/added by the controlled environment. The experiments are simulated using a new fully coupled thermomechanical extension of the constitutive framework developed by this research group. Key features of the framework include the modeling of the reversible A⇄M transformation through a single surface in the deviatoric stress-temperature space that obeys kinematic hardening; with the transformation strain and entropy as the internal variables gov- erned by an associative flow rule; and the inhomogeneous deformation exhibited in tension being modeled as softening. The tube is analyzed in a finite element coupled static displacement transient temperature analysis, and taken through the cool/heat cycle of the experiment. The temperature-strain response is accurately reproduced with the two transformations initiating at essentially the same temperatures as in the experiment and propagating in similar localized banded manners at similar speeds. Reproduction of the complex behavior observed in the ex- periments requires the calibration of the constitutive model, its discretization, and the modeling of the structure and its boundary conditions to work together to near perfection. The simulation also demonstrated that the heat exchange between the structure and the environment, in the present analysis governed by only by convection, requires further enhancement. 
    more » « less
  3. Summary

    Geomaterials such as sand and clay are highly heterogeneous multiphase materials. Nonlocality (or a characteristic length scale) in modeling geomaterials based on the continuum theory can be associated with several factors, for instance, the physical interactions of material points within finite distance, the homogenization or smoothing process of material heterogeneity, and the particle or problem size‐dependent mechanical behavior (eg, the thickness of shear bands) of geomaterials. In this article, we formulate a nonlocal elastoplastic constitutive model for geomaterials by adapting a local elastoplastic model for geomaterials at a constant suction through the constitutive correspondence principle of the state‐based peridynamics theory. We numerically implement this nonlocal constitutive model via the classical return‐mapping algorithm of computational plasticity. We first conduct a one‐dimensional compression test of a soil sample at a constant suction through the numerical model with three different values of the nonlocal variable (horizon)δ. We then present a strain localization analysis of a soil sample under the constant suction and plane strain conditions with different nonlocal variables. The numerical results show that the proposed nonlocal model can be used to simulate the inception and propagation of shear banding as well as to capture the thickness of shear bands in geomaterials at a constant suction.

     
    more » « less
  4. Scale-free phase-field approach (PFA) at large strains and corresponding finite element method (FEM) simulations for multivariant martensitic phase transformation (PT) from cubic Si I to tetragonal Si II in a polycrystalline aggregate are presented. Important features of the model are large and very anisotropic transformation strain tensor εt = {0.1753; 0.1753; −0.447} and stress-tensor dependent athermal dissipative threshold for PT, which produce essential challenges for computations. 3D polycrystals with 55 and 910 stochastically oriented grains are subjected to uniaxial strain- and stress-controlled loadings under periodic boundary conditions and zero averaged lateral strains. Coupled evolution of discrete martensitic microstructure, volume fractions of martensitic variants and Si II, stress and transformation strain tensors, and texture are presented and analyzed. Macroscopic variables effectively representing multivariant transformational behavior are introduced. Macroscopic stress-strain and transformational behavior for 55 and 910 grains are close (less than 10% difference). This allows the determination of macroscopic constitutive equations by treating aggregate with a small number of grains. Large transformation strains and grain boundaries lead to huge internal stresses of tens GPa, which affect microstructure evolution and macroscopic behavior. In contrast to a single crystal, the local mechanical instabilities due to PT and negative local tangent modulus are stabilized at the macroscale by arresting/slowing the growth of Si II regions by the grain boundaries and generating the internal back stresses. This leads to increasing stress during PT. The developed methodology can be used for studying similar PTs with large transformation strains and for further development by including plastic strain and strain-induced PTs. 
    more » « less
  5. Abstract

    A continuum damage model was developed to describe the finite tensile deformation of tough double‐network (DN) hydrogels synthesized by polymerization of a water‐soluble monomer inside a highly crosslinked rigid polyelectrolyte network. Damage evolution in DN hydrogels was characterized by performing loading‐unloading tensile tests and oscillatory shear rheometry on DN hydrogels synthesized from 3‐sulfopropyl acrylate potassium salt (SAPS) and acrylamide (AAm). The model can explain all the mechanical features of finite tensile deformation of DN hydrogels, including idealized Mullins effect and permanent set observed after unloading, qualitatively and quantitatively. The constitutive equation can describe the finite elasto‐plastic tensile behavior of DN hydrogels without resorting to a yield function. It was showed that tensile mechanics of DN hydrogels in the model is controlled by two material parameters which are related to the elastic moduli of first and second networks. In effect, the ratio of these two parameters is a dimensionless number that controls the behavior of material. The model can capture the stable branch of material response during neck propagation where engineering stress becomes constant. Consistent with experimental data, by increasing the elastic modulus of the second network the finite tensile behavior of the DN hydrogel changes from necking to strain hardening.

     
    more » « less