考虑初始不均匀晶粒度的锻造过程数值模拟方法
【技术领域】
[0001]本发明涉及锻造过程数值模拟技术领域，尤其涉及一种考虑初始不均匀晶粒度的锻造过程数值模拟方法。
【背景技术】
[0002]目前，在大型锻件的生产过程中，首先要浇注铸坯，进而在压力机上锻造，利用计算机技术模拟微观组织在金属凝固和塑性变形过程中的演化过程能实现对锻件的成形质量进行预测。阶段的组织演变模拟方法都是基于均匀、统一晶粒度，即假设锻造前晶粒度从表面到芯部都是一致的。但是事实上，金属在凝固过程中会形成由表面到芯部的不均匀晶粒度，在后续的变形过程中，这种不均匀晶粒度会对最后的成形质量产生影响。由于在组织演变模拟方过程中假设了初始均匀晶粒度的锻件，这样就会造成模拟结果的不准确，从而与客观事实不相符。
【发明内容】
[0003](一 )要解决的技术问题
[0004]本发明的目的是提供一种考虑初始不均匀晶粒度的锻造过程数值模拟方法，以克服现有技术中，由于在组织演变模拟方过程中假设了初始均匀晶粒度的锻件，所造成模拟结果的不准确、与客观事实不相符等问题。
[0005]( 二)技术方案
[0006]为了解决上述技术问题，本发明提供了一种考虑初始不均匀晶粒度的锻造过程数值模拟方法，包括以下步骤:
[0007]S1.测量原始坯料的初始晶粒度分布；
[0008]S2.建立标准坯料的有限元模型；
[0009]S3.参照原始坯料的初始晶粒度分布，修改标准坯料的有限元模型的初始晶粒度分布;
[0010]S4.保存修改后的有限元模型。
[0011 ] 优选地，所述步骤SI中，通过procast软件的CAFE模块测得原始还料的初始晶粒度分布。
[0012]或者，所述步骤SI中，通过金相实验和image-Pro plus软件观察记录原始还料的初始晶粒度分布。
[0013]优选地，所述步骤S2中，通过deform软件建立标准还料的塑性成形有限元模型，并生成有限元模型的key文件，Key文件包括:单元的总数、每个单元的节点号、每个节点的空间坐标和每个单元的晶粒度初始值。
[0014]优选地，通过UG7.5软件建立标准坯料的几何模型，然后转换成stl文件并导入到deform软件中，进行网格划分，形成标准坯料的有限元模型。
[0015]优选地，所述步骤S3中，通过节点的空间坐标来确定节点的空间位置，再通过单元的节点位置来确定单元的空间位置，根据每个单元的空间位置的不同，改写每个单元的晶粒度值。
[0016]优选地，根据每个节点的空间坐标先给每个节点设定晶粒度值，再根据每个单元所包含的节点号，取加权平均值，得到该单元的晶粒度值。
[0017]优选地，每个单元包括四个节点。
[0018]优选地，通过编写fortran程序来修改标准还料的有限元模型的初始晶粒度分布，以实现标准坯料的有限元模型的初始不均匀晶粒度分布。
[0019]优选地，所述步骤S4中，保存修改后的key文件。
[0020](三)有益效果
[0021]本发明的考虑初始不均匀晶粒度的锻造过程数值模拟方法根据原始坯料的初始晶粒度分布修改标准坯料的有限元模型的初始晶粒度分布，实现了组织演变模拟过程中坯料的初始晶粒度的不均匀分布。其与现有技术相比具有以下优点:模拟精度高，应用该方法，能够实现初始不均匀晶粒度在deform的分布，在后续的微观组织模拟研究中，其模拟结果比初始均匀晶粒度的模拟结果模拟精度高；使用方便，与手工修改每个单元的晶粒度的值，节约了大量时间和精力，提高了建模的效率；开放性好，该方法的相关程序既能耦合到deform中去，又能独立运行。既能解决特定key文件的晶粒度赋值的要求，又能解决大批量key文件赋值要求；为解决锻造过程中混晶问题提供了思路和技术支持。
【附图说明】
[0022]图1为本发明实施例1中42CrMo棒材的三个截面的晶粒度分布图；
[0023]图2为本发明实施例1中42CrMo棒材的key文件中用户自定义单元变量；
[0024]图3为本发明实施例1中42CrMo棒材的用户自定义单元变量图形显示；
[0025]图4为本发明实施例1中42CrMo棒材的key文件中单元编号；
[0026]图5为本发明实施例1中42CrMo棒材的key文件中节点编号；
[0027]图6为本发明实施例1中42CrMo棒材的key中修改后的用户自定义变量；
[0028]图7为本发明实施例1中42CrMo棒材的修改后的用户自定义变量在deform中的显不；
[0029]图8为本发明实施例2中7050A1法兰的截面的晶粒度分布图；
[0030]图9为本发明实施例2中7050A1法兰的key文件中用户自定义单元变量；
[0031]图10为本发明实施例2中7050A1法兰的用户自定义单元变量图形显示；
[0032]图11为本发明实施例2中7050A1法兰的key文件中单元编号；
[0033]图12为本发明实施例2中7050A1法兰的key文件中节点编号；
[0034]图13为本发明实施例2中7050A1法兰的key中修改后的用户自定义变量；
[0035]图14为本发明实施例2中7050A1法兰的修改后的用户自定义变量在deform中的显示。
【具体实施方式】
[0036]下面结合附图和实施例对本发明的实施方式作进一步详细描述。以下实施例用于说明本发明，但不能用来限制本发明的范围。
[0037]本发明的考虑初始不均匀晶粒度的锻造过程数值模拟方法包括以下步骤:
[0038]S1.测量原始坯料的初始晶粒度分布；
[0039]可以通过procast软件的CAFE模块测得原始坯料的初始晶粒度分布，也可以通过金相实验和image-Pro plus软件观察记录原始还料的初始晶粒度分布。
[0040]S2.建立标准坯料的有限元模型；
[0041]通过UG7.5软件建立标准还料的几何模型，然后转换成stl文件并导入到deform软件中，进行网格划分，通过deform软件建立标准坯料的塑性成形有限元模型，并生成有限元模型的key文件，Key文件包括:单元的总数、每个单元的节点号、每个节点的空间坐标和每个单元的晶粒度初始值。
[0042]S3.参照原始坯料的初始晶粒度分布，修改标准坯料的有限元模型的初始晶粒度分布;
[0043]通过节点的空间坐标来确定节点的空间位置，再通过单元的节点位置来确定单元的空间位置，根据每个单元的空间位置的不同，改写每个单元的晶粒度值，其中，每个单元包括四个节点，根据每个节点的空间坐标先给每个节点设定晶粒度值，再根据每个单元所包含的节点号，取加权平均值，得到该单元的晶粒度值；通过编写fortran程序来修改标准坯料的有限元模型的初始晶粒度分布，以实现标准坯料的有限元模型的初始不均匀晶粒度分布O
[0044]S4.保存修改后的有限元模型；
[0045]即保存修改后的key文件。
[0046]实施例1
[0047]利用procast软件的CAFE模块实现了 42CrMo普通棒材的浇注数值模拟及微观组织模拟结果，得到三个不同截面的晶粒形貌结果，如图1所示。
[0048]利用UG7.5软件建立棒材的几何模型，直径40mm，高200mm，转换成stl文件，导入到deform软件中，进行网格划分，网格数为27437个。
[0049]在deform软件中，打开advanced下element data选择usr添加四个用户自定义变量=USREl为动态再结晶体积分数，其值为0，USRE2为晶粒度，其值为230 μ m，USRE3为动态再结晶晶粒尺寸，其值为0，USRE4为平均晶粒度尺寸，其值为230。
[0050]打开其key文件，如图2所示，USRELM是用户自定义单元变量，I为第一个单元，27437为单元总数，0.0000000E+000为数据格式，4为有四个单元变量，下一行I为第一个单元，后面四个数分别对应着四个用户单元变量的值分别为:USRE1，USRE2，USRE3，USRE4。在deform的图形显示为图3所示。
[0051]通过对key文件的分析，可知每个单元的节点号如图4，ELMCON为单元，27437为单元总数，4为单元节点数，可知是四面体单元，下一行中I为第一个单元，407，406，428，531为第一个单元的四个节点号，以此类推直到第27437个单元。
[0052]Key文件中也可以得到每个节点的空间坐标，如图5所示，图中RZ为节点，6278为节点总数，下一行I为节点号，后面三个数分别对应着这个节点的空间坐标。
[0053]通过key文件可知用户单元变量的数值及储存方式，key文件也清晰的给出了每个单元对应的节点标号，以及每个节点的三维坐标。通过单元的节点位置来确定单元的空间位置，在根据单元的空间位置不同，改写单元的用户自定义变量，比如初始晶粒度USRE2，从而完成初始不均匀晶粒度的分布。又由于key文件的储存严格按照fortran语言的要求，所以编写fortran程序更容易实现初始不均匀晶粒度的分布。
[0054]先简单区分中心区域和表面区域，棒料端面为直径40mm的圆，规定直径30mm为中心区，晶粒度为800 μ m，其余为表面区域，晶粒度为300 μ m。编写fortran程序，实现不均匀晶粒度在deform前处理中的设置。
[0055]在该程序的编写中，定义了整型变量二维数组INTEGER*4ELM(27437，4)，用来储存每个单元的节点号；实型变量二维数组REAL*8XYZ(6278，3)，用来储存每个节点的空间坐标；实型变量二维数组USRE (27437，4)，用来储存每个单元的4个用户定义变量。其编写主要语句如下:
[0056]OPEN(UNIT = 3，FILE = ’ H:\deforming\DEFORM_fforkpiece.KEY’)
[0057]该语句的作用是打开H盘下的key文件，设备号是3。
[0058]READ (3，30) (Al (I)，(XYZ (I, J), J = I, 3), I = I, 6278)
[0059]该语句的作用是把节点的空间坐标以30行号指定的格式储存在XYZ数组中。
[0060]READ (3, 40) (AI(I), (ELM (I, J), J = 1,4),1 = I, 27437)
[0061]该语句的作用是把单元的节点号以40行号指定的格式储存在ELM数组中。
[0062]READ (3，60) (Al (I)，(USRE (I, J), J = I, 4), I = I, 27437)
[0063]该语句的作用是把单元的4个用户单元变量以60行号指定的格式储存在USRE数组中。
[0064]WRITE (3，60) (I, (USRE (I, J), J = I, 4), I = I, 27437)
[0065]把已经改好的四个用户自定义变量存储到key文