Macroscopic yield criteria for plastic anisotropic materials containing spheroidal voids

The combined effects of void shape and matrix anisotropy on the macroscopic response of ductile porous solids is investigated. The Gologanu–Leblond–Devaux’s (GLD) analysis of an rigid-ideal plastic (von Mises) spheroidal volume containing a confocal spheroidal cavity loaded axisymmetrically is extended to the case when the matrix is anisotropic (obeying Hill’s [Hill, R., 1948. A theory of yielding and plastic flow of anisotropic solids. Proc. Roy. Soc. London A 193, 281–297] anisotropic yield criterion) and the representative volume element is subjected to arbitrary deformation. To derive the overall anisotropic yield criterion, a limit analysis approach is used. Conditions of homogeneous boundary strain rate are imposed on every ellipsoidal confocal with the cavity. A two-field trial velocity satisfying these boundary conditions are considered. It is shown that for cylindrical and spherical void geometries, the proposed criterion reduces to existing anisotropic Gurson-like yield criteria. Furthermore, it is shown that for the case when the matrix is considered isotropic, the new results provide a rigorous generalization to the GLD model. Finally, the accuracy of the proposed approximate yield criterion for plastic anisotropic media containing non-spherical voids is assessed through comparison with numerical results.     

Approximate yield criteria for anisotropic metals with prolate or oblate voidsCritères macroscopiques pour des métaux plastiques anisotropes contenant des cavités non sphériques

Following the study of Gologanu et al. (1997) which has extended the well-known approach of Gurson (1975), we propose approximate yield criteria for anisotropic plastic voided metals containing non spherical cavities. The plastic anisotropy of the matrix is described by means of Hill's quadratic criterion. The procedure to establish the closed form expression of approximate macroscopic criteria, in which void shape and plastic anisotropic effects are included, is detailed. The new criteria allow us to recover existing results in the cases of spherical and cylindrical voids in an Hill type plastic matrix. Moreover, they agree with previous criteria for non spherical voids in an isotropic plastic matrix. Finally, for validation purposes, we provide, in the general case of non spherical cavities in the anisotropic matrix, a comparison with the numerical exact two field criteria. To cite this article: V. Monchiet et al., C. R. Mecanique 334 (2006).     

基于非关联流动规律的Gotoh屈服准则的参数确定方法

屈服准则对板料成形过程的理论解析、工艺优化和有限元模拟有着重要的影响. 通过提高屈服准则的各向异性表征能力, 可以确保成形过程的可靠性及实际预测的准确性. 本文基于非关联流动法则, 给出了Gotoh屈服准则一套全新的参数求解方法. 在结合常用屈服准则并考虑流动规律的基础上, 分别以5754O铝合金、DP980先进高强钢和SAPH440结构钢作为研究对象, 进行了不同加载路径下各向异性变形行为的预测. 根据Gotoh屈服准则推导的屈服函数、塑性势函数以及基于关联流动的理论函数计算出屈服应力和各向异性指数$r$值随加载角度的分布趋势, 进而针对平面应力状态的屈服轨迹展开分析, 验证了不同屈服准则和流动规律对各向异性屈服行为的预测精度. 理论与实验数据对比结果表明: 不同屈服准则针对同种板料在流动规律一致的情形下其表征各向异性的能力有显著差异; 相同屈服准则基于不同流动规律其表征能力也具有明显差别. 基于非关联流动的屈服准则能极大地提高精度, 各向异性表征能力显著加强. 相关结果能够为各向异性屈服准则在塑性成形领域的实际应用方案提供重要参考.      

Interpretation of the behavior of metals under large plastic shear deformations: comparison of macroscopic predictions to physically based predictions

A large plastic shear problem is analyzed by application of a macroscopic anisotropic plasticity model (Kuroda, M., 1997. Interpretation of the behavior of metals under large plastic shear deformations: a macroscopic approach. Int. J. Plasticity 13, 359–383), and the results are compared to predictions based on crystal plasticity with the Taylor assumption. It is found that these two different-scale models provide very similar predictions. The interpretations for such similarities are pursued in detail. The present macroscopic model reproduces quite well the change in orientation of anisotropy, which is directly predicted in the crystal plasticity analyses as a macroscopic manifestation of texture development. Consequently, the predictions for the rotation of the yield surface by the different-scale models become very similar. It is clearly shown that, in a macroscopic sense, the rotation of the anisotropic yield surface is a main cause of the axial effects in large plastic shear deformation.     

Forming simulation of aluminum sheets using an anisotropic yield function coupled with crystal plasticity theory

This paper describes the application of a coupled crystal plasticity based microstructural model with an anisotropic yield criterion to compute a 3D yield surface of a textured aluminum sheet (continuous cast AA5754 aluminum sheet). Both the in-plane and out-of-plane deformation characteristics of the sheet material have been generated from the measured initial texture and the uniaxial tensile curve along the rolling direction of the sheet by employing a rate-dependent crystal plasticity model. It is shown that the stress–strain curves and R-value distribution in all orientations of the sheet surface can be modeled accurately by crystal plasticity if a “finite element per grain” unit cell model is used that accounts for non-uniform deformation as well as grain interactions. In particular, the polycrystal calculation using the Bassani and Wu (1991) single crystal hardening law and experimental electron backscatter data as input has been shown to be accurate enough to substitute experimental data by crystal plasticity data for calibration of macroscopic yield functions. The macroscopic anisotropic yield criterion CPB06ex2 (Plunkett et al., 2008) has been calibrated using the results of the polycrystal calculations and the experimental data from mechanical tests. The coupled model is validated by comparing its predictions with the anisotropy in the experimental yield stress ratio and strain ratios at 15% tensile deformation. The biaxial section of the 3D yield surface calculated directly by crystal plasticity model and that predicted by the phenomenological model calibrated with experimental and crystal plasticity data are also compared. The good agreement shows the strength of the approach. Although in this paper, the Plunkett et al. (2008) yield function is used, the proposed methodology is general and can be applied to any yield function. The results presented here represent a robust demonstration of implementing microscale crystal plasticity simulation with measured texture data and hardening laws in macroscale yield criterion simulations in an accurate manner.     

An admissible form of an inelastic constitutive equation—II. Examples of elastoplastic constitutive equation

Some examples of elastoplastic constitutive equation are presented using the general theory reported in the preceding paper (Part I). Some examinations of them are given to show that the theory is self-consistent and useful especially for anisotropic materials or materials with anisotropy resulting from plastic deformation. Mises' and Yoshimura's yield functions and a kind of quadratic function are adopted as the yield function. Formulae of r-value after arbitrary pre-straining are given which are of paramount importance in the field of press-forming of sheet metals. Several examples of stress-strain curves for various loading paths are also given.     

Anisotropic response of high-purity α-titanium: Experimental characterization and constitutive modeling

This paper presents a comprehensive experimental and theoretical investigation of the deformation behavior of high-purity, polycrystalline α-titanium under quasi-static conditions at room temperature. The initial material in this study was a cross-rolled plate with a strong basal texture. To quantify the plastic anisotropy and the tension–compression asymmetry of this material, monotonic tensile and compressive tests were conducted, on samples cut along different directions of the plate. A new anisotropic elastic/plastic model was developed to describe the quasi-static macroscopic response of the aggregate. Key in its formulation is the use of an anisotropic yield criterion that captures strength-differential effects and an anisotropic hardening rule that accounts for texture evolution associated to twinning. A very good agreement between FE simulations using the model developed and uniaxial data was obtained.     

On the non-oscillation criterion for multiplicative anisotropic plasticity at large simple shear deformation

The criterion for non-oscillatory stresses under monotonic large simple shear deformation in the context of multiplicative anisotropic plasticity is discussed. In particular, evolving anisotropy combined with a Hill type of yield criterion is considered. It is shown that a sufficient, but not necessary, criterion for a non-oscillatory stress is ellipticity of the first Piola–Kirchhoff stress. Loss of ellipticity corresponds to a critical value h_cr of the generalized plastic modulus. Similarly, the absence of limit points on the stress–strain relation motivates an alternative criterion in terms of a critical value h_sh  h_cr. Finally, this criterion is demonstrated analytically as well as numerically for an important class of models with evolving anisotropy of the saturation type.     

An approximate lower bound damage-based yield criterion for porous ductile sheet metals

An approximate lower bound damage-based yield criterion is developed for isotropic porous ductile sheet metals. The matrix of the sheet metals is assumed to be elastic-perfectly plastic and obey the von Mises yield criterion with periodically distributed voids. Gurson’s unit-cell model is simplified to characterize the sheet metals. To accommodate biaxial loading, an approach of stress superposition is adopted for the stress analysis. Numerical results were calculated and compared to Gurson’s extended yield criterion and experimental results.     

On the implementation of anisotropic yield functions into finite strain problems of sheet metal forming

In this article the implementation of anisotropic yield functions into finite element investigations of orthotropic sheets with planar anisotropy is discussed within a plane-stress context. Special attention is focused on the proper treatment of the orientation of the anisotropic axes during deformation into the finite-strain range. As an example problem the hydrostatic bulging of a membrane is considered in conjunction with a recently proposed anisotropic yield function. It is shown that the aspects of the plane-stress assumption, which do not come into consideration in isotropic analyses, can play an important role on the accuracy of the solution when the rotation of the orthotropic axes enters the computation directly due to the presence of material anisotropy. When the material anisotropy is considered and when the deformation of the workpiece is not limited to the plane of the undeformed sheet (such as cup drawing, hydrostatic bulging, etc.), the numerical experiments indicate that the only correct formulation is the one based on numerically imposing the requirement that for the plane-stress application, the in-plane material axes have to remain in the plane of the sheet during the deformation.     

Texture development and plastic anisotropy of B.C.C. strain hardening sheet metals

A Taylor-like polycrystal model is adopted here to investigate the plastic behavior of body centered cubic (b.c.c.) sheet metals under plane-strain compression and the subsequent in-plane biaxial stretching conditions. The <111> pencil glide system is chosen for the slip mechanism for b.c.c. sheet metals. The {110} <111> and {112} <111> slip systems are also considered. Plane-strain compression is used to simulate the cold rolling processes of a low-carbon steel sheet. Based on the polycrystal model, pole figures for the sheet metal after plane-strain compression are obtained and compared with the corresponding experimental results. Also, the simulated plane-strain stress—strain relations are compared with the corresponding experimental results. For the sheet metal subjected to the subsequent in-plane biaxial stretching and shear, plastic potential surfaces are determined at a given small amount of plastic work. With the assumption of the equivalence of the plastic potential and the yield function with normality flow, the yield surfaces based on the simulations for the sheet metal are compared with those based on several phenomenological planar anisotropic yield criteria. The effects of the slip system and the magnitude of plastic work on the shape and size of the yield surfaces are shown. The plastic anisotropy of the sheet metal is investigated in terms of the uniaxial yield stresses in different planar orientations and the corresponding values of the anisotropy parameter R, defined as the ratio of the width plastic strain rate to the through-thickness plastic strain rate under in-plane uniaxial tensile loading. The uniaxial yield stresses and the values of R at different planar orientations from the polycrystal model can be fitted well by a yield function recently proposed by Barlat et al. (1997b).     

Constitutive response of idealized granular media under the principal stress axes rotation

This paper is dedicated to the understanding of the phenomena, which give rise to anisotropy and non-coaxiality in granular materials. In achieving three-dimensional numerical simulation under static condition of granular media, granular element method (GEM) is adopted in this study. The method has been incorporated into the so-called mathematical homogenization theory for quasi-static equilibrium problems, which enables us to obtain the macroscopic/phenomenological inelastic deformation response of a representative volume element (RVE). To examine the anisotropic macroscopic deformation properties of the assumed RVE, which is solved by granular element method (GEM), a series of numerical experiments involving the pure rotation of the principal stress axes are carried out, and its results are discussed in relation to induced anisotropy and non-coaxiality.     

A plane stress anisotropic plastic flow theory for orthotropic sheet metals

A phenomenological theory is presented for describing the anisotropic plastic flow of orthotropic polycrystalline aluminum sheet metals under plane stress. The theory uses a stress exponent, a rate-dependent effective flow strength function, and five anisotropic material functions to specify a flow potential, an associated flow rule of plastic strain rates, a flow rule of plastic spin, and an evolution law of isotropic hardening of a sheet metal. Each of the five anisotropic material functions may be represented by a truncated Fourier series based on the orthotropic symmetry of the sheet metal and their Fourier coefficients can be determined using experimental data obtained from uniaxial tension and equal biaxial tension tests. Depending on the number of uniaxial tension tests conducted, three models with various degrees of planar anisotropy are constructed based on the proposed plasticity theory for power-law strain hardening sheet metals. These models are applied successfully to describe the anisotropic plastic flow behavior of 10 commercial aluminum alloy sheet metals reported in the literature.     

材料三维微结构表征及其晶体塑性有限元模拟

材料的力学性能，尤其是在有限变形下所呈现的宏观各向异性，是材料结构设计和服役寿命考虑的关键因素。由于宏观模型不能较好地反映材料微观结构（晶粒的形貌和取向等）对宏观塑性各向异性的影响，因此，本文建立了能实际反映晶粒形貌的三维Voronoi模型，并基于晶体塑性理论对铝合金在有限变形下的响应进行计算。首先，建立反映材料微结构的代表性体积单元RVE模型进行计算，并与实验结果进行对比验证。然后，以单向拉伸为例，分析了有限变形过程中试件的晶粒形貌和取向分布等微观因素对宏观各向异性演化的影响，并从材料和结构两个层面讨论了微观结构对宏观力学性能的影响。结果表明，本文模型能够反映微观结构对宏观力学性能的影响，为实际生产制造领域构件的力学性能提供可靠的预测。     

金属塑性变形极限判据

研究出金属材料塑性变形极限判据，它为断裂力学和失效分析提供理论依据。利用这判据，裂纹尖端的临界应变和相应的临界应力可以确定。在板料成形中，拉应力下颈缩失稳的临界应变值可以定出。节省大量试验工作。