Comparisons of crystal hardening laws in multiple slip

This paper brings together and concisely reviews results from recent analytical investigations on single crystals (variously done alone or with students) in which predictions from different theoretical hardening laws are contrasted and compared with experimental studies. Finitely deforming f.c.c. crystals in both constrained and unconstrained multiple-slip configurations are considered. Four crystal hardening laws are given prominence. Two of these belong to a class of theories in which the physical hardening moduli relating rates-of-change of critical strengths (in the 24 crystallographically equivalent slip systems) to slip-rates are taken as symmetric. These are G.I. Taylor's classic isotropic hardening rule (proposed in 1923), which is almost universally adopted in the metallurgical literature for various approximate analyses of single and poly-crystal deformation, and a 2-parameter modification of Taylor's rule that has an empirical basis in the qualitative features of experimentally determined latent hardening in single slip. The other two hardening laws featured here belong to a class of theories that were introduced in 1977 by this author. This class requires the above moduli to be nonsymmetric and explicity dependent upon the current stress state in such a manner that the following consequences are assured. (1) The deformation-dependent hardening of latent slip systems necessarily develops anisotropically if there is relative rotation of gross material and underlying crystal lattice. (2) The theories admit self-adjoint boundary value problems for crystalline aggregates, hence a variational formulation. (The fact that symmetric physical hardening moduli do not permit variational formulations of polycrystalline problems was shown at the 1972 Warsaw Symposium.) The two members of this class considered here are the original (and simplest p possible) theory of rotation-dependent anisotropy, which was proposed by this author in 1977 and commonly has been referred to as the “simple theory,” and a modification of this theory introduced in 1982 by Peirce, Asaro and Needleman that lies between Taylor's rule and the simple theory in its predictions for finitely deforming f.c.c. crystals. (In a series of five papers during 1977–1979, the simple theory was shown to universally account for the experimental phenomenon of “overshooting” in single slip in both f.c.c. and b.c.c. crystals.) Theoretical results from the various hardening rules are contrasted and compared with finite strain experiments in the metallurgical literature. Both tensile-loaded crystals in 4, 6 and 8-fold symmetry orientations and compressively loaded crystals under conditions of channel die constraint are treated. A postulate of minimum plastic work introduced in 1981 plays a prominent role in the theoretical analyses, in many cases providing a unique solution to the slip system inequalities and deformation constraints (where applicable). The rather remarkable ability of the simple theory to reconcile diverse qualitative features of both constrained and unconstrained finited deformation of f.c.c. crystals is demonstrated. Finally, conditions for total loading (all systems active) in 6-fold symmetry are investigated, and certain concepts regarding the selection of active systems under prescribed straining are critically assessed.     

A cyclic plasticity model for single crystals

The Armstrong–Frederick type kinematic hardening rule was invoked to capture the Bauschinger effect of the cyclic plastic deformation of a single crystal. The yield criterion and flow rule were built on individual slip systems. Material memory was introduced to describe strain range dependent cyclic hardening. The experimental results of copper single crystals were used to evaluate the cyclic plasticity model. It was found that the model was able to accurately describe the cyclic plastic deformation and properly reflect the dislocation substructure evolution. The well-known three distinctive regimes in the cyclic stress–strain curve of the copper single crystals oriented for single slip can be reproduced by using the model. The model can predict the enhanced hardening for crystals oriented for multislip, showing the model's ability to describe anisotropic cyclic plasticity. For a given loading history, the model was able to capture not only the saturated stress–strain response but also the detailed transient stress–strain evolution. The model was used to predict the cyclic plasticity under a high–low loading sequence. Both the stress–strain responses and the microstructural evolution can be appropriately described through the slip system activation.     

On the uniqueness of a rate-independent plasticity model for single crystals

This paper presents the uniqueness and existence conditions for a rate-independent plasticity model for single crystals under a general stress state. The model is based on multiple slips on three-dimensional slip systems. The uniqueness condition for the plastic slips in a single crystal with nonlinear hardening is derived using the implicit function theorem. The uniqueness condition is the non-singularity of a matrix defined by the Schmid tensors, the elasticity, and the hardening rates of the slip systems. When this matrix becomes singular, the limitations on the loading paths that can be accommodated by the active slip systems (i.e., the existence conditions) are also given explicitly. For the compatible loading paths, a particular solution is selected by requiring the solution vector to be orthogonal to the null space of the singular coefficient matrix. The paper also presents a fully implicit algorithm for the plasticity model. Numerical examples of an fcc copper single crystal under cyclic loadings (pure shear and uniaxial strain) are presented to demonstrate the main features of the algorithm.     

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.     

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 study of directionality of latent hardening behavior in aluminium single crystals

By defining the yield stress in a latent hardening test as the first deviation from the elastic straight line, the yielding and hardening behavior on a latent system in the positive and negative slip direction was studied in aluminium single crystals. It is shown that the yield stresses on both the positive and negative latent systems are about equal to or a little lower than the maximum resolved shear stress in the primary test, but much higher than that of the active system. The influence of relative orientation and prestrain on latent hardening and initial work-hardening in the secondary test was also investigated, and it was found that there is a considerable effect on initial work-hardening, but none on latent hardening. With reasonable approximation, a hardening rule for single crystal could be proposed from the experimental results, that is, except for the yield stress on the system negative to the active system that is very low, hardening on the other systems is nearly the same as self-hardening.     

Finite versus small strain discrete dislocation analysis of cantilever bending of single crystals

Plastic size effects in single crystals are investi-gated by using finite strain and small strain discrete dislo-cation plasticity to analyse the response of cantilever beam specimens. Crystals with both one and two active slip sys-tems are analysed, as well as specimens with different beam aspect ratios. Over the range of specimen sizes analysed here, the bending stress versus applied tip displacement response has a strong hardening plastic component. This hardening rate increases with decreasing specimen size. The hardening rates are slightly lower when the finite strain discrete disloca-tion plasticity (DDP) formulation is employed as curving of the slip planes is accounted for in the finite strain formulation. This relaxes the back-stresses in the dislocation pile-ups and thereby reduces the hardening rate. Our calculations show that in line with the pure bending case, the bending stress in cantilever bending displays a plastic size dependence. How-ever, unlike pure bending, the bending flow strength of the larger aspect ratio cantilever beams is appreciably smaller. This is attributed to the fact that for the same applied bend-ing stress, longer beams have lower shear forces acting upon them and this results in a lower density of statistically stored dislocations.     

基于能量极值原理的单晶体塑性滑移的有限变形分析

对FCC单晶体的率无关弹塑性力学响应的本构关系进行了数值模拟。用一个基于能量极值原理的数值计算方法来处理复杂的多面塑性问题，这种算法可以有效地模拟单晶体多滑移系的启动，在这一理论框架下，增量的应力应变关系可以从所构造的能量函数中推导出来。对于滑移系激活情况的判定则可转化为求解活动约束的非线性数学规划问题，通过对时间的离散，此问题又可细化为逐步二次规划问题，并采用有效集法来搜索启动滑移系，进而求得弹塑性本构关系，数值结果表明该方法具有稳定、收敛、可行的特点，在数值计算的基础上研究了单轴拉伸下晶体的不同取向对单晶体硬化程度和滑移系激活情况的影响。     

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

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

On the Torsion of a Class of Nonlinear Fiber Composite Cylinders

This paper presents an analysis of the torsion of a solid or annular circular cylinder consisting of nonlinear material in the form of an elastic matrix with embedded unidirectional elastic fibers parallel to the cylinder axis. The specific class of composite considered is one for which nonlinear fiber-matrix interface slip is captured by uniform cohesive zones of vanishing thickness. Previous work on the effective antiplane shear response of this material leads to a stress–strain relation depending on the interface slip together with an integral equation governing its evolution. Here, we obtain an approximate single mode solution to the integral equation and utilize it to solve the torsion problem. Equations governing the radial distributions of shear stress and interface slip are obtained and formulae for torque–twist rate are presented. The existence of singular surfaces, i.e., surfaces across which the slip and the shear stress experience jump discontinuities are analyzed in detail. Specific results are presented for an interface force law that allows for interface failure in shear.     

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.     

Anisotropic hardening in single crystals and the plasticity of polycrystals

Based on the dislocation structures developed during plastic deformation, an anisotropic hardening law is developed to describe the latent hardening behavior of slip systems under multislip. This theory incorporates the concept of isotropic hardening, kinematic hardening, and the two-parameter representation; it automatically includes the strength differential between the forward and reversed slips and between the acute and obtuse cross slips. The self-hardening modulus of a slip system is found to be “associated” with the latent hardening law involved, and, based on some experimental evidence, two specific sets of self-hardening modulus are suggested. An important feature of this associated modulus is that the slip system with a soft latent hardening (e.g., the reversed system with a Bauschinge effect) will have an enhanced self-hardening modulus. This newly developed hardening law, together with its associated latent hardening moduli, is then applied to examine the strain-hardening behavior of a polycrystal. Although crystals with a stronger latent hardening will, in general, also lead to a stranger strain-hardening for the polycrystal, the stress-strain behavior of the polycrystal using the kinematic hardening law of single crystals is found to be not necessarily softer than that using the isotropic hardening law. Within the range of experimentally measured latent hardening ratio of slip systems, the anisotropic theory is also used to calculate the motion of yield surface of a polycrystal. The general results, employing four selected types of anisotropic hardening, all show the essential features of experimental observations by Phillips and his co-workers. The application is highlighted with a reasonably successful quantitative modeling of initial and subsequent yield surfaces of an aluminum.     

Plane constrained shear of single crystal strip with two active slip systems

Within continuum dislocation theory the plane constrained shear of a single crystal strip with two active slip systems is considered. An analytical solution is found for symmetric double slip which exhibits the energetic and dissipative thresholds for dislocation nucleation, the Bauschinger translational work hardening, and the size effects. Comparison with discrete dislocation simulations shows good agreement between the discrete and continuum approaches. Numerical procedures in the general case of non-symmetric double slip are proposed.     

Theoretical latent hardening in crystals—I. General equations for tension and compression with application to F.C.C. crystals in tension

Equations for latent strengths in single slip, based upon the simple theory of finite distortional crystal hardening introduced by K.S. Havner and A.H. Shalaby (1977), are derived for both tensile and compression tests without restriction as to crystal class. Detailed comparisons between theoretical results and the experiments of P.J. Jackson and Z.S. Basinski (1967) on copper crystals in tension are presented. There is good qualitative agreement between theory and experiment regarding the diversity of anisotropic hardening among slip systems. Moreover, there is satisfactory quantitative agreement between the theory and the extrapolated experimental data in the stage III, large-strain range. It is suggested that further experimental investigation of latent hardening at large prestrains would be desirable.The simple theory predicts anisotropic hardening and the perpetuation of single slip in axial loading of cubic crystals initially oriented for single slip, but predicts symmetric, isotropic hardening of specimens initially oriented in positions of 4, 6 or 8-fold multiple-slip. These predictions are in general accord with experimental observations from tests of f.c.c. and b.c.c. crystals.     

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.     

Shear band bifurcations in ductile single crystals

A shear band bifurcation analysis for a multiply slipping ductile single crystal is presented. The crystal is rate independent, and crystallographic slip is governed by Schmid's law. Attention is focused on double slip in a tensile deformation, so that comparison with previous studies is feasible. The results, which consider the rigid-plastic and elastic-plastic cases separately, show several qualitative features in agreement with experiment and with earlier analytical work based on a plane strain model. The shear band orientations exhibit the expected tendency, but the critical strain hardening values obtained here, while positive, are quite low. The present results also indicate the ready availability of non-uniform slip patterns, or ‘patchy slip’, a phenomenology consistent both with experimental observations and recent numerical studies.     

The effect of biaxial loading on the critical resolved shear stress of zinc single crystals

This paper reports an experimental investigation to test the hypothesis that the critical resolved shear stress is independent of the resolved normal stress on the active slip system plane. The resolved normal stress is varied independently for crystals of constant orientation by applying a uniform biaxial tension stress through slitted rubber strips glued to a flat tension sample. It is found that the critical resolved shear stress is constant until the resolved normal stress on the active slip system becomes twice the resolved shear stress; and then, for further increases in the resolved normal stress, the critical resolved shear stress begins to decrease.     

铝单晶体潜在硬化的指向性行为研究

应用高纯铝单晶体，采用偏离弹性线法定义其在潜在滑移系统在屈服应力，研究了潜在滑移系在正负两个滑移方向上的屈服及应变硬化行为。结果表明，潜在滑移系负行为的差异要远远小于动滑移系的Ｂａｕｓｃｈｉｎｇｅｒ效应，其正负方向的屈服应力一般均等于或略小于预应变时的最大分切应力，大大高于动滑移系的负向屈服应力，潜在滑移系和原始滑移系的相对取向及预应变的大小对单晶体潜在移系在潜在硬化的影响不大，但对其起始过渡区应     

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.