热传导和热辐射大作业报告.doc

WORD资料整理 完美格式可编辑 热传导与热辐射大作业报告 热传导与热辐射大作业报告 WORD资料整理 完美格式可编辑 目录 TOC \o 1-1 \h \u HYPERLINK \l _Toc5749 一、作业题目 PAGEREF _Toc5749 - 1 - HYPERLINK \l _Toc18253 二、作业解答 PAGEREF _Toc18253 - 2 - HYPERLINK \l _Toc13344 个人感想 PAGEREF _Toc13344 - 17 - HYPERLINK \l _Toc11721 附件.计算中所用程序 PAGEREF _Toc11721 - 18 - 完美格式可编辑 一、作业题目 一矩形平板, ，内有均匀恒定热源，在及处绝热，在及处保持温度，初始时刻温度为，如右图1所示： 1、求时，矩形区域内的温度分布的解析表达式； 2、若，，，，，热传导系数 ，热扩散系数。请根据1中所求温度分布用MATLAB软件绘出下列结果，加以详细物理比较和分析： 300s内，在同一图中画出点、、、、（单位：m）温度随时间的变化； 200s内，画出点、、、、（单位：m）处，分别沿x、y方向热流密度值随时间的变化； 画出时刻区域内的等温线； 300s内，在同一图中画出点（单位：m）在分别等于，，情况下的温度变化； 300s内，比较点(9,6) （单位：m）在其它参数不变情况下热导率分别为、和的温度、热流密度变化； 300s内，比较点(9,6) （单位：m）在其它参数不变情况下热扩散系数分别为、和的温度、热流密度变化； 3、运用有限差分法计算2中(b)、(d)和(e)，并与解析解结果进行比较，且需将数值解与解析解的相对误差减小到1‰以下； 4、附上源程序和个人体会； 以报告形式整理上述结果，用A4纸打印上交。 二、作业解答 1、求时，矩形区域内的温度分布的解析表达式； 解答：我们令，则可以得到一个方程和边界条件： （1-1） 将上式分解为一个的稳态问题： （1-2） 和一个的其次问题： （1-3） 其中 则原问题的解根据下式求得： （1-4） 发热强度为常数的特解可从表2-4中查的，则新变量可定义为： （1-5） 将（1-5）带入（1-2）整理得到： （1-6） 若令常数，则上式可以变为： （1-7） 其中 假定可以分离出如下形式： （1-8） 对应于的分离方程为： （1-9） （1-10） 在中特征值问题的解可以直接从表2-2第6条中得到，只需要用a代替L， （1-11） （1-12） 是下面方程的正根： （1-13