Effect of fiber transverse isotropy on effective thermal conductivityof metal matrix composites论文笔记

声明:文中所有观点仅代表个人见解，若有不妥之处欢迎指正。

这篇论文Effect of fiber transverse isotropy on effective thermal conductivity of metal matrix composites reinforced by randomly distributed fibers于2016年被 "Composite Structures"录用。原文链接：Redirecting

本文的主要工作：

1、基于代表性体积元（RVE）的有限元均质法评估导热系数

2、生成RVE：RSA方法的改进

3、比较实验和Hashin–Shtrikman bounds验证预测导热系数的准确

4、结论：横向各向同性的纤维和各向同性的纤维，体积分数，长径比对金属基复合材料导热系数的影响

目录：

1、RVE的生成

2、控制方程

3、温度的周期性边界条件

4、纤维横向各向同性导热系数：局部坐标和全局坐标

5、结论

1、RVE的生成 

任意两纤维不重叠，公式表示：

P1和P2是任意两个纤维；x1和x2是分别位于纤维P1和P2表面上的任意两点；d(x1,x2)是两点的最短距离。纤维不重叠，就是最短距离大于0。

RVE周期性，公式表示：

左侧是平移后的纤维部分；ki取值为-1，0，1，该值取决于纤维x穿透RVE的几个面。二维的示意图如下所示：

详细的算法请参考

Representative volume element for composites reinforced by spatially randomly distributed discontinuous fibers and its applications

2、控制方程

稳态传热公式：

是温度周期性边界条件；是温度梯度；k是导热系数张量。这是abaqus中用到的，一般在编程中不用到。

热通量向量：

q这里是向量，<q>=(<q>1,<q>2,<q>3)T。

q的各分量：

这里，neint和nint分别是单元和RVE的积分点数；W(yI)是单元e中位于积分点yI处的权重；

3、温度的周期性边界条件

分为两部分和，

右侧是常值，故可以写成下面这种形式

rve上的点可以分为三类：顶角处、边上、面上；如图所示

温差通过一个参考点引入周期性边界条件：

4、纤维横向各向同性导热系数：局部坐标和全局坐标

坐标转换公式：G全局坐标，l局部坐标

Rjl是

其中各参数如下：

前两个参数分别是与z轴的夹角，在xoy面上的投影与x轴的夹角。w如图所示，它的作用在于：当纤维的横截面不是对称的时候，例如椭圆，那么纤维旋转就会有一个相应的w。

导热系数的转换：

5、结论

5.1 比较实验和Hashin–Shtrikman bounds验证预测导热系数的准确

Hashin–Shtrikman bounds计算公式：

Hashin–Shtrikman bounds计算

Hashin–Shtrikman bounds计算结果与预测结果：

与实验结果的比较：

5.2 纤维横向各向同性导热系数，体积分数，长径比对金属基复合材料导热系数的影响

5.2.1纤维横向各向同性和各向同性对结果的影响：结论为不显著

横向各向同性的纤维和各向同性的纤维作比较：如下

此外，本文还定义了个参数表示两者的差距

T为横向各向同性，I为各向同性

这个参数为

这个参数的计算结果：

5.2.2体积分数，长径比对金属基复合材料导热系数的影响

结论：

体积分数增加，导热系数线性降低

长径比影响不大

体积分数图：

长径比图：