A soil-bentonite cutoff wall is a type of subsurface vertical barrier constructed by back-filling a trench with a mixture of soil, bentonite, and water. Although soil-bentonite cutoff walls are common, their mechanical behavior is not well understood. Current design procedures do not consider the final stress state of the consolidated soil-bentonite backfill or deformations in adjacent ground. The final stress state in the completed wall is important because it influences the hydraulic conductivity of the cutoff (Barrier 1995), the cutoff's susceptibility to hydraulic fracture, and the magnitude of deformations adjacent to the cutoff wall. Deformations adjacent to the cutoff wall can be significant and can cause damage to adjacent structures. The objectives of this research are to 1) add to the current body of knowledge of the properties of soil-bentonite mixtures, 2) evaluate constitutive models and select a model to represent soil-bentonite, 3) model a soil-bentonite cutoff wall using finite elements, and 4) investigate the influence of several factors on the deformations in adjacent ground. These objectives were met by first summarizing information from the literature on soil-bentonite properties and then performing a laboratory testing program on different soil-bentonite mixtures. Five constitutive models were evaluated to determine how well they match the data from the laboratory testing program. A model referred to as the RS model was chosen to best represent soil-bentonite, and provided a good match of the soil-bentonite behavior. The RS model, which is a special case of a more complicated existing model, is a non-associative Modified Cam Clay type model that has parameters to change the yield surface and plastic potential surface into ellipses of varying shapes. The RS model was implemented into the finite element program SAGE. A finite element model was developed using SAGE to simulate all stages of construction of a soil-bentonite cutoff wall including excavation of a trench under bentonite-water slurry, replacement of the bentonite-water slurry with soil-bentonite backfill, and consolidation of the soil-bentonite backfill. The model was calibrated with a well-documented case history, and predicted deformations in adjacent ground were close to measured deformations. Evaluation of the model indicates that there is good confidence in the prediction of deformations in adjacent ground, but there is lower confidence in the predicted stresses in the consolidated soil-bentonite and settlement of the backfill in the trench. A parametric study was then performed using the finite element model assuming sand sites of varying density and OCR. Deformations in adjacent ground were calculated for various soil conditions, soil-bentonite properties, and trench configurations. A correlation was found between maximum calculated settlement in adjacent ground and factor of safety against trench / Ph. D.

The coastal population is growing, putting extra stress on coastal sediments and protection features, such as beach dunes. Moreover, global warming will increase the frequency of storms, and coastal dunes and other defense infrastructure will be subjected to increased erosion and scouring, endangering the people they are meant to protect. Understanding soil dynamics and fluid interaction is crucial to predict the effects of sand erosion. In particular, the study of wind erosion of sands in coastal dunes is essential due to the protective role these earthen structures have during storm events. One of the challenges about predicting wind erosion in coastal dunes is its extended spatial scale and the associated economic and logistics costs of sampling and characterizing the sediments. Because of this, in-situ testing for sediment characterization is essential. In particular, the usage of free-fall penetrometers (FFP) is appealing due to their portability and robustness. The sediment properties obtained with this type of testing can later be used to assess wind erosion susceptibility by determining, for example, the wind velocity to initiate the erosion process. FFP testing involves dropping an instrumented probe that impacts the soil and measures the kinematics or kinetics during the penetration process. For example, deceleration measurements are used to compute an equivalent quasi-static failure, which is not in line with the dynamic process characteristic of FFP testing. This preassumed failure mechanism is used to back-calculate the sand's geomechanical properties. However, soil behavior is highly complex under rapid loading, and incorporating this behavior into FFP sediment characterization models is challenging. Advanced numerical modeling can improve the understanding of the physics behind FFP testing. This thesis presents various advancements in numerical modeling and erosion models to bridge FFP in-situ testing with predicting the initiation of wind erosion of sands. First, improvements oriented to the Material Point Method (MPM) for modeling in-situ FFP testing are proposed. The numerical results show that the simulation of FFP deployment in sands is affected by strain localization and highlight the importance of considering constitutive models sensitive to different loading rates. Because of the importance of rate effects in soil behavior, the second aspect of this thesis proposes a novel consistency framework. Two constitutive models are adapted to study strain-rate sensitive non-cohesive materials: i) a strain-softening Mohr-Coulomb, and ii) a NorSand model. In addition to increased strength, the proposed framework captures increased dilatation, an early peak deviatoric stress, and relaxation. Finally, a novel sand erosion model is derived using a continuum approximation and limit equilibrium analysis. The erosion law considers geotechnical parameters, the effects of slope, and moisture suction, in a combined manner. The proposed model is theoretically consistent with existing expressions in the literature. It covers a wide range of environmental and geometrical conditions and helps to reconcile the results from FFP testing with the prediction of the initiation of wind erosion. The model was validated in a wind tunnel and is demonstrated to be a viable alternative for predicting sand erosion initiation. This thesis opens up new research prospects, such as improving the soil characterization models or the direct prediction of sand erosion using rapid, reliable, and efficient in-situ testing methods. / Doctor of Philosophy / With global warming and climate change, it is expected that the frequency and intensity of storms will increase. This increment will put extra stress on coastal sediments such as beach sand and coastal dunes, making them prone to erosion. Coastal dunes lose their ability to withstand storms as they erode, potentially making coastal flooding more frequent. In light of this, all stakeholders involved in the protection against coastal disasters must have the tools to predict, prepare for, and mitigate for situations like the ones stated above. An essential aspect of the prediction component is dependent on a successful sediment characterization, for example, determining how much wind the sand can withstand before it erodes. Free-fall penetrometers (FFP) are devices designed to conduct the characterization mentioned above. However, the procedures used to perform this characterization are mainly based on empirical or semi-empirical expressions. Computer models, capable of simulating the physics behind FFP testing, can bring more insight into the process of interaction between FFP devices, sands, and water and can be the basis to improve the characterization methods. The latter results can be utilized for instance to predict wind erosion, including several properties of the sand, such as its mineralogy and shape. This study contributes to developing the computer simulations of FFP deployment and the wind erosion prediction models. Eventually, these developments can help engineers and coastal managers to anticipate and prepare for more frequent coastal hazards.

Constitutive modeling is a method that is useful in providing precise predictions of material response in components subjected to a variety of operating conditions. A process for optimizing the material constants of the Miller constitutive model for uniaxial modeling was developed and implemented in an automated optimization routine. Generally, up to twenty experiments simulating a range of conditions are needed to identify the material parameters for the model. The research sought to minimize the amount of empirical data that is necessary for optimization, aiming to reduce the costs and time necessary to carry out this procedure for more expensive classes of materials. The ultimate goal was to develop a robust method for determining the material constants of a viscoplastic model using a minimum amount of experimental data. An automated optimization routine was implemented into a program, referred to as uSHARP, developed as part of the research to determine constitutive model parameters. Central to the method was the use of more complex stress, strain, and temperature histories than are traditionally used, allowing for the effects of all material parameters to be captured using as few tests as possible. By carrying out successive finite element simulations and comparing the results to simulated experimental test data, the material constants were determined from 75% fewer experiments. By reducing monetary costs and time required, this procedure will allow for a more widespread application of advanced constitutive models in industry, allowing for better life prediction modeling of critical components in high temperature applications.

Rubber-like materials are submitted to cyclic loading conditions in various applications. Fillers are always incorporated within rubber compounds. They improve the mechanical properties but induce a significant stress-softening under cyclic loadings. The physical source of the softening is not yet established and its modeling remains a challenge. For a better understanding of the softening, filled rubbers are submitted to cyclic loadings. In order to quantify the effects of the loading intensity and the number of cycles, original methods are proposed to characterize the softening. To study the influence of the material microstructure on the softening, compounds with various compositions are considered.Non proportional tensile tests including uniaxial and biaxial loading paths are applied in order to highlight the softening induced anisotropy. Such unconventional experimental data are used to provide a general criterion for the softening activation. A constitutive modeling grounded on a thorough analysis of experimental data is proposed. The model is based on a directional approach. The Mullins softening is accounted for by the strain amplification concept and is activated by a directional criterion. The model ability to predict non proportional softened material responses is demonstrated

L’étude présentée porte sur le développement d’un modèle numérique destiné à évaluer les performances balistiques d’une structure multicouche polymère/métal frittée par procédé SPS. Les matériaux sont un alliage d’aluminium 7020 et un polyimide thermoplastique amorphe qui sont ensuite assemblés avec une résine epoxy. Le comportement mécanique de ces trois matériaux a été étudié sur de larges gammes de vitesses de déformations (de 0.0001 /s à 50.000 /s) et de températures (de -70°C à 500°C) correspondant aux conditions extrêmes rencontrées lors d’impacts à hautes vitesses. Afin d’améliorer la précision des résultats, des approches analytiques ont été développées autant pour la modélisation du métal que pour celle les polymères. Après la calibration des modèles, ces derniers ont été implémenté dans ABAQUS®/Explicit (éléments finis) via des subroutines VUMAT en code FORTRAN. Des essais d’impacts de billes à hautes vitesses ont été réalisés sur des cibles monocouches pour valider les modèles numériques. De nombreuses configurations de composites multicouches ont ensuite été étudiées numériquement et leurs performances balistiques ont été comparées. / The present study deals with the development of a numerical model to evaluate the ballistic performance of a polymer/metal multilayer structure sintered by SPS. The materials are an aluminum alloy 7020 and an amorphous thermoplastic polyimide which are then assembled using an epoxy resin. The mechanical behavior of these three materials has been studied over wide ranges of strain rates (from 0.0001 / s to 50,000 / s) and temperatures (from -70 °C to 500 °C) corresponding to the extreme conditions encountered during impacts at high velocities. In order to improve the accuracy of the results, analytical approaches have been developed both for the modeling of the metal and for the polymers. After the calibration of the models, these models were implemented in ABAQUS® / Explicit (finite elements) via VUMAT subroutines in FORTRAN code. Ball impact tests at high velocities were performed on monolayer targets to validate numerical models. Numerous configurations of multilayer composites were then studied numerically and their ballistic performances have been compared.

Ce travail de thèse porte sur une catégorie de matériaux actifs dénommés Elastomères Magnéto-Rhéologiques (EMR). Ces derniers sont composés de particules micrométriques et magnétisables imprégnées dans une matrice élastomère isolante. Il est possible de modifier les propriétés mécaniques de tels matériaux en les soumettant à un champ magnétique externe. Avec pour objectif d’aboutir à une caractérisation couplée (magnéto-mécanique) du comportement des EMRs en grandes déformations et en présence de champs magnétiques élevés, ce travail propose une approche à la fois expérimentale, théorique et numérique.La première partie de ce travail s’intéresse à des aspects expérimentaux où l’influence de la microstructure (isotrope et transverse isotrope) et l’influence de la fraction volumique de particules sont étudiées. Un échantillon dédié est développé afin d’obtenir simultanément des champs mécaniques et magnétiques les plus homogènes possibles dans celui-ci lors d’une caractérisation couplée. La question de l’adhésion interfaciale entre les particules de fer doux et la matrice en silicone est également traitée et il est montré qu’un traitement chimique des particules est nécessaire afin d’éviter toute décohésion avec la matrice lorsque le matériau est soumis à un champ magnétique externe. Avant d’analyser les données obtenues, le dispositif expérimental permettant d’obtenir de manière simultanée une mesure du champ de déformation en trois dimensions et une mesure des champs magnétiques internes, est décrit. Malgré l’ensemble des difficultés expérimentales en grande partie dûes à des phénomènes d’instabilité qui sont omniprésents chez les EMRs, de nombreuses données sont collectées et serviront à la calibration des lois de comportement.La seconde partie de cette thèse couvre la modélisation couplée magnéto-mécanique des EMRs en s’appuyant sur le cadre théorique général des solides magnéto-élastiques proposé par Kankanala, Triantafyllidis et Danas (2004, 2012, 2014). En particulier, la méthode énergétique (qui s’appuie sur l’utilisation d’une fonction d’énergie libre) est préférée et des formulations variationnelles équivalentes (qui diffèrent entre elles simplement par le choix de la variable magnétique indépendante utilisée pour décrire le problème : B, H ou M) sont proposées et implémentées dans des codes numériques 3D s’appuyant sur la méthode des éléments finis. Ces outils numériques sont combinés à la méthode de minimisation des moindres carrés afin d’obtenir l’ensemble des paramètres matériaux du modèle de comportement des EMRs. L’utilisation de simulations numériques est nécessaire car une approche purement analytique ne permettrait pas de modéliser « l’effet de forme » observé expérimentalement. En effet, il est primordial de modéliser ce dernier car dans le cas contraire les paramètres identifiés dépendraient de la forme de l’échantillon expérimental et ne décriraient pas uniquement le matériau.La troisième partie de cette étude décrit en détail l’implémentation numérique des différentes formulations variationnelles proposées précédemment. Dans chacun des cas, il est prouvé que l’utilisation d’éléments isoparamétriques est bien adaptée. De nombreuses difficultés numériques ont été observées dans le cas des formulations variationnelles utilisant le champ de déplacement et le potentiel vecteur magnétique comme variables indépendantes. L’ensemble de ces difficultés (comme par exemple la minimisation de l’énergie potentielle sous la contrainte imposée par la jauge de Coulomb) est surmonté dans ce travail. Avant de décrire les différents problèmes tests utilisés pour s’assurer de la validité et de la précision des codes numériques, les différentes étapes nécessaires à la simulation d’un problème aux limites sont expliquées. Plus précisément, les questions liées aux spécificités des conditions aux limites à appliquer sur le potentiel vecteur magnétique ou encore aux conditions de symétries, sont traitées. / In this thesis, we study a class of “active materials” called Magnetorheological elastomers (MRE) which are ferromagnetic impregnated rubbers whose mechanical properties are altered by the application of external magnetic fields. With the purpose of characterizing the behavior of MREs up to large strains and high magnetic fields, this work brings a completely novel experimental, theoretical and numerical approach.The first part of this study focuses on an experimental investigation of MRE where multiple microstructures (isotropic and transversely isotropic materials) and multiple particles’ volume fraction are tested. A special sample geometry is designed in order to increase the uniformity of internal magnetic and mechanical fields measured during coupled-field experiments. The interfacial adhesion between the iron fillers and the silicone matrix is investigated and we show that when specimens are subjected to external magnetic fields, a silane primer treatment of the particles is needed to prevent debonding at the interface particle/matrix. Then, we present the magneto-mechanical testing setup that allows simultaneous 3D mechanical and magnetic measurements before discussing the results. Even if is found that instabilities are ubiquitous in MREs, lots of useful data are collected and will be used to compute the parameters proposed in the material model.The second part of the thesis is dedicated to the modeling of isotropic MREs. The continuum description proposed by Kankanala, Triantafyllidis and Danas (2004, 2012, 2014) to derive constitutive laws that account for finite strains is used and, in particular, the energetic approach (that requires an energy density function) is chosen. Multiple equivalent variational formulation alternatives (based on different choices of the independent magnetic variable used in the energy function: B, H or M) are given and implemented into 3D finite element (FEM) codes. Based on the use of FEM simulation in combination with least square optimization methods, the previously collected experimental data are fitted and all three energy functions ψB , ψH and ψM are computed. The obtained material model proves to have excellent predictive capabilities when compared to other experiments not used in the fitting process. The use of numerical tools is necessary to make sure that the calculated material parameters are not influenced by the shape of experimental specimens.The last part of this work details the numerical implementation of the different variational formulations. For each one of them, it is found that isoparametric elements are well suited to simulate coupled magneto-mechanical boundary value problems. We show that special care is needed when implementing variational formulations using the displacement vector and the magnetic vector potential as independent variables. Indeed, ensuring the uniqueness of the vector potential requires to numerically enforce the Coulomb gauge, which leads to numerical complications that are addressed in this thesis. Before describing the different patch tests that have been considered to validate the numerical codes, we show which are the valid boundary conditions for the magnetic vector potential and how to use the symmetry properties of a given boundary value problem to reduce its complexity and the computational resources needed to solve it.

The focus of this work is on the modeling and simulation of shock wave propagation in reactive metal powder mixtures. Reactive metal systems are non-explosive, solid-state materials that release chemical energy when subjected to sufficiently strong stimuli. Shock loading experiments have demonstrated that ultra-fast chemical reactions can be achieved in certain micron-sized metal powder mixtures. However, the mechanisms of rapid mixing that drive these chemical reactions are currently unclear. The goal of this research is to gain an understanding of the shock-induced deformation that enables these ultra-fast reactions. The problem is approached using direct numerical simulation. In this work, a finite element (FE) model is developed to simulate shock wave propagation in discrete particle mixtures. This provides explicit particle-level resolution of the thermal and mechanical fields that develop in the shock wave. The Ni/Al powder system has been selected for study. To facilitate mesoscale FE simulation, a new dislocation-based constitutive model has been developed to address the viscoplastic deformation of fcc metals at very high strain rates. Six distinct initial configurations of the Ni/Al powder system have been simulated to quantify the effects of powder configuration (e.g., particle size, phase morphology, and constituent volume fractions) on deformation in the shock wave. Results relevant to the degree of shock-induced mixing in the Ni/Al powders are presented, including specific analysis of the thermodynamic state and microstructure of the Ni/Al interfaces that develop during wave propagation. Finally, it is shown that velocity fluctuations at the Ni/Al interfaces (which arise due to material heterogeneity) may serve to fragment the particles down to the nanoscale, and thus provide an explanation of ultra-fast chemical reactions in these material systems.

Structures like transmission towers, tele-communication masts, dry-docks, tall chimneys, tunnels and burried pipelines under water etc are subjected to considerable uplift forces. The net effect of external loading on the foundations of these structures results in forces that try to pull the foundations out of the ground. Anchors are usually provided to resist such uplift forces. Earlier theoretical research of anchor behavior has focused on elastic response and ultimate pullout capacity. Many investigators have proposed techniques for determining the collapse load of anchors. Essentially the approaches involve the use of limit equilibrium concepts, with some assumptions regarding the shape of the failure surface and/or the influence of the soil above the anchor. The possible effect of dilatency and initial stress state are not considered in these methods. A number of investigators have used the results of small size model anchors to understand the behavior and extrapolated the results for predicting the behavior of large sized anchors. This has lead to unsatisfactory results. It has been clearly shown by Dickin (1989) that the failure displacements and load displacement curve patterns are very different for small and large sized anchors, i.e. they are not just proportional to the size of the anchor. Critical pullout load and the load displacement behavior are required for the complete analysis of anchor foundations. Though, many theories have been proposed to predict the uplift capacity within the limits of accuracy required at engineering level, at present no simple rational method is available for computing deformations. In the present investigation attempts have been made to analyze the load deformation behavior of large size strip anchors in sands, clays and layered soils using two-dimensional explicit finite difference program FLAG (Fast Lagrangian Analysis of Continua), well suited for geomaterials, by assuming soil to be a Mohr-Coulomb material in the case of sands and modified Cam-clay material in the case of clays. It is now well understood that the shearing resistance of a granular soil mass is derived from two factors frictional resistance and the dilatency of the soil. So the peak friction angle can be divided in to two components critical friction angle Фcv and dilation angle Ψ. Critical friction angle is the true friction angle as a result of frictional resistance at interparticle level when the soil is shearing at constant volume. If Фcv for a given soil remains constant, the value of Ψ has to increase with the increase in initial density of soil packing. The dilatency of a soil mass gradually decreases with continued shearing from its initial high value to zero after very large shear strains, when the soil finally reaches a constant, steady volume at critical states. Correspondingly the observed friction angle Ф reduces from its peak value to Фcv at a very large strain. In earlier days, clays used to be characterized by the strength parameters c and Ф. often, under undrained conditions, Ф would be even considered zero. But in the recent developments, it is understood that all the strength of clays is frictional. There is nothing like cohesion. The part of shear strength, which appears to be independent of normal stress, is shown to be the effect of over-consolidation and the resulting dilation. Thus although Cam-clay model uses zero cohesion for all clays, it reflects this component of strength through over-consolidation and in a more realistic way. Hence, it is appropriate to consider the pre-consolidation pressure as parameter in the analysis. More specifically, the various aspects covered in this investigation are as follows. Chapter 1 provides the general introduction. In chapter 2, the existing literature for the analysis of anchors for both experimental and analytical investigations on the pullout capacity of anchors in homogeneous and layered soils and the load deformation behavior of anchors under pullout are briefly reviewed. Chapter 3 deals with the features and the implementation of the two-dimensional explicit finite difference program, Fast Lagrangian Analysis of Continua (FLAC) and the constitutive modeling of soils. It discusses the background and implementation of Strain softening / hardening model. This model is based on the Mohr- Coulomb model with non-associated shear and associated tension flow rules. In this model the cohesion, friction, dilation and tensile strength may harden or soften after the onset of the plastic yield. Further the critical state concepts and implementation of the modified Cam-clay model have been discussed. Cam-clay model originally developed for clays reflects the hydrostatic pressure or density dependent hardening material response. Chapter 4 focuses on the analysis of load deformation behavior of large size anchors in granular soils. Two-dimensional explicit finite difference program (FLAC) is used for the simulations and the soil is modeled as a Mohr-Coulomb strain softening/hardening material In this chapter a series of simulations have been carried out on large size anchor plates, with parametric variation. By analyzing these results, a generalized load deformation relationship for different sizes of anchors and different types of soil have been proposed. The results are presented in the form of influence/design charts which can be used in hand calculations to obtain an estimate of anchor capacity and deformation for a wide range of soil types and size of anchors. Chapter 5 deals with the analysis of the drained and undrained behavior of large size horizontal strip anchors in clays using modified Cam-clay model. Earlier investigators have studied the undrained behavior of anchor plates in clays, but no studies are reported in literature for the drained behavior of anchors in clays. Further it is not clear whether, drained or undrained condition will be critical for an anchor. In this chapter the drained and undrained behavior of large size anchor plates in both normally consolidated and over-consolidated states have been made. It has been found that the undrained pullout capacity of an anchor in a soil of normally consolidated state will always be more than the drained capacity. This is contrast to the usual understanding that undrained behavior is more critical than the drained behavior. In Chapter 6 an attempt has been made to analyze the behavior of large size anchors in two layered sands and in conditions where backfill material has a higher or lower strength than the native soil, for different shape of excavations. Soil is assumed to be a Mohr-coulomb strain softening/hardening material. In Chapter 7 the entire investigation covered in earlier chapters has been synthesized and some specific conclusions have been highlighted.

Asphalt mixtures are complex composites that comprise aggregate, asphalt binder, and air. Several research studies have shown that the mechanical behavior of the asphalt mixture is strongly influenced by the matrix, i.e. the asphalt binder. Therefore, accurate constitutive models for the asphalt binders are critical to ensure accurate performance predictions at a material and structural level. However, researchers who use computational methods to model the micromechanics of asphalt mixtures typically assume that (i) asphalt binders behave linearly in shear, and (ii) either bulk modulus or Poisson’s ratio of asphalt binders is not time dependent. This research develops an approach to measure and model the shear and bulk behavior of asphalt binders at intermediate temperatures. First, this research presents the findings from a systematic investigation into the nature of the linear and nonlinear response of asphalt binders subjected to shear using a Dynamic Shear Rheometer (DSR). The DSR test results showed that under certain conditions a compressive normal force was generated in an axially constrained specimen subjected to cyclic torque histories. This normal force could not be solely attributed to the Poynting effect and was also related to the tendency of the asphalt binder to dilate when subjected to shear loads. The generated normal force changed the state of stress and interacted with the shear behavior of asphalt binder. This effect was considered to be an “interaction nonlinearity” or “three dimensional effect”. A constitutive model was identified to accommodate this effect. The model was successfully validated for several different loading histories. Finally, this study investigated the time-dependence of the bulk modulus of asphalt binders. To this end, poker-chip geometries with high aspect ratios were used. The boundary value problem for the poker-chip geometry under step displacement loading was solved to determine the bulk modulus and Poisson’s ratio of asphalt binders as a function of time. The findings from this research not only improve the understanding of asphaltic materials behavior, but also provide tools required to accurately predict pavement performance. / text

Previous experimental observations have shown that the pseudoelastic response of NiTi shape memory alloys (SMA) is localized in nature and proceeds through nucleation and propagation of localized deformation bands. It has also been observed that the mechanical response of SMAs is strongly affected by loading rate and cyclic degradation. These behaviors significantly limit the accurate modeling of SMA elements used in various devices and applications. The aim of this work is to provide engineers with a constitutive model that can accurately describe the dynamic, unstable pseudoelastic response of SMAs, including their cyclic response, and facilitate the reliable design of SMA elements. A 1-D phenomenological model is developed to simulate the localized phase transformations in NiTi wires during both loading and unloading. In this model, it is assumed that the untransformed particles located close to the transformed regions are less stable than those further away from the transformed regions. By consideration of the thermomechanical coupling among the stress, temperature, and latent heat of transformation, the analysis can account for strain-rate effects. Inspired by the deformation theory of plasticity, the 1-D model is extended to a 3-D macromechanical model of localized unstable pseudoelasticity. An important feature of this model is the reorientation of the transformation strain tensor with changes in stress tensor. Unlike previous modeling efforts, the present model can also capture the propagation of localized deformation during unloading. The constitutive model is implemented within a 2-D finite element framework to allow numerical investigation of the effect of strain rate and boundary conditions on the overall mechanical response and evolution of localized transformation bands in NiTi strips. The model successfully captures the features of the transformation front morphology, and pseudoelastic response of NiTi strip samples observed in previous experiments. The 1-D and 3-D constitutive models are further extended to include the plastic deformation and degradation of material properties as a result of cyclic loading.