大型碳纤维部件VARTM工艺仿真：流动建模与沟流预防

大型CFRP部件的真空辅助树脂传递模塑计算仿真。使用PAM-RTM和RTM-Worx的流动前沿建模、渗透率表征、铺覆仿真以及风电叶片和船舶结构的沟流缓解策略。

大型VARTM制造的挑战

真空辅助树脂传递模塑（VARTM）是大型碳纤维增强聚合物部件——超过80米的风电叶片、船舶壳体段和航空结构工装——的主导制造工艺。与热压罐工艺不同，VARTM依靠真空驱动压差（通常95-98 kPa）将树脂注入干的纤维预成型体。在这些尺度上，即使是纤维渗透率、预成型体压实度或树脂粘度的微小变化也可能产生代价高昂的缺陷：干斑、孔隙和富树脂区，它们会损害力学性能。

挑战因涉及复杂曲率、厚度过渡和夹芯嵌件的零件几何形状而加剧。一个典型的60米风电叶片需要约2.5-3.5吨干碳纤维织物，按40-60层离散铺层排列，具有不同的纤维取向。没有精确的工艺仿真，制造商将面临试错循环，每个迭代仅材料浪费就达5-15万美元。

VARTM流动建模的基础物理

通过纤维多孔介质的树脂流动遵循达西定律，各向异性情况由渗透率张量K描述：

达西速度：u = —(K/μ) · ∇P

其中μ是树脂粘度（环氧体系在注射温度下通常为100-800 mPa·s），∇P是压力梯度。典型双轴和三轴碳纤维无卷曲织物的面内渗透率在主纤维方向（0°）为1×10⁻¹⁰至8×10⁻¹⁰ m²，横向（90°）为0.3×10⁻¹⁰至2×10⁻¹⁰ m²。厚度方向渗透率显著较低，为0.05×10⁻¹⁰至0.3×10⁻¹⁰ m²，对于长宽比超过100:1的薄壳几何体，可近似为二维流动前沿。

仿真软件对比

软件	流动求解器	渗透率输入	铺覆仿真	热耦合	年度许可费
PAM-RTM (ESI)	有限元+控制体积	3D张量K	集成	是（放热固化）	$35,000-55,000
RTM-Worx (Polyworx)	有限差分+VOF	2D/3D K场	外部（CATIA）	是	$25,000-40,000
ANSYS Fluent + 多孔介质	有限体积	自定义UDF	未集成	是（完整CFD）	$40,000-70,000
OpenFOAM（开源）	有限体积+VOF	自定义UDF	未集成	自定义	$0
LIMS (UC Davis)	有限元+控制体积	2D/3D K	未集成	有限	$2,500（学术）

对于大型工业VARTM仿真，PAM-RTM和RTM-Worx占主导地位，因为它们具有专用的预成型体渗透率库、铺覆模块和经验证的沟流模型。两者之间的选择通常取决于现有的CAD/CAM集成以及特定树脂体系是否需要固化动力学耦合。

渗透率表征与预成型体压实

准确的渗透率输入数据是VARTM仿真中最大的误差来源。制造商通常通过以下三种方法之一确定渗透率张量值：

台架注射测试（100×100 mm至300×300 mm）：使用与生产树脂粘度匹配的测试流体进行单向或径向流动实验。成本：每种织物类型$3,000-8,000，包括测试面板制造和数据缩减。
压实耦合渗透率测量：在真空袋压力（95-100 kPa）下用测力传感器测量预成型体厚度的同时测量渗透率。这考虑了真空施加期间发生的15-30%的厚度减少，在典型NCF铺层中纤维体积分数从45%增加到55-58%。
逆参数估计：通过迭代调整渗透率，使简单注射测试的流动前沿位置与仿真匹配，直到模型和实验流动前沿在5% RMS误差阈值内收敛。

纤维体积分数（Vf）是主导变量：根据Carman-Kozeny关系，Vf从50%变为58%可使面内渗透率降低约40-60%。下表显示了在55% Vf下常见碳纤维结构的典型渗透率值：

织物结构	K11（0°方向）×10⁻¹⁰ m²	K22（90°方向）×10⁻¹⁰ m²	K33（厚度方向）×10⁻¹⁰ m²
双轴NCF（±45°）	2.8-4.2	1.8-2.9	0.08-0.15
三轴NCF（0°/+45°/-45°）	3.5-5.5	0.8-1.6	0.06-0.12
单向（0°）	6.0-8.0	0.3-0.6	0.05-0.10
平纹编织（3K, 200 g/m²）	1.5-2.5	1.5-2.5	0.10-0.20
短切毡（450 g/m²）	4.5-7.0	4.0-6.5	0.50-1.00

铺覆仿真与纤维取向

铺覆仿真预测平直织物在双曲面模具表面上的变形方式。对于大型VARTM部件，精确的铺覆至关重要，因为纤维重新取向会改变局部渗透率张量，从而改变树脂流动路径。铺覆仿真的关键输出包括：

局部纤维角度与预期取向的偏差（剪切角，NCF织物起皱前通常为0-30°）
由于织物在曲面区域的压实造成的厚度变化（内半径区域减少5-15%）
双曲率区域相邻织物层之间的间隙形成，这会产生沟流通道
面内剪切刚度和锁定角（织物开始起皱的剪切角，碳纤维NCF通常为25-45°）

PAM-RTM集成的铺覆模块使用带有指数松弛模型的运动学映射算法，在2米级验证部件上实现与实验摄影测量相比±3°以内的纤维角度预测精度。RTM-Worx依赖CATIA复合材料部件设计模块的外部铺覆，该方法基于有限元，精度略高（±1.5°），但需要单独的CATIA许可。

沟流预防策略

沟流——树脂沿预成型体边缘与模具壁之间的间隙、嵌件周围或铺层递减区域优先流动——是大型VARTM部件中干斑缺陷的最常见原因。沟流通道宽度可达0.5-3.0毫米，其渗透率比织物本体渗透率高100-1,000倍。

基于仿真的缓解策略包括：

边缘密封优化：对密封胶带放置（通常3-6毫米宽）进行建模，以阻止周边沟流。仿真根据预成型体厚度和树脂凝胶时间预测最佳密封间距（距净修边线30-80毫米）。
流道介质通道与流道布局：分布介质（高渗透率聚合物网格，K≈1×10⁻⁸ m²）在入口表面建模为离散渗透层。入口流道之间的最佳间距（通常100-400毫米）通过平衡填充时间和分布层中的树脂浪费来确定。
周边排气管路：仿真确定真空排气口位置，防止空气在最后填充区域滞留。高级仿真模拟排气阀顺序（在流动前沿通过时顺序关闭排气口）以最小化树脂排气量，目标排气比为部件重量的5-10%。
嵌件和芯材沟流：夹芯结构中的泡沫和轻木芯材嵌件沿芯材/预成型体界面产生优先流动路径。对此界面处0.2-0.5 mm间隙的仿真预测局部填充时间减少40-60%，需要局部增粘或粘合薄膜屏障。

验证案例：60米风电叶片根部

一家欧洲主要叶片制造商对60米涡轮叶片根部截面（弦长3.2米，宽2.8米，包含单向和三轴碳纤维织物，总计1,200 kg）的VARTM仿真与生产数据进行了验证。PAM-RTM仿真预测总填充时间为38分钟，与实际生产填充时间41分钟吻合（误差在8%以内）。仿真正确识别了后缘泡沫芯嵌件处的干斑风险区，经首个生产部件的超声波C扫描验证。在该位置添加螺旋缠绕真空排气口后，后续部件消除了该缺陷，废品率从12%降至1.5%。

常见问题

VARTM工艺设计的最低可行仿真精度是多少？

对于生产工装验证，填充时间预测与实验值偏差在±15%以内被认为可接受用于初始工艺设计。对于缺陷预测（干斑、孔隙），要求流动前沿形状与RMS误差低于10%的相关性。这些目标通过良好表征的渗透率数据和适当的边界条件建模是可以实现的。

在不知道确切间隙尺寸的情况下，如何模拟沟流？

标准方法是灵敏度分析：模拟三种间隙尺寸情景（0.2 mm、0.5 mm、1.0 mm），覆盖制造公差带。如果流动模式对间隙尺寸变化不敏感（填充时间变化<10%），则沟流不构成风险。如果各情景之间流动时间变化超过30%，则必须使用密封胶带或增粘剂控制间隙。0.5 mm间隙仿真通常与生产数据的最佳匹配。

VARTM仿真能否处理伴随固化的非等温树脂注射？

可以——PAM-RTM和RTM-Worx都支持耦合热-固化分析。树脂粘度使用Kamal-Sourour固化动力学模型建模为温度和固化度的函数。对于注射时间较长（>30分钟）的大型部件，这种耦合至关重要，因为放热反应可使树脂温度升高15-40°C，将粘度降低2-5倍，并可能在薄截面处导致提前凝胶化。