组分比对PTFE-Al-W复合材料动态性能及介观力学响应的影响

《Materials & Design》：Effect of component ratios on dynamic properties and mesoscale mechanical response of PTFE-Al-W composites

　　

在动态压缩下，聚四氟乙烯（PTFE）基体承受大部分塑性变形。较低的PTFE分数在固定全局应变下扩大剪切失稳区，降低复合材料的宏观断裂韧性。开裂和断裂源于动态压缩过程中的剪切失稳（局部变形带）。初始孔隙的存在提供了额外的变形空间，有效延迟了剪切失稳的发生。

**论文解读文章**

**研究背景与问题**

冲击引发含能材料（reactive materials）因其独特的力-化耦合响应而备受关注，在极端动态载荷（如冲击或爆轰）下可快速发生放热反应，为提升常规弹药效能提供新途径。作为改性配方，PTFE-Al-W复合材料在基础PTFE-Al体系中引入钨粉，旨在提高密度以增强穿透能力。组分比、粒径及制备工艺的差异会显著改变其力学行为与化学响应，从而满足不同平台的应用需求。然而，现有研究对组分比如何影响力学性能与失效行为，以及宏观力学性能与微观/介观力学行为之间的关系仍不明确。为此，本文通过霍普金森压杆试验与介观数值模拟相结合的方法，系统研究了组分比对PTFE-Al-W复合材料动态性能的影响，旨在建立宏观与微观行为之间的关联，为材料优化设计提供指导。

**研究内容与结论**

研究人员制备了五种固定PTFE与Al质量比（约2.77:1，对应化学计量反应4Al+3C₂F₄→4AlF₃+6C）的PTFE-Al-W复合材料试样，钨含量从16.3 wt.%递增至84.6 wt.%，理论密度范围2.7–8.9 g/cm³，孔隙率相应变化。通过分离式霍普金森压杆（SHPB）进行动态压缩试验，结合高速相机记录破坏过程，并利用VHX-7000数字显微镜观察微观结构。同时，基于ANSYS参数化设计语言（APDL）建立二维介观模型，采用Johnson-Cook本构模型与Grüneisen状态方程描述材料，利用任意拉格朗日-欧拉（ALE）方法处理PTFE大变形，使用LS-DYNA求解器进行数值模拟，分析剪切失稳与裂纹演化。

研究得出以下结论：组分比显著影响屈服应力、极限强度与极限应变；PTFE基体作为连续网络承担大部分塑性变形；裂纹与断裂源于动态压缩下的剪切失稳；初始孔隙可容纳额外变形，延迟剪切失稳的发生。该研究为理解PTFE-Al-W复合材料宏观性能与介观行为的关系提供了重要依据，论文发表在《Materials》。

**主要技术方法**（不超过250字）

主要关键技术方法包括：1）采用分离式霍普金森压杆（SHPB）进行动态压缩试验，基于一维应力波理论计算应力、应变与应变率（公式(1)–(3)），并使用高速相机（帧率10⁵帧/秒）记录破坏过程；2）介观数值模拟，基于ANSYS平台建立二维模型，采用Johnson-Cook本构模型（公式(6)）和Grüneisen状态方程（公式(7)）描述PTFE、Al、W，利用任意拉格朗日-欧拉（ALE）方法处理PTFE大变形，采用自动单面接触算法，使用LS-DYNA求解器求解，并设置速度边界条件（15 m/s）对应应变率约3200 s^-1。金属颗粒与孔隙假设为圆形，随机分布由APDL生成。试样原料来源：PTFE与Al粉来自贵州贵天诚科技有限公司，W粉来自北京兴荣源科技有限公司。

**研究结果**

**3.1 动态性能（Dynamic properties）**
通过SHPB试验获得真应力-应变曲线。在约3200–3250 s^-1应变率下，D-1试样（14.4% PTFE）的屈服应力最大（34.54 MPa），B-1试样（32.7% PTFE）的极限强度最大（60.05 MPa）。随PTFE比例增加，屈服应力与极限强度先增后减；极限应变在PTFE比例约24%后不再增加。随钨含量增加，极限应变在约60% W后快速下降，表明断裂韧性降低。破坏形态观察：A-1与B-1主要塑性变形，B-1边缘有小裂纹；C-1严重碎裂；D-1与E-1粉碎。断裂韧性随PTFE比例降低而逐渐下降。

**3.2 压缩试样微观结构（Microstructure of the compressed specimens）**
利用VHX-7000数字显微镜观察，A-1表面无裂纹，A-2与B-1表面有大量裂纹；严重破碎的试样（C、D、E）存在贯穿厚度裂纹。更高分辨率下，A-1中PTFE基体出现微裂纹，A-2中裂纹相交，B-2中出现可能引发断裂的宏观裂纹。

**4.1 材料模型与颗粒初始分布（Material model and initial distribution of particles）**
采用Johnson-Cook模型与Grüneisen状态方程描述PTFE（ρ=2.2 g/cm³，A=11 MPa，B=44 MPa，N=0.121，C=1.7×10³，m=1.12，c₀=1.59×10³ m/s，s=1.59，γ₀=0.59）、Al（ρ=2.7 g/cm³，A=265 MPa，B=426 MPa，N=0.34，C=0.015，m=1.295，c₀=5.4×10³ m/s，s=1.34，γ₀=1.97）与W（ρ=19.35 g/cm³，A=1510 MPa，B=117 MPa，N=0.12，C=0.016，m=1.0，c₀=4.0×10³ m/s，s=1.23，γ₀=1.67）。模型中颗粒与孔隙假设为圆形，使用APDL生成随机分布，网格密度收敛性验证确保PTFE失效体积变化<2%。

**4.2 模型有效性分析（Model validity analysis）**
建立与A、B、C三种试样对应的二维模型（孔隙率与组分比对应）。当全局应变达0.29时，模拟的剪切应变分布与实验微观结构吻合：A型仅局部剪切失稳，结构完整；B型剪切失稳区增多并有局部裂纹；C型大量裂纹，接近严重碎裂。验证了模型可有效反映裂纹形态随组分比与孔隙率的变化。D、E因高颗粒含量建模困难未比较。

**4.3 组分比的影响（Effect of component ratios）**
固定孔隙率3%，建立五种体积分数不同的模型（PTFE 75.5%–46.6%）。全局应变0.375时，A型仅少量裂纹，B、C、D型多裂纹，E型完全断裂。剪切应变分布显示，随PTFE比例降低，剪切失稳区扩大。PTFE失效比例历史曲线表明，在大部分压缩过程中，钨含量越高，PTFE失效比例越高，且该比例受全局变形影响。力链分析显示，B型中力链比A型更明显，但随金属比例增加，力链增强而极限强度降低，表明存在最优组分可平衡剪切失稳与强化效果。

**4.4 初始孔隙率的影响（Effect of initial porosity）**
固定PTFE/Al/W体积分数58%/17%/25%，改变孔隙率（1%–9%）。PTFE失效比例-全局应变曲线显示，低应变（<0.22）时各曲线重合，随后孔隙率1%的试样在0.22–0.36应变区间失效比例最高，孔隙率9%的试样在大部分压缩过程中失效比例最小。剪切应变分布表明，孔隙率1%时剪切失效最严重，大孔隙率下剪切失稳位置与路径改变，初始孔隙通过提供额外变形空间延迟了剪切失稳的发生。

**总结讨论与结论**

研究人员指出，当前工作未考虑化学反应效应，实验与数值分析仅限于力学加载与微观结构表征；实验样本数量有限，结果应视为初步证据，未来需更大样本集验证通用性。结论总结如下：
(1) 组分比显著影响屈服应力、极限强度与极限应变。最高屈服应力（34.54 MPa，约3200 s^-1）与极限强度（60.05 MPa）分别由不同配方实现（14.4/5.2/80.4 wt.% 与 32.7/11.8/55.5 wt.%）。
(2) PTFE基体作为嵌入金属颗粒的连续网络，承担大部分塑性变形，传播裂纹的汇合导致试样破裂。
(3) PTFE-Al-W试样的开裂与断裂源于动态压缩下的剪切失稳。剪切带的形成时间、位置与传播路径受组分比与初始孔隙率显著影响。
(4) 随PTFE比例降低，恒定全局应变下的剪切失稳区扩大，导致宏观断裂韧性降低。此外，初始孔隙的存在提供了更多变形空间，从而延迟了剪切失稳的发生。