材料与方法
研究采用牌号为PNK UT3的镍粉（平均粒径约20 µm，纯度99.9%）作为基体材料，以及平均尺寸约40 nm的碳化硅（SiC）纳米颗粒作为改性相。粉末制备过程包括在行星式球磨机中进行机械活化与混合，球料比为5:1，并在氩气保护下进行。为减少铁污染，镍粉预先与钢球进行20小时的“擦拭”处理。随后，将镍粉与预分散的SiC纳米颗粒（含量分别为0.001、0.005和0.015 wt.%）混合20分钟。混合后的粉末在异丙醇（粉液比1:5）和油酸存在下进行分散处理，并辅以超声处理。最终悬浮液在通风橱中室温干燥24小时。
材料 consolidation 通过放电等离子烧结（SPS）在石墨模具中进行，烧结温度为850°C，保温时间20分钟，压力50 MPa，采用直流脉冲模式（脉冲12 ms，间歇2 ms）。烧结后样品在空气中冷却。
材料表征包括通过静水力学称重法测定密度，使用扫描电子显微镜（SEM）和透射电子显微镜（TEM）研究微观结构，通过X射线衍射（XRD）分析相组成。力学性能测试包括维氏显微硬度测量（负载200 gf，保压10秒）、三点弯曲试验（分别在25°C和400°C下进行）以及通过超声波法测定杨氏模量（温度范围25–750°C）。此外，还利用热力学软件包（TERRA, MAAT, IVTANTERMO）和OQMD数据库进行了热力学与动力学建模，计算了Ni-Si-C体系的相图、化学势和碳在镍中的扩散系数。
结果
烧结镍的结构与相组成
TEM分析显示，纯镍烧结后形成超细晶结构，晶粒呈伸长状，纵向尺寸大于0.4 µm，横向尺寸约0.15 µm，表明在非平衡SPS条件下发生了定向再结晶。在晶界处发现了热解碳区域。X射线荧光分析表明，表层碳含量达到1-2 wt.%，而体相碳含量不超过约1 wt.%，显示出碳扩散饱和的梯度特性。XRD图谱证实了面心立方（FCC）镍晶格的保留，晶格参数为3.526–3.527 Å，略高于纯镍的标称值（3.524 Å），表明间隙碳原子存在于固溶体中。根据峰宽估算的相干散射区域尺寸约为390 nm，与亚晶结构吻合。在小角度区域（2θ~10–30°）检测到漫散晕，这是存在X射线非晶相的特征。通过与文献中镍晶格参数随碳浓度变化关系的比较，估算出固溶体中碳浓度约为0.1–0.3 at.%，这意味着相当一部分碳（可达几个原子百分比）以游离碳相（热解碳、石墨）形式存在。热解碳的形成可能与冷却过程中固溶体析出过量碳、石墨工具和垫纸提供的碳源以及油酸的热分解有关。
Ni-SiC复合材料的结构
SEM图像显示了Ni-0.001SiC和Ni-0.005SiC复合材料的微观结构。截面上观察到亚微米和微米级孔隙，晶粒尺寸可达20 µm。经蚀刻后，在Ni-0.005SiC样品中揭示了尺寸达0.6 µm的亚晶，表明存在发达的亚结构，具有Hall-Petch强化的潜力。
热力学建模与化学势分析
利用OQMD数据库计算的Ni-Si-C三元相图在0 K时显示存在SiC、Si2Ni、NiSi、Ni2Si、Ni3Si2、Ni31Si12、Ni3Si和NiSi2等稳定相。在1123 K下，若不考虑固溶体，镍角相图包含纯Ni、石墨和硅化物Ni2Si。若考虑Ni-C和Ni-Si理想固溶体的形成，则会出现宽广的Ni-C固溶体区域，与Ni2Si共存，而Ni-Si固溶体区域极窄。对Ni-C-SiC-NiO-SiO2体系平衡组成的计算表明，随温度升高，会形成NiSi和Ni2Si硅化物。化学势计算揭示了相稳定区对温度的临界依赖性。在300 K时，各相（SiC, Si, Ni, 硅化物）的稳定区边界清晰，SiC区域广阔，表明其在室温下热力学稳定性高。当温度升至1273 K时，声子熵的贡献显著降低了金属和金属间化合物的自由能，导致镍硅化物（尤其是富镍相Ni3Si）的稳定区显著扩大，同时SiC的稳定区变窄。这从热力学上预测了SiC在镍存在下于烧结温度分解的趋势。组分活度分析表明，在Ni-C体系中，镍的活度随碳含量增加先降后升；在Ni-Si-C三元体系中，碳和硅的活度在较宽浓度范围内显著大于1，表明该温度下组分间相互作用弱，硅碳原子间存在显著排斥力。Ni-Si体系中镍的活度随硅浓度增加逐渐降低，表明镍硅原子在固溶体中结合强。Ni-Si-C体系固溶体形成吉布斯自由能的计算显示，在低C和Si含量时，吉布斯自由能等于或大于0；当碳浓度轻微增加（与本工作中报道的相当）时，吉布斯自由能略低于0，这与文献中关于850°C时SiC与镍相互作用分解出游离纳米碳并形成镍硅化物的结果一致。
碳在镍中的扩散动力学
利用Arrhenius方程计算碳在镍中的扩散系数表明，在850°C（1123 K）的烧结温度下，扩散系数（D）约为4 × 10−12m2/s。该计算结果与独立实验数据吻合，表明在SPS过程中，碳从石墨工具可在大约5分钟内渗透入镍深度达100 µm。然而，与量子化学计算（DFT）数据的直接比较揭示了一些差异。DFT计算得出的碳在镍体相中迁移的激活能约为1.62 eV（≈156 kJ/mol），高于本计算所用的值（138.3 kJ/mol）。这种差异可通过实际材料中存在的结构缺陷（晶界、位错）降低有效扩散能垒来解释。晶界扩散的激活能可能显著更低。因此，实验中观察到的高碳饱和速率可能归因于SPS获得的超细晶结构中短路径扩散（晶界扩散）占主导地位。Ni-C体系理想吉布斯自由能、混合吉布斯自由能和混合焓随碳浓度变化的计算表明，混合吉布斯自由能低于理想吉布斯自由能，且随碳浓度增加而降低，表明随温度升高，形成碳在镍中固溶体的可能性增加。对含碳体系（无论碳以游离态还是固溶体形式存在）在298–1198 K温度范围内平衡态的计算结果显示，CO2的形成始于598 K，并在798 K以下活跃进行，之后CO的形成开始占主导，同时伴随NiO的还原和Ni3C的非剧烈形成。
力学性能
含纳米颗粒样品的开孔率平均约为1.3%。力学性能测试结果表明，即使超低含量的SiC纳米颗粒也能显著提高抗弯强度。所有复合材料的相对密度平均超过98.5%，显示了高的烧结质量。实验数据清晰地展示了超低浓度SiC纳米颗粒和SPS方法对镍显微硬度的复杂影响。关键发现是在最小分散相浓度下观察到极端硬化现象，这无法用经典模型解释。具体而言，引入仅0.001 wt.% SiC可使硬度从70 ± 10 HV显著增加至81 ± 3 HV。将添加剂含量增至0.005 wt.%时，硬度进一步达到最大值92 ± 4 HV。然而，将含量增至0.015 wt.%时，硬度下降至82 ± 4 HV，表明最佳值在0.005 wt.%附近。这种性能的极端变化是协同硬化机制的特征。添加0.001 wt.% SiC使20°C下的抗弯强度提升至1130 MPa（相对于纯镍+115%），400°C下提升至976 MPa（+86%）。Ni–0.001SiC复合材料的抗弯强度值（1130 MPa）是超细晶纯镍文献中已知的最高值之一，超过了例如通过剧烈塑性变形方法获得的值，强调了本工作实现的协同强化机制的有效性。随着SiC含量进一步增加，强度有所下降，在0.001 wt.%处形成极值。纯镍的断口呈韧性韧窝状，而含SiC纳米颗粒的样品在弯曲过程中表现出显著的塑性变形后才断裂。
讨论
观察到的镍在0.001 wt.% SiC下的极端强化不能仅用经典的Orowan机制解释，因为颗粒的体积分数极小，计算贡献ΔσOr仅为1–5 MPa。这种在创纪录低分散相含量下观察到的强化极值与现代概念一致，即超低添加物在结构缺陷处偏聚并根本影响其迁移性的关键作用。在本研究中，来自分解的SiC纳米颗粒的碳就是这样的元素。主要强化因素包括：
Hall-Petch强化：源于超细晶亚结构。对于平均晶粒尺寸d ≈ 0.25 µm和Hall-Petch常数kHP≈ 0.1–0.2 MPa·m1/2，计算得到ΔσHP≈ 200–400 MPa。对于具有较大晶粒（d ~20 µm）的区域，贡献降至20–40 MPa，反映了结构的异质性。
固溶强化：由固溶体中的碳（0.1–0.3 at.%）贡献，估算ΔσSS≈ 45–135 MPa。
织构和伸长晶粒形态：SPS在压力下形成的织构和伸长晶粒形态产生额外贡献Δσtexture~50–150 MPa。
SiC分解产物的协同贡献：热力学建模和化学势分析明确证明了Ni-SiC体系在850°C烧结温度下的热力学不稳定性。SiC纳米颗粒的分解导致形成游离碳（晶界处的热解碳）以及假设的镍硅化物（Ni2Si和NiSi）。根据经典Zener理论，抑制晶粒长大的效率由颗粒的体积分数和尺寸共同决定。由SiC受控分解形成的纳米尺度硅化物和热解碳析出物恰好定位于晶界，是异常有效的晶界钉扎中心，这解释了在如此小的总添加量下UFG结构的稳定性。在超低SiC含量（0.001 wt.%）下，形成的硅化物中间层是不连续的，可作为位错的有效障碍，防止形成连续的脆性相。同时，释放的碳部分溶解于镍（固溶强化），部分形成纳米尺度的热解碳析出物，同样抑制位错运动和晶界迁移。
总计算强度增量约为：Δσtotal≈ ΔσHP+ ΔσSS+ Δσtexture≈ 295 至 685 MPa。
烧结温度下碳在镍中的高扩散速率，确保了碳沿晶界的快速重新分布，有助于稳定超细晶结构。因此，0.001 wt.% SiC处的强化极值是以下因素最佳平衡的结果：晶粒细化（Hall-Petch）和织构强化的最大贡献；显著的固溶强化；受控SiC分解产物（热解碳、硅化物）的稳定作用，且不形成大范围脆性相；碳的高扩散迁移性为这些过程提供了有效的动力学保障。当SiC含量增加（>0.001 wt.%）时，分解产物量增加，可能导致形成更连续且可能脆性的硅化物中间层，成为裂纹萌生中心，解释了强度下降的原因。
所获得的超细晶结构具有良好的短期热稳定性。根据Zener理论，位于晶界的纳米尺度热解碳和假设的硅化物颗粒是有效的晶界固定中心，应能防止晶粒合并。然而，镍中的碳固溶体是亚稳态的，热力学建模表明在600–800 K以上温度下，碳从晶格中析出在能量上更有利，因此长时间高温暴露可能导致由于固溶体分解和析出物聚集而引起的强化机制改变。
将本研究的抗弯强度结果（1130 MPa）与通过替代 consolidation 方法硬化的镍进行比较，凸显了具有短处理周期技术的优越性。例如，SPS等快速方法能够最小化晶粒生长，这对于实现Hall-Petch硬化机制至关重要。在本研究中，预机械活化与SPS的结合不仅保留了超细晶结构（~0.25 µm），还观察到了固溶硬化和分散析出物的附加贡献，这是SiC纳米颗粒受控分解的结果。这使得获得的抗弯强度与通过剧烈塑性变形方法获得的最佳值相当，但具有制造更复杂形状产品的潜力。
硬度对晶粒尺寸的依赖关系通常遵循Hall-Petch关系，但在超细晶尺度（d < 1 µm）可能会观察到硬化饱和，这与新变形机制的激活有关。在SPS材料中，除了晶粒尺寸，位错胞、变形孪晶等晶内亚结构也对硬化有关键影响。本研究的显微硬度数据与这些基本规律一致，并展示了晶粒尺寸和分散添加剂的影响。观察到的性能在特定分散相含量下的极值表明，需要仔细选择其浓度，以在晶界硬化和脆化风险之间取得最佳平衡。
结论
通过机械活化和在850°C石墨模具中的放电等离子烧结（SPS），成功制备了Ni–ySiC体系（y = 0.001, 0.005, 0.015 wt.%）复合材料，其力学性能得到改善。
对于Ni–0.001SiC复合材料，其在20°C的抗弯强度相比未添加剂的烧结镍提高了115%（达到1130 MPa），在400°C提高了86%（达到976 MPa）。在更高SiC含量下，强度下降，表明性能极值出现在0.001 wt.%纳米颗粒处。
微观结构分析揭示了烧结镍具有超细晶结构（0.15–0.4 µm），晶粒呈伸长状，晶界存在热解碳，固溶体中碳浓度约为0.1–0.3 at.%。相当一部分碳（约6–7 at.%）以游离相形式存在。
综合热力学（化学势）和动力学（扩散）建模证实了在SPS温度下SiC分解形成游离碳和镍硅化物（Ni2Si, NiSi）的可能性。研究表明，随着温度从300 K升高到1273 K，镍硅化物的稳定区扩大，而SiC的稳定区变窄，这从热力学上有利于界面反应。
0.001 wt.% SiC处的性能极值通过几种强化机制的协同作用得到解释：固溶强化（碳）、晶界强化（超细晶结构的Hall-Petch效应）、位错强化（高位错密度和伸长晶粒形状）以及受控SiC分解产物（热解碳、假设的硅化物）对结构的稳定作用，同时最小化了广泛脆性相的负面影响。
所获得的结果对于开发可用于航空航天工程和动力装置的高强度、热稳定镍基复合材料具有重要意义。
