Drucker-Prager弹塑性本构关系积分:考虑非关联流动与各向同性硬化

考虑非关联流动法则以及各向同性硬化条件,采用广义中点法(Generalized Midpoint Method,GMM)进行Drucker-Prager(DP)弹塑性本构关系数值积分,给出调整后最终应力的解析解.GMM 属于隐式算法,具有良好的计算精度与数值稳定性;最近点投影法(Closest Point Project Method,CPPM)是其特例,具有一阶精度并且无条件稳定.DP 塑性势函数的特殊性质导致上述GMM 解由初始应力状态与应变增量显式确定,无需迭代求解,因此计算效率大幅提高,同时避免了迭代过程的收敛性问题.数值算例证明:当加载偏离角度较大时,GMM(ξ=1/2)的计算精度高于CPPM,可适应更大的加载步长;而对于比例加载,任意GMM 等同于精确解,采用CPPM 可获得最高的计算效率.推导了满足DP 屈服准则厚壁圆筒的弹塑性理论解,对比验证算法精度.采用非关联流动各向同性线性硬化DP 材料模拟厚壁圆筒变形局部化效应.     

RPC200塑性与各向异性损伤耦合本构建模研究

针对活性粉末混凝土 (RPC200:Reactive Powder Concrete 200MPa)的非线性行为,将连续损伤力学和塑性力学结合起来,建立了RPC200的弹塑性各向异性损伤耦合本构。在有效应力空间内,采用非关联流动法则和非均匀等向强化法则建立了RPC200的塑性本构。基于损伤能释放率建立了损伤准则,并由正交法则给出损伤演化法则,同时采用两个不同的损伤硬化法则来描述拉、压性能不同的损伤硬化。证明了该文所建模型与热动力学基本原理的相容性。最后,给出了RPC弹塑性各向异性损伤耦合本构模型数值计算流程,并以简单例证证明了所建模型的预测能力。     

Modeling axisymmetric flow through a converging channel with an arbitrary yield condition

Summary A system of ordinary differential equations describing the velocity, strain and stress fields of an isotropic rigid perfectly plastic material flowing through an axisymmetrical converging channel was derived. The plastic flow behavior was assumed to be described by an arbitrary yield condition and its associated flow rule. The solution of this problem was applied to the case of a yield function suitable for isotropic FCC polycrystals such as aluminum alloys. The singularity of the strain rate field at the channel wall where maximum friction forces occur was discussed and the influence of the yield surface shape on the velocity, strain and stress fields was investigated.     

Effects of pressure-sensitive yielding on stress distributions in crankshaft sections under fillet rolling and bending fatigue tests

In this paper, the evolution equation for the active yield surface during the unloading/reloading process based on the Drucker–Prager yield function and a recently developed anisotropic hardening rule is first presented. A user material subroutine based on the anisotropic hardening rule and the constitutive relation was written and implemented into the commercial finite element program ABAQUS. Computations were first conducted for a simple plane strain finite element model under uniaxial monotonic and cyclic loading conditions. The results indicate that the anisotropic hardening rule with the non-associated flow rule describes well the strength-differential effect and the asymmetric closed hysteresis loops as observed in the uniaxial cyclic loading tests of cast irons. Then, a two-dimensional plane strain finite element analysis of a crankshaft section under fillet rolling and subsequent bending was conducted. For the pressure sensitivity corresponding to the cast iron crankshaft of interest, the critical locations for fatigue crack initiation according to the stress distributions for pressure-sensitive materials agree with the experimental observations in bending fatigue tests of crankshaft sections.     

Experimental study on yielding and failure of an anisotropic clay

The effect of initial and stress-induced anisotropies on yielding and failure of an anisotropic clay has been studied experimentally by loading soil samples along different stress paths under triaxial stress. Five orientations of bedding planes of clay fabric were tested. In all the loading paths, it was observed that the strain increment vectors were not coaxial with the stress vectors. The degree of disassociation depended on initial fabric anisotropy and stress-induced anisotropy. The test results indicated that for the inherently anisotropic material, there is no particular requirement (a) for the yield and plastic potential surfaces to coincide, and (b) for an associated flow rule to hold.     

Experimental mechanical characterization of plastic-bonded explosives

This article deals with the characterization of the static mechanical behavior of an energetic material. Due to its constituents (crystals and a polymeric binder), the behavior is complicated to model. A specific experimental protocol has been proposed in this article. It involves uniaxial tensile and compressive tests, compression under confinement and dynamic mechanical analysis. A constitutive law has been developed. The behavior is described using a Maxwell’s model, in which all the components are influenced by an isotropic damage. The first component takes into account an elasto-plastic behavior. The yield stress evolution is described using a parabolic criterion and an isotropic hardening law. The plastic flow rule is non-associated. A linear visco-elastic behavior is used for the other components. Numerical simulations show that experimental data are quite well reproduced. The last part of the article is devoted to a discussion highlighting the future improvements.     

Viscoplastic modeling of asphalt mixes with the effects of anisotropy, damage and aggregate characteristics

This paper presents an approach for constitutive modeling of the viscoplastic behavior of asphalt mixes. This approach utilizes an anisotropic non-associated flow rule based on the Drucker–Prager yield surface. The selection of this yield surface is motivated by the field stress paths and material properties associated with permanent deformation at high temperatures. The efficacy of the model is demonstrated by analyzing data from compressive triaxial tests conducted at different confining pressures and strain rates for three different mixes. The model parameters are related to the experimental measurements of aggregate shape characteristics, aggregate surface energy, inherent anisotropic distribution of aggregates, and microstructure damage measured using X-ray computed tomography and image analysis techniques. Establishing the relationship between the model parameters and material properties is important in order to optimize the mix properties, and achieve desirable mix performance.     

考虑中间主应力和剪胀特性的深埋圆巷弹塑性应力位移解

平面应变条件下的深埋圆形巷道问题一般忽略中间主应力的影响,但这会与塑性区围岩的实际情况产生较大差异。为了充分考虑中间主应力对深埋圆形巷道的影响,基于平面应变假设与非关联流动法则将Mohr-Coulomb准则精确匹配为Drucker-Prager准则,在此基础上推导了考虑剪胀特性的理想弹塑性材料在塑性阶段的中间主应力表达式,中间主应力与剪胀角密切相关;根据所得的中间主应力表达式结合非关联流动法则,不引入任何假设,得出深埋圆巷塑性区由于剪胀角而发生体变的关系式;进一步推导了考虑中间主应力和剪胀特性的深埋圆形巷道塑性区应力位移解析式,其中径向应力、切向应力及塑性区半径的表达式与卡斯特奈(Kastner)解完全一致,但卡斯特奈(Kastner)解无法得出中间主应力,而新的位移解析式则与以往的文献完全不同;经与以往文献的位移理论解比较分析知,新的位移解答更加合理。因此考虑中间主应力和剪胀特性的解答为深埋圆形巷道的计算与设计提供一定的理论基础。     

A family of bi-potentials describing the non-associated flow rule of pressure-dependent plastic models

In this paper, we discuss in detail the derivation of the bi-potential function used to express the flow rule of pressure-dependent non-associated plastic models. In particular, we show that an infinity of equivalent expressions of the bi-potential can be derived for pressure-dependent models. The method consists in applying a transformation to the plastic strain rate or the stress or both to recover the normality rule. Then, the Fenchel inequality is used to derive the bi-potential function. The construction method itself indicates clearly that a class with infinite equivalent expressions of the bi-potential can be derived.     

Isotropic‐kinematic hardening framework to assess ratcheting response of steel samples undergoing asymmetric loading cycles

The present study intends to characterize ratcheting response of several steel alloys subject to asymmetric loading cycles through coupling the Ahmadzadeh‐Varvani kinematic hardening rule with isotropic hardening rules of Lee and Zavrel, Chaboche, and Kang. The Ahmadzadeh‐Varvani kinematic hardening rule was developed to address ratcheting progress over asymmetric stress cycles with relatively a simple framework and less number of coefficients. Inclusion of isotropic hardening rules to the framework improved ratcheting response of materials mainly over the first stage of ratcheting. Lee and Zavrel model (ISO‐I) developed an exponential function to account for accumulated plastic strain as yield surface is expanded over stage I and early stage II of ratcheting. Isotropic models by Chaboche (ISO‐II) and Kang (ISO‐III) encountered yield surface evolution in the framework by introducing an internal variable that takes into account the prior maximum plastic strain range. The choice of isotropic hardening model coupled to the kinematic hardening model is highly influenced by material softening/hardening response.     

Flow rule for the plastic deformation of particiulate metal-matrix composites

Finite element calculations are performed on models of particulate metal-matrix composites to study the applicability of a quadratic yield function and an associated flow rule. The matrix, here taken to be Al, is described by a J₂ flow rule for an isotropic material and the reinforcing particles, either SiC or TiC, are taken to be elastic. Two different types of three-dimensional models of the composite are considered: (1) a simple cubic lattice of spherical particles and (2) random digital models that are approximately isotropic. Only in the second type of model can the existence of a flow rule be established. The validity of a flow rule in random models is associated with the absence of shear planes that extend throughout the solid without intercepting any particles. Such planes can be drawn in the simple cubic lattice for some shear deformations and particle volume fractions. Localized shear bands occurring on these planes results in a shear response essentially identical to that of the unreinforced matrix material, which precludes the determination of the shear response from the uniaxial deformation of the composite.     

A 3D fourth order fabric tensor approach of anisotropy in granular media

We propose a micromechanical approach for granular media, with a particular account of the texture-induced anisotropy and of the strain localization rule. The approach is mainly based on the consideration of a fourth order fabric tensor able to capture general anisotropy which can be induced by complex distribution of contacts. Incorporation of this fourth order fabric tensor in a suitable homogenization scheme allows to determine the corresponding macroscopic elastic properties of the granular material. For this purpose, in addition to the classical Voigt upper bound, a new kinematics-based localization rule is proposed. It generalizes the one formulated by Cambou et al. [B. Cambou, Ph. Dubujet, F. Emeriault, F. Sidoroff, Eur. J. Mech. A/Solids 14 (1995) 225–276] in the case of an isotropic contact distribution. The results of the complete model compare well to numerical simulations results when available [C.S. Chang, C.L. Liao, Appl. Mech. Rev. 47 (1 Part 2) (1994) 197–207] (case of isotropic distribution of contacts). Finally, the interest of the fourth order fabric tensor based approach combined with the proposed localization rule is shown for different distributions of contacts by comparing its predictions to those given by a second order fabric tensor approach.     

On the influence of kinematic hardening on plastic anisotropy in the context of finite strain plasticity

The paper discusses the application of a newly developed material model for finite anisotropic plasticity to the simulation of earing formation in cylindrical cup drawing. The model incorporates Hill-type plastic anisotropy, nonlinear kinematic and nonlinear isotropic hardening. The constitutive framework is derived in the context of continuum thermodynamics and represents a multiplicative formulation of anisotropic elastoplasticity in the finite strain regime. Plastic anisotropy is described by means of second-order structure tensors which are used as additional tensor-valued arguments in the representation of the yield criterion and the plastic flow rule. The evolution equations are integrated by a form of the exponential map that fullfils plastic incompressibility and preserves the symmetry of the internal variables. The numerical examples investigate the influence of the hardening behaviour on an initially anisotropic yield criterion. In particular, the influence of using the kinematic hardening component of the model in addition to isotropic hardening in the earing simulations is examined. Comparisons with test data for aluminium and steel sheets display a good agreement between the finite element results and the experimental data.     

A logarithmic‐exponential backward‐Euler‐based split of the flow rule for anisotropic inelastic behaviour at small elastic strain

A basic aspect of modern algorithmic formulations for large‐deformation hyperelastic‐based isotropic inelastic material models is the exponential backward‐Euler form of the algorithmic flow rule in the context of the multiplicative decomposition of the deformation gradient. Advantages of this approach in the isotropic context include the exact algorithmic fulfilment of inelastic incompressibility. The purpose of this short work is to show that such an algorithm can be formulated for anisotropic inelastic models as well under assumption of small elastic strain, i.e. for metals. In particular, the current approach works for both phenomenological anisotropy as well as for crystal plasticity. The major difference between the current and previous approaches lies in the fact that the elastic rotation is reduced algorithmically to a dependent internal variable, resulting in a smaller internal variable system. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation of the plastic flow properties of iron single crystals

Summary The present paper deals with the numerical simulation of the plastic flow properties of iron single crystals as well as their influence on the macroscopic elastic-plastic deformation and localization behavior affected by superimposed hydrostatic pressure. Based on experimental observations the onset of plastic yielding on the microscale is described by an extended microscopic yield condition taking into account various microscopic stress components acting on the respective slip systems. In addition, to be able to compute inelastic deformations from a plastic potential, the latter is expressed in terms of workconjugate microscopic stress and strain measures which leads to a non-associated flow rule for the macroscopic plastic strain rate. On the numerical side, generalized functions for constitutive parameters will be used to be able to simulate the single crystal's microscopic deformation behavior observed in experiments. Estimates of the current microscopic stresses and strains are obtained via an efficient and remarkably stable plastic predictor-elastic corrector technique which is incorporated into a nonlinear finite element program. Numerical simulations of uniaxial tests show quantitatively the influence of hydrostatic pressure on current material data. Further numerical studies on the additional constitutive non-Schmid terms elucidate their effect on iron single crystal's macroscopic deformation and localization behavior.     

Anisotropic Gurson‐Tvergaard‐Needleman plasticity and damage model for finite element analysis of elastic‐plastic problems

Implementation and analysis of the anisotropic version of the Gurson‐Tvergaard‐Needleman (GTN) isotropic damage criterion are performed on the basis of Hill's quadratic anisotropic yield theory with the definition of an effective anisotropic coefficient to represent the elastic‐plastic behavior of ductile metals. This study aims to analyze the extension of the GTN model suitable for anisotropic porous metals and to investigate the GTN model extension. An anisotropic damage model is implemented using the user material subroutine in ABAQUS/standard finite element code. The implementation is verified and applied to simulate a uniaxial tensile test on a commercially produced aluminum sheet material for three‐dimensional and plane stress test cases. Spherical and ellipsoidal micro voids are considered in the matrix material, and their effects on the uniaxial stress‐strain response of the material are analyzed. Hill's quadratic anisotropic yield theory predicts substantially large damage evolution and a low stress‐strain curve compared with those predicted by the isotropic model. An approximate model for anisotropic materials is proposed to avoid increased damage evolution. In this approximate model, Hill's anisotropic constants are replaced with an effective anisotropy coefficient. All model‐generated stress‐strain predictions are compared with the experimental stress‐strain curve of AA6016‐T4 alloy.     

From non-planar dislocation cores to non-associated plasticity and strain bursts

In non-close-packed crystalline lattices, e.g. of bcc metals and intermetallic compounds, the stress-state dependence of the Peierls barrier for the motion of a screw dislocation violates Schmid’s law and leads to non-associated plastic flow at the continuum level. Plasticity models based upon distinct yield and flow functions are developed for both single crystals and polycrystalline aggregates that build upon atomic-level simulations of single dislocations. For a random polycrystal, isotropic forms for those functions are proposed and used to study mechanisms of macroscopic deformation. Non-associated flow is shown to have a significant effect on strain localization. Intermittent strain bursts are predicted to arise as a consequence of non-associated flow, particularly for deformations close to the plane strain state and for nearly rate-insensitive response.     

A model for the Mullins effect during multicyclic equibiaxial loading

In this paper, we derive a model to describe the cyclic stress softening of a carbon-filled rubber vulcanizate through multiple stress–strain cycles with increasing values of the maximum strain, specializing to equibiaxial loading. Since the carbon-filled rubber vulcanizate is initially isotropic, we can show that following initial equibiaxial loading the material becomes transversely isotropic with preferred direction orthogonal to the plane defined by the equibiaxial loading. This is an example of strain-induced anisotropy. Accordingly, we derive nonlinear transversely isotropic models for the elastic response, stress relaxation, residual strain and creep of residual strain in order to model accurately the inelastic features associated with cyclic stress softening. These ideas are then combined with a transversely isotropic version of the Arruda–Boyce eight-chain model to develop a constitutive relation for the cyclic stress softening of a carbon-filled rubber vulcanizate. The model developed includes the effects of hysteresis, stress relaxation, residual strain and creep of residual strain. The model is found to compare extremely well with experimental data.     

Stress concentration at irregular surfaces of anisotropic half-spaces

Summary An analytical method is presented to derive the stresses in anisotropic half-spaces with smooth and irregular surface morphologies. The half spaces can be subjected to body forces, surface tractions, and uniform far-field stresses. The general solution is expressed in terms of three analytical functions using the analytical function method of anisotropic elasticity. These three functions are then determined using a numerical conformal mapping technique and an integral equation method. Numerical examples are presented for the stress concentration at irregular surfaces induced by a uniform far-field horizontal stress. The elastic half-spaces are assumed to be transversely isotropic or isotropic, and their surface morphologies are constructed by the superposition of multiple long and symmetric ridges (mounds) and valleys (depressions). For isotropic media, the stress concentration depends only on the half-space surface geometry. It is found here that for anisotropic media, the half-space surface geometry, as well as the orientation of the planes of material anisotropy, have a great effect on the stress concentration. The degree of material anisotropy, on the other hand, has little influence on the stress concentration.     

Effects of plastic anisotropy on crack-tip behaviour

For a crack in a homogeneous material the effect of plastic anisotropy on crack-tip blunting and on the near-tip stress and strain fields is analyzed numerically. The full finite strain analyses are carried out for plane strain under small scale yielding conditions, with purely symmetric mode I loading remote from the crack-tip. In cases where the principal axes of the anisotropy are inclined to the plane of the crack it is found that the plastic zones as well as the stress and strain fields just around the blunted tip of the crack become non-symmetric. In these cases the peak strain on the blunted tip occurs off the center line of the crack, thus indicating that the crack may want to grow in a different direction. When the anisotropic axes are parallel to the crack symmetry is retained, but the plastic zones and the near-tip fields still differ from those predicted by standard isotropic plasticity.