A theory of inelastic behavior of materials Part I. Ideal inelastic materials

Archive for Rational Mechanics and Analysis -     

Experimental study of the inelastic behavior of materials

This paper presents the results of an experimental study of the elastic and viscoplastic properties of alloys using elastoplastic theory. The need to study complex loading of materials, consider the microstructure evolution during deformation, and develop an appropriate experimental design theory and adequate constitutive relations is highlighted.     

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

Archive for Rational Mechanics and Analysis -     

A generalized self-consistent polycrystal model for the yield strength of nanocrystalline materials

Inspired by recent molecular dynamic simulations of nanocrystalline solids, a generalized self-consistent polycrystal model is proposed to study the transition of yield strength of polycrystalline metals as the grain size decreases from the traditional coarse grain to the nanometer scale. These atomic simulations revealed that a significant portion of atoms resides in the grain boundaries and the plastic flow of the grain-boundary region is responsible for the unique characteristics displayed by such materials. The proposed model takes each oriented grain and its immediate grain boundary to form a pair, which in turn is embedded in the infinite effective medium with a property representing the orientational average of all these pairs. We make use of the linear comparison composite to determine the nonlinear behavior of the nanocrystalline polycrystal through the concept of secant moduli. To this end an auxiliary problem of Christensen and Lo (J. Mech. Phys. Solids 27 (1979) 315) superimposed on the eigenstrain field of Luo and Weng (Mech. Mater. 6 (1987) 347) is first considered, and then the nonlinear elastoplastic polycrystal problem is addressed. The plastic flow of each grain is calculated from its crystallographic slips, but the plastic behavior of the grain-boundary phase is modeled as that of an amorphous material. The calculated yield stress for Cu is found to follow the classic Hall-Petch relation initially, but as the gain size decreases it begins to depart from it. The yield strength eventually attains a maximum at a critical grain size and then the Hall-Petch slope turns negative in the nano-range. It is also found that, when the Hall-Petch relation is observed, the plastic behavior of the polycrystal is governed by crystallographic slips in the grains, but when the slope is negative it is governed by the grain boundaries. During the transition both grains and grain boundaries contribute competitively.     

The scale effect on the yield strength of nanocrystalline materials

The microstructure of the nanocrystalline can be divided generally into two parts: grain and interface. When the grain size is about or less than 10 nm, the interface can be divided into grain boundary and triple junctions. The mechanical performance of nanocrystalline materials with complicated microstructures is greatly different from that of the coarse grain materials. In this paper, the nanocrystalline material is considered as a composite with three phases: the grain core, the grain boundaries, and the triple junction. The model analysis for nanocrystalline material deformation is established and the relationship between yield strength and grain size is obtained. The obtained result explains the inverse Hall–Petch relation.     

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.     

Micromechanics modeling of strength for nanocrystalline copper

We present a model in this paper for predicting the inverse Hall–Petch phenomenon in nanocrystalline (NC) materials which are assumed to consist of two phases: grain phase of spherical or spheroidal shapes and grain boundary phase. The deformation of the grain phase has an elasto-viscoplastic behavior, which includes dislocation glide mechanism, Coble creep and Nabarro–Herring creep. However the deformation of grain boundary phase is assumed to be the mechanism of grain boundary diffusion. A Hill self-consistent method is used to describe the behavior of nanocrystalline pure copper subjected to uniaxial tension. Finally, the effects of grain size and its distribution, grain shape and strain rate on the yield strength and stress–strain curve of the pure copper are investigated. The obtained results are compared with relevant experimental data in the literature.     

A porosity-dependent inelastic criterion for engineering materials

Many criteria have been developed to describe the yielding condition, plastic potential, and failure strength of engineering materials. In this paper, the authors first review the characteristics of a few of the more common criteria used for porous materials. It is then shown that the main features of many criteria can be represented by a unique set of recently developed equations. The ensuing multiaxial criterion becomes applicable to a variety of materials and loading states. One of the advantages of the proposed criterion, named MSDP_u, is that it is explicitly porosity-dependent. The validity of this general inelastic criterion is demonstrated using experimental results obtained from various types of materials. A brief discussion follows on the advantages and limitations of the proposed equations.     

Homogenization of inelastic solid materials at finite strains based on incremental minimization principles. Application to the texture analysis of polycrystals

The paper presents new continuous and discrete variational formulations for the homogenization analysis of inelastic solid materials undergoing finite strains. The point of departure is a general internal variable formulation that determines the inelastic response of the constituents of a typical micro-structure as a generalized standard medium in terms of an energy storage and a dissipation function. Consistent with this type of finite inelasticity we develop a new incremental variational formulation of the local constitutive response, where a quasi-hyperelastic micro-stress potential is obtained from a local minimization problem with respect to the internal variables. It is shown that this local minimization problem determines the internal state of the material for finite increments of time. We specify the local variational formulation for a distinct setting of multi-surface inelasticity and develop a numerical solution technique based on a time discretization of the internal variables. The existence of the quasi-hyperelastic stress potential allows the extension of homogenization approaches of finite elasticity to the incremental setting of finite inelasticity. Focussing on macro-deformation-driven micro-structures, we develop a new incremental variational formulation of the global homogenization problem for generalized standard materials at finite strains, where a quasi-hyperelastic macro-stress potential is obtained from a global minimization problem with respect to the fine-scale displacement fluctuation field. It is shown that this global minimization problem determines the state of the micro-structure for finite increments of time. We consider three different settings of the global variational problem for prescribed displacements, non-trivial periodic displacements and prescribed stresses on the boundary of the micro-structure and develop numerical solution methods based on a spatial discretization of the fine-scale displacement fluctuation field. Representative applications of the proposed minimization principles are demonstrated for a constitutive model of crystal plasticity and the homogenization problem of texture analysis in polycrystalline aggregates.     

Higher-order theory for periodic multiphase materials with inelastic phases

An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress and inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model and with finite-element results for transverse loading.     

A variational multiscale constitutive model for nanocrystalline materials

This paper presents a variational multi-scale constitutive model in the finite deformation regime capable of capturing the mechanical behavior of nanocrystalline (nc) fcc metals. The nc-material is modeled as a two-phase material consisting of a grain interior phase and a grain boundary effected zone (GBAZ). A rate-independent isotropic porous plasticity model is employed to describe the GBAZ, whereas a crystal-plasticity model which accounts for the transition from partial dislocation to full dislocation mediated plasticity is employed for the grain interior. The constitutive models of both phases are formulated in a small strain framework and extended to finite deformation by use of logarithmic and exponential mappings. Assuming the rule of mixtures, the overall behavior of a given grain is obtained via volume averaging. The scale transition from a single grain to a polycrystal is achieved by Taylor-type homogenization where a log-normal grain size distribution is assumed. It is shown that the proposed model is able to capture the inverse Hall-Petch effect, i.e., loss of strength with grain size refinement. Finally, the predictive capability of the model is validated against experimental results on nanocrystalline copper and nickel.     

An embedding method for modeling micromechanical behavior and macroscopic properties of composite materials

This paper presents a numerical method for modeling the micromechanical behavior and macroscopic properties of fiber-reinforced composites and perforated materials. The material is modeled by a finite rectangular domain containing multiple circular holes and elastic inclusions. The rectangular domain is assumed to be embedded within a larger circular domain with fictitious boundary loading represented by truncated Fourier series. The analytical solution for the complementary problem of a circular domain containing holes and inclusions is obtained by using a combination of the series expansion technique with a direct boundary integral method. The boundary conditions on the physical external boundary are satisfied by adopting an overspecification technique based on a least squares approximation. All of the integrals arising in the method can be evaluated analytically. As a result, the elastic fields and effective properties can be expressed explicitly in terms of the coefficients in the series expansions. Several numerical experiments are conducted to verify the accuracy and efficiency of the numerical method and to demonstrate its application in determination of the macroscopic properties of composite materials.     

Statistical continuum theory for inelastic behavior of a two-phase medium

A statistical continuum mechanics formulation is presented to predict the inelastic behavior of a medium consisting of two isotropic phases. The phase distribution and morphology are represented by a two-point probability function. The isotropic behavior of the single phase medium is represented by a power law relationship between the strain rate and the resolved local shear stress. It is assumed that the elastic contribution to deformation is negligible. A Green’s function solution to the equations of stress equilibrium is used to obtain the constitutive law for the heterogeneous medium. This relationship links the local velocity gradient to the macroscopic velocity gradient and local viscoplastic modulus. The statistical continuum theory is introduced into the localization relation to obtain a closed form solution. Using a Taylor series expansion an approximate solution is obtained and compared to the Taylor’s upper-bound for the inelastic effective modulus. The model is applied for the two classical cases of spherical and unidirectional discontinuous fiber-reinforced two-phase media with varying size and orientation.     

The application of ideas associated with materials with memory to modeling the inelastic behavior of solid bodies

In this paper, we present a new way to describe the rate-independent inelastic behavior of metals undergoing finite deformations. Experiments indicate that the stress often has a stronger dependence on the inelastic history in the more recent past as compared to that in the distant past. For this reason, an “annihilation” function is used to weight the inelastic history so that less importance is given to the strains in the more distant past. This “annihilation” function does not depend explicitly on time, but instead on the pathlength associated with the history of stress-free or natural configurations in the strain space relative to the current natural configuration. In this formulation the current configuration is adopted as the reference configuration for the kinematic quantities. The constitutive equation for stress is expressed in terms of the strain associated with the current natural configuration relative to the current actual configuration. Equations have been developed to prescribe the change in the natural configuration as the material yields. A general yield function has been defined in terms of the relative natural strain to restrict the manner in which the natural configuration changes. Since the yield conditions are in terms of the relative natural strain, we can account for situations in which a material yields during the process of unloading. For the sake of simplicity, the elastic properties of the material are considered to be constant throughout the deformation with the material remaining isotropic with respect to the natural configuration. With the proposed theory, we have examined the “Bauschinger effect” which is exhibited by metals that are deformed beyond the yield limit before being deformed beyond the yield limit in the reverse direction.     

层状岩体的Cosserat介质均匀化

将层状岩体视为由Cosserat介质组成的复合体,对其进行均匀化等效,在MATLAB平台上进行有限元编程,得到新等效本构关系。不仅能反映层状岩体的弯曲特性,而且可以对节理岩体的力学性能进行宏观和细观描述。通过对洞室算例进行显式节理模型、横观各向同性模型以及均匀化Cosserat等效模型三种方法的数值计算比较,证明理论分析与计算方法可行有效。     

A material force method for inelastic fracture mechanics

A material force method is proposed for evaluating the energy release rate and work rate of dissipation for fracture in inelastic materials. The inelastic material response is characterized by an internal variable model with an explicitly defined free energy density and dissipation potential. Expressions for the global material and dissipation forces are obtained from a global balance of energy-momentum that incorporates dissipation from inelastic material behavior. It is shown that in the special case of steady-state growth, the global dissipation force equals the work rate of dissipation, and the global material force and J-integral methods are equivalent. For implementation in finite element computations, an equivalent domain expression of the global material force is developed from the weak form of the energy-momentum balance. The method is applied to model problems of cohesive fracture in a remote K-field for viscoelasticity and elastoplasticity. The viscoelastic problem is used to compare various element discretizations in combination with different schemes for computing strain gradients. For the elastoplastic problem, the effects of cohesive and bulk properties on the plastic dissipation are examined using calculations of the global dissipation force.     

Prediction of bending rupture strength of non-linear materials with different behavior in tension and compression

Many materials have quite different stress-strain relations in tension and compression. Examples include such diverse materials as rock, cast iron, concrete, tire cord-rubber and soft biological tissues. It is shown by analysis in this paper that DBTC (different behavior in tension and compression) has a profound effect on the flexural strength as predicted by application of fundamental continuum mechanics relations. The theory is applied to a non-linear material model which is shown to be applicable to two widely different materials: concrete, which has more strength degradation in tension than in compression, and steel cord-rubber, which has strength which is enhanced in tension by cord strengthening and degraded in compression by cord microbuckling.     

Dislocation loop nucleation and plastic deformation of nanocrystalline materials

We propose a three-dimensional model of plastic deformation of a mechanically loaded nanocrystalline material by means of heterogeneous nucleation of loops of complete and partial lattice dislocations and by means of grain boundary dislocations on the already formed dislocation loops. We calculate and compare the energy variations characterizing various versions of dislocation loop nucleation. We discover three basic regions of variation of nanocrystalline material grain dimensions which are characterized by their own types of loop nucleation. We also study the role played by loop nucleation in plastic and superplastic deformations of nanocrystalline materials.