凝固过程温度场

合集下载

铸造凝固过程数值模拟

铸造凝固过程数值模拟

铸造凝固过程数值模拟-简介

1.铸造凝固过程数值模拟

1.1 概述

在铸造生产中,铸件凝固过程是最重要的过程之一,大部分铸造缺陷产生于这一过程。凝固过程的数值模拟对优化铸造工艺,预测和控制铸件质量和各种铸造缺陷以及提高生产效率都非常重要。

凝固过程数值模拟可以实现下述目的:

1)预知凝固时间以便预测生产率。

2)预知开箱时间。

3)预测缩孔和缩松。

4)预知铸型的表面温度以及内部的温度分布,以便预测金属型表面熔接情况,方便金属型设计。

5)控制凝固条件。

6)为预测铸应力,微观及宏观偏析,铸件性能等提供必要的依据和分析计算的基础数据。4

铸件凝固过程数值模拟开始于60年代,丹麦FORSUND把有限差分法第一次用于铸件凝固过程的传热计算。之后美国HENZEL和KEUERIAN应用瞬态传热通用程序对汽轮机内缸体铸件进行数值计算,得出了温度场,计算结果与实测结果相当接近。这些尝试的成功,使研究者认识到用计算数值模拟技术研究铸件的凝固过程具有巨大的潜力和广阔的前景。于是世界上许多国家都相继开展了铸件凝固过程数据模拟以及与之相关的研究工作。

1.2 数学模型的建立和程序设计

液态金属浇入铸型,它在型腔内的冷却凝固过程是一个通过铸型向环境散热的过程。在这个过程中,铸件和铸型内部温度分布要随时间变化。从传热方式看,这一散热过程是按导热,对流及辐射三种方式综合进行的。显然,对流和辐射的热流主要发生在边界上。当液态金属充满型腔后,如果不考虑铸件凝固过程中液态金属中发生的对流现象,铸件凝固过程基本上看成是一个不稳定导热过程。因此铸件凝固过程的数学模型正是根据不稳定导热偏微分方程建立的。但还必须考虑铸件凝固过程中的潜热释放。

热型连铸凝固过程微观组织形成的数值模拟

热型连铸凝固过程微观组织形成的数值模拟

热型连铸凝固过程微观组织形成的数值模拟

热型连铸凝固过程微观组织形成的数值模拟

在金属材料的生产过程中,热型连铸是一种重要的凝固方法。通过热型连铸,可以制备出具有优良性能的金属材料。而在热型连铸过程中,凝固微观组织形成的过程对最终材料的性能起着至关重要的作用。为了更好地理解凝固过程中的微观组织形成机制,科学家们采用数值模拟方法进行研究。

热型连铸可以简单地理解为利用铸铁管将金属液注入到铸模中,通过正向或背向冷却,实现金属材料的凝固过程。在凝固过程中,液态金属逐渐转变为固态金属,同时伴随着组织结构的形成。这种组织结构及其形成机制是决定材料性能的重要因素之一。

通过数值模拟方法,研究人员可以在计算机上根据已知的物理建模和模型、热力学和凝固动力学方程,模拟热型连铸凝固过程中微观组织的演变过程。通过模拟计算,可以预测凝固过程中的温度场、相变行为以及组织形态的演变规律,为优化工艺参数和改善材料性能提供理论依据。

在热型连铸凝固过程的数值模拟中,研究人员首先需要构建合适的凝固模型。这个模型通常包含了温度场、相变过程、流动场等因素。热型连铸的凝固模型常常采用二相流动模型,将液相和固相作为两个不可混合的相进行计算。接着,研究人员需要设置合适的边界条件和初始条件,以保证模拟的准确性和可靠性。

在模拟计算中,研究人员常常使用有限元或有限差分等数值方法,将凝固模型中的方程进行离散化,进而求解数值逼近解。通过数值模拟计算,可以得到凝固过程中温度梯度、相变

速率以及组织演变规律等信息。这些信息对于工艺优化和材料性能的改善起到了重要的指导作用。

【材料成型原理--铸造】第4章 液态金属凝固过程中的传热与传质

【材料成型原理--铸造】第4章 液态金属凝固过程中的传热与传质
• 凝固速度:单位时间凝固层增长的厚度。 • 铸件凝固时间的确定方法:试验法、数值模拟法、计算法。 • 1、理论计算法
t
1V1
2b2 A1
L
C1 Ti
T浇 T20
TS
计算温度场有些假设,算出的凝固时间是近似的。 应用较少。
13/33
• 2、经验计算法——平方根定律
t
2
K2
K为凝固系数,ξ为凝固层厚度。
凝固时间与凝固层厚度的平方成正比。
计算结果与实际接近。
适合大平板和结晶间隔小的铸件。
14/33
• 3、“折算厚度”法则
R2 t
K2
R V1 为铸件折算厚度或铸件模数。
A1
由于考虑了铸件的形状因素,更接近实际,是对平 方根定律的修正和发展。
15/33
第二节 凝固过程中的传质
16/33
一、平衡凝固溶质再分配
• 该两式为平衡凝固时溶质再分配的数学模型。
19/33
CS
1
C0k0 f S (1 k0 )
CL
k0
来自百度文库
C0 f L (1 k0 )
• 3、验证 (1)开始凝固时 • 初始条件:fS0,fL1 • 则:CS=k0C0;CL=C0 (2)凝固结束时 • 初始条件:fS1,fL0 • 则:CS=C0;CL=C0/k0
• 1、假设条件: • (1)长度为L的一维体自左至右定向单相凝固; • (2)冷速缓慢; • (3)溶质在固相和液相中充分均匀扩散; • (4)液相温度梯度保持固液界面为平面生长。

钢锭凝固过程温度场数值模拟

钢锭凝固过程温度场数值模拟

Nu me r i c a l Si mu l a t i o n o f Te mp e r a t u r e F i e l d o f I n g o t Du r i n g t h e
SO l i d i f i C a t i O n Pr o c e s swenku.baidu.com
钢 可 采 用 向 浇 注 后 冒 口加 入 2 0 0 m m厚 发 热 剂 增 强 钢 锭凝 固 末 期 钢 液 补 缩 能 力 ,脱 模 时 间 为 浇 注 后 1 2 . 5 h 。
关键 词 :钢锭 ;热物 性参数 ;发热 剂 ;温度场 ;数值模拟 中图分 类号 :T G 2 4 4 + . 1 文献标 识码 :A 文章 编号 :1 0 0 1 — 4 9 7 7( 2 0 1 3 )0 5 — 0 4 1 0 — 0 6
d u r i n g t h e s o l i d i i f c a t i o n p r o c e s s wi t h t h e c o mm e r c i a l s o t f wa r e MS C. Ma r c . a n d t h e t h e r mo p h y s i c a l p r o p e r t i e s o f t h e s t e e l i n mo d e l we r e d e t e r mi n e d b y t h e mi c r o s e g r e g a t i o n mo d e 1 . wh i c h c a n p r e d i c t t h e p h a s e e v o l u t i o n o f

第二章凝固温度场

第二章凝固温度场

(2)有限元法


根据变分原理来求解热传导问题微分方程的一 种数值方法。将连续求解域分割为有限个单元组成 离散化模型,再用变分原理将各单元内的热传导方 程转化为等价的线性方程组,最后求解全域内的总 体合成矩阵。 优点:适合于具有复杂形状的铸件
(一 ) 数学解析法

应用数学方法研究铸件和铸型的传热。 铸件在铸型中的凝固极为复杂: 1.不稳定的传热 2.铸件的传热大多为三维传热 3.释放结晶潜热 4.铸件、铸型的热物理参数随温度而变

T T T T 2 ( 2 2 2 ) a T t c x y z
2 2 2

三、凝固温度场的求解方法
数学解析法


数值方法
差分法
有限元法

数学解析法



主要目的:利用传热学的理论建立表明铸件凝固过 程传热特征的各物理量之间的方程式、即铸件和温 度场数学模型并加以求解。 优点:物理概念、逻辑推理清楚,解的函数表达式 能够清楚表达温度场的各种影响因素,有利于分析 各参数变化对温度高低的影响。 缺点:只适用于简单热传导问题
例:假设具有一个平面的半 无限大铸件在半无限大铸 型冷却。条件如下: 铸件、 铸型 1.铸型和铸件的材质是均质 2.铸型初始温度为t2; 3.设液态金属充满铸型后立 即停止流动且各处温度均 匀及铸件的初始温度为t1 4.坐标原点设在铸型与铸件 接触面上。

金属凝固理论 第4章 液态金属凝固过程中的传热、传质及液体流动

金属凝固理论 第4章 液态金属凝固过程中的传热、传质及液体流动
2020/10/15
T n
Tw Tf
4/56
凝固过程中,热量传递有三种形式:传导、辐射、 对流。
以热传导为主。 热传导过程取决于温度的分布——温度场:温度
随空间和时间的变化。 T = f(x,y,z,t) Fourier热传导方程:
三维非稳态导热:
T t
2T c ( x2
2T y2
小单元,对这些小单元用差分方程式近似代替微分方程式, 给出初始条件和边界条件,然后逐个计算各单元温度。 即使铸件形状很复杂,也只是计算式和程序烦杂而已,原 则上都是可以计算的。 实测法直观、可靠性好,但不方便;解析法适宜简单件, 有许多假设,误差大。 数值模拟法比其它方法准确性高,当单元选得足够小时, 差分方程的离散误差趋于零。 数值模拟法有多种方法,有限差分法应用较多。
凝固过程中铸件与铸型 的温度分布
2020/10/15
7/56
2、铸件凝固温度场
2、凝固传热研究方法
◎解析法:假设一维导热
T 2T t a x2
通解
T A Berf ( x ) 2 at
对铸件:
边界条件
对铸型:
边界条件
T1
Ti
(Ti
T10 )erf ( 2
x) a1t
T2
Ti (T20 Ti )erf ( 2
T f (t) 第一类边界条件: 给定物体表面温度随时

第二章 凝固的温度场

第二章 凝固的温度场
清楚地表达温度场的各种影响因素,有利于直观分析各参
数变化对温度高低的影响。
缺点:通常需要采用多种简化假设,而这些假设往往并不
适合实际情况,这就使解的精确程度受到不同程度的影响。 目前,只有简单的一维温度场(“半无限大”平板、圆柱
体、球体)才可能获得解析解。
第二章 凝固温度场
10
(二) 数值方法
数值方法又叫数值分析法,是用计算机程序来求解
二维传热:
2T T 2T a x 2 y 2 t
一维传热:
T 2T a t x 2
第二章 凝固温度场
6
对具体热场用上述微分方程进行求解时,需要根据具体问
题给出导热体的初始条件与边界条件。
初始条件: 初始条件是指物体开始导热时(即 t
T2
b1T10 b2T20 b1T10 b1T20 x erf 2 a t b1 b2 b1 b2 2
T2 x x0

T20 Ti a2 t
b (T T20 ) T dt时间由铸型导出的热量 为:dQ2 2 2 Adt 2 i Adt t t x 0
T20
由在界面处热流的连续性条件可得:
Ti b1T10 b2 T20 b1 b2
T1
b1 1c1 1
b2 2 c2 2

铸件凝固过程温度场的数值模拟

铸件凝固过程温度场的数值模拟

1 引言
2 前处理 单 元
铸造是一个液态金属充填铸型型腔 ,并在其 中凝 固 和冷却的高温过程 ,包括 了许多对铸件质量产生重要影 响的物理过程和现象。由于缺乏直接考察这一过程的有 效手段 , 导致长期以来 , 对铸造过程的认识及工艺方案 的 设计主要依靠工程技术人员的经验积累和现场试制。这 不可避免地导致铸件质量不稳定 、 工艺准备周期长 、 试制 成本高等许 多不利于铸造行业发展的因素
3 中央处 理 单元
从总体上讲前处理单元所涉及为纯计算机图形学 问
题, 一般采用“ 拿来 主义” 即使用成品软件 , , 而无需另行
开 发。
中央处理单元是模拟系统 的核心。 在有限元 ( 或有限 差分) 网格划分基础上 , 给出凝 固过程 中每个小单元 ( 按
点处理 ) 各个时刻温度 , 以文件形式存贮 , 并对缩孔 、 缩
Ab ta t B s d o ic s in a o t c a im f o se d e t x h n e p o e s d f r n i l q a in a d i nt 1 s c : a e ndsu so b n r me h n s o n t a y h a c a g r c s , i e e t u t n t f i e. n e ae o s i e e n o u i n o e t o d c t e t a d l mp r t r ed i a t g s l i c t n h v e n d r e t l . mc ts l t f a n n t o h c ma h mai l mo e t c f o e e au e f l n c si o i f a i a e b e e i d wi r e i n di o v he v n r b e n s l n n me c l i lt n a e n tai n o r c s h n o e . a t o lms d r u t i u r a mu a i sd mo s t f e a t r k h l s p a e s i s o r o f o si

多晶硅铸造过程温度场模拟仿真

多晶硅铸造过程温度场模拟仿真

多晶硅铸造过程温度场模拟仿真

引言

太阳能电池作为一种清洁能源越来越受到广泛的关注。其光电转换效率很大程度上取决于多晶硅的质量,而多晶硅质量又取决于硅锭定向凝固过程中温度等工艺条件的控制。因此,对多晶硅凝固过程中温度场进行模拟是确定和优化工艺条件的高效、重要技术手段。

目前,国内外已经有一些学者在多晶硅凝固温度场数值模拟方面进行了研究,比如美国的马里兰大学对多晶硅定向凝固炉和热交换炉的温度场进行了模拟比较分析。美国纽约州立大学的郑丽丽博士对太阳能多晶硅定向凝固炉进行了计算模拟。中国有色金属研究总院的刘秋娣等也对多晶硅锭凝固过程的影响因素进行了分析及数值模拟。以往的研究通常假设了特定的边界条件,并且往往缺少实际温度的测量数据。因此,多晶硅铸锭炉温度场模拟过程中边界条件的确定仍然是一个关键问题。

本文提出一种基于PID控制原理对多晶硅铸锭炉边界条件进行反算的方法,并根据反算得出的边界条件对多晶硅定向凝固炉的温度场进行研究。

1多晶硅定向凝固工艺

a)绝热罩温度恒定

b)各固体元件交界界面上无接触热阻

c)忽略炉内气体对流

2.3控制方程

根据多晶硅铸锭炉的传热方式,本文采用FLUENT中的P-1和Rosseland辐射传热模型模拟铸锭炉内的传热。相邻物体之间的导热采用Fourier导热定律,非稳态导热的控制方程:

fluent凝固过程模拟

fluent凝固过程模拟

fluent凝固过程模拟

Fluent凝固过程模拟是通过ANSYS Fluent软件来模拟凝固过程的数

值模拟方法。数值模拟方法通过建立数学模型和计算方法,对实际的自然

界中的流体流动和传热传质等物理现象进行定量描述和预测,从而实现巨

大规模和复杂性的工程系统的分析和设计。

在凝固过程中,液体从高温态到低温态的过程中会发生相变,将液体

转变为固体。对于金属、陶瓷等材料的凝固过程,Fluent可以提供精确

的描述和预测。而对于复杂的非晶态材料,Fluent的凝固模拟同样可以

提供有用的信息,以改进材料制备过程和优化材料性能。

Fluent凝固过程模拟的基本原理是基于多相流动和传热传质理论。

通过数学模型和相应的计算方法,可以描述液态金属在凝固过程中的流动

规律和温度场分布。具体而言,可以采用用于描述流动的Navier-Stokes

方程和用于描述传热传质的能量方程和质量守恒方程。

在Fluent的凝固模拟中,需要定义几个关键的物理和数值参数。首先,需要定义材料的物性参数,如热导率、密度和比热等。这些参数将影

响凝固过程中的能量传递过程。其次,需要定义边界条件,如入口温度、

入口速度和壁面热通量等。这些边界条件将影响材料的凝固速度和凝固结构。

Fluent的凝固模拟可以提供流动场、温度场、凝固率、凝固结构等

信息。通过对这些信息的分析和比较,可以优化凝固过程的设计和操作策略。例如,可以控制壁面热通量的分布,以改变凝固结构的均匀性。此外,还可以通过调整流动条件,如入口速度和剪切力,来控制凝固速度和凝固

结构。

Fluent的凝固模拟应用广泛,包括金属制造、陶瓷制造、半导体制

2.凝固温度场的测定

2.凝固温度场的测定

凝固温度场的测定

一、意义和目的

铸件温度场是指浇注后,某一时刻铸件内部的温度分布规律。在温度场中,向着铸件中心的方向上单位长度的温度变化率称为温度梯度。本实验主要是测定凝固时期各个时刻的铸件温度场。铸件凝固时期的温度场越陡,即温度梯度越大,则铸件冷却得越快,它的凝固速度就越大。

某一时刻铸件温度场中温度相同点所组成的面称为等温面。对于在一个温度范围(结晶间隔)内凝固的合金而言,铸件断面中由达到液相线温度的点所组成的面,称为液相线等温面,或称为液相边界。同样,由达到固相线温度的点组成的面,称为固相线等温面,或称为固相边界。

凝固过程中,铸件断面上液相边界和固相边界之间的区域谓之凝固区域,也就是铸件凝固过程中凝固并存区域。

阐明凝固时期各个时刻的凝固区域大小和它从铸件表面向铸件中心移动规律的曲线,称为凝固动态曲线。测定凝固动态曲线能够比较全面地描绘铸件凝固过程和研究这种过程。

凝固区域的大小即宽度决定了铸件的凝固方式,即逐层凝固方式、糊状凝固方式和中间凝固方式。

凝固区域于狭窄,铸件越是倾向于逐层凝固方式。这种凝固方式的铸件容易形成集中缩孔,便于采取措施(例如用冒口)去除铸件中的集中缩孔;铸件的热裂倾向性小和金属液充型能力较好。

凝固区域越宽,铸件越是倾向于糊状凝固方式。这种凝固方式的铸件容易形成分散性的缩孔即缩松,即使采用冒口亦难以消除这种缩松;铸件的热裂倾向性大和金属液充型能力差。

金属和铸型两方面的各种因素决定了凝固区域的宽窄,也就是决定了铸件的凝固方式。

例如合金的结晶间隔(液相线到固相线之间的温度间隔)越大,铸件的凝固区域就越宽,糊状凝固方式的倾向性就越大。反之,结晶间隔越小,则铸件的凝固区域越窄,逐层凝固方式的倾向越大。

凝固过程数值模拟

凝固过程数值模拟

左边第二项是枝晶间液体流动的热对流,由于潜热释放、两 相区温度梯度及液相率比较小,可忽略不计,式可修改为:
2015/11/13
初始条件:
冲型完毕时金属液和型腔的温度场即是凝固开始的温度场
边界条件:
边界条件是指铸件——铸型交界面或不同种类造型材料之间的 交接面。 界面热流率为:
如果界面换热系数为h,则界面两侧相邻节点之间的等效换热系 数为:
要解决的主要问题:
1.如何计算收缩量? 2.如何进行补缩通道的判断? 3.哪些单元可以进行补缩?
2015/11/13
多热节判断
多热节的判 断就是判断 出补缩通道 的情况,从 而确定出热 节的收缩量 和缩孔位置
2015/11/13
收缩量的计算
每个孤立热节内所有单元在第m时间步长内的体
积收缩量的计算如下:
2015/11/13
补缩单元
当单元的固相率 高于一个临界值 fsc时,就不再具 有流动能力,因 此也就不再具有 补缩能力
2015/11/13
补缩单元
(a)未考虑孤立熔池热节 (b)考虑孤立热节
2015/11/13
温度场模拟软件CAE
模拟的主要步骤:
(1)工艺三维建型 (2)新建工程 (3)前处理网格划分 (4)参数设定与技术分析 (5)后处理数据可视化
铸件凝固过程数值模拟
• 主 讲 人 :贺腾博 • 小组成员:黄萍、徐晓欢、贺腾博、韩丽 梅、刘靓

2011凝固第一章

2011凝固第一章

对于具体给定的铸件和铸型,用迭代法很容易求得φ值。
五、虚拟凝固层厚度s0
当x’=s’=s0时,即凝固刚开始,界面上的热平衡关系可表示为
s ' hi (T f T0 ) s L( ' ) s0 t
2 s 2 s ' 根据(1-24)得: ( ' )s0 t s0
所以
(1-33)
金属浇铸到金属型中,铸件降温很快,铸件和铸 型温度急速变化。对于许多钢锭、永久模铸造和压铸 来说都是金属模,因此有必要了解在金属型中影响凝 固因素的变化。 由于金属型具有很高的导热性能,所以在铸件凝固 过程中热流的限制环节通常不在铸型,而在铸件与铸型 之间的界面,当铸件凝固收缩和铸型受热膨涨而在铸 件铸型间形成气隙时,界面热阻的作用将变得更为突 出。
s m (TM T0 ) L s 型t t
(1-12)
dt 0 Ls ds (TM T0 ) mcm m / 0 t
s t
(1-13)
2 TM T0 S m m cm t s L 铸型
金属
(1-14)
16
对于形状简单的铸型,可以用V/A替代S, 可得
17
四、过热的影响
假定金属无过热,实际不可能。考虑过热,把 过热的热量加到潜热中去,得到总热量,凝固过程传 出的总热量。 H总=L+C (TP-TM) L—结晶潜热;J/kg

铸件凝固过程中热应力场及热裂的数值模拟研究分析.

铸件凝固过程中热应力场及热裂的数值模拟研究分析.

铸件凝固过程中热应力场及热裂的数值模拟研究分析

1 铸件凝固过程数值模拟的意义及概况

自1962年丹麦Fround第1个采用电子计算机模拟铸件凝固过程以来,计算机在铸造工艺研究中得到了广泛的应用,如凝固过程温度场、热应力场的数值模拟,充型过程流速场的数值模拟;组织形态及力学性能的数值模拟等。通过这些单1或复合过程的数值模拟,可以分析铸件中存在的各种缺陷的产生原因,进而采取相应工艺措施来消除缺陷,实现工艺优化,同时可以节省大量的人力、物力和财力,缩短产品从设计到应用的周期,增强产品的市场竞争能力。如今,在芬兰,90%以上的铸造厂在日常中应用铸造模拟软件辅助铸造工艺设计;世界上一些大型的汽车公司的铸造厂,如美国的通用、福特,德国的奔驰等,都把数值模拟软件作为1种日常工具来使用。

近10年来,涌现出了许多优秀的铸造过程数值模拟软件,如美国的ProCast、德国的MAGMASoft、芬兰的CastCAE、西班牙的ForCast、日本的CASTEM、法国的SIMULOR软件等。从功能上看,许多软件可以对砂型铸造、金属型铸造、精密铸造、压力铸造等多种工艺进行温度场、流场、应力场的数值模拟,可以预测铸件的缩孔、缩松、裂纹等缺陷和铸件各部位的组织。国内在经历了10多年的基础研究和发展后,也出现了一些技术水平接近国外商品化的应用软件,可以进行铸钢、铸铁件砂型铸造时的三维温度场模拟及收缩缺陷的预测,以及对铸钢、铝合金件的热应力场进行模拟。总的来说,国外软件的通用性强,能进行铸造全过程的数值模拟,并具有较强的后置处理功能及友好的用户界面。建模方便,易于模型设计和修改,便于用户掌握和使用。其计算精度与运算速度等方面也能满足需要。正因为如此,国外模拟软件已经成为实际生产中的有力工具.国内不少用户趋向于采用大型通用工程软件如:COSMOS、ANSYS、ADINA等进行模拟计算。

凝固过程温度场

凝固过程温度场

(二)数值方法
数值方法又叫数值分析法,是用计算机程序来求解数学模型的近似解(数值 解),又称为数值模拟或计算机模拟。 有限元法:有限元法是根据变分原理来求解热传导问题微分方程的一种数值 计算方法。有限元法的解题步骤是先将连续求解域分割为有限个单元组成的 离散化模型,再用变分原理将各单元内的热传导方程转化为等价的线性方程 组,最后求解全域内的总体合成矩阵。由于有限元法的单元形状可以比较随 意,因此更能适用于具有复杂形状的物体。
数值计算是数值模拟 技术的核心组成部分, 前后处理均为数值计 算服务
而数值计算就是以传 热学为基础来进行的。
后处理
传热学基础
定义:传热学(heat transfer)是研究热量传递规律的科学,是研究由 温差(推动力) 引起的热能传递规律的科学。
基础:热力学第一定律和第二定律 热力学第一定律:热量可以从一个物体传递到另一个物体,也可以与机 械能或其他能量互相转换,但是在转换过程中,能量的总值保持不变。 热力学第二定律:不可能把热从低温物体传到高温物体而不产生其他影 响,或不可能从单一热源取热使之完全转换为有用的功而不产生其他影 响,或不可逆热力过程中熵的微增量总是大于零。又称“熵增定律”, 表明了在自然过程中,一个孤立系统的总混乱度(即“熵”)不会减小。
A:垂直于导热方向的截面积[m2] 。 :导热系数(热导率)[W/mK]
一维无内热源单层平壁中

基于有限元的模拟挤压铸造凝固过程数学模型分析

基于有限元的模拟挤压铸造凝固过程数学模型分析

基于有限元的模拟挤压铸造凝固过程数学模型分析

1.引言

挤压铸造是一种重要的金属加工工艺,通过在金属液体受到一定的挤压后凝固成型,

可以得到具有一定几何形状和性能要求的零件。在挤压铸造过程中,金属液体的凝固过程

会受到多种因素的影响,如温度场、应力场以及固相的凝固行为等。建立准确的数学模型

对挤压铸造过程进行分析和模拟具有重要意义。

有限元方法是一种常用的数值计算方法,可以用于求解复杂的物理问题,包括金属凝

固过程的数学模型。本文将针对挤压铸造凝固过程,利用有限元方法进行数学模型的建立

和分析,对挤压铸造过程中的凝固行为进行深入研究。

2.挤压铸造凝固过程数学模型

挤压铸造凝固过程可以用热传导方程和固态相变方程进行描述。在进行数学建模时,

需要考虑金属液体在凝固过程中温度场和应力场的变化,并结合相变过程对凝固过程进行

描述。下面是对挤压铸造凝固过程进行数学建模的步骤:

2.1 热传导方程

考虑挤压铸造过程中金属材料的热传导行为,可以得到热传导方程如下:

\[ \frac{\partial \rho C_p}{\partial t} = \nabla \cdot (k \nabla T) \]

\( \rho \) 为金属的密度,\( C_p \) 为比热容,\( k \) 为热导率,\( T \) 为温度场。

通过以上两个方程的描述,可以得到挤压铸造过程中金属液体的凝固过程的数学模型,联立热传导方程和固态相变方程可以得到完整的挤压铸造凝固过程数学模型。

3.有限元模拟

利用有限元方法可以求解挤压铸造凝固过程的数学模型。有限元方法是一种离散化的

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

传热学基础
热量传递的三个基本方式
热传导 热对流 热辐射 conduction convection radiation
热传导( thermal conduction )
定义:指温度不同的物体各部分或温度不同的两物体间直接接触时,
依靠分子、原子及自由电子等微观粒子热运动而进行的热量传递现象。 When molecules collide, energy is transferred from the more
最大温度变化率为温度梯度(temperature
gradient).用grad t表示。
温度梯度
温度梯度(temperature gradient)是等温线面法线
方向上的温度变化率。 在温度场中,温度梯度表达了温度在空间上改变的大 小程度,是一个矢量。方向指向温度增大的方向
热流的方向与温度梯度方向相反。

通过单层无限大平壁的稳态导热,可视为一维稳态导热,边界条件可以为第 一类、第三类边界条件。这里仅讨论第一类边界条件。
稳态中单层平壁的导热
1、物理模型及数学模型(第一类边界条件): 其数学描述为:


d dT ( ) qv 0 导热微分方程式: dx dx
边界条件:
T x0 Tw1
energetic(high temperature) molecules to the less energetic
(lower temperature) molecules. 物质的属性:可以在固体、液体、气体中发生 导热的特点: 1. 必须有温差
2. 物体直接接触
3. 依靠分子、原子及自由电子等微观粒子热运动而传递热量
差分法: 差分法是把原来求解物体内随空间、时间连续分布的温度问题, 转化为求在时间领域和空间领域内有限个离散点的温度值问题,再用这些离
散点上的温度值去逼近连续的温度分布。差分法的解题基础是用差商来代替
微商,这样就将热传导微分方程转换为以节点温度为未知量的线性代数方程 组,得到各节点的数值解。
(二)数值方法
导热系数(Heat Conductivity)
非金属材料的导热机理:非金属物质多属于多孔性材料,其内部 孔隙部分充满着空气。其导热机理一般是通过材料的实体和孔隙 空气两部分热量传递综合作用的结果,如果空隙大到一定程度, 也会存在对流换热换热和辐射换热方式。 273K时物质的导热系数
热传导( thermal conduction ) 一、导热微分方程式的表达式
T x Tw2
2、求解方法
无内热源,且导热系数λ 为常数
d 2T 导热微分方程式: 0 dx2
温度分布: T
C1 x C2 其中 C1
(Tw1 Tw 2 )

C2 Tw1
热流通量及热流量: q
dT T T w1 w 2 dx
热阻及热阻分析图:R
第二类边界条件:给出通过物体表面的比热流随时间变化的关系:
第三类边界条件:给出物体周围介质温度以及物体表面与周围介质的换热系数:
稳态中单层平壁的导热
一、通过单层平壁的导热 • 无限大平壁的长度和宽度都远大于其厚度,因而平壁两侧保持均匀边界条件 的稳态导热就可以归纳为一维稳态导热问题。从无限大平壁的结构可分为单 层壁,多层壁和复合壁等类型,如图所示。

缺点:通常需要采用多种简化假设,而这些假设往往并不适合实际情况,
这就使解的精确程度受到不同程度的影响。目前,只有简单的一维温度 场(“半无限大”平板、圆柱体、球体)才可能获得解析解。
(二)数值方法
数值方法又叫数值分析法,是用计算机程序来求解数学模型的近似解(数值
解),又称为数值模拟或计算机模拟。
3分析求解,得出导热物体的温度场。
4利用傅立叶定律和已有的温度场最终确定热流量或热流密度。
对具体热场进行求解时,需要根据具体问题各出导热体的初始条件和边界条件。 初始条件:物体开始导热时(t=0时)的瞬时温度分布。 边界条件:导热体表面与周围介质间的热交换情况。 常见的边界条件有以下三类: 第一类边界条件:给定物体表面温度随时间的变化关系:
凝固过程温度场相关
报告人:陆

温度场
1、基本概念 指某一瞬时物体内各点的温度分布状态。温度是标量,温度场是时间和空间 的函数,也是标量场。 在直角坐标系中: 在柱坐标系中: 在球坐标系中:
根据温度场表达式,可分析出导热过程是几维、稳态或非稳态的现象,温
度场是几维的、稳态的或非稳态的。 例如表示导热过程是二维、稳态的导热现象,温度仅在x、y方向发生变化, 但不随时间变化; 表示导热过程是一维、非稳态的导热现象,温度仅在x方向随时间发生变化。
射能。热辐射的主体与受体是相对的,辐射能的传递是相互往
复发生的,一定时间后双方的辐射速度趋于等同,便出现暂时 的热平衡。 热辐射是物体因自身的温度而具有向外发射能量的本领。热辐射虽然也是热传递
的一种方式,但它和热传导、对流不同。它能不依靠媒质把热量直接从一个系统
传给另一系统。热辐射以电磁辐射的形式发出能量,温度越高,辐射越强。辐射 的波长分布情况也随温度而变,如温度较低时,主要以不可见的红外光进行辐射, 在500℃以至更高的温度时,则顺次发射可见光以至紫外辐射。热辐射是远距离传 热的主要方式,如太阳的热量就是以热辐射的形式,经过宇宙空间再传给地球的。
铸造过程数值模拟 三维实体造型 铸造工艺设计 前处理 而数值计算就是以传 热学为基础来进行的。 数值计算是数值模拟 技术的核心组成部分, 前后处理均为数值计 算服务
网格剖分 数值计算 结果显示
后处理
传热学基础
定义:传热学(heat transfer)是研究热量传递规律的科学,是研究由 温差(推动力) 引起的热能传递规律的科学。
组,最后求解全域内的总体合成矩阵。由于有限元法的单元形状可以比较随 意,因此更能适用于具有复杂形状的物体。
二、导热理论分析方法的基本思路 •导热理论的任务就是要找出任何时刻物体中各处的温度,进而确定热量传递 规律。
1简化分析导热现象,根据几何条件、物理条件简化导热微分方程式。
2确定初始条件及各物体各边界处的边界条件,每一维导热至少有两个边界条 件。从而得到导热现象的完整数学描述,包括:导热微分方程式和单值性条件 (见图)。
数值方法又叫数值分析法,是用计算机程序来求解数学模型的近似解(数值
解),又称为数值模拟或计算机模拟。
有限元法:有限元法是根据变分原理来求解热传导问题微分方程的一种数值 计算方法。有限元法的解题步骤是先将连续求解域分割为有限个单元组成的
离散化模型,再用变分原理将各单元内的热传导方程转化为等价的线性方程
热传导( thermal conduction )
导热的基本定律:
1822年,法国数学家Fourier: 上式称为傅立叶定律(导热基本定律), 是一个一维稳态导热。其中: -: 热量传递的方向与温度梯度方向相反。 Q:热流量,单位时间传递的热量。[W] q:热流密度,单位时间通过单位面积传递的热量[W/ m2]
A:垂直于导热方向的截面积[m2] 。
:导热系数(热导率)[W/mK]
一维无内热源单层平壁中
热传导( thermal conduction )
傅立叶定律(Fourier's Law)
傅里叶定律的文字表述:在导热现象中,单位时间内通过给定截面的热量, 正比例于垂直于该截面方向上的温度变化率和截面面积,而热量传递的方 向则与温度升高的方向相反。 傅立叶定律是热传导的基础。它并不是由热力学第一定律导出的数学表达 式,而是基于实验结果的归纳总结,是一个经验公式。同时,傅立叶定律 是定义材料的一个关键物性,热导率的一个表达式。 另外,如上所述,傅立叶定律是一个向量表达式。热流密度是垂直于等温 面的,并且是沿着温度降低的方向。傅立叶定律适用于所有物质,不管它 处于什么状态(固体、液体或者气体)。
会带着他的热量升上去,凉的气体或液体会降下来替代刚
升上去的气或热的位置,然后受热后继续升上去,刚上去 的丢失了热量后会降下来,这样反复循环,就是对流。就 像平时烧水,先是下面的水受热,然后升上去,上面凉的 水就会降下来然后受热。
热辐射( thermal radiation ) 定义:热辐射是由于物体内部原子振动而发出的一种电磁波的 能量传递。一切自身温度高于0K的物体,都会从表面发射出辐
直角坐标系:
圆柱坐标系:
球坐标系:
热传导( thermal conwk.baidu.comuction ) 一、导热微分方程式的表达式
导热微分方程式一般由导热项、内热源生成项及非稳态项组成。 如图所示。
内能的热增加率 (非稳态项)
导热的净热量 (导热项)
内热源 (内热源生成项)
qv-内热源;-物体的密度; -导热系数;t-温度;τ -时间
热对流( Heat convection )
定义:由流体各质点间的相对位移而引起的热量转移方式称为热对流。
对流包括自然对流和强迫对流。自然对流是由于质点间的温度差或者 密度差引起的浮力流,强迫对流是体系在外力(如机械力、电磁力等)
驱动下产生的质点的相对位移
热对流一般是发生在气体和液体中的,受热的气体或液体
温度场
2、等温面与等温线
三维物体内同一时刻所有温度相同的点的集合称为等 温面(isothermal surface); 一个平面与三维物体等温面相交所得的的曲线线条即 为平面温度场中的等温线(isotherms)。
3、温度梯度
在具有连续温度场的物体内,过任意一点P温度变化 率最大的方向位于等温线的法线方向上。称过点P的
Tw1 Tw2 q
规律:温度分布为直线且斜率大小由导热系数决定;内部各处热流通量 及热流量处处相等;
有内热源,且导热系数λ 为常数 •
2 d 导热微分方程式: T qv dx2

(Tw1 Tw 2 ) qv 2 qv 2 C C2 Tw1 T x C1 x C2;其中 1 温度分布: 2 2
基础:热力学第一定律和第二定律 热力学第一定律:热量可以从一个物体传递到另一个物体,也可以与机 械能或其他能量互相转换,但是在转换过程中,能量的总值保持不变。 热力学第二定律:不可能把热从低温物体传到高温物体而不产生其他影 响,或不可能从单一热源取热使之完全转换为有用的功而不产生其他影 响,或不可逆热力过程中熵的微增量总是大于零。又称“熵增定律”, 表明了在自然过程中,一个孤立系统的总混乱度(即“熵”)不会减小。 热量可以自发地从温度高的物体向温度低的物体进行传递,如果没有能 量转化的途径,热量始终是守恒的。
导热系数(Heat Conductivity) 一、导热系数 • • 定义式:
导热系数在数值上等于单位温度降度(即lK/m)下,在垂直于热流密度的单位 面积上所传导的热流量。导热系数是表征物质导热能力强弱的一个物性参数 。 二、影响因素 • 包括:物质的种类及性质、温度、压力、密度以及湿度 • 各种物质的导热系数相差很大,其根本原因在于不同的物质其导热机理存在 着差异。一般而言,金属的导热系数最大,非金属和液体次之,气体的导热 系数最小。导热系数越大,说明其导热性能越好。由图中可以看出,各类物 质导热系数的一般大小顺序。
导热系数(Heat Conductivity)
现行国家标准(GB 4272—92)规定,平均温度在350℃以下时 导热系数低于0.12时,这种材料称为保温材料。
导热系数(Heat Conductivity)
同一种物质的导热系数也会因其状态的不同而改变,因而导热系数是物质温 度和压力的函数。由于物质温度和压力的高低直接反映物质分子的密集程度 和热运动的强弱程度,直接影响着分子的碰撞、晶格的振动和电子的漂移, 故物质的导热系数与温度和压力密切相关。见下表。
内热源项——结晶潜热 结晶潜热(latent heat of crystallization)是指在温度保持不变的情 况下,单位质量的物质从液态转变到固态时所释放出的热量。
温度
温度
时间
时间
凝固温度场的求解方法
(一) 解析法 (二) 数值方法
(一) 解析法 • • 解析方法是直接应用现有的数学理论和定律去推导和演绎数学方程(或 模型),得到用函数形式表示的解,也就是解析解。 优点:是物理概念及逻辑推理清楚,解的函数表达式能够清楚地表达温 度场的各种影响因素,有利于直观分析各参数变化对温度高低的影响。
相关文档
最新文档