【论文关键词】: 有限元; 数值模拟; 汽车覆盖件; 拉延 

　　【论文摘要】：针对汽车覆盖件冲压的有限元模拟方面的具体问题进行了研究，采用弹塑性有限元的数值模拟及试验研究的方法，对汽车覆盖件拉延过程中的成形进行了研究。针对拉延模拟结果进行应力应变分析，寻找工艺参数的优化方案，改进的工艺方案使破裂情况明显改善。 

　　由于冲压工艺具有生产效率高、质量稳定、成本低以及可加工复杂形状等一系列优点，在机械行业的应用非常广泛，占有十分重要的地位[1]。但是冲压模具的设计主要依据工程师长期积累的经验。对于复杂的成形工艺和模具，设计质量难以得到保证;一些关键性的设计参数要在模具制造出来之后，通过反复的调试、修改才能确定。这样就浪费了大量的人力、物力和时问[2-3]。随着有限元技术和计算机技术的发展，数值模拟已逐渐成为工艺分析及优化设计的有效工具。 
　　 
　　1. 有限元模型的建立和参数设定 
　　 
　　一般汽车覆盖件工艺设计流程具体分析如下: (1) 根据产品图及产品冲压工艺设计，进行详细的车身产品工艺性分析。为了实现拉延或创造良好的拉延条件，必须合理考虑冲压方向、工艺补充部分形状以及压料面形式、拉延筋布置等重要工艺因素。其中包括利用计算机进行的工艺补充面三维设计。(2) 在满足产品使用的前提下，将过剩的质量要求及时反馈给产品设计部门，进行研讨，力争把产品完善到最简单、最合理的工艺要求，以克服产品的过剩质量，减少不必要的工装投入。(3) 利用计算机进行车身产品的冲压工艺性分析，进行图面形状的分析探讨和尺寸公差的分析研究，在充分理解、把握产品使用性能要求的前提下，考虑用户使用和维修，利用塑性加工原理、冲压工艺知识和模具设计结构的有关知识，设计冲压工艺过程图。在设计过程中，同时要分析冲压工艺方案，发现不足之处，进行必要的修正。(4) 模具设计人员按照冲压工艺过程图的基本要求进行模具设计，模具cad设计包括上、下模座，工作部分零件，导向部件，定位零件和进出料装置等设计。数控编程和模型人员按照冲压工艺过程图和模具图进行数控编程和模型制造，最后按照冲压工艺过程模具图要求进行机械加工和模具装配调试，最终调试出合格的产品。 
　　选用某轿车内部地板零件产品图，此零件是一个比较复杂的中小型车身结构件。由于零件拉延深度深，并且具有局部反拉延，因此成形过程估计会出现问题，为了验证问题所在我们利用cae软件进行模拟成形计算。对于复杂冲压零件的成形过程，不但同一时刻不同位置的板坯所承受的变形方式和变形程度不同，而且不同时刻同一位置的板坯所承受的变形方式和变形程度也不同;另外，冲压工艺边界条件的设定对变形路径和各部分的变形程度的影响也非常明显。 
　　一般划分网格时，首先建立一个拓扑结构模型。这一步骤是连接分离的型面，使你可以在网格划分的时候得到连续的网格(两个相连的元素在分界线之问共同享用相同的节点)。系统能通过你所定义的公差自动辨认普通表面之问的分界线，以建立我们所说的拓扑模型。建立好拓扑结构以后，应定义网格划分的参数，并进行网格的自动划分。一般情况下要求用户最少确定四个参数，包括最小元素大小，最大元素大小，两个相连的元素之问的法向夹角，网格的弦高。最小元素的大小影响着网格划分中最小元素的尺寸。当模型的型面比较平坦时它最大元素的大小则受最大元素参数的影响。两个相连的元素之问的法向夹角所起的作用是规定了两个相连元素之问的最大法向夹角，即当两个元素的夹角大于用户给定的值时，这两个元素会分裂为更多的元素，故它影响着倒角和小圆角部分的网格密度，它的值越小网格则越密。例如:一般我们在划分模具网格时，它的拉延圆角最好有五行元素，这时调整法向夹角的参数就可以达到目的。弦高的大小则影响着大网格半径表面上的网格密度，它的值越大，则网格越少。在汽车覆盖件模拟中，板料数据一般都是曲线，因此板料的网格划分与模具的划分不一样。 

　　根据实际需要确定板料特性，应力应变关系=537(0.0102+)0.23mpa，法向各项异性系数为1.8。其他参数如下:扬氏模量2.07e+5 mpa;屈服极限210 mpa;泊松比0.28;板材厚度0.8mm;板料质量密度7.83e-9;r0=1.87, r45=1.27,r90=2.17。由于摩擦系数必须有实验得出，特别是几种常用材料在工业生产中的实际摩擦系数。考虑到汽车覆盖件生产厂家和模具生产厂家的实际，一般不考虑使用润滑油，在拉延前要使用清洁防锈油清理兼润滑。因此我们必须通过试验来得出在几种不同条件下的摩擦系数，例如干摩擦和加清洁防锈油后的摩擦。还有就是拉延筋的拉延阻力在不同形状拉延筋情况下的取值。测定为此我们设计了覆盖件模具的摩擦系数和实际拉延筋拉延阻力的测定的试验，详细试验结果在第六章中。摩擦系数根据测量结果给定0.175 ,拉延筋选取单圆筋，拉延阻力为0.178kn/mm。 
　　 
　　2. 汽车覆盖件冲压的有限元模拟结果分析 
　　 
　　经过计算后，板料的fld如图2所示。在fld图中，红色表示破裂，粉红色表示起皱，而在应变云图中红色表示正应变，深蓝表示副应变。从fld图中我们可以看出四处破裂，分别是大鼓包处，凹坑底部，最下方的小鼓包处，右上方的直壁处。通过主应变和次应变云图可以看出在突起的鼓包顶端处为双向拉应变发生破裂，并目_从板料轮廓的变化发现在有拉延筋的地方板料儿乎没有流动，形成过度胀形，凹坑底部破裂处也同样出现胀形过度问题。而模具拉延直壁处的破裂却是不同形式的，该处的主应变为拉应变，次应变为压应变，为明显的拉深破裂状态。之所以只有这个直壁角破裂是因为这个角离大鼓包最近，并且通过成形过程的模拟我们发现这个直角壁首先成形，从而在凹坑成形前破裂。其它四个角由于拉延高度低并且没有复杂的凸凹变形，都有足够的板料流动量，板料的流动情况良好，所以没有破裂。 
　　 
　　3. 汽车覆盖件冲压工艺改进方案 
　　 
　　在去掉拉延筋，变化压边力后还是无法缓解，于是决定改变模型，我们把拉延直壁消除降低了模具拉延高度;把型面中那一个接近大直角型面过渡改为一个小缓坡，减缓了陡峭程度;由于模具进料困难，所以去掉拉延筋，然后设定压边力为400kn,摩擦系数为0.12，进行模拟后如图4所示。可以看出与未改前的情况有很大的不同，破裂情况明显改善，尤其是右上角直壁处的破裂变得很小，这是由于降低了它的拉延高度。 
　　 
　　4. 结论 
　　 
　　世界上每年的钢材有半数以上被轧制成板料和管料。金属板、管的成形和加工在航空、航天、汽车、船舶及许多民用工业中都占有相当重的比例。因此，提高相应的成形技术和制造水平是一个具有普遍意义的大课题。因此，文章在汽车覆盖件数值模拟和试验研究的基础上，采用有限元的数值模拟及试验研究的方法，对汽车覆盖件拉延过程中的成形进行了数值模拟和试验研究。 
　　 
　　参考文献 
　　[1] 李东升, 黄小明, 胡世光. 汽车覆盖件拉延筋的单元模拟试验研究[d]. 北京航空航天大学学报，1995.21(2):67-71. 
　　[2] 刘建华, 姜全昌, 土志恒. 环形拉深筋成形力的计算及试验研究[j]. 塑性工程学报，2000,7(3):15-18. 
　　[3] 何丹农, 陶宏之, 包向军, 等. 拉深用润滑剂的摩擦系数和压边力的关系研究[j]. 摩擦学学报，2000,20(3):223-225.