1、 西安航空职业学院 毕业论文 基于 PAM-STAMP 保护盖拉深工艺仿真研究 姓 名： 专 业： 航空电子 班 级： 完成日期： 指导教师 : 摘要 ： 拟采用仿真不同模具 R圆角半径、优化拉深次数、润滑边界条件的模拟工艺方案解决零件拉裂的问题，根据仿真得到最优工艺参数 ，采用试验件验证，并根据实际产品情况调整参数，最终形成紧急应对方案，弄清楚该零件的成形机理， 改进 某 批次产品出现的异常状况。 关键词 ： 保护盖 ，拉深，仿真 1 概述 保护盖（见图 1）零件材料为 LF3 铝板，壁厚 2mm，曲线边拉深成形各向异性指数大、易裂、易起皱。特别是出现了一次连续拉裂 6 件的情况。停工 分析
2、原因，并对拉深过程进行仿真模拟，弄清楚零件成型机理。否则反复使用调试料试模或者反复返修模具将消耗大量时间和人力物力。 实施目的与思路：拟采用不同 R 圆角的拉深模及单次多次拉深结合的工艺方案解决零件连续破裂问题，模拟仿 真设计出阳模、阴模、压边圈等边界条件，根据仿真结果得到新工艺方案，采用试验件验证，并根据实际产品情况调整展开参数，最终零件合格率超过 90%。 目前，有限元仿真软件经过几十年的发展和完善，已在国内外广泛应用于工程制造和科学研究领域，成为工程师分析计算的得力助手。因此采用有限元仿真软件模拟、计算，研究零件成形情况，为模具设计提供参考依据，制定工艺方案，可以缩短生产周期，降低制造成
3、本，是一种低成本的工艺研究模式。 图 1 J9445-02 保护盖零件图 2 模型建立 2.1 软件选用 PAM-STAMP 是专长于计算壳 单元的板料成形分析软件，能够分析保护盖的拉深过程，具有宏模块设置功能，设置简单， 软件能分析材料各向异性对拉深和回弹的影响。 2.2 几何模型 保护盖采用双动反拉深工艺，经过理论计算，拉深压边力需要约 10KN，考虑到冲模的装模高度，机械分厂具备此吨位以上的压力机有 40 吨、 63 吨带气垫曲柄机械压力机。 图 2 压力机的工作方式 根据压力机的工作方式（见图 2），拟定了双动拉深工艺方案：解决零件拉裂问题，需先进行仿真计算，观察并对比仿真结果，最后确
4、定工艺方案。 将坯料拉深成形作为仿真计算的焦点，因此建模只需创建坯料 和与其接触的模具，仿真时将压力机的运动与载荷赋予模具，忽略连接件以及压力机部件的作用，从而简化计算（图 3）。 用 UG 建模，所有模型的类型均为片体，保存为 igs 格式文件。 图 3 双动拉深的模型 2.3 材料物理模型 零件材料为 LF3 铝板，仿真计算需要坯料常温下的弹性模量、泊松比、密度、各向异性指数、应力应变曲线、成形极限曲线。我们与北京航空航天大学合作，通过试验，获取了以上性能参数。将 LF3 铝板性能参数输入 PAM-STAMP 软件，创建了 LF3 铝板的材料物理模型（表 1,图 4， 5， 6）。 表 1
5、 LF3 材料参数 密度 (kg/m3) 2.7 103 弹性模量 E(GPa) 70.97 泊松比 0.319 各项异性 r0=0.646， r45=0.927， r90=0.665 本构方程 21221133233222221148i l l 221 NGFHH 图 4 LF3 铝板性能参数设置表 图 5 LF3 铝板应力应变曲线 图 6 LF3 铝板成形极限曲线 2.4 边界条件 分析包括异形压边和拉深两个过程。 将已建立的 LF3 铝板的材料物理模型和料厚 2mm 的属性赋予坯料。双动拉深的坯料为平板，纤维方向为 +X。坯料为弹塑性变形体，屈服准 则为 Hill48，该准则在对复杂应力
6、状态下的坯料成形，可导出线性化的应力应变增量关系，适用于各向异性显著的材料，仿真过程依据此屈服准则建立本构方程。模具均设置为刚体。 拉深过程设置阳模行程为过程变量，将阴模固定不动，坯料与阴模接触的位置设为分析起始点，阳模与压板共同夹紧坯料下行至 30mm 距离设为分析终止点。 3 仿真过程难点分析与前处理设置 3.1 仿真过程难点分析 仿真过程计算的核心是坯料的变形，需要设置网格、运动、接触、载荷等，其中主要难点是坯料的网格设置。在仿真计算中，网格越密，则计算的结果就越精确，对 于拉深仿真，也是如此。模具在仿真计算中为刚体，不会发生变形，网格也保持不变，模具的网格与坯料的初始网格相匹配。要同时
7、考虑到网格单元数量对计算精度和计算效率的影响。本仿真控制初始网格数量 501000，最终网格数量 10000，控制计算时间 4 小时（见表 2）。 坯料网格过粗，单元数少，计算效率高，但计算误差较大，到了计算后期，还有可能会因为某些网格单元变形过大而导致溢出报错，无法得到计算结果；坯料网格过密，单元数多，计算精度高，会造成计算规模几何级数级的增加，大大增加计算时间和计算异常出现的概率，一旦出现错误 ，便无法得到计算结果。 坯料在拉深过程中，并非整个变形区的都一直在变形。变形区中有的区域在拉深前期完成变形，之后便不再变形，有的区域在拉深前期不变形，到了拉深后期才变形，将坯料在拉深过程中会变形的区
8、域网格划分较密，而对其他区域划分较粗，较前两者而言，这是一个软件优化办法（见图 7， 8）。在坯料成形仿真计算中，每步变形都会对每个网格单元进行计算。为了解决这个问题，采用了自适应网格优化方法，该方法适用于计算壳单元，在计算初期，使用较粗的单元网格，随着计算的逐步进行，在某些需要更细密网格描述的区域，将网格按 2n 倍数 进一步细分（见图 9），这样一来，坯料网格该粗的地方粗，该细的地方细，以尽可能小的计算耗费获得更好的计算结果。 表 2 坯料自适应网格优化单元数 初始网格单元数 最终成形后网格单元数 双动压边拉深 76 4216 图 7 坯料的初始网格 图 8 拉深坯料的不变形区与变形区 图
9、 9 坯料自适应网格优化 本项目对保护盖的压边拉深工艺方案进行仿真，坯料都应用自适应网格优化进行了前处理设置，在计算初期，坯料最大网格尺寸为 10mm，随着坯料逐渐变形，网格按角度准则和模具曲率半径准则进一步细分，以此作为仿真的计 算核心，运动、接触、载荷等参数均按照各工艺方案的拉深方式进行设置。 3.2 前处理设置 3.2.1 单动压边 该方案对 1 个工序进行仿真， 共计算 1 次压边（图 10）。 （ 1）网格划分 导入 igs 格式的模型文件，建模时模型类型为片体，导入到 PAM-STAMP 后转化为壳体，进行初始条件下的网格划分，定义网格类型为四边形，最大单元尺寸为 10mm，最小单
10、元尺寸为 0.1mm，最大弦差误差为 0.15mm，曲面相邻网格单元法向最大角度为 15，公差 0.01mm。模型的按以上方式划分网格后的单元数： 坯料： 76 阳模： 496 压边 圈： 945 （ 2）运动设置 阴模固定不动，阳模固定不动， 压边圈只允许沿 Z 轴移动，先移动至与坯料接触，与阳模夹紧坯料后保持压边力 10KN， 坯料不作运动设置。 a）起始 b）终止 图 10 单动压边的运动设置（ +Z=） （ 3）接触设置 所有模具与坯料的接触摩擦系数为 0.12。 （ 4）载荷设置 压边圈由液压垫施压提供压边力，沿 +Z 轴施加 10KN 力，阴模不施加载荷，阳模不施加载荷， 坯料受压
11、边圈的运动而夹紧，故不施加载荷。 （ 5）坯料自适应网格优化： 计算初期坯料最大网格尺寸为 10mm，随着坯料逐渐变形，网 格按角度准则和模具曲率半径准则进一步细分。 角度准则：当坯料上 2 个相邻单元间的夹角大于 10时，这两个单元将被细分。 模具曲率半径准则：当坯料上 1 个单元与临近模具面的曲率相差大于 2mm，这个单元将被细分。 （ 6）分析设置： 阴模行程设为过程变量，计算阴模从分析起始点到分析终止点期间坯料的变形。 （ 7）回弹分析 不计算回弹。 3.2.2 双动压边拉深 该方案对 1 个工序（压紧同时拉深成形）进行仿真， 共计算 1 次拉深（图 11）。 3.3 预成形工序 （
12、1）网格划分 网格类型为四边形，最大单元尺寸为 10mm，最小单元尺寸为 0.1mm，最大弦差误差为0.15mm，曲面相邻网格单元法向最大角度为 15，公差 0.01mm。模型的按以上方式划分网格后的单元数： 坯料： 682 阳模： 496 阴模： 945 压边圈： 2186 （ 2）运动设置 阴模固定不动 ， 阳模只允许沿 Z 轴移动，先移动至与坯料接触，然后与压边圈夹紧坯料继续下行（ -Z），下行距离为 60mm， 压边圈只允许沿 Z 轴移动，随阳模下行而下行 ， 坯料不作运动设置。 a）起始 b）进行中 c）终止 图 11 双动拉深的运动设置 （ +Z=） （ 3）接触设置 所有模具与坯
13、料的接触摩擦系数为 0.12。 （ 4）载荷设置 压边圈由液压垫施压提供压边力，沿 +Z 轴施加 10KN 力，阴模不施加载荷， 阳模向 -Z 轴施加运动速度 V，行程 30mm， 坯料受模具的运动而变形，故不施加载荷。 （ 5）坯料自适应网格 优化 计算初期坯料最大网格尺寸为 10mm，随着坯料逐渐变形，网格按角度准则和模具曲率半径准则进一步细分。 角度准则：当坯料上 2 个相邻单元间的夹角大于 10时，这两个单元将被细分。 模具曲率半径准则：当坯料上 1 个单元与临近模具面的曲率相差大于 2mm，这个 单元将被细分。 （ 6）分析设置 阳模行程设为过程变量，计算阳模从分析起始点到分析终止点期间坯料的变形。 （ 7）回弹分析 不计算回弹。 4 计算结果及分析 4.1 计算结果 4.1.1 单动压边 图 12 单动压边的壁厚云图 图 13 单动压边的应力分析云图 零件单动压边后，从壁厚云图（图 12）来看，零件厚度为 2.0mm，满足图样要求。从应力分析云图（图 13）来看，零件底部无残余应力，应力集中在 R 部位，塑性变形为压边部位。 4.1.2 双动压边拉深 1) 2)
