板料成形静力隐式程序中的接触判断方法.pdfVIP

板料成形静力隐式程序中的接触判断方法 ’徐伟力杨玉英王仲仁+林忠钦 (}上海交通大学汽车工程研究所200030) (哈尔滨_T业大学金属塑性加工专业100050) 摘 要本文给出了一种简单、可靠的接触判断方法，把它应用在一个静力隐式有限元程序的开发中．并对一些 典型的板料成形过程和NUMIsHEET93国际会议的标准考题进行了计算，得到了较好的结果。 关键词静力隐zt板料成形有限元 t+ - ～ 一i ，∈。47 1 引言 近10多年来，板料成形过程的有限元数值模拟一直是世界范围内的研究热点。一方面， 它能增强汽车车身覆盖件模具设计和工艺设计的科学性，从而显著地缩短其制造周期：另一方 面，它能对新设计的车身进行成形性评价，并提出修改的意见，从而尽可能提高其成形性，减 少工序数目和模具套数，降低成形成本。这个车身覆盖件成形领域的“虚拟制造”技术正在美 国、德国、法国、意大利、日本等汽车工业发达国家得到推广和普及[】…。 现有的板料成形有限元程序根据时间积分算法可以分为：静力隐式、静力显式、动力显式 三种。这三种程序各有自己的优点和缺点。为了使它们能真正成为冲压模具设计和冲压工艺设 计的有力工具，迫切需要进一步改进它们的模拟精度和效率。 静力隐式有限元程序由于它的时间积分方案能真正满足板料成形的变形特点，因此它能得 到一个比较稳定的计算结果，并且它能比较有效地计算回弹问题。但它的致命缺陷是收敛性问 题和计算效率问题，国内外学者为此都做了不懈的努力o。7】。 接触判断方法是影响静力隐式有限元程序的主要因素，误判或漏判都将使收敛性和计算效 率下降。为此，本文对它进行了研究。 2接触判断方法 假殴所有的模具单元都是由三角形单元给出，并以任意一个板料节点f和任意一个模具单 元／的接触判断为例(见图1)。 圈l接触判断示例 柞：一个戴荷增量步内，节点f和单元／的运动方程可以分别_蚪j式(1)和式(2)来表示： x=x，+△驴f] y=y，+△V。口r} (1) z=t+△w妒rJ (2) 爿(工一x。一△“df)+曰(y—y。一△Vdf)+c(z—z。一△wdf)=0 式L卜l (x，．”，z，)和(x。，儿，：。)增量加载开始时节点弭口节点d的坐标； (△“p、△v口、△w，)和(△％、△”d、△wd)是节点f和单元／在载荷增量内的位穆增 茸：(爿、口、C)是模具的外法线矢量。 把式(1)代入式(2)中，得到式(3)： mf=仃 (3) 式中 Ⅲ=爿(△“p一△“d)+曰(△V∞一△Vd)+C(△w。一△wd) 门=爿(x，一x。)+B(y，一y。)+C(z，一z。) 式(3)的解町以归结为下面的情况： 若忉l≤e，，o，c和川≤eⅣDrc，那么f有无穷多解： 若l驯≤Prrorc和蚓≥Pr，orc，那么f无解： 若Iml2P，，Drc和H≥P”orc，那么f有唯一解，即：f=”／m。 这儿P，加，c是一个小的正数。 如果r有无穷多解，或者，无解，或者虽然有唯一解，但fO或者f1，则节点f不和单 元J接触，转向节点f和下一个单元／+1的接触判断。否则。求交点七的坐标(把，代入式(1) 中即可)。如果交点々满足式(4)，即它位于单元，的内部，则可以判断交点f接触单元，．否 则转向F一个接触判断。 ls△础6+s△。h+sA船c—s血kl≤PrrD， (4)