立方晶粒正交板材瑞利波速与织构系数的关系

正交板材是由大量微小晶粒组成的正交多晶体材料,而多晶体中晶粒的取向分布(可通过取向分布函数中的织构系数来描述)影响着多晶体材料的力学性能,也必然影响着瑞利波的传播速度。将多晶体材料的织构系数引入到弹性张量中,通过特征值办法,采用线性化处理,推导出立方晶粒正交板材的瑞利波速与织构系数的关系式,在此基础上可通过正交板材瑞利波速的测量获得织构系数,并与通过超声横波纵波测得结果相比,吻合很好。     

A methodology for measuring and modeling crystallographic texture gradients in processed alloys

The explicit representation of internal material structure in alloy processing and in-service performance simulations is becoming increasingly prevalent. This paper presents a methodology for characterizing and representing a spatially-varying orientation distribution function (ODF) that can be used in processing and performance simulations for alloys containing texture gradients. We use thick AA 7050 aluminum plate, which is known to contain texture gradients, as a case study to demonstrate the methodology, which employs a finite element representation of the ODF initialized using individual lattice orientation measurements taken using the electron backscatter pattern (EBSP) technique. As expected, we find that the texture varies significantly through the plate thickness. We use the ODF to examine the effect of the varying texture on the resulting yield strength distribution as embodied by the average Taylor factor. We find that the predicted yield strength anisotropy is different at different locations through the thickness of the plate. We examine the optimal number of orientation measurements necessary for determining the ODF in the presence of this texture gradient. We find that as we increase the number of orientations, the ODF quickly becomes stable but eventually starts to change under the influence of the texture gradient. We also investigate spatial interpolation of the ODF using the finite element representation. We find that, as with finite element representations of other fields, interpolation accuracy depends on the variation of the field variable and the discretization of the domain. In this case, gradients in both physical space and orientation space affect the accuracy of the interpolation. Finally, the effects of the texture gradient on the mechanical response of the material is demonstrated by employing the ODFs taken from various locations through the thickness of the plate in polycrystal plasticity simulations of uniaxial tension and plane strain compression.     

Three-dimensional deformation process simulation with explicit use of polycrystal plasticity models

The combination of massive parallel processing and polycrystal plasticity theory offers the potential for applying detailed microstructural models to macroscopic deformation processes. In this work the finite element method is used to solve for the three-dimensional deformation of a plastic workpiece. The elemental constitutive response is derived from the microstructural response of a polycrystal aggregate situated in the element. Crystal orientations and their respective weighted contributions to the aggregate response are selected to approximate the orientation distribution derived from experimental pole figure measurements. The interaction of the material symmetry adopted in analysis of pole figures and the boundary conditions posed in the plasticity boundary value problem are examined. Through the introduction of distinct aggregates with decreasing crystal to aggregate ratio, an inhomogenous material response is developed where: (1) the orientation distribution becomes well approximated only by a collection of spatially distinct aggregates, and (2) these aggregates experience deformation paths of increasing variation. It is shown that the use of spatially distinct aggregates in a material experiencing local kinematic inhomogeneities throughout its deformation history leads to texture predictions that compare favorably with experimental measurements.     

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

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

Anisotropic yield and plastic flow of polycrystalline solids

Deformation induced anisotropy in polycrystalline solids results mainly from crystallographic slip due to dislocation motion at the grain level and texture development due to grain rotation at the aggregate level. To describe these characteristics, the so-called scale invariance approach is adopted which allows information and constitutive relations pertaining to single slip to be cast in a form of macroscopic constitutive equations. An orientation distribution function (ODF) and a texture tensor are introduced into the earlier version (based on the hypotheses of single slip at the grain level and isotropic distribution of the crystallites at the aggregate level) of the scale invariance framework to describe texture effects in plastically deformed polycrystals. The texture tensor is calculated either directly through the solution of ODF, or indirectly through an appropriate set of evolution equations for the orientation tensors and the use of a closure approximation. Theoretical predictions for anisotropic yield and plastic flow behavior compare well with available experimental data.     

Anisotropic yield criterion for polycrystalline metals using texture and crystal symmetries

An anisotropic yield criterion for polycrystalline metals which uses texture data and takes advantage of crystal symmetries is presented. A linear transformation is developed to map an anisotropic yield surface for a polycrystal to an appropriate isotropic yield surface. The transformation developed reflects the symmetry of the material being modeled. First, the transformation is determined. Then, information regarding the orientation distribution (texture) of the crystals in a polycrystalline aggregate is used to determine, via averaging, the transformation for the polycrystal. The transformation, along with appropriate isotropic yield surface, provides a phenomenological approach to modeling yield, yet accounts for microstructural texture. The approach reduces to the Hill (1950) anisotropic plasticity theory under certain conditions. The yield surfaces and R-values for various face-centered-cubic ( fcc) polycrystalline textures are computed by this method. Results compare favorably with those given by other theories, and with experiment. The method proves to have the computational efficiency of phenomenological approaches to modeling yield, while effectively incorporating the physics of more complex crystallographic approaches.     

织构可调塑性本构模型标定及其应用

提出了利用率相关晶体塑性模型标定织相可调本构模型的求解步骤,得出了一组依赖于晶粒间相互作用假设而独立于具体板材织构的本构相关系数.以此为基础再结合板材织构系数所得出的本构模型系数可避免出现屈服面非外凸的情形.利用所提求解步骤对在不同热处理条件下产生不同织构的AL5052铝合金板的深拉成形过程进行了有限元模拟.结果再现了典型织构在板材成形过程中所出现的塑性各向异性,从而表明求解步骤的可行性.     

Prediction of stored energy in polycrystalline materials during cyclic loading

The effect of initial texture on the stored energy is investigated. Uniaxially loaded polycrystalline materials with initial textures based on the Goss component and the Brass component are analyzed. For reference purposes a single crystal and an initial isotropic crystal orientation distribution are also analyzed. Special attention is directed at the thermomechanical behavior of polycrystalline material during cyclic loading, the temperature evolution and change in stored energy are studied. Cyclic loading of Cook’s membrane is also considered. The simulations are done using a rate-dependent crystal plasticity model for large deformations formulated within a thermodynamic framework. It is shown that incorporation of the latent-hardening into the Helmholtz free energy function and use of evolution laws of appropriate form allows a thermodynamically consistent heat generation due to plastic work.     

STUDY OF TEXTURE EFFECT ON STRAIN LOCALIZATION OF BCC STEEL SHEETS

Using elastic crystalline viscoplastic finite element (FE) annlysis, the formability of BCC steel sheets was assessed. An orientation probability assignment method in the FE modeling procedure, which can be categorized as an inhomogenized material modeling, was newly proposed. In the study, the crystal orientations of three materials, mild steel, dual phase steel and the high strength steel, were obtained by X-ray diffraction and orientation distribution function (ODF) analyses. The measured ODF results have revealed clearly different textures in the sheets, featured by orientation fibers, skeleton lines and selected orientations in Euler angle space, which are closely related to the plastic anisotropy. Then, the crystal orientations were assigned to FE integration points by using this ODF data, individually. The FE analyses of the standard limiting dome height(LDH) test show how the fiber textures affect the extent of strain localization in the forming processes. It was confirmed by comparison with experimental results that this FE code could predict the extreme strain localization and assess the sheet formability. The third author is indebted to the National Natural Science Foundation of China for financial support (Grant No. 59875025) to the research cooperation with OIT, Japan.     

Green's Function for Orthorhombic Aggregates of Weakly Anisotropic Cubic Crystals

Herein a closed but approximate formula of the Green's function is obtained for orthorhombic aggregates of cubic crystallites. This formula, which includes three material constants and three texture coefficients, accounts for the effects of the orientation distribution function (ODF) up to terms linear in the texture coefficients. Thus it is expected that our formula would be applicable to aggregates with weak texture or to materials such as aluminum whose single crystal has weak anisotropy. The approximate formula remains valid and assumes a simpler form when the polycrystal reduces to a weakly anisotropic cubic crystal. Two examples are presented to compare predictions from our formula with those from Nishioka and Lothe's formula and from Synge's contour integral through numerical integration.     

Effects of spatial grain orientation distribution and initial surface topography on sheet metal necking

The finite element method is used to numerically simulate localized necking in AA6111-T4 under stretching. The measured EBSD data (grain orientations and their spatial distributions) are directly incorporated into the finite element model and the constitutive response at an integration point is described by the single crystal plasticity theory. We assume that localized necking is associated with surface instability, the onset of unstable growth in surface roughening. It is demonstrated that such a surface instability/necking is the natural outcome of the present approach, and the artificial initial imperfection necessitated by the macroscopic M–K approach [Marciniak and Kuczynski (1967). Int. J. Mech. Sci. 9, 609–620] is not relevant in the present analysis. The effects of spatial orientation distribution, material strain rate sensitivity, texture evolution, and initial surface topography on necking are discussed. It is found that localized necking depends strongly on both the initial texture and its spatial orientation distribution. It is also demonstrated that the initial surface topography has only a small influence on necking.     

Finite element analysis of localization in FCC polycrystalline sheets under plane stress tension

Localization phenomena in thin sheets subjected to plane stress tension are investigated. The sheet is modelled as a polycrystalline aggregate, and a finite element analysis based on rate-dependent crystal plasticity is developed to simulate large strain behaviour. Accordingly, each material point in the specimen is considered to be a polycrystalline aggregate consisting of a large number of FCC grains. The Taylor model of crystal plasticity theory is assumed. This analysis accounts for initial textures as well as texture evolution during large plastic deformations. The numerical analysis incorporates certain parallel computing features. Simulations have been carried out for an aluminum sheet alloy, and the effects of various parameters on the formation and prediction of localized deformation (in the form of necking and/or in-plane shear bands) are examined.     

FCC晶体中不均匀变形对空洞长大的影响

通过编制率相关有限元用户子程序,采用一个单胞模型研究了FCC晶体中孔洞在单晶及晶界的长大行为,分析了由于晶体取向及变形失配对孔洞长大和聚合的影响。研究结果表明：孔洞的形状和长大方向与晶体取向密切相关;晶界上孔洞的长大速度大于单晶中孔洞的长大速度;晶粒间的变形失配加速了晶界上孔洞的长大趋势,因而使材料易发生沿晶断裂,随着晶粒间取向因子差异的增加,孔洞越易沿着晶界长大。     

A new crystal plasticity scheme for explicit time integration codes to simulate deformation in 3D microstructures: Effects of strain path,strain rate and thermal softening on localized deformation in the aluminum alloy 5754 during simple shear

The predominant deformation mode during material failure is shear. In this paper, a crystal plasticity scheme for explicit time integration codes is developed based on a forward Euler algorithm. The numerical model is incorporated in the UMAT subroutine for implementing rate-dependent crystal plasticity model in LS-DYNA/Explicit. The sheet is modeled as a face centered cubic (FCC) polycrystalline aggregate, and a finite element analysis based on rate-dependent crystal plasticity is implemented to analyze the effects of three different strain paths consisting predominantly of shear. Finite element meshes containing texture data are created with solid elements. The material model can incorporate information obtained from electron backscatter diffraction (EBSD) and apply crystal orientation to each element as well as account for texture evolution. Single elements or multiple elements are used to represent each grain within a microstructure. The three dimensional (3D) polycrystalline microstructure of the aluminum alloy AA5754 is modeled and subjected to three different strain rates for each strain path. The effects of strain paths, strain rates and thermal softening on the formation of localized deformation are investigated. Simulations show that strain path is the most dominant factor in localized deformation and texture evolution.     

Crystal plasticity based analysis of localized necking in aluminum tube under internal pressure

The finite element method is used to numerically simulate localized necking in aluminum alloy tube under internal pressure. The measured electron backscatter diffraction (EBSD) data are directly incorporated into the finite element model and the constitutive response at an integration point is described by the single crystal plasticity theory. The tube is assumed sufficiently long, so that length changes as well as the end effects can be ignored and a plane strain analysis can be performed. Localized necking is assumed to be associated with surface instability, the onset of unstable thinning. It is demonstrated that such a surface instability/necking is the natural outcome of the present approach, and an artificial initial imperfection required by other approaches is not necessary in the present analysis. The effects of spatial grain orientation distribution, material strain rate sensitivity, work hardening, and initial surface topography on necking are discussed. It is found that localized necking depends strongly on both the initial texture and its spatial orientation distribution, while the initial surface topography has a negligible effect on necking.     

Macroscopic theory of evolution of deformation textures

This paper contains the derivation of a new equation for the orientation distribution function (texture evolution) in metallic polycrystals undergoing large plastic deformations in the range of low temperatures. In contrast to other theories, the equation proposed in this paper is hyperbolic. It is shown that solutions of this equation qualitatively reflect properties of textures very well, such as saturation at large plastic deformations. It is also demonstrated that plastic rotations play a particulary important role in time changes of textures. Ann effective analytic solution is constructed for plane textures. A brief discussion of averaging by means of the orientation distribution function shows that the evolution of textures yields entirely different averages of the rate of deformation as those used in crystal plasticity. Simultaneously, the multiplicative decomposition of the deformation gradient, commonly used in theories of large plastic deformations, is shown to be unacceptable in the case of time-dependent textures.     

On the strain hardening and texture evolution in high manganese steels: Experiments and numerical investigation

We present a systematic investigation on the strain hardening and texture evolution in high manganese steels where twinning induced plasticity (TWIP) plays a significant role for the materials' plastic deformation. Motivated by the stress–strain behavior of typical TWIP steels with compositions of Fe, Mn, and C, we develop a mechanistic model to explain the strain-hardening in crystals where deformation twinning dominates the plastic deformation. The classical single crystal plasticity model accounting for both dislocation slip and deformation twinning are then employed to simulate the plastic deformation in polycrystalline TWIP steels. While only deformation twinning is activated for plasticity, the simulations with samples composed of voronoi grains cannot fully capture the texture evolution of the TWIP steel. By including both twinning deformation and dislocation slip, the model is able to capture both the stress–strain behaviors and the texture evolution in Fe–Mn–C TWIP steel in different boundary-value problems. Further analysis on the strain contributions by both mechanisms suggests that deformation twinning plays the dominant role at the initial stage of plasticity in TWIP steels, and dislocation slip becomes increasingly important at large strains.     

PREDICTION OF YIELD FUNCTIONS ON BCC POLYCRYSTALS

By the nonlinear optimization theory, we predict the yield function of single BCC crystals in Hill's criterion form. Then we give a formula on the macroscopic yield function of a BCC polycrystal Ω under Sachs' model, where the volume average of the yield functions of all BCC crystallites in Ω is taken as the macroscopic yield function of the BCC polycrystal. In constructing the formula, we try to find the relationship among the macroscopic yield function, the orientation distribution function (ODF), and the single BCC crystal's plasticity. An expression for the yield stress of a uniaxial tensile problem is derived under Taylor's model in order to compare the expression with that of the macroscopic yield function.     

两种晶体塑性损伤模型在有限变形条件下的比较

对建立在连续损伤力学和内变量理论基础上的两种晶体塑性损伤模型进行了比较,旨在比较两种模型在描述材料物理性能方面的适用性以及由加载而引起变形响应的不同之处.模拟结果显示,两种模型均能反映出塑性各向异性和损伤的演化;由加载而导致有限变形的响应不仅依赖于变形,而且也依赖于晶格取向;尽管两种模型在揭示单晶体的物理性能方面是不同的,但是在预测材料力学性能方面有着相同的预测趋势.