1610.镁合金薄板成形过程有限元模拟——镁合金筒形件拉深成形过程的数值模拟 论文正文.doc

《1610.镁合金薄板成形过程有限元模拟——镁合金筒形件拉深成形过程的数值模拟 论文正文.doc》由会员分享，可在线阅读，更多相关《1610.镁合金薄板成形过程有限元模拟——镁合金筒形件拉深成形过程的数值模拟 论文正文.doc（60页珍藏版）》请在三一文库上搜索。

1、本科生毕业论文 1 毕毕 业业 论论 文文 镁合金薄板成形过程有限元模拟镁合金薄板成形过程有限元模拟镁合镁合 金筒形件拉深成形过程的数值模拟金筒形件拉深成形过程的数值模拟 Magnesium Alloy Sheet Forming process finite element modeling Numerical simulation of Forming process during Deep Drawing of Barrel Type Piece Magnesium Alloy 专专 业：业： 机械设计制造及其自动化机械设计制造及其自动化 班级学号：班级学号： 学生姓名：学生姓名： 指导

2、教师：指导教师： 二二八年八年 六六 月月 本科生毕业论文 2 摘摘 要要 冲压工艺在现代制造业中，特别是汽车、摩托车、轻工、家电和电器仪表 行业应用非常广泛，是目前材料成形及模具技术领域最热门的工艺方法之一。 变形镁合金在以上领域的应用越来越广泛。本文主要研究镁合金筒形件拉深成 形过程的数值模拟。近年来，随着计算机技术和有限元分析技术的飞速发展， 金属板料冲压成型过程的计算机模拟技术正逐渐从实验研究走向实际生产。本 文主要阐述了怎样利用有限元软件 DYNAFORM，对镁合金筒形零件在冲压过程中 出现的缺陷问题进行数值模拟及实验研究，分析了不同参数下的各个变形过程。 分析结果表明，镁合金板料拉

3、深成形的深度越大，越容易拉裂；取合适的压边 力、合适的坯料大小可以有效减少起皱；同时润滑条件对起皱的影响也很大， 有效的润滑可以减少起皱；模具圆角的大小也对板料成形有影响，合理的圆角 大小可以消除起皱、拉裂等现象；坯料相对厚度越小，起皱越严重。通过数值 模拟分析，制定合适的变形参数和有效的润滑条件,可以加强板料成形的合理性。 这对实际生产具有指导意义，从而也可以提高产品质量。 关键词关键词：镁合金筒形零件；冲压；有限元；DYNAFORM 软件 本科生毕业论文 3 引 言 一、 研究工作的目的 镁合金是在工程应用中密度最低的金属结构材料，具有高比强、高比模、 高阻尼、电磁屏蔽以及优异的铸造、切削

4、加工性能和易回收等优点，在汽车、 电子、航空航天、国防等领域具有重要的应用价值和广阔的应用前景，被誉为 “21 世纪绿色工程材料” 。随着很多金属矿产资源的日益枯竭，镁以其资源丰 富而日益受到重视。目前镁及其合金材料的研究已成为世界性热点。我国镁资 源极其丰富，镁工业发展迅速。 本课题基于镁合金的成形机理，利用有限元分析软件 DYNAFORM 来模拟 镁合金筒形件拉深成形过程，分析该过程中产生的起皱现象以及如何优化各个 参数以消除起皱现象，达到完美的拉深效果，更好的应用于实际生产当中，具 有很大的实际应用意义。 二、研究工作的范围以及相关领域 本课题研究的范围广，设计的相关领域颇多，主要包括板

5、料的冲压成型技 术、有限元的数值模拟分析和 DYNAFORM 软件的应用与分析等等。 1、冲压技术发展的特征 近十多年来，随着对发展先进制造技术的重要性获得前所未有的共识，冲 压成形技术无论在深度和广度上都取得了前所未有的进展，其特征是与高新技 术结合，在方法和体系上开始发生很大变化。计算机技术、信息技术、现代测 控技术等冲压领域的渗透与交叉融合，推动了先进冲压成形技术的形成和发展。 纵观上世纪的发展历程可见： 1）冲压性能的研究和改进是与冲压技术的发展相辅相承的。 2）汽车、飞机等工业的飞速发展，以及能源因素都是冲压技术发展的主要 推动力。进入新世纪，环境因素及相关的法律约束日益突出，汽车轻

6、量化设计 和制造成为当前的重要课题。 3）成形过程数字化仿真技术的发展，推动传统冲压技术走向科学化，进入 先进制造技术行列。 4）冲压技术的发展涉及材料、能源、模具、设备等各方面。工艺方法的创 新及其过程的科学分析与控制是技术发展的核心；模具技术是冲压技术发展的 体现，是决定产品制造周期、成本、质量的重要因素。 冲压技术的发展与材料和结构密切相关。预计未来10-15年，环境要求和日益严 本科生毕业论文 4 格的环保法律，将促使汽车材料和结构发生很大变化。为了减少城市的排 2 CO 放量，汽车力求轻量化，其最突出的发展方向是提高所用材料的比强度和比刚 度及发展高效的轻量化结构。现代车身结构中，高

7、强度钢约占25%。目前在继 续开发超高强度钢的同时，结合发展新的“高效结构”和制造技术，争取使车身 重量减少20%以上。但更引人关注的努力方向是扩大铝、镁等低密度合金材料 在汽车上的应用。 2、板料成形模拟技术在工业中的作用 一般说来，板料成形有限元软件可以对薄板成形件的成形结果进行分析、 评价，帮助工程师进行产品工艺可行性研究和模具设计制造，优化板料成形工 艺，优化成形参数，对新材料的评估，对产品的设计开发提出修改建议。有限 元模拟的作用大致可以归纳如下: (1)缩短模具制造周期 从制件的构形设计、模具设计、调试直至投产的整个周期来看，模具制造 是影响整个进度的关键，通过使用有限元模拟技术，

8、能有效提高模具设计，特 别是工艺设计的可靠性，缩短调试周期。通过模拟，找出影响因素(如压边力、 润滑等)的变化趋势，指导模具生产的进行。 (2)节省费用 有限元模拟引人各种材料参数和工艺参数，对冲压工艺进行全面而且有效 的评估，从而减少发生工艺设计不合理的机会，缩短调试周期，减少修模次数， 减少模具设计和制造费用。减少模具报废也是节省费用的一个方面。有时，还 可以取消软模试验。 (3)提高零件的质量和使用性能 保证零件的形状和尺寸精度，有效地控制成形中材料的塑性变形程度，从 而控制材料的塑性硬化程度，改善零件在使用中的力学性能。 (4)降低零件废品率 为零件成形找到合适的生产工艺，保证零件成形

9、的进行。 (5)合理的选用成形材料 根据模拟分析结果，获取生产材料所需的各项力学性能指数，合理的选用成形 材料;用计算机模拟技术进行毛坯尺寸的准确计算，减少不必要的原材料浪费， 从而降低成本。 3、DYNAFORM 软件的功能 DYNAFORM是美国（Engineering Technology Associates Inc.）ETA公司和 LSTC公司联合开发的用于板材成形过程的专用软件。它集成了DYNAFORM本 本科生毕业论文 5 身功能强大的前处理功能和Eta-Post后处理软件以及LSTC公司开发的有限元动 力显示求解器960和970。DYNAFORM的求解器LS-DYNA是采用显隐

10、结合的算 法进行板材成形模拟的最具有代表性的软件，在模拟冲压成形过程中，计算稳 定，效率高。LS-DYNA是目前业界公认的板料成形模拟结果准确性最好的软件 之一。DYNAFORM的后处理功能使用户可以用云图显示板材变形后的应力、 应变信息，材料的厚度分布信息等。作为一款专用的板材成形软件， DYNAFORM可以预测成形过程中板；料的裂纹、起皱、减薄、划痕、回弹， 评估板料的成形性能，从而为板料成形工艺及模具设计工作提供帮助，帮助模 具设计人员显著减少模具的开发设计时间及试模周期。 三、研究工作的主要路线 本课题主要利用 DYNAFORM软件有限元法对板料冲压成形过程进行有限 元法的定量分析和数

11、值模拟，分析镁合金筒形零件成形过程并解决在成形过程 中出现的缺陷。 本科生毕业论文 6 1 文献综述 1.1 冲压技术的现状和发展趋势 近十多年来，随着对发展先进制造技术的重要性获得前所未有的共识，冲 压成形技术无论在深度和广度上都取得了前所未有的进展，其特征是与高新技 术结合，在方法和体系上开始发生很大变化。计算机技术、信息技术、现代测 控技术等冲压领域的渗透与交叉融合，推动了先进冲压成形技术的形成和发展。 冲压技术的真正发展，始于汽车的工业化生产。20 世纪初，美国福特汽车 的工业化生产大大推动了冲术的研究和发展。研究工作基本上在板料成形技术 和成形性两方面同时展开，关键问题是破裂、起皱与

12、回弹，涉及可成形性预估、 成形方法的创新，以及成形过程的分析与控制。但在 20 世纪的大部分时间里， 对冲压技术的掌握基本上是经验型的。分析工具是经典的成形力学理论，能求 解的问题十分有限。研究的重点是板材冲压性能及成形力学，远不能满足汽车 工业的需求。60 年代是冲压技术发展的重要时期，各种新的成形技术相继出现。 尤其是成形极限图（FLD）的提出，推动了板材性能、成形理论、成形工艺和质 量控制的协调发展，成为冲压技术发展史上的一个里程碑。 由于 80 年代有限元方法及 CAD 技术的先期发展，使 90 年代以数值模拟仿 真为中心的和计算机应用技术在冲压领域得以迅速发展并走向实用化，成为材 料

13、变形行为研究和工艺过程设计的有力工具。汽车冲压技术真正进入了分析阶 段，传统的板成形技术开始从经验走向科学化。 冲压成形技术的发展趋势进入 90 年代以来，高新技术全面促进了传统成形 技术的改造及先进成形技术的形成和发展。21 世纪的冲压技术将以更快的速度 持续发展，发展的方向将更加突出“精、省、净”的需求。冲压成形技术将更 加科学化、数字化、可控化。 加入 WTO 以后，中国的汽车工业、航空航天工业等支柱产业必将有大的发 展。我国的冲压行业既充满发展的机遇，又面临进一步以高新技术改造传统技 术的严峻挑战。国民经济和国防建设事业将向冲压成形技术的发展提出更多更 新更高的要求。我国的板料加工领域

14、必须加强力量的联合，加强技术的综合与 集成，加快传统技术从经验向科学化转化的进程。加速人才培养，提升技术创 新能力，提高冲压技术队伍的整体素质和生产企业的竞争力。 1.2 有限元算法和模拟 有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。 本科生毕业论文 7 它是50年代首先在连续体力学领域飞机结构静、动态特性分析中应用的一 种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力 学等连续性问题。 1.1.1 有限元法分析计算的思路和做法可归纳如下： 1、物体离散化 将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。 离散后单元于单元之间利用单元

15、的节点相互连接起来；单元节点的设置、性质、 数目等应视问题的性质，描述变形形态的需要和计算进度而定（一般情况单元 划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大） 。所以有 限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以 一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似 的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。 2、单元特性分析 1）选择位移模式 在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点 力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未 知量时称为混合法。位移法易于

16、实现计算自动化，所以，在有限单元法中位移 法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量 如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用 一些能逼近原函数的近似函数予以描述。通常，有限元法我们就将位移表示为 坐标变量的简单函数。这种函数称为位移模式或位移函数。 2）分析单元的力学性质 根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单 元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹 性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度 矩阵，这是有限元法的基本步骤之一。 3）计算等效节

17、点力 物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是， 对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这 种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去， 也就是用等效的节点力来代替所有作用在单元上得力。 3、单元组集 本科生毕业论文 8 利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来， 形成整体的有限元方程。 有限元模拟算法主要分为：静态隐式和动态显式算法。在实际计算中静 态隐式算法不能保证迭代计算的收敛性而且求解规模巨大，通常用于二维模 拟计算问题。比如本文所用的 CAE 类仿真模拟软件 Dynaform 就是主要以

18、计算二 维、片体等为主而动态显式算法基于动态平衡方程，不需要迭代求解。对于三 维模拟计算问题则更为有效，例如 Deform 软件。模具在成形过程中的模拟分析 一般采用动态显式算法。由于零件的成形过程包含几何材料和接触的非线性问 题。还可以将加载过程的动态显式分析和卸载回弹过程的隐式分析结合起来 进行有效的数值计算。 1.1.2 板料成型过程的有限元模拟： 自从六十年代提出有限元方法，近十年，关于板料成型的数值模拟的研究， 绝大多数采用有限元法。作为一种有效的数值模拟方法，板料成型有限元模拟 的大量研究工作集中在非 线形弹塑性成型问题的定解条件和求解技巧。定解条 件的准确性和求解速度对数值模拟的

19、成败至关重要。有限元方法的灵活性为模 拟各种定解条件提供了可能。 1、材料特性 多数有限元程序处理板料可分为刚塑性与弹塑性两类。还有考虑刚粘塑性 和弹粘塑性材料。塑性变形的本构方程分增量理论、全量理论两种。塑性应力 应变强化模式有线形强化、分段线形强化、指数强化等。 2、单元类型 根据板料冲压成型的力学特点，有限元模拟采用的单元基本上是薄膜板壳 元、弯曲板壳元或两种单元组合。其单元又有二维、三维之分，偶尔也有用实 体元的，模拟对象多为一般轴对称件、半球件、盒形件。 3、边界条件 冲压模具、工件压紧装置的几何形状，外力、摩擦条件、接触条件等边界 条件的准确性是有限元模拟成败的关键之一。 接触问题

20、长期为数学力学工作者所关心，求解途径有二：直接法和数学规 划法。第一种方法在有限元方法中直接设定接触边界条件，通过反复试算获得 解答；第二种途径在给定的接触边界约束下求解势能最小的解，解这种数学规 划问题的方法又分二次规划法、罚函数法、拉格朗日乘子法等几种。 在成形过程中，有的不考虑摩擦问题(模具与夹具间)，有的采用简单的库 仑摩擦定律将摩擦条件转换为接触条件，还有考虑相邻单元之间的摩擦力的。 一般说来，板料成形时冲压模具可认为是刚性的，因此模具与板料间的接 本科生毕业论文 9 触可视为变形体与任意形状刚体之间的接触 。 4、网格划分 有限网格可以在 Dynaform 软件中进行自动的划分。

21、5、求解方法 成功地用于简单模具的方法对复杂问题不一定总是有用，即便解答得到了， 精度还经常成问题。另外，解答的有效性和强度也非常重要，问题甚至不能简 单地用先进的计算机来解决，许多研究是瞄准模拟算法的精度、强度、有效性， 从而创造条件使复杂板料成形模拟变换更可靠。Dynaform 软件中自带的求解器 功能非常强大，较之以前的一些有限元分析软件功能强大了很多，计算结果也 比较接近于实验数值。 6、成形趋势 板料成形的数值模拟就是在定解条件下通过数值求解预测并显示破裂、起 皱、回弹、板料厚度变化等成形趋势。对成形趋势的准确性直接标志数值模拟 的好坏，从塑性成形理论看，预测成形趋势的最重要的因素是

22、应变分布，破裂 判据、起皱判据、回弹、厚度变化等都是在应变分布基础上作出直接或合成的 判断。 1.3 DYNAFORM 的应用 以下结合生产中的实际，举例说明利用 DYNAFORM 软件进行的拉深成形数值 模拟。 对于非规则的拉深件而言，坯料形状对拉深成形中料的流动阻力有很大的 影响，同样的工件，由于坯料几何形状不同，在成形时产生的应力状态有很大 的区 别。所以，如何优化坯料的形状显得尤为重要。下面是某公司生产的 NSW 型手制动机箱壳是非规则的拉伸成形件，见图 2 所示。手制动机壳体采用厚度 为 4mm 的 Q235 材料，在中厚板成形工艺中，它属于非规则的拉伸件。从图 1 可以看出，该件包

23、含盒形拉伸和圆锥形拉伸两种变形，在成形中容易出现断裂 现象。 本科生毕业论文 10 通过实验测得材料性能参数扬氏模量 E=207GPa；泊松比=0.3；屈 服应力 =36483Mpa；硬化指数 n=0.186；板厚=4ram；各向异性系数 s R=0.747；凸凹模与板材之间的摩擦系数 f=0.15，压边圈与板材之间的摩擦 系数 f=0.1。模具部件及板料的几何模型在 ProE 中创建。由 DYNAFORM 的 IGES 接口读入，进行前处理。手制动机壳体其数值模拟有限元模型的单元个数 为 7146 个，节点个数为 7189 个，其中板材的单元个数为 2578 个，节点个数为 2652 个。图

24、 3、图 4 为板料优化前的厚度变化云图和成形极限云图(FLD)；图 5、图 6 为板料优化后的厚度变化云图和成形极限云图。 本科生毕业论文 11 从模拟结果可以知道，板料形状优化前，由于圆锥拉伸部分材料流动变形 阻力较小，而盒形角部进料阻力较大，按成形最小阻力原理，盒形角部板料流 动变形较少，导致底部圆角处开裂，该处厚度变薄达 33，从 FLD 曲线看，该 区域已经超过临界线，即使压边力降到 16t，模拟结果变化不大。板料形状优 化后，圆锥拉伸部分的板料尺寸适当加大，以平衡盒形角部进料阻力，压边力 给定为 26t，从模拟结果看，盒形底部圆角处厚度变薄为 19，从 FLD 曲线看， 本科生毕业

25、论文 12 该区域处于临界线内。 1.4 镁合金 1.4.1 镁合金的特点： 1）重量轻，镁金属是目前世界实际应用中重量最轻的金属结构材料。比重 只有 1.8，分别为铝的1/3和铁的1/4。 2）镁合金相对比强度（强度与质量之比）最高。比刚度（刚度与质量之比） 接近铝合金和钢，远高于工程塑料。其比强度高达133，可以和钛的比强度相媲 美。 3）在弹性范围内，镁合金受到冲击载荷时，吸收的能量比铝合金件大一半， 所以镁合金具有良好的抗震减噪性能。弹性模量小，刚性好，抗震力强，长期 使用不易变形。 4）镁合金熔点比铝合金熔点低，压铸成型性能好。镁合金铸件抗拉强度与 铝合金铸件相当，一般可达250Mp

26、a，最高可达600多Mpa。屈服强度，延伸率与 铝合金也相差不大。 5）镁合金具有良好的压铸成型性能，压铸件壁厚最小可达0.5mm。适应制 造汽车各类压铸件。 6）镁合金件稳定性较高压铸件的铸造行加工尺寸精度高，可进行高精度机 械加工。 7）电磁屏蔽性能，防辐射性能好。对环境无污染。 镁合金作为工业应用中最轻的金属结构材料，由于具有密度低 、比强度和 比刚度高、抗震及减振能力强、电磁屏蔽效果优异、尺寸稳定、导热性好 、机 械加工性能优良及易回收等优点，在电子、电器、通讯 、汽车、交通、航空、 航天等领域具有极其重要的应用价值和广阔的应用前景。但是，镁合金的密排 六方结构，传统上被视为一种塑性加

27、工性能差的金属材料，所以其加工方式局 限于铸造特别是压铸方面，然而，铸造件存在力学性能差，壁厚受到限制、易 产生缺陷等缺点，大大限制了镁合金的应用范围和使用性能。所以，对其塑性 加工技术的研究已成为镁合金研究和发展的方向，并已取得了一定的成果，通 过挤压、轧制等塑性加工工艺加工的镁合金板材具有更高的强度、更好的延展 性和更多样化的力学性能，这就为镁合金冲压工艺的开展奠定了基础。 随着市场对镁产品需求应用领域的不断拓宽，从航空、航天、汽车零部件、 钢铁脱硫、合金压铸件、3C 产品的广泛应用，到民用产品的不断研发，以及镁 合金技术的进一步研究，镁产品的发展愈来愈显现出它独特的不可替代优点。 同时由

28、于镁合金较好的加工性能、较强的机械性能、重量轻、韧性好、屏蔽性 本科生毕业论文 13 好，越来越受到人们的青睐。 1.4.2 镁合金塑性成形的特点 镁合金的塑性加工性能主要与变形温度、合金的成分及组织、变形速率、 成形工艺等影响因素密切相关。镁合金的塑性成形主要有以下几个特点： 1）对加工温度敏感 温度是影响镁合金的塑性变形能力的关键因素，镁合 金的塑性成形温度区间较窄，一般在 300-450的区间内。当成形温度低时， 镁合金的变形能力差，易产生冷裂纹。温度过高时，合金中的某些低熔点相可 能在成形的过程中熔化，在应力作用下易产生热裂纹。此外，成形温度过高， 也会造成零件显微组织的晶粒粗大，直接

29、影响材料服役的性能。 2）合金成分与显微组织的不同，镁合金的塑性成形工艺条件也不同 一般 说来，随着合金元素含量的增加，镁合金的塑性变形能力变差，从而导致其塑 性成形能力下降。原始显微组织粗大也可造成材料的塑性成形能力的降低。 3）加工速率慢 镁合金的塑性变形能力差，在热加工过程中需要动态回复 与动态再结晶来提供其持续变形的能力。而动态回复和动态再结晶的进行程度 与时间有关，这就决定了镁合金的塑性加工速率比较慢，生产效率低。例如， 镁合金的轧制速率一般在 10 m/min 的量级，约是钢的热轧速率的 1/10 左右。 铝合金、低碳钢、铜合金等可以进行高速的锻造，而镁合金的锻造主要是在慢 速的压

30、力机上进行。增大塑性加工变形的速率，极易导致表面裂纹。 4）变形抗力大 在塑性加工的温度范围内，镁合金比铝合金具有更高的屈 服强度。 5）流动性差 镁合金的黏性很大，造成镁合金填充深空模腔的能力要差得 多。对于铝合金要两付模具可以完成的锻件，如果用镁合金制造同样的锻件则 需要三付模具。另外，模具内表面要进行抛光、润滑来促进加工中的金属流动。 6）镁合金的导热性能好 在加工中的热量流失快,如在轧制加工中，镁合 金坯料的热量损失过快，需要加热轧辊及道次间隙加热坯料来保证塑性变形能 力。 7）制件具有力学行为的各向异性 在塑性加工过程中，由于晶内滑移和孪 晶仍起着重要作用，使得再结晶晶粒出现择优取向

31、从而形成结构。镁合金是 hcp 晶体结构，最常见的变成组织结构为0001基面结构。由于结构的存在， 使得镁合金塑性加工后在力学性能上具有各向异性。根据实验结果，可以看出， 对于挤压件来说，在挤压方向上材料强度最高。而在垂直于挤压方向的界面上， 材料的强度较低。另外，材料的强度还与受力状态有关。一般来说，塑性加工 后材料受拉时的屈服强度要比受压时屈服强度高得多。因此，在零件设计时要 本科生毕业论文 14 充分考虑结构的影响。 1.4 镁合金拉深工艺发展趋势 目前，镁合金拉深工艺的研究已经取得了一定的成果，但主要局限在较为 传统的冲压工艺，针对其拉深成形性能进行研究。为了进一步提高镁合金拉深 成形

32、性能，扩大镁合金的应用范围，今后可在以下几个方面展开研究： 1)根据已有的对其它材料拉深工艺的研究，结合镁合金的特点，研究镁合 金板料的化学成分、组织与制造过程等因素所决定的材料性能与拉深成形之间 的关系，为镁合金材料的设计与生产提供技术依据与要求。 2)研究开发新的拉深成形工艺或运用已开发的先进的拉深成形工艺，以此 来提高镁合金的拉深成形能力与生产效率。利用一些先进的拉深技术，如变压 边力控制技术，温度分布控制技术、切分式压边圈辅助摩擦技术、高温液压成 形技术、气胀成形等，并与计算机技术相结合形成的智能化冲压技术已经成为 冲压工艺发展的重要方向与趋势，而由此也必将极大的促进镁合金拉深技术进

33、步与提高。 3)结合有限元(FEM)分析技术，研究镁合金拉深变形行为，预测其成形能力， 为镁合金拉深模具的设计、拉深工艺的制定提供依据和指导，降低生产成本， 提高生产效率。 4)根据镁合金板料特征和零件要求，对拉深模具的几何形状进行优化设计， 以提高其成形能力，并制造出结构更为复杂、精度更高的镁合金拉深件。 5)对镁合金拉深变形过程中组织结构、力学性能等的变化以及各工艺参数 对变形后 镁合金产品组织、性能等的影响进行研究，以明确成形后镁合金零件 的性能是否满足要求并根据产品使用性能要求进行拉深工艺的设计与优化。 2 实验过程 2.1 简单校核模具尺寸 本文选取带凸缘的筒形件为拉深实验件，首先需

34、要确定各个部分尺寸。 给定拉深件直径为 50mm ,材料是牌号为 AZ31B 的镁合金薄板，冲压方向为 90方向轧制，板料厚度为 t=1mm。 2.1.1 拉深模的凸凹模圆角尺寸确定 查冷压模具设计与制造得，带凸缘件拉深凹模圆角半径为=（6-10） d r t，对于本实验=6mm。则凸模半径= (0.61.0)rd=1.0t=6mm。 d r p r 本科生毕业论文 15 857 . 0 7 6 R r m t z 1 . 11 2 1 . 11 . 11 . 11 2 tt z 2.1.2 拉深次数的确定 圆角部位拉深系数 式中，r：拉深件圆角半径，。mm R：圆角部分的坯料半径，。mm 由

35、毛坯相对厚度可查高性能变形镁合金及加工技术书中图 6-36 拉伸 实验得极限拉深系数。因为 ， 所以，只 8 . 077 . 0 1 m mm8 . 077 . 0 1 需要拉深一次。 2.1.3 是否使用压边圈的确定 略估算普通平面凹模拉深时的不起皱条件： 首次拉深： )1 (045 . 0 1 m D t 式中， ：材料厚度，；tmm ：毛坯直径，；Dmm ：首次拉深系数； 1 m 则 所以，02106 . 0 )532 . 0 1 (045 . 0 )1 (045. 001064 . 0 94 0 . 1 1 m D t 在拉深过程中会起皱，需要加压边圈。 2.1.4 拉深模的间隙确定

36、拉深模的凸凹模之间的间隙对拉伸力、零件质量、模具寿命等都有影响。 间隙小，拉伸力大，模具磨损大；但冲压件回弹小，精度高。间隙过小，会使 拉深件严重变薄甚至拉裂。间隙过大，坯料容易起皱，冲压件锥度过大，精度 差。因此，应根据板料厚度及公差、拉深过程板料的增厚情况、拉深次数、拉 深件形状及精度要求等，正确确定拉深间隙。 在该实验条件下，在拉深时会发生起皱现象，所以需用压边圈，查冷冲 压成形工艺与模具设计制造书中表 616 得 式中， ：拉深模的凸凹模间隙，；zmm ： 毛坯厚度，；tmm ：系数，精度要求较高的拉深件取小值，精度要求较低得拉深件取大值。1 . 11 本科生毕业论文 16 由于 t=

37、1.0mm，则 即单边间隙为 1.1mm。 2.1.5 压边力的确定 根据经验公式 qrdDQ d )2( 4 22 其中 ：压边力，QN ：毛坯直径，mmD ：拉深件直径，mmd ：凹模圆角半径，mm d r ：单位压边力，qMpa 则 = 其中，根据简单件单Q44.78370 . 2)6250(94 4 22 Mpa 位压边力一般为 2，所以取 2。MpaqMpa 2.2DYNAFORM 软件模拟 2.2.1 DYNAFORM 软件使用详细步骤： 在应用 DYNAFORM 软件板料成形过程时主要包括三个基本部分，即建立计算 模型、求解和分析计算结果。 具体应用步骤表述如下： 直接在 DYN

38、AFORM 的前处理器中建立模型或在 CAD 软件（如 UG、CATIA、PRO/E 等）中根据拟订的成形方案，建立板料、对应的凸模和凹模 的型面模型以及压边圈等模具零件的面模型，然后存为 IGES、STL 或 DXF 等文 件格式，将上述模型数据导入 DYNAFORM 系统。 利用 DYNAFORM 软件提供的网格划分工具对板料、凸模、凹模、压边圈进 行网格划分，检查并修正网格缺陷（包括单元法矢量、网格边界、负角、重叠 结点和单元等） 。 定义板料、凸模、凹模和压边圈的属性，以及相应的工艺参数（包括接 触类型、摩擦系数、运动速度和压边力曲线等） 。 调整之间的相互位置，观察凸模和凹模之间的相

39、对运动，以确保模具动 作的正确性。 设置好分析计算参数，然后启动 LS-DYNA 求解。 将求解结果读入 DYNAFORM 后处理器中，以云图、等值线和动画等形式显 示数值模拟结果。 分析模拟结果，通过反映的变化规律找到问题的所在。重新定义工具的 本科生毕业论文 17 形状、运动曲线，以及进一步设置毛坯尺寸，变化压边力的大小，调整工具移 动速度和位移等，重新运算直至得到满意的结果。 2.2.2 举例说明模拟起皱现象的主要步骤与分析 已知成形件直径=50mm ,选择凸模直径=24.45mm，凹模直径d p d =25.55mm，凸模圆角=6mm，凹模圆角=6.5mm,拉深高度 h=25mm，则可

40、以 d d p r d r 根据以下公式计算出毛坯的直径。 D drdhdD t 44 . 3 4 2 其中，：凸缘直径，mm；：成形件直径，mm；：拉深深度，mm； ： t ddhr 圆角半径，mm 因为模拟成形的是窄凸缘筒形件，所以根据查冷冲压成形工艺与模具设 计制造选择，因为=50mm，所以=70mm。4 . 1 d dt d t d504 . 1 则， =94mmdrdhdD t 44 . 3 4 2 65044 . 3 25504702 1）模型的建立 在 CAD 软件 UG 中建立模型，包括凸模、凹模和毛坯。 图 2-1 在 UG 中建模 （图示，上端是凹模 PUNCH，中间是毛坯

41、 BLANK，下端是凸模 DIE。因为在 UG 中不容易生成压边圈，所以压边圈在 DYNAFORM 系统中生成。 ） 2）模型的导入和网格划分 本科生毕业论文 18 在 UG 中把 BLANK、DIE、PUNCH 分别生成 IGES 格式的文件，导入 DYNAFORM 系统。之后根据 DYNAFORM 系统中的工具菜单栏中 PreprocessElement 进行 凸、凹模和毛坯的网格划分。之后再进行网格检查，如果划分的网格有漏洞， 则需要修补或重新划分网格。 因为本实验做的是带凸缘的筒形件，属于简单件的模拟，选用的板料毛坯 为圆形，所以使用软件自动划分网格的功能即可。但是网格需要划分多大，是

42、 一个需要解决的问题。因为网格划分原则是越小越好，所以经过前期熟悉软件 的过程，我做了大量试验，结果得出网格划分为 1 是最合理的。所以之后的整 个试验过程，我均选用网格划分为 1。 例如，板料厚度为 1mm，拉深深度为 25mm，坯料直径为 50mm 的镁合金筒形 件拉深数值模拟过程中的网格划分。板料的网格划分，单元个数为 6979，节点 个数为 7056；凹模网格划分，单元个数为 10335，节点个数为 18633；凸模网 格划分，单元个数为 10511，节点个数为 17490；压边圈的网格划分，单元个数 为 3890，节点个数为 29242。 图 2-2 凸、凹模和毛坯的网格划分 3）生

43、成压边圈 利用凸模 PUNCH 的凸缘通过菜单栏上的工具按纽 Preprocess Element Copy 可以生成压边圈。 本科生毕业论文 19 图 2-3 压边圈的生成（绿色部分为生成的压边圈） 4）LS-DYNA 970 运算器运算 合理定位之后，定义毛坯、凸模、凹模和压边圈的属性。测量凸、凹模之 间 的距离、凸模和压边圈之间的距离，之后定义凸模和压边圈的运动。 图 2-4 凹模 DIE 的运动加载曲线 本科生毕业论文 20 图 2-5 压边圈 BINDER 上压边力的加载曲线 凸模运动利用速度控制。实际生产中使用的油压机速度低，凸模实际速度 约为 2mm/s 左右。由于 DYNAFO

44、RM 采用显式算法，对速度不敏感。为加快 运算，可以直接增大时间步长的值，或提高模具的虚拟运动速度。研究表明， 当虚拟速度是实际速度的 1000 倍左右时，模拟结果的对误差较小一。所以这 里的虚拟速度取 2000mm/s。 压边圈运动的控制可以分为两个阶段，第一个阶段是利用速度和力一起控 制，第二个阶段是利用力来控制。结果证明后者对于解决起皱现象是合理的， 因为利用速度和力一起来控制的话，速度控制是几乎不起作用的，凹模与压边 圈之间预留的间隙在模拟过程中无法保证始终如一，所以单用压边力来控制。 压边圈运动利用力控制。最开始做实验时，压边力的大小是根据经验来选 取的，经过查阅大量资料之后，得到了

45、一个经验公式，根据此经验公式算出的 力的大小，经过软件模拟及反求验证，结果证明是可行的，所以压边力的大小 就选取利用经验公式算出的力的大小。 完成以上设置后，保存文件，由板料成形数值模拟软件生成求解器的输入 文件(如 LSDYNA 的 DYNA 控制文件 *Dyn 和 DYNA 模型文件 *mod)， 提交给 DYNAFORM 有限元动力显示求解器(LSDYNA 970 求解器)进行求解。 5）Eta/Post GL 后置处理 LSDYNA 利用*dyn 和*mod 两文件的数据，按照一般有限元运算步骤 迭代计算，输出各个时间步后处理所需的数据文件 d3plot。 Eta/Post GL 读进

46、 d3plot 文件进行后处理，它可用来生成各个时间步的板 料变形图和成形极限图(FLD)；级联菜单 View/Profile(查看/分布)能以表格形 本科生毕业论文 21 式列出多个节点在所选时间步的位移、正(剪)应力和应变、弯矩、扭矩、剪力 或板厚等物理量的数据；菜单 View/List contour value(查看/ytj 出等值线值)可 列出所选节点在某时间步的位移、速度、加速度或其他物理量的数据，画出板 料在整个变形过程中该物理量的时间曲线图，以及找出使它取最大值的节点号 和所处工件的部位；菜单 View /Nodetrace(查看/节点路径)可图示一个或多个 节点当板料连续变形

47、时所经历的路径；而 Eta/Graph 则用曲线显示 d3plot 文件 中各物理量的数据对成形仿真结果进行评估，通过后处理器查看工具的运动、 位移、成形过程中系统的能量变化以及工具与板料之间的接触等等以判断仿真 的可靠性。 6）FLD 成形极限图说明 图 2-4 FLD 成形极限图 如图所示，其显示色系的意义如下： 1)红色区域破裂区域 2)黄色区域有拉裂趋势的区域 3)绿色区域一一正常变形过程（属安全变形） 4)蓝色区域拉压变形区域（安全变形） 5)粉红色区域起皱区域 相邻两色之间的区域在图中表示临界区，对应的曲线则表示临界曲线。依 照如上所述可以预见工艺规程的危险性，如果板料单元实际变形

48、落于临界区， 则说明很危险，废品率较高；如果实际应变值落在临界曲线上，说明有相当的 危险，对各条件须严格控制；如果落在远离临界曲线的地方，则说明过分安全， 板料还有变形的潜力。 将求解后得到的 d3plot 文件利用软件的 Post-Processor 程序进行后续分 析，由图 3-4 最终变形的 FLD 图可知，板料单元进行冲压成形时，其大多数变 形区落在安全区内，一小部分起皱，一小部分有拉裂趋势，有一点破裂。 经 Dynaform 软件对镁合金筒形件拉深成形过程的仿真模拟分析可知，以筒 形件的结构设计相应的拉深模具，是能够得到符合要求的冲压件的，至此，关 本科生毕业论文 22 于零件及其拉

49、深模具的结构论证工作结束。 3 不同试验结果分析比较及讨论 3.1 带凸缘筒形件的拉深变形特点 带凸缘的筒形件，其凸缘未被全部拉入凹模，故其变形程度不能用 d/D 来 反映，根据凸缘直径的大小、拉深高度的大小不同，分以下几种情况。 第一种情况 在拉深过程中，拉深力始终未达到传力区筒壁的极限负荷， 窄凸缘、小高度的工件拉深属于这类。这类凸缘筒形件可一次拉成。 第二种情况 拉深过程在拉深力上升达到传力区筒壁的极限负荷前停止了， 宽凸缘、小高度的工件拉深属于这类。这类带凸缘筒形件可一次拉成。 第三种情况 拉深时拉深力超过传力区筒壁的极限负荷，不能一次拉深，需 要经过多次拉深。窄凸缘、大高度和宽凸缘、大高度两类带凸缘件都可能出现 这种情况。 本文研究的是窄凸缘、小高度的镁合金筒形件的拉深成形数值模拟过程，所 以属于第一种情况，这个结论与