Anisotropic plastic stress fields at a rapidly propagating crack-tip

All the stress components at a rapidly propagating crack-tip in an elastic perfectly-plastic material are the functions of only. Making use of this condition and the equations of steady-state motion, stress-strain relations and Hill anisotropic yield condition, we obtain the general solutions in both the cases of anti-plane and in-plane strain. Applying these two general solutions to propagating Mode III and Mode I cracks, respectively, the anisotropic plastic stress fields at the rapidly propagating tips of Mode III and Mode I cracks are derived.     

Anisotropic plastic stress fields at a slowly propagating crack tip

Under the condition that any perfectly plastic stress components at a crack tip are nothing but the functions of 0 only making use of equilibrium equations. Hill anisotropic yield condition and unloading stress-strain relations, in this paper, we derive the general analytical expressions of anisotropic plastic stress fields at the slowly steady propagating tips of plane and anti-plane strain. Applying these general analytical expressions to the concrete cracks, the analytical expressions of anisotropic plastic stress fields at the-slowly steady propagating tips of Mode I and Mode III cracks are obtained. For the isotropic plastic material, the anisotropic plastic stress fields at a slowly propagating crack tip become the perfectly plastic stress fields.     

A constitutive model for orthotropic elasto-plasticity at large strains

Summary The paper presents a thermodynamically consistent constitutive model for elasto-plastic analysis of orthotropic materials at large strain. The elastic and plastic anisotropies are assumed to be persistent in the material but the anisotropy axes can undergo a rigid rotation due to large plastic deformations. The orthotropic yield function is formulated in terms of the generally nonsymmetric Mandel stress tensor such that its skew-symmetric part is additionally taken into account. Special attention is focused on the convexity of the yield surface resulting in the nine-dimensional stress space. Of particular interest are new convexity conditions which do not appear in the classical theory of anisotropic plasticity. They impose additional constraints on the material constants governing the plastic spin. The role of the plastic spin is further studied in simple shear accompanied by large elastic and large plastic deformations. If the plastic spin is neglected, the shear stress response is characterized by oscillations with an amplitude strictly dependent on the degree of the plastic anisotropy.accepted for publication 2 March 2004     

A finite-deformation theory of plasticity

Based on a multiplicative decomposition of local deformation into elastic and plastic deformations general constitutive equations of elastic-plastic materials are proposed. Two alternative approaches are discussed: one in which the elastic deformation is used as an independent variable, and the other in which the stress is one of the independent variables. The appropriate material symmetries are defined, and it is shown that the plastic spin is absent in the theory of isotropic materials. Analysis of a simple extension is given as an example.     

An analysis of thermally-induced plane waves in elastic-plastic single crystals

A framework for the calculation of thermally-induced plane waves in elastic-plastic single crystals of arbitrary crystallographic symmetry and orientation is presented. Plasticity is described in terms of small strain theory and the available slip-planes which can be arbitrary in number as well as in orientation. The effects of perfect-plasticity modify not only the anisotropic elastic moduli, but also the components of the Grüneisen tensor. The latter effect is a consequence of a non-spherical stress state developed in anisotropic materials during rapid energy-absorption at constant strain. Specific examples of thermally-induced plane waves are presented for both the elastic and plastic response of beryllium and graphite single-crystals.     

Finite-element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is ideally suited to isotropically hardening Prandtl-Reuss materials. Further, the formulation is given in a manner which allows any conventional finite element program, for “small strain” elastic-plastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension.The paper closes with a unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures. Further, a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain, and the inadequacies of some of these are commented upon.     

Anisotropic plastic fields at a rapidly propagating plane-stress crack-tip

Under the condition that all the stress components at a crack-tip are the functions ofθonly,making use of the equations of steady-state motion.Hill anisotropic yield condition and stress-strain relations,we obtain the general solution of anisotropic plastic field at a rapidly propagating plane-stress crack-tip.Applying this general solution to four particular cases of anisctropy,the general solutions of these four particular cases are derived.Finally,we give the anisotropic plastic field at the rapidly propagating plane-stress modeⅠcrack-tip in the case of X=Y=Z     

Elastic-plastic behavior of elastically homogeneous materials with a random field of inclusions

A multiphase material is considered, which consists of a homogeneous elastic-plastic matrix containing a homogeneous statistically uniform random set of ellipsoidal elastic-plastic inclusions. The elastic properties of the matrix and the inclusions are the same, but the so-called “stress-free strains”, i.e. the strain contributions due to temperature loading, phase transformations, and the plastic strains, fluctuate. A general theory of the yielding for arbitrary loading (by the stress and by the temperature) is employed. The realization of an incremental plasticity scheme is based on averaging over each component of the nonlinear yield criterion. Usually, averaged stresses are used inside each component for this purpose. In distinction from this usual practice physically consistent assumptions about the dependence of these functions on the component's values of the second stress moments are applied. The application of the proposed theory to the prediction of the thermomechanical deformation behavior of a model material is shown.     

Theoretical study of acoustoelastic effects caused by plastic anisotropy growth

The general expressions for strain-velocity relations, which are called acoustoelastic effects, of waves propagating in elasto-plastically deformed solids have been derived with aid of Jaumann co-rotational rate. In the previous paper (Kobayashi [1986]), for the slightly orthotropic case, the theoretical results of those effects were presented and it was shown that plastic anisotropy influences the acoustoelastic effects through anisotropic stress increment terms constituting plastic strain increment. In the present paper, those anisotropic stress increment terms are defined by using the unsymmetric yield function proposed by Shriastava, Mróz & Dubey [1973]; the acoustoelastic effects caused by the plastic anistropy growth are analyzed.     

The theory of elastic-plastic deformation at finite strain with induced anisotropy modeled as combined isotropic-kinematic hardening

A Theory of elastic-plastic deformation with strain induced anisotropy based on finite-deformation-valid continuum mechanics is presented. On the foundation of nonlinear kinematics which provides strict uncoupling of elastic and plastic deformation rate terms according to their physical origins, it introduces a basis for the modified plastic rate of deformation D̂^p suggested by G.J. Creus, A.G. Groehs and E.T. Onat in a report entitled “Constitutive Equations for Finite Deformations of Elastic-Plastic Solids,” 1984, in which this variable was suggested in order to give an elegant mathematical structure to the theory. D̂^p is shown to express the resultant rate of deformation in the current configuration of the elastically-plastically deforming material which is envisaged to be generated by the pure plastic flow and the anisotropy-caused spin, both considered to be occurring in the unstressed state.From this basis an elastic-plastic theory is developed in the case where the strain-induced anisotropy takes the form of combined isotropic-kinematic hardening, although the concepts involved also apply to more general anisotropic characteristics. A general evolution equation is adopted for the back stress α, the kinematic-hardening shift of the yield surface, its rate of growth being expressed as a general form-invariant function of α and D̂^p, including a general term expressing the influence of the spin of α because it is embedded in the deforming material.By providing an expression for the total strain rate as the sum of the strain rate D̂^p and an elastic term, linear in the Jaumann derivative of Kirchhoff stress, it is shown that (D̂^p dt) is the residual strain increment following a loading/unloading cycle imposed by a stress increment. By considering materials which obey the normality rule it is also shown that the instantaneous elastic-plastic moduli have the symmetries necessary for generating a rate potential function and hence can be incorporated into Hill's variational principle valid for solving problems involving finite deformation and convenient for finite-element exploitation.     

基于微纳米压入法提取无铅焊料合金弹塑性力学参数的研究

本文通过微纳米压入法结合数值模拟研究了无铅焊料合金SnAg3.5 的弹塑性力学性能，分别采用圆柱形压头及两种不同锥角压头对无铅焊料合金进行压入测试：基于圆柱形压头测试过程中接触面积恒定的特点得到了无铅焊料的弹性模量，进一步采用塑性应变梯度理论对两种锥角压头的测试结果予以修正并通过数值模拟反分析得到相应的特征应力值，同时基于压入特征塑性应变与压头锥角的关系式得到两种不同锥角压头下的特征应变值，在此基础上经求解方程组得到焊料合金的初始屈服应力与应变强化因子，进而得到了焊料合金的幂强化弹塑性本构关系．该方法剔除了压入尺度效应的影响并保证了所得本构关系的唯一性，给出了一种通过原位压入测试表征金属材料弹塑性力学性能的有效方法．     

A framework for numerical integration of crystal elasto-plastic constitutive equations compatible with explicit finite element codes

In this paper we summarize the elements of a numerical integration scheme for elasto-plastic response of single crystals. This is intended to be compatible with large-scale explicit finite element codes and therefore can be used for problems involving multiple crystals and also overall behavior of polycrystalline materials. The steps described here are general for anisotropic elastic and plastic response of crystals. The crystallographic axes of the lattice are explicitly stored and updated at each time step. A plastic predictor–elastic corrector scheme is used to calculate the plastic strain rates on all active slip systems based on a rate-dependent physics-based constitutive model without the need of further auxiliary assumptions. Finally we present the results of numerous calculations using a physics-based rate- and temperature-dependent model of copper and the effect of elastic unloading, elastic crystal anisotropy, and deformation-induced lattice rotation are emphasized.     

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.     

A generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials

In this paper, a generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials is derived. The evolution equation for the active yield surface with reference to the memory yield surface is obtained by considering the continuous expansion of the active yield surface during the unloading/reloading process. The incremental constitutive relation based on the associated flow rule is then derived for a general yield function for pressure insensitive and sensitive materials. Detailed incremental constitutive relations for materials based on the Mises yield function, the Hill quadratic anisotropic yield function and the Drucker–Prager yield function are derived as the special cases. The closed-form solutions for one-dimensional stress–plastic strain curves are also derived and plotted for materials under cyclic loading conditions based on the three yield functions. In addition, the closed-form solutions for one-dimensional stress–plastic strain curves for materials based on the isotropic Cazacu–Barlat yield function under cyclic loading conditions are summarized and presented. For materials based on the Mises and the Hill anisotropic yield functions, the stress–plastic strain curves show closed hysteresis loops under uniaxial cyclic loading conditions and the Masing hypothesis is applicable. For materials based on the Drucker–Prager and Cazacu–Barlat yield functions, the stress–plastic strain curves do not close and show the ratcheting effect under uniaxial cyclic loading conditions. The ratcheting effect is due to different strain ranges for a given stress range for the unloading and reloading processes. With these closed-form solutions, the important effects of the yield surface geometry on the cyclic plastic behavior due to the pressure-sensitive yielding or the unsymmetric behavior in tension and compression can be shown unambiguously. The closed form solutions for the Drucker–Prager and Cazacu–Barlat yield functions with the associated flow rule also suggest that a more general anisotropic hardening theory needs to be developed to address the ratcheting effects for a given stress range.     

VISCOPLASTIC SOLUTION TO FIELD AT STEADILY PROPAGATING CRACK TIP IN LINEAR- HARDENING MATERIALS

An elastic-viscoplastic constitutive model was adopted to analyze asymptotically the tip-field of moving crack in linear-hardening materials under plane strain condition. Under the assumption that the artificial viscosity coefficient was in inverse proportion to power law of the rate of effective plastic strain, it is obtained that stress and strain both possess power law singularity and the singularity exponent is uniquely determined by the power law exponent of the rate of effective plastic strain. Variations of zoning structure according to each material parameter were discussed by means of numerical computation for the tip-field of mode Ⅱ dynamic propagating crack, which show that the structure of crack tip field is dominated by hardening coefficient rather than viscosity coefficient. The secondary plastic zone can be ignored for weak hardening materials while the secondary plastic zone and the secondary elastic zone both have important influence on crack tip field for strong hardening materials. The dynamic solution approaches to the corresponding quasi-static solution when the crack moving speed goes to zero, and further approaches to the HR (Hui-Riedel) solution when the hardening coefficient is equal to zero.     

单晶体塑性滑移有限变形下的应力计算

为了探讨和发展单晶金属材料的非弹性有限变形分析方法,提出一种单晶体各向异性弹塑性分析的计算格式.该方法是一种以初始构形为变形计算参考构形的描述方法,它对单晶体塑性构形的演化用增量计算以跟随加载路径,而在应力计算时在卸载构形的基础上用Hencky对数弹性应变来计算总量的应力以保证计算的稳定和收敛;通过求解满足瞬时屈服条件和应力与弹性应变关系的广义胡克定律的非线性方程组来搜索激活滑移系.     

THE PRECISE AND NEW ANALYSIS FOR A MODE Ⅲ GROWING CRACK IN AN ELASTIC-PERFECTLY PLASTIC SOLID

The near crack line field analysis method has been used to investigate into ModeⅢ quasistatically propagating crack in an elastic-perfectly plastic material.Thesignificance of this paper is that the usual small scale yielding theory has been brokenthrough.By obtaining the general solutions of the stresses and the displacement rate ofthe near crack line plastic region,and by matching the general solutions with theprecise elastic fields(not the usual elastic K-dominant fields)at the elastic-plasticboundary,the precise and new solutions of the stress and deformation fields,the sizeof the plastic region and the unit normal vector of the elastic-plastic boundary havebeen obtained near the crack line.The solutions of this paper are sufficiently precisenear the crack line region because the roughly qualitative assumptions of the smallscale yielding theory have not been used and no other roughly qualitative assumptionshave been taken,either.The analysis of this paper shows that the assumingly“steady-state cas     

I型定常扩展裂纹尖端的弹黏塑性场

考虑材料在扩展裂纹尖端的黏性效应，假设黏性系数与塑性应变率的幂次成反比，对幂硬化材料中平面应变扩展裂纹尖端场进行了弹黏塑性渐近分析，得到了不含间断的连续解，并讨论了I型裂纹数值解的性质随各参数的变化规律. 分析表明应力和应变均具有幂奇异性，并且只有在线性硬化时，尖端场的弹、黏、塑性才可以合理匹配. 对于I型裂纹，裂尖场不含弹性卸载区. 当裂纹扩展速度趋于零时，动态解趋于准静态解，表明准静态解是动态解的特殊形式；如果进一步考虑硬化系数为零的极限情况，便可退化为Hui和Riedel的非线性黏弹性解.