A constitutive relation for compressible plastic materials

Summary In this note a yield criterion and its associated stress-strain increment relations dependent on isotropic pressure, are presented, for a non-hardening elastic-plastic compressible material.     

Damage of ductile materials deforming under multiple plastic or viscoplastic mechanisms

This work presents a model to represent ductile failure (i.e. failure controlled by nucleation, growth and coalescence) of materials whose irreversible deformation is controlled by several plastic or viscoplastic deformation mechanisms. In addition work hardening may result from both isotropic and kinematic hardening. Damage is represented by a single variable representing void volume fraction. The model uses an additive decomposition of the plastic strain rate tensor. The model is developed based on the definition of damage dependant effective scalar stresses. The model is first developed within the generalized standard material framework and expressions for Helmholtz free energy, yield potential and dissipation potential are proposed. In absence of void nucleation, the evolution of the void volume fraction is governed by mass conservation and damage does not need to be represented by state variables. The model is extended to account for void nucleation. It is implemented in a finite element software to perform structural computations. The model is applied to three case studies: (i) failure by void growth and coalescence by internal necking (pipeline steel) where plastic flow is either governed by the Gurson–Tvergaard–Needleman model or the Thomason model, (ii) creep failure (Grade 91 creep resistant steel) where viscoplastic flow is controlled by dislocation creep or diffusional creep and (iii) ductile rupture after pre-compression (aluminum alloy) where kinematic hardening plays an important role.     

A strain localization analysis using a viscoplastic softening model for clay

Strain localization has become an attractive subject in geomechanics during the past decade. Shear bands are well known to develop in clay specimens during the straining process. Strain localization is closely related to plastic instability. In the present paper, a non-linear instability condition for the viscoplastic strain softening model during the creep process is firstly obtained. It is found that the proposed viscoplastic model is capable of describing plastic instability. Secondly, a two-dimensional linear instability analysis is performed and the preferred orientation for the growth of fluctuation and the instability condition are derived. It is worth noting that the two instability conditions are equivalent. Finally, the behavior of the clay is numerically analyzed in undrained plane-strain compression tests by the finite element method, considering a transport of pore water in the material at a quasi-static strain rate. The numerical results show that the model can predict strain localization phenomena, such as shear banding. From the numerical calculations, the effects of strain rate and permeability are discussed.     

Modeling the viscoplastic micromechanical response of two-phase materials using Fast Fourier Transforms

A viscoplastic approach using the Fast Fourier Transform (FFT) method for obtaining local mechanical response is utilized to study microstructure-property relationships in composite materials. Specifically, three-dimensional, two-phase digital materials containing isotropically coarsened particles surrounded by a matrix phase, generated through a Kinetic Monte Carlo Potts model for Ostwald ripening, are used as instantiations in order to calculate the stress and strain-rate fields under uniaxial tension. The effects of the morphology of the matrix phase, the volume fraction and the contiguity of particles, and the polycrystallinity of matrix phase, on the stress and strain-rate fields under uniaxial tension are examined. It is found that the first moments of the stress and strain-rate fields have a different dependence on the particle volume fraction and the particle contiguity from their second moments. The average stresses and average strain-rates of both phases and of the overall composite have rather simple relationships with the particle volume fraction whereas their standard deviations vary strongly, especially when the particle volume fraction is high, and the contiguity of particles has a noticeable effect on the mechanical response. It is also found that the shape of stress distribution in the BCC hard particle phase evolves as the volume fraction of particles in the composite varies, such that it agrees with the stress field in the BCC polycrystal as the volume of particles approaches unity. Finally, it is observed that the stress and strain-rate fields in the microstructures with a polycrystalline matrix are less sensitive to changes in volume fraction and contiguity of particles.     

A viscoplastic strain gradient analysis of materials with voids or inclusions

A finite strain viscoplastic nonlocal plasticity model is formulated and implemented numerically within a finite element framework. The model is a viscoplastic generalisation of the finite strain generalisation by Niordson and Redanz (2004) [Journal of the Mechanics and Physics of Solids 52, 2431–2454] of the strain gradient plasticity theory proposed by Fleck and Hutchinson (2001) [Journal of the Mechanics and Physics of Solids 49, 2245–2271]. The formulation is based on a viscoplastic potential that enables the formulation of the model so that it reduces to the strain gradient plasticity theory in the absence of viscous effects. The numerical implementation uses increments of the effective plastic strain rate as degrees of freedom in addition to increments of displacement. To illustrate predictions of the model, results are presented for materials containing either voids or rigid inclusions. It is shown how the model predicts increased overall yield strength, as compared to conventional predictions, when voids or inclusions are in the micron range. Furthermore, it is illustrated how the higher order boundary conditions at the interface between inclusions and matrix material are important to the overall yield strength as well as the material hardening.     

Transient squeeze flow of viscoplastic materials

The transient, axisymmetric squeezing of viscoplastic materials under creeping flow conditions is examined. The flow of the material even outside the disks is followed. Both cases of the disks moving with constant velocity or under constant force are studied. This time-dependent simulation of squeeze flow is performed for such materials in order to determine very accurately the evolution of the force or the velocity, respectively, and the distinct differences between these two experiments, the highly deforming shape and position of all the interfaces, the effect of possible slip on the disk surface, especially when the slip coefficient is not constant, and the effect of gravity. All these are impossible under the quasi-steady state condition used up to now. The exponential constitutive model, suggested by Papanastasiou, is employed. The governing equations are solved numerically by coupling the mixed finite element method with a quasi-elliptic mesh generation scheme in order to follow the large deformations of the free surface of the fluid. As the Bingham number increases, large departures from the corresponding Newtonian solution are found. When the disks are moving with constant velocity, unyielded material arises only around the two centers of the disks verifying previous works in which quasi-steady state conditions were assumed. The size of the unyielded region increases with the Bingham number, but decreases as time passes and the two disks approach each other. Their size also decreases as the slip velocity or the slip length along the disk wall increase. The force that must be applied on the disks in order to maintain their constant velocity increases significantly with the Bingham number and time and provides a first method to calculate the yield stress. On the other hand, when a constant force is applied on the disks, they slow down until they finally stop, because all the material between them becomes unyielded. The final location of the disk and the time when it stops provide another, probably easier, method to deduce the yield stress of the fluid.     

Crossflow of viscoplastic materials through tube bundles

The flow of viscoplastic materials through staggered arrays of tubes is analyzed. The mechanical behavior of the materials is assumed to obey the generalized Newtonian liquid (GNL) model, with a viscosity function given by the biviscosity law. The governing equations of this flow are solved numerically using a finite-volume method with a non-orthogonal mesh. For a representative range of the relevant parameters, results are presented in the form of velocity, pressure and viscosity fields. The pressure drop is also given as a function of rheological and geometric parameters.     

Upper yield points of non-linear viscoplastic materials

An approximate solution to the non-linear differential equation governing the behavior of a non-linear, viscoplastic material in tension is given. This solution is used to generate an accurate estimate for the upper yield point for the material. The results are compared to those previously obtained by other means.     

On the elastic–viscoplastic behavior of nanocrystalline materials

A new constitutive law is introduced to quantify the macroscopic effect of grain boundary dislocation emission on the behavior of pure face center cubic nanocrystalline materials. It is postulated that an emitted dislocation ends its trajectory in the grain boundary opposite to the source causing mass transfer. Dislocation emission by grain boundary ledges, considered here as the primary grain-boundary sources, is modeled as a thermally activated mechanism and the penetration of an emitted dislocation is assimilated as a soft collision. The macroscopic behavior of the material is retrieved via the use of a secant self-consistent scheme. The material is seen as a two-phase composite where the inclusion phase represents grain cores, their behavior is driven by dislocation glide, and where the matrix phase, governed by the newly introduced dislocation emission and penetration mechanism, represents both grain boundaries and triple junctions. The long range stress field arising from the presence of grain boundaries is taken into account in the critical glide resistance stress at 0 K in the inclusion phase. The model is applied to polycrystal copper and results in pure tension and creep are compared to experiments. Good agreements between the experimental measurements and the model predictions are observed.     

Study of plastic/viscoplastic models with various inelastic mechanisms

This paper describes a general framework for the development of plastic or viscoplastic constitutive equations. As the applications are focused on cyclic loadings, only small strains are considered, with an additive decomposition of the total strain into a thermo-elastic part, and several inelastic parts, the evolution of which is determined by several plastic or viscoplastic criteria. Quadratic or linear (crystallographic) criteria could be used, so that the approach is able to describe the contribution of several physical levels, or deformation mechanism, to the inelastic behavior. The present work is restricted to the case of quadratic criteria, and specially to the study of the various interactions which can be introduced between the mechanisms. The most important case is the coupling between kinematic hardening variables which allows to describe: (1) either normal rate sensitivity or inverse rate sensitivity; (2) plasticity-creep interaction; (3) ratcheting for high mean stress but either adaptation or plastic shakedown for lower mean stress.     

黏塑性本构计算的稳定性分析

提出本构方程计算方法的稳定性问题,针对黏塑性本构计算的显式精确算法的稳定性进行分析,发现该算法并非无条件稳定,使用小扰动方法给出了其计算稳定的必要条件,稳定性条件对数值计算中的时间步长提出限制要求。通过有限元算例验证了分析的正确性,计算结果也表明理论推导得到的稳定性公式能够准确预测满足计算稳定性条件要求的最大时间步长与各参数之间关系。     

A semi-analytical model for the behavior of saturated viscoplastic materials containing two populations of voids of different sizes

This paper presents a micromechanical model for a porous viscoplastic material containing two populations of pressurized voids of different sizes. Three scales are distinguished: the microscopic scale (corresponding to the size of the small voids), the mesoscopic scale (corresponding to the size of the large voids) and the macroscopic scale. It is assumed that the first homogenization step is performed at the microscopic scale, and, at the mesoscopic scale, the matrix is taken to be homogeneous and compressible. At the mesoscopic scale, the second homogenization step, on which the present study focuses, is based on a simplified representative volume element: a hollow sphere containing a pressurized void surrounded by a nonlinear viscoplastic compressible matrix. The nonlinear behavior of the matrix, which is expressed using the results obtained in the first homogenization step, is approached using a modified secant linearization procedure involving the discretization of the hollow sphere into concentric layers. Each layer has uniform secant moduli. The predictions of the model are compared with the more accurate numerical results obtained using the finite element method. Good agreement is found to exist with all the macroscopic stress triaxialities and all the porosity and nonlinearity values studied.     

A theory of inelastic behavior of materials Part II. Inelastic materials

Archive for Rational Mechanics and Analysis -     

Constitutive model for plastic deformation of nanocrystalline materials with shear band

A phase mixture model was used to study the plastic deformation behaviors in hardening stage of nanocrystalline materials. The material was considered as a composite of grain interior phase and grain boundary (GB) phase. The constitutive equations of the two phases were determined in term of their main deformation mechanisms. In softening stage, a shear band deformation mechanism was presented and the corresponding constitutive relation was established. Numerical simulations have shown that the predications fit well with experimental data. The investigation using the finite-element method (FEM) provided a direct insight into quantifying shear localization effect in nanocrystalline materials.     

Thermomechanical model of elastic–plastic materials with defect structures

A complete thermomechanical model for elastic–plastic materials that account for the interaction of defects in structures is proposed. The characteristics of affine-metric space can be used as thermodynamic variables for constructing the model. Applications to describe the pattern of fracture for underground mines are made.