塑性成形的刚塑性有限元方法概述.doc

塑性成形的刚塑性有限元方法概述
徐小波
摘要：总结了国外有关塑性成形的刚塑性有限元方法的研究现状以及刚塑性有限元法的概述和基本理论。指出三维成形有限元模拟在工业设计生产中具有广泛的应用前景。
关键词：有限元法塑性成形数值模拟
一、引言
21世纪的塑性加工产品向着轻量化、高强度、高精度、低消耗的方向发展。塑性精密成形技术对于提高产品精度、缩短产品交货期、减少或免除切削加工、降低成本、节省原材料、降低能耗，当前的生产的发展，除了要求锻件具有较高的精度外，更迫切地是要解决复杂形状地成形问题，同时还要不断提高锻件地质量、减少原料的消耗、提高模具寿命，促使降低锻件成本、提高产品的竞争能力。
塑性加工问题的研究方法主要有三种：理论解析研究方法、试验研究方法和数值模拟研究方法。这几种方法中，理论解析法突出的优点是求解直接，能给出力学量与参数间的函数全局关系，对揭示变形的力学本质和指导实践有重要意义。但这种方法只能求解简单的或经过简化的问题，对于复杂问题，求解复杂、难度大。
试验研究方法在理论解析与数学手段尚不完善的情况下，是一种不可缺少的研究方法；结果可靠，常作为理论解析与数值模拟的验证或对比数据；此外，试验研究可以发现新现象、新规律。然而，试验研究的局限性在于对复杂成形过程的研究有时试验手段与试验方法无法实现或难以达到要求；另外，耗资大、周期长、工作量大，为此，试验方法的应用存在严重的局限性，并且优化显得特别重要。
数值模拟的方法，克服了理论解析法求解复杂问题的困难，能减少试验工作量，近年来得到很大的发展。特别是基于变分原理的有限元法，由于其单元形状的多样性与方法本身的特点，原则上可以运用于分析任何材料模型、任意边界条件、任意形状的零件的塑性成形过程，得到广泛的应用。
二、刚塑性/刚粘塑性有限元法概述
塑性有限元法可以分为流动性塑性有限元（包括刚塑性有限元和刚粘塑性有限元法）和固体塑性有限元（包括小变形弹塑性有限元和大变形弹塑性有限元）两大类。
弹塑性有限元法是1967年由Marcal和King首先提出的，1968年山田嘉昭根据屈服准则的微分形和法向流动法则，推导出弹塑性应力一应变矩阵。弹塑性有限元的主要优点是考虑弹性变形和塑性变形的相互联系，不仅可以计算工件的变形、应力和应变的分布以及变形力等信息，而且可以有效的处理卸载问题，计算残余应力、残余应变和回弹。因此，弹塑性有限元宜于处理板料成形等问题。但是为了保证计算精度和解的收敛性，每一次加载步不能过大，以便只有很少的单元达到屈服状态。这种以小变形为理论基础的弹塑性有限元处理变形较大的塑性问题时，所需要的计算时间很长，并且随着变形的逐步增大会出现明显的误差。多年来国外学者一直致力于大变形的弹塑性有限元研究，Hibbit等在1970年首先提出了建立在有限元变形理论基础上的大变形有限元列式，Osias和McMeeking等在20世纪70年代中期分别采用Euler描述法建立了大变形有限元列式。此后大变形弹塑性有限元法不断发展，解决了塑性加工中的许多实际问题，但是仍有许多工作有待进一步完善。
针对弹塑性有限元法存在的问题，Lee和Kobayashi于1973年提出了刚塑性有限元法。该方法采用刚塑性材料模型，忽略了材料的弹性变形部分。虽然该方法也基于小应变的本构关系，但它不是像标准的弹塑性有限元法那样采用应力、应变增量形式来求解，因此，每一载荷步长可以取得大一些，从而减少计算时间，这样就克服了弹塑性有限元的不足。但是，刚塑性有限元法不能确定刚性区的应力、应变分布，也不能处理卸载问题。在锻造、挤压和拉拔的体积成形问题中，坯料的塑性变形量很大，弹性变形量相对比较小，同时回弹对锻件的精度影响不大。因此，用刚塑性有限元法模拟体积成形过程是比较适宜的。金属在高温或某些应变速率敏感的材料在常温条件下表现出的材料粘性，对材料塑性变形规律有较大影响，在有限元模拟中必须加以考虑。Zienkiewicz于1972年提出了粘塑性材料有限元法，进一步扩展了有限元方法在金属体积成形过程中的模拟作用。
三、刚塑性/刚粘塑性有限元法的基本理论
1．基本假设
在研究大变形金属塑性成形问题时，将金属坯料视为刚塑体，即把变形中的某些过程理想化，以便于数学上进行处理。并将坯料金属材料作如下假设：
（1）不考虑材料的弹性变形；
（2）材料均质且各向同性；
（3）变形过程中，材料体积保持不变；
（4）不考虑体积力；
（5）材料的变形服从Levy-Mises流动理论；
（6）加载条件给出刚性区与塑性区的界限。
2．塑性变形的边界值问题
在固体力学中塑性变形问题的求解可以看成一个边界值问题进行处理。刚塑性/刚粘塑性变形的边界值问题可以描述如下：
式中，SF为给定力面，为力面上的给定面力，SU为给定速度面，为速度面上的给定速度。
3．变分原理
为