基于非傅里叶热传导层状复合陶瓷结构有限元分析.pdfVIP

基于非傅里叶热传导的层状复合陶瓷结构有限元分析 基于非傅里叶热传导的层状复合陶瓷结构有限元分析· 王小贞1 张士元2郑百林3贺鹏飞3 1福州大学物理与信息学院福建福州，350002 2福建工程学院土木工程系福建福州，350007 3同济大学航空航天与力学学院，上海，200092 摘要层状复合陶瓷结构的热涂层和夹层厚度达到了微、纳米尺度，其温度分析模型应采 用微尺度下的非傅里叶热传导模型，因此其热应力与宏观尺度下的傅里叶热传导定 律计算的结果有一定差异。本文将非傅里叶热传导模型与傅罩叶热传导模型相结合 所得的层状复合结构的温度场作为热载荷，运用有限元方法对热应力进行分析。研 究表明，由于非傅里叶热传导模型考虑了热的波动性，其热涂层的热应力比按傅里 叶热传导定律计算的结果高，夹层可以有效地减缓热涂层的拉应力，同时夹层热物 理性能参数(如：松弛时间、声速)对热涂层应力场有明显影响。 关键词层状复合结构，微、纳米尺度，非傅里叶热传导，热应力，有限元 为了满足微电子、航空、航天以及高温应用等技术的要求，抵抗热冲击的陶瓷材料往 往制作成层状复合结构形式。在宏观尺度下，忽略了某些微尺度(包括时间和空19)因素的 影响下，基于傅里叶热传导定律建立的宏观热传导，成功地解决了许多自然界和工程领域 所遇到的热传导问题，人们对此已做了一定研究Il邡筇l。但随着人们对微尺度领域的深入 探索和研究，特别是纳米技术(空间微尺度)和高频短脉冲加热技术(时间微尺度)的出现，经 典傅里叶热传导定律所隐含的热量传播速度为无穷大的假设已不再成立，傅里叶热传导定 律已不再适用，应该采用统计力学的观念和量子力学来解决微、纳尺度问题。 中热传播的波动现象，并定量测定了热波传播速度比声速小一个数量级。近几年来随着 微尺度研究的深入和制作工艺的提高，许多学者建立了各种时间和空间微尺度条件下的 导热模型：Qiu和Tienl8，9岐寸超短脉冲激光(不产生电流)加热金属薄膜的研究中提出了微 绝缘体和半导体空间微尺度的迟滞行为时也得出同样的热传导方程；M旬umdar(1993) 114,15,161在研究半导体声学薄膜时提出了纯声子辐射模型。 本文针对纳米尺度层状陶瓷复合结构(热涂层和夹层为纳米尺度的陶瓷层)，将非傅 里叶热传导模型(用于热涂层和夹层)和宏观经典傅里叶热传导模型(用于结构层)相结合 求解温度场，用有限元方法分析其热应力，并与完全采用傅里叶热传导模型(用于热涂层、 ‘福建省科技厅(JB06035)、国家自然基金资助项目． 电话 ·200· 数学、力学、物理学、高新技术研究进展——2008(12)卷 夹层和结构层)所得的热应力进行了比较，并研究了夹层不同物理性能(如：松弛时间、 声速)对热应力的影响。 1．温度数学模型和求解 对于纳米尺度层状复合陶瓷结构(见图1)，由于热涂层和夹层是绝缘并且其尺寸接近 绝缘声子的平均自由程(10～500nm)，可采用纯声子散射模型，结构层厚度达到宏观尺寸， 采用经典傅里叶热传导定理计算温度场。 对于陶瓷热涂层和夹层，双曲型纯声子散射方程(1)如下。 图1超薄陶瓷热涂层夹层结构 V2h9r5N 其中c是声子的平均速度，fⅣ是声子系统中动量守恒的一般过程中的松弛时间，h是 声子系统中发生动量损失的碰撞过程中的松弛时间。 0、 对于两个结构层，采用经典傅里叶热传导方程(2)如下： a2下 a下 兰—}=口号：；=-…………(6Xb+h，2b+hx2b+2办)(2) 撒‘ Ot 其中口为热扩散系数。 热涂层和夹层厚度为b，结构层厚度都为h；初始条件为，当时间为零时，在热涂 层、夹层和结构层内的温度为瓦(3000c)并均匀分布，并且与外界无热对流和交换，即保 持绝热；热冲击作用下，边界条件为戈=0处温度突然从瓦降低到 处始终保持绝热；热涂层和结构层之间良好热交换接触，即保持温度和热流的相等