1. 研究背景：热塑性复合材料的“新宠”地位

玻璃纤维增强聚丙烯（GF/PP）等连续纤维增强热塑性复合材料，正成为汽车工业的明星材料。它们具备：

• 轻量化

• 快速成型

• 可回收利用

• 优异的抗冲击性能

但问题也随之而来：热塑性基体（如PP）表现出明显的压力依赖性和速率依赖性塑性行为，这与传统热固性复合材料（如碳纤维/环氧）截然不同。直接套用为热固材料开发的损伤模型，往往会导致预测失败。

因此，建立一个适用于热塑性复合材料的微观力学模型，成为工程与学术界的迫切需求。

2. 研究目的：让计算机成为“虚拟试验场”

本研究的目标是：

建立一个高保真的微观力学模型，能够预测GF/PP复合材料在横向拉伸下的损伤行为，并识别影响其性能的关键设计参数。

具体包括：

• 系统确定统计意义上的代表体积元（RVE）

• 建立压力相关塑性 + 损伤 + 界面脱粘的本构模型

• 通过实验验证模型有效性

• 作为材料设计的虚拟测试工具

3. 研究方法：三步走，步步为营

3.1. RVE的确定——不只是“画个框”

方法一：两点概率函数（几何路径）

• 计算纤维相的空间相关性

• 确定相关长度，保证RVE大于该长度

• 发现GF/PP具有层状纤维簇分布，RVE必须覆盖整个单层厚度（≈255μm）

方法二：Hill-Mandel均匀化（力学路径）

• 施加均匀应变边界条件

• 计算均匀化弹性张量的特征值

• 特征值稳定时的窗口尺寸即为RVE最小尺寸

✅ 最终RVE尺寸：300×255μm2

3.2. 有限元模型构建（Abaqus）

• 建立3D RVE模型

• 施加周期性边界条件

• 对边界上的纤维进行移除与再生成，保证纤维体积分数不变

• 最小纤维间距限制为0.4μm，确保网格质量

3.3. 本构模型三大模块

组分	模型	关键特性
玻璃纤维	线弹性	各向同性
PP基体	Drucker-Prager塑性 + 韧性损伤	压力依赖性、速率依赖性（简化）、损伤演化
纤维/基体界面	双线性内聚力模型	最大应力判据 + 断裂能控制

4. 研究过程与重要图分析

图1 & 图2（原文Page 3-4）：微观结构 + 两点概率函数

图1 展示了GF/PP的典型层状纤维分布：中部密集，两侧富树脂。

图2 显示不同方向的S2(r)：

0∘（沿纤维层）波动小

90∘（穿过纤维层）波动大，波长≈15倍纤维直径

✅ 结论：RVE必须包含整个单层厚度，且宽度≥300μm。

图3（Page 4）：RVE尺寸判定

图3a：S2(r) 随窗口尺寸增大趋于稳定
图3b：标准差同步减小

✅ 几何RVE最小尺寸：300×255μm2

图4 & 图5（Page 5）：均匀化弹性张量特征值

图4：不同窗口尺寸的3D FEM模型
图5a：6个特征值在窗口≥125×255μm后稳定
图5b：标准差收敛

✅ 力学RVE最小尺寸：125×255μm2，最终取300μm宽为安全尺寸。

图6（Page 7）：最终RVE模型

• 包含纤维、基体、界面（内聚力单元）

• 满足周期性边界条件

• 网格以六面体为主

图7（Page 8）：基体应变率分布 + 速率依赖性简化处理

图7a：基体局部应变率比宏观高48.9%
图7b：插值得到该应变率下的PP应力-应变曲线

✅ 创新：无需复杂粘塑性模型，简化且稳定。

图8（Page 8）：基体损伤模型

图8a：弹塑性损伤模型示意
图8b：损伤变量随等效塑性位移演化

• 损伤起始：等效塑性应变

• 损伤演化：断裂能 + 特征长度 → 避免网格依赖

图9（Page 8）：界面内聚力模型

• 双线性牵引-分离法则

• 损伤起始：最大应力判据

• 损伤演化：线性软化，由断裂能控制

图10（Page 9）：模型验证——关键图

图10 展示模型预测与实验结果对比：

• 弹性模量：预测 4.41GPa，实验 4.43GPa

• 非线性段与失效点均在实验散布范围内

✅ 模型可作为虚拟测试工具

图11（Page 10）：损伤演化过程

• 损伤先在纤维密排区萌生

• 最终形成横向贯穿裂纹

• 裂纹穿过富树脂区，提供额外阻力

✅ 结论：纤维间距是控制损伤的关键。

5. 研究重难点

难点	解决方案
热塑性基体压力依赖性塑性	双曲线Drucker-Prager模型
速率依赖性	简化方法：平均应变率 + 插值曲线
纤维成层分布	两点概率函数 + 层厚方向分析
界面脱粘	内聚力单元 + 最大应力判据
网格依赖性	特征长度 + 断裂能
周期性边界下的纤维截断	移除 + 再生成边界纤维

6. 研究结论：谁才是真正的“控制者”？

模型能力

• 成功预测GF/PP在横向拉伸下的应力-应变曲线、损伤起始、裂纹扩展

• 与实验高度吻合

关键设计参数影响排序

参数	影响程度
纤维间距	★★★★★（最关键）
界面强度	★★★★☆
界面断裂能	★★★★☆
基体强度/塑性	★★☆☆☆
基体韧性	★☆☆☆☆（影响不大）

重要设计建议

• ❌ 提高基体韧性不是最优策略

• ✅ 优化纤维分散性（减少聚集） 更有效

• ✅ 提高界面断裂能 有助于延迟失效

---

7. 未来展望

尽管模型已具备较强的预测能力，仍可向以下方向扩展：

1. 多轴载荷：压缩、剪切、组合载荷

2. 多尺度：将RVE响应嵌入层合板模型

3. 动态加载：更精确的速率依赖性模型

4. 工艺影响：热残余应力、纤维取向

5. 材料数据库：构建不同参数下的性能图谱

---

附：原文图片位置索引

内容	原文页码
微观结构与两点概率函数	Page 3-4, Figures 1–2
RVE尺寸判定	Page 4, Figure 3
均匀化弹性张量	Page 5, Figures 4–5
最终RVE模型	Page 7, Figure 6
应变率分布与速率处理	Page 8, Figure 7
基体损伤模型	Page 8, Figure 8
界面内聚力模型	Page 8, Figure 9
模型验证（核心）	Page 9, Figure 10
损伤演化过程	Page 10, Figure 11

注：更多关于复合材料多尺度有限元rve建模的前沿知识小编之前有推荐，可以详查置顶文章：数智融合：FE²、PINN、能量法与AI代理模型在结构工程中的应用

如果您觉得文章不错，欢迎点赞、关注、收藏及转发~