A constitutive model of the deformation behaviour of twinning induced plasticity (TWIP) steel at different temperatures

The mechanical behaviour of Fe-18Mn-0.6C-1Al (wt%) TWIP steel was modelled in the temperature range from room temperature to 400°C. The proposed constitutive model was based on the Kocks-Mecking-Estrin (KME) model. The model parameters were determined using extensive experimental measurements of the physical parameters such as the dislocation mean free path and the volume fraction of twinned grains. More than 100 grains with a total area of ~300μm2 were examined at different strain levels over the entire stress-strain curve. Uniaxial tensile deformation of the TWIP steel was modelled for different deformation temperatures using a modelling approach which considers two distinct populations of grains: twinned and twin-free ones. A key point of the work was a meticulous experimental determination of the evolution of the volume fraction of twinned grains during uniaxial tensile deformation. This information was implemented in a phase-mixture model that yielded a very good agreement with the experimental tensile behaviour for the tested range of deformation temperatures. © 2014 Elsevier B.V.

On the strength of dislocation interactions and their effect on latent hardening in pure Magnesium

This study is dedicated to the quantification of latent hardening and its effect on the plasticity of pure hexagonal magnesium. To this end, discrete dislocation dynamics simulations are used to (1) extract latent hardening parameters coupling different slip systems, and to (2) assess the validity of two existing constitutive models linking slip system strength to dislocation densities on all slip systems. As hexagonal materials deform via activation of different slip modes, each with different mobilities and lattice friction stress, the effects of the latter on latent hardening evolution are also investigated. It is found that the multi-slip formulation proposed by Franciosi and Zaoui gives accurate predictions when multiple interactions are involved while the formulation suggested by Lavrentev and Pokhil systematically overestimates the flow stress. Similar to FCC materials, it is also found that collinear interactions potentially contribute the most to latent hardening. Basal/pyramidal c + a interactions are found to be very strong, while interactions involving second-order pyramidal c + a primary dislocations appear to be the weakest ones. Finally, the latent hardening parameters, extracted from the discrete dislocation dynamics simulations, are used in polycrystal simulations and the impact of finely accounting for latent hardening on predictions of the macroscopic anisotropic response is shown to be of significant importance.

A microstructure based analytical model for tensile twinning in a rod textured Mg alloy

Abstract A model for tensile twinning during the compression of rod textured magnesium is developed based on the idea that these twins nucleate at grain boundaries and that when the twin number density per grain is low these twins readily give rise to the formation of other 'interaction' twins in adjacent grains. Experimental observations of twin aspect ratios measured at a single grain size and twin number densities measured over four grain sizes were used to determine model material parameters. Using these, the model provides reasonable predictions for the observed magnitudes and trends for the following observations:Effect of grain size and stress on twin volume fraction, fractional twin length and the fraction of twin contact.Effect of grain size on the yield stress.Effect of grain size on the general shape of the stress-strain curve at low strains. A parametric study shows the model to be quite robust but that it is particularly sensitive to the value of the exponent assumed for the twin nucleation rate law. It is seen that preventing the formation of interaction twins provides an important avenue for hardening and that the flow stress is also particularly sensitive to the relaxation of the twin back stresses. The model shows the importance of taking microstructure into account when modelling twinning.

Constitutive analysis of hot deformation behavior of a Ti6Al4V alloy using physical based model

The effect of deformation parameters on the flow behavior of a Ti6Al4V alloy has been studied to understand the deformation mechanisms during hot compression. Cylindrical samples with partially equiaxed grains were deformed in the α+β phase region at different thermo-mechanical conditions. To develop components with tailored properties, the physically based Estrin and Mecking (EM) model for the work hardening/dynamic recovery combined with the Avrami equation for dynamic recrystallization was used to predict the flow stress at varying process conditions. The EM model revealed good predictability up to the peak strain, however, at strain rates below 0.01s^-1, a higher B value was observed due to the reduced density of dislocation tangles. In contrast, the flow softening model revealed higher value of constants a and b at high strain rates due to the reduction in the volume fraction of dynamic recrystallization and larger peak strain. The predicted flow stress using the combined EM+Avrami model revealed good agreement with the measured flow stress resulted in very low average absolute relative error value. The microstructural analysis of the samples suggests the formation of coarse equiaxed grains together with the increased β phase fraction at low strain rate leads to a higher flow softening.

A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic - Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Latent hardening size effect in small-scale plasticity

We aim at understanding the multislip behaviour of metals subject to irreversible deformations at small-scales. By focusing on the simple shear of a constrained single-crystal strip, we show that discrete Dislocation Dynamics (DD) simulations predict a strong latent hardening size effect, with smaller being stronger in the range [1.5 µm, 6 µm] for the strip height. We attempt to represent the DD pseudo-experimental results by developing a flow theory of Strain Gradient Crystal Plasticity (SGCP), involving both energetic and dissipative higher-order terms and, as a main novelty, a strain gradient extension of the conventional latent hardening. In order to discuss the capability of the SGCP theory proposed, we implement it into a Finite Element (FE) code and set its material parameters on the basis of the DD results. The SGCP FE code is specifically developed for the boundary value problem under study so that we can implement a fully implicit (Backward Euler) consistent algorithm. Special emphasis is placed on the discussion of the role of the material length scales involved in the SGCP model, from both the mechanical and numerical points of view.

Latent hardening size effect in small-scale plasticity

We aim at understanding the multislip behaviour of metals subject to irreversible deformations at small-scales. By focusing on the simple shear of a constrained single-crystal strip, we show that discrete Dislocation Dynamics (DD) simulations predict a strong latent hardening size effect, with smaller being stronger in the range [1.5 µm, 6 µm] for the strip height. We attempt to represent the DD pseudo-experimental results by developing a flow theory of Strain Gradient Crystal Plasticity (SGCP), involving both energetic and dissipative higher-order terms and, as a main novelty, a strain gradient extension of the conventional latent hardening. In order to discuss the capability of the SGCP theory proposed, we implement it into a Finite Element (FE) code and set its material parameters on the basis of the DD results. The SGCP FE code is specifically developed for the boundary value problem under study so that we can implement a fully implicit (Backward Euler) consistent algorithm. Special emphasis is placed on the discussion of the role of the material length scales involved in the SGCP model, from both the mechanical and numerical points of view.

Properties of granular analogue model materials: A community wide survey

We report the material properties of 26 granular analogue materials used in 14 analogue modelling laboratories. We determined physical characteristics such as bulk density, grain size distribution, and grain shape, and performed ring shear tests to determine friction angles and cohesion, and uniaxial compression tests to evaluate the compaction behaviour. Mean grain size of the materials varied between c. 100 and 400 μm. Analysis of grain shape factors shows that the four different classes of granular materials (14 quartz sands, 5 dyed quartz sands, 4 heavy mineral sands and 3 size fractions of glass beads) can be broadly divided into two groups consisting of 12 angular and 14 rounded materials. Grain shape has an influence on friction angles, with most angular materials having higher internal friction angles (between c. 35° and 40°) than rounded materials, whereas well-rounded glass beads have the lowest internal friction angles (between c. 25° and 30°). We interpret this as an effect of intergranular sliding versus rolling. Most angular materials have also higher basal friction angles (tested for a specific foil) than more rounded materials, suggesting that angular grains scratch and wear the foil. Most materials have an internal cohesion in the order of 20–100 Pa except for well-rounded glass beads, which show a trend towards a quasi-cohesionless (C < 20 Pa) Coulomb-type material. The uniaxial confined compression tests reveal that rounded grains generally show less compaction than angular grains. We interpret this to be related to the initial packing density after sifting, which is higher for rounded grains than for angular grains. Ring-shear test data show that angular grains undergo a longer strain-hardening phase than more rounded materials. This might explain why analogue models consisting of angular grains accommodate deformation in a more distributed manner prior to strain localisation than models consisting of rounded grains.

Modeling of hardening and softening processes in Mg alloys

The tensile deformation behavior of a range of supersaturated Mg-Al solid solutions and an as-cast magnesium alloy AM60 has been studied. The Mg-Al alloys were tested at room temperature while the alloy AM60 was tested in the temperature range 293-573 K. The differences in the deformation behavior of the alloys is discussed in terms of hardening and softening processes. In order to identify which processes were active, the stress dependence of the strain-hardening coefficient was assessed using Lukac and Balik's model of hardening and softening. The analysis indicates that hardening involves solid solution hardening and interaction with forest dislocations and non-dislocation obstacles such as second phase particles. Cross slip is not a significant recovery process in the temperature range 293-423 K. At temperatures between 473 and 523 K the analysis suggests that softening is controlled by cross slip and climb of dislocations. At temperatures above 523 K softening seems to be controlled by dynamic recrystallisation. (C) 2004 Elsevier B.V. All rights reserved.

Aggregation and caking processes of granular materials:continuum model and simulation with application to sugar

Aggregation and caking of particles are common severe problems in many operations and processing of granular materials, where granulated sugar is an important example. Prevention of aggregation and caking of granular materials requires a good understanding of moisture migration and caking mechanisms. In this paper, the modeling of solid bridge formation between particles is introduced, based on moisture migration of atmospheric moisture into containers packed with granular materials through vapor evaporation and condensation. A model for the caking process is then developed, based on the growth of liquid bridges (during condensation), and their hardening and subsequent creation of solid bridges (during evaporation). The predicted caking strengths agree well with some available experimental data on granulated sugar under storage conditions.

A mechanistic description of combined hardening and size effects

Mechanistic models based on geometrically necessary dislocations are re-examined in light of recent experiments exhibiting the indentation size effect. A simple method is developed to combine work hardening, solid solution hardening, radiation hardening and size effects. The model is verified by experiments in ionic salt crystals. © 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

A Finite Element to Perform 3D Shakedown Analysis for Limited Kinematic Hardening Materials

The use of the Design by Analysis (DBA) route is a modern trend in pressure vessel and piping international codes in mechanical engineering. However, to apply the DBA to structures under variable mechanical and thermal loads, it is necessary to assure that the plastic collapse modes, alternate plasticity and incremental collapse (with instantaneous plastic collapse as a particular case), be precluded. The tool available to achieve this target is the shakedown theory. Unfortunately, the practical numerical applications of the shakedown theory result in very large nonlinear optimization problems with nonlinear constraints. Precise, robust and efficient algorithms and finite elements to solve this problem in finite dimension has been a more recent achievements. However, to solve real problems in an industrial level, it is necessary also to consider more realistic material properties as well as to accomplish 3D analysis. Limited kinematic hardening, is a typical property of the usual steels and it should be considered in realistic applications. In this paper, a new finite element with internal thermodynamical variables to model kinematic hardening materials is developed and tested. This element is a mixed ten nodes tetrahedron and through an appropriate change of variables is possible to embed it in a shakedown analysis software developed by Zouain and co-workers for elastic ideally-plastic materials, and then use it to perform 3D shakedown analysis in cases with limited kinematic hardening materials

Shakedown Analysis Considering Limited Kinematic Hardening Materials

In design or safety assessment of mechanical structures, the use of the Design by Analysis (DBA) route is a modern trend. However, for making possible to apply DBA to structures under variable loads, two basic failure modes considered by ASME or European Standards must be precluded. Those modes are the alternate plasticity and incremental collapse (with instantaneous plastic collapse as a particular case). Shakedown theory is a tool that permit us to assure that those kinds of failures will be avoided. However, in practical applications, very large nonlinear optimization problems are generated. Due to this facts, only in recent years have been possible to obtain algorithms sufficiently accurate, robust and efficient, for dealing with this class of problems. In this paper, one of these shakedown algorithms, developed for dealing with elastic ideally-plastic structures, is enhanced to include limited kinematic hardening, a more realistic material behavior. This is done in the continuous model by using internal thermodynamic variables. A corresponding discrete model is obtained using an axisymmetric mixed finite element with an internal variable. A thick wall sphere, under variable thermal and pressure loads, is used in an example to show the importance of considering the limited kinematic hardening in the shakedown calculations