筒形件拉深过程的数值模拟.pdfVIP

维普资讯 j 一 l 第 l8卷 第3期 太 原 重 型 机 械 学 院 学 报 Vo1．18 №3 1997年 9月 JOURNALOFTAIYUANHEAVYMACHINERY INSTITUTE Sep．1997 筒形件拉深过程的数值模拟 叮臼谬 太太(原原重重型型机机械学院压力加工系系．．太太原原030024)(深深圳圳黄黄金金灯饰集团，深圳518028)lP弓(J(‘ 摘 要 本文研究了板抖成彤模拟中减速系数选取和压边圈简化模型．嫡制了相应 的有限元数值模拟程序SPID。井模拟研究了筒彤件拉深过程中板厚变化规律及摩擦系 数的影响，计算结果与试验值比较吻合，表明该分析程序的可靠性。 关键词 筒形件；拉深；有限元法；数值模拟 中圈分类号 TG386 带压边的拉深成形过程，如图1所示。凸模底部与板料接触部分是以D1为直径的圆， 法兰变形区的内缘圆周半径为D2，形成了以D1一D2的环形塑性变形区域a．当凸模下行， 其底部首先与板料接触，并带动环形区域内的板料发生变形， 逐渐包容到凸模圆角上。变形的主要特征是板料厚度的减 薄。在法兰部分的板料开始进入凹模之前的变形阶段，称为 ＼j J f IJ S 拉深的初始阶段(图1中I所示)。随着凸模下行，拉深力增 加，法兰部分的板料被拉人凹模，形成拉深件的筒壁部分，为 ． 拉深的成形段(图1中Ⅱ所示)。 1 有限元模拟程序相关问题的技术处理 [1-] 本文的有限元模拟程序是用Fortran语言编制，在 In． tel_486微机上开发的，简称SPID。程序流程图略。 1．1 衰减因子口的选取 Ⅱ 由于塑性变形时的材料非线性效应，刚塑性有限元法基 于结点速度向量求解时，通常采用Newton—Rophson法线性 比并迭代求解，速度迭代修正公式为 图1 筒形件拉深变形过程 收稿 日期：1997—03—01 本文第一作者：男．1940年l2月生，压力加工系教授，摸具技术研究所所长．研究方向为冲压新拄 术蔑模具CAD 维普资讯 206 太 原 重 型 机 械 学 院 学 报 1997正 !一 ! !一 (1) 式中 ——迭代次数；y一1——第( 一1)迭代步的近似解；Av 第 迭代步的修正量； 卢一 减速因子，0 ≤l 减速因子 的取值直接影响到速度迭代求解的收敛性的效率，优其迭代计算 口的初值 。 更为重要。在拉深过程的初始阶段，毛坯与模具接触边界变化较为剧烈，有限元模拟时速 度解收敛性较差，因而采用 的选取与当前速度相对误差范数相关的方法。速度相对误差 范数定义为： = ll△ ll／ll 一1ll (2) 而 的选取可见表 l。 表 l 减速系数初值Po的选取 o．001lo．001~0．01lo．01--0．1lo．1～o．5 O．5 o．8 lo．5 Io．3 Io．1 O．O5 1．2 压边圈简化模型 压边I器产生的