Effects of the stress state on plasticity and ductile failure of an aluminum 5083 alloy

The experimental and numerical work presented in this paper reveals that stress state has strong effects on both the plastic response and the ductile fracture behavior of an aluminum 5083 alloy. As a result, the hydrostatic stress and the third invariant of the stress deviator (which is related to the Lode angle) need to be incorporated in the material modeling. These findings challenge the classical J₂ plasticity theory and provide a blueprint for the establishment of the stress state dependent plasticity and ductile fracture models for aluminum structural reliability assessments. Further investigations are planned to advance, calibrate and validate the new plasticity and ductile fracture models.     

一种适用于岩土的扩展强度及屈服准则

土壤材料是一种典型的摩擦型材料,然而天然岩石却具有一定的凝聚力,而金属材料则完全是凝聚型材料. 在分析三种典型的材料强度准则表达式基础上,即SMP,Lade-Duncan以及广义Von-Mises准则,通过利用应力张量的不变量表达形式,提出了一种扩展准则即VML准则,该准则能够分别退化为上述3种典型准则. 在偏平面上,新准则能够描述从曲边三角形到圆形在内的多种开口形态;在子午面上,采用幂函数作为破坏准则公式,能够描述静水压力对于强度特性影响的非线性性质. 而对于土壤的屈服性质,岩土材料具有典型的压剪耦合特性,因此,为了描述剪切与等方向压缩两种路径下的体积耦合现象,采用水滴型屈服面作为屈服准则. 对于偏平面上的截面形状,讨论了给定球应力下偏应力强度值的分布形式及特点,讨论了应力罗德角对于偏平面上强度曲线的凹凸性的影响. 最后,通过多种材料的破坏与屈服试验成果,用所提新准则进行了验证. 通过强度以及屈服特性测试对比,验证了所提VML准则的合理性.      

Damage accumulation and fracture initiation in uncracked ductile solids subject to triaxial loading

A damage plasticity model for ductile fracture is proposed. This model is established on the cylindrical coordinate system of principal stress space. Experimental results show that fracture initiation in uncracked ductile solids is sensitive to the hydrostatic pressure and dependent on the Lode angle. The joint effects of pressure and Lode angle define a fracture envelope in principal stress space. Plastic deformation induced damage is calculated by an integral of the damage rate measured at current loading and deformation status with respect to the fracture envelope. A power law damage rule is proposed to characterize the nonlinearity in damage accumulation. A damage-related weakening factor is adopted to describe the material deterioration. The material parameters are calibrated from standard laboratory tests. The proposed model is numerically implemented. Four simulations with emphasis on crack path prediction are presented.     

Multiaxial constitutive model accounting for the strength-differential in Inconel 718

The nickel-base alloy Inconel 718 exhibits a strength-differential, that is, a different plastic flow behavior in uniaxial tension and uniaxial compression. A phenomenological viscoplastic model founded on thermodynamics has been extended for material behavior that deviates from classical metal plasticity by including all three stress invariants in the threshold function. The model can predict plastic flow in isotropic materials with or without a flow stress asymmetry as well as with or without pressure dependence. Viscoplastic material parameters have been fit to pure shear, uniaxial tension, and uniaxial compression experimental results at 650°°C. Threshold function material parameters have been fit to the strength-differential. Four classes of threshold functions have been considered and nonproportional loading of hollow tubes, such as shear strain followed by axial strain, has been used to select the most applicable class of threshold function for the multiaxial model as applied to Inconel 718 at 650 °C. These nonproportional load paths containing corners provide a rigorous test of a plasticity model, whether it is time-dependent or not. A J₂J₃ class model, where J₂ and J₃ are the second and third effective deviatoric stress invariants, was found to agree the best with the experimental results.     

On stress-state dependent plasticity modeling: Significance of the hydrostatic stress, the third invariant of stress deviator and the non-associated flow rule

It has been shown that the plastic response of many materials, including some metallic alloys, depends on the stress state. In this paper, we describe a plasticity model for isotropic materials, which is a function of the hydrostatic stress as well as the second and third invariants of the stress deviator, and present its finite element implementation, including integration of the constitutive equations using the backward Euler method and formulation of the consistent tangent moduli. Special attention is paid for the adoption of the non-associated flow rule. As an application, this model is calibrated and verified for a 5083 aluminum alloy. Furthermore, the Gurson-Tvergaard-Needleman porous plasticity model, which is widely used to simulate the void growth process of ductile fracture, is extended to include the effects of hydrostatic stress and the third invariant of stress deviator on the matrix material.     

New true-triaxial rock strength criteria considering intrinsic material characteristics

A reasonable strength criterion should reflect the hydrostatic pressure effect, minimum principal stress effect,and intermediate principal stress effect. The former two effects can be described by the meridian curves, and the last one mainly depends on the Lode angle dependence function. Among three conventional strength criteria, i.e.Mohr–Coulomb(MC), Hoek–Brown(HB), and Exponent(EP) criteria, the difference between generalized compression and extension strength of EP criterion experience a firstly increase then decrease process, and tends to be zero when hydrostatic pressure is big enough. This is in accordance with intrinsic rock strength characterization. Moreover, the critical hydrostatic pressure I_c corresponding to the maximum difference of between generalized compression and extension strength can be easily adjusted by minimum principal stress influence parameter K. So, the exponent function is a more reasonable meridian curves, which well reflects the hydrostatic pressure effect and is employed to describe the generalized compression and extension strength.Meanwhile, three Lode angle dependence functions of L_(MN),L_(WW), and L_(YMH), which unconditionally satisfy the convexity and differential requirements, are employed to represent the intermediate principal stress effect. Realizing the actual strength surface should be located between the generalized compression and extension surface, new true-triaxial criteria are proposed by combining the two states of EP criterion by Lode angle dependence function with a same lode angle. The proposed new true-triaxial criteria have the same strength parameters as EP criterion. Finally, 14 groups of triaxial test data are employed to validate the proposed criteria. The results show that the three new true-triaxial exponent criteria,especially the Exponent Willam-Warnke criterion(EPWW)criterion, give much lower misfits, which illustrates that the EP criterion and L_(WW) have more reasonable meridian and deviatoric function form, respectively. The proposed new true-triaxial strength criteria can provide theoretical foundation for stability analysis and optimization of support design of rock engineering.     

Tension–torsion fracture experiments—Part I: Experiments and a procedure to evaluate the equivalent plastic strain

Ductile failure experiments on a double notched tube (DNT) specimen subjected to a combination of tensile load and torque that was applied at a fixed ratio is presented. The experimental results extend those in Barsoum and Faleskog (2007a) down to zero stress triaxiality. A new and robust evaluation procedure for such tests is proposed, and a simple relation for the equivalent plastic strain at failure for combined normal and shear deformation, respectively, is developed. Tests were carried out on the medium strength medium hardening steel Weldox 420, and the high strength low hardening steel Weldox 960. The experimental results unanimously show that ductile failure not only depends on stress triaxiality, but is also strongly affected by the type of deviatoric stress state that prevails, which can be quantified by a stress invariant that discriminates between axisymmetric stressing and shear dominated stressing, e.g., the Lode parameter. Additional experiments on round notch bar (RNB) specimens are recapitulated in order to give a comprehensive account on how ductile failure depends on stress triaxiality, ranging from zero to more than 1.6, and the type of stress state for the two materials tested. This provides an extensive experimental data base that will be used to explore an extension of the Gurson model that incorporates damage development in shear presented in Xue et al. (2013) (Part II).     

Modification of the Gurson Model for shear failure

Recent experimental evidence points to limitations in characterizing the critical strain in ductile fracture solely on the basis of stress triaxiality. A second measure of stress state, such as the Lode parameter, is required to discriminate between axisymmetric and shear-dominated stress states. This is brought into the sharpest relief by the fact that many structural metals have a fracture strain in shear, at zero stress triaxiality, that can be well below fracture strains under axisymmetric stressing at significantly higher triaxiality. Moreover, recent theoretical studies of void growth reveal that triaxiality alone is insufficient to characterize important growth and coalescence features. As currently formulated, the Gurson Model of metal plasticity predicts no damage change with strain under zero mean stress, except when voids are nucleated. Consequently, the model excludes shear softening due to void distortion and inter-void linking. As it stands, the model effectively excludes the possibility of shear localization and fracture under conditions of low triaxiality if void nucleation is not invoked. In this paper, an extension of the Gurson model is proposed that incorporates damage growth under low triaxiality straining for shear-dominated states. The extension retains the isotropy of the original Gurson Model by making use of the third invariant of stress to distinguish shear dominated states. The importance of the extension is illustrated by a study of shear localization over the complete range of applied stress states, clarifying recently reported experimental trends. The extension opens the possibility for computational fracture approaches based on the Gurson Model to be extended to shear-dominated failures such as projectile penetration and shear-off phenomena under impulsive loadings.     

On some features of the effective behaviour of porous solids with J2- and J3-dependent yielding matrix behaviourSur quelques effets de l'angle de Lode sur les caractéristiques macroscopiques de matériaux en rupture ductile

Some features od the constitutive behaviour of voided materials taking into account possible effects of the Lode angle in the yielding behaviour of the matrix are discussed. The Gurson approach is used to this end. After providing a parametric representation of the effective behaviour of such materials, some closed-form results are given for pure shear stress states and also at very high stress triaxialities. In the former case corresponding to a zero macroscopic mean stress, the contour of the yield domain in the π-plane has exactly the shape of the yield surface of the matrix in the deviatoric plane, but a size reduced by a factor $1?f$, with f the porosity of the voided material. In the latter, effective yield stresses for the voided material are slightly different from the Gurson result and found to be set by the yield stress at a microscopic stress Lode angle $\frac{\pi }{3}$ for very high positive triaxiality and by the yield stress at a microscopic stress Lode angle 0 for very high negative triaxiality. This last result is extended for porous materials with yielding depending further on the hydrostatic stress, fully exhibiting the interaction between volumetric and shear interactions on the yielding behaviour of isotropic porous materials. Applications to many usual yielding criteria for the matrix are also provided.     

Constitutive modeling of pressure dependent plasticity and fracture in solder joints

A new constitutive model describing the pressure dependence of plasticity and fracture in solder joints is proposed. The pressure dependence of the solders flow stress is obtained from the Peierls stress necessary to move mobile dislocations. Conservation of volume during plastic deformation is achieved through a deviatoric associated flow rule. Fracture is considered as a result of damage accumulation depending on the amount of plastic strain and triaxiality of stress the sample has experienced. The constitutive model was implemented in the FEM code Abaqus using the user subroutine UMAT. The model predictions were compared to experimental results of tensile tests performed with Sn-3.5Ag solder joints of various thicknesses. Ultimate tensile strength and ductility showed a strong dependency on the solder joints geometry. It is demonstrated that the model can nicely be fitted to experimental results.     

弹体贯穿混凝土数值模拟的改进材料模型

研究混凝土结构在冲击载荷下的力学特性对武器以及防护结构的设计和评估具有重要意义,而合适的材料模型可以更准确地预测混凝土结构的力学行为和破坏模式。因此,本文中提出了一种改进的混凝土塑性损伤材料模型来描述其在冲击载荷下的力学响应。该改进模型考虑了压力-体积应变关系、应变率效应、洛德角效应和塑性损伤累积对混凝土材料力学特性的影响,并引入了一个与损伤相关的硬化/软化函数来描述压缩状态下的应变硬化和软化行为。随后,通过对3个独立的强度面进行线性插值得到了该改进模型的破坏强度面,并采用部分关联流动法则考虑了混凝土材料的体积膨胀特性。最后,开展了单个单元在不同加载条件下和弹体贯穿钢筋混凝土靶的数值模拟,验证了该改进模型的可行性、准确性以及预测性能提升。     

Numerical simulation of fracture mode transition in ductile plates

Fracture mode of ductile solids can vary depending on the history of stress state the material experienced. For example, ductile plates under remote in-plane loading are often found to rupture in mode I or mixed mode I/III. The distinct crack patterns are observed in many different metals and alloys, but until now the underlying physical principles, though highly debated, remain unresolved. Here we show that the existing theories are not capable of capturing the mixed mode I/III due to a missing ingredient in the constitutive equations. We introduce an azimuthal dependent fracture envelope and illustrate that two competing fracture mechanisms, governed by the pressure and the Lode angle of the stress tensor, respectively, exist ahead of the crack tip. Using the continuum damage plasticity model, we demonstrate that the distinctive features of the two crack propagation modes in ductile plates can be reproduced using three dimensional finite element simulations. The magnitude of the tunneling effect and the apparent crack growth resistance are calculated and agree with experimental observations. The finite element mesh size dependences of the fracture mode and the apparent crack growth resistance are also investigated.     

基于应力三轴度和罗德参数的6061和7075铝合金材料断裂失效分析

分别对6061铝合金和7075铝合金材料的缺口圆棒试件和凹槽平板试件进行准静态拉伸试验,并采用ABAQUS软件对拉伸过程进行数值模拟。模拟结果与试验测试结果吻合很好,验证了有限元模型的合理性和可靠性。通过有限元模拟,分别给出了不同试件的应力三轴度和罗德参数随等效塑性应变的变化曲线以及两种材料的失效轨迹,并对它们进行了分析讨论。结果表明:形状相同、材料性质不同的试件,应力三轴度的演化规律不同;材料的失效应变受应力三轴度和罗德参数的影响,并且不同性质的材料对罗德参数的敏感性不同。     

Effects of superimposed hydrostatic pressure on sheet metal formability

The effect of superimposed hydrostatic pressure on sheet metal formability is studied analytically and numerically. A tensile sample of power-law hardening material under superimposed hydrostatic pressure is first analyzed using the classical isotropic plasticity theory. It is demonstrated that the superimposed hydrostatic pressure p lowers the true tensile stress level at yielding by the amount of p, while material work-hardening is independent of p. It is showed, based on the Considère criterion, that the superimposed hydrostatic pressure increases the uniform strain. The effect of superimposed hydrostatic pressure on sheet metal formability is further assessed by constructing the Forming Limit Diagram (FLD) based on the M-K approach. It is found that the superimposed pressure delays the initiation of necking for any strain path. The difference in predicted FLDs between the superimposed hydrostatic pressure and the stress component normal to the sheet plane is demonstrated. Finally, the effect of superimposed hydrostatic pressure on fracture in round bars under tension is studied numerically using the finite element method, based on the Gurson damage model. The experimentally observed transition of the fracture surface, from the cup-cone mode under atmospheric pressure to a slant structure under high pressure, is numerically reproduced.     

钨合金弹体对混凝土靶的超高速侵彻机理

为研究钨合金弹体超高速侵彻混凝土靶的相关机理,构建了适用于超高速撞击的金属强度模型、失效模型和混凝土的本构模型,对93钨合金弹体超高速撞击混凝土靶问题进行了数值模拟。开展了钨合金弹体超高速撞击混凝土靶实验,分析了靶板成坑特性,研究了侵彻总深度和残余弹体长度随撞击速度的变化规律,理论分析了长杆钨弹超高速撞击混凝土的侵彻模型和混凝土靶内的应力波传播。得到以下主要结论:（1）利用金属及混凝土的新本构模型获得的超高速撞击混凝土靶的破坏形貌数值模拟结果与实验结果一致;（2）超高速撞击条件下混凝土靶成坑为“弹坑＋弹洞”形,成坑体积与弹体动能近似成正比;（3）超高速撞击条件下,侵彻深度随弹速提高呈现先增大后减小的现象,高速段侵深降低是弹体经历销蚀侵彻后“刚体侵彻阶段”减少造成的;（4）建立的钨合金超高速撞击混凝土侵彻分析模型,可用来预估侵彻深度、残余弹长、蘑菇头直径等参数;（5）采用建立的超高速撞击混凝土靶内应力波传播理论模型得到的计算结果与实验结果吻合较好。     

A homogenization-based constitutive model for isotropic viscoplastic porous media

An approximate model based on the “second-order” nonlinear homogenization method is proposed to estimate the effective behavior of isotropic, viscoplastic, porous materials. The model is constructed in such a way that it reproduces exactly the behavior of a “composite-sphere assemblage” in the limit of hydrostatic loadings, and therefore coincides with the hydrostatic limit of Gurson’s criterion in the special case of ideal plasticity. As a consequence, the new model improves on earlier homogenization estimates, which have been found to be quite accurate for low triaxialities but overly stiff for sufficiently high triaxialities and nonlinearities. Additionally, the estimates delivered by the model exhibit a dependence on the third invariant of the macroscopic stress tensor, which has a nontrivial effect on the effective response of the material at moderate triaxialities. The proposed model is compared with exact results obtained for a special class of porous materials with sequentially laminated microstructures. The agreement is found to be quite good for the entire range of stress triaxialities, and all values of the porosity and nonlinearity considered.     

复杂应力状态下岩土材料非线性本构模型研究

基于统一强度理论,将洛德参数引入统一强度理论,推导出材料的统一强度参数,进而分析中间主应力以及主剪应力系数对材料统一强度参数的影响。在此基础上,对邓肯张双曲线模型进行改进,使得该模型能够反映复杂应力状态下的应力应变关系。采用粘土的平面应变试验进行验证,结果表明:在平面应变情况下,当b=0.4时,模型能够较好地反映粘土的应力应变关系,及其强度参数。分析了洛德参数以及中间主应力系数对模型的影响,进一步说明岩土材料存在主应力效应。该模型能够反映不同材料在复杂应力状态下的应力应变关系,有其更为广泛的适用性。