模拟退huo算法的特点_剖析亚洲锻造会议感受（三）锻造工艺及模拟软件的应用...

1、概述   

金属成形过程就是一个由变化的温度场和微观组织场耦合的变形过程。这一过程可由一组微分方程来描述。这组微分方程包括：应力平衡方程，描述应变～位移／应变率～速度关系的几何方程，描述材料应力～应变、应变率、温度、微观组织关系的本构方程。描述微观组织变化与温度、应力、应变、应变率及其他类型微观组织变化关系的微观组织演化方程，以及变形体的一组力学和热学边界条件和包括初始微观组织在内的初始条件。

金属成形工艺过程的有限元模拟实质上就是在已知工件坯料几何形状、边界条件、初始条件及工件材料的所有一切参数的条件下用有限元方法求解这一组微分方程。通常以变形体的节点速度和温度为求解变量。考虑成形过程中的某一时刻、当变形体的速度场和温度场解出以后，通过积分可以得到变形体的位移场及变形体现时的各点坐标。据此由几何方程可进一步计算出变形体的应变率。应变：再用材料的本构方程由初始微观组织、温度、应变、应变率计算出应力；用微观组织的演化方程由初始微观组织、应变、应变率和应力计算出现时的微观组织变化。由边界的应力可以求得模具所受到的压力以及所需要的压机载荷。如果计算中将模具和锻件坯料都算作变形体，则模具的温度和变形可同时求得。如果在计算中加进去材料的破坏准则，在计算应力和应变时可以用破坏准则去判断现时的应力应变是否达到了破坏的程度以及发生何种破坏。

对于模锻，在合模后由工件的坐标和模具的位置可以知道是否有模具未充满和折叠缺陷。可见，锻造工艺有限元模拟分为两步：①用户输入要模拟的对象：工件模具的几何信息材料参数初始状态和边界条件，②模拟软件根据所输入的数据求解微分方程组，计算出所需要的各种物理量，并将这些计算结果输出给用户。这就是说锻造工艺模拟可以在不作任何试验的情况下就能使技术人员知道他所设计的工艺、模具和锻件坯料是否合理，如果不合理，他可以修改设计重新输人数据再模拟一次直到设计满意为止。后者就是锻造工艺CAD。可以看出，应用这项技术可以最大限度地减少试验次数，使工艺优化和新产品试制降低成本、缩短周期。     

当然，这一切是建立在模拟软件和计算机硬件均为理想的情况下。但是现时的情况并不理想：在硬件上，计算机的内存不够大，计算速度不够快。在软件上由于相关学科发展的限制，与具体材料有关的方程还不能完全精确地描述材料变化的实际过程。软件所设定的接触的边界条件也还不能完全精确地描述工件与模具之间实际发生的物理变化。

由于实际工况的复杂性及现实试验条件的限制，用户所输人的数据，特别是材料数据和接触边界数据，也会与实际情况有差别。另外为了适应现在计算机的速度和内存的限制，实际计算中引入了一些必要的简化。这一切都是造成模拟结果与实际锻造工艺有一定误差的原因。

当然，经过科研人员的多年的努力，在很多情况下这个误差在工业上是可以接受的，至少目前的模拟已经可以给出正确的趋势性判断。正因为如此，锻造工艺模拟及工艺CAD技术已被工业界所接受商业化软件的出现使这种技术的工业应用发展更加迅速。      

为了使这项技术更大地发挥作用，人们以提高模拟精度、计算速度和节省计算内存为目的的研究取得了很大的进展。

2、锻造工艺模拟技术的研究进展   

(1)关于材料的本构关系     

当有限元算法确定后，模拟软件中所使用的材料本构关系与实际模拟的材料的真实性能的差别大小就是影响模拟精度的关键因素。在金属成形模拟软件中通常使用传统的弹塑性或刚塑性本构关系。大多数软件都考虑了有限变形影响，这对实际变形很大的金属成形过程是十分必要的。

弹塑性模型计算精确，但由于要判断屈服和卸载加大了计算量。刚塑性模型相对计算简单但不能模拟回弹过程。在这些模型中通常只考虑应变硬化。至多考虑了动态回复，因此能满足冷锻和温锻工艺模拟的常规要求。但是对于那些成形过程中伴随有微观组织变化的情况(如热锻、超塑性成形)使用这种传统的本构关系便会产生很大的误差。例如热推锻中动态再结晶的发生将引起应力软化，因此使用传统的本构关系模拟热锻甚至会引起计算出的应力变化趋势错误。在这种情况下由于微观组织与宏观变形产生强烈的耦合。必须使用考虑微观组织与宏观变形耦合的本构关系，并要给出微观组织的演化方程。

近年来这方面的研究取得了明显的进展。例如，为了准确模拟热锻中的应力并预报其中的微观组织变化，德国Aachen大学的R.Kopp教授及其合作者将动态再结晶过程和晶粒长大过程用一组经验方程描述并插入到有限元程序中成功模拟了高温镦粗试样的晶粒度变化。长期以来，可用于金属热变形有限元模拟的考虑微观组织变化的本构理论研究始终没有进展。

90年代初，北京机电研究所的研究人员应用不可逆热力学的内变量理论和细观力学的方法根据材料科学给出的微观组织变化的机理分别建立了考虑三种变形机理和多种微观组织变化的超塑性本构关系以及专门针对热锻变形的考虑动态再结晶过程的热粘塑性本构关系。后者已被插入到有限元软件中，多个算例均得到计算与实验相符合的满意结果。目前进一步的研究工作仍在进行。

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

有了正确的本构方程之后，进一步的问题是如何通过试验测得这些包括在本构方程中的材料参数。对于传统的弹塑性/刚塑性本构关系，通常采用简单应力状态试验(拉伸，压缩，扭转等)来测试其中的材料参数。使用这种方法的基本要求是试样内应力、应变和温度均匀，不均匀性越大测出的材料参数误差越大。

在高温条件下当考虑微观组织变化时，要做到试样内微观组织完全均匀是很困难的，但是由于微观组织变化与宏观变形之间的非线形关系，微观组织空间分布的很小的差异会引起宏观应力应变的很大差别，同样宏观变形的不均匀性也会引起微观组织更大的不均匀。因为这种情况下试样的变形已不是简单应力状态，而应看作为一个复杂结构了。因此我们面对的问题是如何从一个复杂结构变形的试验结果反算出材料的本构参数。对干线性材料已有人建立了一种逆有很无法，但对于物理非统线和几何非线性很强的高温锻造过程，这种方法便无能为力了。

最近我国学者提出了一种试验与有限元相结合的分步迭代法。这种方法的收敛判据是有限元计算结果与试验结果之差小于某个小的常数。这种方法是首先用传统的方法由试验结果计算出材料参数的初始值。然后代入有限元程序中模拟试验过程，并根据计算结果与试验结果之差去修正这些参数。与普通的数值法代法相比，由于各参数对试验结果各物理量的影响关系十分复杂，因此修正公式要复杂得多。目前这种迭代方法的理论研究工作正在进行之中。      

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

与一般的结构分析相比，锻造工艺模拟的特点是，工件与模具的接触边界是随时间变化的。这种接触边界的处理和计算涉及到摩擦机理、接触与脱离搜索方法及判断准则、法向接触力计算方法等几个方面。

虽然目前这些问题已有了不少解决方法并已用于各种金属成形模拟的软件中，但是由于金属成形模具形状的复杂性，现有方法还有很多需要改进之处，所以至今接触问题算法仍然是当前金属成形模拟领域的研究热点之一。

例如在三维模拟软件中，通常模具表面的几何形状用很多平面网格逼近，这种方法虽然接触判断计算简单，但是因模具形状描述越精确要求网格越密，因此接触搜索所需要的计算机时也就越多。对于复杂形状的工件，工件的有限元网格和模具网格都很多时，在零件成形的后期，甚至出现处理接触边界的机时超过每一个时间步长所需总机时的一半。为此提出用参数曲面来逼近模具表面，例如MARC公司的AUTOFORGE软件使用的是B样条曲面。北京机电所正在进行的研究工作中用自然曲面和B样条曲面共同描述模具表面。这种描述方法，不仅减少了描述模具表面的网格数，缩减了接触搜索的时间，并且使法线连续变化消除了接触锁住现象。

在搜索技术方面，近年来也提出了如全局搜索、局域搜索等很多新的方法。摩擦接触力算法也有罚函数法、Lagrange乘子法等多种算法。在摩擦机理研究方面，很多人设计了专用装置对不同条件下的金属成形的摩擦规律进行了大量的实验研究，发现实际情况与库仑摩擦定律有明显差别并对库仑定律提出了分段描述的修正方法。      

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

有限元的自身特点决定了变形体网格的质量对计算精度影响很大。因此在金属成形模拟的整个过程中应保持网格质量不至于太差。但是由于锻件形状的多样性和复杂性，以及金属成形的大变形特征，研究工件初始网格生成和变形过程中对畸变过大的网格进行重新划分的方法就成了金属成形有限元模拟领域的另一个研究热点，尤其是对于三维体积成形模拟问题，这一问题的研究更显得重要，其难度也就更大。

当前大多数体积成形模拟软件都使用四面体单元，这种单元的特点是网格发生和重划算法比较简单，对复杂形状边界表面协调性好，但是计算精度低。因此使用六面体网格是人们的努力目标。为此国际上对此问题进行了多年的研究，但至今还是尚未成熟。目前存在多种六面体网格生成算法，如立方体填充法，由四面体网格到六面体网格的单元转换法，和由实体表面向实体内部逐层生成法。填充法能实现自动化且效率高，但边界上的单元质量差。单元转换法也可实现自动化，缺点也是网格质量差。第三种方法生成的网格质量最好，但是因为要与复杂边界已生成的网格协调，因此内部各网格的点与点、面与面关系非常复杂，实现这种方法难度很大。     

3、关于锻造工艺模拟软件的应用

由于锻造生产发展的需要和计算机技术水平的不断提高，锻造工艺模拟技术已经基本成熟并已经走向工业应用。应该指出，几乎所有的现行商业软件的预报能力都与实际的金属成形工艺存在一定的差距。面对在这种情况，注意对现有软件的缺陷进行补偿，使所得到的数值模拟结果充分发挥作用，是用户值得注意的问题。      

(1)重视输人数据的正确性      

用户所输入的原始数据对最终模拟结果影响很大。为了减少输入数据所引起的误差，用户必须十分重视这些基础数据的准确性。例如，为了获得准确的摩擦边界数据，应该针对所研究的项目进行专门的摩擦实验测量摩擦系数。为了给数值模拟中的开裂准则准备原始数据，应该进行材料成形性能实验。近期的研究表明材料的传热系数与材料所受的压力有关，因此这个参数应该通过加压试验测量出来。如上所述，对于考虑微观组织变化的热锻模拟，应该特别注意使用正确的测试方法。     

(2)实验补偿法      

应该准确了解模拟软件的功能，并据此实事求是地确定数值模拟的目的。对于软件无能为力的问题配以实验研究。例如，在模拟一个复杂的，特别是多道次的锻造工艺之前，应事先专门设计一个简化了的单道次工艺并对它用实验方法和数值方法同时进行模拟，以此确定模拟计算的误差并找出合适的输人参数。最后再用有限元模拟方法进行实际多道次工艺模拟。      

(3)理论分析补偿法      

这种方法要求用户根据所模拟的工艺过程、模拟软件的基本原理和主要计算方法、通过理论分析估算哪类数据可能产生误差，产生多大误差。以便在使用计算结果时扣除误差的影响。还可以分析出误差产生的原因并通过调整输人数据或使用用户子程序来消除计算误差。

常用软件

1、DEFORM

DEFORM作为基于有限元的工艺仿真系统,现在由美国SFTC公司推出井商品化,用于分析金属体积成形及其相关的各种成形工艺和热处理工艺,通过在计算机上模拟整个加工过程,帮助工程师和设计人员设计工具和产品的工艺流程,减少昂贵的现场试验成个:提高模具设计效率,降低生产和材料成本:缩短新产品的研究开发周期. DEFORM是专门针对金属成形模拟仿真的,其友好的图形用户界面便于工程设计人员可以很方便地进行前期的数据准备工作,这样使设计人员可以把精力集中在成形的模拟分析上,另外DEFORM的自动网格重划分优化技术特别适合分析大型复杂的成形工程问题。

DE-FORM具有如下突出特色:

(1) 模拟范围广: 材料流动,锻造负荷,模具应力,晶粒演变,缺陷成因等;

(2)适用工艺广:D冷、热、温过程;2锻造、拉拔、挤压、镦头、镦粗、轧制、摆辗、机加工车削、冲切等;3正火、退火、淬火、回火、时效、渗碳、蠕变、硬化等; 卦鼹霎瓤,猎关回得

(3) 适用设备多:液压机、锻锤、摩擦压力机、机械压力机、轧机、摆辗机等; 溪赐命天

(4) 材料库丰富:200余种专用材料数据库。

2、QFORM 

QFORM锻造模拟软件是基于莫斯科钢铁合金学院几十年的研究成果,2002年推出具有三维模拟功能的QFORM3D,2007年10月推出QFORM2D/3D 4.3 中文版。QFORM2D/3D 4.3中文版是一款专为锻造开发的有限元模拟软件,其功能目前只适用于锻造行业的大部分锻造工艺,目前在德国、俄罗斯,中国、瑞士等国使用较多.QFORM软件适合于以下情况的模拟:冷却和加热过程中的热弹性变形、多孔材料或粉末锻造的3D模拟、多模具模拟(多达20个,可在任何方向移动及旋转)、开坯锻造模拟、冷热模锻及切边、多工位辊锻制坯模拟、电镦、连续挤压,微观组织演变。

1.QFORM 的特点

(1)模具和工件表面三角单元网格生成是基于二次曲面逼近,内部四面体单元网格生成是基于表面三角单元拓展,这种二次曲面逼近生成三角单元的方法能够进行扭曲和弯曲,这样就可以高精确的配合模具表面。

(2)程序完全自动生成网格。网格密度由模具、工件形状,解决方案等其它参数确定·由程序自动监测生成。但使用者仍可以修改基本的网格生成控制参数,允许根据模拟精度要求修改网格生成的精密程度,改变有限单元节点总量。模拟时程序可以自动生成非均匀的有限单元,在临界区域自动生成更小的有限元网格以保证模拟精度。

(3) 全自动调节时间步距而实现了自控增量程序。

(4) 变形材料为刚-粘-塑(3D)性和弹-粘-塑(2D)性材料,模具材料为弹-塑性材料。

(5)工件材料流动应力取决于应变、应变率和温度,其它参数(工件及弹性材料的密度、比热容和热等特性,模具材料的典型强度)都与温度有关。

(6) 兼容四核处理器计算机,迭代解算器和直接解算器并行模式运行。 直接和迭代解算器允许在双CPU及四核计算机进行并行计算,当程序检测到多CPU后就自动运行,也可以根据需求对解算器进行人工选择。

(7)模拟过程中,随着模拟的进行,模拟结果会自动记录在案。

(8)程序含有设备数据库,润滑剂数据库和材料数据库,用户自己可对数据库进行修改和扩展。

3、FORGE

FORGE软件由法国CEMEF(材料成形研究中心)研究开发,CEMEF是由 Ecole desMines de Paris和ARMINES联合建立的研究中心.。FORGE软件基于有限元方法建立,被用于模拟热、温和冷锻金属成形工艺。FORGE软件可以进行3D模拟,也可以进行2D模拟(2D模拟应用于长零件成形,研究横断面的变形过程(平面变形);或者应用于轴对称零件成形,研究径向切面变形过程(轴对称变形)。

该软件模拟热锻成形时使用热一粘塑性模型,而模拟温锻和冷锻成形时使用热一弹一塑性模型,从而可以评估残余应力和成形后的几何形状变化。FORGE软件可应用于不同的领域,汽车、航空航天、手表制造、标准件等行业。稳定的自动网格生成和再生技术,在2D和3D模拟情况下可以模拟几何形状非常复杂的零件。

4、ABAQUS

ABAQUS 是一套功能强大的工程模拟有限元软件,其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。ABAQUS包括一个可模拟任意几何形状的单元库,并拥有各种类型的材料模型库。可以模拟典型工程材料的性能,其中包括金属、橡胶、高分子材料、复合材料、钢筋混凝土、可压缩超弹性泡沫材料以及土壤和岩石等地质材料。作为通用的模拟工具,ABAQUS除了能解决大量结构问题,还可以模拟其它工程领域的许多问题,例如热传导、质量扩散、热电耦合分析、声学分析、岩土力学分析(流体渗透/应力耦合分析)及压电介质分析。

一个完整的ABAQUS分析过程,通常由三个明确的步骤组成:前处理、模拟计算和后处理。 

在前处理阶段需定义物理问题的模型并生成一个ABAQUS输人文件。通常的做法是使用ABAQUS/CAE或其它前处理模块,在图形环境下生成模型。而一个简单问题也可直接用文件编辑器来生成 ABAQUS输入文件。模拟计算阶段用 ABAQUS/Standard 求解模型所定义的数值问题,它在正常情况下是作为后台进程处理的。

完成一个求解过程所需的时间可以从几秒到几天不等,这取决于所分析问题的复杂程度和计算机的运算能力。一旦完成了模拟计算得到位移、应变、应力或其它基本变量,就可以对计算结果进行分析评估,即后处理。通常,后处理是使用 ABAQUS/CAE 或其它后处理软件中的可视化模块在图形环境下交互式地进行,读人核心二进制输出数据库文件后,可视化模块有多种方法显示结果,包括彩色等值线图,彩色渲染图和z一y平面曲线图等。