魔神 发表于 2009-11-23 22:53:20

计算机仿真技术在磨削温度场中的应用

&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class=px14>&nbsp;<H1><FONT id=FontSizeSettings5><FONT size=2>1 引言</FONT></FONT></H1><DD><FONT id=FontSizeSettings5>磨削是一种应用广泛的精密加工方法。在对磨削加工的研究中,由于对加工机理的认识有限,因此对磨削加工过程的实际调整多数是靠试凑法(即凭操作者所积累的大量经验知识)来进行,特别是有关磨削温度分析模型,多是通过单因素获得的。随着计算机性能的日益提高,仿真技术在工业中的应用越来越广泛,给磨削理论的研究带来了新的思路,使我们有可能克服传统研究方法的局限性,深入研究磨削过程中磨削温度的变化,建立系统的磨削温度场理论模型。 </FONT><DD><FONT id=FontSizeSettings5>&nbsp; </FONT><DD><FONT id=FontSizeSettings5>仿真就是在模拟环境下实现和预测产品在真实环境下的性能和特征(动态和静态),它包含了从建模、施加负载和约束到预测产品在真实情况下的响应等一系列步骤。通过对仿真试验过程的观察和估计,得到被仿真系统的仿真输出参数和基本特性,由此来估计和推断实际系统的真实参数和真实性能。仿真技术借助计算机,可以在复杂的磨削加工过程中得到不同输入参数下的各种磨削温度场的变化,从而为深入研究磨削加工机理创造了条件。 </FONT><H1><FONT id=FontSizeSettings5><FONT size=2>2 建立磨削温度场的数学模型</FONT></FONT></H1><DD><FONT id=FontSizeSettings5>采用有限元法建立磨削温度场的数学模型。由于整个磨削温度场满足能量守恒定律,所以磨削温度场的热传导方程式为:<BR><TABLE width=550><TBODY><TR><TD align=middle><TABLE><TBODY><TR align=middle><TD rowSpan=3><FONT size=2>PC</FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>-</FONT></TD><TD><FONT size=2>∂</FONT></TD><TD rowSpan=3><FONT size=2>(k<SUB>x</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>)-</FONT></TD><TD><FONT size=2>∂</FONT></TD><TD rowSpan=3><FONT size=2>(k<SUB>y</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>)-</FONT></TD><TD><FONT size=2>∂</FONT></TD><TD rowSpan=3><FONT size=2>(k<SUB>z</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>)-<FONT face=symbol>r</FONT>Q=0</FONT></TD></TR><TR bgColor=#006600 height=1><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD></TR><TR align=middle><TD><FONT size=2>∂t</FONT></TD><TD><FONT size=2>∂x</FONT></TD><TD><FONT size=2>∂x</FONT></TD><TD><FONT size=2>∂y</FONT></TD><TD><FONT size=2>∂y</FONT></TD><TD><FONT size=2>∂z</FONT></TD><TD><FONT size=2>∂z</FONT></TD></TR></TBODY></TABLE></TD><TD width=80><FONT size=2>(在<FONT face=symbol>W</FONT>内)</FONT></TD></TR></TBODY></TABLE></FONT><DD><FONT id=FontSizeSettings5>其中<FONT face=symbol>W</FONT> 为整个域,它由三类边界条件组成,即: <TABLE width=550><TBODY><TR><TD align=middle><FONT size=2><FONT face=symbol>q</FONT>=<FONT face=symbol>q</FONT>( l,t)</FONT></TD><TD width=80><FONT size=2>(在l<SUB>1</SUB>边界上)</FONT></TD></TR><TR><TD align=middle><TABLE><TBODY><TR align=middle><TD rowSpan=3><FONT size=2>k<SUB>x</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>n<SUB>x</SUB>+k<SUB>y</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>n<SUB>y</SUB>+k<SUB>z</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD><FONT size=2></FONT></TD><TD rowSpan=3><FONT size=2>n<SUB>z</SUB>=q</FONT></TD></TR><TR bgColor=#006600 height=1><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD></TR><TR align=middle><TD><FONT size=2>∂x</FONT></TD><TD><FONT size=2>∂y</FONT></TD><TD><FONT size=2>∂z</FONT></TD></TR></TBODY></TABLE></TD><TD width=80><FONT size=2>(在<FONT face=symbol>G</FONT><SUB>2</SUB>边界上)</FONT></TD></TR><TR><TD align=middle><TABLE><TBODY><TR align=middle><TD rowSpan=3><FONT size=2>k<SUB>x</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>n<SUB>x</SUB>+k<SUB>y</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD rowSpan=3><FONT size=2>n<SUB>y</SUB>+k<SUB>z</SUB></FONT></TD><TD><FONT size=2>∂<FONT face=symbol>q</FONT></FONT></TD><TD><FONT size=2></FONT></TD><TD rowSpan=3><FONT size=2>n<SUB>z</SUB>=a(<FONT face=symbol>q</FONT><SUB>a</SUB> -<FONT face=symbol>q</FONT>)</FONT></TD></TR><TR bgColor=#006600 height=1><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD><TD><FONT size=2></FONT></TD></TR><TR align=middle><TD><FONT size=2>∂x</FONT></TD><TD><FONT size=2>∂y</FONT></TD><TD><FONT size=2>∂z</FONT></TD></TR></TBODY></TABLE></TD><TD width=80><FONT size=2>(在<FONT face=symbol>G</FONT><SUB>3</SUB>边界上)</FONT></TD></TR></TBODY></TABLE></FONT><DD><FONT id=FontSizeSettings5>根据有限单元法的离散化原理,把工件划分成有限个单元,并把磨削过程中的热负载加于边界上的单元,即把整体温度载荷离散为与实际载荷等效的节点载荷,代入上述三类边界条件即可得到磨削温度场的有限元模型。 </FONT><DD><FONT id=FontSizeSettings5>磨削温度场的场函数温度既是空间域<FONT face=symbol>W</FONT>的函数,又是时间域t的函数。在仿真模型中,采用循环迭代的方法,把热载荷的加载过程离散为有限个极短的过程(在单个过程中认为时间域和空间域是耦合的)。以干磨为例,在极短的时间内在某一磨削区加载一固定热流,在下一时间段内移至另一区域加载固定热源,并把上一次所得的计算结果作为本次分析的初始条件,在这一极短的时间内,认为场函数温度与时间和空间有关,但两者并不耦合,这样就可采用部分离散的方法得出如下所示磨削温度场的数学模型:<BR></FONT><CENTER><P align=center><FONT id=FontSizeSettings5><IMG style="BORDER-RIGHT: rgb(0,0,0) 1px solid; BORDER-TOP: rgb(0,0,0) 1px solid; BORDER-LEFT: rgb(0,0,0) 1px solid; BORDER-BOTTOM: rgb(0,0,0) 1px solid" alt=1 src="http://www.chmcw.com/upload_files/article/20/1_lqlwyp20071023072815232.GIF" border=0></FONT></P></CENTER><FONT id=FontSizeSettings5>式中:<FONT face=symbol>r</FONT>——材料密度 </FONT><DD><FONT id=FontSizeSettings5>c——材料比热 </FONT><DD><FONT id=FontSizeSettings5>t——时间 </FONT><DD><FONT id=FontSizeSettings5>k<SUB>x</SUB>,k<SUB>y</SUB>,k<SUB>z</SUB>——材料沿x,y,z 方向的热传导系数 </FONT><DD><FONT id=FontSizeSettings5><FONT face=symbol>q</FONT>——<FONT face=symbol>G</FONT><SUB>1</SUB>边界上的给定温度,<FONT face=symbol>q</FONT>=<FONT face=symbol>q</FONT>(<FONT face=symbol>G</FONT>,t) </FONT><DD><FONT id=FontSizeSettings5>q——<FONT face=symbol>G</FONT><SUB>2</SUB>边界上的给定热流量,q=q(<FONT face=symbol>G</FONT>,t) </FONT><DD><FONT id=FontSizeSettings5><FONT face=symbol>a</FONT>——对流换热系数<FONT face=symbol>q</FONT><SUB><FONT face=symbol>a</FONT></SUB>——在自然对流条件下,为外界环境温度;在强迫对流条件下,为边界层的绝热壁温度。<FONT face=symbol>q</FONT><SUB><FONT face=symbol>a</FONT></SUB>=<FONT face=symbol>q</FONT><SUB><FONT face=symbol>a</FONT></SUB>(l,t) </FONT><DD><FONT id=FontSizeSettings5>上式将时间域和空间域的偏微分方程问题在空间域内离散为n个节点为温度<FONT face=symbol>q</FONT>(t)的常微分方程的初值问题。求解这n个线性常微分方程,得到节点温度列阵,即可得出整个工件的温度分布。该公式适用于干磨和湿磨两种情况。 </FONT><H1><FONT id=FontSizeSettings5><FONT size=2>3 磨削温度场的仿真过程及结果</FONT></FONT></H1><DD><FONT id=FontSizeSettings5>以平面磨削温度场为例,建立磨削温度场的仿真模型。如图1 所示,假设磨削接触弧区为ABB'A',面热源长度l≈(d<SUB>s</SUB>a<SUB>p</SUB>)<SUP>&frac12;</SUP>且平行于工件运动方向。 <TABLE align=center><TBODY><TR><TD align=middle width=300><FONT size=2></FONT><P align=center><FONT size=2><IMG style="BORDER-RIGHT: rgb(0,0,0) 1px solid; BORDER-TOP: rgb(0,0,0) 1px solid; BORDER-LEFT: rgb(0,0,0) 1px solid; BORDER-BOTTOM: rgb(0,0,0) 1px solid" alt=1 src="http://www.chmcw.com/upload_files/article/20/1_kar33s20071023072815767.GIF" border=0><BR><B>图1 工件表面面热源示意图</B></FONT></P></TD><TD align=middle width=300 rowSpan=3><FONT size=2></FONT><P align=center><FONT size=2><IMG style="BORDER-RIGHT: rgb(0,0,0) 1px solid; BORDER-TOP: rgb(0,0,0) 1px solid; BORDER-LEFT: rgb(0,0,0) 1px solid; BORDER-BOTTOM: rgb(0,0,0) 1px solid" alt=q src="http://www.chmcw.com/upload_files/article/20/1_umvrfs20071023072815866.GIF" border=0><BR><B>图3 磨削温度场的仿真过程示意图</B></FONT></P></TD></TR><TR><TD align=middle width=300><FONT size=2></FONT><P align=center><FONT size=2><IMG style="BORDER-RIGHT: rgb(0,0,0) 1px solid; BORDER-TOP: rgb(0,0,0) 1px solid; BORDER-LEFT: rgb(0,0,0) 1px solid; BORDER-BOTTOM: rgb(0,0,0) 1px solid" alt=q src="http://www.chmcw.com/upload_files/article/20/1_gt8tb820071023072816179.JPG" border=0><BR><B>图2 干磨时TC4 工件的磨削温度场示意图</B></FONT></P></TD></TR><TR><TD align=middle width=300><FONT size=2></FONT><P align=center><FONT size=2><IMG style="BORDER-RIGHT: rgb(0,0,0) 1px solid; BORDER-TOP: rgb(0,0,0) 1px solid; BORDER-LEFT: rgb(0,0,0) 1px solid; BORDER-BOTTOM: rgb(0,0,0) 1px solid" alt=1 src="http://www.chmcw.com/upload_files/article/20/1_zjaaa520071023072816239.GIF" border=0><BR><B>图4 表面层下不同深度处干磨和湿磨仿真与试验结果的比较</B></FONT></P></TD></TR></TBODY></TABLE></FONT><DD><FONT id=FontSizeSettings5>输入各项磨削加工参数,并通过仿真模型的运算最后得出磨削温度场的三维等温图。通过改变各项参数的大小来比较输出温度场的变化,可以清楚地了解到参数变化对温度场的影响。 </FONT><DD><FONT id=FontSizeSettings5>以干式磨削难加工材料TC4钛合金工件为例,磨削用量如下:砂轮转速n<SUB>s</SUB>= 143r/min,砂轮直径d<SUB>s</SUB>=245.2mm,砂轮速度v<SUB>s</SUB>=18.46m/ s,工件速度v<SUB>w</SUB>=14m/min=23.33cm/s,磨削深度a<SUB>p</SUB>=0.01mm,切入磨削。试件宽度b=10mm、长度l=170mm。平均切向磨削力F<SUB>t</SUB>=66.001N。该工件的温度场示意图如图2所示。 </FONT><DD><FONT id=FontSizeSettings5>同样,在加载了对流换热的边界条件后就可以得到湿式磨削温度场的示意图,经过与实测温度相比较,误差在10%以内。图3为磨削温度场的仿真过程示意图,图4为表面层下不同深度处干磨和湿磨仿真与试验结果的比较。 </FONT><DD><FONT id=FontSizeSettings5>由图4 可以看出,在相同的磨削条件下,干磨与湿磨的最高温度相差很大(一般在200~300℃之间),尤其是工件表层温度梯度较大,干磨时极易引起烧伤,而湿磨时的温度远远低于工件的烧伤温度。运用仿真模型所得分析结果与实测值相差不大,误差在10%以内(仿真分析所得结果比实测值约高10%)。这是由于测温试件本身具有误差,而且实测温度是距工件表层一定深度下的温度值,所以仿真分析所得结果更接近工件实际温度值。 </FONT><H1><FONT id=FontSizeSettings5><FONT size=2>4 结论</FONT></FONT></H1><DD><FONT id=FontSizeSettings5>在磨削温度场的研究中,应用计算机仿真技术可以准确地得出整体磨削温度分布图,从而减少了使用测试仪器所带来的误差,所得等温图更为简便和直观。同时,由于磨削温度场的复杂性,采用仿真模型更易于分析不同加工参数对磨削弧区温度的影响,从而发现磨削温度场的变化规律。应用仿真技术还可以对磨削温度场进行优化:在改变加工参数的条件下,使磨削温度场的温度变化趋向于合理,从而减少磨削烧伤的产生。 </FONT><DD><FONT id=FontSizeSettings5>磨削温度场的计算机仿真技术不仅适用于目前流行的各种高效磨削技术,同时也为磨削全过程的仿真打下了很好的基础。</FONT></DD>
页: [1]
查看完整版本: 计算机仿真技术在磨削温度场中的应用

中国磨床技术论坛
论 坛 声 明 郑重声明:本论坛属技术交流,非盈利性论坛。本论坛言论纯属发表者个人意见,与“中国磨削技术论坛”立场无关。 涉及政治言论一律删除,请所有会员注意.论坛资源由会员从网上收集整理所得,版权属于原作者. 论坛所有资源是进行学习和科研测试之用,请在下载后24小时删除, 本站出于学习和科研的目的进行交流和讨论,如有侵犯原作者的版权, 请来信告知,我们将立即做出整改,并给予相应的答复,谢谢合作!

中国磨削网