参数求解方法对屈服准则的各向异性行为描述能力的影响

详细分析基于应力各向异性和变形各向异性两种求解Hill48屈服准则参数的方法。在给出两种各向异性参数求解表达式的基础上,具体分析Hill48屈服准则本身的局限性。以5754O铝合金板为研究对象,进行不同方向的单向拉伸试验。采用两种各向异性参数求解方法,基于Hill48屈服准则推导不同方向拉伸过程中的理论应力-应变曲线和拉伸过程中的变形规律。通过对比理论与试验结果具体分析参数求解方法对屈服准则精度的影响。基于两种参数求解方法,进行5754O铝合金板拉深试验的有限元模拟,讨论不同求解方法对凸耳现象的描述精度。得出结论：当对应力各向异性为主的问题进行分析时,应采用应力各向异性法求解;当对变形各向异性为主的问题进行分析时,则应采用变形各向异性法求解。研究结果对屈服准则在板料成形方面的合理应用具有重要的参考价值。     

比例/非比例加载路径下5754O汽车用铝合金板各向异性变形行为的精确解析

随着越来越高的汽车轻量化需求,铝合金板在现代汽车工业中的应用越来越广。在不同加载路径下,包括比例和非比例加载,5754O铝合金板在塑性成形过程中具有复杂的各向异性规律。试验表明5754O铝合金板的各向异性规律随变形量的增加会发生改变,因此在常参数屈服准则理论框架下,基于传统的单一曲线假设难以对5754O铝合金板在整个塑性变形过程中的各向异性行为进行精确描述。鉴于上述问题,并同时考虑到大变形过程中材料变形的稳定性,对Yld2000-2d屈服准则进行改进。基于改进的Yld2000-2d屈服准则和单一曲线假设推导不同方向的单向拉伸应力应变曲线,并与试验结果进行了对比。结果表明,与原始的Yld2000-2d屈服准则不同,基于改进的Yld2000-2d屈服准则,传统的单一曲线假设仍然适用于5754O铝合金板各向异性问题。给出不同强化方式在比例加载路径下的统一性和非比例加载路径下的分散性证明。基于改进的Yld2000-2d屈服准则和等向强化和混合两种强化方式,推导非比例加载路径下板料的应力应变曲线。基于试验结果,验证了推导的理论曲线的精度。实现了5754O铝合金板在比例和非比例加载路径下变形行为的精确描述,为其工业应用提供了重要的理论支撑。     

Finite element method for sheet forming based on an anisotropic strain-rate potential and the convected coordinate system

A variational formulation and the associated finite element (FE) equations have been derived for general three-dimensional deformation of a planar anisotropic rigid-plastic sheet metal which obeys the strain-rate potential proposed by [11.]. By using the natural convected coordinate system, the effect of geometric change and the rotation of planar anisotropic axes were efficiently considered. In order to check the validity of the present formulation, a cylindrical cup deep drawing test was modeled for a 2008-T4 aluminum alloy sheet sample. Eating simulations were performed and planar anisotropic material properties were experimentally determined. Even though quantitative agreement was not fully achieved, reasonably good agreement was found between the FE simulation and the experiment in thickness strain distribution and caring. No numerical difficulty due to planar anisotropy was encountered, and the computational procedure was found to be very stable, requiring only moderate computational time. The results have shown that the present formulation for planar anisotropic deformation can provide a good basis for the analysis of sheet metal forming processes for planar anisotropic materials, especially for aluminum alloy sheets.     

提高铝合金板料成形极限的新型可控压边力技术

提出了一种可随位置和行程同步优化压边力的新型控制策略，并将其用于不等深铝合金阶梯盒零件的冲压，在零件成形过程中按照优化策略在不同压边圈位置随冲压行程施加最优压边力，使其冲压深度从恒定压边力下的45mm提高到60mm。在多点变压边力压机上的实验验证了这种优化策略的有效性，从而提高了可控压边力技术的鲁棒性和实用性。     

Earing predictions for strongly textured aluminum sheets

Metallic alloy sheets develop crystallographic texture and plastic anisotropy during rolling. Deep drawing of a cylindrical cup from a rolled sheet is one of the typical forming operations where the effect of this anisotropy is most evident. Generally, in the finite element analyses of this process, the evolution of anisotropy during forming is neglected. In this paper, results of an experimental program carried out to quantify the anisotropy of aluminum alloy AA5042-H2 are reported. In addition to tensile tests along seven directions in the plane of the sheet, cup-drawing tests were conducted. It was observed that the material displays eight ears. The effects of the evolution in anisotropy and the directionality in hardening on the predictions of the earing profile for this material are investigated using a new methodology that incorporates multiple hardening curves corresponding to uniaxial tension along several orientations with respect to the rolling direction, and to biaxial tension. Yielding is described using the anisotropic yield function Yld2000-2D [1] and a form of CPB06ex2 yield function [2], which is tailored for metals with no tension–compression asymmetry. It is shown that even if distortional hardening is neglected, the latter yield function predicts a cup with eight ears as was observed experimentally. Consideration of distortional only leads to improved accuracy in prediction of the non-uniformity of the cup height profile.     

不同硬化模型对铝合金板冲压成形模拟结果的影响

研究不同塑性变形硬化模型对汽车5182-O铝合金板材冲压成形模拟结果的影响。采用材料单向拉伸试验得到应力应变关系曲线,基于Hollomom、Krupskowsky与Power方程对曲线进行拟合,建立材料室温下塑性变形硬化模型,对厚度为1.5 mm和0.85 mm的5182板材进行冲压试验和有限元模拟分析,对比分析冲压试验与模拟结果。试验与模拟结果显示,当板料厚度为1.5 mm时,板料冲压试验的成形力最大为42.95 kN,板料拉深深度为30.58 mm,基于Power方程计算得到的最大成形力为41.5kN与试验结果比较接近,Hollomom方程计算得到的拉深深度为30.546 mm,板材成形厚度分布与试验结果比较接近;当板料厚度为0.85 mm时,板料冲压试验的成形力最大为34.47kN,板料拉深深度为33.792 mm,基于Power方程计算得到的最大成形力为34.27 kN与试验结果比较接近,Hollomom方程计算得到的拉深深度为33.636 mm,板材成形厚度分布与试验结果比较接近。基于三种硬化模型铝合金冲压成形过程的计算模拟分析结果,并通过与试验对比得到不同硬化模型对铝合金板材冲压成形计算模拟的影响,进一步为汽车铝合金覆盖件在成形工艺的研究分析提供理论指导。     

圆筒形件拉深失稳及各因素影响分析

对板料成形中圆筒形件拉深的破裂失稳及产生破裂失稳的临界压边力进行研究.由于凸、凹模圆角及其间隙的存在,圆筒形件拉深的筒壁区实际为凸、凹模圆角之间的公切线部分.根据Mises-Hill屈服函数及Tresca准则求出凸缘变形区、凹模圆角区和筒壁区的应力分布,得到危险断面处的应力表达式,从而求出不产生破裂失稳的临界压边力的解析表达式,并进一步分析获得拉深比、硬化指数、厚向异性系数、摩擦因数以及径向推力等因素对临界压边力的影响规律.采用液压压边与周缘加径向推力的拉深模具对08Al板料进行拉深试验,试验结果与理论计算结果具有很好的一致性.     

圆筒件拉深变形过程中应力分布的数值模拟

拉深是拉深类零件的一个基本成型工序,为考察冲压拉深变形过程中材料的应力分布状态,利用ANSYS有限元软件对紫铜板拉深过程进行数值模拟。结果表明,拉深变形过程中变形主要集中在处于凹模端面上的凸缘部分,凹模圆角部变形比较复杂,易发生应力集中现象而导致材料被拉断,合理的压边可有效防止拉深起皱。     

An elastoplastic analysis of flange wrinkling in deep drawing process

The onset of flange wrinkling of a deep drawing cup is analyzed as an elastoplastic bifurcation problem. The flange is modeled as an elastoplastic annular plate subject to axisymmetric radial tension along its inner edge. As observed in the laboratory as well as practical industrial applications, aluminum alloy sheets usually wrinkle in the plastic range. Therefore, the critical condition governing the onset of elastoplastic wrinkling is formulated within the context of the general bifurcation theory. A closed-form solution for the critical drawing stress is developed based on an assumed nonlinear plastic stress field and the deformation theory of plasticity. The theory properly accounts for the plastic anisotropy of the aluminum sheets and the critical drawing stress at the onset of wrinkling is also compared against the one employing the flow theory of plasticity. The predicted critical bifurcation stress and the wave numbers are compared to those obtained by Senior's one-dimensional theory. It is demonstrated that there is a strong dependency of the critical bifurcated stress at the onset of wrinkling on the shear stress induced on the flange. The effects of flange width, drawing ratios, material properties, strain hardening on the onset of wrinkling are investigated. The differences between the present theoretical approach and Senior's theory are emphasized.     

Influence of constitutive model in springback prediction using the split-ring test

The aim of this paper is to compare several plastic yield criteria to show their relevance on the prediction of springback behavior for a AA5754-0 aluminum alloy. An experimental test similar to the Demeri Benchmark Test [Demeri MY, Lou M, Saran MJ. A benchmark test for springback simulation in sheet metal forming. In: Society of Automotive Engineers, Inc., vol. 01-2657, 2000] has been developed. This test consists in cutting a ring specimen from a full drawn cup, the ring being then split longitudinally along a radial plan. The difference between the ring diameters, before and after splitting, gives a direct measure of the springback phenomenon, and indirectly, of the amount of residual stresses in the cup. The whole deep drawing process of a semi-blank and numerical splitting of the ring are performed using the finite element code Abaqus. Several material models are analyzed, all considering isotropic and kinematic hardening combined with one of the following plasticity criteria: von Mises, Hill’48 and Barlat’91. This last yield criterion has been implemented in Abaqus. Main observed data are force-displacement curves during forming, cup thickness according to material orientations and ring gap after splitting. The stress distributions in the cup, at the end of the drawing stage, and in the ring, after springback, are analyzed and some explanations concerning their influence on springback mechanisms are given.     

利用激光局部硬化效应提高2B06铝合金板材拉深成形性研究

针对2B06铝合金复杂零件成形困难问题,提出了利用激光热处理局部硬化提高板材成形性的思路。在通过激光热处理试验研究了铝合金板的激光硬化效应的基础上,采用数值模拟计算了铝合金板激光热处理过程中激光光斑路径和其周边热影响区域的峰值温度场,并通过实际测温验证了其准确性。提出并构建了耦合性能梯度的差性坯料模型,对激光局部硬化的杯形件拉深成形性进行了模拟和试验研究。结果表明,激光扫描方式可以形成稳定的梯度温度场并对周边非加热区影响较小,且可以通过多道次扫描方式设计不同宽度范围的大梯度差性板材坯料。力学性能试验表明激光热处理可以有效地提高铝合金的局部加工硬化能力,这种效应可以有效抑制杯形件拉深时圆角大变形区的减薄,从而提升了板材的拉深性能。在上述基础上,将激光局部热处理用于2B06铝合金航空复杂口框零件的成形,通过设计激光处理路径和参数,获得合理的局部硬化区域,可有效地避免在加强筋处出现过度减薄导致的破裂,大大提高复杂零件的成形性。     

1060纯铝微杯形件微拉深成形研究

微拉深成形以其高效率、低成本、高质量等优点成为一种可行的制造开口、空心、薄壁微型零件的微细加工方法。然而,由于微型化的影响,使得微拉深要比宏观拉深更加困难。研制精密微型落料-拉深复合模,成形出内径仅为1 mm的微杯形件。试验材料采用厚度为50 μm的1060纯铝,在真空氛围下经300 ℃分别退火1 h和8 h。微拉深试验在SANS CMT-5504电子万能试验机上进行,试验速度为0.05 mm/s,润滑剂为聚乙烯薄膜。研究退火工艺和压边力对拉深力、拉深比、制耳和起皱的影响。试验结果表明,顺利成形出质量良好、拉深比为2.1的微杯形件;经退火处理的坯料,拉深力降低,拉深比增大,制耳减弱;随着压边力的增加,起皱减弱,但当压边力超过一定数值后则导致拉裂。     

A rigid-plastic finite-element formulation for the analysis of general deformation of planar anisotropic sheet metals and its applications

General three-dimensional deformation of a planar anisotropic rigid-plastic sheet metal which obeys Hill's quadratic yield function is considered in the present analysis. The explicit expressions of the yield criterion and the related equations are given in the surface curvilinear co-ordinate system. A convective co-ordinate system is used in order to take into account the effect of geometric change during one step in the incremental analysis. The expression of the effective strain increment during one step is obtained in closed form by introducing assumptions on the deformation path and by proper consideration of the rotation of the axes of anisotropy during deformation. Considering the equilibrium at the deformed state, a variational formulation is derived to determine the deformation during a step. The corresponding finite-element equations are found in order to analyze the general deformation of planar anisotropic sheet metals. Two computational examples are chosen and computed by using the developed FEM program in order to verify the present formulation. The flange deformation in deep drawing of a circular diaphragm of planar anisotropy is analyzed and compared with the existing solution. The effects of planar anisotropy in the circular diaphragm bulge test are also investigated.     

数值模拟技术及其在冲压加工中的应用

应用Autoform数值模拟软件的板料成形仿真技术,以某筒形件为例,探讨了从CAD造形到CAE拉深成形数值模拟的全过程,预测了板料在成形过程中可能出现的缺陷,为模具制造过程中及时修改和优化冲压工艺方案以及有效解决模具形状参数化调整问题节省了大量时间,同时为模具设计提供了可靠性保证。     

Experimental and finite element investigation on wrinkling of circular single layer and two-layer sheet metals in deep drawing process

This paper deals with the comparison of wrinkling behavior of both single and two-layer sheets in the deep drawing process. Unfortunately, due to the significant difference in the material properties of the two layers, it is very challenging for the finite element method to predict the location and relative possibility of wrinkling in the multilayer sheets. Blank holder force (BHF) has a significant effect on the failure mode of sheet metals. By decreasing BHF, wrinkling turns out as the principal failure mode, whereas its increase results in splitting and fracture. Thus, this paper investigates the effect of blank holder forces on wrinkling and fractures in drawing process of single and two layers. Moreover, the effect of material and arrangement of layers on wrinkling, fracture, and required deformation forces in the cylindrical deep drawing will be discussed. Results show that the optimum blank holder force for two-layer sheets is affected by the material and arrangements of lay-ups.     

润滑条件下铝合金板成形模拟中摩擦模型的研究

对润滑条件下铝合金薄板筒形件拉深成形过程进行了工艺分析，在作出一些简化和假设的基础上，建立了基于流体润滑的筒形件拉深成形过程摩擦模型；将该模型运用到板料成形过程有限元分析中，计算出筒形件成形时摩擦系数动态变化的数值；开发出基于探针测试的铝合金板温成形动态过程摩擦测试系统，并对筒形件成形时的摩擦系数进行了测试实验。计算结果与实验数据表明，板料成形过程中的摩擦系数不是一个常数，随着凸模行程的增加，摩擦系数具有增加的趋势。     

Modeling of anisotropic plastic behavior of ferritic stainless steel sheet

Modeling of anisotropic plastic behavior of ferritic stainless steel sheet (Type 409) was investigated using the three yield functions of Hill [A theory of the yielding and plastic flow of anisotropic metals. Proceedings of Royal Society of London, Series A 1948;193:281–97.], Barlat and Lian [Plastic behavior and stretchability of sheet metals. Part I: A yield function for orthotropic sheets under plane stress conditions. International Journal of Plasticity 1989;5:51–66] and Barlat et al. [Plane stress yield function for aluminum alloy sheet. Part I: Theory. International Journal of Plasticity 2003;19:1297–319.] (referred to as Yld2000-2d) criteria. Mechanical behaviors were characterized based on uniaxial tension, balanced biaxial bulge, and disk compression tests. Directionalities of yield stresses and r values were predicted from the three criteria and compared with experimental results. In order to verify the modeling accuracy of the three functions under complex loading conditions, cylindrical cup drawing and limiting dome height tests were carried out numerically and experimentally. It has been demonstrated that the result from Yld2000-2d criterion exhibits good agreement with experimental data. The effects of anisotropic hardening on earing and necking were also investigated based on the different levels of plastic work.     

圆锥形凹模径向分块压边拉深工艺实验研究

采用圆锥形凹模拉深工艺可以提高成形极限,但需要用压边圈将板坯先压成与凹模面吻合的形状,当变形程度较大时,板坯很容易起皱。为了克服这一缺点,提出了将圆锥形凹模与径向分块压边方法结合的工艺,该工艺可有效改善压边圈与板坯的约束状态,从而达到抑制起皱的目的。对圆筒形件的拉深成形,采用了刚柔复合的径向分块压边圈结构,设计了圆锥形凹模径向分块多压边圈拉深模,取不同凹模半锥角的圆锥形凹模进行了圆筒形件的拉深成形实验。实验表明,新的压边方法能有效克服初始成形过程的起皱,可与锥度较小的凹模一起使用。采用凹模半锥角为45°的凹模,得到AA5754、AA6061和08Al三种板材的极限拉深系数分别为0.410、0.431、0.373,显著提高了成形极限。对圆锥形凹模的拉深成形,给出了理论计算成形极限的方法,理论与实验结果非常接近。     

板料拉延过程的数值分析

本文采用有限应变理论中增量形式的total Lagrangian列式和J 2流动法则,用薄膜壳模型对板料拉延过程进行了数值分析。板料的应力-应变关系为指数形式,考虑了材料厚向异性的影响,设板料与模具的接触表面适用库仑摩擦定律。分析中考虑了压边的作用,引入了凸缘起皱和简部破裂的判断准则。计算结果与实验数据进行了比较。     

平面应力条件下轴对称拉深成形法兰区起皱和破裂分析

分别采用平面应力和平面应变假设条件,对轴对称拉深成形法兰区的应力分布进行了分析比较,两种情况下的径向应力计算值相差较小,但周向应力计算值相差较大。有限元模拟表明,平面应力条件下得到的解析结果与模拟值非常接近,表明平面应力假设条件比平面应变假设条件更接近于实际情况。在平面应力条件下,建立了轴对称成形法兰区起皱失稳条件和圆筒形件破裂失稳条件,导出了临界压边力的计算式。对于某一具体的拉深成形问题,计算得到了临界压边力与拉深位置的关系曲线。分析结果有助于进一步认识拉深成形各变形区的变形特点,建立更符合实际的起皱、破裂准则,预测成形缺陷及临界压边力。