3铸造成形_温度场模拟2008.ppt

第三讲 铸造成形温度场数值模拟 热传递的基础 热量传递方式可以通过三种方式进行：热传导，热对流，热辐射。 这三种传热方式在铸造成形过程都存在。 (1) 热传导 物体各部分之间不发生相对位移时，依靠分子、原子及自由电子等微观粒子的热运动进行的热量传递称为热传导，简称导热。 只有在物体处于不同温度时，热量才能从一个物体传递到另一个物体，或从物体的某一部分传递到物体的另一部分。热总是从温度高的地方流向温度低的地方，铸件凝固冷却时，铸件内部的温度高于外界，因此铸件内部向其外侧以及铸型传递热量。 在三维迪卡尔坐标系统，连续介质各点在同一时刻的温度分布叫做温度场，温度场的一般可表达为：T=?(x,y,z,t)。 若温度场不随时间变化，则称做稳定温度场，由此产生的导热为稳定导热；若温度场随时间改变，则称做不稳定温度场，不稳定温度场的导热为不稳定导热。 导热的基本定律Fourier定律 在归纳大量实验基础上，J.B.Fourier于1882年指出，单位时间内由热传导而通过单位面积的热量（比热流量）与温度梯度成正比。 在一维空间，Fourier定律表示成： 三维热传导微分方程 如物体中无内热源，根据能量守恒定律： 流入热量－流出热量=微元体内蓄热量的增加 而单位时间内微元体蓄热量增量： 整理得到： 令 ， 为导温系数（m2/s），得： (2) 热对流 热对流是指流体中温度不同的各部分相互混合的宏观运动引起热量传递的现象。 对流换热可以用Newton冷却定律来描述，即 式中，q为热流密度； 为对流换热系数； 为流体的特征温度； 为固体边界温度。 (3) 热辐射 物体通过电磁波传递能量的方式称为辐射。物体会因各种原因发出辐射能，其中因热的原因发出辐射能的现象称为热辐射。 辐射换热可以用Stefen-Boltsman定律来描述，即 式中，q为热流密度， 为表面的绝对温度； 为辐射黑度； 为Stefen-Boltsman常数。 数值模拟（CAE）软件组成 三大基本模块： 前处理：离散，构造迭代计算的数组空间 计算分析：计算，迭代计算产生解数值空间 后处理：图像显示，数字图形，可视化动态化 华铸CAE前处理－网格划分 STL文件 STL(Stereolithography)是为快速原型制造RPM (Rapid Prototype Manufacture)服务的文件格式，是由美国3D SYSTEMS公司于1988年制定的一个接口协议，被工业界认为是目前的准标准。 华铸CAE前处理－网格划分 RPM技术问世后，各商用CAD软件公司纷纷推出了带有STL文件输出格式的软件包，包括基于工作站与微机的软件系统。这种文件格式类似于有限元的网格划分，它将物体表面划分成很多小三角形，用这些空间三角形小平面来逼近原CAD实体。文件的数据结构非常简单，而且独立于CAD系统。STL有文本格式及二进制两种形式。 在STL文件中，保存有三角形小平面的矢量信息，用于表达物体表面的内外指向，每个小平面必须有一个单位矢量从实体内指向实体外部。“右手法则”规定物体的内腔，它确定每个小平面的顶点的次序，并检验法矢的指向。 网格剖分原理 1 切片法 有限差分法(FDM)的网格单元多为正交六面体。网格剖分的任务就是将铸件、铸型等目标对象的三维实体在三个坐标轴方向上分别剖片，其结果将得到由众多小六面体“堆砌”而成的目标对象。 2 射线穿透法 基于STL的射线穿透法原理是：用平行于一坐标轴的一组射线来穿透STL文件记录的三维实体，由于STL文件保存的是很多三角形小平面，而射线和任意三角形小平面的关系只有三种：射线落在三角形小平面内，射线与小平面相交，射线与小平面不相交。求出射线与三角形小平面的交点坐标，再将坐标排序，根据交点所对应的三角形小平面的法向矢量来判断射线上某一段和整个实体的关系，从而判断此段上的网格的材质属性。 数值模拟的总思路 1 设计工艺方案 2 求解数理方程 3 揭示行为细节 4 推断质量缺陷 5 改进工艺方案 返回第2步循环 实用求解方法 数值求解要素 有限差分求解的性质 数学模型 温度场 Fourier 方程 Fourier 方程 潜热的处理 等效比热法 温度回升法 热焓法 * * 材料科学与工程学院凝固模拟研究室材料成形与模具技术国家重点实验室 廖敦明 陈立亮 2008级硕士研究生材料成形原理及数值模拟 实际铸件边初值条件复杂，难求分析解 借电子计算机高速大容量，可求数值解 空间离散化 方程差分化 分离出时间变量 按时间步长