三维板料成形有限元数值模拟.pdfVIP

三维板料成形有限元数值模拟+ 张晓静孔永明周贤宾李敏 (北京航空航天大学机械工程及自动化学院板料冲压成形研究中心100083) 摘要本文讨论了板料成形动态显式有限元模拟的基本理论，介绍了成形模拟中接触摩擦基本算法。并对起皱、 破裂、回弹的模拟判断方法及研究现状进行了探讨。最后，结合NuMIsHEET’99三组标准考题的模拟研究．讨 论，有限元模拟技术在成形模拟中的应用。 ] 一l， I ’‘ 关键词板料成形有限元接触摩擦起皱破裂回弹 7 ’‘ 。， 1引言 板料冲压成形是现代工业中的一种重要的加工方法，在航空航天、汽车等各个领域均有广泛 的应用。然而如何保证板料塑性成形零件的成形质量，降低废品率，减少模具的返工或报废，缩 短模具设计与加工周期一直是板料塑性成形领域的一个难题。 板材成形，从力学角度而言，是一个同时包括几何非线性、材料非线性、边界非线性等问题 的非常复杂的力学过程。以往人们在分析这类问题时多采用基于塑性理论的解析方法。由于其复 杂性，人们不得不作出较多的简化和假设，这就使得理论分析的结果局限性较大，难以得出整个 变形过程的解答，特别是在计算回弹、起皱和残余应力时困难更大，迫使人们不得不突破解析法 的限制，寻求新的分析手段。随着计算机的应用和发展以及有限元方法的成熟，近年来发展了板 材成形数值模拟技术。该技术在减少甚至取消试模过程，缩短产品开发周期，降低产品开发成本 方面发挥着越来越重要的作用，已逐渐成为模具设计及优化的重要手段。 本文针对动态显式算法，论述了有限元数值模拟技术的基本理论和几个关键技术。并以法国 户s，公司开发的扳料成形有限元数值模拟专业软件——丹!^街zt^俨和美国觑66以K静b，DH＆ ∞地邶P”公司开发的通用非线性有限元分析软件—_觚叫Q己西为基础，对|Ⅳ明以衄r’99的三组 标准考题进行了模拟计算。 2基本理论 由连续介质力学理论，采用∞l咖应力，当前状态下的运动方程可以表示为 口¨+，一户“．一v“，=o (1) 其中，是体力，p是质量密度，v是阻尼系数。 将运动平衡方程和边界条件代入虚功原理基本方程，并引入G胁n公式，可以得到动力问题 的虚功方程。 ￡P蚱6虬艘+￡VqJ虬m2￡／：d％拯+￡i=，d坼打+￡霄J蚱dr一￡口。国，扪 … 然后将退化的埘珏珈1壳单元位移方程 ’航宁基金资助项目(96H51099)，长谷基金资助项目 (3) (其中，Ⅳ：为插值函数、u；为节点位移。) 及H圮￡正交各向异性屈服函数 ，2j志[F(口：：一口”)2+G(J。，一口n)2+H(口t，一口n)2 (4) +2￡口：，2+2腑，。2+2Ⅳ口．：2卜；孑2 j (其中，F、G、H、￡、M、Ⅳ分别是材料的各向异性参数，盯为材料的屈服应力) 代八虚功方程，得到动力问题的有限元方程如下 M子+cd：尸一F (5) 式中，Ⅳ为一致性质量矩阵，c为一致性阻尼矩阵，尸为外力向量，，为内力向量。 引入时间中心差分，并将质量矩阵M与阻尼矩阵c简化为对角阵，得到相互独立的方程组