Strain gradient crystal plasticity effects on flow localization

In metal grains one of the most important failure mechanisms involves shear band localization. As the band width is small, the deformations are affected by material length scales. To study localization in single grains a rate-dependent crystal plasticity formulation for finite strains is presented for metals described by the reformulated Fleck–Hutchinson strain gradient plasticity theory. The theory is implemented numerically within a finite element framework using slip rate increments and displacement increments as state variables. The formulation reduces to the classical crystal plasticity theory in the absence of strain gradients. The model is used to study the effect of an internal material length scale on the localization of plastic flow in shear bands in a single crystal under plane strain tension. It is shown that the mesh sensitivity is removed when using the nonlocal material model considered. Furthermore, it is illustrated how different hardening functions affect the formation of shear bands.     

Effects of texture on shear band formation in plane strain tension/compression and bending

In this study, effects of typical texture components observed in rolled aluminum alloy sheets on shear band formation in plane strain tension/compression and bending are systematically studied. The material response is described by a generalized Taylor-type polycrystal model, in which each grain is characterized in terms of an elastic–viscoplastic continuum slip constitutive relation. First, a simple model analysis in which the shear band is assumed to occur in a weaker thin slice of material is performed. From this simple model analysis, two important quantities regarding shear band formation are obtained: i.e. the critical strain at the onset of shear banding and the corresponding orientation of shear band. Second, the shear band development in plane strain tension/compression is analyzed by the finite element method. Predictability of the finite element analysis is compared to that of the simple model analysis. Third, shear band developments in plane strain pure bending of a sheet specimen with the typical textures are studied. Regions near the surfaces in a bent sheet specimen are approximately subjected to plane strain tension or compression. From this viewpoint, the bendability of a sheet specimen may be evaluated, using the knowledge regarding shear band formation in plane strain tension/compression. To confirm this and to encompass overall deformation of a bent sheet specimen, including shear bands, finite element analyses of plane strain pure bending are carried out, and the predicted shear band formation in bent specimens is compared to that in the tension/compression problem. Finally, the present results are compared to previous related studies, and the efficiency of the present method for materials design in future is discussed.     

动态载荷下剪切变形局部化、微结构演化与剪切断裂研究进展

总结和评述了近年来金属与合金变形局部化的形成、微结构演化与剪切断裂方面作者和相关的研究工作成果. 材料包括低碳钢, SS304不锈钢, Fe-15%Ni-15%Cr单晶, Al-Li合金,α-Ti和Ti-6Al-4V, Al/SiCp复合材料等.综述内容主要包括:采用改进的Hopkinson扭杆装置,对剪切变形局部化形成、发展和演化过程进行了实验观察与数值模拟;采用"侧剖"与"对接"等定点方法制备电子显微镜薄膜试样,对剪切带内相变与再结晶、非晶转变、旋涡结构等动态变形现象,以及与宏观动态力学行为对应的位错胞的形成、发展和坍塌等微结构特征进行了观测;提出了应变和应变率同时作为剪切带形成的两个必要条件的直接实验证据;在剪切带内发现了α'$-马氏体相变现象,以及相变产物与母体之间的晶体学关系;通过位错单滑移或交滑移等微观剪切最后发展成为宏观剪切的机制;对剪切带内再结晶结构的观测和对再结晶动力学本构关系的定量描述;对剪切带特别是``白色'腐蚀带(或相变带)的形成机制的分析和新的解释,指出 ``白色'是带内亚结构取向趋于一致,其在光学或扫描显微镜下很难辨认这些亚结构的取向差所致,并非表明剪切带内一定发生了相变;通过截断实验和实时跟踪观测发现,剪切带内微裂纹的萌生与聚合是材料承载能力骤然下降并导致最后断裂的主控因素.此外,本文对近年来在准静态和循环加载下材料的局部化形变与剪切断裂的实验结果予以简要评述,指出其微观机制与动态载荷下的截然不同, 是由位错的平面滑移所控制的,与热效应无关的等温变形.      

Elastic/crystalline viscoplastic finite element analyses of single- and poly-crystal sheet deformations and their experimental verification

The elastic/crystalline viscoplastic constitutive equation, based on a newly proposed hardening-softening evolution equation, is introduced into the dynamic-explicit finite element code “Itas-Dynamic.” In the softening evolution equation, the effective distance and the angle between each slip system of a crystal are introduced to elucidate the interaction between the slip systems, which causes a decrease of dislocation density. The polycrystal sheet is modeled by Voronoi polygons, which correspond to the crystal grains; and by the selected orientations, which can relate to the texture, they are assigned to the integration points of the finite elements. We propose a direct crystal orientation assignment method, which means that each integration point of finite element has an assigned orientation, and its orientation can be rotated independently. Therefore, this inhomogeneous polycrystal model can consider the plastic induced texture development and subsequent anisotropy evolution. The parameters of the constitutive equation are identified by uni-axial tension tests carried out on single crystal sheets. Numerical results obtained for sheet tensions are compared with experimental ones to confirm the validity of our finite element code. Further, we investigate the following subjects: (1) how the initial orientation of single crystal affects slip band formation and strain localization; (2) how the grain size and particular orientations of the grain affect the strain localization in case of a polycrystal sheet. It is confirmed that the orientation of a single crystal can be related to the primary slip system and the deformation induced activation of that system, which in turn can be related to the slip band formation of the single crystal sheet. Further, in case of a polycrystal sheet, the larger the grain size, the more the strain localizes at a specific crystal, which has the particular orientation. It is confirmed through comparisons with experiments that our finite element code can predict the localization of strain in sheets and consequently can estimate the formability of sheet metals.     

Finite element simulation of PSB macroband nucleation and propagation in single crystal nickel cycled at low plastic strain amplitudes

Substructure models for vein matrix and persistent slip band (PSB) structures are extracted from a uniaxial mixtures model that was developed to simulate cyclic loading experiments on nickel single crystals oriented for single slip. Reverse magnetostriction is included as well. These substructure models are implanted in a single crystal plasticity framework with fully anisotropic elasticity. The resulting constitutive models are incorporated in finite element models to simulate the process of PSB macroband formation and propagation. Perturbation elements (PEs), elements assigned with PSB properties, are used as the loci for PSB macroband nucleation. Transition of elements with vein matrix properties to elements with PSB properties is triggered at integration points by a shear stress criterion applied on slip systems. The resulting finite element models successfully demonstrate the process of PSB formation and propagation, and plastic strain amplitude partitioning between vein matrix and PSB macrobands. The effect of model boundary constraints, strain increment dependence, mesh sensitivity, PE distribution, specimen axis misorientation, and PSB volume fraction generated is examined.     

Studies of scale dependent crystal viscoplasticity models

A scale dependent crystal viscoplasticity model with a second strain gradient effect is introduced, as a simple extension of the conventional crystal plasticity theory. We confine attention to a single crystal undergoing slip on a single slip system under small strain conditions. Connections between this model and other existing theories are investigated in some detail. Furthermore, some basic predictions of the model, due to the second gradients and the material viscosity, are illustrated, using a constrained simple shear problem for a thin strip bounded by two rigid walls. The effect of viscosity on evolution of the boundary layer is examined, as well as the behavior of the thin strip undergoing reverse/cyclic shear loading, and the ability to predict plastic flow localization.     

Rheology of highly concentrated anionic surfactants

Surfactant solution flow behavior is of great importance to both the chemical and consumer product industries. Most studies on the flow behavior of surfactant solutions, however, have focused on the dilute regime. Seldom reported is rheology in the highly concentrated regime where typically these surfactants are processed and delivered. First, we present here the phase diagram for the ternary system: water and two anionic surfactants (sodium salt of lauric and oleic acid) at different temperatures. Then, we present both linear viscoelastic and steady shear flow results in the high (70 to 90%) surfactant regime. We find that high values of the shear modulus are directly dependent on the quantity of surfactant crystals and that the formation of a lamellar liquid crystal phase at 45°C affects both modulus and flow of the system. Lamellar crystals create a stiff network resulting in wall slip at large shear strain. Using serrated plates removes slip at the wall and we find a shear rate where microfractures localize in a preferential plane and the material flows. This behavior is reversible.     

A Dedicated DIC Methodology for Characterizing Plastic Deformation in Single Crystals

This paper focuses on the development of an appropriate Digital Image Correlation (DIC) methodology based on Image Registration and dedicated for characterizing the plastic deformation in single crystals. A pure nickel single crystal specimen is plastically deformed in tension and investigated by DIC technique. Based on the measured kinematic fields, the proposed method enables to identify the slip activity on the crystal surface and to locate precisely the slip band interfaces at microscale which behave as kinematic discontinuities. The computed displacement data are projected on a well-defined physical basis containing slip details, then the strain fields can be derived directly from a set of analytical functions. The possible errors in displacement induced by this projection approach are evaluated. Finally, some results of the evaluated strain fields are presented. It demonstrates that the developed DIC methodology allows quantitative characterization of a heterogeneous deformation process and promotes further relationships to be established between slip activity and strain field evolution in single crystals.     

基于有限变形单晶塑性滑移的微观热-力耦合模型及其有限元计算

基于率相关的晶体塑性滑移理论,论文考虑晶体内部塑性变形产生的热以及快速热冲击作用下温度急剧变化产生热应力的热-力双向耦合效应,建立起微观单晶的瞬态热-弹-塑性耦合模型,推导出与温度有关的剪应变率和弹塑性切模量公式.根据论文建立的模型,对ABAQUS软件进行二次开发[1],数值模拟出<001>/{100}单晶Cu在单轴拉伸状态下的应力、应变与温度的关系和弹性模量的变化,结果如下:轴向应力随温度升高先呈线性增加再呈非线性减小,轴向应变随温度增加而增加;弹性模量随塑性变形的增加而降低,与分子动力学模拟的趋势[2]是一致的.数值实验表明,论文建立的模型和算法是正确合理的,且计算量远远小于分子动力学模拟.     

Active slip-band separation and the energetics of slip in single crystals

This research supports recent efforts to provide an energetic approach to the prediction of stress–strain relations for single crystals undergoing single slip and to give precise formulations of experimentally observed connections between hardening of single crystals and separation of active slip-bands. Non–classical, structured deformations in the form of two-level shears permit the formulation of new measures of the active slip-band separation and of the number of lattice cells traversed during slip. A formula is obtained for the Helmholtz free energy per unit volume as a function of the shear without slip, the shear due to slip, and the relative separation of active slip-bands in a single crystal. This formula is the basis for a model, under preparation by the authors, of hardening of single crystals in single slip that is consistent with the Portevin-Le Chatelier effect and the existence of a critical resolved shear stress.     

Cylindrical void in a rigid-ideally plastic single crystal III: Hexagonal close-packed crystal

The analytical solution is derived for the plane strain stress field around a cylindrical void in a hexagonal close-packed single crystal with three in-plane slip systems oriented at the angle π/3 with respect to one another. The critical resolved shear stress on each slip system is assumed to be equal. The crystal is loaded by both internal pressure and a far-field equibiaxial compressive stress. The deformation field takes the form of angular sectors, called slip sectors, within which only one slip system is active; the boundaries between different sectors are radial lines. The stress fields are derived by enforcing equilibrium and a rigid, ideally plastic constitutive relationship, in the spirit of anisotropic slip line theory. The results show that each slip sector is divided into smaller regions denoted as stress sectors and the stress state valid within each stress sector is derived. It is shown that stresses are unique and are continuous within stress sectors and across stress sector boundaries, but the gradient of stresses is not continuous across the boundaries between stress sectors. The solution shows self-similarity in that the stresses over the entire domain can be determined from the stresses within a small region adjacent to the void by invoking certain scaling and symmetry properties. In addition, the stress state exhibits periodicity along logarithmic spirals which emanate from the void. The results predict that the mean value of in-plane pressure required to activate plastic deformation around a void in a single crystal can be higher than that necessary for a void in an isotropic material and is sensitive to the orientation of the slip systems relative to the void.     

A mesoscale investigation of strain rate effect on dynamic deformation of single-crystal copper

A combined finite element (FE) simulation and discrete dislocation dynamics (DD) approach has been developed in this paper to investigate the dynamic deformation of single-crystal copper at mesoscale. The DD code yields the plastic strain based on the slip of dislocations and serves as a substitute for the 3D constitutive form used in the usual FE computation, which is implemented into ABAQUS/Standard with a user-defined material subroutine. On the other hand, the FE code computes the displacement and stress field during the dynamic deformation. The loading rate effects on the yield stress and the deformation patterning of single-crystal copper are investigated. With the increasing of strain rate, the yield stress of single-crystal copper increases rapidly. A critical strain rate exists in each single-crystal copper block for the given size and dislocation sources, below which the yield stress is relatively insensitive to the strain rate. The dislocation patterning changes from non-uniform to uniform under high-strain-rate. The shear stresses in the bands are higher than that in the neighboring regions, which are formed shear bands in the crystal. The band width increases with the strain rate, which often take places where the damage occurs.     

Crystal plasticity analysis of dislocation emission from micro voids

Slip deformation in the vicinity of a micro void in metal crystals is analyzed by a crystal plasticity technique, and the geometrically necessary dislocations, which accompany the gradient of plastic shear strain on slip systems, are evaluated. Aggregates of dislocation segments on pairs of slip systems that have the same slip directions but different slip planes exhibit a rhombus-shaped structure, and the structure is shown to be equivalent to prismatic dislocation loops of the interstitial type. Material transport and growth of voids are discussed in terms of GN dislocations.     

Adiabatic shear band localization of inelastic single crystals in symmetric double-slip process

Summary The main objective of the present paper is the development of a viscoplastic regularization procedure valid for an adiabatic dynamic process for multi-slips of single crystals. The next objective is to focus attention on the investigation of instability criteria, and particularly on shear band localization conditions.To achieve this aim, an analysis of acceleration waves is given, and advantage is taken of the notion of the instantaneous adiabatic acoustic tensor. If zero is an eigenvalue of the acoustic tensor, then the associated discontinuity does not propagate, and one speaks of a stationary discontinuity. This situation is referred to as the strain localization condition, and corresponds to a loss of hyperbolicity of the dynamical equations. It has been proved that for an, adiabatic process of rate-dependent (elastic-viscoplastic) crystal, the wave speed of discontinuity surface always remains real and different from zero. It means that for this case the initial-value problem is well-posed. However, for an adiabatic process of rate-independent(elastic-plastic) crystal, the wave speed of discontinuity surface can be equal zero. Then the necessary condition for a localized plastic deformation along the shear band to be formed is as follows: the determinant of the instantaneous adiabatic acoustic tensor is equal to zero. This condition for localization is equivalent to that obtained by using the standard bifurcation method. Based on this idea, the conditions for adiabatic shear band localization of plastic deformation have been investigated for single crystals. Particular attention has been focused on the discussion of the influence of thermal expansion, thermal plastic, softening and spatial covariance effects on shear band localization criteria for a planar model of an f.c.c. crystal undergoing symmetric primary-conjugate double slip. The results obtained have been compared with available experimental observations.Finally, it is noteworthy that the viscoplasticity regularization procedure can be used in the developing of an unconditionally stable numerical integration algorithm for simulation of adiabatic inelastic flow processes in ductile single crystals, cf. [21].The paper has been prepared within research programme sponsored by the Committee of Scientific Research under Grant 3 P404 031 07.     

A peridynamic implementation of crystal plasticity

This paper presents the first application of peridynamics theory for crystal plasticity simulations. A state-based theory of peridynamics is used (Silling et al., 2007) where the forces in the bonds between particles are computed from stress tensors obtained from crystal plasticity. The stress tensor at a particle, in turn, is computed from strains calculated by tracking the motion of surrounding particles. We have developed a quasi-static implementation of the peridynamics theory. The code employs an implicit iterative solution procedure similar to a non-linear finite element implementation. Peridynamics results are compared with crystal plasticity finite element (CPFE) analysis for the problem of plane strain compression of a planar polycrystal. The stress, strain field distribution and the texture formation predicted by CPFE and peridynamics were found to compare well. One particular feature of peridynamics is its ability to model fine shear bands that occur naturally in deforming polycrystalline aggregates. Peridynamics simulations are used to study the origin and evolution of these shear bands as a function of strain and slip geometry.     

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.     

A STUDY OF THE ONSET OF KINK BAND

This paper suggests the asymmetric-plastic theory of crystalline solidsconcerning the plastic rotation of crystal.The deformation of crystalline solidsundergoes three microprocesses:lattice distortion,slip over active crystallographicplanes and the rotation of the crystal.In the asymmetric-plastic theory of crystalsuggested,the corotational rates of symmetric Cauchy stress and moment stresscorrespond to the rates of elastic strain and the gradient of the rotational rate of thecrystal respectively.The Schmid yield condition and Prager's consistency conditionincorporating antisymmetric stress are formulated.Then the asymmetric-plastic modelof crystalline solids is applied for the investigation of the onset of kink band by astandard stability analysis.The orientation of the kink band is perpendicular to theprimary slip system.The width of the kink band is the function of the“characteristiclength”of the microstructure of metal materials.     

Reverse shear on inclined planes at the tip of a sharp crack

Plane strain plastic yielding at a crack tip has been represented by edge dislocations with Burgers vectors parallel to symmetrical planes inclined at 70° and 45° to the plane of the crack. The plastic displacement and the stresses near the crack tip were calculated by a numerical method and the effect of a reduction in applied stress was determined. Removal of the whole or a part of the initial load produces reverse shear in regions of the slip band nearest the crack tip. The amount of reverse shear depends only on the reduction in the load and not on its initial value. The reverse shear is associated with the presence of negative dislocations and the stresses near the crack tip may become compressive even though the applied (remote) stress is still tensile. The degree and extent of compression depends on the reduction in applied stress and on its original value. It is argued that the residual compressive stresses produced under fluctuating loads may produce crack closure and crack arrest. The effect of residual plasticity in a slip band left behind a growing crack has been estimated. It is shown that after an overload the excess residual plasticity opposing crack opening rises to a maximum value when the crack tip has advanced some distance from the point where the overload was applied.     

晶体塑性变形离散滑移模型及有限元分析

基于韧性单晶体实验现象，建立了描述晶体塑性变形的离散滑移模型．该模型的主要特点是：晶体滑移变形在宏观上是不均匀的，滑移带的分布是离散的．利用晶体塑性理论对模型进行了有限变形有限元分析，计算结果揭示了晶体滑移的离散行为，模拟的应力 应变曲线与实验曲线相吻合     

部分剪力连接钢—砼组合梁的滑移及曲率分析

对受均布荷载作用的简支部分剪力连接钢－砼组合梁的滑移应变和滑移的分布规律进行了理论分析,得出部分剪力连接钢－砼组合梁的剪力连接程度系数与滑移应变和滑移的关系。分析了考虑滑移效应的梁的曲率及其与剪力连接程度系数之间的关系。结果与结构在弹性阶段的实验相当吻合,对组合梁挠度的精确计算有指导意义。