一种基于晶粒形貌预测确定增材制造成形工艺参数的方法与流程

1.本发明属于增材制造技术领域，特别涉及一种基于晶粒形貌预测确定增材制造成形工艺参数的方法。


背景技术：

2.增材制造(又名3d打印)，是基于离散、堆积原理，根据三维模型采用材料逐层堆积的形式实现零件直接制造的技术。相比于传统的减材制造(如车削、铣削加工等)，该技术是一种原理性的变革，具有无需模具、快速响应、材料利用率高、可成形任意复杂构件等优点，其中，金属材料增材制造技术具有巨大的发展潜力和应用前景，是我国增材制造技术发展的重要战略方向，在航空航天、汽车及医疗等领域得到广泛应用。
3.根据所用热源的不同，金属增材制造主要分为激光、电子束和电弧增材制造技术。电弧增材制造是采用电弧作为热源，按照设定的轨迹熔化同步送进的金属丝材，实现金属零件的逐层制造。相比于激光和电子束增材制造技术，电弧增材制造具有成形效率高、成本低、尺寸大等技术优势，然而，电弧增材制造金属构件成形时存在较大熔池，成形零件通常具有粗大的晶粒组织，在合金凝固过程中，固
‑
液界面的生长形貌是影响合金组织性能以及缺陷形成的重要因素。因此，实现对增材制造零件晶粒形貌的准确预测有助于高性能金属结构件的制造。
4.金属增材制造中，零件的温度分布和冷却速度分布是极不均匀的，使得零件微观组织分布存在极大差异性。目前，获取增材制造零件晶粒形貌的方法主要包括物理实验法和数值模拟方法。
5.实验方法是对增材制造成形零件微观结构进行三维电子背散射衍射(ebsd)分析，需耗费大量时间和人力，且是基于成形零件的破坏性分析。材料微观组织数值模拟方法主要有相场模型和元胞自动机法和动力学蒙特卡洛法等。相场模型可以用于模拟材料凝固过程中枝晶形成过程、粗化现象以及固态相变等微观组织演化过程，在定量描述枝晶形貌方面具有一定优势。然而由于其计算量巨大，主要用来模拟求解二维问题，且仅能模拟较小区域尺度。元胞自动机法可以直接模拟枝晶形貌轮廓的生长过程，是应用较为广泛的材料微观组织介观尺度模拟方法，然而该方法是在已知瞬态温度场分布结果的基础上进行求解计算，同时模拟材料微观组织的求解网格较为精细，因此，仍需较长的计算时间。动力学蒙特卡洛法是一种基于晶界自由能最小原理的随机算法，需另外引入模型以建立模拟过程中的温度场和时间与真实温度信息和时间的关联。以上增材制造零件的微观组织模拟方法普遍存在计算量大、计算时间长等问题。


技术实现要素：

6.本发明所要解决的技术问题是针对现有技术的不足，提供一种基于温度场分布的增材制造零件晶粒形貌预测结果，确定增材制造成形工艺参数的方法。
7.为了实现上述目的，本发明采取的技术方案如下：
8.一种基于晶粒形貌预测确定增材制造成形工艺参数的方法，包括如下步骤：
9.(1)基于增材制造成形初始工艺参数，以及成形零件和基板的尺寸，建立有限元热模型；
10.(2)仿真增材制造过程，获得成形零件的瞬态温度场分布结果；
11.(3)基于步骤(2)零件的瞬态温度场分布结果，提取观测点瞬态温度变化结果；
12.(4)计算观测点的冷却速度、温度梯度和凝固速度；
13.(5)将步骤(4)计算得到的观测点的温度梯度和凝固速度值绘制到凝固图上，实现对该组工艺参数下成形零件晶粒形貌的预测；
14.(6)判断该组工艺参数下成形零件晶粒形貌是否满足设计要求，通过不断调整工艺参数，使得优化工艺参数下成形零件达到预期晶粒形貌，最终确定最佳工艺参数。
15.具体地，步骤(1)中，有限元热模型为增材制造零件的三维几何模型，包括基板和成形零件两部分，所述基板为三维块体结构，所述零件为单道多层或多道多层自下而上依次堆叠而成的结构件；对三维几何模型进行网络划分，模型采用全六面体网格进行划分，基板上靠近零件区域划分较为致密的网格，远离零件区域划分较疏网格；确定基板和成形零件的热物理性能参数及成形工艺参数，计算所需的热物理性能参数包括密度、熔化潜热、固相温度和液相温度、以及随温度变化的热传导率和比热等；设置热源模型参数，并且设置初始温度条件和换热边界条件。
16.具体地，步骤(2)中，增材制造实验中，高能束按照设定的路径持续熔化送入的金属丝材来实现零件的逐点成形。仿真增材制造过程采用“生死单元”方法仿真增材制造过程，按照高能束移动路径有序地激活一组离散单元，从而实现金属材料的逐点熔覆；在完成以上模型参数设置后，通过求解热分析控制方程，获得成形零件的瞬态温度场分布结果。
17.具体地，步骤(3)中，分析成形零件的温度场分布结果，选取成形零件内部关键位置作为观测点，提取观测点的瞬态温度变化结果。所述的观测点选取成形零件内部关键位置，提取观测点的瞬态温度变化结果；所述的关键位置为一处以上，包括零件体心、零件上表面中心点等位置。
18.具体地，步骤(4)中，基于观测点的瞬态温度变化结果，获得观测点处金属材料的凝固冷却速度、温度梯度和凝固速度；
19.其中，金属材料的凝固冷却速度由下式定义：
[0020][0021]
其中，t
s
和t
l
分别表示金属材料的固相温度和液相温度；t
s
和t
l
分别表示金属材料达到固相温度和液相温度的时间；当金属材料经历多个热循环时，选取高于固相线温度的最后一条热循环曲线数据计算金属材料的凝固冷却速度；
[0022]
成形件在t＝t
l
时的温度梯度g由fourier定律得出：
[0023][0024]
式中q表示热流矢量，即单位时间通过单位面积的热能，该值从有限元模拟结果中计算得出；k为液相温度t
l
下金属材料的热导率；凝固速度r由凝固时的冷却速率和温度梯
度计算得出：
[0025][0026]
具体地，步骤(5)中，所述金属材料的凝固图通过查阅该种金属材料晶粒形貌从柱状晶向等轴晶转变的判据公式计算绘制得出，为本领域公知方法。
[0027]
具体地，步骤(5)中，根据观测点处计算得到的温度梯度和凝固速度值，将其绘制到金属凝固图上，根据该点所在凝固图上的微观组织区域，实现对零件晶粒形貌预测。
[0028]
具体地，步骤(6)中，若初始工艺参数下成形的零件达到预期晶粒形貌，则采用该组工艺参数，完成零件的成形；否则，调整零件成形时的功率和扫描速度等工艺参数，返回步骤(1)重新进行建模以及晶粒形貌预测工作，直至达到预期晶粒形貌，最终确定最佳工艺参数。
[0029]
有益效果：
[0030]
本发明基于有限元模型技术，采用“生死单元”方法仿真零件增材制造过程，获得成形零件的温度场分布结果，通过分析观测点的瞬态温度变化结果，计算得出观测点处金属材料的凝固冷却速度、温度梯度和凝固速度，在对比分析金属凝固图后实现对零件晶粒形貌的预测，实现零件成形时工艺参数优化，获得预期晶粒形貌的零件。该方法相比于其他增材制造晶粒形貌模拟方法计算量较小，更为快捷方便，从而对实现增材制造工艺参数优化，细化成形零件晶粒具有重要意义。
附图说明
[0031]
下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述和/或其他方面的优点将会变得更加清楚。
[0032]
图1是本发明的流程示意图。
[0033]
图2是本发明实施例中增材制造三道三层零件的示意图。
[0034]
图3是本发明实施例中增材制造三道三层零件的全六面体网格划分模型示意图。
[0035]
图4是金属材料熔覆示意图。
[0036]
图5是三道三层成形零件a点处的温度随时间变化曲线。
[0037]
图6是ti
‑
6al
‑
4v钛合金凝固图。
具体实施方式
[0038]
根据下述实施例，可以更好地理解本发明。
[0039]
需要说明的是，本发明提出的增材制造零件晶粒形貌预测方法适用于多种高能束金属增材制造技术，包括激光、电子束或电弧增材制造技术，本实施例以模拟电弧增材制造技术为例。
[0040]
晶粒形貌预测方法可以适用于模拟多种金属材料，包括钛合金、镍基合金、不锈钢、铝合金、镁合金、铁基合金或高熵合金等。在本实施例中，基板材料为ti
‑
6al
‑
4v钛合金，成形零件材料为ti
‑
6al
‑
4v钛合金，利用abaqus有限元软件模拟电弧增材制造过程。
[0041]
本实施例提供的多层多道电弧增材制造金属零件晶粒形貌预测方法主要步骤如图1所示，包括以下几个部分：
[0042]
步骤一，建立有限元热模型：
[0043]
建立增材制造零件的三维几何模型，如图2所示，包括基板和成形零件两部分，所述基板为ti
‑
6al
‑
4v钛合金板，尺寸为长200mm，宽100mm，高5mm，所述成形零件为三道三层依次堆叠而成的ti
‑
6al
‑
4v钛合金零件，每条熔道长80mm。
[0044]
对三维几何模型进行网络划分，模型采用八节点六面体线性热传导单元(dc3d8)，由于零件及其周围区域在增材制造过程中存在较高的温度梯度和冷却速度，为兼顾模型精度和计算效率，采用渐变网格，基板上靠近零件区域划分较为致密的网格，远离零件区域划分较疏网格，如图3所示。确定基板和成形零件的热物理性能参数，包括密度、熔化潜热、固相温度(1878k)和液相温度(1928k)、以及随温度变化的热传导率和比热等。同时确定零件成形工艺参数，如表1所示；选取双椭球热源模型作为模拟电弧增材制造过程的热源模型，同时为基板设置初始温度25℃以及对流和辐射边界条件。
[0045]
表1零件成形工艺参数
[0046][0047]
步骤二，仿真增材制造过程：
[0048]
为实现增材成形时金属材料的逐点增加过程，模拟中采用“生死单元”方法仿真增材制造过程，即按照高能束移动路径有序地激活一组离散单元，从而实现金属的逐点熔覆。图4为金属增材制造过程的模拟示意图，底部的灰色单元为基板，顶层的透明单元表示被“杀死”的单元，蓝色单元为已沉积金属。随着热源的逐步移动，热源所到之处的对应单元逐渐被“激活”，如图中红色单元(set
‑
1)表示电弧正在熔覆该处金属。完成以上模型参数设置后，可生成用于数值计算的热分析模型文件，并采用abaqus命令提示窗口执行计算，通过求解热分析控制方程，获得成形零件的瞬态温度场分布结果和热流矢量分布结果。
[0049]
步骤三，提取观测点瞬态温度变化结果：
[0050]
基于零件的瞬态温度场分布结果，选取第一层中间熔道末端位置作为观测点，如图3中a点。分析该点的温度随时间变化曲线，如图5所示，该点经历了八次热循环，当沉积第一层熔道，该点经历了第1～2次峰值温度；沉积第二层熔道时，该点经历了第3～5次峰值温度；沉积第三层熔道时，该点经历了第6～8次峰值温度。第一次峰值温度高于熔点且为八次峰值温度中的最高点，表明当热源经过该点时，金属完全熔化为液态；当第四个峰值温度来临时，热源移动到第二层中间熔道，该点的温度超过固相线温度，金属重熔，随后的峰值温度均小于固相线温度，表明后续熔道沉积时，该点只经历后热作用。
[0051]
步骤四，计算观测点的冷却速度、温度梯度和凝固速度：
[0052]
通过分析观测点a点的温度随时间变化曲线，选取第四个峰值温度下降沿区段对a点的冷却速度、温度梯度和凝固速度进行计算。为更清晰地解释a点冷却速度、温度梯度和凝固速度的计算过程，将第四个峰值温度曲线单独提炼出来进行解释说明，如图5b所示。ti
‑
6al
‑
4v钛合金的固相温度t
s
为1878k，液相温度t
l
为1928k，a点在经历完第四个峰值温度后冷却到液相温度t
l
和固相温度t
s
所用的时间分别为480.252s(t
l
)和480.422s(t
s
)。因此，观测点a点金属材料的凝固冷却速率为：
[0053][0054]
基于步骤二中热模型求解获得的热流矢量分布结果，提取观测点a点在t＝t
l
时的热流矢量q，该值为10081510w/m2，因此，a点金属材料在t＝t
l
时的温度梯度g为：
[0055][0056]
k为液相温度t
l
下金属材料的热导率，对于ti
‑
6al
‑
4v钛合金而言，该值为34w m
‑1k
‑1。根据观测点a点凝固时的冷却速率和温度梯度结果，可以求解出凝固速度r为：
[0057][0058]
步骤五，将观测点的温度梯度和凝固速度值绘制到凝固图上，预测晶粒形貌：
[0059]
根据步骤四中计算结果，可知观测点a点的温度梯度g为2965k/cm，凝固速度r为0.08226cm/s，将该点的温度梯度和凝固速度值绘制到ti
‑
6al
‑
4v钛合金的凝固图上，如图6所示，可以发现该点位于柱状晶区域，表明在表1所示的工艺参数下，电弧增材制造成形零件在观测点a点处的晶粒形貌为柱状晶。通过以上温度场模拟、参数计算和凝固图绘制过程，实现对电弧增材制造零件晶粒形貌的预测。
[0060]
步骤六，判断该组工艺参数下成形零件晶粒形貌是否满足设计要求，完成零件成形实验：
[0061]
根据步骤五中ti
‑
6al
‑
4v钛合金的凝固图绘制结果，可知表1所示工艺参数下成形的零件具有柱状晶形貌，满足柱状晶粒的预期晶粒形貌要求，因此，采用表1所示工艺参数完成ti
‑
6al
‑
4v钛合金零件的成形。
[0062]
本发明提供了一种基于晶粒形貌预测确定增材制造成形工艺参数的方法的思路及方法，具体实现该技术方案的方法和途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。