外缘内曲翻边工艺与数值模拟分析

1、外缘内曲翻边工艺及数值模拟研究1绪论1.1 外缘翻边及分类翻边是将金属平板坯料或半成品工序件的一部分，沿其一定的轮廓线使其内法兰部分变大、成为有竖边边缘的冲压成形方法。当然，也有不变成竖边，只把坯料中某一部分的孔径加以扩大的情况有的书上称为扩孔)。翻边主要用于制造出与其他零件的装配部位或者为了提高零件的刚度而用来加工出特定的形状，在大型板金成形 时，也可作为控制破裂或折皱的手段。在汽车覆盖件中，多数零件都要经过翻边加工， 翻边部分主要用于各覆盖件之间的连接(焊接、铆接、粘接等)，有时翻边也是产品流线或者美观方面的需要。圆孔翻边、内曲外缘翻边等属于伸长类翻边；外缘外曲翻边等属于压缩类翻边， 而非

2、圆孔翻边通常是伸长类翻边、压缩类翻边和弯曲成形的组合形式。当翻边的变形区边缘为一直线时，翻边成形就转变为弯曲成形，可以说，弯曲是翻边的一种特殊情况。但翻边时，金属板料变形状态比弯曲时要复杂的多。外缘翻边分为内曲翻边和外曲翻边两种。外曲翻边属于压缩类翻边，其变形性质和应力状态类似于不用压边圈的浅拉深。内曲翻边属于伸长类翻边，与孔的翻边相似。当翻边变形面为平面时为平面翻边；翻边变形面为曲面，称为曲面翻边。这里主要讨论内曲外缘翻边的情况1-3。沿不封闭的内凹曲线进行的翻边成形称为内曲外缘翻边，如图 1.1 所示，代表了伸长类翻边成形的一般形式，其变形区的应力状态和变形特点都与圆孔翻边相似。两者不同之

3、处在于内曲翻边时，竖边根部各处应力、应变分布不均匀，曲率半径不同的部位应力应变状态也各不相同。因此，在制备成形坯料时，必须考虑到这一变形的不均匀性，使坯料外轮廓的凸凹变化与材料变形程度大小取得一致。内曲翻边变形程度可用参数来表示：即KA=f/(R-f)(1)KA 值越大，内曲翻边的变形程度越大。变形区的切向拉应力和拉应变在竖边边缘达到最大；并在接近翻边变形结束时刻，变形发展到最为严重程度。外曲翻边成形时的变形区材料承受切向压应力的作用，大到一定量值后，会使变形区材料发生压缩失稳，因此，变形区材料起皱是外曲翻边成形顺利进行的主要障碍。在实际生产中，通常都须对不参与变形的坯料突缘平面施加足够大的压

4、边力。图1.1 内外缘翻边示意图主要特点：圆孔翻边及沿内凹曲线的翻边为典型的伸长类翻边。翻边时，带有圆孔的环状坯料周边被压边圈压紧，位于凹模孔口以内的坯料部分得到塑性变形。在变形过程中，随着凸模下降，变形的坯料单元体在凸模和凹模的圆角处发生弯曲，孔的直径增大，而坯料单元体相对凸模作径向移动，逐渐靠近凹模腔内壁，在此移动过程中，坯料单元体在径向平面内，由凸模下方平直形状移到凸模圆角时产生弯曲，在随后离开凸模圆角时，再次由弯曲变直4-8。1.2 翻边工艺的发展简史及现状我国在板料成形数值模拟方面起步较晚，与发达国家(美、日等)相比晚了十几年。经过多年的发展，我国在板料成形数值模拟方面已经取得了很大

5、进展，但主要集中在部分高校中。华中理工大学针对不完全对称盒形件的成形特点，开发了有限变形弹塑性薄膜有限元程序对其进行分析研究。吉林工业大学采用更新的 Lagrange 法以及有限元变形虚功率增量型原理的弹塑性大变形有限元法,研究了金属板料成形的塑性流动规律以及成形过程中发生的起皱、裂纹等现象，首次提出了多点成形时非连续接触边界约束的处理方法，建立了基于 Mindlin 壳理论三维金属板料成形过程分析的有限元模型，编制了用于板料多点成形分析的有限元专用软件，成功分析了多点成形时的金属流动规律。哈尔滨工业大学采用刚粘塑性本构关系，开发了粘塑性板壳成形有限元分析程序，并对方盒件的成形过程进行了分析，

6、对板料粘性介质胀形过程中的应变速率变化也进行了模拟研究。北京航空航天大学则对板料成形过程中的接触摩擦和悬空区起皱进行数值模拟。上海铁道学院的李尧臣用有限元法模拟了金属板材的冲压成形过程，分析了金属板材在冲压过程中的屈曲现象，建立了增量形式的变分原理，跟踪了板料起皱的发展、折叠、衰减的全过程。上海交通大学对板料成形的回弹进行了较为系统的研究，提出在板料回弹模拟中采用修正的拉格朗日法较为合适9-10。卢险峰教授通过模拟以及实验验证，得出将翻边变形过程可分为具有不同特征的4 个阶段。即变形区材料从变形开始出现斜角段斜角与直壁段圆角与直壁段要求高度的直壁段的几个具体阶段。这使得对翻边变形过程的认识从现

7、有的宏观定性分析进入到细观定量分析。得出采用那些用翻边变形中料厚不变的假设条件所进行的理论分析是不甚合理和不够精确的。数值模拟对深入认识翻边变形过程及其变形特点, 对翻边工艺的理论研究和工程实践具有指导意义11-14。对于翻边的研究，国外学者 Wang N.M.和 Wang C T.曾经采用完全解析的手段， 分别提出了内曲翻边和外曲翻边成形的变形分析计算模型，可以用来分析翻边零件的变形程度，却不便于事先确定修边线的位置，公式推导繁琐，而且没有考虑材料各向异性对翻边的影响。对于翻边回弹方面的研究有 H.Livatyali 等利用 CAE 技术对直翻边中的回弹现象进行了预测，并提出了相应的消除方法

8、。M.J.Finn 将 LS-DYNA 和LS-NIKE3D 结合起来模拟了福特某轿车前翼子板的成形回弹问题，为求解大型复杂外缘成形回弹问题提供了一种简便有效的方法。国内学者曹颖、李大永等在 Wang N.M. 和 Wang C T.的基础上根据全量理论和膜应变假设建立了各自的翻边解析理论模型， 用来确定翻边零件的坯料尺寸，但他们的模型只能用来确定简单地翻边零件的毛坯尺寸15-16。毛坯形状和尺寸的确定一直是板料成形技术中的一大难题，由于板料变形的复杂性，想获得理想的毛坯形状并不容易。毛坯外形的确定问题于 20 世纪 50 年代就已提出，早期提出的有查图和查表法、实验逐次逼近法等，后来又提出了

9、拼合法、滑移线法、物理模拟法、几何映射法和有限元法等。目前应用较多的方法有滑移线法、几何映射法、基于有限元的试错法、反向法、理想成形法、回溯法、增减体积法、类比法、约束优化法、迭代法、人工网络技术等等。1.3 板料成型工艺研究的主要方法简介1.3.1 工艺实验法工艺实验法基本过程是: 首先根据实验目的进行工艺设计，然后在工艺设计的基础上设计并加工制造出模具，并用生产用实际坯料材料来进行实验，然后根据实验情况，检验所要求的变形能否达到零件的设计要求，如果不能，则需反复试验直至达到要求为止。另外，根据实验目的的不同，实验过程不尽相同。该方法的优点是所得的结果准确可靠，缺点是成本较高，且实验周期较长

10、。通常，在科学研究中，工艺实验法是一种极其重要的研究方法，经常被采用， 并且这种方法是检验其他方法所得结果是否合理的重要依据。例如，在翻边工艺研究中，我们可以通过它来进行理论方面的研究，诸如坯料设计计算问题、翻边模具设计问题、应力应变分析以及翻边力的计算等一系列问题都可以通过这种方法来进行研究17-19。1.3.2 理论计算与分析方法理论计算是塑性力学解析的最精确方法。这种方法是将平衡微分方程和塑性条件进行联解，以求出物体塑性变形时的应力分布和应变状态，进而求出变形力。在联解过程中，积分常数根据自由表面和接触表面上的边界条件确定。必要时还须利用应力与应变的关系式和变形连续方程。在一般情况下，共

11、有 16 个未知数，需联列求解16 个方程。但实际上同时联立求解 16 个微分方程是相当困难的，通常只能求解一些简单的二维问题。为此，不少学者在如何简化求解方面作了各种努力，导出了许多近似的求解方程和实验计算法(如主应力法20-21滑移线法、变形功法22-23和极限分析法24-25等)。1.3.3 数值模拟法在工程技术领域内，对于力学问题或其他场问题，已经得到了基本微分方程和相应的边界条件。一般说来，微分方程的边值问题只是在方程的性质比较简单，问题的求解域的几何形状十分规则的情况下，或是对问题进行充分简化的情况下，才能求得解析解。而实际的材料成形问题求解域往往是十分复杂的，而且场方程往往相互祸

12、合， 因此无法求得解析解，而对问题进行过多简化所得到的近似解可能误差很大，甚至是错误的。因此，复杂工程问题的求解必须采用数值方法。有限元法是一种新的现代数值方法。它将连续的求解域离散为有限个单元组成的组合体。这样的组合体能用来模拟和逼近求解域。有限元法另一重要步骤是利用在每一单元内假设的近似函数来表示全求解域上未知场函数。单元的近似函数通常由未知场函数在各个单元节点上的函数值以及单元插值函数表达。因此，在一个问题的有限元分析中，未知场函数的节点值就成为新的未知量，从而使一个连续的无限自由度问题化为离散的有限自由度问题。一经求出这些节点未知量，就可以利用插值函数确定单元组合体上的场函数。显然，随

13、着单元数目的增加，即单元尺寸的缩小，解答的近似程度将不断改进。如果单元满足收敛条件，得到的近似解最后将收敛于精确解26-27。有限元法的基本思想是把连续体视为离散单元的集合体来考虑。在应用有限元法分析问题时，首先采用 “化整为零”的办法，将连续体分解为有限个性态比较简单的“单元”，对这些单元分别进行分析;然后采用 “积零为整”的办法，将各单元重新组合为原来的连续体的简化了的模型，通过求解这个模型(例如位移)在若干离散点上的数值解;最后，根据得到的数值解再回到各个单元中计算其他物理量(例如应变、应力) 28。有限元法起源于 20 世纪 40 年代提出的结构力学中的矩阵算法。它起初是作为一种力学分

14、析的数值计算方法，后来发展成为求解偏微分方程边值、初值问题的一种一般的离散化方法。有限元法是目前进行非线性分析的最强有力的工具29-30。在工程实际问题的分析、解决中，有限元法将发挥着其日益重要的作用31-32。在塑性成形过程的有限元模拟中，根据材料本构关系的不同可将有限元法分为小变形弹塑性有限法、有限应变弹塑性有限方法、刚塑有限元法和弹塑性有限元法33。1.4 课题来源和研究意义1.4.1 课题来源前已述及，翻边是拉伸类冲压成形的基本工序之一，在生产实际中，有很多要用到翻边工艺方法制造的零件。然而，现今的研究中，关于翻边工艺方面的研究相对比较少，与翻边有关的文献资料更为有限。因此，有必要对翻

15、边工艺进行更多、更全面和更深入的研究。本课题的研究内容和方向就是在翻边的发展历史和研究现状的基础上，针对在翻边工艺中研究不够深入全面的地方，或者说比较薄弱的地方而提出的。本课题来源是导师研究课题，也是江西省 “十五”重点学科建设项目 材料加工工程建设项目的内容之一。涉及内容主要有翻边变形中材料厚度变薄规律、冷变形硬化规律、基于Pro/E 模型分析的翻边预加工小孔孔径的设计计算以及开口翻边工艺研究等。这些问题中有的内容前人己有了一定的研究但不够深入，有的内容甚至属于空白。因此，探讨和研究这些内容都是具有较大的理论价值和实践意义，这些内容的预期研究结果必将对翻边工艺乃至冲压工艺理论的完善和对冲压实

16、践水平的提高产生积极的影响34-36。1.4.2 研究意义传统的外缘翻边成形坯料的确定主要依赖于经验和直觉，诸多因素是通过生产的多次反复试验来调整的，结果很难达到现代生产的要求，产品的质量也不易保证。随着计算机技术的不断进步以及有限元技术的不断发展。近年来发展了采用有限元法对板材成形过程进行计算机模拟和分析的新技术，在工程实际中得到迅速地应用和发展。近年来，在板料冲压成形的相关行业中，采用各种仿真软件进行外缘翻边成形的坯料求解和优化，有着极为广泛的应用，这些表明： 借助于国内外成熟的有限元模拟软件与试验结合进行技术研究在技术和方法上是可 行的37-3。81.4.3 研究内容本课题首先系统地归纳

17、和总结了冲压成形中主要的模拟方法，并以典型的外缘翻边件为研究对象，进行模拟比较，分析坯料相关的因素，总结出坯料求解与成形结果之间的关系；从而优化相关参数，改进成形的工艺，使实际生产中成形可能出现的问题在成形前得到预测并给予解决，继而推广到其它零件，给类似的外缘翻边成形提供有价值的参考。本课题以翻边成形中的坯料变形和工艺参数影响为研究对象，主要研究内容如下： 外缘翻边零件的坯料研究。通过数值模拟软件 Dynaform 对坯料的相关因素进行分析，从模拟结果总结出规律，并进行坯料优化。 坯料变形过程中的 FLD 图分析和厚度变形图分析。 坯料形状、内孔半径和凹模圆角等工艺参数对变形结果的影响。本章小

18、结本章主要叙述了外缘翻边的基本原理和特点，基于此基础上提出本课题的研究，概述了当前外缘翻边的国内外现状确定了本文的主要研究内容。2 计算机模拟与 DYNAFORM 软件2.1 计算机模拟方法众所周知，传统的模具设计方法在很大程度上依赖于设计者的经验和技巧。模具设计费时、加工周期长，对设计中产生的问题只能靠试模后的修正，返工率高且花费较大。近年来，随着计算机硬件、软件技术的飞速发展和对材料成形过程物理规律研究的深入，材料成形过程计算机模拟技术取得了很大的进展。在冲模的设计与制造中， 引进 CA 工/CAM 技术，可以减少冲模设计中繁重的劳动，使复杂的计算能够自动化及优化，把数控加工和设计构形紧密

19、联为一体，缩短了制造周期。数值模拟即是通过数值计算得到用微分方程边值问题来描述的具体材料成形问 题中制件和模具的速度场(位移场)、应变场、应力场、温度场等等，据此预测制件组织性能的变化以及可能出现的缺陷。模拟式设计是对冲压变形过程进行描述，并在计算机图形终端上显示，以揭示金属变形规律，研究各工艺参数对成形过程的影响，利用计算机图形技术可将冲压成形过程、应力-应变分布、载荷一行程曲线等显示在屏幕上，直观地、动态地呈现在研究设计人员面前。这样，设计者在制模、组装、试模以前，能通过这个虚拟的材料加工过程检验制件的最终形状、尺寸、性能等是否符合设计要求，并可以在计算机上对一些工艺条件， 如变形速度、毛

20、坯尺寸、模具形状等进行反复调整与修改，也有助于正确选用设备和模具材料。正因为模拟方法的这些种种优点，如今以模拟式设计方法代替传统的经验式设计方法，已成为国际上成形工艺及模具设计发展的一个新方向。计算机辅助设计，主要用于降低设计劳动量和方便数控加工。二十世纪中后期， 国外一些公司研制的用于冲压成形过程分析、冲模设计、制造、分析的软件相继进入市场，如 ABAQUS, DYNAFORM, AUTOFORM 以及 Pro/E 等。采用这些软件的主要目的是: 优化设计过程，缩短模具设计和制造周期;通过各种力学模拟，提高冲压件质量等。最终目的是建立塑性变形理论基础上的分析设计系统。塑性成形过程模拟(CAE

21、)与塑性成形工艺及模具计算机辅助设计(CAD)技术的集成，与可视化技术和人工智能的结合，将形成一种方便高效的智能化设计和研究手段。它不仅能用于检验和优化设计， 也可用于探索新的塑形成性工艺和材料。采用计算机软件来模拟冲压变形过程，例如在计算机上首先根据经验准则设计冲压毛坯，然后通过计算机模拟在冲压模具上 “冲压”这个被假设的“冲压毛坯”。这个“试冲”的结果也可以在计算机图形终端上显示出来。如果模拟结果表明被选择的冲压毛坯形状尺寸不能满足冲压产品的设计要求，例如，翻边时所得零件高度达不到翻边件的设计要求:或者在成形过程中有冲压缺陷产生，那么就可以修改“冲压毛坯” 形状尺寸和采用必要的工艺措施，然

22、后重复在计算机上模拟及 “试冲气。计算机模拟现已广泛的应用于生产实际中，已被实践验证其可靠性、正确性，发挥着越来越重要的作用。2.2 数值模拟软件 DYNAFORM 的介绍2.2.1 DYNAFORM 软件简介DYNAFORM 软件是美国ETA 公司和LSTC 公司联合开发的用于板料成形数值模拟的专用软件，是 LS-DYNA 求解器与 ETA/FEMB 前后处理器的完美结合，是当今流行的板料成形与模具设计的 CAE 工具之一。在其前处理器 (Preprocessor) 上可以完成产品仿真模型的生成和输入文件的准备工作。求解器(LS-DYNA) 采用的是世界上最著名的通用显示动力为主、隐式为辅的

23、有限元分析程序，能够真实模拟板料成形中各种复杂问题。后处理器 (Postprocessor) 通过 CAD 技术生成形象的图形输出，可以直观的动态显示各种分析结果。 Dynaform 软件基于有限元方法建立 , 被用于模拟钣金成形工艺。 Dynaform 软件包含 BSE DFE Formability 三个大模块，几乎涵盖冲压模模面设计的所有要素，包括：定最佳冲压方向、坯料的设计、工艺补充面的设计、拉延筋的设计、凸凹模圆角设计、冲压速度的设置、压边力的设计、摩擦系数、切边线的求解、压力机吨位等。Dynaform 软件可应用于不同的领域，汽车、航空航天、家电、厨房卫生等行业。可以预测成形过程中

24、板料的裂纹、起皱、减薄、划痕、回弹、成形刚度、表面质量，评估板料的成形性能，从而为板成形工艺及模具设计提供帮助。2.2.2 DYNAFORM 的主要特色1集成操作环境，无需数据转换完备的前后处理功能，实现无文本编辑操作，所有操作在同一界面下进行2求解器采用业界著名、 功能最强的 LS-DYNA ，是动态非线性显示分析技术的创始和领导者，解决最复杂的金属成形问题。3工艺化的分析过程囊括影响冲压工艺的 60 余个因素以 DFE 为代表的多种工艺分析模块有好的工艺界面，易学易用4固化丰富的实际工程经验2.2.3 DYNAFORM 的功能介绍1. 基本模块DYNAFORM 提供了良好的与 CAD 软件

25、的 IGES 、VDADXF ，UG 和 CATIA 等接口, 以及与 NASTRAN, IDEAS, MOLDFLOW 等 CAE 软件的专用接口，以及方便的几何模型修补功能。IGES 模型转入自动消除各种孔 DYNAFORM 的模具网格自动划分与自动修补功能强大， 用最少的单元最大程度地逼近模具型面。比通常用于模具网格划分的时间减少了 99%！初始板料网格自动生成器，可以根据模具最小圆角尺寸自动确定最佳的板料网格尺寸，并尽量采用四边形单元，以确保计算的准确性。Quick Set-up, 能够帮助用户快速地完成分析模型的设置，大大提高了前处理的效率。与冲压工艺相对应的方便易用的流水线式的模拟

26、参数定义，包括模具自动定位、自动接触描述、压边力预测、模具加载描述、边界条件定义等等。用等效拉延筋代替实际的拉延筋，大大节省计算时间，并可以很方便地在有限元模型上修改拉延筋的尺寸及布置方式。多工步成形过程模拟 :网格自适应细分，可以在不显著增加计算时间的前提下提高计算精度。显、隐式无缝转换， eta/DYNAFORM 允许用户在求解不同的物理行为时在显、隐式求解器之间进行无缝转换，如在拉延过程中应用显式求解，在后续回弹分析当中则切换到隐式求解。三维动态等值线和云图显示应力应变、工件厚度变化、成形过程等，在成形极限图上动态显示各单元的成形情况，如起皱，拉裂等2. BSE(板料尺寸计算)模块采用一

27、步法求解器，可以方便地将产品展开，从而得到合理的落料尺寸。3. DFE（模面设计）模块DYNAFORM 的 DFE 模块可以从零件的几何形状进行模具设计，包括压料面与工艺补充。 DFE 模块中包含了一系列基于曲面的自动工具，如冲裁填补功能、冲压方向调整功能以及压料面与工艺补充生成功能等，可以帮助模具设计工程师 进行模具设计。基于几何曲面所有的功能都是基于 NURB 曲面的。所有的曲面都可以输出用于模具的最终设计。导角单元导角功能使用户对设计零件上的尖角根据用户指定的半径快速进行导角，以满足分析的要求。冲裁填补功能根据成形的需要，自动填补零件上不完整的形状。能在填补区同时生成网格与曲面。拉延深度

28、与负角检查图形显示零件的拉延深度与负角情况。冲压方向调整功能自动将零件从产品的设计坐标系调整到冲压的坐标系。压料面生成功能可以根据零件的形状自动生成四种压料面。生成的压料面可以根据用户的输入参数进行编辑与变形以满足设计要求。工艺补充面生成功能可以根据产品的大小、深度及材料生成一系列轮廓线。然后将这些轮廓线生成曲面并划分网格形成完整的工艺补充部分。还可以对生成的轮廓线进行交 MORPHINGDFE 模块中提供了线、曲面及网格的变形功能，可以很容易地处理POL、冲裁填补、工艺补充设计以及压料面设计。2.2.4 DYNAFORM 软件在板料成形过程中的分析流程在应用DYNAFORM 软件分析板料成型

29、过程时包括三个基本部分，即建立计算模型、求解和分析计算结果，其流程如图 2.1 所示。图 2.1 DYNAFORM 流程图直接在 DYNAFORM 的前处理器中建立模型，或在 CAD 软件（如 UG、CATIA 等） 中根据拟定或初定的成形方案，建立板料、对应的凸模和凹模的型面模型以及压边圈等模具零件的面模型，然后存为 IGES、STL 或 DXF 等文件格式，将上述模型数据导入DYNAFORM 系统。利用 DYNAFORM 软件提供的网格划分工具对板料、凸模、凹模、压边圈进行网格划分，检查并修正网格缺陷（包括单元法矢量、网格边界、负角、重叠结点和单元等）。定义板料、凸模、凹模和压边圈的属性，

30、以及相应的工艺参数（包括接触类型、摩擦系数、运动速度和压边力曲线） 。调整板料、凸模、凹模和压边圈之间的相互位置，观察凸模和凹模之间的相对运动，以确保模具动作的正确性。设置好分析计算参数，然后启动 LS-DYNA 求解。将求解结果读入 DYNAFORM 后处理器中，以云图、等值线和动画等形式显示数值模拟结果。分析模拟结果，通过反映的变化规律找到问题的所在。重新定义工具的形状、运动曲线，以及进一步设置毛坯尺寸，变化压边力的大小，调整工具移动速度和位移等，重新运算直至得到满意的结果。本章小结本章主要介绍了 DYNAFORM 的基本功能以及特点，了解了为什么要用 DYNAFORM 进行数值模拟分析及

31、它的分析过程。3 外缘翻边的数值模拟的基本步骤3.1 数据库操作导入模型启 动 eta/DYNAFORM,选择菜单 File-Import，将上面所建立的“*.igs”模型文件导入到数据库中，如图 3.1 所示。图3.1 导入模型编辑数据库中的零件层选择菜单 Parts-Edit，显示出编辑零件层的对话框。如图 3.2 所示。所有的已经定义的零件层都显示在列表中，零件层用名字和识别号标示出来。可以改变零件层的名字和标识号，同时，也可以从数据库中删除零件层。编辑修改各零件层的名称、编号（注意编号不能重复）和颜色，将毛坯层命名为 BLANK，将凸模命名为 PUNCH， 将压边圈命名为 BINDER

32、，凹模命名为 DIE，单击 OK 按钮确定。图3.2 编辑零件层的对话框图3.3Element工具栏工具网格划分设定当前零件层为 PUNCH 层，选择 Preprocesss 中的 Element 菜单项，弹出如图3.3 所示的 Element 工具栏。单击图中的 Surface Mesh 按钮，弹出如图 3.4 所示的Surface Mesh 对话框，设定 MAX.Size(最大单元值为 1)，其他各项采用默认值。单击 Select Surfaces 按钮，选择需要划分的曲面，分别如图3.5 和图 3.6 所示，最后得到的单元网格如图 3.7 所示。 网格检查为了防止划分所得到的网格存在一些

33、影响分析结果的潜在缺陷，需要对得到的网格单元进行检查。选择Preprocess 中 Model Check/Repair 选项，弹出如图 3.8 所示Model Check/Repair 对话框，最常用的检查为以下两项。在图 3.8 中，单击“自动翻转单元法向”（Auto Plate Normol）工具按钮，弹出如图 3.9 所示的 CONTROL KEYS 对话框，单击 CURSOR PICK PART 选项，选择工具面确定法线的方向，如图 3.10 所示。在图 3.8 中，单击“边界线显示”（Boundary Display）工具按钮，所得结果如图 3.11 所示。在观察边界线显示结果时，

34、为了更好地观察结果中存在的缺陷，可将曲线、曲面、单元和节点都不显示。.图 3.4 Surface Mesh 对话框图 3.5 选择划分网格的曲面图 3.6 选择PUNCh 层的曲面划分网格图 3.7 PUNCH 划分单元网格图 3.8 自动旋转单元法向检查图3.9单元的选取方式图3.11边界线显示项检查结果图3.10单元法向的选择3.2 传统设置定义工具选择 Tools Define tools 按钮.，弹出如图4.12所示的Define Tools对话框。Tool Name下拉表框中有Die、Punch、Binder三个选项。选择Binder选项，单击图4.12中得Add按钮，弹出如图4.1

35、3所示的Select Part对话框，选择需要定义的Binder零件层后，单击OK按钮确定，分别选择图4.12中得Die、Punch,依次定义Die、Punch工具。图3.12 Define Tools对话框图3.13Binder工具的指定定义毛坯设置工艺参数选择 Tools Define blank 按钮，弹出如图3.14所示的DiefineBlank对话框。单Add按钮添加BLANK零件层到Include Parts List列表框中。需要设定的参数有Material（材料）和Property(属性)两项。单击Material选项中的None 按钮，弹出如图3.15所示的Material对

36、话框。从Material Library中选择所需的材料类型CQ添加到Material列表框中作为模拟过程中成形材料的属性。单击Property选项的None按钮，弹出如图3.16所示的Property对话框，单击NEW按钮出现如图3.17所示的Proerty参数设定的BELYTSCHKOTSAY对话框。设定毛坯厚度值（UNIFORM THICKNESS）为1.0mm，其他参数采用系统默认值。最终毛坯工艺参数的设定值如图3.18所示。图3.14 Define Blank对话框图3.16 Property对话框图3.15 Material对话框图3.18 毛坯工艺参数设定自动定位工具图3.17

37、BELYTSCHKOTSAY对话框选择ToolsPositionAuto Position菜单项，弹出如图3.19所示的Auto Position Tools对话框。设定Master Tools为BLANK，Slave Tools为PUNCH、DIE和BINDER，接触间隙值设定为1.0mm。最后得到的各个工具的位置图如图3.20所示。图3.19Auto Position Tools图3.20 工具自动定位结果测量 PUNCH 的行程距离选择ToolsPosition tools Min.Distance菜单项，弹出如图3.21所示的Min.Dis between Tools 对话框。选择所要

38、测量的PUNCH和DIE两个工具，得到距离值为46.9mm， 这是下一步设定PUNCH运动曲线的设定依据。定义 PUNCH 运动曲线选择 ToolsDefine Tools 菜单项，弹出如图 3.22 所示的 Define Tools 对话框， 在 Tool Name 下拉列表中选择所要设定运动曲线的工具 PUNCH，单击 Define Contact 按钮，弹出如图 3.23 所示的 Tools Contact 对话框，对PUNCH 的接触参数进行设定， 此处采用系统的默认值。单击Define Load Curve 按钮，弹出如图 3.24 所示的 Tools Load Curve 对话框，

39、选择曲线类型为 Motion，单击 Auto 按钮弹出如图 3.25 所示的Motion Curve 对话框。 设定 PUNCH 的运动速度为 2000mm/s ， 运动位移为46.9-1.1=45.8mm。最终设定的 PUNCH 运动曲线如图 3.26 所示。图3.21Min.Dis between Tools图3.22Define Tools图3.23Tools Contact 对话框对话框对话框图3.24 tool load curve对话框图3.25 设置对话框图3.26 PUNCH的运动位移曲线定义压边圈（BINDER）的压力曲线在图 3.22 中下拉列表中选择 BINDER 选项，

40、同上所述设定 BINDER 的接触参数，采用系统默认值。单击该图中的 Define Load Curve 按钮，弹出如图 3.27 所示的 Tool Load Curve 对话框，选择曲线类型为 Force，单击 Auto 按钮弹出如图 3.28 所示的Force/Time Curve 对话框，设定 BINDER 的压边力为 20000N。最终设定 BINDER 的压力载荷曲线如图 3.29 所示。图3.27 Tool Load Curve 对话框图3.28 Force/Time Curve 对话框预览工具的运动图4.29 BINDER的压力载荷曲线3.3 设置分析参数及求解计算选择Tools

41、Animate菜单项，弹出如图3.30所示的Animate对话框。单击Play按钮可以观看工具的模拟运动。选择AnalysisLS-Dyna菜单项，弹出如图3.31所示的Analysis对话框。保存文件，再单击Analysis对话框的OK按钮。选择file submit dynafrom input deck 按钮，选择刚刚保存的文件打开，出现计算程序，如图3.32所示。图3.30 Animate对话框图3.31 Analysis对话框图3.32 计算程序3.4 后置处理绘制变形过程单击菜单栏中的 PostProcess 进入 DYNAFORM 后处理程序，即通过此接口转入到EtaPost 后

42、处理界面。选择 FileOpen 菜单项，浏览到保存结果文件的目录，选择正确的文件格式，然后选择 d3plot 文件，单击 Open 按钮读入结果文件。系统默认的绘制状态是绘制变形过程（Deformation），可在帧下拉列表中选择 All Frames 选项，然后单击“播放”按钮用来以动画的形式显示过程的变化，也可选择单帧对过程中的某步进行观察，如图 4.33 所示。最终得到的零件外形如图 3.34 所示。绘制厚度变化过程，成形极限图单击如图3.35所示的两个工具按钮Thickness和FLD，可绘制成形过程中得毛坯厚度的变化过程（如图3.36所示）和零件的成形极限图（如图3.37所示），同

43、样在“帧” 的下拉列表框中选择All Frames选项，然后单击play按钮以动画的形式显示过程的变化，也可以选择单帧对过程中的某步进行观察，根据计算数据分析成形结果是否满足工艺要求。图3.33 变形过程的绘制图3.35“成形过程控制” 工具按钮图3.34 零件的最终外形图图3.36 最终零件的壁厚分布情况图3.37 最终零件的FLD本章小结本章主要介绍了 DYNAFORM 进行外缘翻边模拟的基本步骤，如何得到最终的模拟结果图。熟悉后处理操作，利于进行数值模拟分析。4 外缘内曲翻边过程的数值模拟分析4.1 变形过程分析4.1.1 板料厚度变化分析在 YDNAFORM 的后置处理过程中，绘制了零

44、件的最终零件厚度变化图和最终零件的 FLD 图。而读懂这两幅图对外缘翻边的数值模拟分析至关重要。研究随着模拟进行厚度的变化过程。如图4.1-图4.4所示。第一阶段（第1步）：压边圈压住位置变薄，凸模还未接触板料，其他不变化。图 4.1 板料厚度变化过程（一）第二阶段（第10步）：凸模接触板料，挤压板料，被挤压部分受拉应力变薄，其他位置无明显变化。图 4.2 板料厚度变化过程（二）第三阶段（第11步）：边缘受拉应力大，变薄，底部被拉伸轻微变薄。图 4.3 板料厚度变化过程（三）第四阶段（第16步）：边缘位置得到金属补充破裂被补变厚，底部受拉伸轻微破裂。图 4.4 板料厚度变化过程（四）4.1.2

45、 板料 FLD 图分析由 FLD 图可以看到，随着模拟从开始到结束坯料形状不断变化，下面逐步观察和分析，如图 4.5-图 4.9 所示。第一阶段（第1步）：原始坯料。图4.5 板料FLD图变化过程（一）第二阶段（第3 步）：压边圈压住坯料，被压住的坯料部分已变形，其余部分因金属流动发生小的变形。图4.6 板料FLD图变化过程（二）第三阶段（第10步）：凸模接触坯料，接触部分发生大的变形，压边圈压住部分无变形，底部待变形。图4.7 板料FLD图变化过程（三）第四阶段（第11步）：凸模运动，边缘部分发生大变形，底部待变形，压边圈压住部分未变形。图4.8 板料FLD图变化过程（四）第五阶段（第17步

46、）：凸模运动到极限位置，底部发生大变形，边缘未变形图4.9 板料FLD图变化过程（五）4.2 各工艺参数对外缘翻边的影响分析4.2.1 板料内孔半径的影响为了分析板料内孔半径对外缘内曲翻边的影响，设计了如表 4-1 的若干组试验。仅改变板料半径，其他工艺参数不变序号板料中心角（度）板料内孔半径(mm)凹模圆角半径（mm）压边力(104N)418042205180622061808220工艺参数表 4-1（1）基于零件厚度变化分布图的翻边分析运用 DYNAFORM 软件对表 4-1 的各组参数进行数值模拟，得出各组试验的最终零件的壁厚分布图如图 4.10图 4.12 所示。图 4.10 试验 1

47、 的最终零件的壁厚分布图图 4.11 试验 2 的最终零件的壁厚分布图图 4.12 试验 3 的最终零件的壁厚分布图根据零件壁厚分布图分析不同板料内孔半径对外缘内曲翻边的影响，为了更好的分析，我们引入一个板料最大厚度变化量x（原板料厚度与翻边后最小厚度的差值）。由DYNAFORM模拟计算可得，板料内孔半径为4的x为0.268；板料内孔半径为6的x 为0.25；板料内孔半径为8的x为0.231。为了更直观的观察x随内孔半径的变化现制作如图4.13的坐标图。由图我们可以得出结论；在一定范围内，板料内孔半径越小， 板料厚度变化越明显，所以板料内孔半径越大越有利于翻边工艺。图 4.13 坐标图（2）基

48、于零件 FLD 图的翻边分析运用 DYNAFORM 软件模拟出 3 组试验的 FLD 图，如图 4.14-4.16 所示。图 4.14 试验 1 的零件最终FLD 图图 4.15 试验 2 的零件最终FLD 图图 4.16 试验 3 的零件最终FLD 图4.2.2 板料凹模圆角半径的影响由图4.14-4.16观察可知，试验1的法兰处发生了些许破裂；试验2的法兰处些许地方发生较严重起皱现象，底面发生轻微起皱；而试验3法兰处和底面都只发生了轻微的起皱。从而得出结论；在一定范围内，板料内孔半径越大翻边过程中越不容易破裂，板料内孔半径越小越容易破裂。工艺参数序号为了分析凹模圆角半径对外缘内曲翻边的影响

49、，设计了如表 4-2 的若干组试验。仅改变板料半径，其他工艺参数不变板料中心角（度）板料内孔半径(mm)凹模圆角半径（mm）压边力(104N)418061205180622061806表 4-2420（1）基于零件厚度变化分布图的翻边分析运用DYNAFORM软件对表4-2的各组参数进行数值模拟，得出各组试验的最终零件的壁厚分布图如图4.17图4.19所示。图 4.17 试验 4 的最终零件的壁厚分布图图 4.18 试验 5 的最终零件的壁厚分布图图4.19 试验6的最终零件的壁厚分布图由DYNAFORM模拟计算可得，凹模圆角为1的的x为0.258；凹模半径为2的x为0.25；凹模半径为4的x为

50、0.238，制作坐标图如图4.20。由图得出结论，在一定范围内，凹模圆角越大，板料厚度变化越明显，所以凹模圆角半径越大越有利于翻边工艺。（2）基于零件FLD图的翻边分析图4.20运用 DYNAFORM 软件模拟出 3 组试验的 FLD 图，如图 4.21-4.23 所示。图 4.21 试验 4 的零件最终FLD 图图 4.22 试验 5 的零件最终FLD 图图 4.23 试验 6 的零件最终FLD 图由图4.21-4.23观察得知，凹模圆角半径为1的板料要比凹模圆角半径为2和4的板料在法兰处和底面更容易破裂，起皱；而凹模圆角半径为2的板料又比凹模圆角半径为4的板料更容易起皱。所以得出结论；在一定范围内，凹模圆角半径越大，板料越不容易发生破裂，更有利于翻边工艺。4.2.3 板料形状的影响为了分析板料形状对翻边工艺的影响，故设计两组矩形板和圆形板的外缘内曲翻边试验，而其他工艺参数相同。运用DYNAFORM软件对两组试验进行数值模拟，得出两组组试验的最终零件图如图4.24，图4.25所示。图 4.24 矩形板外缘内曲翻边图 4.25 圆形板外缘内曲