基于非关联流动规律的Gotoh屈服准则的参数确定方法

屈服准则对板料成形过程的理论解析、工艺优化和有限元模拟有着重要的影响.通过提高屈服准则的各向异性表征能力,可以确保成形过程的可靠性及实际预测的准确性.本文基于非关联流动法则,给出了Gotoh屈服准则一套全新的参数求解方法.在结合常用屈服准则并考虑流动规律的基础上,分别以5754O铝合金、DP980先进高强钢和SAPH440结构钢作为研究对象,进行了不同加载路径下各向异性变形行为的预测.根据Gotoh屈服准则推导的屈服函数、塑性势函数以及基于关联流动的理论函数计算出屈服应力和各向异性指数r值随加载角度的分布趋势,进而针对平面应力状态的屈服轨迹展开分析,验证了不同屈服准则和流动规律对各向异性屈服行为的预测精度.理论与实验数据对比结果表明:不同屈服准则针对同种板料在流动规律一致的情形下其表征各向异性的能力有显著差异;相同屈服准则基于不同流动规律其表征能力也具有明显差别.基于非关联流动的屈服准则能极大地提高精度,各向异性表征能力显著加强.相关结果能够为各向异性屈服准则在塑性成形领域的实际应用方案提供重要参考.     

各向异性屈服准则的发展及实验验证综述

鉴于材料的屈服行为对板料成形的重要性, 人们对各向异性屈服准则进行了长期研究. 本文对各向异性屈服准则的发展进行了较为全面的回顾, 对Hill系列、Hosford系列和Drucker系列3类屈服准则分别进行归纳. 重点介绍不同类型屈服准则的适用范围及缺陷, 总结目前国内外所采用的不同的实验验证方法, 最后指出各向异性屈服准则在数值模拟中应用的难点及今后的研究方向.      

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

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

塑性各向异性板料本构关系及冲压成形的数值模拟

推导了具有一般屈服函数形式的弹塑性速率型本构关系;给出了用于板料成形的Hill塑性各向异性屈服模型下本构关系的具体形式;用有限元动力显式计算程序MSC/DYTRAN模拟了金属板料的冲压成形;通过算例分析,考察了塑性各向异性对凸耳形成和大小以及对成形模拟结果准确性的影响;数值结果和实验结果表明：各向（厚向）异性本构模型比各向同性本构模型更真实地反映了板料的成形性。     

汽车薄钢板应力应变曲线及屈服轨迹的研究

采用十字形双向拉伸的实验方法对两种汽车用薄钢板BH220和SPEN进行了不同 加载路径下的双向拉伸试验，得到了不同应力状态下的应力应变关系曲线，同时，根据单位 体积塑性功相等的原则，确定了两种钢板等效塑性应变从0.2\%$\sim$2\%的实验屈服轨迹. 结果分析表明：不同加载路径下板料的应力应变关系不同，随着加载比例由单拉到等双拉状 态，板料的硬化指数逐步增大；实验屈服轨迹呈外凸性，且以等双拉为界的上下部分屈服轨 迹不对称，随着变形程度的增加，屈服轨迹向外扩大，但单拉时强化程度最小，而等双拉 时最大. 对BH220和SPEN钢板的实验屈服轨迹与几种常用理论屈服轨迹的比较发现，Hosford各向 异性屈服准则的理论轨迹与实验结果最为接近，Hill48准则与实验结果相差最大，此外一 向被视为只适用于各向同性材料的Mises准则与实验结果也较为接近，其他几个屈服准则的 理论屈服轨迹与实验点相差较大.     

弹塑性耦合和广义正交法则

用塑性增量理论表述岩石、土、混凝土及某些复合材料的本构性质时,必须考虑这些介质的如下特点:(1)应变软化性质;(2)弹性系数随塑性变形的发展而变化,这称为弹塑性耦合;(3)塑性势函数与屈服(或加载)函数不同,通常称之为非关联流动规律。对上述的岩土类介质,在刚性试验机上做试验,可以得到包括强化阶段和软化阶段的全部应力应     

各向异性本构关系在板料成形数值模拟中的应用

对几种能表达面内各向异性的屈服准则Hill、Barlat-Lian、Barlat进行了比较。以弹性变形服从各向同性广义虎克定律的情况下，给出了基于张量算法推导的弹塑性本构关系的一般表达式，并由此导出了相应屈服准则的弹塑性本构关系的显式表达。借助ABAQUS软件本构模块用户子程序接口，分别实现了这些屈服准则在ABAQUS的嵌入。以模拟方形盒的拉延过程为例，分析了不同的屈服准则在板料成形过程数值模拟中的应用。模拟结果表明，基于弹塑性本构关系一般表达所列出的相应屈服准则的显式表达式是正确的；在采用壳元来模拟板料成形时，采用Barlat准则的模拟结果和采用Barlat-Lian准则的结果差别不大。     

基于SMP破坏准则的饱和砂土弹塑性本构模型

利用土体的塑性流动理论，提出了用于描述饱和砂土在单调荷载作用下的应力一应变反应性质的弹塑性本构模型。土体总的变形由三部分组成：即弹性应变、与体积屈服机制相关的塑性应变和与剪切屈服机制相关的塑性应变，其中与剪切屈服机制相关的塑性应变的得出是基于SMP破坏准则。通过将模型预测的结果与试验结果进行对比，表明该模型能够较为准确地描述饱和砂土在单调加载条件下的反应性质。     

巷道围岩塑性状态判定分析方法

本文依照莫尔-库仑屈服条件和非关联流动准则, 从地应力场及地下巷道开挖工程特点出发, 提出对围岩塑性屈服状态进行判定的方法。对无支护巷道和水泥锚杆支护巷道围岩进行有限元分析表明, 巷道开挖以后及时支设锚杆, 对出现塑性屈服的围岩起到了等效提高其刚度和强度的作用, 从而增加了围岩稳定性。与现场实测结果对比说明, 非关联流动准则比关联流动准则更符合实际。     

室温下镁合金板循环塑性随动强化本构模型研究

本文在具有各向异性屈服强度和拉压不对称的CPB06屈服准则的基础上,建立了基于随动强化的循环塑性本构模型．通过引入滑移、孪晶以及去孪等不同变形模式下的背应力演化方程,对室温下镁合金板材异常循环硬化行为进行了模拟．选取了AZ31B-O和AZ31B两种镁合金板材,通过拉伸-压缩-拉伸(T-C-T)和压缩-拉伸(C-T)等不同加载路径下的部分实验曲线确定模型的参数,采用三次插值多项式建立了背应力参数与上一变形模式中累积的等效塑性应变（即预应变）之间的函数关系．使用本模型对剩下的实验曲线进行了预测,发现预测结果与实验结果有良好的一致性,说明了当前模型的正确性．     

A finite element formulation based on non-associated plasticity for sheet metal forming

In the present paper, a finite element formulation based on non-associated plasticity is developed. In the constitutive formulation, isotropic hardening is assumed and an evolution equation for the hardening parameter consistent with the principle of plastic work equivalence is introduced. The yield function and plastic potential function are considered as two different functions with functional form as the yield function of Hill [Hill, R., 1948. Theory of yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. A 193, 281–297] or Karafillis–Boyce associated model [Karafillis, A.P. Boyce, M., 1993. A general anisotropic yield criterion using bounds and a transformation weighting tensor. J. Mech. Phys. Solids 41, 1859–1886]. Algorithmic formulations of constitutive models that utilize associated or non-associated flow rule coupled with Hill or Karafillis–Boyce stress functions are derived by application of implicit return mapping procedure. Capabilities in predicting planar anisotropy of the Hill and Karafillis–Boyce stress functions are investigated considering material data of Al2008-T4 and Al2090-T3 sheet samples. The accuracy of the derived stress integration procedures is investigated by calculating iso-error maps.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

Anisotropic parameter identification using inhomogeneous tensile test

In this contribution, an inverse identification strategy of constitutive laws for elastoplastic behaviour is presented. The proposed inverse algorithm is composed on an appropriate finite element calculation combined with an optimisation procedure. It is applied to identify material anisotropic coefficients using a set up of easy performed laboratory tests. The used experimental data are the plane tensile test and the off axes tensile tests. The identified behaviour models are mainly based on Hill's quadratic yield criterion. Two cases of this yield criterion have been considered: the transverse isotropic and the orthotropic one under an associated and non-associated flow rule assumptions for each case. The yield surface has been assumed to expand isotropically (isotropic strain hardening law) as a function of the plastic work.In order to better describe anisotropic plastic properties of the studied materials, a recently planar anisotropic yield function is used. It is a non-quadratic yield criterion which takes account of anisotropic yield stresses as well as anisotropic strain ratios. It is subsequently shown that the agreement between inverse identification results and experimental measurements were improved.We prove also that the presented strategy is a good alternative to the simplified homogeneous tests assumption, especially for the plane tensile test.     

Modified Burzynski criterion with non-associated flow rule for anisotropic asymmetric metals in plane stress problems

The Burzynski criterion is developed for anisotropic asymmetric metals with the non-associated flow rule (NAFR) for plane stress problems. The presented pressure depending on the yield criterion can be calibrated with ten experimental data, i.e., the tensile yield stresses at 0°, 45°, and 90°, the compressive yield stresses at 0°, 15°, 30°, 45°, 75°, and 90° from the rolling direction, and the biaxial tensile yield stress. The corresponding pressure independent plastic potential function can be calibrated with six experimental data, i.e., the tensile R-values at 0°, 15°, 45°, 75°, and 90° from the rolling direction and the tensile biaxial R-value. The downhill simplex method is used to solve these ten and six high nonlinear equations for the yield and plastic potential functions, re- spectively. The results show that the presented new criterion is appropriate for anisotropic asymmetric metals.     

A non-associated flow rule for sheet metal forming

An improved model of material behavior is proposed that shows good agreement with experimental data for both yield and plastic strain ratios in uniaxial, equi-biaxial, and plane-strain tension under proportional loading for steel, aluminum and possibly other alloys. This model is based on a non-associated flow rule in which the plastic potential and yield surface functions are defined by quadratic functions of the stress tensor. The plastic potential aspect of the model is identical to that proposed by Hill for a quadratic anisotropic plastic potential defined in terms of measured r values. The new model differs in that the yield surface, although also defined by a quadratic function of the stress tensor, is defined independently of the plastic potential in terms of measured yield stresses. The model is developed and implemented in an FEM code that is based on a convected coordinate system. Since the associated flow rule, which assumes equivalency between the plastic potential and yield functions, is commonly accepted as a valid law in the theory of plastic deformation of most metals, the arguments for the associated flow rule are also discussed.     

A pressure-sensitive yield criterion under a non-associated flow rule for sheet metal forming

Spitzig and Richmond [Acta Metall. 32 (1984) 457] proposed that plastic yielding of both polycrystalline and single crystals of steel and aluminum alloys shows a significant sensitivity to hydrostatic pressure. They further showed that under the associated flow rule, this pressure sensitivity leads to a plastic dilatancy, i.e. permanent volume change, that is at least an order of magnitude larger than observed. Indeed, the plastic dilatancy for most materials is on the order of the measurement error and must be zero in the absence of phase change and significant void nucleation during plastic deformation. A non-associated flow rule based on a pressure sensitive yield criterion with isotropic hardening is proposed in this paper that is consistent with the Spitzig and Richmond data and analysis. The significance of this work is that the model distorts the shape of the yield function in tension and compression, fully accounting for the strength differential effect (SDE). This capability is important because the SDE is sometimes described through kinematic hardening models using only pressure insensitive yield criteria.     

Evaluation of advanced anisotropic models with mixed hardening for general associated and non-associated flow metal plasticity

The main objective of this paper is to develop a generalized finite element formulation of stress integration method for non-quadratic yield functions and potentials with mixed nonlinear hardening under non-associated flow rule. Different approaches to analyze the anisotropic behavior of sheet materials were compared in this paper. The first model was based on a non-associated formulation with both quadratic yield and potential functions in the form of Hill’s (1948). The anisotropy coefficients in the yield and potential functions were determined from the yield stresses and r-values in different orientations, respectively. The second model was an associated non-quadratic model (Yld2000-2d) proposed by Barlat et al. (2003). The anisotropy in this model was introduced by using two linear transformations on the stress tensor. The third model was a non-quadratic non-associated model in which the yield function was defined based on Yld91 proposed by Barlat et al. (1991) and the potential function was defined based on Yld89 proposed by Barlat and Lian (1989). Anisotropy coefficients of Yld91 and Yld89 functions were determined by yield stresses and r-values, respectively. The formulations for the three models were derived for the mixed isotropic-nonlinear kinematic hardening framework that is more suitable for cyclic loadings (though it can easily be derived for pure isotropic hardening). After developing a general non-associated mixed hardening numerical stress integration algorithm based on backward-Euler method, all models were implemented in the commercial finite element code ABAQUS as user-defined material subroutines. Different sheet metal forming simulations were performed with these anisotropic models: cup drawing processes and springback of channel draw processes with different drawbead penetrations. The earing profiles and the springback results obtained from simulations with the three different models were compared with experimental results, while the computational costs were compared. Also, in-plane cyclic tension–compression tests for the extraction of the mixed hardening parameters used in the springback simulations were performed for two sheet materials.