0 引 言

随着生活水平与科学技术的日益发展，人们对汽车的质量与外观上的要求与日俱增，进而随之而来的，涉及金属板材的应用领域也愈来愈广泛，而且对于复杂件成形的生产要求亦越来越高。汽车覆盖件[1-2]是汽车组成的重要部分，而且它与一般冲压件相比，具有薄壁结构，形状复杂，结构尺寸大，表面质量，以及尺寸精度要求高，以及生产成本高等特点。而其中翼子板[3-5]是汽车覆盖件中成形困难的典型零件，由于它不仅具有外观要求质量高，尺寸精度要求高，模具设计与制造成本高等特点，而且需要经过拉延，切边，翻边和整形等多道工序，因此在成形过程中很容易出现起皱，破裂与变形不足等缺陷。传统的冲压成形[6-9]由于自身工艺特点所限，零件缺陷较多，同时所需费用较高，而充液拉深[10-13] (hydromechanical deep drawing, HDD)具有拉深比大、成形零件表面质量高、尺寸精度高、以及能够进行复杂零件的成形等优点[14-16]，所以充液成形是翼子板成形的有效途径，能够很好地符合翼子板的成形需求。

小狗包弟的故事发生在“文革”期间，巴金写作此文是“文革”结束之后。“文革”对学生来讲只是个概念，是历史教科书上短短的几行字；学生对此既没有感性认知，也没有理性认识。时代背景对学生来说无疑是一个盲点，也是理解文本的一个难点。如果不了解时代背景，就不可能很好地理解文本，通过阅读文本来认识巴金反思的价值和意义也就无从谈起。文本本身并不难，而且文本中就隐含着时代背景的大量信息，因此，我就采用了通过对问题的思考、对文本的解读来了解时代背景的方式。

以汽车前翼子板为研究对象，在零件成形工艺性分析的基础上，采用有限元分析软件对翼子板进行充液拉深模拟，以零件的成形极限图(FLD)为分析标准，分析了拉延阻力系数、液室压力以及压边力对翼子板成形质量的影响，再选取翼子板一截面，以壁厚分布为分析依据，进一步研究压边力对翼子板成形性能的影响规律。为覆盖件充液成形工艺的进一步研究提供了理论基础和生产依据。

1 零件成形工艺性分析

图1是某汽车前翼子板的三维模型，从图中可以看出翼子板与引擎盖配合处的形状较为复杂且存在台阶结构；与门板配合处拉延深度较浅，形状平坦，拉延时容易产生拉延不充分现象，同时在轮罩处存在凸起和翻边。按照传统加工方法很难实现一次拉深成形，需要经过拉延、切边、翻边和整形等多道工序，会耗费时间和材料，同时存在成形精度不足等问题，而采用充液拉深则只需一道工序，可以大大减少工作量，提高成形精度和成形极限。

  

图1 汽车前翼子板三维模型

2 模型的建立

板料选取Stl4ZF [17]深拉深冷轧钢板，厚度为1 mm，其性能参数如表1所示。(E为杨氏模量，Rm为屈服强度，Re为抗拉强度，r为各向异性指数，n为硬化指数，v为泊松比。)材料模型为36号各向异性弹塑性材料，板坯单元采用弹塑性材料。板坯单元采Belytschko-Tsay [18]薄壳单元。毛坯与模具之间的接触类型为单向面-面接触。

 

表1 材料性能参数

  

ERmRernv20717532018023028

有限元模型模拟的准确性高低在很大程度上取决于所采用的材料模型能否真实准确地反映材料变形时的屈服特征。由于三参数Balart模型能够合理地描述具有较强各向异性结构金属板材的屈服行为，故而它是专门针对金属薄板成形而建立起来的材料屈服准则模型[19]。该模型是目前用于汽车覆盖件冲压成形分析时最理想的模型。三参数Balart材料模型所采用如下形式的各向异性屈服准则：

F=a|k1+k2|m+a|k1-k2|m+[2-a)|2k2|m=2m

式中：；；；

汽车翼子板结构复杂，拉深成形时各处板料进料不一致，成形件局部出现起皱和拉裂等缺陷，拉延筋的设置可以有效地改善板料的拉延情况，提高板料的变形均匀性。根据翼子板各处进料规律，将翼子板设置为8段拉延筋，拉延筋采用半圆形，深度为5 mm，如图4所示，由于直壁区进料过快，采用3段拉延筋，而弯曲处进料情况复杂，将拉延筋进行分段设计。

为了使汽车翼子板充分成形，同时减起皱与破裂缺陷等的产生，在汽车翼子板充液成形过程中，需要将翼子板空洞补充完整，使翼子板更好地成形，提高板料的整体利用率。经过分析研究，对其主截面线采取多端段设计，同时附加一段斜壁，为了使各段补充面光滑地连接起来，还需要手动进行局部补充面设计，最后确定完整的工艺补充面如图2所示。

2.1 工艺面补充

8）脚手架处于顶层连墙点之上的自由高度不得大于6m。当作业层高出其下连墙件2步或4m以上、且其上尚无连墙件时，应采取适当的临时撑拉措施。

  

图2 工艺补充面

确定翼子板模型后，由于BSE板料估算坯料尺寸快速且准确，故而利用Dynaform里面的BSE模块进行坯料尺寸估算，零件以数据类型IGS格式导入有限元模拟软件中，1号线为有限元软件估算出的坯料尺寸，但是由于翼子板的复杂性以及需留出拉深余量，实际生产需将1号线扩大化与规则化为2号线，如图3(a)所示。模型导入仿真软件后，通过偏置和复制等工序，生成模拟所用的压边圈和充液凹模，得到其有限元模型如图3(b)所示。

  

图3 坯料尺寸及有限元模型

3 拉延筋参数的设定及优化

式中的x，y分别表示板材的纤维方向和横向；a，h，p是描述塑性各向异性的材料系数；是等效应力；对于体心立方晶体材料，指数m取6；对于面心立方晶体材料，m取8；r0，r45，r90为各向异性材料在轧制方向成 0，45°，90°时的各向异性值。三参数Balart 对 应 Dynaform的36号材料模型[20-21]。

  

图4 拉延筋设置

拉延筋所提供的阻力严重影响翼子板成形性能，根据大量的数值模拟，提出三种拉延筋研究方案，如表2所示。

按照表2的三组拉延阻力设定方案，采用Dynaform对翼子板成形性能进行模拟，分别按照3种研究方案进行加载，三组拉延方案翼子板成形极限图如图5所示。当采用方案1时，从图5(a)中可以看出翼子板局部发生严重开裂，翼子板整体拉延不充分，这是由于阻力系数比较小，拉延筋没有很好的发挥其作用，在翼子板成形时，各处板料发生严重进料流动不均匀；当采用方案3时，从图5(c)中可以看出，翼子板还没有完全成形就已经发生了局部破裂的现象，这是由于阻力系数设置过大，阻碍翼子板成形时板料进料，导致板料发生局部应力集中而破裂；当采用方案2时，从图5(b)中可以发现，板料流动性好，翼子板可以均匀成形，所以经过如上的对比分析确定方案2为拉延筋最佳加载方式。

 

表2 拉延筋参数

  

拉延筋编号方案1方案2方案3阻力百分数/%偏离距离/mm阻力百分数/%偏离距离/mm阻力百分数/%偏离距离/mm12345678510105555525252525353535401530301525151520252525253535354060906060606060602525252535353540

  

图5 翼子板成形极限图

4 加载路径的设定及优化

汽车翼子板充液拉深成形过程中，液室压力是翼子板成形性能的重要工艺参数，液室压力的加载时间与大小是液室压力的研究重点，根据翼子板的成形特点以及数值模拟研究，设计了三条液室压力加载曲线，如图6所示。

  

图6 液室压力加载曲线

加载曲线1加载方式为在凸模开始之时就施加于板材提供背压，并随着加载时间的增大而增大，液室压力加载曲线3是在0.028 s开始施加液压，在0.033 s时加载到6 MPa，保持0.002 s后增加到14 MPa， 再保持0.002 s，最后在0.04 s时达到25 MPa。而液室压力加载曲线2是在0.028 s开始施加液压，在0.033 s时加载到2 MPa，保持0.002 s后增加到8 MPa，再保持0.002 s，最后在0.04 s时达到25 MPa。

该小型撬装式LNG气化站的规格为:在竖直方向上还能够对换热器进行布置，并且不会改变现有水平方向上的空间大小。优化方案确定为:

采用加载曲线1时，液室压力加载过早，板料成形时出现大面积的凹陷而无法正常成形，故而不能满足翼子板成形需求，如图7(a)所示。采用加载曲线3时，虽然翼子板能够成形，但是其多处出现破裂现象，如图7(c)所示，因此也不能满足翼子板需求。采用液室压力加载曲线2时，翼子板可以成形而且没有破裂，如图7(b)所示。经过对比分析，确定方案2为最优加载路径。

  

图7 翼子板成形模拟结果

5 压边力的设定及优化

为了进一步分析压边力对翼子板成形质量的影响规律，在其壁厚变化较大的沟槽区选取截面A-A，顺时针选取15个点，如图9(a)所示，分别测量每个点在不同压力值下的壁厚，得到各个点在不同压边力作用下壁厚分布曲线如图9(b)所示。从壁厚分布曲线中可以看出，翼子板在弯曲处6点的减薄量最大，最容易发生破裂缺陷；在1点、14点、15点处发生了起皱，壁厚值已经超过了板料的初始值1 mm。通过对比分析三种不同压边力作用下的壁厚分布曲线，发现当压边力为2 000 kN时，各个典型点的减薄量最大，在点6处最小壁厚达到了0.66 mm；当压边力为75 kN时，各处减薄量最小，点6处的最小壁厚为0.78 mm；当压边力为1 000 kN时，点6处最小壁厚为0.74 mm。结合不同压边力的成形极限图和壁厚分布曲线，发现当压边力为1 000 kN时既可以减少起皱又可以满足翼子板的成形需求，所以将压边力定为1 000 kN。

1)根据汽车前翼子板成形的复杂性，成形时拉延筋需要进行多处单独设置，保证板料成形时的进料均匀性，根据多次模拟分析，选择的拉延筋为8段，同时设计了三组拉延阻力系数方案，以翼子板成形极限图作为评定标准，对比分析最优方案为：方案二；

当压边力为750 kN时，局部区域出现了严重起皱和拉延不充分的缺陷，如图8(a)所示，这是由于压边力过小，造成板材的流动不充分。当压边力为2 000 kN时，虽然翼子板没有出现严重起皱现象，但是其局部出现破裂，如图8(c)所示，这是由于压边力过大，所以阻碍了板料的正常成形，导致板料进料不充足而严重减薄。当压边力为1 000 kN时，翼子板成形质量较好，同时局部起皱得到了很好的控制，如图8(b)所示。

  

图8 翼子板成形极限图

由于压边力的大小会影响板料的流动，因此压边力是翼子板充液成形时的一个重要工艺参数。

  

图9 翼子板壁厚分布

6 结 论

经过大量的数值模拟的研究，设定三组翼子板压边力的加载方案，且分别为750 kN，1 000 kN与2 000 kN。

片面认知：提到畜禽养殖，人们就会把它和污染联系在一起，养殖和污染似乎已经划上了等号。而实际上种养结合的养殖场，对畜禽废弃物利用是非常充分的，它和通常意义上的环境排放污染物有本质区别。

2)针对翼子板的充液成形，进行大量的数值模拟，从而设定三组液室压力路径加载方案，并以翼子板的成形极限图(FLD)作为评定标准，对比分析得出的最优加载方案是在0.028 s开始施加液压，在0.033 s时加载2 MPa，保持0.002 s后增加到8 MPa， 再保0.002 s，最后在0.04 s时达到25 MPa；

采用SPSS 18.0软件对数据进行分析处理，计量资料以（均数±标准差）表示，采用t检验；计数资料以（n，%）表示，以P＜0.05表示差异具有统计学意义。

3)根据多次模拟研究，采用三组不同的压边力，并以翼子板的成形极限图(FLD)作为分析依据；选取翼子板沟槽区的一个截面，进一步分析压边力对其壁厚的影响情况，以截面的壁厚分布图作为分析标准，最终经过对比分析选取最优的压边力为1 000 kN。

参 考 文 献:

[1] 刘晓晶，刘博，陈晓晓，等. 汽车复杂梁形件冲压成形及回弹数值模拟[J]. 哈尔滨理工大学学报，2015，20(1)：76-79.

[2] 刘晓晶，王祥刚，李峰，等. 基于UG的汽车覆盖件模具的标准件库开发[J]. 哈尔滨理工大学学报，2011，16(2)：19-21.

[3] 吴文娟. 浅谈汽车前翼子板的结构设计[J]. 机电工程技术，2011，40(8)：196-206

[4] 张德超. 汽车车身造型与结构分析[D].湖北：武汉理工大学，2006：1-69.

[5] 武振锋，贾凡. 基于CATIA的轿车前翼子板曲面造型与质量分析[J].汽车技术，2011(1)：54-57.

[6] 肖冰娥，温彤，于建明. 大高径比H68铜套的多道次深拉深工艺与模具设计[J]. 热加工工艺，2013，42(9)：87-89.

[7] 陈文琳，路遥，彭李静，等. 多道次拉深复合成形工艺研究[J]. 精密成形工程，2011，3(1)：23-26.

[8] 陶宏之. 板料多道次拉深成形模拟技术研究与工艺优化[M]. 上海：上海交通大学，2000：45-74.

[9] 汪建敏，钱春苗，崔会杰. 数值模拟技术在板料多次拉深成形中的应用[J]. 热加工工艺，2012，5：136-139.

[10] 孙志莹，郎利辉，孔德帅. 铝合金马鞍形件充液成形工艺模拟分析[J]. 精密成形工程，2015，7(1)：46-65.

[11] XU Y C, CHENG Y, YUAN S J. Influences of Loading Paths on Thickness of Aluminum Alloy Cup with Hydromechanical Deep Drawing [J]. International Journal of Materials and Product Technology, 2010(38):173-183.

[12] GRIZELJ B, CUMIN J, ERGIC T. Application of Hydroforming Process in Sheet Metal Formation [J]. Journal of Materials Processing Technology, 2009, 51(4): 333-340.

[13] LEE M Y, SOHN S M, KANG C Y, et al. Study on the Hydroforming Process for Automobile Radiator Support Members[J]. Journal of Materials Processing Technology, 2002(130): 115-120.

[14] KEUM Y T, LEE K B. Section Finite Element Analysis of Forming Processes for Aluminum-alloy Sheet Metals[J]. International Journal of Mechanical Sciences, 2000(42): 1913-1918.

[15] GUO J D, ZHAO C C, CAO M Y. Process of Back Pressure Deep Drawing with Solid Granule Medium on Sheet Metal[J]. Journal of Central South University, 2014, 21(7): 2617-2626.

[16] 赵军，易际明，倪颋. 基于数值模拟的深筒形件多次拉深工艺与实验研究[J]. 热加工工艺， 2012(5)：94-101.

[17] 张顺才. 超深冲冷轧钢板St14ZF退火组织与力学性能关系的探讨[J]. 宝钢技术，1994(1)：41-51.

[18] VIGNJEVIC R, CAVALCANTI M J. Numerical Simulation of the Lynx Helicopter Main Lift-frame Component Collapse [J]. International Journal of Crashworthiness，1997，2(1)：25-38.

[19] 陈正阳. 一种板料成形屈服准则的建立与应用[D].北京：北方工业大学，2015.

[20] 赵海鸥. LS-DYNA 动力分析指南拉延筋技术[M].北京：兵器工业出版社，2003.

[21] 许允超. 高强度冷热轧汽车板材成形性能数值模拟及实验研究[D].内蒙古：内蒙古科技大学，2015：1-24.