Dynamics of strings made of phase-transforming materials

This paper presents a theory to describe the dynamical behavior of a string made of a phase-transforming material like a shape-memory alloy. The study of phase boundaries, the driving force acting on them and the kinetic relation governing their propagation is of central concern. The paper proposes a qualitative experimental test of the notion of a kinetic relation, as well as a simple experimental method for measuring it quantitatively. It presents a numerical method for studying general initial and boundary value problems in strings, and concludes by exploring the use of phase transforming strings to generate motion at very small scales.     

Experimental study of the phase transformation plasticity of 16MND5 low carbon steel under multiaxial loading

This paper is concerned with the experimental behaviour of a 16MND5 steel (french vessel steel) under complex loading. A particular attention is paid to plasticity induced by phase transformation. We present an experimental set-up to apply thermo-mechanical loads under tension-torsion. This apparatus enables us to reach temperature of 1200 °C at a maximum heating rate of 60 °C s⁻¹ and a high cooling rate of −30 °C s⁻¹. A series of tests is performed in order to show the rule of loading on transformation plasticity.     

A large strain thermoviscoplastic formulation for the solidification of S.G. cast iron in a green sand mould

This paper presents a large strain thermoviscoplastic formulation for the analysis of the solidification process of spheroidal graphite (S.G.) cast iron in a green sand mould. This formulation includes two different non-associate constitutive models in order to describe the thermomechanical behaviour of each of such materials during the whole process. The performance of these models is evaluated in the analysis of a solidification test.     

A model for plasticity kinetics and its role in simulating the dynamic behavior of Fe at high strain rates

The recent diagnostic capability of the Omega laser to study solid-solid phase transitions at pressures greater than 10 GPa and at strain rates exceeding 10⁷ s⁻¹ has also provided valuable information on the dynamic elastic-plastic behavior of materials. We have found, for example, that plasticity kinetics modifies the effective loading and thermodynamic paths of the material. In this paper we derive a kinetics equation for the time-dependent plastic response of the material to dynamic loading, and describe the model’s implementation in a radiation-hydrodynamics computer code. This model for plasticity kinetics incorporates the Gilman model for dislocation multiplication and saturation. We discuss the application of this model to the simulation of experimental velocity interferometry data for experiments on Omega in which Fe was shock compressed to pressures beyond the α-to-ε phase transition pressure. The kinetics model is shown to fit the data reasonably well in this high strain rate regime and further allows quantification of the relative contributions of dislocation multiplication and drag. The sensitivity of the observed signatures to the kinetics model parameters is presented.     

A multi-axial, multimechanism based constitutive model for the comprehensive representation of the evolutionary response of SMAs under general thermomechanical loading conditions

We present a fully general, three dimensional, constitutive model for Shape Memory Alloys (SMAs), aimed at describing all of the salient features of SMA evolutionary response under complex thermomechanical loading conditions. In this, we utilize the mathematical formulation we have constructed, along with a single set of the model’s material parameters, to demonstrate the capturing of numerous responses that are experimentally observed in the available SMA literature. This includes uniaxial, multi-axial, proportional, non-proportional, monotonic, cyclic, as well as other complex thermomechanical loading conditions, in conjunction with a wide range of temperature variations. The success of the presented model is mainly attributed to the following two main factors. First, we use multiple inelastic mechanisms to organize the exchange between the energy stored and energy dissipated during the deformation history. Second, we adhere strictly to the well established mathematical and thermodynamical requirements of convexity, associativity, normality, etc. in formulating the evolution equations governing the model behavior, written in terms of the generalized internal stress/strain tensorial variables associated with the individual inelastic mechanisms. This has led to two important advantages: (a) it directly enabled us to obtain the limiting/critical transformation surfaces in the spaces of both stress and strain, as importantly required in capturing SMA behavior; (b) as a byproduct, this also led, naturally, to the exhibition of the apparent deviation from normality, when the transformation strain rate vectors are plotted together with the surfaces in the space of external/global stresses, that has been demonstrated in some recent multi-axial, non-proportional experiments.     

The dynamic growth of a single void in a viscoplastic material under transient hydrostatic loading

We have examined the problem of the dynamic growth of a single spherical void in an elastic-viscoplastic medium, with a view towards addressing a number of problems that arise during the dynamic failure of metals. Particular attention is paid to inertial, thermal and rate-dependent effects, which have not previously been thoroughly studied in a combined setting. It is shown that the critical stress for unstable growth of the void in the quasistatic case is strongly affected by the thermal softening of the material (in adiabatic calculations). Thermal softening has the effect of lowering the critical stress, and has a stronger influence at high strain hardening exponents. It is shown that the thermally diffusive case for quasistatic void growth in rate-dependent materials is strongly affected by the initial void size, because of the length scale introduced by the thermal diffusion. The effects of inertia are quantified, and it is demonstrated that inertial effects are small in the early stages of void growth and are strongly dependent on the initial size of the void and the rate of loading. Under supercritical loading for the inertial problem, voids of all sizes achieve a constant absolute void growth rate in the long term. Inertia first impedes, but finally promotes dynamic void growth under a subcritical loading. For dynamic void growth, the effect of rate-hardening is to reduce the rate of void growth in comparison to the rate-independent case, and to reduce the final relative void growth achieved.     

Mechanical wave momentum from the first principles

Axial momentum carried by waves in a uniform waveguide is considered based on the conservation laws and a kind of the causality principle. Specifically, we examine (without resorting to constitutive data) steady-state waves of an arbitrary shape, periodic waves which speed differs from the speed of its form and binary waves carrying self-equilibrated momentum. The approach allows us to represent momentum as a product of the wave mass and the wave speed. The propagating wave mass, positive or negative, is the excess of that in the wave over its initial value. This general representation is valid for mechanical waves of arbitrary nature and intensity. The finite-amplitude longitudinal and periodic transverse waves are examined in more detail. It is shown in particular, that the transverse excitation of a string or an elastic beam results in the binary wave. The closed-form expressions for the drift in these waves functionally reduce to the Stokes’ drift in surface water waves (a half the latter by the value). Besides, based on the general representation an energy–momentum relation is discussed and the physical meaning of the so-called “wave momentum” is clarified.     

Dissipative interface waves and the transient response of a three-dimensional sliding interface with Coulomb friction

We investigate the linearized response of two elastic half-spaces sliding past one another with constant Coulomb friction to small three-dimensional perturbations. Starting with the assumption that friction always opposes slip velocity, we derive a set of linearized boundary conditions relating perturbations of shear traction to slip velocity. Friction introduces an effective viscosity transverse to the direction of the original sliding, but offers no additional resistance to slip aligned with the original sliding direction. The amplitude of transverse slip depends on a nondimensional parameter η=c_sτ₀/μv₀, where τ₀ is the initial shear stress, 2v₀ is the initial slip velocity, μ is the shear modulus, and c_s is the shear wave speed. As η→0, the transverse shear traction becomes negligible, and we find an azimuthally symmetric Rayleigh wave trapped along the interface. As η→∞, the inplane and antiplane wavesystems frictionally couple into an interface wave with a velocity that is directionally dependent, increasing from the Rayleigh speed in the direction of initial sliding up to the shear wave speed in the transverse direction. Except in these frictional limits and the specialization to two-dimensional inplane geometry, the interface waves are dissipative. In addition to forward and backward propagating interface waves, we find that for η>1, a third solution to the dispersion relation appears, corresponding to a damped standing wave mode. For large-amplitude perturbations, the interface becomes isotropically dissipative. The behavior resembles the frictionless response in the extremely strong perturbation limit, except that the waves are damped. We extend the linearized analysis by presenting analytical solutions for the transient response of the medium to both line and point sources on the interface. The resulting self-similar slip pulses consist of the interface waves and head waves, and help explain the transmission of forces across fracture surfaces. Furthermore, we suggest that the η→∞ limit describes the sliding interface behind the crack edge for shear fracture problems in which the absolute level of sliding friction is much larger than any interfacial stress changes.     

Onset of necking in electro-magnetically formed rings

The electromagnetic forming (EMF) process is an attractive manufacturing technique, which uses electromagnetic (Lorentz) body forces to fabricate metallic parts. One of the many advantages of EMF is the considerable ductility increase observed in several metals, with aluminum featuring prominently among them. The quantitative explanation of this phenomenon is the primary motivation of this work.To study the ductility increase due to EMF, an aluminum ring is placed outside a fixed coil, which is connected to a capacitor. Upon the capacitor's discharge, the time varying current in the coil induces a larger current in the ring specimen and the resulting Lorentz forces make it expand. The coupled coil-ring electromagnetic and thermomechanical problem is solved, using an experimentally obtained constitutive model for a particular aluminum alloy. Our results show that for realistic imperfections, the EMF ring starts necking at strains about six times larger than its static counterpart, as observed experimentally. This study establishes the importance of solving the fully coupled electromagnetic and thermomechanical problem and provides insight on how different constitutive parameters influence ductility in an EMF process.     

A lattice-based model of the kinetics of twin boundary motion

In this paper we derive a macroscopic kinetic law for twin boundary motion from a lattice dynamical model. The model is developed for compound and type-1 twins and it is explicitly illustrated for a Cu-Al-Ni shape memory alloy. The governing multiple-well energy is calculated using an effective interatomic potential; a Frenkel-Kontorowa type model is developed for the dynamics at the lattice scale; and a quasi-continuum approximation is used to determine the continuum-scale kinetics. The model predicts that compound twins in the Cu-Al-Ni shape memory alloy are an order of magnitude more mobile than type-1 twins which is consistent with experimental observations.     

The plastic deformation of non-homogeneous polycrystals

To describe the work hardening process of polycrystals processed using various thermomechanical cycles with isochronal annealing from 500 to 900 °C, a dislocation based strain hardening model constructed in the basis of the so-called Kocks–Mecking model is proposed. The time and temperature dependence of flow stress is accounted via grain boundary migration, and the migration is related to annihilation of extrinsic grain boundary dislocations (EGBD’s) by climb via lattice diffusion of vacancies at the triple points. Recovery of yield stress is associated with changes in the total dislocation density term ρ^T. A sequence of deformation and annealing steps generally result in reduction of flow stress via the annihilation of the total dislocation density ρ^T defined as the sum of geometrically necessary dislocations ρ^G and statistically stored dislocations ρ^S. The predicted variation of yield stress with annealing temperature and cold working stages is in agreement with experimental observations. An attempt is made to determine the mathematical expressions which best describe the deformation behaviour of polycrystals in uniaxial deformation.     

A computational model for the indentation and phase transformation of a martensitic thin film

We propose a computational model for a stress-induced martensitic phase transformation of a single-crystal thin film by indentation and its reverse transformation to austenite by heating. Our model utilizes a surface energy that allows sharp interfaces with finite energy and a penalty that forces the film to lie above the indenter and undergo a stress-induced austenite-to-martensite phase transformation. We introduce a method to nucleate the martensite-to-austenite phase transformation since in our model the film would otherwise remain in the martensitic phase in a local minimum of the energy.     

An analytical model of rumpling in thermal barrier coatings

Multilayer thermal barrier coatings (TBCs) deposited on superalloy turbine blades provide protection from combustion temperatures in excess of 1500 °C. One of the dominant failure modes comprises cracking from undulation growth, or rumpling, of the highly compressed oxide layer that grows between the ceramic top coat and the intermetallic bond coat. In this paper, a mechanistic model providing an analytical approximation of undulation growth is presented for realistic cyclic thermal histories. Thickening, lateral growth straining and high temperature yielding of the oxide layer are taken into account. Undulation growth in TBC systems is highly nonlinear and characterized by more than 20 material and geometric parameters, highlighting the importance of a robust yet computationally efficient model. At temperatures above 600 °C, the bond coat creeps. Thermal expansion mismatch occurs between the superalloy substrate and the oxide layer and, in some systems, the bond coat. In addition, some bond coats, such as PtNiAl, exhibit a martensitic phase transformation accompanied by nearly a 1% linear expansion, giving rise to a large effective mismatch. These two mismatches promote undulation growth. Nonlinear interaction between the stress in the bond coat induced by the constraining effect of the thick substrate and normal tractions applied at the surface of the bond coat by the compressed, undulating oxide layer produces an increment of undulation growth during each thermal cycle, before the stress decays by creep. A series of problems for systems without the ceramic top coat are used to elucidate the mechanics of undulation growth and to replicate trends observed in a series of experiments and in prior finite-element simulations. The model is employed to study for the first time the effect on undulation growth of a shift in the temperature range over which the transformation occurs, as well as the relative importance of the transformation compared to thermal expansion mismatch. The role of the top coat and other viable ways of reducing undulation growth are considered.     

An experimental investigation of shock wave propagation in periodically layered composites

In heterogeneous media, scattering due to interfaces/microstructure between dissimilar materials could play an important role in shock wave dissipation and dispersion. In this work, the influence of interface scattering on finite-amplitude shock waves was experimentally investigated by impacting flyer plates onto periodically layered polycarbonate/6061 aluminum, polycarbonate/304 stainless steel and polycarbonate/glass composites. Experimental results (obtained using velocity interferometer and stress gage) show that these periodically layered composites can support steady structured shock waves. Due to interface scattering, the effective shock viscosity increases with the increase of interface impedance mismatch, and decreases with the increase of interface density (interface area per unit volume) and loading amplitude. For the composites studied here, the strain rate within the shock front is roughly proportional to the square of the shock stress. This indicates that layered composites have much larger shock viscosity due to the interface/microstructure scattering in comparison with the increase of shock strain rate by the fourth power of the shock stress for homogeneous metals. Experimental results also show that due to the scattering effects, shock propagation in the layered composites is dramatically slowed down and the shock speed in composites can be lower than that either of its components.     

Multiaxial constitutive model accounting for the strength-differential in Inconel 718

The nickel-base alloy Inconel 718 exhibits a strength-differential, that is, a different plastic flow behavior in uniaxial tension and uniaxial compression. A phenomenological viscoplastic model founded on thermodynamics has been extended for material behavior that deviates from classical metal plasticity by including all three stress invariants in the threshold function. The model can predict plastic flow in isotropic materials with or without a flow stress asymmetry as well as with or without pressure dependence. Viscoplastic material parameters have been fit to pure shear, uniaxial tension, and uniaxial compression experimental results at 650°°C. Threshold function material parameters have been fit to the strength-differential. Four classes of threshold functions have been considered and nonproportional loading of hollow tubes, such as shear strain followed by axial strain, has been used to select the most applicable class of threshold function for the multiaxial model as applied to Inconel 718 at 650 °C. These nonproportional load paths containing corners provide a rigorous test of a plasticity model, whether it is time-dependent or not. A J₂J₃ class model, where J₂ and J₃ are the second and third effective deviatoric stress invariants, was found to agree the best with the experimental results.     

On probabilistic aspects in the dynamic degradation of ductile materials

Dynamic loadings produce high stress waves leading to the spallation of ductile materials such as aluminum, copper, magnesium or tantalum. The main mechanism used herein to explain the change of the number of cavities with the stress rate is nucleation inhibition, as induced by the growth of already nucleated cavities. The dependence of the spall strength and critical time with the loading rate is investigated in the framework of a probabilistic model. The present approach, which explains previous experimental findings on the strain-rate dependence of the spall strength, is applied to analyze experimental data on tantalum.     

Rippling and a phase-transforming mesoscopic model for multiwalled carbon nanotubes

We propose to model thick multiwalled carbon nanotubes as beams with non-convex curvature energy. Such models develop stressed phase mixtures composed of smoothly bent sections and rippled sections. This model is motivated by experimental observations and large-scale atomistic-based simulations. The model is analyzed, validated against large-scale simulations, and exercised in examples of interest. It is shown that modelling MWCNTs as linear elastic beams can result in poor approximations that overestimate the elastic restoring force considerably, particularly for thick tubes. In contrast, the proposed model produces very accurate predictions both of the restoring force and of the phase pattern. The size effect in the bending response of MWCNTs is also discussed.