Influence of Anisotropic Yield Model on Simulation Results of Deep Drawing Forming of Pickling Steel Plate
-
摘要:
采用Autoform软件分别基于Hill' 48,Barlat' 89和BBC-2005各向异性屈服模型对SRTC1酸洗钢板拉深成形进行有限元模拟,通过制耳高度和厚度分布数值模拟与试验结果的对比,研究了屈服模型对酸洗钢板拉深成形模拟的影响,获得了最佳屈服模型。结果表明:采用Hill' 48,Barlat' 89和BBC-2005屈服模型模拟的钢板拉深成形筒形件制耳高度差与试验结果的相对误差分别为16.20%,11.27%,7.74%,厚度最大差值模拟结果与试验结果的相对误差分别为8.28%,5.25%,2.30%。BBC-2005屈服模型综合考虑了材料应力和变形各向异性参量,在制耳高度和厚度分布方面均最有最高的预测精度。
Abstract:The finite element simulation of SRTC1 pickling steel plate deep drawing was carried out based on Hill' 48, Barlat' 89 and BBC-2005 anisotropic yield models by using Autoform software, respectively. And the numerical simulation of ear height and thickness distribution were compared with the test results. The influence of yield model on the simulation of deep drawing forming of pickling steel plate was studied, The optimal yield model was obtained. The results show that the relative errors between ear height difference of the deep drawing forming steel plate cylindrical part simulated by Hill' 48, Barlat' 89 and BBC-2005 yield models and test results were 16.20%, 11.27% and 7.74%, respectively. The relative errors of the maximum thickness difference between the simulated results and the test results were 8.28%, 5.25% and 2.30%, respectively. The BBC-2005 yield model had the highest prediction accuracy in both the ear height and thickness distribution because of the consideration of the stress and deformation anisotropy parameters of the material.
-
Keywords:
- deep drawing /
- anisotropy /
- yield model /
- numerical simulation
-
0. 引言
热轧酸洗钢板(简称酸洗板)的表面质量和使用性能介于热轧钢板和冷轧钢板之间,因其性价比较高,广泛应用于汽车和家电领域。但酸洗板的各向异性显著,拉深制耳问题作为其常见缺陷,成为热轧酸洗板工业应用的重要障碍[1]。近年来,诸多研究人员以典型筒形件为例,对拉深成形过程中的金属流动规律及其影响因素等进行了研究[2-4]。林艳丽等[5]建立了制耳形状与板坯面内各向异性系数分布规律的定量解析关系,对满足Hill48屈服模型的板材进行了分析,预测了制耳出现的位置和数量。胡成武等[6]以DC04冷轧钢板为研究对象,采用Hill´ 48屈服模型研究了筒形件变间隙拉深过程中压边力对拉深制耳的影响,获得了变间隙拉深过程中的压边力-制耳率曲线。孙雅倩[7]基于三参数Barlat´ 89屈服模型建立了SPCC冷轧钢板筒形件拉深仿真模型,发现制耳率与板料直径呈正相关,与压边力呈负相关。目前针对冷轧钢板拉深变形行为,学者们大多选择固定的屈服模型,通过有限元仿真分析研究冲压工艺参数、模具参数等对拉深成形的影响,对于热轧酸洗钢板的拉深变形行为,相关研究鲜有报道。
不同屈服模型对于板材拉深成形金属流动规律的预测精度不尽相同,选择合适的屈服模型是描述板材拉深变形行为的前提条件。为此,作者针对某钢厂生产的SRTC1酸洗钢板开展了拉深成形数值模拟与试验研究,采用不同各向异性屈服模型模拟了制耳高度和厚度分布,通过与试验结果对比确定了最佳屈服模型,以期为实际零件的成形工艺研究分析提供理论参考。
1. 典型的各向异性屈服模型
Autoform有限元分析软件[8]被广泛应用于汽车行业中的冲压成形仿真,其内含多个用于描述钢材力学行为的各向异性屈服模型,如常用的Hill´ 48[9],Barlat´ 89[10]以及BBC-2005[11]等。
1.1 Hill' 48屈服模型
1948年,英国学者Hill首次将材料各向异性参数引入到屈服准则中,提出了正交各向异性材料的屈服模型Hill´ 48,为板料塑性变形各向异性理论的建立奠定了基础。Hill´ 48屈服模型的函数表达式为
(1) 式中:τ为剪切应力;
为等效应力;x,y,z为材料各向异性主轴;σxx为x轴应力,以此类推;F,G,H,L,M,N为彼此独立的各向异性系数,与材料力学性能密切相关,由试验确定。 通常情况下,板料在成形过程中大多处于平面应力状态,即σzz,σyz,σzx均为0,则式(1)可简化为
(2) Hill´ 48屈服模型充分考虑了材料的各向异性且囊括剪切应力元素,适用于二维和三维问题,其函数形式简单,参数求解便捷,因此被广泛应用于解决实际工程问题。许多有限元软件都包含该模型,如Dynaform、Abaqus等软件。
1.2 Batlat' 89屈服模型
1989年,法国学者Barlat等为了更好地研究板料的各向异性,提出了各向异性屈服模型Barlat´ 89。Barlat´ 89屈服模型在平面应力条件下建立,函数表达式为
(3) (4) (5) 式中:a1,h1,c1,p1为相互独立的参数;m1为非二次屈服函数的指数;K1,K2代指公式(4),(5);x和y为材料各向异性主轴。
对于面心立方材料,m1=8;对于体心立方材料,m1=6;当m1=2时,Barlat´ 89屈服模型转化为Hill´ 48屈服模型。Barlat´ 89屈服模型能合理地描述具有较强织构的各向异性金属板料的屈服行为,有效模拟板料拉深成形过程中法兰的塑性流动规律,但由于该屈服模型参数由3个各向异性系数确定,对屈服应力的描述不是很准确,且只能用于平面应力状态。
1.3 BBC-2005屈服模型
德国学者Banabic等在上述屈服模型的基础上提出了一系列BBC屈服模型,其中较为典型的是BBC-2005屈服模型,其数学表达式为
(6) (7) (8) (9) 式中:λ为可调数,可根据材料的实际情况进行调节;Λ,Γ,Ψ为应力张量平面应力分量的函数;a2,R,S,T,P,Q为相互独立的各向异性参数,可由材料3个单向拉伸状态下的应力(σ0,σ45,σ90;0,45,90代表与轧制方向所成角度),塑性应变比r(r0,r45,r90)以及等双拉状态下的应力σb和塑性应变比rb通过数值计算的方法获得。
BBC-2005屈服模型具有较大的柔性,适用于钢、铝等材料的塑性成形。
2. 试验方法
2.1 拉伸试验
试验材料为首钢生产的厚度1.65 mm的SRTC1酸洗钢板。根据GB/T 228—2021,在试验钢板上分别沿轧制方向(0°方向)、与轧制方向成45°角(45°方向)和垂直轧制方向(90°方向)制取标距长度为80 mm的“哑铃型”试样,采用Zwick-Z100型拉伸试验机进行准静态拉伸试验,应变速率为0.001 s−1。由拉伸试验得到如图1所示的真应力-真应变曲线,对曲线进行拟合计算,获取硬化指数,计算公式为
(10) 式中:σl为流动应力;εp为塑性应变;k为材料参数;n为硬化指数。
根据GB/T 5027—2016,基于体积不变原理,通过测试试样长度和宽度方向的变形量,计算塑性应变比,公式为
(11) 式中:b0,b为试样拉伸前后的宽度;L0,L为试样拉伸前后的标距长度。
计算得到拉伸性能参数见表1,可知SRTC1酸洗钢板3个方向的力学性能存在明显差异,各向异性较为显著。
表 1 SRTC1酸洗钢板不同方向的拉伸性能参数Table 1. Tensile performance parameters of SRTC1 pickling steel plate in different directions方向/(°) 屈服强度/MPa 抗拉强度/MPa n r 0 335 444 0.15 0.84 45 360 447 0.13 1.08 90 365 454 0.14 0.99 2.2 拉深成形试验
采用BSC-400型高强汽车板综合成形试验机进行室温拉深成形试验,试样为直径90 mm的圆板料,凸模直径为50 mm,压边圈外径为70 mm,压边圈内径为54 mm,凹模内径为54 mm,凸模圆角半径为5 mm,凹模圆角半径为10 mm。将板料涂抹好润滑剂后放置于压边圈之上,压边圈带动板料向上运动与凹模压紧,然后凸模以0.5 mm·s−1的速度向下运动,压边力恒定为100 kN。
采用千分尺沿着与轧制方向成0°至90°方向每隔10°测量顶部到底部距离获取制耳高度,将试样沿径向切开,从筒壁边缘向筒底中心取点,测试径向厚度分布。
3. 拉深成形数值模拟
采用Autoform有限元软件对拉深成形过程进行仿真模拟,模型如图2所示,网格单元类型设置为EPS-11,仿真工艺参数与拉深试验保持一致,凸模运动速度为0.5 mm·s−1,压边力为100 kN,摩擦因数为0.12。根据SRTC1酸洗钢板拉伸性能参数制作仿真数据卡片,分别应用Hill´ 48、Barlat´ 89和BBC-2005等3种屈服模型进行模拟。
采用Autoform软件后处理模块中“Vector Measurement”功能,选取模型变形后底部网格节点与顶部边缘网格节点,读取垂直距离获取制耳高度;获取试样径向厚度分布,从筒壁边缘向筒底中心提取节点厚度数据。
4. 拉深成形试验与模拟结果对比
4.1 制耳高度对比
图3(d)中的RD表示轧制方向。由图3可知,3种屈服模型都预测得到SRTC1酸洗钢板拉深成形后产生了4个制耳,与试验结果相似。由图4可见:基于Hill´ 48和Barlat´ 89屈服模型预测的制耳高度低于试验结果,且预测的制耳峰(制耳高度峰)位置与试验结果有明显偏差,基于BBC-2005屈服模型预测的制耳高度及峰值位置与试验吻合度较高。Hill´ 48,Barlat´ 89,BBC-2005模型预测得到的制耳高度差(制耳峰与制耳谷垂直高度差)分别为1.19,1.26,1.31 mm,与试验结果(1.42 mm)的相对误差分别为16.20%,11.27%,7.74%。说明BBC-2005模型对材料流动状态的预测更符合实际结果,这是因为该屈服模型同时考虑了材料的应力以及变形各向异性参量。
4.2 厚度分布对比
由图5可见:3种屈服模型模拟的径向厚度分布趋势与试验结果保持一致,由筒壁至筒底,酸洗钢板先增厚再减薄再增厚,最薄位置为筒形件底部圆角处,此处模拟厚度值稍高于试验结果。Hill´ 48,Barlat´ 89,BBC-2005模型预测得到的筒形件厚度最大差值(厚度最大值与最小值的差)分别为0.588,0.571,0.555 mm,与试验结果(0.543 mm)的相对误差分别为8.28%,5.25%,2.30%。对比可知BBC-2005模型对于材料流动预测更为准确,更能全面描述材料的力学行为。结合前文不同屈服模型对于拉深制耳预测结果来看,SRTC1酸洗钢板成形模拟选择BBC-2005屈服模型更合适。
5. 结论
(1)单向拉伸试验测得SRTC1酸洗钢板沿轧制方向、与轧制方向呈45°和垂直于轧制方向的屈服强度和塑性应变比存在明显差异,表明SRTC1酸洗钢板具有显著的各向异性。
(2)Hill´ 48,Barlat´ 89和BBC-2005屈服模型模拟的SRTC1酸洗钢板拉深成形筒形件的制耳高度差与试验结果的相对误差分别为16.20%,11.27%,7.74%,径向厚度最大差值模拟结果与试验结果的相对误差分别为8.28%,5.25%,2.30%。相比之下,BBC-2005屈服模型综合考虑了材料应力和应变各向异性参量,在制耳高度和厚度分布方面均有最高的预测精度。
-
表 1 SRTC1酸洗钢板不同方向的拉伸性能参数
Table 1 Tensile performance parameters of SRTC1 pickling steel plate in different directions
方向/(°) 屈服强度/MPa 抗拉强度/MPa n r 0 335 444 0.15 0.84 45 360 447 0.13 1.08 90 365 454 0.14 0.99 -
[1] 董学涛,金世铨,王亚东,等. SPHE压缩机壳体冲压制耳缺陷分析[J]. 金属世界,2021(3):45-48. DONG X T ,JIN S Q ,WANG Y D ,et al. Failure analysis of earing defect of SPHE compressor steel[J]. Metal World,2021(3):45-48.
[2] 仇建桐,邓沛然,邵威,等. 基于Dynaform的铝合金筒形件拉深成形[J]. 锻压技术,2020,45(5):49-55. QIU J T ,DENG P R ,SHAO W ,et al. Deep drawing of cylindrical parts for aluminum alloy based on Dynaform[J]. Forging and Stamping Technology,2020,45(5):49-55.
[3] 陈旭升AZ31镁合金筒形件拉深成形规律的研究秦皇岛燕山大学2021陈旭升. AZ31镁合金筒形件拉深成形规律的研究[D]. 秦皇岛:燕山大学,2021. CHEN X SStudy on drawing forming law of AZ31 magnesium alloy cylindrical partsQinhuangdaoYanshan University2021CHEN X S. Study on drawing forming law of AZ31 magnesium alloy cylindrical parts[D]. Qinhuangdao:Yanshan University,2021.
[4] 范豇宇,刘帅,占乐,等. 高强度钢筒形件拉深工艺研究[J]. 燕山大学学报,2021,45(3):202-207. FAN J Y ,LIU S ,ZHAN L ,et al. Study on cylinder deep drawing of high strength steel[J]. Journal of Yanshan University,2021,45(3):202-207.
[5] 林艳丽,何祝斌,袁杭,等. 金属板材拉深制耳形状与各向异性系数的解析关系模型[J]. 塑性工程学报,2021,28(6):19-25. LIN Y L ,HE Z B ,YUAN H ,et al. Analytical relationship model of ear profile and anisotropy coefficient of sheet metal deep drawing[J]. Journal of Plasticity Engineering,2021,28(6):19-25.
[6] 胡成武,李光,毛远征,等. 压边力对变间隙拉深制耳的影响[J]. 塑性工程学报,2019,26(2):132-137. HU C W ,LI G ,MAO Y Z ,et al. Influence of blank holder force on variable gap deep drawing earing[J]. Journal of Plasticity Engineering,2019,26(2):132-137.
[7] 孙雅倩. 圆筒形件拉深过程中凸耳行为研究[J]. 锻压装备与制造技术,2019,54(5):73-77. SUN Y Q. Study on the behavior of earring during deep drawing process of cylindrical part[J]. China Metalforming Equipment and Manufacturing Technology,2019,54(5):73-77.
[8] 金磊,王宝川,韩龙帅,等. 基于AutoForm的材料仿真数据影响因素研究[J]. 锻压技术,2019,44(9):180-186. JIN L ,WANG B C ,HAN L S ,et al. Research on influencing factors of material simulation data based on AutoForm[J]. Forging and Stamping Technology,2019,44(9):180-186.
[9] HILL R. A theory of the yielding and plastic flow of anisotropic metals[J]. Proceedings of the Royal Society of London Series A:Mathematical and Physical Sciences,1948,193(1033):281-297. [10] BARLAT F ,LIAN K. Plastic behavior and stretchability of sheet metals. part I:A yield function for orthotropic sheets under plane stress conditions[J]. International Journal of Plasticity,1989,5(1):51-66. [11] BANABIC D ,ARETZ H ,COMSA D S ,et al. An improved analytical description of orthotropy in metallic sheets[J]. International Journal of Plasticity,2005,21(3):493-512.