温度场在铸造工艺CAD中应用

《温度场在铸造工艺CAD中应用》由会员分享，可在线阅读，更多相关《温度场在铸造工艺CAD中应用（49页珍藏版）》请在金锄头文库上搜索。

1、材料科学与工程学院 温度场在铸造工艺CAD中应用 报告人：陆皓 1 材料科学与工程学院 温度场在铸造工艺CAD中应用 一、差分法解热传导方程 二、结晶潜热的处理 三、缩孔预测算法 四、缩松预测算法 五、课题相关 2 材料科学与工程学院 差分格式的建立 3 材料科学与工程学院 向前差商 向后差商 中心差商 差商法：用差商代替微商，用微商代替微分 一阶差商 差分格式的建立 4 材料科学与工程学院 二阶差商 差分格式的建立 5 材料科学与工程学院 差分法解热传导方程 差分格式的建立 6 材料科学与工程学院 Q1Q2 X-XXX+X 单元自身内能变化Q： 流入热量Q1： 流入热量Q2： Q=Q1+Q2

2、 A 差分格式的建立 7 材料科学与工程学院 i,ji-1,ji+1,j i,j+1 i,j-1 单元自身内能变化Q： 流入热量Q1： 流入热量Q2： Q1Q2 Q3 Q4 流入热量Q3： 流入热量Q4： 差分格式的建立 8 材料科学与工程学院 差分格式的建立 同理Q=Q1+Q2+Q3+Q4 当 9 材料科学与工程学院 差分格式的稳定性与收敛性 求解之前所得的差分格式方程式时必须要考虑差分格式 的稳定性和收敛性： 当空间步长x0和时间步长t0时，差分方程的解 是否逼近偏微分方程的解，即差分方程收敛性问题。 差分方程在计算过程中所产生的误差是不断扩大，还是 可以控制，即差分方程稳定性问题。 在适

3、当的条件下，稳定性是收敛性的充分必要条件，所 以在这里我们只考虑稳定性条件。稳定性就是：如果初始条 件和边界条件有微小的变化时，若解的最后变化是微小的， 则称解是稳定的，否则是不稳定的。 10 材料科学与工程学院 差分格式的稳定性与收敛性 一维稳定性条件 二维稳定性条件 三维稳定性条件 对于一维显示格式的差分方程式的稳定性条件为 项的系 数大于等于零，即： 同理，对于二维和三维中（ ） 对于 的情况可选择其中最小的一个代入稳定 性条件进行计算。 11 材料科学与工程学院 铸件凝固过程热传导差分格式的求解条件 在数值模拟过程中，除了上述差分格式及收敛条件外，还需 要考虑一些其他的求解条件。 1.

4、几何条件：铸件和铸型的形状及网格参数； 2.物理条件：铸件、铸型和大气的物理参数 导热系数，热容，密度等 3.初始条件 初始条件就是指在t=0时，被求对象各处的值。此时，应为各单 元的初始温度。一般设定，铸件温度为浇注温度，铸型和冷铁 为室温或者预热后温度。忽略浇注过程的热量传递，初始时刻 为液体完全充满金属铸型的状态。 12 材料科学与工程学院 4.边界条件 在铸件凝固过程中，根据材料种类分类的边界有：冒口和大气 、铸件与铸型、冷铁与铸型、铸型和地面、铸型与大气等。根 据传热的特性分类为：绝热、给定热流密度、定温、热传导、 热对流和热辐射六种类型的边界条件。 4.1对称面（绝热） 如果铸件和

5、铸型是对称的，则可以将铸件和铸型从对称面处分 开，只计算其中一半或者更少。此时对称面设置为绝热边界， 可在绝热面上虚设一层单元，参与计算，但是由于绝热，对系 统温度场无影响。 铸件凝固过程热传导差分格式的求解条件 13 材料科学与工程学院 4.2铸件-铸型界面（热传导） 铸件和铸型由于导热系数不同，界面换热系数可分为理想和非 理想状态（利用串联热阻叠加特性）： 非理想状态：存在界面热阻h 铸件凝固过程热传导差分格式的求解条件 14 材料科学与工程学院 铸件凝固过程热传导差分格式的求解条件 4.3铸型-空气界面（对流和辐射的综合换热） h=hc+hr 经验公式 P21 4.4铸型-地面（定温）

6、地面的温度可以认为是恒温，则此边界是定温边界。在计算它们 之间的热传导系数时可把它处理成为理想接触的情况。 15 材料科学与工程学院 结晶潜热的处理 金属凝固过程中，结晶潜热的释放是金属凝固过程区别于一般 导热过程的显著特点。 其中结晶潜热我们可以通过下式得到： 其中，fs为固相率。这里我们考虑固相率和当前温度与液固相 线是线性关系。 16 材料科学与工程学院 结晶潜热的处理 对于不同的合金采用不同的方法，一般分为两种： 对于有一定结晶温度范围的合金采用等价比热法； 对于窄结晶温度范围的合金或者纯金属采用温度回升法。 等价比热法：对于具有一定凝固温度范围的合金，我们采用“ 液相线以下的领域中，

7、由于固相率的增加所释放出来的潜热相 当于散失的热量”这一观点来处理潜热的释放，即物体的比热 可以理解为单位质量的物体降低单位温度时释放的热量，同理 ，单位质量具有凝固温度范围的合金在凝固阶段降低单位温度 时释放的热量也可以认为等于其比热。 17 材料科学与工程学院 结晶潜热的处理 此时释放的比热由两部分组成，为物体的真比热和凝固潜热， 这种比热我们称为等价比热。等价比热计算公式为： 等价比热， 真比热， 单位质量金属降低单位温度时释放的潜热 18 材料科学与工程学院 结晶潜热的处理 这种方法假定在液固相线间的凝固阶段，在任何温度下降低单 位温度释放的潜热量是均等的。但实际上金属在凝固阶段，降

8、低单位温度凝固时所析出的固相量，随着降温而减少。 所以合金的比热为： 式中 液态比热 固态比热， 19 材料科学与工程学院 结晶潜热的处理 CpL C Cp Q0 TLTs Cps Ts TL 时间轴上的潜热处理情况 Q0单位质量金属降低单位温度时释放的潜热，kj/kg。 20 材料科学与工程学院 结晶潜热的处理 温度回升法 21 材料科学与工程学院 结晶潜热的处理 温度回升法 22 V1 V2 V3 材料科学与工程学院 缩孔预测算法 缩孔产生的机理： 阶段:从浇注温度到液相线温度，金属液无相变，温度下降引起体积收缩 。 阶段:从液相线到固相线温度区间，属于凝固阶段，相变引起体积收缩。 阶段

9、:从固相线到室温，金属液已完全凝固，温度下降引起体积收缩。 23 材料科学与工程学院 缩孔预测算法 缩孔预测过程： 计算温度场 计算固相率场 孤立熔池划分 确定熔池顶部单元，计算顶部单元体积 计算熔池收缩量 比较熔池收缩量和顶部单元体积 修改单元属性 24 材料科学与工程学院 缩孔预测算法 缩孔预测的步骤： 1.压力场的计算 压力场的计算主要是为了判断熔池（所有相互连通的可流动单元组成的区 域。）顶部单元的位置和数量，如果当前单元为可流动单元（固相率小于 临界固相率的单元，即可以流动补缩的单元），则其静压力从其底部取值 。 压力计算分为两部分： 1.1 单元自身重量产生的静压力 其中 自身重量

10、产生的静压力； 当前计算单元Y方向上的步长 1.2 当前单元的顶部可流动单元对该单元的压力P 其中 顶部可流动单元对该单元的压力； 顶部可流动单元与 当前单元的高度差。 25 材料科学与工程学院 缩孔预测算法 缩孔预测的步骤： 2.熔池顶部单元的确定 由上述对两种单元所受压力的物理意义我们可以判断出，顶部单元为其顶 部没有其余液态单元的部分，也就是Ps =P的单元即为熔池顶部单元。其除 了受到自身的静压力之外，不再受到任何其余压力。 3.收缩量的计算 总的收缩量是由单元的收缩量之和组成，在时刻t与t+ t之间，熔池的总 收缩量为所有熔池顶部单元的收缩量。收缩包括三种收缩，分别为，液态 收缩，液

11、固收缩和固态收缩。 26 材料科学与工程学院 缩孔预测算法 缩孔预测的步骤： 4.缩孔的形成计算 判断出凝固进程中t到t+ t之间熔池的顶部单元之后，将其体积与熔池总 得收缩量的累积进行比较，当总得收缩量大于熔池顶部单元时，便将熔池 顶部单元设置为缩孔 缩孔的形成是个连续的过程，所以我们在模拟的过程中，不断的累加熔池 收缩量，并与熔池顶部单元体积进行比较。直到符合上式时，便将顶部单 元设置为缩孔，并在收缩量中减去缩孔的体积，继续向下一步进行累积。 27 材料科学与工程学院 缩孔预测算法 熔池的划分 28 材料科学与工程学院 缩孔预测算法 熔池的划分 29 材料科学与工程学院 缩孔预测算法 熔池

12、的划分 30 材料科学与工程学院 缩松预测算法 缩松产生的机理： 结晶温度宽的合金其凝固方式通常体积凝固为主，液态金属中的 小晶核容易发展成发达的树枝晶。 当树枝晶形成骨架时，尚未凝固的金属液被分割成一个个不连通的 小熔池。 小熔池的金属液发生液态收缩和液固收缩，已经凝固的金属则发生 固态收缩，两者之间的体积差会在枝晶间产生间隙。 当枝晶间隙没有足够的金属液进行补缩时，会产生细小的孔洞，这 就是缩松。 31 材料科学与工程学院 缩松预测算法 缩松的形成判据：改进新山判据 32 材料科学与工程学院 缩松预测算法 式中 t f 凝固时间；V 被计算铸件某部分的体积；A 被计 算铸件某部分的散热面积

13、 缩松的形成判据：改进新山判据 推导过程 另一方面，由凝固收缩而引起的流动速度u为： 式中 凝固收缩率 33 材料科学与工程学院 缩松预测算法 缩松的形成判据：改进新山判据 推导过程 由（2）（3）式可得 因为冷却速度 又因为金属液在枝晶间的流动属于多孔介质间的流动, 符合达尔西公式, 即 , 所以有 流路压力差, 即压力损失 34 材料科学与工程学院 缩松预测算法 缩松的形成判据：改进新山判据 推导过程 由公式(7)可知:G/ R 越小, 就越大,金属液流动所需的压力差就越 大, 说明金属液流动的阻力就越大, 对凝固收缩的补缩越困难, 则容易 形成缩松。如果取某一临界值Kc (为常数), 使

14、G/ RK c 为产生缩松 的判据。这就是G/ R 方法的由来。显然, G/ R 只考虑了压力损失的影 响,而忽略了压力P的影响。如果铸件的某个部位在凝固过程中, 对此部 位进行补缩的金属液在其流路上的压力损失较大, 但此部位有足够的压 力使金属液能够克服这一压力损失, 到达需要补缩的凝固部位,则此部 位不会产生缩松, 否则将产生缩松。 在凝固过程中的铸件内部, 某一部位的凝固收缩会导致未凝固部位的金 属液向此部位流动进行补缩,这一流动是在一定的压力作用下进行的。 如果在这一流动的流程上阻力相当大, 即流程压力损失相当大,金属液 还没有流到被补缩部位其压力就已经完全损失, 凝固部位得不到补缩而

15、 产生缩孔。如果金属液的压力足够大, 能够克服流程压力损失, 凝固部 位可以得到金属液的补缩, 则不会产生缩松。 35 材料科学与工程学院 缩松预测算法 缩松的形成判据：改进新山判据 推导过程 P 当前计算单元的压力；P0与当前计算单元相邻单元的压力 显然，当P0时，金属液流动的压力全部损失殆尽，该单元将得不到相邻 单元的补缩而产生缩松，所以满足下式将会产生缩松。 由 可得 36 材料科学与工程学院 缩松预测算法 式中：G 温度梯度，/cm；R冷却速度，/min；K为缩松产生的 临界值 为两单元的距离，上下左右为一个空 间步长，对角线为 在二维中，需要分别计算次点与周围八个 单元中可补缩单元的

16、温度梯度，然后在这 八个值中取最大值来判断。即，当只要存 在可以进行补缩的单元时，我们便认定他 不会产生缩松。 37 材料科学与工程学院 缩松预测算法 式中：G 温度梯度，/cm；R冷却速度，/min ；K为缩松产生 的临界值 式中： 单元当前时刻和一个时间步长后的温度 时间步长 K为一个实验测定的数值，在铸钢件中，缩松产生的临界值是 38 材料科学与工程学院 目前工作概况动态网格技术 动态网格的提出 在模拟计算过程中，一个非常重要的问题就是模拟过程的计算速度和精度 。对于现有的设备来说，计算一个几千万网格的铸件，大概需要24个小时 （工作站上，8核32个线程）。 对于一个比较大的件进行计算的时候，将 铸件分成几千万个网格时，其网格尺寸也可能达到几个毫米甚至几十个毫 米，在一些边缘地带，网格将无法准确的描述出铸件的轮廓，在