三维板料成形过程的显式有限元分析ξ.pdfVIP

第 13 卷第 3 期 计算结构力学及其应用 V o l. 13 N o. 3 1996 年 8 月 COM PU TA T IONAL STRU CTU RAL M ECHAN ICS AND A PPL ICA T ION S A ugust 1996 三维板料成形过程的显式有限元分析 李光耀 (湖南大学工程软件研究所, 长沙, 410082) 摘　要　本文采用基于随动坐标系的假定应变域壳单元及显式有限元格式求解三维板料成形问题。 板材料服从H ill 各向异性弹塑性准则, 板料与模具之间的接触界面由主仆面接触搜寻法处理, 接触 力由罚参数法计算。文中给出了几个三维成形过程的计算实例。数值算例表明, 本文方法具有较高的 计算精度和计算效率, 可在微机上分析中等复杂程度的成形过程。 关键词　壳体; 接触; 有限元 分类号　 242 21; 38 O T G 1　引　　言 板金属冲压成形是现代工业中一种十分重要的加工方法, 在汽车、家用电器、仪器仪表等 领域得到广泛应用。从力学的观点来看, 板料成形是一个包括大变形、大转动、材料非线性在内 的动态摩擦接触问题。由于板料在成形过程中经受非常复杂的塑性大变形, 在传统的模具设计 中, 一些成形缺陷如破裂、起皱、回弹只能靠不断的试模和修模来消除, 这就导致了很低的生产 效率, 并且要求设计人员具有丰富的生产和设计经验。 随着计算机和大规模计算技术的发展, 数值模拟技术已越来越广泛地用于模具设计。通过 对成形过程进行高精度的数值仿真, 设计人员可在计算机内观察设计参数对成形性能的影响, 〔1, 2, 3, 4, 5 〕 并可方便地调整设计参数直至得到满意的冲压计算结果 。近几年来, 国外特别是日本 和西德的各大汽车公司投入了大量的人力和财力研制三维板金成形过程的有限元分析软件。 在我国, 虽然几个汽车制造厂对该技术深感兴趣, 但研究过程进展非常缓慢, 更谈不上实际应 用。 本文采用基于随动坐标系的假定应变域壳单元及显式有限元格式求解三维板料成形问 题。接触界面由主仆面接触搜寻法处理, 板料服从H ill 各向异性弹塑性屈服准则, 模具做刚体 考虑, 摩擦接触力由库仑定律及罚参数法计算。数值算例表明, 本文方法能较准确地模拟三维 板料成形过程, 并且具有较高的计算效率, 可在微机上分析中小型成形件的成形过程。 2　基本公式 2 1　壳单元 本文采用 4 节点壳单元。对每一个单元, 定义局部坐标系如下: 国家自然科学基金资助。 　　　本文于 1995 年 12 月 25 日收到。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 254 　　　　　　　　　　　　　　　 计算结构力学及其应用　　　　　　　　　　　13 卷 e3 = n n , 　　n = r 13 ×r24 ( 1) e1 = gg , (2) g = r 12 - (r 12 e3 ) e3　　e2 = e3 ×e1 (3) 其中 = - , 为节点 的位置矢量, 如图 1 所示。该局部坐 r ij rj r i r i i 标系与单元一起运动, 所有的物理量都在局部坐标系中定义。 由等参插值, 得单元中的速度为