探究锻造工艺的模拟技术

　　1.关于材料的本构关系

　　当有限元算法确定后，模拟软件中所使用的材料本构关系与实际模拟的材料的真实性能的差别大小就是影响模拟精度的关键因素。在金属成形模拟软件中通常使用传统的弹塑性或刚塑性本构关系。大多数软件都考虑了有限变形影响，这对实际变形很大的金属成形过程是十分必要的。弹塑性模型计算精确，但由于要判断屈服和卸载加大了计算量。刚塑性模型相对计算简单但不能模拟回弹过程。在这些模型中通常只考虑应变硬化，至多考虑了动态回复，因此能满足冷锻和温锻工艺模拟的常规要求。但是对于那些成形过程中伴随有微观组织变化的情况(如热锻、超塑性成形)使用这种传统的本构关系便会产生很大的误差。例如热锻中动态再结晶的发生将引起应力软化，因此使用传统本构关系模拟热锻甚至会引起计算出的应力变化趋势错误。在这种情况下由于微观组织与宏观变形产生强烈的耦合，必须使用考虑微观组织与宏观变形耦合的本构关系，并要给出微观组织的演化方程。近年来这方面的研究取得了明显的进展。例如，为了准确模拟热锻中的应力并预报其中的微观组织变化，德国Aachen大学的R.Kopp教授及其合作者将动态再结晶过程和晶粒长大过程用一组经验方程描述并插入到有限元程序中成功模拟了高温镦粗试样的晶粒度变化。长期以来，可用于金属热变形有限元模拟的考虑微观组织变化的本构理论研究始终没有进展、90年代初，北京机电研究所的研究人员应用不可逆热力学的内变量理论和细观力学的方法根据材料科学给出的微观组织变化的机理分别建立了考虑三种变形机理和多种微观组织变化的超塑性本构关系以及专门针对热锻变形的考虑动态再结晶过程的热粘塑性本构关系。后者已被插入到有限元软件中，多个算例均得到计算与实验相符合的满意结果。目前进一步的研究工作仍在进行。

　　2.关于材料本构参数的测试方法

　　有了正确的本构方程之后，进一步的问题是如何通过试验测得这些包括在本构方程中的材料参数。对于传统的弹塑性/刚塑性本构关系，通常采用简单应力状态试验(拉伸，压缩，扭转等)来测试其中的材料参数。使用这种方法的基本要求是试样内应力、应变和温度均匀。不均匀性越大测出的材料参数误差越大。在高温条件下当考虑微观组织变化时，要做到试样内微观组织完全均匀是很困难的，但是由于微观组织变化与宏观变形之间的非线形关系，微观组织空间分布的很小的差异会引起宏观应力应变的很大差别，同样宏观变形的不均匀性也会引起微观组织更大的不均匀。因为这种情况下试样的变形已不是简单应力状态，而应看作为一个复杂结构了。因此我们面对的问题是如何从一个复杂结构变形的试验结果及算出材料的本构参数。对于线性材料已有人建立了一种逆有限元法，但对于物理非线性和几何非线性很强的高温锻造过程，这种方法便无能为力了。最近我国学者提出了一种试验与有限元相结合的分步迭代法。这种方法的收敛判据是有限元计算结果与试验结果之差小于某个小的常数。这种方法是首先用传统的方法由试验结果计算出材料参数的初始值，然后代入有限元程序中模拟试验过程，并根据计算结果与试验结果之差去修正这些参数。与普通的数值迭代法相比，由于各参数对试验结果各物理景的影响关系十分复杂，因此修正公式要复杂得多。目前这种迭代方法的理论研究工作正在进行之中。

　　3.关于接触边界的处理和计算方法

　　与一般的结构分析相比，锻造工艺模拟的特点是，工件与模具的接触边界是随时间变化的。这种接触边界的处理和计算涉及到摩擦机理、接触与脱离搜索方法及判断准则、法向接触力计算方法等几个方面。虽然目前这些问题已有了不少解决方法并已用于各种金属成形模拟的软件中，但是由于金属成形模具形状的复杂性，现有方法还有很多需要改进之处，所以至今接触问题算法仍然是当前金属成形模拟领域的研究热点之一。例如在三维模拟软件中，通常模具表面的几何形状用很多平面网格逼近，这种方法虽然接触判断计算简单，但是因模具形状描述越精确要求网格越密，因此接触搜索所需耍的计算机时也就越多。对于复杂形状的工件，工件的有限元网格和模具网格都很多时，在零件成形的后期，甚至出现处理接触边界的机时超过每一个时间步长所需总机时的一半。为此提出用参数曲面来逼近模具表面，例如MARC公司的AUTOFORGE软件使用的是B样条曲面。北京机电所正在进行的研究工作中用自然曲面和B样条曲面共同描述模具表面。这种描述方法，不仅减少了描述模具表面的网格数，缩减了接触搜索的时间，并且使法线连续变化消除了接触锁住现象。在搜索技术方面，近年来也提出了如全局搜索、局域搜索等很多新的力法。摩擦接触力算法也有罚函数法。Lagrange乘子法等多种算法。在摩擦机理研究方面，很多人设计了专用装置对不同条件下的金属成形的摩擦规律进行了大量的实验研究，发现实际情况与库仑摩擦定律有明显差别并对库仑定律提出了分段描述的修正方法。

　　4.关于网格生成和重划分算法

　　有限元的自身特点决定了变形体网格的质量对计算精度影响很大。因此在金属成形模拟的整个过程中应保持网格质量不至于太差。但是由于锻件形状的多样性和复杂性，以及金属成形的大变形特征，研究工件初始网格生成和变形过程中对畸变过大的网格进行重新划分的方法就成了金属成形有限元模拟领域的另一个研究热点，尤其是对于三维体积成形模拟问题，这一问题的研究更显得重要，其难度也就更大。当前大多数体积成形模拟软件都使用四面体单元，这种单元的特点是网格发生和重划算法比较简单，对复杂形状边界表面协调性好，但是计算精度低。因此使用六面体网格是人们的努力目标。为此国际上对此问题进行了多年的研究，但至今还是尚未成熟。目前存在多种六面体网格生成算法，如立方体填充法，由四面体网格到六面体网格的单元转换法，和由实体表面向实体内部逐层生成法。填充法能实现自动化且效率高，但边界上的单元质量差。单元转换法也可实现自动化，缺点也是网格质量差。第三种方法生成的网格质量最好，但是因为要与复杂边界已生成的网格协调，囚此内部各网格的点与点、面与面关系非常复杂，实现这种方法难度很大。