复合材料的热传导
复合材料热传导的理论基础
复合材料的热传导特性
CFRP和GFRP等复合材料的热传导与均质材料有什么区别？
纤维增强复合材料具有强烈的异向性。纤维方向的热导率主要由纤维贡献，垂直方向由基体（树脂）控制。结果是不同方向的k值可能相差10倍以上。
有效热导率的理论
纤维方向（并联模型）和垂直方向（串联模型）的有效热导率如下。
这里 $V_f$ 是纤维体积分率，$k_f$ 是纤维的热导率，$k_m$ 是基体的热导率。
具体数值的差异有多大？
对于PAN系碳纤维（$k_f=10$ W/(mK)）、环氧树脂（$k_m=0.2$ W/(mK)）、$V_f=0.6$ 的情况：
- $k_{\parallel} = 0.6 \times 10 + 0.4 \times 0.2 = 6.08$ W/(mK)
- $k_{\perp} \approx 0.45$ W/(mK)
相差超过13倍。忽视这一点用各向同性来分析会导致温度分布完全不同。
Hashin-Shtrikman极限
更精确的评估使用Hashin-Shtrikman的上下界。
实测值位于这个上下界之间。使用Halpin-Tsai模型或有限元均质化（RVE分析）可以进行更精确的预测。
纤维配向随机的短纤维材料会怎样？
随机配向的情况趋向各向同性，但注射成形品在流动方向会出现纤维配向，导致部分异向性。从Moldflow等注射成形仿真中导出纤维配向张量，并将其映射到热导率张量的方法已经得到实际应用。
复合材料法则，19世纪50年代至今
Maxwell（1873年）首次理论化了分散球状粒子的复合材料的等效热导率。他的公式至今仍用于铜粒子混合聚合物基底的TIM（热界面材料）设计，形式为λeff ≈ λm(λp+2λm+2φ(λp−λm))/(λp+2λm−φ(λp−λm))，保留在教科书中。
复合材料热传导的数值计算方法
RVE分析均质化
需要建模到纤维级别并进行分析吗？
使用代表体积单元（RVE）进行多尺度分析是标准方法。建立纤维直径7μm的碳纤维六角排列的RVE，在各个方向施加温度差，求取有效热导率张量。
RVE的大小如何确定？
纤维直径的10～20倍是经验值。如果增大RVE尺寸结果不变，则足够。COMSOL和Digimat支持自动RVE生成和参数化均质化。
积层板的建模
对于积层板（例：[0/90/45/-45]s），需要将每层的热导率张量按积层方向进行旋转并叠加。在Abaqus中，通过ORIENTATION和SHELL SECTION定义各层的配向角。
| 积层构成 | 面内k [W/(mK)] | 面外k [W/(mK)] |
|---|---|---|
| UD [0]8 | 6.0 / 0.45 | 0.45 |
| 交叉 [0/90]2s | 3.2 / 3.2 | 0.45 |
| 准各向同性 [0/45/90/-45]s | 3.2 / 3.2 | 0.45 |
交叉层和准各向同性的面内方向变得均匀了。
是的。但无论哪种结构，面外方向都由基体控制，热导率仍然较低。确保面外方向的热路径是CFRP结构热设计的最大课题。通过Z销或碳纳米管的面外引入来改善的研究正在进行。
并联法则和串联法则的精度差异
纤维增强树脂（CFRP）的面内热导率可以用并联法则（混合法则）预测，误差在±5%以内，但厚度方向即使用串联法则误差也常超过20%。1990年代NASA Langley在碳纤维/环氧树脂中进行实测，建议厚度方向使用Hashin-Shtrikman上下限模型（NASA TM-4756）。
复合材料热传导的实务应用
实务中的注意事项
复合材料热分析中最需要注意什么？
最重要的是材料坐标系的正确定义。由于每层的纤维方向不同，需要准确设置每个单元的坐标系。
Ansys ACP的连接
使用Ansys Composite PrepPost（ACP）可以自动化复合材料积层定义到热分析模型的转换。
1. 在ACP中定义积层构成（层顺序、配向角、厚度）
2. 自动生成材料坐标系
3. 传输数据到稳态热分析
4. 自动应用每层的各向异性k
不用手动设置坐标系是大优点。
Abaqus的Abaqus/CAE复合材料铺层功能可实现等效功能。COMSOL的复合材料模块提供对应支持。
与实测值的比较
复合材料的热导率因测量方法而异。
| 测量方法 | 适用范围 | 精度 |
|---|---|---|
| 激光闪光法 (LFA) | 面外方向 | ±5% |
| 稳态法（保护热板法） | 面外方向 | ±3% |
| Angstrom法 | 面内方向 | ±10% |
面内和面外的测量方法不同。
激光闪光法测量面外方向的热扩散率，通过比热和密度换算为 $k = \alpha \rho c_p$。面内方向的测量样品加工困难，精度也较差。比较分析和实测值时需考虑这些不确定性。
智能手机散热材料的演进
2010年代初期的智能手机只是贴铜箔石墨片（λ≈400 W/m·K），但2019年之后的高端机（三星Galaxy S10及以后）采用气相沉积石墨（VGCFs）复合材料，面内热导率超过1500 W/m·K。
复合材料热传导的软件比较
专用工具的比较
有专门针对复合材料热传导的工具吗？
下面比较同时支持均质化（微观）和宏观分析的工具。
| 工具 | 均质化 | 宏观分析 | 特点 |
|---|---|---|---|
| Digimat (e-Xstream/MSC) | RVE自动生成 | FEA连接 | 注射成形材料纤维配向连接 |
| COMSOL Composite | 内置RVE | FEM | 多物理场耦合 |
| Ansys ACP + Mechanical | Ansys Materials均质化 | FEM | Workbench工作流 |
| Abaqus + Micromechanics插件 | RVE分析 | FEM | Python自动化强大 |
Digimat在注射成形连接上很强。
Digimat可以从Moldflow或Moldex3D中读入纤维配向张量，自动映射各单元的各向异性k。对于短纤维强化PA66和PPS等树脂成形品的热分析，实际上是标准工具。
数据库的应用
复合材料热物性值波动大，使用可信的数据库很重要。
- AGATE/NCAMP: 航空航天用预浸料认证数据
- Granta (Ansys): 材料数据库，包含复合材料特性
- CAMPUS: 热塑性树脂物性数据库
不用自己测量，可以从数据库查询。
航空航天领域的认证预浸料数据充实，但新材料必须进行实际测量。数据库的值仅作参考，在设计初期阶段使用。
Ansys Composite PrepPost的问世
Ansys（原ANSYS）于2010年发布ACP（Ansys Composite PrepPost），实现了按积层角度自动转换热导率张量的功能。此前用户通常手写Python脚本进行坐标变换矩阵，输入错误导致的分析事故时有发生。
复合材料热传导的先端研究
CNT/石墨烯强化复合材料
用碳纳米管或石墨烯提高热导率的研究进展如何？
单根CNT的热导率超过3000 W/(mK)，但复合材料中的界面热阻（Kapitza阻力）限制了改善。目前面外k约从0.5改善到1.5 W/(mK)。
$G_{interface}$ 是CNT-基体界面的热导，约为10～100 MW/(m2K)。
多功能复合材料
多功能复合材料同时具有结构和热功能，其设计使用多尺度热分析。例如EV电池壳体的CFRP化，需同时优化结构强度和散热性。
同时满足结构和热的设计很难。
使用多目标优化（帕累托前沿）可视化结构重量和最大温度的权衡，决定积层构成和厚度。使用Ansys optiSLang或modeFRONTIER。
3D打印复合材料
连续纤维3D打印（Markforged等）中，打印路径决定了纤维方向。从打印路径数据直接映射热导率张量的工作流已经开发。
制造过程和分析直接连接。
设计制造一体化（Design-for-AM）与热分析的融合是未来大趋势。可以利用打印方向的异向性积极优化散热路径。
CNT复合材热导率的极限
碳纳米管（CNT）单体热导率超过3000 W/m·K，但CNT/环氧树脂复合材料中界面热阻（Kapitza阻力）占支配地位，实测值仅5～10 W/m·K。2004年Philip Kim等人的纳米级热输运界面阻力定量化论文（Nano Lett.）成为转折点。
复合材料热传导的故障排查
常见问题
复合材料热分析中最容易失败的是什么？
1. 材料坐标系的不匹配
问题: 假设纤维方向沿全局X方向输入k值，但实际零件形状中纤维沿着涂装线弯曲。
对策: 对于曲面零件，使用Ansys ACP或Abaqus复合材料铺层功能自动沿表面生成坐标系。平面投影定义坐标系会在角部产生偏差。
2. 各向同性假设的错误
"各向异性很麻烦，就用各向同性的平均值吧"可以吗？
完全不可以。面内k=3 W/(mK)、面外k=0.45 W/(mK)的平均值1.7来输入，会导致面外方向的温度降低被低估3.8倍。对于面外热阻占主导的电子基板等问题是致命误差。
3. 铺层厚度变化的热影响
厚度过渡部（铺层厚度变化处）的横截面变化会产生局部温度梯度集中。网格不能追随铺层厚度变化形状会导致温度场不准确。
对策: 手动细分铺层厚度变化区域的网格。使用ACP的涂装功能准确反映铺层厚度变化位置。
4. 纤维体积分率的波动
实际零件中 $V_f$ 与设计值相差约±5%。$V_f$ 变化5%会导致 $k_{\parallel}$ 变化约5%，$k_{\perp}$ 变化约15%。面外方向敏感性高，需考虑波动建立设计安全裕度。
面外方向总是瓶颈。
是的。复合材料热设计中"如何降低面外方向热阻力"是永恒课题。
忽视各向异性的设计失误
2008年某卫星制造商用各向同性（λ=5 W/m·K）对CFRP构体进行热分析，忽视了厚度方向实际的λ≈0.6 W/m·K，导致面板温度低估30℃。发射后温度超限，需在轨道上重新编程加热器功率，这一真实案例已被报告。
帮助
更多
报告