基于ANSYS软件的钢球淬火过程分析.docx
- 文档编号:2366538
- 上传时间:2023-05-03
- 格式:DOCX
- 页数:13
- 大小:250.07KB
基于ANSYS软件的钢球淬火过程分析.docx
《基于ANSYS软件的钢球淬火过程分析.docx》由会员分享,可在线阅读,更多相关《基于ANSYS软件的钢球淬火过程分析.docx(13页珍藏版)》请在冰点文库上搜索。
基于ANSYS软件的钢球淬火过程分析
H1NGXIAUNIVERSITY
举业诈文
(2015届)
题目基于ANSYS软件的钢球淬
火过程分析
学院数学计算机学院
专业信息与计算科学
年级201]级
学生学号
学生姓名杨坤指导教师马文涛
基于ANSYS软件的钢球淬火过程分析
数学计算机学院信息与计算科学专业2015届杨坤
摘要:
钢球广泛应用于工业生产中的多个领域,钢球的淬火处理是保证钢球质量的关键生产工艺之一.ANSYS有限元软件在温度场的模拟过程中,很好地结合了材料变温过程材料热物性参数的变化.特别适用于钢球淬火过程温度场的准确计算.本文利用ANSYS有限元分析软件对钢球淬火过程温度场进行有限元模拟分析,得到了钢球温度随淬火时间变化的分布关系和钢球体内不同位置的温度分布情况.
关键词:
钢球淬火;温度场;ANSYS软件;模拟
中图分类号:
AnalysisoftheSteelBallQuenchingProcessBasedonANSYSSoftware
Abstract:
Steelballiswidelyusedintheindustrialproduction,andhardenedsteelballisthekeytoensurethequalityoftheproductionprocessofsteelball・NonlinearthermophysicsparametersareconsideredinANSYSfiniteelementsoftware,soitcanbeusedtocomputetemperatureinthequenching・InthispapermakesuseofANSYSfiniteelementsoftwaretosimulatethetemperaturefieldinquenchingprocessofsteelball,andobtainstherelationshipbetweentemperatureofsteelballandquenchingtimeandtemperaturedistributionofdifferentposition.
Keywords:
steelballquenching;temperaturefield;ANSYSsoftware;simulation
1
1
1
2
2
2
3
4
4
4
5
6
6
6
6
8
10
错误!
未定义书签。
错误!
未定义书签。
基于ANSYS软件的钢球淬火过程分析
钢球催火的基本原理
钢球淬火的模拟方法
2淬火过程数学模型
温度场控制方程
淬火过程导热微分方程
导热微分方程的定解条件
钢球淬火过程有限分析
淬火模型求解
模型算例验证
有限元分析
3基于AXSYS软件的温度场分析
非稳态温度场实例
钢球淬火的ANSYS模型
前处理阶段
后处理阶段
4结术语
参考文献
致谢
基于ANSYS软件的钢球淬火过程分析
1引言
材料及软件介绍
(1)钢球主要用来制造轴承和粉碎物料等,有着较广泛的用途.钢球一般采用铸造或者山钢锭锻造的方式生产.生产过程中,钢球的淬火处理往往是其生产过程中的关键工艺之一.钢球的淬火处理一般采用把常温的钢球重新加热快速冷却的方式进行.重新加热淬火时,淬火温度容易控制,钢球的部分碳化物溶解,碳向基体扩散,可以有效地提高钢球的表面硬度.传统粗放式生产方式常采用对一批(筐)钢球进行整体加热后一起倒入水中冷却的方式进行.
(2)ANSYS软件是一个功能强大而灵活的大型通用有限元分析软件.John-Swanson博士与1970年开发的,经过儿年对原始软件修改,LI前,ANSYS软件的功能更加强大,使用更加便捷.新版本ANSYS软件中可以用来结构静力学分析、结构动力学分析、结构屈曲分析、热力学分析、电磁场分析、声场分析、压电分析、流体动力学分析.同时进行产品分析时,使用ANSYS软件对产品性能进行仿真分析,可以发现产品的问题,降低产品成本,缩短设计周期,提高设计的成功率.
钢球催火的基本原理
淬火⑴是把钢球加热到临界值以上,保温一定的时间,然后以大于临界冷却速度,从而获以马氏体为主的不平衡组织的一种热处理方法.
在淬火冷却过程中,因为钢球内部温度分布不均匀,组织转变过程的不均匀而形成热应力和相变应力,这些应力的存在将直接影响钢球的组织性能和使用寿命.如果热处理不当,将会造成钢球组织性能达不到预定要求,甚至会产生过量变形或开裂而报废.
钢球淬火的模拟方法
U前普遍采用比较成熟的有限元法来模拟淬火过程的变化.所使用的有限元软件有ADINA,NSYS,MARC等.大多数的模拟结果都得到了实验测试数据的支持,取得了令人满意的效果.汁算机模拟技术是发展材料学与工程技术的桥梁.热处理计算机模拟有赖于讣算传热学、相变动力学、热弹塑性力学、数值计算等学科的发展.它不仅为热处理的计算机模拟研究与应用提供了坚实的基础,同时它乂促进了热处理模拟技术的迅速发展.热处理模拟技术不仅有强大的计算功能,而且还可以显示任何时刻工件内任意截面上的温度场、组织场、应力场,能使操作者观察到各种场量随时间推移的情况,也可以显示出用户所感兴趣的任何点上的温度、应力、组织随时间的变化曲线.
2淬火过程数学模型
温度场控制方程
固体热传导的控制方程是Fourier导热方程,最基本的Fourier定律在稳态条件下可描
述为:
dT
4=-^—
ox
其中,你是x方向上单位表面积的热流率,称为比热流率或热流密度是材料的
导热系数(w/(加•©);<是x方向上的温度梯度「c/加);“-”号表示传热的方向和温
Ox
度梯度的方向相反.
淬火过程导热微分方程
在淬火过程中,钢球内部的热量是以热传导方式向外传递的,在钢球表面,主要是通过淬火介质与钢球表面的热对流而散发热量根据能量守恒定律和傅里叶定律,建立导热物体中温度场的数学表达式[2】,式
(2)为三维的钢球热传导偏微分方程:
oT
山于钢球是中心对称物体,则钢球沿■儿z方向的导热系数几是相同的,所以
(2)式可写
18(/7)
Idt
为:
d(dT+———
dz\dz,))2
方程式中,等式左边项的为钢球体内升温需要的热量;等式右边项分别为■y和z方向流入钢球体内的热量;最后一项为钢球体内产生的热量.其中,°是材料的密度(畑/亦);c为材料的比热容(丿/(畑・K));f为时间($);心人虫分别为材料沿匚儿z方向的热导率(W/(〃?
・K));0=°(x,y,乙/)为材料内部的热源密度(VV/kg).
导热微分方程的定解条件
导热微分方程的定解条件包括初始条件和边界条件【2】.
(1)初始条件.初始条件是指求解问题的初始温度场,是计算的出发点,也就是在零时刻温度场分布.其各点的温度值已知或者遵从某种函数关系,如式(4)所示.
T\ls0=T0(W)
山于淬火冷却过程模拟的初始条件一般都忽略加热过程,即假定淬火前的加热是均匀透烧,初始温度场为均匀温度场,初始的应力为零.所以上述函数关系式可以变为式(5).
TU=7;=900oC
其中,7;为已知常温,是常数.
(2)边界条件.边界条件是求解温度场最要的条件,也是指物体表面或者边界与周圉环境的热交情况,通常有三类边界条件.不同边界条件下的导热情况及表达式如下:
1第一类边界条件是指物体边界上的温度函数已知,表示为
2笫二类边界条件是指物体边界上的温度函数已知,表示为
式中,"为物体边界的外法线方向.
3第三类边界乂称为对流边界条件,是指物体与周围环境介质间的对流转热系数和介质的温度已知,表示为
式中,T是物体表面温度;7;为介质温度;k为物体表面换热系数.
根据边界条件可知,钢球淬火一般用第三类边界条件,山给出的数据上式可以写为:
-A^-=k(T-Tf)=1OO(T-2O)
钢球淬火过程有限分析
淬火模型求解
影响钢球淬火的物理学因素很多.比如:
物体密度.物体比热容、物体热导率以及初始温度等.山于钢球的体积相对较小,球内各点的温度变化可以忽略不计.根据能量守恒定律,我们可以得到温度随时间变化的表达式:
此外,山于钢球体和周用流体之间存在热对流,因而满足固体热传导的控制方程.而钢球冷却时流体的温度不随时间变化,那么(10)式中的7;为常数,并且方程是一个一阶微分方程,然后分离变量并积分.
钢球达到一定温度所需要的时间为:
hAT-Tf
钢球在某一时间内的温度为:
模型算例验证
首先,计算毕奥数看是否可以应用集中容量法求解,钢球的体积与表面积之比,既
有:
v-nr
Vqr
L=—==—=0.05//Z
”A4力・23
Bi-仏=
solid
=(嘶加•K)x(0.05叭。
“429
7OW/m-k
因为显然可以使用集中容量法求解.现利用方程确定30$后钢球的温度为:
T—Tf
—•Qvr\
(hA'f
r-20
——avr\
/(100W/加2・k)x(50“、
CAU 900-20r 、(7800畑/m)x(448J/畑•K)x(0.05加)) =>r=741.17°C 说明: 相关数据第三节给出. 有限元分析 根据钢球淬火热传导的控制函数和偏微分方程,以及已知的边界条件和初始条件可知,通常称它为初边值问题.这类问题的求解特点是在空间域内用有限单元网格划分,而在时间域内则用有限差分网格划分,实质上是有限单元和有限差分的混合解法.山时间两点差分(时间向前差分,空间中心差分)可得: 昭-27;: +粘T=-2T: +T二7;二-27;+尤=丄巧-打x2y2z22t 化简,有: 仃二+T二"二+兀J二+T”J 4. \—7t2 fiiiY 4-亠I •>十十•> I2)r疋丿 1 1厶fZu 2十r十2 l疋: r匸丿」 T: 于是,有: 环)+1-2/2 1 1 亠 卜丄) * 令Fo=tX 得: 矿=局(尤+為)+[1-2尸叩;: 稳定收敛标准为: 3基于ANSYS软件的温度场分析 非稳态温度场实例 有一个直径为的钢球,球的初始温度为900°C,将其突然放入温度为20°C的水中,表面传热系数100VW(m2-°C)的流体介质中,放入水中时间为Is.钢球材料性能如下: 密度: 7800Kg/m\热导率: 70网(m・°C)、比热容: 448J/(Kg-°C). 计算: (i)第Is,10s,30s和第50s的整个钢球的温度分布. (ii)钢球边界上任一点的温度在50s内的变化曲线图. 钢球淬火的ANSYS模型 前处理阶段 (1)定义单元类型.在ANSYS单元库中有一百多种不同的单元类型,每一个单元类型都具有一个惟一的编号和一个标识单元类别的前缀,如BEAM4、PLANE77、SOL工D96等[3] 1确定自山度是否相容.根据自山度的不同可选择的单元有线、面或三维实体的单元类型. 2在淬火冷却温度场分析中,首先,单元必须要能够进行热分析;其次,所选的单元还应能够进行瞬态非线性分析;最后还应能进行工件内部应力分析. 根据这些要求,本文选取三维的SOLIDS.SOL工DS具有3维磁场、热场、电场、压电场和结构场分析能力,并能在各场之间实现有限的藕合. (2)定义材料性能参数.钢球的密度与温度变化关系不是很大,按常数处理,对最后结果影响不大.对比热容、热导率、换热系数等这些温度影响较大的热物性参数,必须根据其对应关系进行输入,以确保讣算结果的准确性. (3)建立儿何模型.有限元分析是对物理现象的模拟,是对真实情况的数值近似.通过对分析对象划分网格,求解有限个数值来近似模拟真是环境的无限个未知数. ANSYS中的模型是由ANSYS自带的建模工具进行建模.一般来说,比较复杂的模型可以用其它的CAD软件较为方便.本文的研究对象为钢球,结构简单,用ANSYS直接建模.钢球截面如图3・1所示. (4)网格划分.ANSYS中划分网格的方式主要有两种,分别为自山网格划分和映射网格划分.自山划分的网格对单元形状没有任何限制,用这种方式划分的网格排列不规则;而映射网格划分对包含的单元形状有限制,用映射网格划分方式得到的网格具有规则的儿何形状,而且它对载荷的施加和收敛的控制相当有利.因而,在实际应用中一般优先选用映射网格划分.本文釆用自山网格划分,建立实体模型并划分网格后如图3・2・ 图3-2网格划分图 (5)加载荷载求解.ANSYS软件的分析类型分为: ・Static静态分析 •Modal模态分析 •Harmonic简谐激励分析 •Transient瞬态激励分析 •Spectrum响应谱分析 •EigenBuekling线性屈曲分析 由于淬冷过程中有相变发生,所以选取瞬态激励分析,即Transient.另在模型上加载初始温度为900°C(淬火开始时温度),模型边界上施加温度载荷20°C(淬火介质的温度).再设置载荷包括载荷步长、子步等后,设置记录ANSYS运算过程中每一步迭代运算结果,以便后处理时查看淬火中间过程.设置完后,就可以运算求解了. 后处理阶段 (1)III设定的求解条件可以得到淬火过程50s内零件各部分的温度场分布,并可模拟温度场随时间的动态变化过程.下图所示为零件表面及心部在淬火Is,10s,30s,50s时的温度场的分布.从图中可以看出,在淬火Is时,表面温度迅速降低,而心部温度仍然维持高温,因而表面淬透性很大,这与实际情况相符合.随着淬火时间的延长,表面温度持续下降,但是温度降低速率显着下降叫 图3・3t=\s和/=10$时钢球体内温度分布图 POST26 TEMP^3 TEX匸4 图3・4t=305和/=50$时钢球体内温度分布图 (2)根据节的差分格式计算,在同一时刻各节点的温度分布情况与ANSYS有限元模拟的数据结果大致相同,说明模拟结果具有可靠性.如图3-5所示对6、3、4结点进行研究,考察其温度随时间变化情况.6节点山于处于轮轴外边缘交线上,因而在刚开始淬火时,温度下降极为迅速.而3节点由于处在钢球的中间部位,由于内部热源的补充使其温降速度下降缓慢.同样的道理,3节点与4节点的温度时间关系也存在较大的差别.同样可以对于其感兴趣的点,可利用该方法求出其在任何淬火时它任何刻的准确温度值. AN APR132015 09: 33: 00 STEELBALLA^CiLSISOFAWHEEL 图3・5表示为温度随时间的变化曲线 4结术语 淬火过程是一个复杂的物理变化过程,在这一变化过程中,钢球的温度、组织和应力都会发生变化,并且他们之间相互作用、相互祸合.本文以大型有限元软件ANSYS为平台,钢球以水火介质的淬冷过程的温度场、热应力场及残余应力场进行了数值模拟计算,主要得出以下结论: (1)对于儿何外形复朵的钢球,应用ANSYS求解其淬火过程温度场,计算速度较快,并且计算结果可靠. (2)对于热物性参数为非线性的材料,特别适合利用ANSYS进行求解分析. (3)利用ANSYS进行淬火过程温度场分析,可以为淬火过程热应力,残余应力的计算准备温度条件. (4)钢球在水为介质的淬冷条件下,不管是高温区还是在低温区冷却速度都很大,这对钢球淬冷后性能特别是抗拉强度不利: 但在低温区冷却速度比较低,这有利于工件淬冷后抗拉强度的提高. 参考文献 [1]商跃进.有限元原理与ANSYS应用指南|M].北京: 清华大学出版.2005. [2]胡红军.ANSYS10.0料工程有限元分析实例教程[M].北京机械工^出版.200& ⑶王建江胡仁B.ANSYS结构与热力学有限元分析实例指导教程[M].北京机械工业出版.200& [4]赖宏•刘天模•45刚零件淬火过程温度场的ANSYS模拟[J|.重庆大学学报.2003,26(3): 82-84. 致谢 本文是在导师马文涛副教授的悉心指导下完成的.导师渊博的学识、严谨的治学态度和认真负责的工作作风,使我获益匪浅,并将受益终生.在此论文即将完成之际,向我的导师马老师致以崇高的敬意和衷心的感谢. 感谢宇夏大学数学计算机学院各位老师在整个课题的进行过程中,给予我的亲切关 怀和指导,在此向你们表示深深的谢意和诚挚的祝福.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 ANSYS 软件 淬火 过程 分析