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

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

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

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

单轴载荷下X80钢的包申格效应研究

本文通过单轴拉伸和压缩试验研究了X80管线钢的包申格效应(BE)。采用正向与反向加载方法研究材料变形历史特性。测定了X80钢的简单拉伸试验曲线,其应力-应变关系表明,该材料具有理想弹塑性特点。为了得到X80钢的BE,在不同预变形下对几个试件分别进行加载,并当给定的预应变值分别达到0.63%,0.67%,0.95%,1.27%和1.55%时就卸载。随后再进行反向加载实验,并记录应力应变曲线。该钢材反向加载时出现加工硬化,且屈服强度比正向加载时要低。正反向加载之间的屈服强度差值随着预应变增加而增大;当预应变超过0.95%时,反向屈服强度达到恒量。实验表明,X80钢的反向加载特性可用Remberg-Osgood关系拟合。最后给出了屈服强度降和预塑性应变之间的经验公式。     

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

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

软化Mohr-Coulomb材料强度参数的测定方法

为了得到试件的粘聚力和内摩擦角随轴向塑性压应变变化的曲线提出本方法。试件的弹塑性本构关系遵循相关联的Mohr-Coulomb强度准则;对常规三轴试验,试件受力进入塑性状态后,处在棱椎状屈服面的棱上,加载过程遵循Koiter流动法则。按经典塑性力学理论,推导得到轴向塑性压应变与轴向应力与轴向应变的关系;在常规三轴试验机上获得不同围压下试件的全程应力-应变曲线,进而可得到各自围压下轴向塑性压应变随加载过程的变化曲线;把来自不同围压下对应同一轴向塑性压应变的应力分别代入屈服面方程,即可求得对应的粘聚力和内摩擦角。结果表明,Mohr-Coulomb材料的两个强度参数的变化由轴向塑性压应变确定。轴向塑性压应变可以作为塑性变形的状态参数,它和试件的受力过程可以唯一确定试件的变形过程。     

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

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

含瓦斯煤的破坏和流动法则

讨论了含瓦斯煤破坏形式及相应的破坏条件,以应变分量表示应力分量建立了关联于Drucker-Prager准则和Lagrange准则的正则流动法则,考虑了剪切屈服和拉伸破坏两种状态之间的转换。剪切屈服后的流动是双向的,而拉伸破坏后流动是单向的;对于剪切屈服和拉伸破坏后的流动,应变变化路径不同,应变增量与应力增量之间的关系也不同;应力增量与应变增量的关系,不能由几个方程简捷地表示,需要用流程才能表达。     

完整岩石拉压应力状态下的张裂破坏准则

一点的应力状态可由3个主应力$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示, 当规定主应力以压为正时, 沿最大主应力$\sigma_{1}$方向将产生收缩变形, 若中间主应力$\sigma_{2}$和最小主应力$\sigma_{3}$都远小于$\sigma_{1}$, 则沿$\sigma_{2}$和$\sigma_{3}$方向会产生横向扩张变形, 当横向扩张变形达到一定极限时, 将会在平行于$\sigma _{1}$的方向产生张裂破坏. 如何建立这种张裂破坏的强度准则目前尚缺乏研究, 最大拉应变理论(第二强度理论)有时被用来解释张裂破坏, 但最大拉应变理论难以应用于三向受力状态. 本文分别用$\varepsilon_{1}$, $\varepsilon_{2}$表示最大张应变和次大张应变, 则最大拉应变理论认为当$\varepsilon_{1}$达到单向拉伸屈服应变时, 材料将产生破坏. 而本文将根据$\varepsilon_{1}+\varepsilon_{2}$之和达到极限值$\varepsilon_u$来建立张裂破坏准则. 可以证明$\varepsilon_{1} +\varepsilon_{2}$所表示的是$\sigma_{1}$主平面的面积增长率. 当$\sigma_{3}<\sigma_{2} \ll \sigma_{1}$时, 大部分岩石都具有脆性破坏的特点, 所以可将破坏前的岩石视为满足广义胡克定律的线弹性材料, 这样用$\varepsilon_{1}$, $\varepsilon_{2}$表示的强度准则可通过$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示. 在这个过程中还可考虑岩石在拉伸和压缩时具有不同弹性参数和强度的特点, 并可通过单向拉伸和单向压缩的破坏状态来确定$\varepsilon_u$. 不管$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$是压应力, 还是拉应力, 或者$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$中有拉有压的情形, 基于$\varepsilon_{1} +\varepsilon_{2} =\varepsilon_u$都可建立相应的强度准则. 所建立的准则可以反映中间应力$\sigma_{2}$对强度的影响规律, 通过建立的强度准则还可以证明: 静水拉力能引起屈服, 而静水压力不能产生屈服; 压缩破坏能使塑性体积增大, 其结果比Mohr-Coulomb准则更能反映实际情形. 并通过拉压应力状态下的试验数据验证了所建立的强度准则, 所得理论计算结果和已有的试验数据吻合得很好. 通过提出的强度准则和圆盘劈裂的试验结果, 可获得更为可靠的岩石单轴抗拉强度.      

砂土的应力路径本构模型

将微元应力路径线性逼近，转变成与其充分接近且易于计算应变的等平均应力微元和等应力比微元，计算任意加荷应力路径所产生的塑性应变，建立了双屈服面的砂土应力路径本构模型．模型体现了岩土塑性理论分量屈服和非关联流动法则的要求，在p，q平面内根据双线性的屈服线确定了加卸载准则．结合广义非线性强度理论采用变换应力三维化方法简单、合理地使模型实现三维化．通过试验数据的验证表明，砂土应力路径本构模型可以合理地描述各种应力路径下砂土的变形和强度特性。     

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

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

Fracture prediction in stretch forming using finite element simulation combined with ductile fracture criterion

Summary  A criterion for ductile fracture is introduced in the finite element simulation of sheet metal forming. From the calculated histories of stress and strain in each element, the fracture initiation site and the critical stroke are predicted by means of the ductile fracture criterion. The calculations are carried out for axisymmetric stretch forming of various aluminium alloy sheets and their laminates clad by mild steel sheets. The predictions so obtained are compared with experimental observations. The results show that the combination of the finite element simulation and the ductile fracture criterion enables the prediction of forming limit in a wide range of sheet metals. Accepted for publication 11 August 1996     

A New Experimental Technique for the Multi-axial Testing of Advanced High Strength Steel Sheets

This paper deals with the development of a new experimental technique for the multi-axial testing of flat sheets and its application to advanced high strength steels. In close analogy with the traditional tension-torsion test for bulk materials, the sheet material is subject to combined tension and shear loading. Using a custom-made dual actuator hydraulic testing machine, combinations of normal and tangential loading are applied to the boundaries of a flat sheet metal specimen. The specimen shape is optimized to provide uniform stress and strain fields within its gage section. Finite element simulations are carried out to verify the approximate formulas for the shear and normal stress components at the specimen center. The corresponding strain fields are determined from digital image correlation. Two test series are performed on a TRIP-assisted steel sheet. The experimental results demonstrate that this new experimental technique can be used to investigate the large deformation behavior of advanced high strength steel sheets. The evolution of the yield surface of the TRIP700 steel is determined for both radial and non-proportional loading paths.     

Forming simulation of aluminum sheets using an anisotropic yield function coupled with crystal plasticity theory

This paper describes the application of a coupled crystal plasticity based microstructural model with an anisotropic yield criterion to compute a 3D yield surface of a textured aluminum sheet (continuous cast AA5754 aluminum sheet). Both the in-plane and out-of-plane deformation characteristics of the sheet material have been generated from the measured initial texture and the uniaxial tensile curve along the rolling direction of the sheet by employing a rate-dependent crystal plasticity model. It is shown that the stress–strain curves and R-value distribution in all orientations of the sheet surface can be modeled accurately by crystal plasticity if a “finite element per grain” unit cell model is used that accounts for non-uniform deformation as well as grain interactions. In particular, the polycrystal calculation using the Bassani and Wu (1991) single crystal hardening law and experimental electron backscatter data as input has been shown to be accurate enough to substitute experimental data by crystal plasticity data for calibration of macroscopic yield functions. The macroscopic anisotropic yield criterion CPB06ex2 (Plunkett et al., 2008) has been calibrated using the results of the polycrystal calculations and the experimental data from mechanical tests. The coupled model is validated by comparing its predictions with the anisotropy in the experimental yield stress ratio and strain ratios at 15% tensile deformation. The biaxial section of the 3D yield surface calculated directly by crystal plasticity model and that predicted by the phenomenological model calibrated with experimental and crystal plasticity data are also compared. The good agreement shows the strength of the approach. Although in this paper, the Plunkett et al. (2008) yield function is used, the proposed methodology is general and can be applied to any yield function. The results presented here represent a robust demonstration of implementing microscale crystal plasticity simulation with measured texture data and hardening laws in macroscale yield criterion simulations in an accurate manner.     

An approximate lower bound damage-based yield criterion for porous ductile sheet metals

An approximate lower bound damage-based yield criterion is developed for isotropic porous ductile sheet metals. The matrix of the sheet metals is assumed to be elastic-perfectly plastic and obey the von Mises yield criterion with periodically distributed voids. Gurson’s unit-cell model is simplified to characterize the sheet metals. To accommodate biaxial loading, an approach of stress superposition is adopted for the stress analysis. Numerical results were calculated and compared to Gurson’s extended yield criterion and experimental results.     

Failure prediction in anisotropic sheet metals under forming operations with consideration of rotating principal stretch directions

An approximate macroscopic yield criterion for anisotropic porous sheet metals is adopted to develop a failure prediction methodology that can be used to investigate the failure of sheet metals under forming operations. Hill's quadratic anisotropic yield criterion is used to describe the matrix normal anisotropy and planar isotropy. The approximate macroscopic anisotropic yield criterion is a function of the anisotropy parameter R, defined as the ratio of the transverse plastic strain rate to the through-thickness plastic strain rate under in-plane uniaxial loading conditions. The Marciniak–Kuczynski approach is employed here to predict failure/plastic localization by assuming a slightly higher void volume fraction inside randomly oriented imperfection bands in a material element of interest. The effects of the anisotropy parameter R, the material/geometric inhomogeneities, and the potential surface curvature on failure/plastic localization are first investigated. Then, a non-proportional deformation history including relative rotation of principal stretch directions is identified in a critical element of a mild steel sheet under a fender forming operation given as a benchmark problem in the 1993 NUMISHEET conference. Based on the failure prediction methodology, the failure of the critical sheet element is investigated under the non-proportional deformation history. The results show that the gradual rotation of principal stretch directions lowers the failure strains of the critical element under the given non-proportional deformation history.     

Elastic–plastic behavior of steel sheets under in-plane cyclic tension–compression at large strain

Elastic–plastic behavior of two types of steel sheets for press-forming (an aluminum-killed mild steel and a dual-phase high strength steel of 590 MPa ultimate tensile strength) under in-plane cyclic tension–compression at large strain (up to 25% strain for mild steel and 13% for high strength steel) have been investigated. From the experiments, it was found that the cyclic hardening is strongly influenced by cyclic strain range and mean strain. Transient softening and workhardening stagnation due to the Bauschinger effect, as well as the decrease in Young's moduli with increasing prestrain, were also observed during stress reversals. Some important points in constitutive modeling for such large-strain cyclic elasto-plasticity are discussed by comparing the stress–strain responses calculated by typical constitutive models of mixed isotropic–kinematic hardening with the corresponding experimental observations.     

Advances in experiments on metal sheets and tubes in support of constitutive modeling and forming simulations

This work is a review of experimental methods for observing and modeling the anisotropic plastic behavior of metal sheets and tubes under a variety of loading paths, such as biaxial compression tests; biaxial tension tests on metal sheets and tubes using closed-loop electrohydraulic testing machines; the abrupt strain path change method for detecting a yield vertex and subsequent yield loci without unloading; in-plane stress reversal tests on metal sheets; and multistage tension tests. Observed material responses are compared with the predictions of phenomenological plasticity models. Special attention is paid to the plastic deformation behavior of materials commonly used in industry, and to verifying the validity of conventional anisotropic yield criteria for those materials and associated flow rules at large plastic strains. The effects of using appropriate anisotropic yield criteria on the accuracy of simulations of forming defects, such as large springback and fracture, are also presented to highlight the importance of accurate material testing and modeling.