计算机在材料科学中应用.ppt

计算机在材料科学中的应用 讲授内容 基本原理 基本方法 软件的使用及应用 只能在基本层次上进行 要 求 听课 记录 上网查询资料 写出课程小结 几大内容 数学模型 数值分析方法 物理场的数值模拟 加工过程控制 材料检测 数据与图像处理 Internet的应用 材料分类 按照组成和结构分 金属、无机非金属、有机高分子、复合材料 按照性能特征和作用分 结构、功能材料 按照用途分 建筑、能源、电子、耐火、医用、耐蚀 研究内容 材料的组成 材料的结构 材料的性能 材料的制备工艺 材料的使用性能 材料的各方面之间的相互关系 计算机研究材料 新材料的设计 计算机模拟 工艺过程的优化及自动控制 数据和图像处理 信息检索 第三章 材料科学研究中的主要物理场的数值模拟 第一节 温度场的计算 各种材料的加工、成型过程中与加热、冷却等传热过程有着密切的联系，所以利用计算机解决传热问题是极为有力的。 一.导热方程 固体在稳态条件下的Fourier 导热方程 qx 是x 方向的热流密度(W/m2);λ是材料的热导率(W/m.K) 偏导数为x 方向上的温度梯度，负号表示传热的方向与温度梯度的方向相反。 三维状态时，瞬态温度场的场变量T在有方程： ρ为材料密度(kg/m3), c为材料的比热容(J/kg.K),t为时间，三个λ分别为沿x,y,z方向的导热率（w/m.K），Q是物体内部的热源密度(W/kg). 这是一个热量平衡方程，即体元升温所需的热量应等于流入体元的热量与体元内产生的热量的总和。 若边界条件和内部热源密度Q不随时间变化则经过一定时间后物体内部各点的温度将达到平衡，即有稳态热传导方程： 此方程必须附加初始条件和边界条件(即定解条件)才能得到唯一解。 二.初始条件与边界条件 1.初始条件 T|t=0 =T0 或 T|t=0 = T0(x,y,z) 2.三类边界条件 物体表面或边界与周围环境的热交换通常有三种类情况，既有三类边界条件。 （1）第一类边界条件 已知物体边界上的温度分布函数 T|s =Tw 或 T|s = Tw(x,y,z,t) （2）第二类边界条件 已知物体边界上的热流密度 这实际上就是Fourier定律 （3）第三类边界条件 已知物体与其周围环境介质间的对流传热系数k和介质温度Tf 当k取不同的值时，上述三类边界条件均可以统一用第三类边界条件式表达，以方便计算机编程计算。 三、平面温度场的有限差分求解 1、二维稳态导热问题的求解 四条边界有四种不同的边界条件：三类边界条件和绝热条件。 (1)划分网格 △x＝xi+1- xi ， △y＝yi+1- yi (2)建立差分方程 二维各向同性、无内热源的稳态热传导微分方程有 应用四个边界条件： 对流传热边界条件 热流边界条件 绝热边界条件 给定温度边界条件 用Ti,j 表示结点(i,j)处的温度，以差商代替微商，则有： 代入并舍去误差，令△x=△y后则有 四个边界条件的差分形式如下： 对流传热边界条件 热流边界条件 绝热边界条件 给定温度边界条件 这样一个差分方程加四个边界条件式共同构成方程组，从而通过求解得到结果。 3、简单算例 P68图3－6，已知所有边界点的温度且无内热源，利用稳态导热求解线性方程组。 四、有限元法求解 1、基本方程 试探函数：构造的函数，满足边界条件 加权余量法：用构造函数计算会有余量，则对余量的加权积分使其为零。 型函数Ni KT＝P （热传导矩阵*结点温度阵＝温度载荷阵） 2、平面温度场有限元法求解 (1)单元划分：三结点三角形单元 三角形的每个顶点称为结点，按逆时针方向对3个结点进行编号i,j,m (2)试探函数：在单元中寻找满足单元的近似函数很方便，也就将试探函数称为插值函数，可以假设每个单元上的温度分布T是坐标的线性函数，有 T＝a1+a2x+a3y 它满足在三个结点上为相应结点的温度值，代入得到 T＝NiTi+ NjTj+ NmTm 其中：Nk＝(ak+bkx+cky)/(2A) (k=i,j,m) A是三角形的面积，Nk 称为型函数 (3)单元积分计算 (4)单元总体合成 (5)求解 KT＝P 3、有限元分析实例 工业烟囱壁中的温度分布情况 4、采用ANSYS程序计算潜水艇壳体的温度分布情况 第四章 材料科学与行为工艺的计算机模拟 材料行为工艺是通过调整材料在加工过程中的组织性能来改善其使用性能，利用计算机模拟材料可以部分代替传统的真实试验，提高了效率、节省费用。 第三节 相图计算发展历程和计算软件介绍 相图是描述相平衡系统的