A finite strain constitutive model for metal powder compaction using a unique and convex single surface yield function

Constitutive modelling of metal powder compaction processes is a challenge in view of realistic simulations. To this end, the article under consideration has two objectives: the first goal is to present a new unique and convex single surface yield function for pressure dependent materials, which is also applicable to other areas of granular materials such as soils or concrete. The flexibility is shown at various materials. The yield function is based on a log-interpolation of two known simple yield functions. A convexity proof of the new yield function is provided. The second objective is to propose a new rate-independent finite strain plasticity model for metal powder compaction, which is based on the multiplicative decomposition of the deformation gradient into an elastic and a plastic part with evolution equations for internal variables representing the basic behaviour of powder materials under compaction conditions. These variables are used for the evolution of the yield function in order to represent the compressible hardening behaviour of powder materials. On the basis of the constitutive model, the material parameters are identified at experimental data of copper powder.     

An elasto-plastic damage constitutive theory and its prediction of evolution of subsequent yield surfaces and elastic constants

Based on pair functional potentials, Cauchy-Born rule and slip mechanism, a material model assembling with spring-bundle components, a cubage component and slip components is established to describe the elasto-plastic damage constitutive relation under finite deformation. The expansion/shrink, translation and distortion of yield surfaces can be calculated based on the hardening rule and Bauschinger effect defined on the slip component level. Both kinematic and isotropic hardening are included. Numerical simulations and predictions under tension, torsion, and combined tension-torsion proportional/non-proportional loading are performed to obtain the evolution of subsequent yield surfaces and elastic constants and compare with two sets of experimental data in literature, one for a very low work hardening aluminum alloy Al 6061-T6511, and another for a very high work hardening aluminum alloy annealed 1100 Al. The feature of the yield surface in shape change, which presents a sharp front accompanied by a blunt rear under proportional loading, is described by the latent hardening and Bauschinger effect of slip components. Further, the evolution law of subsequent yield surfaces under different proportional loading paths is investigated in terms of their equivalence. The numerical simulations under non-proportional loading conditions for annealed 1100 Al are performed, and the subsequent yield surfaces exhibit mixed cross effect because the kinematic hardening and isotropic hardening follow different evolution tendency when loading path changes. The results of non-proportional loading demonstrate that the present model has the ability to address the issue of complex loading due to the introduction of state variables on slip components. Moreover, as an elasto-plastic damage constitutive model, the present model can also reflect the variation of elastic constants through damage defined on the spring-bundle components.     

花岗岩破坏过程能量演化机制与能量屈服准则

为了明确岩石破坏的能量演化特性，结合单轴实验和颗粒流程序获得花岗岩的细观力学参数，进行不同应力状态的花岗岩实验，研究不同围压下花岗岩破坏过程的能量演化机理并推导能量屈服准则。获得以下主要结论:花岗岩破坏过程中低围压下内部损伤出现较早而高围压较晚，表明低围压花岗岩内部损伤是渐进发展过程，而高围压下内部损伤一旦出现便快速发展破坏；高围压花岗岩峰值前一定应变范围弹性应变能基本保持不变，吸收的能量全部转化为耗散能，表明高围压破坏时花岗岩内部损伤程度严重；弹性应变能经历不断积累并达到弹性储能极限而后减小的变化过程，而弹性储能极限与围压之间存在线性变化规律，因此高围压下岩体开挖卸荷时极易诱发大量弹性应变能的急剧释放，引起围岩失稳甚至发生岩爆；花岗岩峰值破坏时的能量比与围压无关，为一定值；基于能量原理导出了能量屈服准则，该准则包含岩性参数和所有主应力，能够综合反映岩石破坏影响因素。     

A weakened form of Ilyushin's postulate and the structure of self-consistent Eulerian finite elastoplasticity

From a general standpoint in terms of internal variables, we formulate a general theory of self-consistent Eulerian finite elastoplasticity based on the additive decomposition of the Eulerian strain rate, i.e., D=D^e+D^p, as well as two consistency criteria. In this theory, the elastic behaviour is characterized by an exactly integrable elastic rate equation for D^e with a general form of complementary elastic potential. It is assumed that the yield function depends in a general manner on the Kirchhoff stress and the internal variables. Moreover, the plastic rate equation for D^p and the evolution equation for each internal variable are allowed to assume general forms relying on the just-mentioned variables and the stress rate. It is indicated that two consistency criteria, i.e., the self-consistency for the elastic rate equation and Prager's yielding stationarity, lead to the unique choice of objective rates, i.e., the logarithmic rate.The structure of the above theory is further studied and examined by virtue of a weakened form of Ilyushin's postulate. In a spinning frame defining the logarithmic rate, we introduce the notion of standard elastoplastic strain cycle, which starts at a point not on but inside a yield surface and incorporates only one infinitesimal plastic subpath. We show that this type of strain cycle is always possible. Then, by ruling out strain cycles starting at points on yield surfaces we propose a weakened form of Ilyushin's postulate, which says that the changing rate of the stress work done along every standard strain cycle should be non-negative, whenever the incorporated plastic subpath tends to vanish. By virtue of simple, rigorous procedures, we demonstrate that this weakened form of Ilyushin's postulate is adequate to ensure direct results concerning the normality rule and the convexity of the yield surface in the context of the foregoing Eulerian finite elastoplasticity theory. Specifically, with an exactly integrable elastic rate equation defining D^e, we prove that, in the space of the Kirchhoff stresses, the difference (D−D^e) is just the gradient of the yield function multiplied by a plastic multiplier, and thus bears the very kinematical and physical feature of plastic strain rate. Furthermore, we prove that, in the space of the Kirchhoff stresses, the elastic domain bounded by each yield surface should be convex. The main results are derived in a self-contained manner within the context of an Eulerian theory of finite elastoplasticity, without involving issues concerning how to define intermediate stress-free states and plastic strains, etc.     

Prediction of forming limit curves using an anisotropic yield function with prestrain induced backstress

The Forming Limit Diagram (FLD), a plot of the maximum major principal strains that can be sustained by sheet materials prior to the onset of localized necking, is a useful concept for characterizing the formability of sheet metal. Both experimental and numerical results in the literature have shown that the level of the FLD is strongly strain path dependent and the prediction of FLD depends on the shape of the initial yield function and its evolution. In this work, to improve the accuracy of FLD prediction under nonlinear strain paths for a given material, the evolution of the yield function is proposed in terms of the changes of its center and its curvature. The center of the subsequent yield surface after preloading and unloading will be determined via a backstress tensor, and the curvature change will be reflected by changing the exponent in the yield function. Both parameters are functions of the effective plastic strain and will be determined using the forming limit strains obtained from two nonlinear tests. Using this approach, a combination of Marciniak–Kuczynski (M–K) analysis (Marciniak, Z., Kuczynski, K. 1967. Limit strains in the processes of stretch-forming sheet metal. Int. J. Mech. Sci. 9, 609.) and a general anisotropic yield criterion developed by Karafillis and Boyce (Karafillis, A.P., Boyce, M.C. 1993. A general anisotropic yield criterion using bounds and transformation weighting tensor, J. Mech. Phys. Solids, 41, 1859) is used to predict nonlinear FLDs of both Al2008-T4 and Al6111-T4. Excellent agreements were obtained between computed FLDs with the experimental data of Graf and Hosford (Graf, A., Hosford, W.F. 1993a. Calculations of forming limit diagrams for changing strain paths. Metall. Trans. A. 24, 2497; Graf, A., Hosford, W.F. 1993b. Effect of changing strain paths on forming limit diagrams of Al 2008-T4. Metall. Trans. A. 24, 2503; Graf, A., Hosford, W.F. 1994. The influence of strain path changes on forming limit diagrams of Al 6111-T4. Int. J. Mech. Sci. 36, 897). This prediction capability provides a powerful tool in the design and optimization process of 3D sheet metal forming where strain path changes are inevitable.     

A FE formulation for elasto-plastic materials with planar anisotropic yield functions and plastic spin

In this article a stress integration algorithm for shell problems with planar anisotropic yield functions is derived. The evolution of the anisotropy directions is determined on the basis of the plastic and material spin. It is assumed that the strains inducing the anisotropy of the pre-existing preferred orientation are much larger than subsequent strains due to further deformations. The change of the locally preferred orientations to each other during further deformations is considered to be neglectable. Sheet forming processes are typical applications for such material assumptions. Thus the shape of the yield function remains unchanged. The size of the yield locus and its orientation is described with isotropic hardening and plastic and material spin.The numerical treatment is derived from the multiplicative decomposition of the deformation gradient and thermodynamic considerations in the intermediate configuration. A common formulation of the plastic spin completes the governing equations in the intermediate configuration. These equations are then pushed forward into the current configuration and the elastic deformation is restricted to small strains to obtain a simple set of constitutive equations. Based on these equations the algorithmic treatment is derived for planar anisotropic shell formulations incorporating large rotations and finite strains. The numerical approach is completed by generalizing the Return Mapping algorithm to problems with plastic spin applying Hill’s anisotropic yield function. Results of numerical simulations are presented to assess the proposed approach and the significance of the plastic spin in the deformation process.     

Analysis of circular hole expansion with generalized yield criteria

Expansion of a circular hole, embedded in an infinite elastoplastic sheet, is studied within the framework of large strain plane-stress plasticity. Material response is modeled by deformation type theories with two families of generalized isotropic yield criteria. Two distinct problems are examined in detail: hole enlargement under internal pressure and hole expansion under remote tension. Strain hardening and elastic compressibility are fully accounted for.Numerical illustrations reveal constitutive sensitivity of stress and deformation profiles. For the internally pressurized hole the specific power needed to create a new volume unit reaches an asymptotic level practically independent of yield criteria. The specific cavitation power is used to derive a simple relation for the ballistic limit in quasi-static plate perforation, showing good agreement with experimental results. Under remote tension the hole expands spontaneously when external stress approaches a limit which is found to be highly sensitive to the yield criteria.     

Mesoscopic natural element analysis of elastic moduli, yield stress and fracture of solids containing a number of voids

A two-dimensional mesoscale simulation method based on the natural element method, which is a kind of meshless method, is developed and applied to the analysis of overall elastic moduli, macro yield stress and void-linking fracture. The calculated results are compared with the theoretical solutions for overall elastic moduli, the experimental macro yield stress for aluminum and brass as well as the improved Gurson’s yield function, and the experimental void linking fracture to show the validity of the proposed method.     

Theoretical considerations upon the MK model for limit strains prediction: The plane strain case,strain-rate effects,yield surface influence,and material heterogeneity

The paper presents a study of the Marciniak and Kuczynski (MK for short) model for the prediction of limit strains of orthotropic sheet metal under in-plane proportional biaxial stretching. In two particular cases analytical results can be obtained if the groove of the MK model is oriented along one of the in-plane symmetry axes. The first case is the plane strain loading mode. Necessary and sufficient conditions are derived for the MK-predicted plane strain limit strain to match exactly the experimentally measured limit strain. An example of material, the AA5182-O aluminum alloy, that does not satisfy these conditions is discussed. It is shown then that if a power-law strain rate sensitivity is included in the hardening law then the MK-model can match exactly any target plane strain limit strain. The second case is the non-hardening case for positive strain ratios. This case allows for an insight into the way the MK-predicted limit strains depend upon the yield function. Based on the theory developed for the plane strain case, material heterogeneity as a possible cause for unstable plastic flow is further discussed. It is shown that such heterogeneities can be modeled by perturbing the rate of deformation with an eigenstrain. This allows for an extension of the MK-model to sheets of uniform thickness.     

Dynamic Brittle Fracture as a Small Horizon Limit of Peridynamics

We consider the nonlocal formulation of continuum mechanics described by peridynamics. We provide a link between peridynamic evolution and brittle fracture evolution for a broad class of peridynamic potentials associated with unstable peridynamic constitutive laws. Distinguished limits of peridynamic evolutions are identified that correspond to vanishing peridynamic horizon. The limit evolution has both bounded linear elastic energy and Griffith surface energy. The limit evolution corresponds to the simultaneous evolution of elastic displacement and fracture. For points in spacetime not on the crack set the displacement field evolves according to the linear elastic wave equation. The wave equation provides the dynamic coupling between elastic waves and the evolving fracture path inside the media. The elastic moduli, wave speed and energy release rate for the evolution are explicitly determined by moments of the peridynamic influence function and the peridynamic potential energy.     

Plastic behavior of some yield stress fluids: from creep to long-time yield

The yielding of several reversible yield stress fluids is studied during scissometric-like creep experiments. The temporal evolution of the apparent deformation is recorded for applied stresses close and below the usual yield stress. Similarly to solids, three main creep regimes are observed. First, a primary creep regime displaying a temporal power law evolution of the deformation rate occurs, followed by a temporal minimum, which leads to an apparent flow of the material. This local minimum, defined as the “transition time,” and the subsequent fluidization can be observed at long times. The evolution of this time as a function of the applied stress appears to follow a universal law reminiscent of fracture behavior in hard solids.     

圆管截面齐次广义屈服函数与结构极限承载力

为克服圆管截面广义屈服准则不满足比例加载条件,导致采用弹性模量调整法求解该类结构极限承载力时存在计算结果受荷载初值影响、计算精度受损等问题,利用回归分析和最小二乘法研究建立了圆管截面广义屈服函数的齐次多项式,通过误差分析确定了齐次化多项式的阶次;据此定义了圆管截面薄壁构件的单元承载比、承载比均匀度和基准承载比,为高承载比薄壁单元的判别及其弹性模量调整提供了动态判据,进而依据能量守恒准则建立了以单元承载比为基本参数的模量调整公式,结合下限原理提出了圆管截面薄壁结构极限承载力分析的弹性模量缩减法。研究表明,选取齐次化多项式的广义屈服函数能更加准确地考虑各项内力对结构极限承载力的综合影响,具有良好的计算精度和效率,可应用于复杂圆管截面薄壁结构的极限承载力分析中。     

Where is the elastic limit?

The elastic limit provides a convenient concept for the design of mechanical and structural parts which should exhibit no permanent deformation after loading and subsequent unloading. The offset tensile yield strength of a material for small offsets, such as the 0.01-percent offset, is considered to be a good approximation of the elastic limit. TheWT-bend tester provides an alternative method to the tension test of determining the offset yield strength of materials. The specimens are subjected to cyclic bending and energy dissipation is used as a yield criterion. The stress as a function of the offset has been determined for a number of alloys. For some of the materials investigated cyclic stresses at levels considerably below the 0.01-percent offset yield strength caused significant changes in mechanical properties. Furthermore, for some highly cold-worked materials substantial cyclic softening could be observed. This raises the question: where is the elastic limit? It is hypothesized that no true elastic limit exists and that it would be possible only to determine an anelastic limit.     

Prediction of localized thinning in sheet metal using a general anisotropic yield criterion

The forming limit diagram (FLD) is a useful concept for characterizing the formability of sheet metal. The ability to accurately predict the FLD for a given material has been shown to depend on the shape of the selected yield function. In addition, both experimental and numerical results have shown that the level of the FLD is strongly strain path dependent. In this work, a combination of Marciniak–Kuczynski (M–K) analysis and a general anisotropic yield criterion developed by Karafillis and Boyce (Karafillis, A.P., Boyce, M.C., 1993. A general anisotropic yield criterion using bounds and transformation weighting tensor. J. Mech. Phys. Solids 41, 1859) is used to predict localized thinning of sheet metal alloys for linear and nonlinear strain paths. A new method for determining the constants in the yield criterion is proposed. The optimal values are obtained by fitting the initial yield stresses and calculated FLD under linear strain paths with the experimental measurement. Using this approach, accurate yield functions can be defined for both Al2008-T4 and Al6111-T4. Comparisons of computed FLDs with the experimental data of Graf and Hosford (Graf, A., Hosford, W.F., 1993b. Effect of changing strain paths on forming limit diagrams of Al 2008-T4. Metall. Trans. A. 24, 2503; Graf, A., Hosford, W.F., 1994. The influence of strain path changes on forming limit diagrams of Al 6111-T4. Int. J. Mech. Sci. 36, 897) show good agreements.     

Statistical aspects of the identification of material parameters for elasto-plastic models

Summary The presented method to identify material parameters for inelastic deformation laws is based on the numerical analysis of inhomogeneous stress and strain fields received from suitable experiments. Tensile and bending tests were carried out to obtain elastic and hardening parameters. The deformation law for small elasto-plastic strains is presented as a system of nonlinear differential and algebraic equations (DAE) consisting of the stress–strain relation, evolution equations for the internal variables and the yield condition. Different rules for the evolution equations of isotropic, kinematic and distorsional hardening are proposed. The DAE are discretized using an implicit Euler method, and the resulting system of nonlinear algebraic equations is solved using the Newton method. Deterministic optimization procedures are preferred to identify material parameters from a least-squares functional of numerical and measured comparative quantities. The gradient of the objective function was calculated using a semianalytical sensitivity analysis. Due to measurement errors, the optimal sets of material parameters are non unique. The approximate estimation of confidence regions and the calculation of correlation coefficients is presented. The results of several optimization processes for material parameters of elasto-plastic deformation laws show a good agreement between measured and calculated values, but they show also problems which may occur if systematic errors will not be recognized and deleted. Received 30 September 1999; accepted for publication 8 March 2000     

A coupled temperature and strain rate dependent yield function for dynamic deformations of bcc metals

A coupled temperature and strain rate microstructure physically based yield function is proposed in this work. It is incorporated along with the Clausius–Duhem inequality and an appropriate free energy definition in a general thermodynamic framework for deriving a three-dimensional kinematical model for thermo-viscoplastic deformations of body centered cubic (bcc) metals. The evolution equations are expressed in terms of the material time derivatives of the elastic strain, accumulated plastic strain (isotropic hardening), and the back stress conjugate tensor (kinematic hardening). The viscoplastic multipliers are obtained using both the Consistency and Perzyna viscoplasticity models. The athermal yield function is employed instead of the static yield function in the case of the Perzyna viscoplasticity model. It is found that the static strain rate value, at which the material shows rate-independent behavior, varies with the material deformation temperature. Computational aspects of the proposed model are addressed through the finite element implementation with an implicit stress integration algorithm. Finite element simulations are performed by implementing the proposed viscoplasticity constitutive models in the commercial finite element program ABAQUS/Explicit [ABAQUS, 2003. User Manual, Version 6.3. Habbitt, Karlsson and Sorensen Inc., Providence, RI] via the user material subroutine coded as VUMAT. Numerical implementation for a simple compression problem meshed with one element is used to validate the proposed model implementation with applications to tantalum, niobium, and vanadium at low and high strain rates and temperatures. The analysis of a tensile shear banding is also investigated to show the effectiveness and the performance of the proposed framework in describing the strain localizations at high velocity impact. Results show mesh independency as a result of the viscoplastic regularization used in the proposed formulation.     

Orientation tensors in simple flows of dilute suspensions of non-Brownian rigid ellipsoids,comparison of analytical and approximate solutions

General analytical solutions are obtained for the planar orientation structure of rigid ellipsoid of revolutions subjected to an arbitrary homogeneous flow in a Newtonian fluid. Both finite and infinite aspect ratio particles are considered. The orientation structure is described in terms of two-dimensional, time-dependent tensors that are commonly employed in constitutive equations for anisotropic fluids such as fiber suspensions. The effect of particle aspect ratio on the evolution of orientation structure is studied in simple shear and planar elongational flows. With the availability of analytical solutions, accuracies of quadratic closure approximations used for nonhomogeneous flows are analyzed, avoiding numerical integration of orientation distribution function. In general, fourth-order orientation evolution equations with sixth-order quadratic closure approximations yield more accurate representations compared to the commonly used second-order evolution equations with fourth-order quadratic closure approximations. However, quadratic closure approximations of any order are found to give correct maximum orientation angle (i.e., preferred direction) results for all particle aspect ratios and flow cases.     

Complete Integrability of Shock Clustering and Burgers Turbulence

We consider scalar conservation laws with convex flux and random initial data. The Hopf–Lax formula induces a deterministic evolution of the law of the initial data. In a recent article, we derived a kinetic theory and Lax equations to describe the evolution of the law under the assumption that the initial datum is a spectrally negative Markov process. Here we show that: (i) the Lax equations are Hamiltonian and describe a principle of least action on the Markov group; (ii) the Lax equations are completely integrable and linearized via a loop-group factorization of operators; (iii) the associated zero-curvature equations can be solved via inverse scattering. Our results are rigorous for N-dimensional approximations of the Lax equations, and yield formulas for the limit N → ∞. The main observation is that the Lax equations and zero-curvature equations are a Markovian analog of known integrable systems (geodesic flow on Lie groups and the N-wave model respectively). This allows us to introduce a variety of methods from the theory of integrable systems.