On the modelling of anisotropic material behaviour in viscoplasticity

A uniaxial viscoplastic deformation is motivated as a discrete sequence of stable and unstable equilibrium states and approximated by a smooth family of stable states of equilibrium depending on the history of the mechanical process. Three-dimensional crystal viscoplasticity starts from the assumption that inelastic shearings take place on slip systems, which are known from the particular geometric structure of the crystal. A constitutive model for the behaviour of a single crystal is developed, based on a free energy, which decomposes into an elastic and an inelastic part. The elastic part, the isothermal strain energy, depends on the elastic Green strain and allows for the initial anisotropy, known from the special type of the crystal lattice. Additionally, the strain energy function contains an orthogonal tensor-valued internal variable representing the orientation of the anisotropy axes. This orientation develops according to an evolution equation, which satisfies the postulate of full invariance in the sense that it is an observer-invariant relation. The inelastic part of the free energy is a quadratic function of the integrated shear rates and corresponding internal variables being equivalent to backstresses in order to consider kinematic hardening phenomena on the slip system level. The evolution equations for the shears, backstresses and crystallographic orientations are thermomechanically consistent in the sense that they are compatible with the entropy inequality. While the general theory applies to all types of lattices, specific test calculations refer to cubic symmetry (fcc) and small elastic strains. The simulations of simple tension and compression processes of a single crystal illustrates the development of the crystallographic axes according to the proposed evolution equation. In order to simulate the behaviour of a polycrystal the initial orientations of the anisotropy axes are assumed to be space-dependent but piecewise constant, where each region of a constant orientation corresponds to a grain. The results of the calculation show that the initially isotropic distribution of the orientation changes in a physically reasonable manner and that the intensity of this process-induced texture depends on the specific choice of the material constants.     

Integration of self-consistent polycrystal plasticity with dislocation density based hardening laws within an implicit finite element framework: Application to low-symmetry metals

We present an implementation of the viscoplastic self-consistent (VPSC) polycrystalline model in an implicit finite element (FE) framework, which accounts for a dislocation-based hardening law for multiple slip and twinning modes at the micro-scale grain level. The model is applied to simulate the macro-scale mechanical response of a highly anisotropic low-symmetry (orthorhombic) crystal structure. In this approach, a finite element integration point represents a polycrystalline material point and the meso-scale mechanical response is obtained by the mean-field VPSC homogenization scheme. We demonstrate the accuracy of the FE-VPSC model by analyzing the mechanical response and microstructure evolution of α-uranium samples under simple compression/tension and four-point bending tests. Predictions of the FE-VPSC simulations compare favorably with experimental measurements of geometrical changes and microstructure evolution. Specifically, the model captures accurately the tension–compression asymmetry of the material associated with twinning, as well as the rigidity of the material response along the hard-to-deform crystallographic orientations.     

On the coupling of plastic slip and deformation-induced twinning in magnesium: A variationally consistent approach based on energy minimization

The present paper is concerned with the analysis of the deformation systems in single crystal magnesium at the micro-scale and with the resulting texture evolution in a polycrystal representing the macroscopic mechanical response. For that purpose, a variationally consistent approach based on energy minimization is proposed. It is suitable for the modeling of crystal plasticity at finite strains including the phase transition associated with deformation-induced twinning. The method relies strongly on the variational structure of crystal plasticity theory, i.e., an incremental minimization principle can be derived which allows to determine the unknown slip rates by computing the stationarity conditions of a (pseudo) potential. Phase transition associated with twinning is modeled in a similar fashion. More precisely, a solid-solid phase transition corresponding to twinning is assumed, if this is energetically favorable. Mathematically speaking, the aforementioned transition can be interpreted as a certain rank-one convexification. Since such a scheme is computationally very expensive and thus, it cannot be applied to the analysis of a polycrystal, a computationally more efficient approximation is elaborated. Within this approximation, the deformation induced by twinning is decomposed into the reorientation of the crystal lattice and simple shear. The latter is assumed to be governed by means of a standard Schmid-type plasticity law (pseudo-dislocation), while the reorientation of the crystal lattice is considered, when the respective plastic shear strain reaches a certain threshold value. The underlying idea is in line with experimental observations, where dislocation slip within the twinned domain is most frequently seen, if the twin laminate reaches a critical volume. The resulting model predicts a stress-strain response in good agreement with that of a rank-one convexification method, while showing the same numerical efficiency as a classical Taylor-type approximation. Consequently, it combines the advantages of both limiting cases. The model is calibrated for single crystal magnesium by means of the channel die test and finally applied to the analysis of texture evolution in a polycrystal. Comparisons of the predicted numerical results to their experimental counterparts show that the novel model is able to capture the characteristic mechanical response of magnesium very well.     

A combined experimental-numerical approach for elasto-plastic fracture of individual grain boundaries

The parameters for a crystal plasticity finite element constitutive law were calibrated for the aluminum–lithium alloy 2198 using micro-column compression testing on single crystalline volumes. The calibrated material model was applied to simulations of micro-cantilever deflection tests designed for micro-fracture experiments on single grain boundaries. It was shown that the load–displacement response and the local deformation of the grains, which was measured by digital image correlation, were predicted by the simulations. The fracture properties of individual grain boundaries were then determined in terms of a traction–separation-law associated with a cohesive zone. This combination of experiments and crystal plasticity finite element simulations allows the investigation of the fracture behavior of individual grain boundaries in plastically deforming metals.     

The viscoplastic behavior of hastelloy-X single crystal

A viscoplastic constitutive model for Hastelloy-X single crystal material is developed based on crystallographic slip theory. The constitutive model was constructed for use in a viscoplastic self-consistent model for isotropic Hastelloy-X polycrystalline material, which has been described in a recent publication. It is found that, by using the slip geometry known from the metallurgical literature, the anisotropic response can be accurately predicted. The model was verified by using tension and torsion data taken at 982°C (1800°F). The constitutive model used on each slip system is a simple unified visoplastic power law model in which weak latent interaction effects are taken into account. The drag stress evolution equations for the octahedral system are written in a hardening/recovery format in which both hardening and recovery depend on separate latent interaction effects between the octahedral crystallographic slip systems. The strain rate behavior of the single crystal material is well correlated by the constitutive model in uniaxial and torsion tests, but it is necessary to include latent information effects between the octahedral slip systems in order to obtain the best possible representation of biaxial cyclic strain rate behavior. Finally, it was observed that the single crystal exhibited dynamic strain aging at 871°C (1600°F). Similar dynamic strain aging occurs at 649°C (1200°F) in the polycrystalline version of the alloy.     

铁电单晶Pb(Mg1/3Nb2/3)O3-0.32PbTiO3的实验及本构模型研究

对压力作用下沿[001]晶向极化的Pb(Mg弛豫型铁电单晶;极化旋转(相变);黏塑性模型;本构;细观力学国家自然科学基金,教育部全国优秀博士学位论文作者专项基金2006-11-06对压力作用下沿[001]晶向极化的Pb(Mg1/3Nb2/3)O3-0.32PbTiO3(PMN0.32PT)弛豫型铁电单晶的应力应变行为进行了实验研究,实验结果表明铁电单晶〈001〉晶向的应力应变行为和铁电多晶有本质的不同,是传统的电畴翻转机理所难以解释的,所提出的极化旋转(相变)模型合理地解释了实验中观察到的现象;基于提出的极化旋转(相变)模型,采用细观力学方法建立了铁电单晶的细观本构模型.在模型中采用黏塑性公式描述铁电单晶可能的8个相变系统的相变行为.为了验证模型的可靠性,用该模型模拟了铁电单晶〈001〉晶向的应力应变实验曲线.计算表明,该模型能较好地模拟铁电单晶〈001〉晶向的相变行为.     

Forming limit comparisons for FCC and BCC sheets

In this paper, numerical simulations of forming limit diagrams (FLDs) are performed based on a rate-sensitive polycrystal plasticity model together with the Marciniak–Kuczynski (M–K) approach. Sheet necking is initiated from an initial imperfection in terms of a narrow band. The deformations inside and outside the band are assumed to be homogeneous, and conditions of compatibility and equilibrium are enforced across the band interfaces. Thus, the polycrystal model need only be applied to two polycrystal aggregates, one inside and one outside the band. Both FCC and BCC crystals are considered with 12 distinct slip systems for an FCC crystal and 24 distinct slip systems for a BCC crystal. The response of an aggregate comprised of many grains is based on an elastic–viscoplastic Taylor-type polycrystal model. With this formulation, the effects of initial imperfection intensity and orientation, crystal elasticity, strain-rate sensitivity, single slip hardening, and latent hardening on the FLD can be assessed. Identical initial textures are considered for both FCC and BCC polycrystals and the predicted FLDs are compared with each other.     

Finite element modeling of crystal plasticity with grains shaped as truncated octahedrons

Different modeling strategies are tested for the prediction of texture development and microscopic strain heterogeneity in cold-rolled ULC steel and in multiphase steel under uniaxial tension. The polycrystalline aggregate is represented by a finite element mesh that is loaded under periodic boundary conditions. Grains are shaped as cubes or as truncated octahedrons, defining three levels of mesh refinement. Simulations rely on a simplified implementation of crystal plasticity, in which elastic strains are considered infinitesimal. The constitutive law is integrated fully implicitly in a reference frame bounded to the crystal lattice. Accuracy of the time-integration procedure is assessed by referring to some predictions of the crystal plasticity algorithm developed by Kalidindi et al. [Crystallographic texture evolution in bulk deformation processing of FCC metals. J. Mech. Phys. Solid 40 (1992) 537–569]. Then, results of the micro–macro modeling are compared to experimental data. It is found that the simulations with truncated octahedral grains yield improved predictions compared to those with cuboidal grains.     

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.     

Incorporation of twinning into a crystal plasticity finite element model: Evolution of lattice strains and texture in Zircaloy-2

A crystal plasticity finite element code is developed to model lattice strains and texture evolution of HCP crystals. The code is implemented to model elastic and plastic deformation considering slip and twinning based plastic deformation. The model accounts for twinning reorientation and growth. Twinning, as well as slip, is considered to follow a rate dependent formulation. The results of the simulations are compared to previously published in situ neutron diffraction data. Experimental results of the evolution of the texture and lattice strains under uniaxial tension/compression loading along the rolling, transverse, and normal direction of a piece of rolled Zircaloy-2 are compared with model predictions. The rate dependent formulation introduced is capable of correctly capturing the influence of slip and twinning deformation on lattice strains as well as texture evolution.     

Microstructural effects on yield surface evolution in cubic metals using the viscoplastic ?-model

Polycrystalline yield surfaces of metals are a good way to characterize the anisotropy of plastic deformation. The evolution of these surfaces is impossible to accurately reproduce without taking into account the evolution of the material microstructure such as texture development. In this paper, a numerical computation of yield surfaces using the viscoplastic ?-model is proposed. Results concerning face-centered cubic metals subjected to a plane strain compression test are presented. The influence of several mechanical parameters (strain hardening, strain rate sensitivity coefficient and accumulated deformation) on subsequent yield surfaces evolution is studied. The analysis of the change in the shape and size of the yield surfaces shows that the results depend strongly on the parameter ? which controls the strength of the interactions in the polycrystal. In addition, the predictions are compared to the widely used viscoplastic self-consistent model as well as to experimental yield loci taken from the literature for various aluminum alloy sheets. A fairly good qualitative agreement between our ?-model results and the experimental ones is found. The probable links between the parameter ? and the microstructural features such as the stacking fault energy and the grain size of the polycrystal are also briefly discussed.     

Evolution of plastic zones in dynamically loaded plates using different elastic–viscoplastic laws

In the framework of viscoplastic theory many different laws were developed, accounting for material behaviors like creep, relaxation or evolution of overstresses. Though each model is able to predict in uni-axial material tests the values of stresses depending on plastic strains and plastic strain rates the question is if solutions of simulations are still realistic if the viscoplastic law is applied on structural deformations. In the present study strain rate sensitive metal plates are subjected to shock waves. The purpose is to compare simulation results obtained with different elastic–viscoplastic laws to experiments in order to determine the most appropriate material model. By subjecting circular metal plates experimentally to shock wave loadings realistic deformation histories are measured. The measurements are compared to simulation results obtained with different viscoplastic laws. The aim is to find out the accuracy of each model concerning the predictions of displacements, shape formings, spread of plastic zones and evolutions of inner bending moments.     

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.     

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.     

单晶有限变形的热力耦合计算模型

以弹性变形梯度作为基本变量,结合热力学理论构造了单晶有限变形的热、力耦合计算模型。该模型考虑了温度、变温速率以及塑性耗散等条件对单晶有限变形的影响,相对于传统的以弹性变形梯度为基本变量的晶体塑性模型,算法能够体现温度效应的影响。采用隐式的积分方法对建立的控制方程进行计算以保证求解过程的稳定。以1100Al单晶为例计算了不同升温、降温速率,以及不同应变率影响下的材料应力-应变的响应。结果表明,模型能较好地反映变温过程中,单晶各向异性性质的演化以及应力、应变之间关系的变化。     

CuAlNi单晶形状记忆合金压缩特性的实验研究

具有相变伪弹性特性的CuAlNi单晶是目前应用最广泛的形状记忆合金之一。这种材料被广泛应用在工程、生物和医学科学等领域。由于单晶是各向异性，没有多晶中晶粒之间的相互作用，因而在特定晶向上的力学性能稳定。但是这种材料的一些基本性质，如压缩状态下马氏体的发生、生长和传播等还没有人详细研究。本文主要研究CuAlNi单晶在特定晶向上的变形过程，它的应力—应变特性，并利用特殊的显微成像系统首次获得了沿[110]方向压缩时二维马氏体的发生、生长和传播过程。     

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.     

Crystal plasticity model with enhanced hardening by geometrically necessary dislocation accumulation

A strain gradient dependent crystal plasticity approach is used to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. Material points are considered as aggregates of grains, subdivided into several fictitious grain fractions: a single crystal volume element stands for the grain interior whereas grain boundaries are represented by bi-crystal volume elements, each having the crystallographic lattice orientations of its adjacent crystals. A relaxed Taylor-like interaction law is used for the transition from the local to the global scale. It is relaxed with respect to the bi-crystals, providing compatibility and stress equilibrium at their internal interface. During loading, the bi-crystal boundaries deform dissimilar to the associated grain interior. Arising from this heterogeneity, a geometrically necessary dislocation (GND) density can be computed, which is required to restore compatibility of the crystallographic lattice. This effect provides a physically based method to account for the additional hardening as introduced by the GNDs, the magnitude of which is related to the grain size. Hence, a scale-dependent response is obtained, for which the numerical simulations predict a mechanical behaviour corresponding to the Hall-Petch effect. Compared to a full-scale finite element model reported in the literature, the present polycrystalline crystal plasticity model is of equal quality yet much more efficient from a computational point of view for simulating uniaxial tension experiments with various grain sizes.