单晶体和双晶体微观层次变形行为的有限元分析

从微观层次上研究金属材料的变形行为,将位错引入到本构关系中,用硬化函数描述材料的硬化规律,考虑了变形的率相关性,采用三维模型用大变形有限单元法对单晶体在单向拉伸载荷和循环载荷作用下的变形行为、双晶体在单向拉伸作用下滑移系的开动进行了模拟计算,得到了与实验一致的计算结果。     

A general kinematical analysis of double slip

General kinematic solutions for double slip in fcc and bcc crystals are presented which are free from constitutive assumptions: that is, the analysis does not presuppose equal amounts of slipping on equallystressed slip systems, in contrast to the standard solutions (wherein Taylor hardening is implicitly assumed). The axis rotations and limiting positions on a stereographic projection are illustrated for several different slip-system combinations, initial axis positions (none on a symmetry line), and proportional slip ratios in both fcc and bcc crystals. It is suggested that the solutions have particular application to the experimental study of double slip in the tensile test of prestrained crystals.     

Kinematic,stress, and hardening analysis in finite double slip

A kinematic, stress, and hardening analysis of finite double slip in fcc crystals under axial loading is presented. The relative amounts of slip in classic 1925 experiments by Taylor and Elam (determined analytically in [Int. J. Plasticity 9 (1993) 159–179] by comparing theoretical and experimental cones of unextended directions), together with load-extension data and other measurements, are used to calculate resolved shear stress vs. slip curves and assess predictions of several finite distortional hardening theories. In particular, a new hardening rule is introduced that gives very close agreement with the anisotropic experimental results in double slip yet is consistent with the axisymmetric deformation, lattice stability, and isotropic hardening that are characteristic of fcc and bcc crystals in high symmetry axial-load orientations involving 6- and 8-fold slip.     

Theoretical latent hardening in crystals—II. BCC crystals in tension and compression

The general latent hardening law of single slip derived in the first paper of this series (Havner, Baker and Vause, 1979) is applied to an analysis of “overshooting” phenomena in bcc crystals in tension and compression. This new law, which predicts anisotropic hardening of latent slip systems, is based upon the simple theory of finite distortional crystal hardening introduced by Havner and Shalaby (1977).Because of historical ambiguities regarding identification of the slip plane in bcc metals, parallel analyses are presented corresponding to two separate criteria: (i) slip on {110}, {112} and {123} crystallographic planes only; and (ii) slip on the plane of maximum resolved shear stress containing a 〈111〉 direction. It is established that the new hardening law is a theory of “overshooting” in bcc crystals according to either identification of the slip plane.A qualitative comparison between theoretical results and two experimental papers on Fe crystals is included. The general difficulties in making comparisons with the experimental literature on finite distortional latent hardening are briefly discussed.     

The role of texture in tension–compression asymmetry in polycrystalline NiTi

The purpose of the present study is to thoroughly understand the stress–strain behavior of polycrystalline NiTi deformed under tension versus compression. To do this, a micro-mechanical model is used which incorporates single crystal constitutive relationships and experimentally measured polycrystalline texture into the self-consistent formulation. For the first time it is quantitatively demonstrated that texture measurements coupled with a micro-mechanical model can accurately predict tension/compression asymmetry in NiTi shape memory alloys. The predicted critical transformation stress levels and transformation stress–strain slopes under both tensile and compressive loading are consistent with experimental results. For textured polycrystalline NiTi deformed under tension it is demonstrated that the martensite evolution is very abrupt, consistent with the Luders type deformation experimentally observed. The abrupt transformation under tension is attributed to the fact that the majority of the grains are oriented along the [111] crystallographic direction, which is soft under tensile loading. Since single crystals of the [111] orientation are hard under compression it is also demonstrated that under compression the martensite in textured polycrystalline NiTi evolves relatively slower.     

Self-consistent modeling of large plastic deformation, texture and morphology evolution in semi-crystalline polymers

A self-consistent model for semi-crystalline polymers is proposed to study their constitutive behavior, texture and morphology evolution during large plastic deformation. The material is considered as an aggregate of composite inclusions, each representing a stack of crystalline lamellae with their adjacent amorphous layers. The deformation within the inclusions is volume-averaged over the phases. The interlamellar shear is modeled as an additional slip system with a slip direction depending on the inclusion's stress. Hardening of the amorphous phase due to molecular orientation and, eventually, coarse slip, is introduced via Arruda-Boyce hardening law for the corresponding plastic resistance. The morphology evolution is accounted for through the change of shape of the inclusions under the applied deformation gradient. The overall behavior is obtained via a viscoplastic tangent self-consistent scheme. The model is applied to high density polyethylene (HDPE). The stress-strain response, texture and morphology changes are simulated under different modes of straining and compared to experimental data as well as to the predictions of other models.     

A modified model for simulating latent hardening during the plastic deformation of rate-dependent FCC polycrystals

A strain hardening model for the plastic deformation of rate-dependent FCC crystals is proposed based on experimental observations previously reported for single crystals. This model, which is an extension of that employed by et al. [1983], includes both the self-hardening and latent hardening of the slip systems. The differential hardening of the latent systems is assumed to arise from the interaction between glide dislocations and forests. With this hardening model and a rate-sensitive crystal plasticity theory, the deformation behavior of FCC polycrystals can be predicted from the deformation response of the constituent single crystals. As examples, the uniaxial tensile behaviour of pure aluminum and copper polycrystals is simulated using the extended model, and the results are compared with published experimental data. The effects of latent hardening on polycrystal deformation, especially on flow stress and the formation of tensile textures, are discussed.     

Strain localization in ductile single crystals

This paper is concerned with an analysis of strain localization in ductile crystals deforming by single slip. The plastic flow is modelled as rate-insensitive, and localization, viewed as a bifurcation from a homogeneous deformation mode to one which is concentrated in a narrow ‘shear band’, is found to be possible only when the plastic hardening modulus for the slip system has fallen to a certain critical value h_cr, sensitive to the precise form of the constitutive law governing incremental shear. We develop the general form of this constitutive law, incorporating within it the possibility of deviations from the Schmid rule of a critical resolved shear stress, and we show that h_cr may in fact be positive when there are deviations from the Schmid rule. It is suggested that micromechanical processes such as ‘cross-slip’ in crystals provide specific cases for which stresses other than the Schmid stress may influence plastic response and, further, there is an experimental association of localization with the onset of large amounts of cross-slip. Thus, we give the specific form of h_cr for a constitutive model that corresponds to non-Schmid effects in cross-slip, and we develop a dislocation model of the process from which we estimate the magnitude of the parameters involved. The work supports the notion that localization can occur with positive strain-hardening, h_cr > 0, and the often invoked notions of the attainment of an ideally-plastic or strain-softening state for localization may be unnecessary.     

Nonlocal continuum crystal plasticity with internal residual stresses

We derive a three-dimensional constitutive theory accounting for length-scale dependent internal residual stresses in crystalline materials that develop due to a non-homogeneous spatial distribution of the excess dislocation (edge and screw) density. The second-order internal stress tensor is derived using the Beltrami stress function tensor φ that is related to the Nye dislocation density tensor. The formulation is derived explicitly in a three-dimensional continuum setting for elastically isotropic materials. The internal stresses appear as additional resolved shear stresses in the crystallographic visco-plastic constitutive law for individual slip systems. Using this formulation, we investigate two boundary value problems involving single crystals under symmetric double slip. In the first problem, the response of a geometrically imperfect specimen subjected to monotonic and cyclic loading is investigated. The internal stresses affect the overall strengthening and hardening under monotonic loading, which is mediated by the severity of initial imperfections. Such imperfections are common in miniaturized specimens in the form of tapered surfaces, fillets, fabrication induced damage, etc., which may produce strong gradients in an otherwise nominally homogeneous loading condition. Under cyclic loading the asymmetry in the tensile and compressive strengths due to this internal stress is also strongly influenced by the degree of imperfection. In the second example, we consider simple shear of a single crystalline lamella from a layered specimen. The lamella exhibits strengthening with decreasing thickness and increasing lattice incompatibility with shearing direction. However, as the thickness to internal length-scale ratio becomes small the strengthening saturates due to the saturation of the internal stress.Finally, we present the extension of this approach for crystalline materials exhibiting elastic anisotropy, which essentially depends on the appropriate Green function within φ.     

A semi-implicit integration scheme for rate independent finite crystal plasticity

An efficient new numerical integration scheme is presented for rate independent crystal plasticity theory. A key feature of this approach is the ability to identify active slip systems prior to determining their shearing rate. Options are described for various cases of slip system hardening, including self hardening and latent hardening. Alternatives for the constitutive update are explored, including hyperelasticity based on the multiplicative decomposition of the deformation gradient as well as application of the consistency condition in a much more efficient hypoelastic formulation. Several conclusions are drawn concerning the influences of elastic and plastic properties on the activation of slip systems and their subsequent shearing rates. Key among these is the fact that, once activated, shearing rates are independent of the levels of shear flow resistance on the slip systems, provided that the plastic hardening moduli are much less in magnitude than the elastic moduli, as is usually the case. Determination of active slip systems and their shearing rates depend on the degree of elastic anisotropy of the crystal, but not on the magnitude of elastic stiffness.     

有限变形下多晶晶体塑性模型算法及应用

用Sanna和Zacharia^[1]所提出的延性单晶本构模型的积分算法和Taylor多晶模型假设研究了时间步长和硬化模型的选取对多晶集合体的应力应变响应和织构演化的影响。该算法是利用变形梯度乘法分解获得弹性变形梯度演化方程，用增量迭代法积分该方程，显式更新各滑移上的临界分切剪应力。算例的结果表明该算法具有时间步大，计算效率高的特点，另外，不同硬化模型的选取对多晶集合体应力应变响应的预测有明显的影响但对织构演化的预测影响不大。     

A crystallographic model for the study of local deformation processes in polycrystals

Most engineering materials possess a polycrystalline structure. Under load the anisotropy of the constituent grains leads to strong inhomogeneities of stresses and strains on the grain level. In order to investigate the local deformation processes, a new crystallographic model for pure fcc metals in the low temperature range has been developed. It is based on the framework of crystal plasticity and uses the finite element method (FEM). The rate dependent constitutive equations consider isotropic as well as kinematic hardening, whereby the mutual interactions of dislocation processes on the different slip systems are taken into account. Comprehensive calculations show that the essential features of both single crystals—which serve as a test object for the constitutive equations—and polycrystals are reproduced correctly. Moreover the simulations allow a deeper understanding of the mechanisms that control the local deformation behaviour of metals, especially of the mutual interactions of slip system activity, local hardening and resulting local strain. Furthermore, the model may serve as a physically motivated base for a later inclusion of damage terms which allow investigations of damage and fatigue on the local scale.     

Analysis of dynamic shear bands in an FCC single crystal

We study plane strain dynamic thermomechanical deformations of an fcc single crystal compressed along the crystallographic direction [010] at an average strain rate of 1000 sec⁻¹. Two cases are studied; one in which the plane of deformation is parallel tothe plane (001) of the single crystal, and another one with deformation occuring in the plane (101̄) of the single crystal. In each case, the 12 slip systems are aligned symmetrically about the two centroidal axes. We assume that the elastic and plastic deformations of the crystal are symmetrical about these two axes. The crystal material is presumed to exhibit strain hardening, strain-rate hardening, and thermal softening. A simple combined isotropic-kinematic hardening expression for the critical resolved shear stress, proposed by Weng, is modified to account for the affine thermal softening of the material. When the deformation is in the plane (001) of the single crystal, four slip systems (111)[11̄0], (111̄)[11̄0], (11̄;1̄;)[110], and (11̄1)[110] are active in the sense that significant plastic deformations occur along these slip systems. However, when the plane of deformation is parallel to the plane (101̄;) of the single crystal, slip systems (11̄;1)[110], (11̄1)[011], (111)[11̄0], and (111)[01̄1] are more active than the other eight slip systems. At an average strain of 0.0108, the maximum angle of rotation of a slip system within a shear band, about an axis perpendicular to the plane of deformation, is found to be 20.3° in the former case, and 22.9° in the latter.     

A dislocation-based constitutive law for pure Zr including temperature effects

In this work, a single crystal constitutive law for multiple slip and twinning modes in single phase hcp materials is developed. For each slip mode, a dislocation population is evolved explicitly as a function of temperature and strain rate through thermally-activated recovery and debris formation and the associated hardening includes stage IV. A stress-based hardening law for twin activation accounts for temperature effects through its interaction with slip dislocations. For model validation against macroscopic measurement, this single crystal law is implemented into a visco-plastic-self-consistent (VPSC) polycrystal model which accounts for texture evolution and contains a subgrain micromechanical model for twin reorientation and morphology. Slip and twinning dislocations interact with the twin boundaries through a directional Hall–Petch mechanism. The model is adjusted to predict the plastic anisotropy of clock-rolled pure Zr for three different deformation paths and at four temperatures ranging from 76 K to 450 K (at a quasi-static rate of 10⁻³ 1/s). The model captures the transition from slip-dominated to twinning-dominated deformation as temperature decreases, and identifies microstructural mechanisms, such as twin nucleation and twin–slip interactions, where future characterization is needed.     

A dislocation density-based single crystal constitutive equation

Single crystal constitutive equations based on dislocation density (SCCE-D) were developed from Orowan’s strengthening equation and simple geometric relationships of the operating slip systems. The flow resistance on a slip plane was computed using the Burger’s vector, line direction, and density of the dislocations on all other slip planes, with no adjustable parameters. That is, the latent/self-hardening matrix was determined by the crystallography of the slip systems alone. The multiplication of dislocations on each slip system incorporated standard 3-parameter dislocation density evolution equations applied to each slip system independently; this is the only phenomenological aspect of the SCCE-D model. In contrast, the most widely used single crystal constitutive equations for texture analysis (SCCE-T) feature 4 or more adjustable parameters that are usually back-fit from a polycrystal flow curve. In order to compare the accuracy of the two approaches to reproduce single crystal behavior, tensile tests of single crystals oriented for single slip were simulated using crystal plasticity finite element modeling. Best-fit parameters (3 for SCCE-D, 4 for SCCE-T) were determined using either multiple or single slip stress–strain curves for copper and iron from the literature. Both approaches reproduced the data used for fitting accurately. Tensile tests of copper and iron single crystals oriented to favor the remaining combinations of slip systems were then simulated using each model (i.e. multiple slip cases for equations fit to single slip, and vice versa). In spite of fewer fit parameters, the SCCE-D predicted the flow stresses with a standard deviation of 14 MPa, less than one half that for the SCCE-T conventional equations: 31 MPa. Polycrystalline texture simulations were conducted to compare predictions of the two models. The predicted polycrystal flow curves differed considerably, but the differences in texture evolution were insensitive to the type of constitutive equations. The SCCE-D method provides an improved representation of single-crystal plastic response with fewer adjustable parameters, better accuracy, and better predictivity than the constitutive equations most widely used for texture analysis (SCCE-T).     

Constitutive analysis of elastic-plastic crystals at arbitrary strain

A time-independent constitutive framework is constructed for crystals capable of crystallographic shearing. No restrictions are placed on their elasticity or on the amounts of slip. In deriving properly-objective relations due account is taken of the separate motions of the material and lattice. A precise specification of the structure of the flow rule leads to conditions for the existence of plastic potentials, as well as to an exact statement of the slip criterion required for normality of the plastic strain-increment in conjugate variables. Objective hardening laws are also examined and sufficiency conditions obtained for uniqueness of the slip magnitudes.     

Rate-independent crystalline and polycrystalline plasticity,application to FCC materials

This paper deals with the simulation of the mechanical response and texture evolution of cubic crystals and polycrystals for a rate-independent elastic–plastic constitutive law. No viscous effects are considered. An algorithm is introduced to treat the difficult case of multi-surface plasticity. This algorithm allows the computation of the mechanical response of a single crystal. The corresponding yield surface is made of the intersection of several hyper-planes in the stress space. The problem of the multiplicity of the slip systems is solved thanks to a pseudo-inversion method. Self and latent hardening are taken into account. In order to compute the response of a polycrystal, a Taylor homogenization scheme is used. The stress–strain response of single crystals and polycrystals is computed for various loading cases. The texture evolution predicted for compression, plane strain compression and simple shear are compared with the results given by a visco-plastic polycrystalline model.     

有限弹塑性变形的三维组集式本构模型

本文将文[1]中提出的三维组集式弹塑性本构模型推广应用于有限变形分析,导出了全量型和增量型本构关系在初始构形上的拉格朗日(Total Lagrange)形式和瞬时构形上的拉格朗日(Updated Lagrange)形式。文中对晶体单轴拉伸中的宏观剪切带进行了分析。预测结果与实验吻合。从而说明这种本构模型能够模拟有限变形中的几何非线性效应和晶体材料塑性变形中的宏观力学行为。