声明:文中所有观点仅代表个人见解,若有不妥之处欢迎指正。
这篇论文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的生成
任意两纤维不重叠,公式表示:
RVE周期性,公式表示:
![](https://img-blog.csdnimg.cn/direct/05f5f95b10eb427a8ac690777aa4d77e.png)
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验证预测导热系数的准确
![](https://img-blog.csdnimg.cn/direct/fec682f91d12484c8ace2e4c03dbb337.png)
![](https://img-blog.csdnimg.cn/direct/94b74e2413374da791ea595d96536909.png)
![](https://img-blog.csdnimg.cn/direct/54946dcb91314952bf21442df3d31b86.png)
5.2 纤维横向各向同性导热系数,体积分数,长径比对金属基复合材料导热系数的影响
5.2.1纤维横向各向同性和各向同性对结果的影响:结论为不显著
横向各向同性的纤维和各向同性的纤维作比较:如下
此外,本文还定义了个参数表示两者的差距
T为横向各向同性,I为各向同性
这个参数为
这个参数的计算结果:
5.2.2体积分数,长径比对金属基复合材料导热系数的影响
结论:
体积分数增加,导热系数线性降低
长径比影响不大
体积分数图:
长径比图: