增材制造过程的热—力学模拟

增材制造过程的热—力学模拟

Abaqus/Standard为热机械模拟提供了一个通用框架，允许您定义适当的边界条件、载荷、相互作用、约束和温度相关的材料属性。本节概述了可用于但不限于添加制造工艺热力学模拟的特殊技术。

增材制造工艺的热应力分析:

• 允许在时间和空间上精确指定加工条件，如材料沉积顺序、激光扫描路径、激光特性和加工环境；
• 允许精确控制不同模拟需求的解决方案的保真度；
• 随后可以分析热处理应用、支座移除和/或机械性能测试等；
• 并且通常*比本征应变分析更精确*，但是计算成本更高。

本节讨论：

• 增材制造过程的热—力模拟。
• 热力学模拟中的渐进元素激活过程。
• 通过对流和辐射进行渐进式冷却。
• 通过移动热源进行渐进式加热。
• 控制模拟的规模和解决方案的保真度。
• 解决收敛困难
• 输入文件模板

增材制造过程的热—力模拟

增材制造过程的连续热应力分析包括由该过程在印刷零件上引入的热负荷的瞬态热传递分析，随后是由热分析的温度场驱动的静态结构分析。该分析还可以包括支撑结构(如果构建需要)和基底，在基底上构建零件和支撑，以考虑它们对热传导、零件变形和残余应力的影响。

对于热传递分析，必须在模拟中包括制造过程的以下方面：

• 渐进式材料沉积：每个增材制造过程随着时间增加新材料。
• 沉积材料的逐步加热：对于一些增材制造工艺，新沉积的材料层被加热到熔点，导致材料熔合到下面的层或基底上。
• 打印部件的渐进冷却：在打印部件时，其冷却表面会随着时间不断变化。

对于应力分析，

• 来自传热分析的温度场用于驱动应力分析。
• 可以使用与传热分析类似的渐进材料沉积技术。
• 与温度相关的材料属性可用于获得精确的应力结果。

热力学模拟中的渐进元素激活过程

在Abaqus/Standard中，可以使用渐进元素激活来模拟增材制造过程中的渐进材料沉积(参见渐进元素激活)。部件和支撑结构可以被逐步激活。在分析过程中，任何被定义为渐进激活的元素都被材料填充(完全或部分)或保持为空(非激活)。在使用用户子程序UEPACTIVATIONVOL进行分析的过程中，可以控制每个增量中每个元素的激活和材料体积分数。

材料沉积序列在事件序列中定义(例如，以j挤出机事件序列的时间和空间坐标的形式；更多细节见Event series)，并使用用户子例程UEPACTIVATIONVOL进行处理，以定义在给定时间增量内向元素添加了多少材料。为了简化此过程，您可以从UEPACTIVATIONVOL中访问刀具路径—网格相交实用程序，以计算每个元素中材料的体积分数的变化和时间增量，同时考虑有限元网格与指定的材料沉积事件系列数据的相交(请参阅刀具路径-网格相交模块)。或者，可以更新材料方向，使其与轨迹对齐。

此外，Abaqus/Standard为一些最常见的增材制造工艺提供了许多简化的解决方案，不需要您编写用户子程序。这些技术在增材制造专用技术中有所描述。

初始温度

在零件的增材制造过程中，新沉积的材料在给定的初始温度下进入。您可以通过为要打印的零件的网格节点指定初始条件来定义该初始温度值(请参见定义初始温度)。ABAQUS将初始温度应用于这些元素内部的材料，在这些元素第一次被激活的地方。

使用相对粗糙的时间增量和网格的零件级热机械仿真的应力分析是上述陈述的例外。在这种情况下，您应该指定一个代表松弛温度的初始温度，超过该温度时，热应变引起的热应力可以忽略不计。该初始温度通常不同于在制造过程中沉积的新材料的初始温度。有关详细信息，请参见控制模拟的比例和解决方案保真度。

应力分析中的位移输出

在Abaqus/Standard中使用渐进式元素激活时，可以控制非活动元素的行为以跟随或不跟随模型中激活元素的变形(请参见控制非活动元素的行为)。除了位移和旋转(U、UT和UR)之外，这两种行为预计在小变形的限度内会产生相似的结果。

跟随变形的非活动元素，也称为*“安静”元素，总是出现在模型中并参与求解，但它对整体响应的贡献可以忽略不计。在这种情况下，附着到非活动元素的节点在其任何附着的元素被激活之前可能会经历非零位移。节点输出变量U*、UT和UR代表**从分析开始测量的位移和旋转**，包含节点非激活和激活期间的位移贡献。Abaqus还提供了节点输出变量UACT、UTACT和URACT，这些变量对应于从连接到节点的单元第一次被激活开始测量的位移和旋转。

不跟随变形的非活动元素不会影响模型的刚度，也不会参与求解。任何附加到非活动元素的节点都保持在它们的初始位置。在这种情况下，节点输出变量UACT、UTACT和URACT分别与输出变量U、UT和UR相同。

不管为非活动元素选择的行为如何，激活时元素的配置通常不同于原始配置，因为活动元素和非活动元素共享的节点会发生位移(请参见初始配置)。激活元素时，激活时的配置是后续元素计算的参考。因此，输出变量E代表从单元被激活开始测量的应变。

时间增量

在包含渐进元素激活的分析中使用的时间增量可能会影响最终结果。假设执行两个分析来模拟材料珠的材料挤出过程，该材料珠被建模为一行元素，在两个不同的时间增量下：小的时间增量每个增量激活一个元素，而大的时间增量每个增量激活两个元素。在两种应力分析中，每隔一个单元的初始配置是不同的，导致残余应力和变形的结果不同。通过执行时间步进收敛研究，可以在热力学分析中选择适当的时间增量。

通过对流和辐射进行渐进式冷却

在增材制造过程中，随着新材料的沉积，先前暴露的材料表面被覆盖，随着时间的推移，新的自由表面被创建。您可以在连续变化的自由表面上定义表面对流和辐射(请参见在Abaqus/Standard中的元素演化面上指定基于元素的膜条件和在Abaqus/Standard中的元素演化面上指定基于元素的辐射条件)。Abaqus/Standard持续跟踪在构建过程中任何给定时间反映打印零件当前形状的演变自由表面，并仅对这些表面应用膜和辐射载荷。

通过移动热源进行渐进式加热

许多类型的增材制造工艺涉及移动热源(激光、电子束、焊枪等)的应用来熔合原材料。您可以在用户子程序UMDFLUX中定义一个或多个非均匀分布的热源(请参见在用户子程序UMDFLUX中定义运动或静止的非均匀热流)。

移动热源的扫描路径在Event series中定义(例如，以热源的时间、空间坐标和功率的形式；更多细节见event series)，并使用用户子程序UMDFLUX进行处理，以定义在给定时间增量内在给定单元上发生的每个加热事件。为了简化此过程，您可以从UMDFLUX中访问刀具路径—网格相交实用程序，以计算每个单元在每个增量中接收的热通量(请参阅刀具路径-网格相交模块)。

此外，ABAQUS/Standard提供了许多简化的解决方案，以指定与增材制造工艺中一些最常见的热源相对应的热事件，这些热源不需要您编写用户子程序。这些技术在增材制造专用技术中进行了描述。

控制模拟的规模和解决方案的保真度

增材制造过程仿真是一个时间和空间上的多尺度问题。通过选择适当的时间增量和网格大小，可以控制模拟的规模和解决方案的保真度。一般来说，您可以考虑保真度范围两种极端类型的模拟:过程级模拟(高保真)与零件级模拟(低保真度)。

通过使用小的时间增量和精细的网格来执行详细的工艺级模拟，每个印刷层厚度至少有一个单元，并且在发生主动熔化或熔合的作用区域内有数个单元。该模拟捕获了通常在作用区域内和附近发现的快速变化的温度和高温梯度，因此可以准确预测残余应力和变形。

可以指定与温度相关的热和力学材料属性。在热传递分析中，可以使用潜热定义对熔化和凝固过程中的热能释放和吸收进行模拟(请参见潜热)。对于使用金属材料的工艺，也可以在应力分析中模拟退火效应或在特定温度下的熔化(请参见退火或熔化)。此外，温度历史的精确模拟有助于对印刷过程中的材料相变和微结构演变进行建模，这可以通过用户子程序(如UMAT、UMATHT或USDFLD)来实现。例如，可将不同材料相(如粉末、液体或不同固相)的材料属性指定为代表材料相的场变量或溶液相关变量的函数，并在USDFLD中作为温度历史的函数进行评估(请参阅指定场变量相关性和将材料数据指定为溶液相关变量的函数)。工艺级模拟可以对增材制造过程的详细物理特性进行建模并提供准确的结果，但它通常具有较高的计算成本，并且由于在快速变化的温度条件下使用与温度相关的非线性材料属性，可能会受到收敛困难的影响。

通过对事件的时间顺序进行适当的平均(集中)并使用较粗的网格来执行零件级模拟。例如，模型可以具有比物理打印层厚度大几倍的单元尺寸，并且仅使用一个或几个时间增量来打印一个单元层。只要渐进加热和冷却的热能平衡被正确建模，热传递分析通常可以捕获远场温度演变(远离作用区)。然而，模拟可能无法正确捕捉局部快速温度变化，因为指定的集中、快速移动的热源序列在时间和空间上都集中在一起。换句话说，温度结果通常不包含准确的熔化和凝固历史。在这种情况下，为了正确模拟应力分析中的退火或熔化效应，您必须指定一个初始温度，该温度表示松弛温度，高于该温度时热应变在打印部件中引起的热应力可以忽略不计。在单元激活时，松弛温度是发生初始热收缩的温度。该弛豫温度可以通过实验或详细的工艺级模拟来校准。零件级模拟在计算上有效地预测打印零件中的变形和应力，并具有合理的精度。

解决收敛困难

由于与温度相关的材料行为、渐进式元件激活或其他因素，增材制造过程的热—力学模拟通常是非线性的。下面讨论了常见的收敛问题以及如何解决这些问题的建议。

在传热分析中，收敛困难通常是由材料性质的非线性引起的。

• 如果为材料指定了与温度相关的电导率，则方程组是不对称的。在这种情况下，使用非对称矩阵存储和求解方案（参见 Abaqus/Standard 中的矩阵存储和求解方案）可能有助于提高收敛性。

• 当考虑潜热效应时，它们会在分析中引入严重的非线性并对收敛产生负面影响。 Abaqus 提供了通过指定平滑过渡或延长固相线—液相线温度区间来减少非线性的选项（请参阅定义平滑潜热过渡）

在应力分析中，在激活单元的过程中会遇到收敛困难。

• 单元在被激活之前可能会过度扭曲从而导致收敛困难。在这种情况下，您应该指定非活动单元跟随变形以防止单元过度变形（请参阅控制非活动单元的行为）。
• 单元在激活时经历的初始热应变（由于当前温度与指定初始温度之间的差异）可能导致收敛困难，而这无法通过减小时间增量来解决。为了克服这个问题，Abaqus 提供了一个在一段时间内增加热应变的选项，而不是在激活时立即应用热应变（请参阅应用初始热应变）。渐变热应变会影响分析结果的准确性。例如，如果在激活下一层单元时，一层中的单元的热应变还在变化，则新激活的单元的无应变配置与热应变完全应用时的情况不同。应该使用小于处理一层所需时间的渐变时间常数。
• 如果材料定义包括塑性，则分析可能会由于用于加快求解速度的外推方案而过度迭代。可以通过关闭外推来防止此问题（请参阅 Abaqus/Standard 中的增量）。

输入文件模板

以下模板显示了传热分析的输入 heat.inp:

HEADING
…
** Choose the uncoupled heat transfer element type
*ELEMENT, TYPE=DC3D8
** Define material deposition motion and/or heat source scan path
**EVENT SERIES TYPE
*EVENT SERIES
** Define additional information for material deposition and/or moving heat sources
*PARAMETER TABLE TYPE
*PROPERTY TABLE TYPE
*TABLE COLLECTION
*PARAMETER TABLE
*PROPERTY TABLE
** Define elements that can be activated during an analysis
*ELEMENT PROGRESSIVE ACTIVATION
…
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=θZθZ, STEFAN BOLTZMANN=σσ
*INITIAL CONDITIONS, TYPE=TEMPERATURE
…
*STEP, UNSYMM=YES
*HEAT TRANSFER
…
** Turn on progressive element activation in the step
*ACTIVATE ELEMENTS
*table collection for material deposition (if defined)*
** Apply thermal loads and boundary conditions
*BOUNDARY
*Data line to define nonzero boundary conditions on temperature degree of freedom*
*FILM
*element number or element set name, FFS or FFSNU, θ0θ0, h*
*RADIATE
*element number or element set name, RFS, θ0θ0, ϵϵ
*DFLUX
*element set name, MBFNU, blank entry, table collection for a heat source (if defined)*
…
** Write all nodal temperatures to the output database file, heat.odb
*OUTPUT, FIELD
*NODE OUTPUT
NT
*END STEP

以下模板显示了后续静态结构分析的输入：

HEADING
…
** Choose the continuum element type compatible with the heat transfer element type used
ELEMENT, TYPE=C3D8
…
** Define material deposition motion
EVENT SERIES TYPE
EVENT SERIES
** Define additional information for material deposition 
PARAMETER TABLE TYPE
PROPERTY TABLE TYPE
TABLE COLLECTION
PARAMETER TABLE
PROPERTY TABLE
** Define elements that can be activated during an analysis
ELEMENT PROGRESSIVE ACTIVATION
…
INITIAL CONDITIONS, TYPE=TEMPERATURE
BOUNDARY
Data lines to specify zero-valued boundary conditions on displacement degree of freedom
…
STEP, EXTRAPOLATION=NO
STATIC
…
** Read in all nodal temperatures from the output database file, heat.odb
TEMPERATURE, FILE=heat.odb
** Turn on progressive element activation in the step
ACTIVATE ELEMENTS
table collection for material deposition (if defined)
** Apply structural loads and boundary conditions
…
END STEP