变压边力冲压成形工艺设计

本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 I 页 共页目 录1 绪论 .11.1 引言 .11.2 方盒件变压边力冲压成形的主要研究方向 .21.2.1 冲压成形工艺 .21.2.2 变压边力 .21.2.3 板料冲压成形 CAE 分析方法 .31.3 主要研究目标及内容 .41.4 论文的组织结构 .52 方盒件冲压成形工艺理论分析 .72.1 方盒件成形的拉深成形原理 .72.2 方盒形件变形特点 .82.3 方盒形件主要拉深失效形式 .93 方盒形件成形过程的变压边力分析与研究 .123.1 临界压边力影响因素分析 .123.2 板料成形安全区域研究 .143.3 变压边力中的控制曲线 .153.4 仿真与分析步骤 .154 基于 CAD/CAE 的方盒件变压边力冲压成形仿真与建模 .174.1 DYNAFORM 软件概述.174.2 基于 CAD/CAE 的方盒件模型建立 .184.3 方盒件冲压工艺模拟与分析前处理 .194.4 方盒件冲压工艺模拟与后处理分析 .214.5 不同压边力控制曲线仿真结果及分析 .224.6 同一压边力不同控制曲线结果分析 .294.7 小结 .32结束语 .33致谢 .34参考文献 .35 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 1 页 共 36 页1 绪论1.1 引言拉深是板料冲压成形中的主要成形方式,在汽车、航空航天、石油化工等诸多等领域均有广泛的应用。例如,汽车车身覆盖件、大油罐等的成形均采用拉深工艺。然而,如何保证板料塑性成形零件的质量,降低废品率,减少模具的返工,缩短模具设计周期,一直是板料塑性成形领域的一大难点。板料成形过程包括了非常复杂的物理现象,涉及力学中的三大非线性问题:几何非线性、物理非线形和边界非线性 1。因此,难以用传统的弹塑性理论的解析法进行研究。利用传统的方法进行冲压件工艺分析时，为了避免材料成形过程中出现断裂、起皱、颈缩等不良影响，必须反复修改成形加工的某些参数或修改模具形状，使得工艺过程耗资大、产品开发周期长，已不能适应激烈的场竞争和现代工业的发展要求。随着计算机技术的迅速发展及有限元方法的成熟,特别是商用有限元软件的不断完善,促进了板料成形中数值模拟技术的发展。该技术减少了试模过程,缩短了产品开发周期,降低了成本，对提高制成品的制成质量提供了良好的保障。方盒形件广泛应用于电子部件和汽车部件中，如电池盒、半导体盒和汽车反光镜、汽车油箱等。在实际生产中拉深件多为非轴对称形状，而方盒形件是典型的非轴对称形的拉深件，在方盒形件拉深过程中材料的变形较为复杂，因此以方盒形件为研究对象对拉深工艺作进一步的探讨更具有实际意义 2。传统的压力机在冲压过程中压边力一般为恒定值，这制约了材料的成形性能，也影响了冲压件的尺寸精度和稳定性。为了获得更好的成形质量，变压边力控制曲线对成形性能的影响和变压边力预测成为目前冲压技术研究的关键技术之一。因此，选用非轴对称件方盒件，借助 CAE 对不同类型压边力控制曲线的板料成型性能的仿真，寻找合适的压边力控制曲线，将弹塑性力学理论和人工智能技术结合起来，提出预测板料拉伸过程中压边力控制曲线的办法，大幅度减少设计和试验的工作量，降低成本费用，提高板料的成形性能，最终奠定对复杂形状冲压件压边力的智能预测研究的基础，对智能数控冲压机床的研究和设计都具有重要的现实意义。本课题针对方盒件冲压成形对压边力控制的要求，采用变压边力方式，首先选取了简单的非轴对称件方盒形件作为研究起点，借助 CAD/CAE 手段通过对不同类型压边力控制曲线的板料成形性能仿真，寻找合适的压边力控制曲线，探讨方盒 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 2 页 共 36 页件冲压成形工艺。1.2 方盒件变压边力冲压成形的主要研究方向盒形件是金属薄板拉深成形中较为典型的冲压件，其变形规律具有一定的典型性，研究这类件的成形规律，不仅对这类件成形工艺参数和工艺步骤的确定是至关重要的，同时也是进一步认识复杂件成形规律的基础。因此国内外对此展开了大量的研究。主要集中在冲压成形工艺、变压边力和板料冲压成形 CAE 分析方法等方向上。1.2.1 冲压成形工艺板料成形加工是金属加工的一种重要工艺方法，它不仅生产效率高、原材料消耗少，而且可以有效地改善金属材料的力学性能和组织，因而在国民经济中得到了极为广泛的应用，特别是在航空、宇航、汽车、造船、电器、五金等工业部门。板料冲压成形的方式有很多，通常可用四种基本变形方式来认识，即膨胀、拉延、伸长类翻边以及弯曲变形板料冲压成形的方式有很多，通常可用四种基本变形方式来认识，即膨胀、拉延、伸长类翻边以及弯曲变形 3。一般来说，板料成形性能依赖于压力、拉力、拉伸速率、温度这些与金属抗伸长断裂有关的因素，金属材料的尺寸、形状、第二相微粒的分布状况等对其性能影响也很大 4。冲压是一种高效率的生产方式，适用于大批量的生产，其材料的成本占生产成本的 60%80%，材料利用率的高低主要取决于排样方案的优劣。根据不同的加工工艺与加工方式，有普通单排、普通双排、对头单排、双头双排 4 种排样方式可供选择。在金属材料成形中，材料的成形过程受到材料的参数，工艺参数，材料尺寸等诸多因素的影响，其中材料参数主要包括材料的抗疲劳强度，屈顺强度，应变硬化指数，各向异性指数等。工艺方面参数主要是指冲压速度，压边力的大小和分布，成形润滑情况等等 5,6。1.2.2 变压边力变压边力技术研究的关键在于：如何在保证成形质量前提下改善材料的成形极限。变压边力可以分为随时间变化的变压边力、随位置变化的变压边力、随时间和位置综合变化的变压边力。(1)随时间变化的变压边力对成形性能的影响：随时间变化的变化的变压边力是指在薄板成形过程中，在压边圈采用整体压边 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 3 页 共 36 页圈的前提下，压边力仅随凸模行程变化的变压边力，采用这种压边力可以有效改善材料的成形性能。意大利的巴勒莫大学的L.Fratini和F.Micari进行了大量的模拟和实验，对不同材料和拉深比的轴对称件在不同压边力变化方式下进行成形性能研究，得出选择恰当的渐增的压边力变化方式可以使拉深深度最大的结论。新加坡的S.Thiruvarudchelvan以铝板、铜板杯形件拉深为例，分析了压边力与冲压力成正比（约34）拉深过程，在实验的基础上得出采用变压边力形式将有利于提高圆筒形件的成形质量，降低了拉延开始时的拉裂危险性，消除起皱时所用的最大压边力比采用恒压边力所需最大压边力要小，并且杯壁部分的减薄量小于恒定压边力作用时的减薄量。Ahmetoglu等人对渐减式压边力曲线进行了进一步的实验研究，实验和有限元分析的结果表明，从较大的压边力逐步减少到仅能消除起皱的恒定压边力的变化曲线能显著减小法兰起皱的高度，并且可以避免较大的恒定压边力所造成的断裂。(2)随位置变化的变压边力对成形性能的影响：随位置变化的变压边力是指在薄板成形过程中,在冲压件的法兰区不同位置施加大小不同但恒定的压边力，采用这种变压边力控制技术,不仅可以提高薄板的成形性能,减少和消除成形过程中起皱和断裂等缺陷,而且能提高冲压件的尺寸精度和冲压过程的稳定性。(3)随时间和位置综合变化的变压边力对成形性能的影响：随时间和位置综合变化的变压边力是指在薄板成形过程中,在冲压件的不同位置施加大小不同的压边力。这种变压边力技术是综合随时间变化的变压边力技术和随位置变化的变压边力技术，其更能提升材料成形性能。特别是对非对称件，这种变压边力更能改善材料的成形极限和质量。目前,国内外对变压边力改善板料成形性能的内在机理仍不够十分了解，研究结果不能直接用于指导实际生产，同时，最佳变压边力的确定仍依赖于手工“试错” ，并且缺乏有效的确定最佳变压边力优化设计方法。研究发现，通过数值模拟与正交优化的有效结合，对方盒件成形进行研究也表明了正交设计优化方法对于方盒件冲压成形参数的优化是一种有效可行的途径。在冲压成形过程中，压边力的大小是影响冲压件起皱和破裂是否出现的主要因素 7,8。适当的控制压边力在拉深加工中的数值，能够更好的控制制品成形的质量。 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 4 页 共 36 页1.2.3 板料冲压成形 CAE 分析方法板料成形数值模拟技术涉及到塑性力学、有限变形理论、有限元理论、板壳理论、计算数学、计算机辅助设计、计算机图形学和可视化技术等多方面的理论和技术。一般来说，一个成功的金属板料冲压成形 CAE 系统从功能上讲都是由前置处理、有限元求解器和后置处理三部分组成 9。板料成形 CAE 分析的一般过程为前处理(网格划分、定义约束条件、施加条件和边界条件等)，求解(数值计算) 和后处理分析结果(应力、应变、材料厚度分布、能量等历史曲线、云图及动画、切取截面显示、回弹、FLD 图等)。金属板料冲压成形 CAE 技术主要是通过对结果的分析，使产品设计、模具结构更合理、工艺参数更满足要求。板料成形 CAE 分析的一般步骤如图 1.1 所示。DYNAFORM 是目前最常使用的 CAE软件。DYNAFORM 的主要功能包括分析拉伸、成形、回弹、翻边、切边等板料成形过程中的不同工序，也可以进行多步成形（或多工序加工）分析 10。通过用户已定义好的冲压工艺及模具曲面形状来预测成形状态，其中包括减薄拉裂、起皱、回弹等各种问题，同时可以对成形力、压边力、拉延筋、模具磨损等各种工艺问题进行分析，以便优化工艺和模具设计。对于容易出现局部开裂的现象，拟定利用 DYNAFORM在计算机上基于实际的边界条件进行数值模拟计算，验证和优化工艺参数。图1.1 板料成形CAE分析的一般步骤1.3 主要研究目标及内容本项课题以非轴对称件方盒形件为研究对象。通过对简单的方盒件进行有限元数值模拟。了解板料在拉深成形过程中冲压成形状况。通过比较几种不同情况下冲压成形效果较好的一种。并进一步通过变压边力法来优化其成形工艺。 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 5 页 共 36 页本课题的主要内容归纳如下：（1）分析国内外的研究现状，确定本课题的研究内容；（2）分析方盒形件拉深过程中所产生的缺陷，确定影响成形质量的因素，并初步提解决方案；（3）介绍板料成形的有限元数值模拟基本理论及 Dynaform 软件的数值模拟过程，确定研究方盒件成形工艺的具体仿真实验方案；（4）选择板料及模具的参数，使用 Pro/E 建立 CAD 模型，Dynaform 建立冲压模具、压边圈、板料的有限元模型，建立方盒件模型，选定合适的试验参数 11,12；（5）采用不同的压边力对方盒件来冲压成形性能进行比较，得出较优解。本论文将对仿真模型进行以下分析：（1）拉深过程施加不同的压边力，找出压边力对成形质量的影响特点，从而确定合适的压边力大小；（2）修改变压边力中控制曲线，分析不同压边力下同一加载方式的不同冲压成形效果；（3）在已寻找到的合适压边力下，分析同一压边力下 6 种不同加载方式的不同冲压成形效果，最后确定较好的加载方式。1.4 论文的组织结构本文以板料拉深成形中比较有代表意义的方盒形件为主要研究对象，基于CAD/CAE 的方盒件变压边力冲压成形工艺研究使用板料成形模拟仿真专用软件包DYNAFORM 对方盒形件进行仿真和分析，获取合适的压边力控制曲线，并寻找方盒形件成形的合适的压边力和合适的压边力下最佳的加载方式。全文分为五章，各章内容如下：第 1 章：简要介绍了课题研究的背景和意义，分析了目前方盒形件变压边力冲压成形相关技术的主要研究方向，确立了课题的研究内容和目标，阐述了论文的结构安排。第 2 章：分析方盒形件拉深过程的变形特点，并简要分析了板料拉深成形失效形式。第 3 章：介绍了方盒形件成形过程的压边力分析与研究的理论和方法，对方盒形件成形时破裂和起皱压边力理论进行了研究，分析了影响板料临界压边力的因素，并讨论了板料成形的安全区域。在板料成形安全区域，使用变压边力技术，且在变 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 6 页 共 36 页压边力过程中加载不同的控制曲线，以控制施加在压边圈上不同部位压边力的大小并简要列出 6 种不同变压边力的加载方式。最后给出了方盒形件成形过程的具体研究方法及步骤。第 4 章：基于 CAD/CAE 的方盒件变压边力冲压成形进行工艺研究，使用模拟仿真专用软件 DYNAFORM 对 6 种不同形式的压边力曲线控制下的方盒形件成形过程进行了仿真分析，确定了合适的压边力并在较优的成形效果下寻找控制曲线中较好的加载方式。最后简要总结了本课题研究所取得的成果与不足，并对今后的工作进行展望。 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 7 页 共 36 页2 方盒件冲压成形工艺理论分析板料拉深是在模具的作用下，板平面内产生切向压应力和径向拉应力，坯料通过拉深凹模向直壁流动，使板料成形为空心零件，或浅的空心毛坯成形为更深的空心零件的冲压工序 13,14。拉深变形总是与弯曲的胀形等其他方式的变形同时发生。因此需要从方盒形件的拉深成形工艺着手，了解其主要成形缺陷机理以及板料成形失效中破裂和起皱的评价标准从而寻找到合适的压边力变化的规律和范围 15。2.1 方盒件成形的拉深成形原理方盒形件拉深成形是属于盒形件一种基本的拉深成形方式。如图 2.1 所示，宽度为 L 的板料在冲压力、压边力和圆角半径为 的凸模、圆角半径为 、转角半径prdr为 、截面宽度为 B 的凹模共同作用下，被拉入凸模和凹模之间的间隙里，成为开cr口的方盒形件。方盒形件上高度为 h 的直壁部分是由毛坯的法兰区部分转化而成的。图 2.1 方盒形件拉深原理方盒形件可以看作由直边部分及圆角部分组成。方盒形件拉深变形时，圆角部分近似圆筒形件的拉深，直边部分近似板料弯曲。因此方盒形件的拉深成形是圆角部分拉深、直边部分弯曲两种方式的组合 16。但是，方盒形件的圆角及直边是联系在一起的整体，因而变形时必然又相互的作用及影响，以至圆角不是简单的圆筒件拉深，直边部分不是单纯的平板弯曲。在直边部分横向间距缩小，愈靠角部缩小愈多；纵向间距增大，愈靠盒形转角处增大愈大。这是由于在拉深时，圆角部分的材料要向直边部分流动，使直边部分材料受压。其结果使横向间距变小，纵向间距变大，愈近角部变化愈大。圆角部分并不是单纯的弯曲变形。由于圆角部分的材料要向直边部分流动，因而直边部分也增加了横向压缩、纵向伸长的变形。而圆角部分， 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 8 页 共 36 页由于直边区存在金属的流动，使得圆角部分的变形程度大为减小 17,18。2.2 方盒形件变形特点因此，盒形件得变形特点可归纳为以下几点：（1）盒形件拉深时，角部变形基本上与圆筒形件拉深变形相似，只是由于金属向直边流动，使得径向应力 在角部的分布是不均匀的，拉应力 在圆角中间处最rr大，直边中间处最小，而且平均拉应力比相应圆筒形件的拉应力小的多，于是减小了危险断面的载荷，不易被拉断，故盒形件的成形极限较高。（2）拉深时，压应力 的分布也如同 分布不均匀，在圆角中间处最大，直r边中间处最小，如果直边较长，中间部分 近似为零。（3）盒形件的最大拉应力 和最大压应力 均出现在角部。因而破裂、起皱r等现象也多在角部发生。在远离角部的直边部分一般不会产生起皱 19,20。（4）由于变形速度的不一样，从而导致应力分布不均匀。因此在方盒形件的圆角区和直边区的过渡区域，存在剪切应力 。在过渡区域的中间处剪切应力 最大， 然后向两边逐渐减小。由此可见，毛坯的变形主要集中在凹模表面的凸缘上，拉深的过程就是使凸缘部分逐渐收缩，转化为方盒形件壁的过程。拉深过程中凸模圆角以上的材料为主要变形区，而底部通常是不参加变形的不变形区。盒形件拉深时，不仅在直边部分的凹模圆角处产生弯曲变形，而且在法兰区产生切向收缩、径向伸长的变形与圆角区的变形性质基本相同。通过分析盒形件成形机理指出，盒形件变形时，变形区圆角处与直边处的金属由于位移不同而诱发的剪应力能降低圆角处传力区的轴向拉应力；由于变形区直边处产生拉深变形，使圆角处变形区的变形抗力得到降低。方盒形件是一种非轴对称的零件，它的筒壁是由直边部分与圆角部分构成。拉深变形时，应力、应变在变形区内沿周边的分布很不均匀，而且随着零件的几何参数、毛坯形状以及拉深条件的变化，其不均匀程度也会变化 21,22。此外还有板材本身的各向异性等，这些构成了盒形件拉深的一系列特点，也是盒形件拉深的难点所在。因此，人们对于盒形件的成形机理，缺陷的分析，成形极限和材料各向异向性对拉深的影响等作了大量的研究。由于盒形件形状的非轴对称性和自身结构的复杂 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 9 页 共 36 页性，决定了盒形件成形时不同于圆筒形件最大特点就是不均匀性，这种不均匀性主要体现在以下几个方面：（1）受力分布不均匀性盒形件拉深变形过程中，在法兰变形区，圆周部分的正应力沿切向和径向的分布、直边部分正应力沿横向和纵向的分布都是不均匀的。在已变形的直壁部分，直边区和圆周区所承受的单位载荷也是不同的。（2）变形分布不均匀性盒形件成形过程中变形主要在法兰区进行，而冲头的底部材料基本不变形或有少量的胀形部分。研究表明:法兰区平面内圆周区和直边区的变形程度是不同的，圆周区的变形程度大于直边区。（3）变形速度不均匀性盒形件拉深成形时，直边部分的流动速度大于圆周部分的流动速度。但是，直边区和圆周区是一连续的整体，因而其位移速度场的分布也是连续的 23，不存在间断点，两区之间也不存在明确的分界线。2.3 方盒形件主要拉深失效形式由于盒形件的成形情况较为复杂，既包括圆角部分的圆筒形拉深,又存在直边上的弯曲和拉延,加之盒形件的应力应变状态也很难分析,所以它的成形工艺设计和控制变形非常困难，其成形过程常常产生许多缺陷，如起皱、拉裂、回弹、内凹、弹复补偿不当、滑痕、形状歪扭等。方盒形件在拉深时的主要失效形式有破裂、起皱、回弹、内凹等。（1）起皱所谓起皱是指在原为光滑的平面或曲面上出现褶皱的现象。在板料冲压成形时，为了使板料产生塑性变形，模具要对板料施加一定的外力，在板料内部产生复杂的应力状态。由于板厚方向尺寸与其他两个方向的尺寸相比很小很小，因此变形过程中厚度方向是不稳定的。故当外力在板平面内引起的压应力使板厚方向达到失稳极限时便产生失稳起皱。起皱是板料拉深成形中的一种有害现象，轻微的起皱会影响产品的形状精度、尺寸精度以及表面质量等，而严重的起皱可能会妨碍和阻止加工过程的顺利进行。因此，对起皱问题进行研究，深入了解其产生机理，科学地掌握其发生的规律，对板料成形技术的进步具有十分重要的意义。在板料冲压过程中产生的起皱是多种多样的，大致可以分为压缩力、剪切力、 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 10 页 共 36 页不均匀拉深力以及平面内弯曲力等四类。方盒件成形过程中，一般会在法兰区、壁部、过度区以及凸模头部部分产生皱褶。法兰区的皱褶也是其它部位产生皱褶的原因。法兰区的皱褶是由于料在成形过程中受到切向压应力和径向拉应力的作用，板料开始顺着凹模口的法线方向凹模腔延伸，法兰随着收缩变小，应力发生变化，多余的板料则相互挤压，切应力过大，引起法兰失稳，产生皱褶。盒形件圆周处垂直壁上产生褶皱的原因是在于法兰上先发生皱褶引起的，或是拉深凹模圆角 的磨损。dr在凸模底面上出现的皱褶是因板料在型腔内被凸模头部与凹模底部所压形成再拉深状态。冲程结束后，在凸模面及凸模肩的拐角处产生皱褶，是因为拐角处的材料在拉深过程中变薄，多余的材料相互挤压产生应力所致。图 2.2 板料起皱类型（2）破裂破裂是拉深失稳在薄板冲压成形中的主要表现之一。在以拉为主的变形方式中，板料往往过度变薄甚至拉断，这种现象实质上和起皱一样，也是变形不能稳定进行的结果。不同的是，拉伸失稳只可能发生在材料的塑性变形阶段。（a）底裂也称为拉深破裂，多出现在冲头底部的圆角区，严重时直边部分也会出现，如图 2.3 a 所示，这种破裂首先发生在底角圆弧与直壁相接的过渡区，其原因是这个部位在变形过程中发生两次弯曲，变形程度较大，材料变薄严重，当变形力超过其承载能力时，就产生破裂，多在最大变形力出现之前。（b）壁裂发生在筒壁的圆周和直边区。一般有 V 型壁裂和横向壁裂两种情况。 V 型壁裂，是一种比较常见的破裂。多发生在直壁圆角靠近凹模肩部的地方(图 2.3 c)，一般呈“V”型，无论变形力是否达到最大力，都可能产生这种破裂。关于这种破裂机理的基本观点是:法兰变形区金属流动不均匀 8，圆周部分变形程度大，硬化严重；圆周 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 11 页 共 36 页部分材料增厚，压边力集中于此，形成直壁承载时的应力集中；凹模圆角处金属材料经历了强烈的纵向弯曲、反弯曲和横向弯曲，金属严重变薄，致使破裂。图 2.3 常见的几中不同破裂情况（c）横向壁裂发生在盒形件直壁区直边上的横向裂纹。多发生在尺寸比较大的盒形件，特别是有拉深筋或在凹模肩部受到胀形及反复弯曲的直壁部位更容易发生。一般观点认为横向壁裂的原因主要有：拉深深度过大、凹模圆角半径过小、压边力过大及拉深筋的存在、润滑条件不好等。（3）回弹在板料成形结束阶段，随着变形力的释放或消失，成形过程中存储的弹性变形能要释放出来，引发内应力重组，进而导致零件整体形状改变，这种现象称为回弹。产生回弹的原因主要由两个：一是因为当板料内外缘表面纤维进入塑性状态，而板料中心仍处在弹性状态，这时当凸模上升去除外载后板料产生弹性回跳；二是金属塑性变形时总是伴有弹性变形，所以板料弯曲时，即使内外层纤维全部进入塑性状态，当凸模上升去处外力后，弹性变形消失，也会出现弹性回弹。回弹是板料成形中不可避免的现象，它的存在影响零件成形精度，增加试、修模以及成形后校形的工作量，故在生产实际中迫切需要采取行之有效的措施。零件的最后回弹是整个成形历史的累积效应，而板料成形过程中与模具几何形状、材料特性、摩擦接触等诸多因素密切相关，所以板料成形的回弹问题非常复杂。 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 12 页 共 36 页3 方盒形件成形过程的变压边力分析与研究在冲压成形过程中，压边力（简称BHF）主要用于产生摩擦阻力，以增加板料中的拉应力，控制材料流动，使拉深过程中不起皱或不破裂。若压边力过小，会产生起皱现象；压边力过大，则会使板料产生破裂现象。故压边力的大小是影响冲压件起皱和破裂是否出现的主要因素。变压边力技术研究的关键在于：如何在保证成形质量前提下改善材料的成形极限。由于工作条件的差异、仿真软件的差异、仿真模型建立的差异、起皱区域的划分、摩擦条件的不确定性等因素，导致合适的压边力曲线还没有明确的结论。因而，进一步从理论上解析板料拉深成形过程，并结合实际生产研究出最有利于提高板料成形性能的压边力曲线是很有必要的。3.1 临界压边力影响因素分析从理论上讲，板料起皱和断裂的临界压边力并不是唯一的。下面就在前面方盒形件破裂和起皱的临界压边力理论推导的基础上具体分析模具几何形状、板料厚度 、t毛坯宽度L、摩擦系数 、厚向异性指数 、强度系数K、应变强化指数n、材料强度r、材料屈服极限 等因素对方盒形件成形质量的影响 24。bs（1）模具的几何形状凸模圆角半径 和凹模圆角半径 主要影响破裂临界压边力曲线，对起皱临界prdr压边力曲线影响比较小。不难想象，较大的凸模和凹模圆角有利于增大成形安全区域，凹模圆角半径对破裂临界压边力的影响大于凸模圆角对破裂临界压边力的影响。其中凸模圆角半径过小，成形阻力大，凸模圆角半径过大，容易造成板料冲压时，底部材料的破裂。同时合理的转角半径 ，对于提高成形质量也是有很大影响，如cr果转角半径过小，板料很容易发生破裂。同时凹模截面尺寸B，对于拉深成形也有很大的影响，凹模截面尺寸越小，材料越不容易进入模腔，材料容易发生起皱。（2）板料厚度 t板料的厚度越大，起皱临界压边力越大。当板料的厚度增加时，材料越不容易产生破裂失效，故其破裂临界压边力也增加。但是厚度较大会引起冲压过程中板料的内凹，经过仿真发现在 =1.8mm时，材料易发生内凹。t 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 13 页 共 36 页图3.1 板料厚度对成形影响（3）毛坯外形尺寸L合理的毛坯外形不仅可以减小加工用料，而且能减少甚至免除后续的修边工序，大大降低了生产成本。另一方面，合理的毛坯外形能改善成形过程中材料变形，减轻局部变薄，因而有利于提高产品质量。（4）摩擦系数 摩擦系数的大小对拉深成形过程的影响具有两面性。从减少板料在拉深过程中的摩擦损耗，减少方盒形件壁传力区的负担来讲，凹模与压边圈和板料接触的工作表面应比较光滑，尽量减小摩擦系数，有利于改善盒形件金属流动条件，有利于材料的拉深成形。另一方面，凸模与方盒形件壁之间的摩擦可以增大拉深能力。根据前面的结论可知，随着摩擦系数 的增大，破裂临界压边力减小。当摩擦系数 增 大时，则有板料法兰变形区的径向拉应力增大，相应地切向拉应力减小。因而 增大时，拉深中的起皱临界压边力则减小，有利于降低起皱的可能性。（5）厚向异性指数 r由于板料在生产过程中要经历轧制与退火等工艺，使得板料塑性在各个方向上存在一定的差异，即板料的塑性各向异性。厚向异性指数反映了板料在板平面方向和厚度方向上由于各向异性而引起应变能力不一致的情况。反映了板料在板平面内承受拉力或压力时抵抗变薄或变厚的能力。 值越大，表示板料在宽度方向上易于产生变形，可以减少起皱的可能性；在拉r应力作用下， 值越大，板料在受拉处不易变薄，可以减少被拉裂的可能性。当r值增大时，破裂临界压边力也增大，而起皱临界压边力几乎不受r值的影响，因此安全区域随着厚向异性指数r值的增大而增大。（6）强度系数K强度系数K表示材料的强度性能，其数值取决于材料的种类。当K值增大时，起皱临界压边力增大，由前面的分析可知，破裂失效时的极限载荷随着K值增加而增加， 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 14 页 共 36 页故破裂临界压边力也增加。（7）应变强化指数n应变强化指数n值也称为加工硬化指数或硬化系数，它表示板料在冷变形过程中材料的变形抗力(强度)随变形程度增大而增加的性质。n值越大的材料，在相同的变形程度下，应力增加得越多。故加工硬化越强的材料在拉深过程中越不容易产生细颈，因而危险断面的变薄甚至拉断现象也就可以相应的延迟很多。因而可想而知，破裂临界压边力随n值的增大而增大，而随着n值的增大，起皱临界压边力减小 25。所以安全区域随着应变强化指数n值的增大而增大。（8）材料的强度 b由前面的分析可知，当材料强度 增大时，破裂临界压边力增大，而其对起皱b临界压边力影响不大。（9）屈服极限 s屈服极限 越小，材料越容易发生塑性变形，材料在拉深成形时所需要的拉深s力就越小。（10）成形速度 和温度vT成形速度 、温度 对板料成形性能的影响因成形方式及材料而异。3.2 板料成形安全区域研究传统压力机压边力一般都是恒定值，不随凸模行程而变化，从而制约了板料的成形性能，同时也影响了冲压件的质量。实际上在拉深过程中，毛坯发生断裂和起皱失效的趋势在不同的阶段是不同的，合理的压边力的取值范围应该是由凸模行程内的起皱极限曲线和断裂极限曲线所围成的区域所决定的 26。该区域即是板料成形的“安全区域”如图 3.2 所示。安全区域的形状是受模具的几何参数、板料的力学性能以及工艺参数等多种因素影响的，合理的压边力变化曲线必须考虑这些因素。图 3.2 板料成形安全区域 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 15 页 共 36 页3.3 变压边力中的控制曲线在安全区域确定的情况下，压边力控制曲线的形式对成形极限的影响可以通过施加恒定压边力控制曲线、渐减压边力控制曲线、渐增压边力控制曲线等多种压边力控制曲线得以体现。凸模行程/mm压边力KN1234651-渐增 2-渐减 3-先增后增减型 4-先减后增型 5-恒定型 6-先增恒定后减型图 3.3 压边力加载模式 虽然在板料成形安全区域中可以采用多种压边力控制曲线，但是哪一种压边力曲线更能提高成性极限和成形质量，这是需要通过具体仿真分析的方法确定合适的压边力控制曲线。3.4 仿真与分析步骤为了板料制备简便也为了便于比较，本课题采用比较具有典型性的方盒件作为研究对象，从板料拉深成形过程的力学变形特点出发，并用板料成形模拟仿真专用软件包 DYNAFORM 对方盒形件进行仿真和分析，获取合适的压边力控制曲线，并对方盒形件拉深成形最佳压边力控制曲线进行实验，以寻求合适的的加载方式。（1）首先，使用 Pro/E 建立方盒件模型，对生成的实体模型单击文件，保存副本，文件类型选择*igs 格式。（2）导入 DYNAFORM，对模型进行单元网格化处理。定义板料的材料与属性，模具间距。（3）定义凸模的运动和压边力，对模型进行有限元计算，后处理。（4）分析成形极限图（FLD） ，并截下 FLD 图，回到前处理，重新定义压边力载荷曲线作用力取值范围为 100KN300KN。找出冲压成型效果合适的的压边力。（5）在合适的的压边力下，按照 6 种不同的加载方式（渐增，渐减，先增后减， 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 16 页 共 36 页先减后增，恒定，先增恒定后减）修改加载控制曲线，然后仿真分析、后处理。比较应力图、应变图、薄厚图，找出最佳的加载方式。 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 17 页 共 36 页4 基于 CAD/CAE 的方盒件变压边力冲压成形仿真与建模方盒形件是薄板金属冲压中较难成形的一类零件，并且在其成形过程中的变形特点具有一定的典型意义，因此很有必要对其进行数值模拟。在方盒形件拉深成形过程中，板材不同部位的受力状态、变形方式以及变形性质存在较大差异 27，材料的性能参数、模具几何参数和压边力等因素，都影响着方盒形件的成形规律和拉深性能。本章应用 DYNAFORM 有限元软件，方盒件变压边力采用动力显式算法模拟了方盒形件拉深成形过程，基于 CAD/CAE 进行冲压成形工艺研究。4.1 DYNAFORM 软件概述计算机技术的飞速发展和有限元技术的日趋成熟，采用 CAE 仿真分析金属在塑性成形过程中的变形规律在实际生产中得到了越来越广泛的应用。CAE 技术的成功应用不仅大大缩短了模具和新产品的开发周期，降低了成本，提高了企业的市场竞争力，而且还有利于有限元分析法和传统的实验方法相结合，从而推动现代制造业的快速发展。本文采用了DYNAFORM来对方盒形件拉深过程仿真研究，DYNAFORM作为一款专用的板料成形软件，能够模拟压边、拉深、翻边、弯曲、多工步成形等工艺，可以预测成形过程中板料的破裂、起皱、减薄、划痕、回弹，评估板料的成形性能。通过对这几个比较流行的软件进行对比，发现DYNAFORM无论在易操作性、前处理、后处理、模拟的准确度都具有很大的优势。与其他的有限元分析软件相比，DYNAFROM 软件有如下几个方面的技术特点：（1）DYNAFORM具有功能丰富的前处理器。它具有强大的图形文件导入功能，能够方便而无数据丢失地读入IGES格式文件以及UG、Pro-E、Catia等主流CAD软件的图形文件，同时用户也可以在DYNAFORM中很方便地创建点、线、面等几何模型。（2）DYNAFORM的求解器采用了业界非常著名的非线性动力显式有限元软件Ls-dyna。Ls-dyna是采用显隐结合的算法进行板料成形模拟的最具有代表性的软件，计算稳定，效率高，模拟结果准确性很好。虽然Ls-dyna也能够进行仿真，但是其材料库相对DYNAFORM而言就有一些欠缺，软件的易操作性也较差。（3）DYNAFORM 具有强大的后处理功能。eta-Post 是 ETA 公司开发的一款专门针对 DYNAFORM 的后处理软件，它可以方便用户直观地得到求解结果。在 eta-Post 中新增加的 GRAPH 模块中，用户可以利用曲线图表功能来显示拉深过程中各种参数随 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 18 页 共 36 页时间变化的曲线。（4）支持从个人机、工作站到巨型机的所有硬件平台。可兼容个人机、工作站、大型机及巨型机等硬件平台上的全部数据文件。在个人机、工作站、大型机及巨型机等硬件平台上具有统一的用户界面。可与大多数的 CAD 软件集成并有接口。具有智能网格划分。良好的用户开发环境。4.2 基于 CAD/CAE 的方盒件模型建立方盒形件成形应考虑到冲压件受工件几何参数、模具几何参数、材料性能、摩擦系数、材料厚度、冲压力、拉深极限、变压边力等因素的影响，因此要获得良好的实验效果，首先应建立合理工艺模型，步骤如图 4.1：用 P r o / E 建 立 方 盒 件 板 料 及 模 具 的 几 何 模 型以 I G E S 格 式 将 零 件 模 型 导 入 D Y N A F O R M零 件 层 编 辑 改 名模 具 有 限 元 网 格 划 分模 具 网 格 检 查 及 修 补定 义 成 形 工 具 （ 凸 模 、 凹 模 、 压 边 圈 ）定 义 材 料 属 性 ， 选 择 材 料 模 型 几 单 元 公 式划 分 毛 坯 网 格设 置 成 形 参 数 ( 如 压 边 力 、 冲 压 速 度 、 修 改 加 载 曲线 形 状 等 )建 立 方 盒 件 模 型求 解 器 仿 真 计 算应 力 应 变 、 厚 度 分 析成 形 极 限 图 （ F L D ）材 料 流 动设 计 结 果 是 否 满 意后处理前处理不 满 意满 意开 始结 束修改成形参数图 4.1 方盒件建模流程图 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 19 页 共 36 页（1）确定方盒形件的尺寸在 PRO/E 中建立凸模、凹模、板料的模型，在草绘状态下绘制凸模、凹模、板料的尺寸。对于凸模和凹模草绘后进行拉深，拉深深度为 20mm；对板料草绘后填充。（2）根据工件尺寸利用 Pro-e 建立方盒形件的*.igs 格式文件对生成的实体模型单击文件，保存副本，文件类型选择*igs，单击确定，在输出 IGES 窗口，选取曲面，坐标缺省，单击确定。这样就能把相应的模型转换为*igs文件，供 Dynaform 调用。同时将相应的实体模型保存以便在后面改变模具参数重新建模所用。4.3 方盒件冲压工艺模拟与分析前处理（1）将*.igs 数据导入 DYNAFORM 并对模型进行单元网格化处理。打开 Dynaform 软件，单击 file，选择 input 项，找到刚才保存的*igs 文件，依次将凸模、板料、凹模导入，凹模导入两次，因为要将其中一个凹模文件做成压边圈。单击 Parts，选择 Edit，在 Edit part 窗口为导入文件重新命名。单击 Preprocess（前处理） ，选择 surface 对导入的模型进行处理，删除不必要的表面，在 Surface 窗口选择 ，在 selectByCursor 选择 ，然后在窗口中选中不需要的表面，单击 OK 即可。单元网格化处理，单击 Preprocess（前处理） ，选择 Element，选中 ，相应的单元的参数设置如下图，Max.Size 为 30.000，Min.Size 为 0.5，Chordal Dev为 0.150，Angle 为 20.000,Gap tol 为 2.500，Ignore Hole Size 为 0.000。然后点击 Select Surfaces，点击 ，分别对凸模、凹模、压边圈进行网格化。对如板料的网格化 Tools，选择 Blank Generator，单击 SURFACE，选中板料，在 Mesh Size 窗口，Tool Radius 为 1.750000，单击确定，板料单元网格化完成。 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 20 页 共 36 页图 4.2 网格化后的模型（2）定义板料的材料、属性与模具间距板料的材料和属性设定，单击 Tools，选择 Define Blank，单击 Add 添加板料B 为毛胚，然后点击 Material，单击 Material Library，进入材料库窗口选择低碳钢 Mid 中 DQSK，Type37 所对应的材料。相应的中文参数如表 4.1：表 4.1 材料性能参数厚度t/ mm宽度L/ mm屈服极限 /MPas强度系数K/ MPa厚向异性指数 r应变强化指数n1 200 154.30 512.2 1.65 0.23单击 Property，在 Property 窗口单击 New 按钮，新建板料属性，UNIFORM THICKNESS（板厚为 1mm）为 1.000000E+000。模具间隙的确定，首先定义模具，单击 Tools，Define Tools，选中 User Define Tools，然后在 User Define Tools Neme 中打击 New 按钮新建用户工具，并在 Include Parts List 中选择相应的模型件。其中 t 对应凸模，b 对应板料，y 对应压边圈，a 对应凹模。定义模具间的距离，单击 Tools，Position Tools 中 Move Tools，选择要移动的模具，在 Distance 中输入移动距离，移动方向为 Z Translation。使得 t 在 Z 方向移动 0.5mm，y 在 Z 方向移动 0.5mm,a 在 Z 方向移动-0.5mm。（3）定义凸模的运动和压边力单击 Tools，选择 Define Tools，选择 User Define Tools，选中 t，然后单击Define Load Curve 按钮，Curve Type 选择 Motion，Z 方向，选择 Set Death Time,单击 Auto 按钮，定义凸模的运动，Veloctity（速度），Stroke Dist（冲压行程）负号表示 Z 的负方向，然后单击 OK 即可。定义压边力，选中 y，单击 Define Load Curve 按钮，Curve Type 选择 本 科 毕 业 设 计 说 明 书 （ 论 文 ） 第 21 页 共 36 页Force，Z 方向，选择 Set Death Time,单击 Auto 按钮,定义作用在压边圈上的压边力，Force 输入压边力的大小，打击 OK 即可。（4）对模型进行有限元计算单击 Analysis，选择 LS_Dyna，为了在后处理中能够较好的观察成形过程，一般设定 STEP=20。求解器采用 Full Run Dyna,求解器精度采用单精度。同时在计算机内存较大时，为了加快运算速度，可以适当提高 DYNAFORM 运算器的内存值，即可进行有限元的计算。图 4.3 Analysis 的参数设置4.4 方盒件冲压工艺模拟与后处理分析在仿真结束后， 可以进入 DYNAFORM 的 PostProcess 的后处理环境，进行一系列后处理。（1）等值线对单元应力和相关的结果等进行显示。在同一等值线上是以同一颜色显示的，相应的等值线值在图形窗口右边的颜色柱显示(图4.4(a)。（2）矢量图标将结果以矢量形式显示(图 4.4(b)，能够对材料的流动性进行显示。（3）成形极限图用来评价板料的可成形性(安全和失效区域)。图中每一点的 X 坐标和 Y 坐标代 本 科 毕 业 设