关于1561铝合金曲面件的多点成形工艺
付文智1, 刘晓东1, 王洪波2, 闫德俊2, 刘晓莉2, 李明哲1, 董玉其2, 曾振华2, 刘桂彬2
1.吉林大学 无模成形技术开发中心,长春 130022
2.中船黄埔文冲船舶有限公司,广州 511464

作者简介:付文智(1960-),男,教授,博士.研究方向:无模成形及数字化制造.E-mail:fwz@jlu.edu.cn

摘要

为了掌握新型舰船用材料1561铝合金板类件的成形规律,通过拉伸试验获得了1561铝合金的机械性能特性曲线。通过有限元软件模拟了1561铝合金的多点成形和卸载过程,结果表明:1561铝合金板料在成形过程中易起皱,具有较大的回弹量。采用回弹补偿法和多道次成形工序对某船用曲面件进行多点成形试验,并通过三维扫描和Geomagic软件对成形精度进行了分析,结果表明:采用回弹补偿法和多道次成形方法可以很好地解决曲线的回弹和起皱问题。

关键词: 塑性加工工艺与设备; 1561铝合金; 数值模拟; 回弹补偿; 多点成形; 多道次成形
中图分类号:TG386 文献标志码:A 文章编号:1671-5497(2017)06-1822-07
Multi-point forming process of 1561 aluminum alloy surfaces
FU Wen-zhi1, LIU Xiao-dong1, WANG Hong-bo2, YAN De-jun2, LIU Xiao-li2, LI Ming-zhe1, DONG Yu-qi2, ZENG Zhen-hua2, LIU Gui-bin2
1.Dieless Forming Technology Center, Jilin University, Changchun 130022, China
2.CSSC Huangpu Wenchong Shipbuilding Company Limited,Guangzhou 511464, China
Abstract

In order to know the formability of 1561 aluminum alloy sheet as a new ship material,its mechanical properties were obtained by tensile test. The multi-point forming and unloading process of aluminum alloy sheet was simulated by finite element software ABAQUS. Simulation results showed that 1561 aluminum alloy material was easy to wrinkle during the forming process and had a large springback. Three-dimensional scanning and Geomagic analysis were carried out, and the results verify the feasibility of using the springback compensation to compensate the springback and multi-steps forming method to solve the winkling in the multi-point forming of 1561 aluminum alloy. At the same time, it is shown that simulation results are reliable.

Keyword: plastic processing technology and equipment; 1561 aluminum alloy; numerical simulation; springback compensation; multi-point forming; multi-pass forming
0 引 言

随着船舶制造向轻量化、高速化、大型化方向发展, 具有小密度、大比强度、较好稳定性的铝镁合金成为替代船用钢材的优选材料[1, 2]。目前在船舶制造中5083铝合金使用得较多, 但焊接后的矫形处理对制造精度的影响较大[3]。1561铝合金(俄罗斯牌号)[4]相比5083铝合金具体更高的强度, 采用1561铝合金作为船用材料将会提高舰船的质量和性能。多点成形[5, 6, 7]是一种适用于复杂件成形的柔性成形技术[8, 9, 10, 11], 它是由一系列紧密排列、高度可调的基本体单元组成[12]。相比于传统模具成形, 一套离散化的多点模具通过数字化控制即可成形多种形状和尺寸的曲面件, 省去了传统模具制造、调试等成本, 真正实现了“ 一机多用” 的构想, 适用于曲面件的板料成形。

目前, 国内、外对1561铝合金成形性的研究鲜有报道, 为了迎合新产品开发的需要, 本文利用有限元软件对某船用典型曲面件进行了多点成形数值模拟, 研究了1561铝合金的回弹补偿, 掌握了1561铝合金成形件的回弹规律, 并对模拟结果进行了试验验证, 结合3DSS三维扫描和Geomagic软件分析对成形件的成形精度进行了分析。

1 1561铝合金的拉伸试验

为了使多点成形过程模拟计算得到的结果更精确, 使计算结果能够指导实际生产, 按照《GB/T228.1— 2010》[13]设计试样, 如图1所示, 在电子万能试验机上进行室温拉伸断裂试验。表1为1561铝合金的拉伸试验数据, 可见, 1561铝合金的抗拉强度和屈服强度的平均值分别为345.9和266.67 MPa; 弹性模量E为67.516 GPa; 延伸率平均值为21.8%。1561铝合金的真实应力-应变曲线如图2所示。

图1 拉伸试样Fig.1 Tensile test specimen

表1 1561铝合金力学性能试验数据 Table 1 1561 aluminum alloy mechanical properties test data

图2 1561铝合金的应力-应变曲线Fig.2 Stress-strain curve of 1561 aluminum alloy

2 多点成形有限元模型

为了更好地分析1561铝合金板料的多点成形过程, 利用有限元分析方法中的动态显式算法对1561铝合金的加载成形过程进行计算, 同时结合静态隐式算法对卸载回弹过程进行数值模拟。图3为多点成形有限元模型, 由板料、聚氨酯弹性垫以及上、下多点模具组成。成形板料选用厚度为3 mm的1561铝合金。多点模具形面尺寸为500 mm× 700 mm; 基本体单元为10 mm× 10 mm; 冲头半径为8 mm, 板料形状如图3所示。与板料变形量相比, 冲头变形量非常小, 故冲头设为刚体单元。同时为了简化计算模型, 假设板料为各向同性。上、下多点冲头采用R3D4单元划分网格, 网格单元尺寸为2 mm。板料采用S4R壳单元划分网格, 网格单元尺寸为3 mm。聚氨酯弹性垫采用C3D8R单元划分网格, 网格单元尺寸为3 mm。聚氨酯弹性垫与板料和多点模具之间的摩擦因数取0.3。

图3 多点成形有限元模型Fig.3 Multi-point forming finite element model

3 数值模拟结果
3.1 成形情况分析

图4为本文开发成形的某船外板零件所需板料。为了叙述方便, 设板料宽度中心长度方向为X方向; 长度中心宽度方向为Y方向; 板料平面法线方向为Z方向, 如图4所示。

图4 成形目标件所需板料Fig.4 Sheet required to form target

图5为成形得到的某船外板多点成形件回弹前、后的轮廓对比图, 其中网格图为回弹前工件轮廓图。图6为成形件卸载前、后沿XY方向的高度值, 可以发现, 1561铝合金在卸载前、后回弹较大, 前、后形状差别较大, 其中Y方向形状波动起伏不光滑, 发生了起皱, X方向曲面过渡平顺。由此可得, 零件在Y方向较难成形, 皱纹沿着X方向。

图5 无回弹补偿量工件回弹前后比照图Fig.5 Sheet required to form target

图6 成形件的高度Fig.6 Height of forming part

由此可见, 1561铝合金的成形件具有较大的回弹量, 且与目标件在XY方向的回弹量不同, 因此应采用差异补偿法, 令XY方向的形面补偿系数不同。回弹补偿法是通过对模具形面进行预修正, 对模具施加与回弹方向相反的补偿量, 使卸载后的曲面接近目标形状。

对于双曲率零件, 由于零件各个方向曲率不一致, 导致XY方向的回弹量不同, 因此XY方向的补偿量不同。补偿面上任一位置高度Z的确定方法为:假设目标曲面件上任一点的高度值Z'h, X方向型面补偿系数为m, Y方向型面补偿系数为n, 则零件的回弹补偿量计算公式为:

Z=(1+m)h, 位置在X轴上(1+n)h, 位置在Y轴上(1+m+n)h, 位置既不在 X轴也不在Y轴上1

3.2 成形路径对成形件的影响

合理的变形路径能够使变形板料受力均匀, 变形均匀, 而且能够避免成形缺陷的产生, 获得成形效果良好的零件。板料沿理想路径发生变形, 任意时刻t的形变为[14]:

U(X, t)=Q(t)· Diag[λ α (t)(t)]· QT(t) (2)

式中:Q(t)为正交旋转张量; Diag[λ α (t)(t)]为Hencky应变的主值组成的对角张量, λ α (t)(t)为U(X, t)的主值, α (t)为变形参量。

通过建立泛函方程, 采用数值方法构造的中间构形为:

tx* =argminψ(tx)=argmintΩ12(U-U^):(U-U^)dtΩ3

式中:“ :” 为张量计算符号。

构造的中间构形的单元节点坐标计算公式为:

ψ(Xnode)Xnode=etΩeU-U^):uXnodedtΩe=0(4)

式中:tΩ et时刻构形的体积; Xnode为初始构形中单元的节点坐标; U、分别为数值方法计算出的和理想路径的中间构形变形张量。

通过式(4)可得到t时刻的中间构形。因初始构形及中间构形都是按照理想路径设计, 因此通过式(4)可得到最优成形路径。通过选取变形量不同的最优成形路径上的中间构形, 并通过数值仿真计算, 可得到多道次成形的优化成形路径, 最终确定的目标件成形路径依次为X0.5Y1.2、X0.7Y1.2、X0.9Y1.4、X1.0Y1.6。

采用差异补偿法补偿多点型面, 并进行一次成形和四道次成形, 其成形对比如图7所示; 零件卸载前、后XY方向的高度值如图8所示。图7(a)所示的一次成形零件中部发生起皱, 而图7(b)所示的四道次成形零件过渡平顺, 质量良好, 皱纹完全消失。这是由于多点模具成形采用一次成形方法时, 材料一次性变形较大, 增加回弹补偿量后, 成形深度增加, 容易超过材料一次成形所允许的最大变形量, 在无压边力作用下成形较易起皱。利用多点模具高度可调的优点, 采用多道次成形方法, 每次成形深度较小, 使板料各部分变形在变形过程中最大限度地均匀化, 最终达到目标件尺寸, 可得到质量较高的目标件。由图8可以看出:

图7 不同路径下成形件光照效果图Fig.7 Shaping effect of forming parts under different paths

图8 不同成形路径下成形件沿中线的高度值Fig.8 Height values of forming part along center line under different forming paths

成形件在卸载回弹后在Y方向的回弹量较大, 在X方向回弹量较小, 说明, Y方向比X方向较难成形, 需要在Y方向给予更多的补偿量。采用多道次成形方法, 板料每次变形相对一次变形较小, 相对弯曲半径R/t较小, 且变形充分、均匀, 相对于一次成形具有较小的回弹量。

图9为根据模拟结果多次补偿模具型面选取的成形路径依次为X0.5Y1.2、X0.7Y1.2、X0.9Y1.4、X1.0Y1.6四道次成形方法成形零件得到的X方向和X1.0Y1.6四道次成形方法成形零件得到的X方向和Y方向高度值。可以看出, 对于1561铝合金这种回弹较大的材料, 当成形件的变形量较小时, X方向几乎没有变形, 而Y方向轻微变形。采用多道次逐渐过渡的方法成形, 由于每次变形量较小不超过板料成形所允许最大变形量, 成形件表面过渡光顺, 无起皱现象发生, 可得到质量良好的成形件。

图9 四道次成形件沿中线的高度值变化Fig.9 Change of height value along center line in multi-pass forming

4 试验验证

对材料为1561铝合金的某船用外板零件在多点成形压力机上进行了多道次成形, 成形件公差要求为± 1 mm。多点成形试验的现场照片如图10所示。

图10 多点成形试验Fig.10 Multi-point forming experiment

图11 使用样板对成形件进行测量Fig.11 Molded parts were measured using template

板料尺寸、基本体尺寸与数值模拟相同。成形结果如图11所示, 可以看出, 采用多道次成形工序, 得到的成形件无压痕和起皱等缺陷, 表面连续光滑, 与数值模拟结果吻合。通过样板对成形件检测, 可以看出:形件在X方向和Y方向与目标尺寸吻合较好, 误差较小。以上验证了在1561铝合金多点成形时, 采用回弹补偿法补偿回弹和多道次工艺解决形件起皱的可行性, 也验证了有限元分析结果对生产实践指导的可靠性。

使用3DSS三维扫描仪对一道次和四道次多点成形件进行光学扫描, 获取成形件表面的点云数据; 使用Geomagic分析软件对成形件点云数据与目标件进行比较, 可得到成形件的误差分布规律。成形件的形状误差云图和误差分布图分别如图12和图13所示。由图可以看出, 一次成形件误差80%分布在-1~1 mm, 95%分布在-2~2 mm; 而四道次成形件误差91%分布在-1~1 mm, 96%分布在-2~2 mm。相较于一次成形, 四道次成形件误差较小, 成形质量高, 符合产品要求, 进一步验证了数值模拟的可靠性和采用多道次工序成形1561铝合金的可行性。

图12 形状误差云图Fig.12 Shape error cloud

图13 成形件误差分布Fig.13 Error distribution of forming parts

5 结 论

(1)新型材料1561铝合金具有较高的屈服强度和较小的弹性模量, 在多点成形后回弹较大, 小变形量在卸载后可基本完全回弹, 因此, 1561铝合金成在成形过程中需给予较大的回弹补偿量。

(2)由数值模拟结果可以看出, 当成形件曲率较大时, 1561铝合金在一次成形时易起皱, 而采用多道次成形方法后, 解决了成形件中心部位起皱的问题, 成形误差± 1 mm的区域由80%提高到91%, 误差减小, 成形精度提高。

(3)通过试验对1561铝合金多点成形的可行性及数值模拟的可靠性进行了验证, 试验结果对1561铝合金零件的实际成形有一定的指导意义。

The authors have declared that no competing interests exist.

参考文献
[1] 闫德俊, 韩端锋, 王毅, . 1561铝合金双面双弧TIG焊接接头的组织和力学性能[J]. 中国有色金属学报, 2016, 26(10): 2065-2070.
Yan De-jun, Han Duan-feng, Wang Yi, et al. Microstructure and mechanical properties of double-sided double-arc TIG welded joint of 1561 aluminum alloy[J]. Chinese Journal of Nonferrous Metals, 2016, 26(10): 2065-2070. [本文引用:1]
[2] 闫德俊, 韩端锋, 王毅, . 1561铝合金TIG深熔焊接头组织与力学性能[J]. 焊接学报, 2016, 37(12): 29-32.
Yan De-jun, Han Duan-feng, Wang Yi, et al. Microstructure and mechanical properties of TIG Deep fused joints of 1561 aluminum alloy[J]. Transactions of the China Welding Institution, 2016, 37(12): 29-32. [本文引用:1]
[3] 韩善果, 闫德俊, 刘晓莉, . 1561铝合金TIG焊接接头组织与力学性能分析[J]. 焊接技术, 2015, 44(1): 18-21.
Han Shan-guo, Yan De-jun, Liu Xiao-li, et al. Microstructure and mechanical properties of TIG welded joints of 1561 aluminum alloy[J]. Welding Technology, 2015, 44(1): 18-21. [本文引用:1]
[4] 韩善果, 郑世达, 蔡得涛, . 双丝CMT焊接参数对1561铝合金焊接接头力学性能和组织的影响[J]. 焊接, 2016(5): 31-34.
Han Shan-guo, Zheng Shi-da, Cai De-tao, et al. Effect of double wire CMT welding parameters on mechanical properties and microstructure of1561 aluminum alloy welded joints[J]. Welding & Joining, 2016(5): 31-34. [本文引用:1]
[5] Gong X P, Li M Z, Hu Z Q, et al. Research on continuous multi-point forming technology for three-dimensional sheet metal[J]. International Journal of Materials and Product Technology, 2010, 38(2/3): 210-222. [本文引用:1]
[6] 谭富星, 李明哲, 钱直睿, . 钛合金颅骨修复体的多点成形数值模拟[J]. 吉林大学学报: 工学版, 2006, 36(6): 851-855.
Tan Fu-xing, Li Ming-zhe, Qian Zhi-rui, et al. Numerical simulation on process of multi-point forming for titanium alloy cranial prosthesis[J]. Journal of Jilin University (Engineering and Technology Edition), 2006, 36(6): 851-855. [本文引用:1]
[7] Li M Z, Liu Y H, Su S Z, et al. Multi-point forming: a flexible manufacturing method for a 3-D surface sheet[J]. Journal of Material Processing Technology, 1999, 87(1-3): 277-280. [本文引用:1]
[8] Qiu Ning-jia, Li Ming-zhe, Li Ren-jia. The shape analysis of three-dimensional flexible rolling method[J]. Journal of Engineering Manufacture, 2016, 230(4): 618-628. [本文引用:1]
[9] Shim D S, Yang D Y, Chung S W, et al. Optimization of forming steps in the incremental forming of twisted shapes using a line array roll set (LARS)[J]. International Journal of Precision Engineering and Manufacturing, 2010, 11(5): 715-723. [本文引用:1]
[10] Li R J, Li M Z, Qiu N J, et al. Surface flexible rolling for three-dimensional sheet metal parts[J]. Journal of Materials Processing Technology, 2014, 214(2): 380-389. [本文引用:1]
[11] 付文智, 王妤, 李明哲, . 柔性轧制成形中曲面件起皱的数值分析与控制方法[J]. 光学精密工程, 2016, 24(10): 407-414.
Fu Wen-zhi, Wang Yu, Li Ming-zhe, et al. Numerical analysis and control method of curved surface wrinkles in flexible rolling forming[J]. Optics and Precision Engineering, 2016, 24(10): 407-414. [本文引用:1]
[12] 龚学鹏, 李明哲, 卢启鹏, . 连续多点成形中的成形载荷分析[J]. 光学精密工程, 2012, 20(6): 1288-1295.
Gong Xue-peng, Li Ming-zhe, Lu Qi-peng, et al. Analysis on forming load in continuous multi-point forming[J]. Optics and Precision Engineering, 2012, 20(6): 1288-1295. [本文引用:1]
[13] GB/T228. 1-2010金属材料拉伸试验第1部分: 室温试验方法[S]. [本文引用:1]
[14] 裴永生, 李明哲, 蔡忠义, . 板材变路径多点成形的理论分析与实现[J]. 农业机械学报, 2003, 34(2): 114-116.
Pei Yong-sheng, Li Ming-zhe, Cai Zhong-yi, et al. Theoretical analysis and realization of multi-point forming of plate-changing path[J]. Journal of Agricultural Mechanization, 2003, 34(2): 114-116. [本文引用:1]