叠加型A-F类随动强化模型塑性应变的数值计算法

以Chaboche随动强化模型为例,在M isses屈服准则及正交流动准则的前提下,推导了叠加型A rm-strong-F rederick(A-F)类随动强化模型塑性应变的数值计算法,联合利用四阶龙格-库塔法与径向返回法实现数值计算中的内部平衡迭代。同时推导了统一切向矩阵以便确定每一平衡迭代后的试算应变。利用AN SY S提供的U PF s将算法嵌入到AN SY S有限元程序,实现了叠加型A-F类随动强化模型塑性应变的数值计算,并利用四边形单元模拟了单轴循环加载时的棘轮应变,计算结果能够很好地与实验值吻合。     

Meshless analysis of plasticity with application to crack growth problems

The governing equations for classical rate-independent plasticity are formulated in the framework of meshless method. The special J₂ flow theory for three-dimensional, two-dimensional plane strain and plane stress problems are presented. The numerical procedures, including return mapping algorithm, to obtain the solutions of boundary-value problems in computational plasticity are outlined. For meshless analysis the special treatment of the presence of barriers and mirror symmetries is formulated. The crack growth process in elastic–plastic solid under plane strain and plane stress conditions is analyzed. Numerical results are presented and discussed.     

An implicit consistent algorithm for the integration of thermoviscoplastic constitutive equations in adiabatic conditions and finite deformations

The so-called viscoplastic consistency model, proposed by Wang, Sluys and de Borst, is extended here to the integration of a thermoviscoplastic constitutive equation for J₂ plasticity and adiabatic conditions. The consistency condition in this case includes not only strain rate but also the effect of temperature on the yield function. Using the backward Euler integration scheme to integrate the constitutive equations, an implicit algorithm is proposed, leading to generalized expressions of the classical return mapping algorithm for J₂ plasticity, both for the iterative calculation of the plastic multiplier increment and for the consistent tangent operator when strain rate and temperature are considered also as state variables of the hardening equation. The model was implemented in a commercial finite element code and its performance is demonstrated with the numerical simulation of four Taylor impact tests.     

Material force in micromorphic plasticity

Micromorphic theory envisions a material body as a continuous collection of deformable particles with finite size and inner structure. It is considered as the most successful top-down formulation of a two-level continuum model, in which the deformation is expressed as a sum of macroscopic continuous deformation and internal microscopic deformation of the inner structure. In this work, we revisit the original micromorphic theory and further construct a mathematical theory of micromorphic plasticity with generalized strain-based return mapping algorithm. The concept of material forces, which may also be referred as Eshelbian mechanics, was first derived for micromorphic thermo-visco-elastic solid, and, now in this work, it is extended to the micromorphic plasticity. The balance law of pseudo-momentum is formulated. The detailed expressions of Eshelby stress tensor, pseudo-momentum, and material forces are derived. Following this formulation, the failure mechanisms of micromorphic thermo-visco-elastoplastic materials can be further investigated.     

A multi-field incremental variational framework for gradient-extended standard dissipative solids

The paper presents a constitutive framework for solids with dissipative micro-structures based on compact variational statements. It develops incremental minimization and saddle point principles for a class of gradient-type dissipative materials which incorporate micro-structural fields (micro-displacements, order parameters, or generalized internal variables), whose gradients enter the energy storage and dissipation functions. In contrast to classical local continuum approaches to inelastic solids based on locally evolving internal variables, these global micro-structural fields are governed by additional balance equations including micro-structural boundary conditions. They describe changes of the substructure of the material which evolve relatively to the material as a whole. Typical examples are theories of phase field evolution, gradient damage, or strain gradient plasticity. Such models incorporate non-local effects based on length scales, which reflect properties of the material micro-structure. We outline a unified framework for the broad class of first-order gradient-type standard dissipative solids. Particular emphasis is put on alternative multi-field representations, where both the microstructural variable itself as well as its dual driving force are present. These three-field settings are suitable for models with threshold- or yield-functions formulated in the space of the driving forces. It is shown that the coupled macro- and micro-balances follow in a natural way as the Euler equations of minimization and saddle point principles, which are based on properly defined incremental potentials. These multi-field potential functionals are outlined in both a continuous rate formulation and a time-space-discrete incremental setting. The inherent symmetry of the proposed multi-field formulations is an attractive feature with regard to their numerical implementation. The unified character of the framework is demonstrated by a spectrum of model problems, which covers phase field models and formulations of gradient damage and plasticity.     

A study of microbend test by strain gradient plasticity

Metallic materials display strong size effect when the characteristic length associated with plastic deformation is on the order of microns. This size effect cannot be explained by classical plasticity theories since their constitutive relations do not have an intrinsic material length. Strain gradient plasticity has been developed to extend continuum plasticity to the micron or submicron regime. One major issue in strain gradient plasticity is the determination of the intrinsic material length that scales with strain gradients, and several microbend test specimens have been designed for this purpose. We have studied different microbend test specimens using the theory of strain gradient plasticity. The pure bending specimen, cantilever beam, and the microbend test specimen developed by Stolken and Evans (Stolken, J.S., Evans, A.G., 1998. A microbend test method for measuring the plasticity length scale Acta Mater. 46, 5109–5115) are found suitable for the determination of intrinsic material length in strain gradient plasticity. However, the double cantilever beam (both ends clamped) is unsuitable since its deformation is dominated by axial stretching. The strain gradient effects significantly increase the bending stiffness of a microbend test specimen. The deflection of a 10-μm thick beam is only a few percent of that estimated by classical plasticity.     

Design sensitivity analysis for parameters affecting geometry,elastic–viscoplastic material constant and boundary condition by consistent tangent operator-based boundary element method

This paper presents a design sensitivity analysis method by the consistent tangent operator concept-based boundary element implicit algorithm. The design variables for sensitivity analysis include geometry parameters, elastic–viscoplastic material parameters and boundary condition parameters. Based on small strain theory, Perzyna’s elastic–viscoplastic material constitutive relation with a mixed hardening model and two flow functions is considered in the sensitivity analysis. The related elastic–viscoplastic radial return algorithm and the formula of elastic–viscoplastic consistent tangent operator are derived and discussed. Based on the direct differentiation approach, the incremental boundary integral equations and related algorithms for both geometric and elastic–viscoplastic sensitivity analysis are developed. A 2D boundary element program for geometry sensitivity, elastic–viscoplastic material constant sensitivity and boundary condition sensitivity has been developed. Comparison and discussion with the results of this paper, analytical solution and finite element code ANSYS for four plane strain numerical examples are presented finally.     

Interfacial energy effects within the framework of strain gradient plasticity

In the framework of strain gradient plasticity, a solid body with boundary surface playing the role of a dissipative boundary layer endowed with surface tension and surface energy, is addressed. Using the so-called residual-based gradient plasticity theory, the state equations and the higher order boundary conditions are derived quite naturally for both the bulk material and the boundary layer. A phenomenological constitutive model is envisioned, in which the bulk material and the boundary layer obey (rate independent associative) coupled plasticity evolution laws, with kinematic hardening laws of differential nature for the bulk material, but of nondifferential nature for the layer. A combined global maximum dissipation principle is shown to hold. The higher order boundary conditions are discussed in details and categorized in relation to some peculiar features of the boundary surface, and their basic role in the coupling of the bulk/layer plasticity evolution laws is pointed out. The case of an internal interface is also studied. An illustrative example relating to a shear model exhibiting energetic size effects is presented. The theory provides a unified view on gradient plasticity with interfacial energy effects.     

现代统一塑性理论

经典的本构理论不能精确描述材料表现出的率相关性和非线性本构关系,导致现代统一塑性理论的出现.各派现代统一塑性理论形式不相同,但最根本的观点接近:①材料变形由弹性变形与非弹性变形两部分组成,不再将非弹性变形区分为塑性变形与蠕变;②存在着两种物理参量:背应力X与粘滞系数K_0 σ—X构成非弹性形变的动力,K是非弹性应变的阻力,因而ε~(in)=f((σ-X)/K);③原始材料X= 0,K=K_0,K_0为较大的标量,所以加载初始表现出线弹性,续继加载X与K按演化规律变化导致非线性形变规律的出现.各派统一塑性理论的主要差异在于X与K的演化公式不同.各有其长亦有其短.在大量试验基础上取各派之长,舍各派之短,融合为一科学的理论是今后发展的方向.      

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.     

Bammann-Chiesa-Johnson粘塑性本构模型材料参数的一种识别方法

粘塑性本构模型能否成功模拟金属高应变率大应变变形过程依赖于材料参数识别结果的好坏。由于BCJ模型考虑了应变率、温度与材料硬化之间的耦合效应以及应变率、温度历史效应,同时模型中包含了多个材料参数,因此很难通过试验直接识别模型的材料参数。本文针对BCJ模型中的耦合效应和历史效应,基于对模型中材料参数物理涵义的界定,给出了一种对材料参数解耦、分离并进行估计的方法,获得了模型材料参数估计公式,估计了材料参数的取值范围。在此基础上,编制了BCJ模型应力积分径向返回算法和粒子群优化算法的计算程序,应用重新设计了BCJ模型耦合效应和历史效应的反分析方法,在参数取值范围内对材料参数进行了优化识别。以OFHC Cu为例,应用提出的识别方法对BCJ模型的材料参数进行了识别,计算结果和试验结果符合较好。     

Strain gradient effects on microscopic strain field in a metal matrix composite

The purpose of this paper is to demonstrate the improved modeling accuracy of a finite-deformation strain gradient crystal plasticity formulation over its classical counterpart by conducting a joint experimental and numerical investigation of the microscopic details of the deformation of a whisker-reinforced metal-matrix composite. The lattice rotation distribution around whiskers is obtained in thin foils using a TEM technique and is then correlated with numerical predictions based on finite element analyses of a unit-cell of a single crystal matrix containing a rigid whisker. The matrix material is first characterized by a classical, scale-independent crystal plasticity theory. It is found that the classical theory predicts a lattice rotation distribution with a spatial gradient much higher than experimentally measured. A strain gradient crystal plasticity formulation is then applied to model the matrix. The strain gradient formulation accounts for both strain hardening and strain gradient hardening. The deformation thus predicted exhibits a strong dependence on the size of the whisker. For a constitutive length scale comparable to the whisker diameter, the spatial gradient of the lattice rotation is several times lower than that predicted by the classical theory, and hence correlates significantly better with the experimental results.     

塑性材料的统一屈服条件

从一个新的角度对经典屈服条件进行研究，并以最大剪应力是使塑性屈服发生的主要是否 是唯一因素为出发点，提出了包含材料特性以参数的统一屈服条件。     

A large strain plasticity model for anisotropic materials — composite material application

In this work a generalized anisotropic model in large strains based on the classical isotropic plasticity theory is presented. The anisotropic theory is based on the concept of mapped tensors from the anisotropic real space to the isotropic fictitious one. In classical orthotropy theories it is necessary to use a special constitutive law for each material. The proposed theory is a generalization of classical theories and allows the use of models and algorithms developed for isotropic materials. It is based on establishing a one-to-one relationship between the behavior of an anisotropic real material and that of an isotropic fictitious one. Therefore, the problem is solved in the isotropic fictious space and the results are transported to the real field. This theory is applied to simulate the behavior of each material in the composite. The whole behavior of the composite is modeled by incorporating the anisotropic model within a model based on a modified mixing theory.     

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

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

两种弹塑性算法的比较

分析了弹塑性一致切线模量算法与径向返回算法，发现一致切线模量算法并不如有关文献报道的那样，比径向返回算法有快得多的收敛速度，两个算例表明了这一点。这使对一致切线模量算法有比较全面的看法。     

A simple theory of strain gradient plasticity based on stress-induced anisotropy of defect diffusion

A new strain gradient plasticity theory is formulated to accommodate more than one material length parameter. The theory is an extension of the classical J₂ flow theory of metal plasticity to the micron scale. Distinctive features of the proposed theory as compared to other existing theories are the simplicities of mathematical formulation, numerical implementation and physical interpretation.     

A mathematical basis for strain-gradient plasticity theory. Part II: Tensorial plastic multiplier

A phenomenological, flow theory version of gradient plasticity for isotropic and anisotropic solids is constructed along the lines of Gudmundson [Gudmundson, P., 2004. A unified treatment of strain-gradient plasticity. J. Mech. Phys. Solids 52, 1379-1406]. Both energetic and dissipative stresses are considered in order to develop a kinematic hardening theory, which in the absence of gradient terms reduces to conventional J₂ flow theory with kinematic hardening. The dissipative stress measures, work-conjugate to plastic strain and its gradient, satisfy a yield condition with associated plastic flow. The theory includes interfacial terms: elastic energy is stored and plastic work is dissipated at internal interfaces, and a yield surface is postulated for the work-conjugate stress quantities at the interface. Uniqueness and extremum principles are constructed for the solution of boundary value problems, for both the rate-dependent and the rate-independent cases. In the absence of strain gradient and interface effects, the minimum principles reduce to the classical extremum principles for a kinematically hardening elasto-plastic solid. A rigid-hardening version of the theory is also stated and the resulting theory gives rise to an extension to the classical limit load theorems. This has particular appeal as previous trial fields for limit load analysis can be used to generate immediately size-dependent bounds on limit loads.     

蠕变－弹塑性－损伤耦合响应的有限元分析

本文处理各向异性非线性材料的蠕变。弹塑性－损伤耦合响应的数值计算。建议了一个计算应力的三级向后欧拉积分算法。导出了一个利用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ迭代的一般的直接应力返回映射算法。同时求解应力向量和蠕变、塑性、损伤的内状态变量。也导出了用于全局Ｎｅｗｔｏｎ－Ｒａｐｈ－ｓｏｎ迭代过程的一致性切线矩阵公式。给出的数值例题结果表明所提出的算法和公式在模拟耦合本构行为上的能力和可靠性。     

Wedge indentation of thin films modelled by strain gradient plasticity

A plane strain study of wedge indentation of a thin film on a substrate is performed. The film is modelled with the strain gradient plasticity theory by Gudmundson [Gudmundson, P., 2004. A unified treatment of strain gradient plasticity. Journal of the Mechanics and Physics of Solids 52, 1379–1406] and analysed using finite element simulations. Several trends that have been experimentally observed elsewhere are captured in the predictions of the mechanical behaviour of the thin film. Such trends include increased hardness at shallow depths due to gradient effects as well as increased hardness at larger depths due to the influence of the substrate. In between, a plateau is found which is observed to scale linearly with the material length scale parameter. It is shown that the degree of hardening of the material has a strong influence on the substrate effect, where a high hardening modulus gives a larger impact on this effect. Furthermore, pile-up deformation dominated by plasticity at small values of the internal length scale parameter is turned into sink-in deformation where plasticity is suppressed for larger values of the length scale parameter. Finally, it is demonstrated that the effect of substrate compliance has a significant effect on the hardness predictions if the effective stiffness of the substrate is of the same order as the stiffness of the film.