材料工程基础课程大作业.doc.docVIP

安 徽 建 筑 大 学 材料与化学工程学院 《工程基础》 大作业 班级姓名 2017 年 月 1 计算模拟 3 1.1 导热问题数值解法 3 1.1.1概述 3 1.1.2 数值解法原理 3 2数值模拟案例 9 2.1MALTAB PDE-tool工具箱 9 2.2计算案例 9 2.3模拟具体步骤 10 2.3.1在MATLAB命令窗口，输入命令 pedtool 10 2.3.2 建立几何模型 11 2.3.3 设置边界条件 12 2.3.4 PDE方程的选择 13 2.3.5 网格划分 14 2.3.6 绘温度场分布图 15 1 计算模拟 1.1 导热问题数值解法 1.1.1概述 随着计算机的普及应用和性能的不断改善，以及相关的数值计算方法的发展和应用程序的开发，传热学数值计算方法作为数值求解传热问题的有效工具也得到了相应的发展。利用计算机求解传热学问题愈来愈受到人们的普遍重视，而且在计算复杂传热问题中显示出它的优越性，因而成为传热学的一个重要的分支。数值传热的相关内容也很自然地成为工程类学生学习传热学课程的不可缺少的部分。 1.1.2 数值解法原理 (1)建立控制方程及定解条件 对于一些更为复杂的导热问题，如复杂的几何形状和边界条件以及物性变化较大的情况，分析求解往往很复杂或者根本不可能。此时求解问题的唯一途径是利用数值分析的办法获得数值解 数值求解通常是对微分方程直接进行数值积分或者把微分方程转化为一组代数方程组再求解。这里要介绍的是后一种方法。如何实现从微分方程到代数方程的转化又可以采用不同的数学方法，如有限差分法、有限元法和边界元法等。 简要地介绍用有限差分析方法从微分方程确立代数方程的处理过程。有限差分法的基本思想是把原来在时间和空间坐标中连续变化的物理量，如温度、压力、速度和热流等。用有限个离散点上的数值集合来近似表示。有限差分的数学基础是用差商代替微商，导数。而几何意义是用函数在某区域内的平均变化率代替函数的真实变化率。 在右图中可以看出有限差分 表示的温度场与真实温度场的区别。图中用T0、T1、T2…表示连续的温度场T；Δx为步长，它将区域的x方向划分为有限个数的区域，Δx0、Δx1、Δx2…，它们可以相等，也可以不相等。当Δx相等时，T1处的真实变化率a可以用平均变化率b、c或d来表示，其中b、c和d分别表示三种不同差分格式下的温度随时间的变化率，即: b为向后差分格式： c为向前差分格式: d为中心差分格式: 这种差分格式也可以推广到高阶微商的情形。对于二阶微商的差分格式可以在一阶差分格式的基础上得出: (2)区域离散化 这里我们以一个矩形长柱体的非稳态导热过程为例来讨论区域离散化问题。如果不考虑矩形长柱体长度方向上的温度变化，那么它是一个二维非稳态导热问题，如图， 上图表示了长柱体矩形截面上区域离散化的情况。图中可见，对于给定的 空间区域，在x方向上的步长为Δx，在y方向上的步长为Δy，用它们作为空间尺度可以将矩形区域划分成纵横交错的网格，交点称为节点。然后以节点为中心，在两个节点的中心处划分界限，定出节点的控制面积，对于三维情况则为控制体积或控制容积，因而常在一般意义上称之为节点的控制体。 控制体的形状是随着坐标系的不同而改变的，这里的控制体是一个个的矩形面积。网格的步长在每一个方向上可以均匀划分，也可以不均匀的划分。因此，选用不同的步长和不同的划分方法，可以将同一区域划分出不同大小、不同数目的控制区域，以及不同数目的节点数。显然，随着步长的不断减小，节点数目的不断增加，由节点温度表示的离散的温度场就会更加接近连续的温度场，但计算工作量也会随之增加。在时间方向上离散化的步长常用Δη来表示，Δη的选取也是可大可小的，也可以随时间的进程而变化。显然，无限小的时间步长Δη亦会使得离散温度变化接近连续的温度改变，但随之而来的是相应的计算工作量的增加。 （3）建立节点物理量的代数方程 控制体热平衡法建立节点方程的过程是将能量守恒方程应用于控制体，建立该节点与周围节点之间的能量平衡关系式，再利用傅立叶的导热定律，最后获得控制体节点温度与周围节点温度之间的关系式。 考察右图中的节点P及其控制体，由能量平衡关系应有 式中，QW、QE、QS和QN分别为邻近节点W、E、S和N通过传导方式传给节点P的热流量；QV为单位时间控制体内热源的发热量，ΔΕ为控制体单位时间内热能的增加量。由导