A finite strain elastic-viscoplastic self-consistent model for polycrystalline materials

A large strain elastic-viscoplastic self-consistent (EVPSC) model for polycrystalline materials is developed. At single crystal level, both the rate sensitive slip and twinning are included as the plastic deformation mechanisms, while elastic anisotropy is accounted for in the elastic moduli. The transition from single crystal plasticity to polycrystal plasticity is based on a completely self-consistent approach. It is shown that the differences in the predicted stress-strain curves and texture evolutions based on the EVPSC and the viscoplastic self-consistent (VPSC) model proposed by Lebensohn and Tomé (1993) are negligible at large strains for monotonic loadings. For the deformations involving unloading and strain path changes, the EVPSC predicts a smooth elasto-plastic transition, while the VPSC model gives a discontinuous response due to lack of elastic deformation. It is also demonstrated that the EVPSC model can capture some important experimental features which cannot be simulated by using the VPSC model.     

Effect of secondary orientation on notch-tip plasticity in superalloy single crystals

Numerical and experimental evolutions of slip fields in notched Ni-Base Single Crystal superalloy tensile specimens are presented as a function of secondary crystallographic orientation. The numerical predictions based on three-dimensional anisotropic elasticity and crystal plasticity are compared with experimental observations. The results illustrate the strong dependence of the slip patterns and the plastic zone size and shape on the secondary orientation of notches, which can have important consequences on crack initiation. Specific orientations or non-symmetric notch geometries lead to non-symmetric patterns on both sides of the sample. The computations show that strongly different plastic zones are expected in the core of the sample and at free surfaces. The ability of the anisotropic elastic model to anticipate the plastic domains, based on identifying dominant slip systems, is confirmed by the crystal plasticity computations, at low load levels. An important observation is that kink shear banding is a real deformation mode operating at crack tips and notches in high strength nickel-based single crystal superalloys for specific orientations.     

Numerical simulation of shape memory alloys structures using interior-point methods

This paper is concerned with the numerical simulation of quasi-static evolutions in shape memory alloys structures. In micromechanical models of shape memory alloys, the phase transformation is described by an internal variable representing the volume fractions of the different variants of martensite. In contrast with the framework of standard plasticity, that internal variable is physically constrained to satisfy a set of inequalities at each point. From a structural simulation point of view, handling such local constraints in the evolution problem is not obvious. Robustness and computational time are particularly sensitive issues. In this paper is proposed a new approach which essentially consists in reformulating the non-linear incremental problem as a linear complementarity problem. That approach notably calls for an easy implementation in FEM codes. Comparison with some experimental and numerical results from literature shows the relevancy of the proposed approach.     

孪生诱发塑性对V-5Cr-5Ti合金应变硬化行为的影响

钒合金(V-Cr-Ti)作为潜在重要的聚变反应堆用结构材料, 近年来受到广泛的关注. 为了研究 V-5Cr-5Ti 合金不同应变率压缩下的应变硬化行为, 特别是孪生对塑性变形的影响, 以位错密度和孪晶演化为基础, 建立了该合金的应变硬化模型. 模型中考虑了孪晶中的位错滑移对材料塑性应变的贡献. 模拟结果表明, 由于孪生诱发塑性, 从而使动态压缩时的位错密度小于准静态加载时的, 这使得 V-5Cr-5Ti 合金在动态压缩时的应变硬化率比准静态加载时的小. 当孪晶形成后, 位错滑移引起的塑性应变率随应变增大而增大, 并逐渐接近加载应变率, 而孪生引起的塑性应变率则随应变增大而减小.      

Modeling of deformation and rotation bands and of deformation induced grain boundaries in IF steel aggregate during large plane strain compression

A computation using crystal plasticity modeling of an actual IF steel aggregate plane strain compression deformation, underlines the formation of different deformation bands morphologies and grain splitting occurrence, already experimentally observed by different authors. The model based on dislocation densities as internal variables, developed in the framework of finite deformation and implemented in the Finite Element Method, is able to capture the main characteristics of different inhomogeneities and to analyze their formation and further development with strain, from the determination of the active and latent slip systems, and also from the quantification of their dislocation densities and corresponding glide rates evolutions. The respective boundary conditions and material properties effects are discussed.     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

Non-local crystal plasticity model with intrinsic SSD and GND effects

A strain gradient-dependent crystal plasticity approach is presented to model the constitutive behaviour of polycrystal FCC metals under large plastic deformation. In order to be capable of predicting scale dependence, the heterogeneous deformation-induced evolution and distribution of geometrically necessary dislocations (GNDs) are incorporated into the phenomenological continuum theory of crystal plasticity. Consequently, the resulting boundary value problem accommodates, in addition to the ordinary stress equilibrium condition, a condition which sets the additional nodal degrees of freedom, the edge and screw GND densities, proportional (in a weak sense) to the gradients of crystalline slip. Next to this direct coupling between microstructural dislocation evolutions and macroscopic gradients of plastic slip, another characteristic of the presented crystal plasticity model is the incorporation of the GND-effect, which leads to an essentially different constitutive behaviour than the statistically stored dislocation (SSD) densities. The GNDs, by their geometrical nature of locally similar signs, are expected to influence the plastic flow through a non-local back-stress measure, counteracting the resolved shear stress on the slip systems in the undeformed situation and providing a kinematic hardening contribution. Furthermore, the interactions between both SSD and GND densities are subject to the formation of slip system obstacle densities and accompanying hardening, accountable for slip resistance. As an example problem and without loss of generality, the model is applied to predict the formation of boundary layers and the accompanying size effect of a constrained strip under simple shear deformation, for symmetric double-slip conditions.     

2D vortex interaction in a non-uniform flow

In a two-dimensional incompressible fluid, we study the interaction of two like-signed Rankine vortices embedded in a steady shear/strain flow. The numerical results of vortex evolutions are compared with the analytical results for point vortices. We show the existence of vortex equilibria, and of merger for initial distances larger than those without external flow. The evolutions depend on the initial orientation of the vortices in the external flow.     

Non-linear and linear evolution of perturbation in stochastic basic flows

In this paper, we examine the non-linear and linear evolutions of perturbation in stochastic basic flows with two-dimensional quasi-geostrophic equations on a sphere. As the analytic solutions for the considered quasi-geostrophic equations are not available, the Fourier finite volume element method is used to perform numerical simulation. It is found that, the non-linear and linear evolutions of perturbation in stochastic basic flow will be consistent for a short period of time and small stochastic fluctuations when they are consistent in the deterministic basic flow. However, the tangent linear model will fail to approximate the original non-linear model when the time period is considerably long and stochastic fluctuation becomes large. Moreover, the global energy decays faster for stochastic basic flow with stronger fluctuations.     

Texture gradient simulations for extrusion and reversible rolling of FCC metals

A texture simulation method is described for some complex plane strain deformation paths during hot shaping of FCC metals. The method employs both finite element calculations and a polycrystal plasticity model based on the Relaxed-Constraints (RC) Taylor hypothesis and a viscoplastic constitutive law. We have considered the {111}<110> slip systems and the {100}, {110}, {112} <110> non-octahedral slip systems. The finite element codes simulate the strain paths of material flow during a shaping process. The local velocity gradients, expressed in the macroscopic reference coordinates, are rewritten in the local flow line coordinates using a kinematic analysis for steady-state flow. Secondly, for the different deformation paths, the RC polycrystal plasticity model is used to numerically simulate the local deformation texture evolutions as a function of depth. Texture simulations are carried out for two deformation processes combining hot compression and shear: extrusion and reversible rolling. For extrusion, the simulated pole figures and ODFs show the typical texture variations through the thickness of an extruded 6082 aluminium alloys, i.e. (β-fibre in the centre and a TD rotated copper component near the surface. It is shown that hot reversible rolling should develop a strong pure shear texture {001}<110> near the surface.     

基于细观拓扑结构演化的颗粒材料剪胀性分析

剪胀性是包括岩土材料在内的摩擦性颗粒材料的重要特征之一,其形成机制与颗粒体系内部拓扑结构的演化有关.基于颗粒体系细观数据,可对颗粒体系内部的拓扑结构特征及演化进行分析,进而建立拓扑演化与宏观剪胀变形之间的联系.采用离散单元法,根据密实、中密和松散摩擦性颗粒材料双轴试验的宏微观数据,从拓扑参量演化及接触网络拓扑变化所引起...     

Strain gradient plasticity by internal-variable approach with normality structure

This paper presents a constitutive formulation for materials with strain gradient effects by internal-variable approach with normality structure. Specific micro-structural rearrangements are assumed to account for the inelasticity deformations for this class of materials, and enter the constitutive formulations in form of internal variables. It is further assumed that the kinetic evolution of any specific micro-structural rearrangement may be fully determined by the thermodynamic forces associated with that micro-structural rearrangement, by normality relations via a flow potential. Macroscopic gradient-enhanced inelastic behaviours may then be predicted in terms of the microscopic internal variables and their conjugate forces, and thus a micro–macro bridging formulation is available for strain-gradient-characterised materials. The obtained formulations are first applied to crystallographic materials, and a crystal gradient plasticity model is developed to account for the influence of microscopic slip rearrangements on the macroscopic gradient-dependent mechanical behaviour for this class of materials. Micro-cracked geomaterials are also treated with these formulations and a gradient-enhanced damage constitutive model is developed to address the impacts of the evolutions of micro-cracks on the macroscopic inelastic deformations with strain gradient effects for these materials. The available formulations are further compared with other thermodynamic approaches of constitutive developing.     

Generalized plasticity and some models for geomechanics

In this brief note, we have (a) introduced a very general form of plasticity definition; (b) indicated several types of models used for monotonic and transient loading of soil-like materials, which do not depend on time effects. The details of each form are given elsewhere but we believe that the new framework provides an easier and more general interpretation of a variety of behaviour form.     

Computationally efficient database and spectral interpolation for fully plastic Taylor-type crystal plasticity calculations of face-centered cubic polycrystals

A new computationally efficient database approach to fully plastic Taylor-type crystal plasticity calculations is presented in this paper. In particular, we explore strategies that circumvent the need to repeatedly solve sets of highly non-linear, extremely stiff, algebraic equations with poor convergence characteristics that are inherent to these calculations. The suggested strategies consist of computing only once all of the needed variables in crystal plasticity calculations, storing them, and retrieving the values of interest according to the need in any specific simulation. An algorithm is presented here that facilitates this approach, and involves local spectral interpolation using discrete fourier transform (DFT) methods. The approach described here results in major improvements in the computational time over both the conventional crystal plasticity calculations and our previously developed spectral approach using generalized spherical harmonics (GSH). Details of this new approach are described and validated in this paper through a few example case studies.     

解Drucker-Prager塑性问题的二阶锥互补法

基于经典弹塑性理论中多数屈服准则具有凸锥数学结构的事实,将在大规模计算中更具潜力的锥规划法引入弹塑性分析。考虑到弹塑性流动理论有关联与非关联之分,本文提出利用锥型互补法求解弹塑性问题。具体以Drucker-Prager弹塑性模型为例,首先利用最大塑性功耗散原理和圆锥对偶理论等工具,建立了弹塑性本构方程的等价二阶锥互补模型;然后,基于参变量变分原理和有限元技术,建立了弹塑性增量分析的二阶锥线性互补模型;最后,利用一类半光滑Newton算法求解。数值算例表明了本文方法的有效性。     

An electromechanical atomic-scale finite element method for simulating evolutions of ferroelectric nanodomains

In this paper, a novel atomic-level computational method of perovskite ferroelectrics is established by combining the shell model and atomic-scale finite element method (AFEM). Its applicability is carefully testified for both bulk and nanoscale ferroelectrics, by comparing the calculated structural parameters and polarizations with the molecular dynamics (MD) simulations, first-principles calculations and experiment results. A comparison of the CPU time demonstrates that the developed method has a computational speed about 10 times over that of shell model MD method and its advantage becomes more evident as the computational scale becomes larger. Moreover, two effective calculation skills of long-range Coulomb force are introduced which can further enhance the computational efficiency by about 10 times. Using the developed atomic-level method, we investigate the various patterns of nanoscale domain structures in BaTiO₃ and their evolutions under electrical loadings. A domain structure with coexistence of vortex and streamline polarization patterns is revealed. Furthermore, the simulations of domain evolutions not only reproduce well the two-step 90° domain switching process observed in experiments on single domain under an anti-parallel electric field, but also provide a full evolution diagram among different domain patterns under various electric fields. A quantitative analysis indicates that the direction-dependent coercive field of multi-domain structure can be well described by that of single domain based on a simple analytical model. This study on domain patterns and evolutions may help us understand the behaviors of ferroelectrics from the atomic level.     

An effective computational algorithm for rate-independent crystal plasticity based on a single crystal yield surface with an application to tube hydroforming

Up to now, several computational methods have been proposed for crystal plasticity models. The main objective of these computational methods has been to overcome the problem with the non-uniqueness of active slip systems during the plastic deformation of a single crystal. Crystal plasticity models based on a single crystal yield function have been proposed as alternative algorithms to overcome this problem. But the problem with these models is that they use a highly non-linear yield function for the crystal, which makes them computationally expensive. In this paper, a computational method is proposed that would modify a single crystal yield function in order to make it computationally efficient. Also to better capture experimental data, a new parameter is introduced into the single crystal yield function to make it more flexible. For verification, this crystal plasticity model was directly applied for the simulation of hydroforming of an extruded aluminum tube under complex strain paths. It was found that the current model is considerably faster than the previous crystal plasticity model based on a power-law type single crystal yield surface. Due to its computational efficiency, the current crystal plasticity model can also be used to calculate the anisotropy coefficients of phenomenological yield functions.     

Constitutive model for uniaxial transformation ratchetting of super-elastic NiTi shape memory alloy at room temperature

The transformation ratchetting of super-elastic NiTi shape memory alloy was observed by the uniaxial stress-controlled cyclic tests [Kang, G.Z., Kan, Q.H., Qian, L.M., Liu, Y.J, 2009a. Ratchetting deformation of super-elastic and shape memory NiTi Alloys. Mech. Mater. 41, 139–153]. It is concluded that the NiTi alloy presents apparent ratchetting behaviour, and the ratchetting is collectively caused by the cyclic accumulation of residual induced-martensite and the transformation-induced plastic deformation (i.e., namely transformation ratchetting). Based on the experimental results, a cyclic constitutive model was constructed in the framework of generalized plasticity [Lubliner, J., Auricchio, F., 1996. Generalized plasticity and shape memory alloys. Int. J. Solids Struct. 33, 991–1003] to describe the transformation ratchetting of super-elastic NiTi alloy. The proposed model simultaneously accounts for the evolutions of residual induced-martensite and transformation-induced plastic strain during the stress-controlled cyclic loading by introducing an internal variable z_c, i.e., cumulated induced-martensite volume fraction. The dependence of transformation ratchetting on the applied stress levels and the phase transformation hardening behaviour of the NiTi alloy are also considered in the developed model. The anisotropic phase transformation behaviours of the alloy presented in the tension and compression cases are described by employing a Drucker–Prager-typed transformation surface. It is shown that the simulated results of transformation ratchetting obtained by the proposed model are in good agreement with the corresponding experiments, since the typical features of transformation ratchetting are reasonably captured by the proposed model.