基于球型填充相复合材料有效导热系数的计算.pdf
《基于球型填充相复合材料有效导热系数的计算.pdf》由会员分享,可在线阅读,更多相关《基于球型填充相复合材料有效导热系数的计算.pdf(6页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第3 1 卷第9 期2 0 0 8 年9 月合肥工业大学学报(自然科学版)J O U R N A LO FH E F E IU N I V E R S I T YO FT E C H N O L O G YV o L3 1N o 9S e p t 2 0 0 8基于球型填充相复合材料有效导热系数的计算刘祥宽1,M e h a r iS a l o m o n1 2,胡献国1(1 合肥工业大学机械与汽车工程学院,安徽合肥2 3 0 0 0 9;2 德国达姆斯达特应用大学,达姆斯达特D-6 4 2 9 5)摘要:文章运用数值单元体法建立了计算填充二硫化钼聚甲醛基复合材料有效导热系数的模型,提出了用
2、有限元法和用求单元体模型的等效热阻来确定该复合材料等效导热系数2 种计算方法,通过计算表明2 种方法的计算结果比较接近,并运用上述方法研究了二硫化钼含量对聚甲醛基复合材料导热系数的影响规律。关键词:复合材料;等效导热系数;有限元中图分类号:T Q 3 1 4 2 5 8文献标识码:A文章编号:1 0 0 3 5 0 6 0(2 0 0 8)0 9 1 3 7 8 0 4E f f e c t i v et h e r m a lc o n d u c t i v i t yc a l c u l a t i o nm o d e l so ft h ec o m p o s i t e sw
3、i t hs p h e r i c a lf i l l i n g sL I UX i a n g-k u a n1,M e h a r iS a l o m o n1”,H UX i a n-g u o1(1 S c h o o lo fM a c h i n e r ya n dA u t o m o b i l eE n g i n e e r i n g,H e f e iU n i v e r s i t yo fT e c h n o l o g y,H e f e i2 3 0 0 0 9,C K l l a;2 H o c h s e h u l eD a m m t a
4、d t,U n i v e r s i t yo fA p p l i e dS c i e n c e s,H a a r d t r i n g1 0 0,D-5 4 2 9 5D a r m s t a d t,G e m m n y)A b s t r a c t:T h en u m e r i c a le l e m e n tm e t h o di Su s e dt oe s t a b l i s han u m e r i c a lm o d e lt Op r e d i c tt h ee f f e c t i v et h e r m a lc o n d u
5、 c t i v i t y(E T C)o ft h ep o l y o x y m e t h y l e n e(P O M)c o m p o s i t e sf i l l e dw i t hm o l y b d e n u md i s u l f i d e(M o S 2)p a r t i c l e s Af i n i t ee l e m e n tm e t h o da n da ne f f e c t i v eh e a tr e s i s t a n c em e t h o da r ea d o p t e dt Oc a l c u l a
6、t et h ee f f e c t i v et h e r m a lc o n d u c t i v i t yo ft h ec o m p l e xm a t e r i a l T h ec a l c u l a t i o nr e s u l t so ft h et w om e t h o d sa r ec l o s et oe a c ho t h e r,a n dt h ee f f e c to ft h ef i l l e rc o n t e n to nt h ee f f e c t i v et h e r m a lc o n-d u c
7、t i v i t yo ft h ec o m p l e xm a t e r i a l i Sa l s os t u d i e db ym e a n so ft h et w om e t h o d sa b o v e-m e n t i o n e d K e yw o r d s:c o m p l e xm a t e r i a l;e f f e c t i v et h e r m a lc o n d u c t i v i t y;f i n i t ee l e m e n t0 引言复合材料构件具有较好的自润滑性和耐磨性,且综合性能优异,已经广泛应用于航空
8、、汽车、精密仪器及电子电器等行业。但复合材料表面摩擦温升对其摩擦磨损性能具有重要的影响,大部分聚合物复合材料的导热性能差、熔点不高,若摩擦产生的热量不能及时地排出,将会影响构件的正常工作,导热系数是表征复合材料热性能的主要参数。?国内外对复合材料传热特性进行了一系列的研究,目前已就两相的复合材料的导热与物理性能建立了多种预测模型j 如B r u g g e n m e n 模型及C h e n g-V o c h o n 模型等 1 。对于填充聚合物复合材料,M a x w e l l 等通过求解电场能量的方程推导出基体球形粒子复合材料的电导率计算公式。E u e k e n 将电导率换为导热
9、系数,从而得到M a x w e l l E u c k e n 方程 2 3。计算复合材料的导热系数时,不仅需要考虑填料的含量,还需考虑各相的尺寸、形状、排列方式等微观结构形式 3“。此外,一些研究者认为复合材料的导热系数计算可归结为单元体等效导热系数的计算,而单元体的等效导热系数可以通过数值计算或有限元计算来获得。采用数值计算法分析以塑料相为连续相的复合材料时,先计算各通道的并联热阻,再将各通道并联热阻串联得到等效热阻 5 3;有限元法是运用收稿日期:2 0 0 7-0 9 1 2;修改日期:2 0 0 7 1 1 1 2基金项目:国家自然科学基金资助项目(5 0 4 7 5 0 7 1)
10、;安徽省自然科学基金资助项目(0 7 0 4 1 4 1 5 2)作者简介:刘祥宽(1 9 8 2-),男江西占安人,合肥工业大学硕上生;胡献N 1(1 9 6 3-),男,安徽安庆人,博士,合肥工业大学教授,博士生导师 第9 期刘祥宽,等:基于球型填充相复合材料有效导热系数的计算1 3 7 9单元体法建立计算其有效导热系数的模型 6,7 3。本文通过对复合材料的微观结构分析,提出了一种导热系数的计算模型。以这种模型为研究对象,设计了2 种不同的计算方法,并研究了分散相颗粒含量对塑料基复合材料导热系数的影响,2种计算方法的计算结果比较接近。1复合材料导热系数的计算方法聚甲醛基复合材料是一种二相
11、结构材料,聚甲醛(P O M)为基础相,二硫化钼(M o S 2)为分散相,鉴于这两相之间的混容性很小,导热系数并不等于各组分热物理性质的简单平均值,它与复合材料中各相组分的导热系数、容积比、分散相在复合材料中的形状、排列以及取向等因素有关 8 。1-1P O M 基复合材料基本单元的计算模型根据最小热阻力法则和比等效导热系数相等法则,只要复合材料的单元体与总体有相等的比等效热阻,不论单元体的尺度大小,这种单元体与总体的导热系数相等 9 ,这样,研究复合材料的等效导热系数,只要研究与之相等的单元体的导热系数即可。在传热模型中将整块复合材料看成是由大量芯部包含一个球形颗粒的正方体组成的密堆体,热
12、流从顶面进入,如图1 所示。设正方体的边长为z,颗粒半径为r,体积分数为口,可求得颗粒体积分数7 与参量r 的关系为r=(等)专圈1 单元体分析模型1 2 复合材料等效导热系数的有限元计算方法根据以上对复合材料的基本单元计算模型的假设,可以采用数值单元法求复合材料的导热系数。选取其中一个立方体作为研究对象,分散相M o S z 位于该立方体中心,采用有限元分析软件模拟单元体传热,在顶面施加第一类边界条件,即施加温度载荷;底面施加第三类边界条件,即施加对流载荷,左右两边界设为绝缘边界条件。设单元体的有效导热系数为A。,由傅里叶定律知A。=啬(2)其中,T 为单元体上下表面的温差;q 7 为通过单
13、元体的热流密度;L 为单元体厚度。温差与通过单元体的热流密度可由有限元计算结果得到。1 3 确定复合材料等效热阻方法热流在材料中传递过程就象电流在导体中传导一样,在传递过程中材料会对热流有阻力,好比电流在导体中受到电阻,这种阻力称为热阻。借用电路中的一些计算原则和方法,对于均质材料的热阻,由下式得R L(A A)(3)其中,A 为材料的导热系数;L 为热流通道长度;A为热流通道面积。复合材料的等效热阻计算可以归结为单元体热阻计算,假设单元体尺寸为单位长1,把单元体划分成若干个微元块,形成热阻的网格图,然后根据(3)式计算出微元块的热阻,而微元块之间的连接等效为它们热阻的连接,这样整个单元体的热
14、阻就是由各微元体热阻之间的串、并联得到。为了简化模型,把沿着热流方向划分的微元体归结为一个通道,把在热流垂直方向划分的微元体归结为一层,如图2 所示。(b)圈2 单元体模型的热阻网格划分 1 3 8 0合肥工业大学学报(自然科学版)第3 1 卷当分散相为球体,球体的半径为r,在计算其热阻时,把热流通道分作两大系统:系统表示为在热流主方向上只有连续相的通道;系统表示为在热流主方向上遇到球分散相的通道。在系统中的热阻为R 7 一L(钡1)一(1 2 r)(A 1)(4)在系统中的热阻为以球心为原点(见图2 a),在Y 处取一微元体厚曲,在该层的热阻可以看成是截面积为7 c(产一Y 2)的分散相热阻
15、与截面积为1 7 c(r 2-y 2)的连续相热阻的并联。微元体层的热阻为】d R2 丽F 了丽a l y i 万 了甄(5)其中,天。为连续相的导热系数;A z 为分散相的导热系数。在整个含有分散相层中的热阻,由一r到r 的积分得到R J 二丽顽瓦而 芳可i 驴i 弼(6)单元体总的热阻就是由系统与系统组成,由(3)、(4)、(6)式可以得到单元体的有效导热系数的计算公式为址串=州 L 而赢丽正翟芊F 而+(1 2 r)(7)当A:J:l。时,数学处理得A并一串=A,“-n三;手享;弓三手萝萎享垩琴需ccA:A,一,一,+1-2r,(8)2 导热系数的计算结果与分析2 1有限元法模拟的单元体
16、温度分布在计算单元体等效导热模型中,假设单元体为单位长1。已知聚甲醛基体的导热系数为0 2 3 0 4W(1 T I K),二硫化钼的导热系数为2 4 8 5W(m K)。单元体顶面施加2 5 的温度载荷,底面施加第三类边界条件,施加对流载荷,环境温度为2 0,对流换热系数为1 2 5W(m 2 K),单元体周围绝热,使用A N S Y S 模拟的温度云图,如图3 所示。从图3 可以看出,由于聚甲醛与M o S z的导热系数不同,而且为不相容相,在温度云图上所表现出温度梯度不是规则的递减,而是与M o S 2 分散相在聚甲醛基体的含量与形状有关。图4 表示了沿单元体中心路径的温度分布情况。由图
17、4 可以看出,单元温度不是随着路径线性的下降,在通过分散相路径时温度变化比较小,几乎成一条水平线。另外分散相质量分数对温度分布影响大,分散相M o S 2 质量分数越高水平线越长,质量分数为0 8 的单元体中心路径温度变化比较快,水平线最短,上下两面的温度下降最快。质量分数为1 0 的单元体中心路径温度变化最缓慢,水平线也最长,上下两面的温度下降最小,微球体的热阻就比较小,这表明在聚甲醛中添加M o S 2 可以改善材料的导热性能 1 0 。2 4 1 8 5 一-、2 j 6 4彳藩Il12 2、I 一一刍。-r,一_-二A 2 1 4 6 7-蝴4 丽F一2 0图3 单元体温度云图00 2
18、0 40 60 81 O中心距离岬圈4 单元体中沿路径温度变化曲线2 2 分散相质量分数的影响以分散相M o S 2 颗粒质量分数分别为0 8、1 o、3 O、5 o、8 0、1 0 0 为研究对象,由(1)式可以得到在不同质量分数下的分散相M o S 2 微球半径。用同样的方法做出复合材料的计算模型,分别采用有限元模拟方法和确定复合材料等效热阻方法来求其等效导热系数,得到分散相M o S 2 质 第9 期刘祥宽,等:基于球型填充相复合材料有效导热系数的计算1 3 8 1量分数对复合材料的导热系数的影响。从表1 中可以看到随着分散相M o S 2 质量分数的增加,复合材料的导热系数增加,并且添
19、加分散相M o S 2 对整个复合材料有效导热系数的影响很大,添加3 o、5 o 的M o s 2 后,复合材料的导热系数提高了2 5 4、3 7 9,说明在聚甲醛基中添加M o S 2 分散相能有效地改善其导热性能,这对熔点不是很高的复合材料来说是非常重要的。表1 不同情况下复合材料导热系数w(m K)由图5 可以发现,在分散相质量分数为1 0 以内,2 种方法计算结果吻合比较好,特别是在M o S z 为3 和5 时,发现2 种计算结果仅相差0 0 0 12 和0 0 0 52。另外M o S z 的质量分数过高,就不能均匀的分散在聚甲醛基体里,会发生团聚,使复合材料的整体性能下降。,釜白
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 填充 复合材料 有效 导热 系数 计算
限制150内