增材模拟论文精读

背景：为什么做

• 金属增材制造涉及大量相互耦合的工艺和材料参数，妨碍参数优化、材料开发以及实时监测与控制。主文 p1-p2
• 匙孔由蒸发诱导反冲压力形成，其形貌会改变激光能量耦合和熔池流动；不稳定匙孔还与孔隙、飞溅、球化和异常组织相联系。主文 p2
• 已有研究识别了未熔合、凹陷区失稳、挥发性元素蒸发和析氢等成孔机制，但仍缺少能定量区分各机制贡献并预测孔隙类型和幅度的模型。主文 p2
• 研究目标是用量纲分析寻找跨材料、工艺条件和设备仍成立的紧凑标度律，并用原位数据和多物理场模型解释其物理来源。主文 p1-p2

研究方法

模型

• 三维模型耦合气相与基体相体积分数、不可压 Navier-Stokes 动量方程、焓形式能量方程及金属蒸气组分输运。SI p2-p4
• 模型包含激光多次反射、传热传质、Marangoni 流、表面张力、蒸发反冲压力、熔化凝固和金属蒸气羽流；激光吸收用 Fresnel 关系和射线追踪计算。SI p2-p4
• 模型不包含等离子体生成，也不考虑蒸气羽流或等离子体造成的激光衰减。主文 p3
• 能量收支逐时计算反射、对流、辐射、蒸发和传导；飞溅功率因模型不能准确给出，改由 X 射线图像中液滴体积近似。SI p5-p6

数值方法

• 控制方程用 ANSYS FLUENT 2020 R1 的有限体积法求解，采用非迭代 PISO、二阶隐式时间离散、PRESTO 压力格式和二阶迎风动量离散。SI p4
• 界面采用显式 VOF 与 CICSAM 格式；初始六面体网格边长 8 μm，高于固相线温度的区域做两级局部加密，最小边长 2 μm。SI p4
• 时间步在 1×10^-10 至 1×10^-8 s 之间自动变化，以保持全局 Courant 数小于 1。SI p4
• 共做 30 个不同功率和扫描速度的多物理场算例，其中 Al6061 为 10 例，Ti-6Al-4V 为 20 例。SI p5；主文 p5
• 作者从归一化控制方程和边界条件做量纲分析，建立 Keyhole number、归一化扩散长度等低维参数，再用实验和模拟数据拟合并检验标度关系。主文 p2；SI p6-p9

实验

• 作者在 Al6061 和 SS316 裸板上做单道激光熔化，激光功率 208 至 520 W，扫描速度 0.3 至 1.2 m/s，并用氩气保护。主文 p7
• 原位成像在 Argonne Advanced Photon Source 32-ID-B 光束线完成，标称空间分辨率 1.93 μm/pixel，帧率 20,000 至 50,000 fps，曝光时间 1 至 40 μs。主文 p7
• 匙孔和熔池轮廓经 ImageJ 处理后逐帧人工检查与测量，记录各图像栈的最大值、最小值、均值和标准差。主文 p7
• 孔隙率标度律使用六项既有独立研究的数据，覆盖不同材料、工艺参数和粉末床设备；这部分不是本文新做的孔隙实验。主文 p5-p6

结果：算了什么、做了什么实验

• 匙孔纵横比 e*=e/r0 与 Keyhole number 呈近线性关系 e*=0.4(Ke-1.4)。作者据此划分导热模式 Ke<1.4、过渡模式 1.4≤Ke≤6.0 和匙孔模式 Ke>6.0。主文 p2
• 匙孔纵横比波动满足 Δe*=0.36Ke^0.86；Ke<16 时波动小、几乎无孔且少见飞溅，Ke>30 时匙孔表现出混沌和失稳。主文 p3
• 30 个多物理场算例给出吸收率关系 η=0.7[1-exp(-0.6KemLd*)]，并与文献中的实验吸收率数据比较。主文 p5
• 能量分解显示吸收和反射占主导，对流、蒸发、辐射和飞溅的归一化功率之和约为 0.1，辐射损失可忽略。主文 p5
• 孔隙率由归一化能量密度 NED 与 KemLd* 两个量控制，二维标度律对汇总数据的决定系数 R² 为 0.92；NED 约高于 5 后，未熔合孔隙可避免，而匙孔孔隙项开始主导。主文 p5
• 作者把匙孔稳定性与孔隙形成定量联系起来，并指出 Keyhole number 出现在匙孔形貌和孔隙率两个独立问题中。主文 p5、p7

讨论：围绕什么展开

• 作者认为方法可推广到激光焊接、熔覆及其他能束制造，但电子束的吸收机制和真空环境不同，需要修改标度律。主文 p6
• 本文 X 射线实验使用裸板。作者依据既有研究认为 50 至 100 μm 粉层对匙孔深度或匙孔孔隙边界影响较小，但承认粉层可能使拟合常数需要调整。主文 p6；审稿 p5、p9-p10
• 匙孔标度律只适用于充分发展的准稳态熔池，不适用于扫描道起止位置和激光转向处的瞬态熔池。主文 p6-p7；审稿 p8-p9、p16-p17
• 基板温度和环境压力在实验中固定，因此本文没有量化它们对匙孔形成与稳定性的影响。主文 p2；审稿 p7
• 审稿讨论的主要焦点是裸板结果对粉末床、多道多层和大尺度打印的适用性，以及准稳态限制、电子束推广和蒸发参数在模型与标度律中的不同作用。审稿 p1-p3、p8-p14、p16-p17

创新点

• 提出 Keyhole number，用单一无量纲量统一描述三种材料及不同功率、速度和光斑条件下的匙孔纵横比与稳定性。主文 p1-p3
• 用超高速原位 X 射线测量传统事后截面无法观察的匙孔波动，并给出波动幅度与 Ke 的定量标度律。主文 p3
• 以包含射线追踪和蒸气羽流的高保真模型分解瞬态能量收支，为吸收率标度和变量筛选提供物理解释。主文 p3-p5；SI p2-p5
• 把原先由十个参数控制的孔隙率问题压缩为 NED 与 KemLd* 的二维关系，并定量连接未熔合孔隙与匙孔孔隙两个区域。主文 p5-p6

背景：为什么做

• 增材制造成形件中的孔隙会降低极限强度、疲劳强度和断裂强度，匙孔模式下的孔隙是激光粉末床熔融与激光焊接中的常见缺陷。主文 p1
• 既有离位实验能看到孔的形貌与位置，却不能直接观察形成过程；原位X射线实验已经辨认出不同类型的匙孔孔隙。主文 p1
• 作者认为局部反冲压力分布、能量分布与匙孔波动的关系、以及孔隙形成时的熔池流动仍缺少机制解释。主文 p1

研究方法

模型

• 模型求解不可压缩牛顿层流的质量、动量和能量守恒，并用Boussinesq近似处理浮力。主文 p9
• 模型包括传热、熔池流动、金属蒸发产生的反冲压力、Marangoni效应、表面张力和糊状区Darcy阻力。主文 p1、p9至10
• 激光多次反射由射线追踪和Fresnel方程计算，采用镜面反射假设；封闭气泡按绝热理想气体处理。主文 p9至10

数值方法

• 采用有限体积法，熔池自由表面用体积分数法捕捉，糊状区阻力采用Blake-Kozeny模型。主文 p9
• 所有主文算例网格尺寸为4 μm，物理时间为2000 μs；具体软件或代码名、计算域和求解器时间步均未给出。主文 p2、p9
• 主文设置Case 1至5，补充材料另做一个关闭大部分Darcy阻力的Case 2对照算例。SI p1至2

实验

• 本研究未报告新做的实验；仿真参数除低压算例外取自Zhao等人的原位X射线实验，并用其图像校验气泡形成与捕获过程。主文 p2；SI p1
• 所有仿真使用裸板以排除随机铺粉影响，作者明确指出粉床影响仍需另行研究。主文 p1至2

结果：算了什么、做了什么实验

• 匙孔孔隙形成包含两个阶段：不稳定匙孔塌陷产生瞬时气泡，随后气泡被凝固前沿捕获。主文 p2、p8
• 匙孔后壁局部反冲压力减弱和反向流促使塌陷；气泡下方的高速流造成较低局部压力，阻碍气泡上浮并促使其近水平移动。主文 p2、p4
• 在相同激光功率下，扫描速度升高时孔隙尺寸显著减小，孔形趋于球形并沿熔道底部水平分布。主文 p4至6、p8
• 低环境压力算例中没有形成匙孔孔隙；后壁反冲压力增大，使匙孔形状更稳定。主文 p6、p8
• 低压下匙孔深度波动反而更大，因此深度波动不能单独表征稳定性；能量吸收率波动对形状变化更敏感。主文 p8

讨论：围绕什么展开

• 当前模型不能模拟声波驱动的孔隙机制和相爆引起的激波，论文的机制讨论集中在拖曳力相关的第二类孔隙。主文 p2
• 仿真与实验在临界工况下的最大孔尺寸和匙孔深度存在偏差，作者归因于过程随机性以及模型未包含离焦、等离子体和羽流等影响。主文 p2、p4至5
• 作者把匙孔深度、深度波动和能量吸收率变化讨论为孔隙风险指标，其中后壁吸收率的波动与匙孔形状变化联系最直接。主文 p5至6、p8
• 现有Blake-Kozeny阻力是各向同性的，而增材制造糊状区多为柱状晶；作者认为后续模型应考虑枝晶形貌导致的阻力各向异性。SI p3
• 低压能减少匙孔孔隙，但会加剧铝、镁、锰等元素的蒸发损失，原料成分可能需要相应调整。主文 p8

创新点

• 论文把匙孔塌陷后的瞬时气泡生成与气泡最终被凝固前沿捕获明确分成两个阶段，并分别解释驱动力。主文 p2、p8
• 模型把局部反冲压力、糊状区Darcy阻力和气泡周围的速度压差放入同一形成机制，并与原位X射线图像逐阶段比较。主文 p2至4；SI p1至3
• 论文提出用匙孔表面能量吸收率的波动辅助判断匙孔形状稳定性，而不只看匙孔深度波动。主文 p6、p8
• 论文用低压算例说明后壁反冲压力分布改变可以抑制匙孔塌陷，并同时讨论低压造成的元素挥发代价。主文 p6、p8

背景：为什么做

• DED可制造或修复高价值构件，但孔隙会降低力学性能，尤其是疲劳寿命，限制其安全关键应用。主文 p1
• 金相和离线X射线CT只能看到凝固后的孔，不能直接捕捉孔形成、长大和迁移的动态过程。主文 p1
• 既有原位研究主要针对LPBF，已有DED研究又多处于高能量密度的匙孔模式；工业相关的传导模式DED具有更大的熔池和粉末撞击，孔演化物理并不相同。主文 p2
• 已有DED多物理场模型主要研究沉积、铺层和流场，没有系统包含孔的形成、迁移、合并和逸出，也缺少内部X射线数据验证。主文 p2

研究方法

模型

• 模型把质量、动量和温度方程全耦合求解（fully coupled），包含相变、糊状区Darcy阻力、重力、表面张力与Marangoni效应、导热、潜热和辐射；激光由射线追踪施加，沉积质量按质量守恒加入。主文 p8、p12；SI Supplementary Method 1
• 自由表面用CLSVOF捕捉；小气泡按Lagrangian点粒子追踪，大气泡用Eulerian level-set显式解析气液界面、形状和合并。点粒子适用性由Reynolds数、Weber数和Stokes松弛时间判定。主文 p12；SI Supplementary Method 1、Supplementary Fig. 15；审稿 p26至p27
• 粉末撞击用两种方式建模：强制表面速度与温度扰动，以及直接粉末撞击；前者施加驻波，后者模拟颗粒撞击、局部冷却和流胞重排。主文 p9至p12

数值方法

• 数值格式基于CIP三次插值。基础流场和三处气泡合并算例用16 µm网格，长时间气泡轨迹、推挤、逸出和粉末撞击算例用32 µm网格。主文 p12；SI Supplementary Method 1
• 文中没有给求解器或代码名称、完整计算域、单元总数、实际时间步、收敛准则和机时。主文 p12；SI Supplementary Method 1；主文、SI、审稿往来全文检索
• 实验图像在ImageJ和MATLAB中处理，采用平场校正、Otsu分割和分水岭分离粉末或孔。主文 p11至p12；SI Supplementary Fig. 2

实验

• 实验在ESRF ID19束线上用BAMPR-II装置进行，材料是气雾化RR1000镍基高温合金粉末，D50为90 µm。主文 p2、p11
• 研究使用传导模式，实用参数覆盖激光功率100至200 W、基板移动速度1至2 mm/s、送粉率1.8至2.7 g/min；实测光斑约390 µm。主文 p11；SI Supplementary Discussion 1、Supplementary Table 2
• 高速成像的空间分辨率为4 µm、帧率为20 kHz；较长时段成像的空间分辨率为3.7 µm、帧率为1 kHz。作者还用断层成像检查整条沉积道及道尾大孔。主文 p4、p6、p11；SI Supplementary Fig. 4

结果：算了什么、做了什么实验

• 作者把孔演化归纳为五个环节：粉末或前层带入气泡、气泡合并长大、凝固界面推挤、Marangoni流维持大泡滞留，以及最终逸出或凝固捕获。主文 p2至p4
• 在所用气雾化粉末和传导模式下，粉末内部氩孔是主源，其输入孔数约为其他来源的2至4倍；作者没有观察到匙孔孔或颗粒入池卷气。主文 p4至p5
• 小泡在回流中合并成大泡。移动速度为较低设置时，大泡最大直径约180 µm；较高设置时临界直径约120 µm。摘要给出多至约 70 个直径 20–50 µm 的小泡合并成一个大孔；正文给的是形成最大气泡平均要 30 个以上（2 mm/s）与 70 个以上（1 mm/s）。主文 p4、p6、p10
• 定量计数中，逸出、捕获和合并速率分别约为1.09、0.05和0.07个每毫秒，说明在该工况下多数小泡从表面逸出，少量被凝固前沿捕获或并入大泡。主文 p6、p8；SI Supplementary Fig. 7
• 计算再现了前、后及两侧的四个回流胞、三处气泡合并、Marangoni剪切流压住大泡及大泡逸出，并显示粉末撞击可打乱原流场、形成小流胞并改变气泡轨迹。主文 p8至p10；SI Supplementary Fig. 8至13

讨论：围绕什么展开

• 大泡临界尺寸与熔池大小相关：较低移动速度、较高激光功率和较大送粉率产生更大的熔池，也允许大泡继续合并长大；层数影响不明显。主文 p6；SI Supplementary Fig. 5至6
• 气泡去向由尺寸、位置和局部流结构共同决定。小泡可随回流循环、逸出或被界面捕获；大泡会压缩其上方Marangoni流，压力和剪切合力可向下抵消浮力，直到气泡长大或流场被扰动。主文 p6至p9；SI Supplementary Discussion 2、Supplementary Table 3
• 与先前匙孔模式DED研究不同，本研究把粉末内孔识别为主源。作者认为差别来自气雾化粉末、传导模式和更大的熔池，而不是一种普适的新孔源。主文 p5、p10；审稿 p16
• 审稿讨论集中在工业代表性、二维投影能否证明合并、受力解释是否完整，以及这些现象相对既有LPBF、焊接和凝固研究的新意。作者据此补充了参数对照、受力分解、点粒子适用范围和模型限制。审稿 p2至p4、p8至p11、p17至p28
• 模型没有直接研究不同尺寸气泡在凝固界面处的推挤与捕获阈值；二维X射线也存在厚度方向重叠歧义。这两点限制了对界面机制和单次合并事件的定量判定。审稿 p9、p11

创新点

• 在 DED 的大熔池内直接追踪到小泡反复合并成大泡并周期性逸出的过程，并量化出多至约 70 个直径 20–50 µm 的小泡合并成一个约 180 µm 的大孔。主文 p6、p10至p11；审稿 p7、p25
• 研究给出了DED中凝固界面推挤大泡的原位证据，并区分稳态推挤与停激光后大小气泡共同被捕获的情形。主文 p3至p4、p10至p11
• 原位高速X射线和三维多物理场模型互相约束，使气泡迁移、合并、受压滞留和逸出不只停留在表面观测或定性示意。主文 p2、p8至p10；SI Supplementary Discussion 2
• 粉末撞击被拆成驻波扰动和直接撞击两类算例，用来解释随机入射下仍出现较有组织的气泡运动；不过审稿人对这一解释和整体创新幅度仍保留意见。主文 p9至p10、p12；审稿 p21至p22、p25、p28至p29

背景：为什么做

• PBF工艺参数会改变成形组织，组织又影响最终力学性能；参数众多使试错实验成本高、耗时长，因此作者把数值模拟作为研究工艺与组织关系的手段。主文 p1
• 既有增材制造晶粒模拟以元胞自动机和二维相场为主，三维相场工作较少，且主要限于单道。主文 p1-2
• 既有模型很少处理循环热作用下热影响区的晶粒粗化，也通常把离散粉床连续化；把离散多晶粉末颗粒纳入晶粒演化仍有缺口。主文 p2
• 作者以316L不锈钢为对象，目标是在三层三道PBF中同时处理形核、外延与竞争生长、搭接区重熔再生长及热影响区粗化，并用实验和文献中的纳米颗粒细化结果作比较。主文 p2

研究方法

模型

• 粉末尺度热流体流动模型解析单个粉末颗粒，求解质量、动量和能量守恒，并考虑导热、熔化潜热、表面张力、Marangoni效应、黏性、蒸发和反冲压力。主文 p9
• 相场模型用非守恒变量ζ区分液相与固相，用η_i区分不同取向晶粒，并通过时变Ginzburg-Landau方程描述液固转变、晶粒生长和粗化。主文 p9-10
• 晶界能随晶界相对晶体取向呈各向异性，晶界迁移率随温度变化，使晶粒在温度场中发生取向相关的竞争生长。主文 p10
• 形核模型把经典形核理论给出的异质形核率接入Poisson seed方法；模型按局部相场变量自动选取液固界面，在网格单元中随机引入晶核，无需预设形核位置和临界过冷度。主文 p10-11
• 热流体模型与相场模型采用单向耦合：先输出若干时刻的温度场，再在空间和时间上线性插值后输入更细的相场网格。主文 p11
• 模型忽略成分场与液固两相间的成分扩散，并令液固转变由温度控制。主文 p9

数值方法

• 主算例的热流体模型使用商业求解器Flow3D，计算域为2 × 0.4 × 0.2 mm³，网格尺寸4 μm，文中报告约300万个单元。主文 p2、p11；SI p5
• 主相场算例使用自研GPU并行代码，计算域为350 × 350 × 146 μm³，网格尺寸1 μm，约1788万个单元，时间步长1 μs。主文 p2、p11；SI p7
• 相场温度历史覆盖三层、每层三道，共九道扫描；相场物理时间约0.009 s。主文 p11
• 初始组织只使用36个随机晶体取向；单个粉末颗粒按尺寸放入约5至10个小晶粒。主文 p2；SI p1-2

实验

• 作者用商业DiMetal-100设备制备10 × 10 × 30 mm³的316L矩形样品，粉末粒径20至40 μm，采用相邻层旋转90°的扫描路径，截面经抛光和化学腐蚀后用SEM及EBSD观察。主文 p9
• 验证实验Case I采用190 W激光功率、1000 mm·s⁻¹扫描速度、80 μm道间距和30 μm层厚。SI p2
• 纳米颗粒细化算例采用文献实验Case II的参数：100 W、83.33 mm·s⁻¹、120 μm道间距和50 μm层厚；该实验并非本文作者新做。主文 p8、p11；SI p2

结果：算了什么、做了什么实验

• 三层三道模拟再现了熔池内形核和生长、粉末与基板及前一道晶粒的外延生长、搭接区重熔再生长、竞争生长，以及热影响区粗化。主文 p2-5
• 第三层第三道的模拟熔池平均宽度和深度为111.9 μm与71.8 μm；顶层实验值为116.5 μm与79.2 μm，作者据此认为熔池形貌一致。主文 p3
• 全部扫描道的晶粒中，83.4%的长短轴比大于1.5，按文中判据属于柱状晶；长短轴比分布与实验结果相符。主文 p6
• 同层内平均晶粒体积由第一道到第三道增大，层间则由下层向上层增大；作者分别归因于热积累和下层晶粒向上层的竞争延伸。主文 p6
• 第一层第三道的形核数略多于第二层第三道；作者把这一差别解释为温度梯度与冷却速率的共同作用，并认为本算例中冷却速率变化的影响更大。主文 p6-8
• 加入TiB₂纳米颗粒的简化单道模拟定性再现了熔池内部大量细晶和边界处拉长晶粒，但模拟细化晶粒大于文献实验值。主文 p8

讨论：围绕什么展开

• 层间扫描方向旋转90°会同步改变热传导方向和优先生长方向，形成跨层折线状晶粒，并降低形成强织构的可能性。主文 p5
• 基板、部分熔化粉末和未重熔下层中的晶粒粗化，来自后续扫描造成的升温、晶界迁移率提高及原位热处理。主文 p5
• 作者讨论了温度梯度和冷却速率对连续冷却形核的耦合作用，并提醒恒定拉速实验中改变温度梯度也会改变冷却速率。主文 p7-8
• 纳米颗粒既可提供异质形核位置，也可能钉扎晶界；当前模型只含前一作用，未包含Zener钉扎，因此只能定性复现细化且高估晶粒尺寸。主文 p8
• 作者把成分扩散缺失、可模拟层道数有限和体素表面不光滑列为局限，并讨论用机器学习温度场、自适应网格、多GPU或表面重构扩展模型的可能性。主文 p8-9

创新点

• 把经典形核理论和Poisson随机事件接入三维相场模型，并由局部相场变量自动确定异质形核区域，避免预设形核位置和临界过冷度。主文 p8、p10-11
• 在同一三维三层三道算例中联合处理液固转变、形核、竞争与外延生长、搭接重熔再生长及热影响区粗化。主文 p2、p8
• 初始粉床保留离散粉末颗粒，并把每个颗粒表示为多晶而非单晶，从而能够模拟粉末内部晶粒的外延和粗化。主文 p2
• 模型除与本文制备样品的熔池和晶粒形貌比较外，还被用于纳米颗粒诱导细化的独立场景，展示了模型的扩展用途，同时暴露出未含钉扎效应的边界。主文 p3-8

背景：为什么做

• NIST AM Bench 2022 Challenge 3 要求盲预测 IN718 激光粉末床熔融中单道和多道的固相冷却速率、液相冷却速率、熔化温度以上停留时间及熔池几何。主文 p1-p2
• 高保真热流体模型需要处理 Marangoni 流、汽化、反冲压力、自由表面演化和激光多次反射，扩展到零件尺度时计算代价高。主文 p2
• 作者要解决的问题是：能否用经标定的简化圆柱热源吸收部分缺失物理，在激光功率、扫描速度和光斑直径变化时仍给出定量预测。主文 p2

研究方法

模型

• 模型是三维瞬态热流体模型，求解连续性方程、Navier-Stokes 方程和能量方程，并包含 Marangoni 驱动流动、糊状区阻力、辐射、对流、汽化热损失和相变潜热。主文 p9-p10
• 激光以径向为高斯分布的圆柱体积热源表示，热源深度、吸收率和半径由启发式、基于小孔标度律、物理引导三种方案标定。主文 p10-p12；SI Supplementary Methods 1-2
• 物理引导方案把体积能量密度与热源深度、吸收率及半径联系起来，其中半径关系经过修改，使热源纵横比随体积能量密度变化。主文 p12；SI Supplementary Method 2
• 模型没有求解自由表面演化、激光多次反射或保护气流，因此圆柱热源决定了熔池截面的基本形状。主文 p6、p9

数值方法

• 作者使用西北大学自研的 C++ 有限体积 AM-CFD 求解器，并以 Higher-Order Proper Generalized Decomposition，简称 HOPGD，构建非侵入式降阶代理模型来辨识三个热源参数。主文 p2、p9-p12；SI Supplementary Method 5
• AM-CFD 采用随激光移动的局部高分辨率自适应网格；网格位置变化时，状态变量通过线性插值传递到下一时间步。主文 p10
• 单道计算没有沿用实验的 10 mm 全轨迹，而是按熔池达到稳态所需的距离缩短计算域；最高和最低体积能量密度工况分别约在 3 ms 和 1 ms 达到稳态。主文 p3；SI Supplementary Method 3

实验

• 验证对象是 NIST AM Bench 2022 Challenge 3 的 IN718 裸板实验，包括七种单道激光参数和基准参数下的 X、Y 两种多道扫描。主文 p2、p8-p9；SI Supplementary Method 6
• 原位热成像给出固相与液相冷却速率和熔化温度以上停留时间，打印后截面测量给出熔池几何。主文 p1、p9；SI Supplementary Method 6
• 固相冷却速率由 1260 到 1150 摄氏度区间的线性拟合斜率求得，液相冷却速率使用 1400 到 1336 摄氏度区间，熔上时间以 1298 摄氏度为阈值。主文 p3-p5；SI Supplementary Method 6

结果：算了什么、做了什么实验

• 物理引导方案总体改善了七个单道工况的冷却速率、熔上时间和熔池几何预测，尤其改善了高体积能量密度下启发式方案的偏差。主文 p3-p7
• 正文称物理引导方案把熔上时间差异降到 20% 以内，但表 2 中仍有多个工况超过 20%，因此这不是对全部七个工况都成立的界限。主文 p5、p7
• 正文称物理引导方案把熔池几何预测带到实验值的 5% 以内；表 2 显示宽度符合这一表述，但深度有一个工况超过 5%。主文 p5、p7
• 多道计算的冷却速率、熔上时间和尺寸多数与实验相差 10% 以内，但熔池重叠深度和宽度的形貌差异约为 55%。主文 p6、p8-p9
• 作者指出，在高能量密度条件下忽略激光反射、吸收和金属蒸气中的传播，可使冷却速率和熔上时间高估约 100%。主文 p9

讨论：围绕什么展开

• 讨论的主线是简化热源能否代替昂贵的自由表面、光线追迹和多物理场计算。定量标量可以通过标定改善，但熔池截面形貌仍受热源先验形状限制。主文 p2、p6、p9
• 高体积能量密度下，小孔和汽化改变了真实热输入；启发式方案为维持熔池尺寸而输入过多能量，从而抬高峰值温度和冷却速率。主文 p3-p4
• Y 向多道截面出现交替深浅。作者比较了达到稳态所需时间和负 Y 向保护气流两种解释，并因模型未含气流而把后者列为假说，并非已验证结论。主文 p6；SI Supplementary Method 6
• 作者把该方法定位为走向多道、多层和零件尺度计算的近似方案，而不是完整的高保真熔池动力学模型。主文 p9

创新点

• 作者提出基于小孔标度律的圆柱热源标定方案，以及把体积能量密度与熔池纵横比相关联的物理引导标定方案。主文 p2、p11-p12；SI Supplementary Methods 1-2
• 物理引导方案不新增自由表面方程，而是修改热源参数随体积能量密度的关系，用较简单的有限体积模型覆盖不同功率、速度和光斑。主文 p9、p12
• HOPGD 代理模型用于从已有 AM-CFD 数据和文献熔池数据中辨识热源参数，减少反复调用全模型的标定成本。主文 p2、p10-p11；SI Supplementary Method 5
• 论文把三种标定方案放到 NIST 盲测框架下，同时比较单道热历史、熔池尺寸和多道结果。主文 p1-p9

背景：为什么做

• 激光粉末床熔融逐层扫描时会随机积累孔洞等缺陷，使零件难以稳定通过资格认证；即便致密度超过99%，少量微孔也会影响力学性能。主文 p1
• 蒸气射流和气流会卷吸周围粉末并产生飞溅、去粉区与再沉积，而现有超高速X射线仍可能缺少捕捉熔池形貌和微孔快速变化所需的时空分辨率。主文 p1
• 研究目标是把激光参数、粉末、飞溅和熔池的相互依赖关系连接到缺陷形成，并建立可控制瞬态过程的判据。主文 p1，p6

研究方法

模型

• 作者使用ALE3D建立热传导与流体动力学耦合模型，并加入基于非结构三维网格的完整激光射线追踪和温度相关材料性质，使吸收率成为模型输出而非常数输入。SI p2，p29；主文 p2
• 模型不求解外部气体动力学，也不求解孔内气体；蒸气反冲压力作为液面边界条件，粉末初始化为固定的预烧结状态，飞溅场景由数字实验人工加入。主文 p2；SI p2
• 作者以Weber数比较惯性压力和表面张力，并用比例积分微分控制器调节功率，使激光下方示踪点保持在熔点和给定深度轨迹上，由此得到功率映射。主文 p3-p5；SI p6-p7

数值方法

• 模型输出三维温度和速度等场量，空间尺度约为2 µm；传热与流体动力学在纳秒尺度耦合，但文中没有给出网格拓扑或数值时间步长。主文 p1
• 起始段孔洞机制的可视化从共6000条射线中随机显示100条；SI还提到仿真以0.5 µs尺度观察波动，但没有把它定义为求解时间步长。SI p8
• 比例积分微分控制器用Ziegler-Nichols方法分阶段调节P、PI和PID项；训练算例数量没有给出。SI p7

实验

• 原位透射X射线实验在Advanced Photon Source的32-ID和35-ID束线开展，配合连续波光纤激光、振镜扫描和高速相机，记录单道熔化、飞溅、遮光与起始段回溅。SI p2-p4，p7
• 作者还用EOS M280打印Ti-6Al-4V样品，通过背散射电子显微镜连续截面、Fiji/ImageJ配准、DREAM.3D重建和ParaView显示终点冻结凹坑。主文 p5；SI p9

结果：算了什么、做了什么实验

• 完整射线追踪再现了裸板从传导模式向近钥孔模式的吸收率变化；200 W以上时，带粉末与裸板的吸收率曲线重合，粉末作用减弱。主文 p1-p2；SI p2
• 飞溅数字实验中，300 W激光能将大团聚体侧向排出，150 W则会把团聚体压入熔池并造成轨迹不连续；正面碰撞还会把一个团聚体裂解为多个新飞溅点。主文 p2-p3
• 当飞溅尺寸接近光斑并遮住高强度中心时，快速冷却会形成微孔；对SS316L和50 µm光斑，作者推导的排出阈值为排出阈值 P_threshold 约等于 157 乘以扫描速度（主文原话 power is equal to ~157 times the scan speed），其中功率单位为 W、扫描速度单位为 m/s。主文 p3；SI p5-p6
• Ti-6Al-4V和SS316L的临界流速分别为6.1和3.8 m/s，对应的熔深变化率上限分别约为1.5和0.9 m/s；超过该尺度的起始段高速回流会产生大尺寸回溅。主文 p3-p4；SI p6
• 功率映射把起始段熔深变化率降到0.13 m/s，抑制回溅和起始孔洞；缓慢降功率还避免了约43 µm的终点冻结凹坑，并在转弯时压制钥孔和孔洞。主文 p5；SI p8-p10

讨论：围绕什么展开

• 作者把工艺窗口解释为两类约束的夹区：功率必须高到足以排出遮光飞溅，又不能高到引发钥孔、过强瞬态回流或大尺寸回溅。主文 p6；SI p6
• 预烧结或多束预处理可减少松散粉末运动和飞溅，但作者明确指出还要处理粉末可回收性和氧化问题。主文 p2-p3
• 稳定性判据给的是非理想几何下的参考尺度，并非精确闭式预测；控制示踪点轨迹的方法被用于起始、终止和转弯扫描。主文 p3-p5；SI p6-p7
• 模型省略气体流动并固定粉末，这降低了计算代价，却也意味着气体卷吸、孔内气体和最终孔形不能由这些算例完整给出。主文 p2；SI p2，p8

创新点

• 把完整激光射线追踪并入热流体模型，让吸收率随粉末和熔池形貌自动变化，覆盖传导到近钥孔的功率范围。主文 p1-p2；SI p2
• 提出并用仿真与原位X射线证据连接了激光排出、自复制飞溅和飞溅遮光致孔三种机制。主文 p2-p3；SI p3-p6
• 用Weber数稳定性尺度和反馈控制生成位置随时间变化的功率映射，同时处理扫描起始、终止与转弯处的瞬态缺陷。主文 p3-p6；SI p6-p10

背景：为什么做

• 激光粉末床熔融可直接制造复杂金属零件，但孔隙会损害耐腐蚀性和疲劳寿命；深窄匙孔的波动与塌陷是孔隙来源之一。主文 p.1
• 功率和扫描速度工艺图长期偏经验化，匙孔孔隙边界是否清晰、可复现，以及粉末层会把边界移动多少，此前缺少定量答案。主文 p.1
• 孔洞形成跨越微米和亚微秒尺度，常规光学观测难以同时满足穿透、空间分辨率与帧率要求。主文 p.1；补充材料 p.2、p.6

研究方法

模型

• 作者把匙孔形成归因于金属蒸发反冲压力，并用表面张力、热毛细力、黏性阻力和声波惯性作用解释孔洞的形成与迁移。主文 p.1、p.4、p.6；补充材料 p.7
• 微射流冲击压力用水锤关系估算，匙孔底温度由 Bernoulli 方程和饱和蒸气压估算，黏性阻力系数由热毛细力与黏性阻力平衡估算。补充材料 p.4-p.5

数值方法

• ImageJ 用于背景校正、反差翻转和几何及速度测量；匙孔深度再按激光功率与扫描速度的解析关系做回归。补充材料 p.3、p.5
• 三维射线追迹模型计算给定匙孔壁上的吸收强度，物理处理为多次反射与 Fresnel 吸收；图 S15D 比较三角形底部与针状底部，输入为光斑 100 μm、功率 205 W。具体建模细节转引到 Ref. 9，本文未复述。补充材料 p.4、p.23

实验

• 材料为 Ti-6Al-4V 裸板和粉末床试样，粉末层厚度为 105±8 μm；连续高斯激光沿直线扫描，工艺范围为 84-560 W、150-700 mm/s。补充材料 p.2
• Advanced Photon Source 的原位同步辐射成像使用两台高速相机：工艺图测量采用约 2.0 μm/pixel、50 kHz，细节成像采用约 3.0 μm/pixel、1.087-5.0 MHz 和约 100 ps 单脉冲曝光。补充材料 p.2
• 兆赫兹孔洞形成实验选择靠近孔隙边界的 205 W、500 mm/s 条件；独立的高帧率事件记录见 Movies S11 和 S12。补充材料 p.2、p.32

结果：算了什么、做了什么实验

• 有粉和无粉条件下，稳定熔化区与匙孔孔隙区之间都出现平滑而尖锐的边界，粉末只使边界小幅移动。主文 p.1、p.6
• 边界上的临界匙孔深度满足 dc = a(tanθ - 1)，匙孔孔隙常数 a 约为 83 μm，大于稳定区有效光斑约 62 μm；水平截距接近 tanθ = 1。主文 p.4（文字层里 µ 丢成了 m，单位 µm 已对照 PDF 页面核过）
• 近边界孔洞依次经历微小匙孔形成、瞬时孔洞分离、针状匙孔底部形成，以及孔洞快速远离并被凝固前沿钉扎。主文 p.4
• 匙孔塌陷速度约 13 m/s，微射流以约 16 m/s 撞击孔壁，对应约 140 MPa 的水锤压力；随后孔洞可在不足 1 μs 内由近乎静止加速到约 10 m/s。补充材料 p.4、p.24
• 在所分析条件下，强温度梯度区边缘距匙孔壁 42±2 μm，热毛细力与黏性阻力平衡位置为 28±3 μm，估得阻力系数约 6。主文 p.6；补充材料 p.5

讨论：围绕什么展开

• 主导机制是匙孔尖端失稳产生孔洞，声波及其诱导液体惯性把孔洞推离强温度梯度区；若能量不足，热毛细力会把孔洞重新吸回匙孔。主文 p.4、p.6；补充材料 p.7
• 低功率区域还存在次级机制：粉末遮光加剧匙孔深度波动，孔洞在等待凝固前沿到达时也可能受局部熔体阻力带离匙孔。主文 p.2、p.6；补充材料 p.7
• 作者把针状匙孔底部的爆发行为解释为相爆，但措辞明确是推测；支持证据是从内部开始膨胀、壁面粗糙起皱及孔洞先塌陷后膨胀。主文 p.4
• 作者保守地使用 acoustic waves，因为实验没有直接确认波的幅值和持续时间；后续可把实时被动声学或超高速 X 射线衍射与成像联用。补充材料 p.7-p.8
• 研究也解释了小球形孔的来源歧义：它们既可能继承自粉末，也可能由边缘稳定的熔化条件生成。主文 p.6

创新点

• 用原位高速 X 射线直接画出有粉与无粉条件下清晰、可复现的匙孔孔隙边界，并把粉末影响限制在可比较的工艺图中。主文 p.1、p.6
• 提出并实测匙孔深宽比相关的线性边界关系，用统一的匙孔孔隙常数连接高功率速度段与低功率速度段。主文 p.3-p.4
• 以兆赫兹成像捕捉匙孔塌陷、微射流、孔洞分裂与反弹、针状底部形成和孔洞迁移的连续事件链。主文 p.4-p.5；补充材料 p.16-p.30
• 把近边界孔隙形成分成声波驱动的主导机制和低功率波动主导的次级机制，并给出两者可由孔洞加速快慢区分的判据。主文 p.6；补充材料 p.7

背景：为什么做

• 金属增材制造希望通过局部组织控制获得按位置设计的性能，但 L-PBF 中常见的粗大柱状晶使这一目标难以实现。主文 p1-p2
• 柱状晶与等轴晶的转变受温度梯度和凝固界面速度的时空变化控制，光斑强度分布提供了除扫描策略以外的组织调节手段。主文 p2
• 实验试错难以直接捕捉粉末尺度熔池流动、蒸发、快速凝固和热循环共同作用下的三维晶粒演化，因此作者转向高保真热流体和晶粒生长耦合模拟。主文 p2

研究方法

模型

• ALE3D 求解三维瞬态质量、动量和能量守恒方程，并纳入 Marangoni 对流、蒸发反冲压力、蒸发冷却及熔池自由表面内的多次激光反射。主文 p2-p3
• 晶粒模型采用三维元胞自动机。每个晶粒用沿 FCC 晶体 <100> 方向生长的八面体包络表示，模型显式处理熔化、Poisson 随机形核、枝晶尖端生长和相邻元胞捕获。主文 p3-p4
• 元胞自动机的温度驱动力来自 ALE3D 的瞬态温度场，形核密度服从以过冷度为自变量的高斯分布，生长速度采用过冷度的幂律。主文 p3

数值方法

• ALE3D 在非结构网格上采用混合有限元和有限体积格式，并用完整光线追迹描述随熔池表面变化的激光能量沉积。主文 p2-p3
• CA 计算域为 0.86 mm × 0.30 mm × 0.30 mm，规则网格为 430 × 150 × 150，立方元胞边长 2 μm，共 9,675,000 个元胞。主文 p4；总数按三轴格点数相乘（430 × 150 × 150 = 9 675 000，导师复算）
• 作者采用自适应时间步，限制枝晶生长前沿单步推进距离不超过一个元胞；数值时间步长和总步数没有给出。主文 p4
• DREAM.3D 根据 EBSD 统计生成等效基板晶粒，并对 CA 输出的连续切片进行重构、分割和统计；晶体学织构的 ODF 与织构指数由 MTEX 计算。主文 p4；SI p1-p2

实验

• 作者用 EBSD 表征 316L 不锈钢基板的晶粒尺寸分布和晶体学织构，作为 DREAM.3D 合成初始组织的输入。实测等效球直径的拟合均值为 12.37 μm，标准差为 2.22 μm。主文 p4；SI p1
• 论文将模拟趋势与文献 [12] 的单道实验横截面对照。该实验中 CG、LE、TE 光斑对应的等轴晶面积比例分别为 2%、28%、77%。主文 p19
• 本文没有报告新开展的完整多道或多层打印实验，文中的光斑形状实验对照来自文献 [12]。主文 p19、p22 参考文献 [12]

结果：算了什么、做了什么实验

• 在相同功率和扫描速度下，TE 光斑形成更短、更宽、更浅的熔池，LE 光斑形成更长、更窄、更深的熔池，CG 介于二者之间。主文 p8
• 300 W、1.8 m/s 条件下，CG、TE、LE 在开光阶段形成的晶核数分别为 37、88、29；计入关光后阶段后，总数为 65、99、64。主文 p10、p16、p19
• 350 W、0.75 m/s 条件下，CG、TE、LE 在开光阶段形成的晶核数分别为 901、1355、822，TE 仍最早开始形核且形核增长最快。主文 p18-p19
• 不开启形核时，CG、TE、LE 的织构指数为 0.421、0.307、0.549；开启形核后变为 0.407、0.292、0.478，LE 的柱状晶织构最强，TE 最弱。主文 p10
• 关光后趋势发生反转。CG 和 LE 的熔池更深，冷却至形核过冷度所需时间更长，因此关光后形成的晶核多于 TE；第二组参数下 CG、TE、LE 的熔池深度分别为 54、36、80 μm。主文 p18、p20
• 二维切面上看似等轴的晶粒可能是从切面外延伸进来的柱状晶。三维模拟因此指出，只靠单个二维截面会高估等轴晶比例。主文 p9-p10、p20-p21

讨论：围绕什么展开

• 作者用 Hunt 柱状晶到等轴晶转变图解释形核差异，将局部温度梯度 G 和固液界面速度 R 投影到柱状、混合和等轴区域。主文 p12、p14
• 在 83 μs 时，CG、TE、LE 落入混合区的 G-R 点比例为 0.0%、6.9%、0.5%；在 243 μs 时为 2.3%、5.6%、1.9%，支持 TE 更早、更频繁形核。主文 p14-p15
• 开光阶段的控制量是熔池尾部过冷区宽度。TE 的尾部温度梯度更低且界面速度更高，因而具有更宽的过冷区；关光后没有持续热输入，熔池深度转而决定可积累的过冷度。主文 p16、p18
• 作者人为采用很大的形核密度以促成足够的形核事件。第一组和第二组参数分别使用 10^15 m^-3 和 10^13 m^-3，因此晶核数主要用于比较光斑形状，不是对实验绝对晶核数的预测。主文 p20-p21
• 模型只计算短单道，作者明确说明这是高保真模拟成本所限；Hunt 模型又假定固化区内温度梯度恒定，所以结论的外推对象主要是几何机制和相对趋势。主文 p21

创新点

• 将完整激光光线追迹和粉末尺度瞬态熔池流动与三维元胞自动机耦合，使任意光斑形状下的热历史可直接驱动晶粒形核和生长。主文 p2-p3
• 把光斑形状对组织的作用拆成开光与关光两个机制：开光时形核倾向随熔池宽度增加，关光后随熔池深度增加。主文 p18、p21
• 用三维晶粒连通性揭示二维截面中的伪等轴形貌，为实验 EBSD 截面的解释给出具体警示。主文 p9-p10、p20-p21
• 论文不仅比较稳态附近的三种光斑，还分析启停瞬态、另一组功率和速度、形核过冷度以及两类基板组织扰动。主文 p14、p18-p19；SI p10-p11

背景：为什么做

• 粉床选区激光熔化会出现孔隙、开裂、残余应力、变形、飞溅和表面粗糙等缺陷，实验参数研究通常依赖试错，成本高且不够系统。主文 p1-p2
• 该过程同时包含多相流、颗粒运动、传热和热致相变，因而难以用单一尺度的数值方法完整描述。主文 p2
• 连续介质 CFD 难以直接解析颗粒运动，未解析 CFD-DEM 的网格通常大于颗粒，也不能解析颗粒周围流场。主文 p2-p4
• 作者的目标是建立可解析粉末、熔池及其相互作用的方法，同时覆盖重复铺粉和熔化过程。主文 p2-p3

研究方法

模型

• 铺粉阶段使用未解析 CFD-DEM，熔化阶段使用半耦合解析 CFD-DEM，两阶段分开求解后组成完整流程。主文 p3-p5
• 浸入边界法解析固体颗粒周围流体；颗粒所占位置另建虚拟 CFD 域，使颗粒内部温度梯度和局部熔化也由多相 CFD 求解。主文 p5-p6
• 动量方程纳入 Darcy 阻力、表面张力、反冲压力和 Marangoni 流，能量方程纳入传导、耗散、熔化、对流、辐射和蒸发潜热项。主文 p6-p7
• 激光吸收采用 Beer-Lambert 形式，并按粉末未熔和熔池形成后的不同吸收状态切换衰减系数。主文 p7-p8

数值方法

• 模型在 OpenFOAM 与 LIGGGHTS 中实现，以 CFDEM 作为耦合引擎。主文 p8
• 自由表面用 iso-Advector 重构，速度和压力用 PISO-IB 迭代耦合。主文 p6、p8-p9
• 材料为 Ti-6Al-4V；正文把粉末尺度写为9至28 µm，Fig. 1色标则标为Particle Radius，随后给出平均直径约36 µm；基板厚80 µm、粉层厚100 µm。主文 p9-p11；半径与直径的表述口径不完全清楚
• CFD时间步为1×10^-7 s，DEM时间步为1×10^-8 s；主网格收敛性比较5、4、3.5 µm，后续采用4 µm。主文 p10、p13

实验

• 作者没有报告新做的实验，验证使用已发表文献中的原位观测、表面形貌、粗糙度趋势、高速成像和SEM图像。主文 p13-p25
• 单轨验证改变激光功率与能量分布，多轨验证改变道间距，并将数值现象与文献实验作定性或趋势比较。主文 p13-p16

结果：算了什么、做了什么实验

• 5、4、3.5 µm这3种网格所得熔道轮廓相近；5与4 µm轮廓最大差约4 µm，温度最大差约130 K，作者随后选4 µm。主文 p11-p13
• 在60 cm/s扫描速度下，150、200、300 W算例再现了随功率变化的波纹与堆积道头；较低功率还可能因熔化不足产生内部孔隙。主文 p13-p15、p22
• 200 W时，均匀能量分布形成的熔道比Gaussian分布更平、更规则，平均粗糙度也更小。主文 p14、p17-p18
• 模型捕捉到颗粒落入熔池或被扫离熔池的运动，但作者说明蒸气可能是颗粒运动的主导因素，本研究没有考虑它（模型本身含环境气体）。主文 p14-p15、p19
• 道间距从1.25倍激光直径减到0.75倍时，粘附颗粒和不连续逐步减少；25%搭接同时减弱这两类缺陷。主文 p10、p16、p21
• 多层模型再现了边缘粗糙和波纹，作者把粗糙归因于粘附粉末、部分熔化粉末和球化造成的不连续。主文 p17-p19、p25

讨论：围绕什么展开

• 作者围绕激光功率引起的熔池尺寸和温度梯度变化，解释波纹密度、尺寸和表面质量的变化。主文 p13-p14
• Gaussian能量集中在光斑中心，造成更高温度梯度和更尖锐波纹；均匀分布则降低熔道起伏。主文 p14
• 多轨部分把球化与台阶边缘不连续联系起来，并讨论用轨道搭接改善表面质量。主文 p15-p16
• 对于孔隙植入体等工作区很小的结构，轨道搭接可能受空间限制，缺陷不能简单靠同一策略消除。主文 p17-p19
• 论文把未考虑高速金属蒸气和计算成本列为两项主要缺点；加入蒸气可能还要缩小网格。主文 p22-p24

创新点

• 半耦合做法用解析 CFD-DEM求颗粒受力，同时把颗粒位置映射为虚拟 CFD 相来求内部温度梯度和部分熔化，绕开一个DEM颗粒只有一个温度的限制。主文 p5、p20
• 作者把解析熔化与未解析铺粉接成统一流程，使模型能够循环处理熔化、凝固表面替换、平衡和滚筒铺粉。主文 p16-p17、p19-p22
• 模型把iso-Advector、浸入边界法和温度相关激光吸收模型组合起来，同时处理自由表面、颗粒-流体作用和激光穿透。主文 p6-p8、p20
• 作者称，熔化时粉末运动此前尚未由已有研究尝试，而该耦合模型能够直接捕捉这种运动。主文 p14

背景：为什么做

• L-PBF中的复杂热史、大温度梯度和高冷却速率会产生高位错密度及特殊位错结构，但这些结构的起源以及后续受载时的稳定性仍有争议。主文 p1
• 显微方法通常只能看到制造后的最终位错，制样还可能改变结构；原位电子显微难以用于受激光照射的L-PBF样品，同步辐射XRD也不能分辨单个位错结构。主文 p2
• 作者集中回答热应力是否足以形成高密度位错结构、结构能否稳定、结构是否只在特定位置晶粒形成，以及多层制造温度场怎样影响位错演化。主文 p2
• 分子动力学和离散位错动力学计算代价高且时间尺度有限；作者因此采用能描述位错曲率与胞状图样的高维连续体位错动力学。主文 p2

研究方法

模型

• 模型把有限应变热-力晶体塑性与温度相关CDD完全耦合；每个滑移系用总位错密度、刃型与螺型几何必要位错密度、总曲率密度四个状态变量。主文 p3-p4
• 塑性滑移服从Orowan关系，位错速度由分解剪应力与临界分解剪应力控制；临界应力和弹性常数随温度变化。主文 p3-p4
• 模型加入总位错密度的统计湮灭项，湮灭率与位错速度绝对值、总位错密度平方及特征湮灭距离成正比。主文 p6-p7
• 当前CDD只处理单一活跃滑移系，不含多滑移系相互作用；位错增殖只由曲线位错长度增长表示。主文 p3-p4
• 多尺度链条由热流动温度场、相场晶粒结构、现象学多晶CPFEM和晶粒内CDD组成；本文把既有大尺度结果中的位移、投影剪切应变与平均温度传给小尺度模型。主文 p7-p9

数值方法

• 模型在有限元框架MOOSE中实现；应力平衡用标准连续Galerkin法，CDD状态变量用不连续Galerkin法和隐式时间积分。主文 p3、p5
• 界面通量把位错密度分解成正/负、刃/螺分量，并按迎风方向选取通量，因此能处理位错偶极等高度不连续密度场。主文 p2、p5-p6
• 外边界对位错通量使用周期边界条件，四个CDD方程都加入很小的人工扩散项用于数值稳定。主文 p6
• 主薄滑移面域为40 × 40 × 0.4 μm，完整晶粒域为8 × 8 × 12 μm，单元边长均为0.4 μm；对应单元数分别为10000和12000。主文 p9；单元数按域尺寸除以单元边长算（导师复算）

实验

• 本文没有开展新的制造或显微实验；L-PBF工艺参数取自Wang等的316L实验，晶粒取向取自同材料的EBSD实验。主文 p7
• 作为CDD输入的多晶CPFEM此前已通过特定取向晶粒族的模拟晶格应变与XRD测量对比进行标定。主文 p7
• 本文的验证使用已发表的TEM密度区间、DDD结果和文献显微图像；作者明确称直接、全面的实验比较仍受数据不足限制。主文 p15-p17、p21-p22

结果：算了什么、做了什么实验

• Grain 1在加热后形成明显胞状位错结构，约从670 μs稳定到1170 μs；高、低密度区分别达到约4.5和2.5 μm^-2。主文 p9-p10
• Grain 1最终高低密度结构的特征周期约12.6 μm，胞状高密度区平均直径约6.3 μm。主文 p10-p11
• 不同位置晶粒给出不同终态：熔池外且经历非单调热载荷的Grain 1和4结构更清楚，凝固后单调受载的Grain 2和3最终位错密度更高。主文 p11-p12
• 增大湮灭距离会降低总体位错密度，拉大高低密度区之比，并形成密度接近零的小胞区；减小湮灭距离则总体密度升高，低密度区更连通，结构的特征尺度增大。主文 p12-p13、p15-p17
• 冷却后的单滑移结构在分解剪应力约16.5 MPa以内大体稳定，超过后会继续演化以适应更大的塑性变形。主文 p13
• 4层循环使胞状结构逐轮变大；含湮灭的模拟形成了包围胞区的低密度区，且平均位错密度在首次凝固冷却后最高。主文 p13-p15
• FFT给出单层湮灭算例的特征峰约6.7-10 μm，4层循环算例约12 μm；作者据此认为循环热载荷可把特征尺度放大两倍或更多。主文 p15-p16

讨论：围绕什么展开

• 即使不引入微偏析和析出，热膨胀、收缩所致的塑性变形以及CDD中的输运、增殖机制也能形成位错胞结构，支持热应力是其来源之一。主文 p15
• 熔池内晶粒在凝固后的热应力下积累更多位错，顶部自由表面附近的首次热循环增量最大；后续循环不一定继续增加密度，湮灭可使其下降。主文 p15
• 湮灭机制控制密度水平、高低密度对比和结构尺度，但当前用单一中间湮灭距离代替刃位错与螺位错各自的距离，因此仍需分别建模。主文 p15-p17
• 单元边长0.4 μm不能分辨更小结构，且本文只分析一组工艺参数；与文献实验的尺度比较属于量级验证。主文 p16-p17
• 模型假定凝固后的初始状态，最大剪切应变约4.7%；位错密度定义在参考构形中，较大变形时还需重写输运方程。主文 p17
• 单滑移模型不能直接重现实验中的多滑移位错网络；作者指出多滑移扩展同时受计算能力和矢量位错密度多极展开策略限制。主文 p17

创新点

• 提出温度相关的四状态变量CDD，并与热-力晶体塑性完全耦合，用于L-PBF热史下的位错输运。主文 p1、p3-p4
• 构造连续Galerkin与不连续Galerkin混合求解，并用正/负、刃/螺分解定义守恒界面通量，以处理强不连续位错密度场。主文 p2、p5-p6
• 把多晶L-PBF计算得到的晶粒级真实温度和变形历史传给单晶CDD，比较熔池内外不同位置晶粒及多层循环。主文 p2、p7-p9、p13-p14
• 在同一框架中分析湮灭距离、冷却后机械加载和多层循环对结构稳定性及特征尺度的影响。主文 p12-p16

背景：为什么做

• 增材制造中的非均匀热循环会产生热应力，热裂纹是常见且破坏性很强的缺陷；其形成跨越多物理过程、时间尺度和长度尺度，数值研究困难。主文 p1-p2
• 已有工作多从显微组织表征或多尺度建模研究热裂，但作者认为直接面向热裂纹的多物理数值模型仍很有限，机制也未充分澄清。主文 p2
• 作者指出，完全耦合的内聚相场热断裂模型在热力学一致性、裂纹对弹性能和导热的反馈、断裂性质随温度变化等方面仍缺少系统研究，且此前尚未用于增材制造热裂纹。主文 p3

研究方法

模型

• 模型采用小变形、弹性静力和准静态脆性断裂假设，主变量为位移、损伤相场和温度。自由能分为弹性、断裂和热能三部分，并由 Clausius-Duhem 不等式推导本构与耦合关系。主文 p3-p4
• 弹性应变能只让拉伸部分随损伤退化；裂纹不可逆性用历史变量和 KKT 条件约束。作者比较 AT2 与内聚相场 CPF，并强调 CPF 的损伤起始阈值和对相场长度尺度不敏感。主文 p5-p7
• 断裂韧度 Gc 写成温度函数；导热率按损伤退化，使完全开裂区不再出现非物理的穿裂纹导热。当前实现忽略应变率和损伤率对应的瞬态耦合项。主文 p7-p8
• PBF 热裂模拟除热应变外还加入凝固收缩，并用有效热膨胀系数表示固相线到液相线温区内的相变收缩。主文 p15-p16

数值方法

• 控制方程用二维四节点四边形有限元离散；每个时间步先单体求解温度和位移，再固定更新后的温度、位移求损伤，循环至全局残量变化低于容差。主文 p9-p10
• 代码基于 MOOSE，使用自动微分生成 Jacobian，并用 MultiApp 在主应用与损伤子应用之间传递场变量。主文 p10
• 验证算例包括单边缺口拉伸和陶瓷板淬冷；前者检查 CPF 的长度尺度敏感性与裂纹导热退化，后者与已发表的淬冷实验图样比较。主文 p10-p15
• PBF 热裂先使用椭圆熔池的解析温度场，再使用三维现象学 PBF 热模型求得的稳态截面温度和相态作为二维平面应变热裂计算的初始条件。主文 p16-p20

实验

• 本文未报告作者新做的物理实验。淬冷验证使用 Jiang 等人的实验图像和裂纹统计，PBF 热裂比较使用 Stopyra 等人的单道实验截面图。主文 p14-p15、p18
• 实验对比用于检查裂纹形貌和趋势；本文的 PBF 参数研究被作者定位为定性检查，不是材料或工艺参数的定量标定。主文 p21、附录B p24

结果：算了什么、做了什么实验

• 单边缺口算例中，CPF 的全局力学响应对长度尺度几乎不敏感，而 AT2 的峰值载荷随长度尺度增加而下降；CPF 的损伤带也更紧凑。主文 p11
• 裂纹使导热率退化后，温度在裂纹两侧出现跳变，热流绕过裂尖并在裂尖集中；若保持导热率不变，热流仍穿过完全开裂区。主文 p12-p14
• 淬冷算例再现了裂纹起始、竞争生长、停裂和层级形貌；作者报告其裂纹数量、间距及图样与引用实验相符。主文 p14-p15
• 解析温度场算例显示，凝固收缩主要形成沿固液相变区的环向裂纹；提高人为温度梯度参数后才出现较弥散的中心损伤，暴露了解析温度场的局限。主文 p16-p19
• 数值 PBF 温度场下，提高激光功率或降低扫描速度会增强中心裂纹并使形貌转向 keyhole 模式；降低功率或提高速度则更接近只有环向裂纹的 conduction 模式。主文 p20-p21
• 附录把弹性模量改为温度相关后，具体裂纹结果发生定量变化，但功率增大时由 conduction 模式转向 keyhole 模式的趋势不变。附录B p23-p24

讨论：围绕什么展开

• 作者围绕热、力、裂纹三场如何互相反馈展开讨论，重点是温度依赖断裂韧度、裂纹引起的导热退化和凝固收缩应变。主文 p3、p7-p8、p15-p17
• 对热裂形貌的解释分成两条机制：凝固收缩主导环向裂纹，熔池附近的高温度梯度主导中心裂纹；功率和速度通过线能量密度改变温度梯度。主文 p17、p20-p21
• 作者明确承认当前弹性脆性框架在温度超过液相线时并不成立，常数材料性质只是初步近似；后续需要加入塑性、延性断裂及更准确的相演化。主文 p3-p4、p21、附录B p23
• 解析椭圆温度场不能准确重现中心裂纹，现象学 PBF 温度场虽改善形貌，但网格、时间离散和完整热历史未在本文交代，其他参数还转引到 Yang 等人的论文。主文 p17-p18、p21

创新点

• 作者给出一个由热力学基本原理推导的热力耦合相场断裂框架，同时包含损伤对弹性能和导热率的反馈以及断裂韧度的温度依赖。主文 p3-p8、p21
• 文章把长度尺度不敏感的 CPF 放入多物理热断裂问题，并通过单边缺口算例与 AT2 比较。主文 p5-p7、p10-p11
• 按作者的表述，这是将该相场热断裂方法用于增材制造热裂纹的首次尝试；模型把凝固收缩、解析熔池温度场和数值 PBF 温度场纳入同一套热裂分析。主文 p3、p15-p21
• 参数研究不只给出裂纹形貌变化，还把环向裂纹和中心裂纹分别关联到凝固收缩与温度梯度，为 conduction 和 keyhole 两类热裂模式提供了机制解释。主文 p17、p20-p21