发明名称 一种确定立铣刀铣削过程切削热分配比例的方法
摘要 本发明公开了一种确定立铣刀铣削过程切削热分配比例的方法,步骤是:计算刀屑接触区的表面热流密度;利用ABAQUS软件对铣刀三维模型进行有限元仿真分析得到铣刀的温度场结果;进行刀具测温实验,利用OMEGA软件得到实验中刀杆上部多个测温点的温度:根据实验中多个测温点的位置将仿真结果中对应点的温度提取出来与刀杆上多个测温点的温度进行对比,根据对比结果调整热分配系数A%,当仿真结果与测温实验结果相差±15%时,即可得出此工况下的具体的热分配系数值A%。本发明方法具有实验时间较短,计算准确的优点,工程实际应用效果好。可以方便的使用有限元分析软件进行模拟仿真,并通过实验数据与仿真结果进行比较、调整,从而得到较为精确的热分配比例结果。
申请公布号 CN106484961A 申请公布日期 2017.03.08
申请号 CN201610837004.4 申请日期 2016.09.20
申请人 天津大学 发明人 任成祖;杨新鹏;陈光;邹云鹤;靳新民
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 李丽萍
主权项 一种确定立铣刀铣削过程切削热分配比例的方法,其特征在于,包括以下步骤:步骤一、计算刀屑接触区的表面热流密度g,包括:首先,根据式(1)求出铣刀切削热功率P,<maths num="0001"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mfrac><mrow><mi>V</mi><mo>&times;</mo><msub><mi>dF</mi><mi>R</mi></msub><mo>&times;</mo><msub><mi>sin&phi;</mi><mi>n</mi></msub><mo>&times;</mo><mi>s</mi><mi>i</mi><mi>n</mi><mi>&lambda;</mi></mrow><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>&phi;</mi><mi>n</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001117173850000011.GIF" wi="1510" he="88" /></maths>式(1)中:V为切削速度,α<sub>n</sub>为法向前角,α<sub>n</sub>=tan<sup>‑1</sup>(tanα<sub>r</sub>cosα<sub>h</sub>),α<sub>r</sub>为径向前角,α<sub>h</sub>为螺旋升角;φ<sub>n</sub>为法向剪切角,<img file="FDA0001117173850000012.GIF" wi="378" he="72" />λ为刀屑摩擦角,λ=19.1+0.29α<sub>n</sub>(°);<maths num="0002"><math><![CDATA[<mrow><msub><mi>dF</mi><mi>R</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>dF</mi><mi>s</mi></msub><msub><mi>cos&eta;</mi><mi>s</mi></msub></mrow><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mi>&lambda;</mi><mrow><mo>(</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo>(</mo><mrow><msub><mi>&phi;</mi><mi>n</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>n</mi></msub></mrow><mo>)</mo><mo>-</mo><msub><mi>tan&lambda;cos&eta;</mi><mi>c</mi></msub><mi>s</mi><mi>i</mi><mi>n</mi><mo>(</mo><mrow><msub><mi>&phi;</mi><mi>n</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>n</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001117173850000013.GIF" wi="1519" he="84" /></maths>式(2)中,dF<sub>s</sub>为微元的剪切力,η<sub>c</sub>为切屑流出速度,令η<sub>c</sub>=α<sub>h</sub>,η<sub>s</sub>为剪切带流动方向角,<maths num="0003"><math><![CDATA[<mrow><msub><mi>&eta;</mi><mi>s</mi></msub><mo>=</mo><msup><mi>tan</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>tan&eta;</mi><mi>c</mi></msub><msub><mi>sin&phi;</mi><mi>n</mi></msub><mo>-</mo><msub><mi>tan&lambda;</mi><mi>s</mi></msub><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><msub><mi>&phi;</mi><mi>n</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>cos&alpha;</mi><mi>n</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001117173850000014.GIF" wi="1518" he="87" /></maths>式(3)中,α<sub>s</sub>为刃倾角等于α<sub>h</sub>,<maths num="0004"><math><![CDATA[<mrow><msub><mi>dF</mi><mi>s</mi></msub><mo>=</mo><mfrac><mrow><mi>&tau;</mi><mi>h</mi></mrow><mrow><msub><mi>cos&lambda;</mi><mi>S</mi></msub><msub><mi>sin&phi;</mi><mi>n</mi></msub></mrow></mfrac><mi>d</mi><mi>z</mi><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001117173850000015.GIF" wi="1517" he="90" /></maths>式(4)中,h为每齿进给量,dz为微元长度,τ为平均剪切应力,<maths num="0005"><math><![CDATA[<mrow><mi>&tau;</mi><mo>=</mo><mfrac><mn>1</mn><msqrt><mn>3</mn></msqrt></mfrac><mo>&lsqb;</mo><mi>A</mi><mo>+</mo><mi>B</mi><msup><mrow><mo>(</mo><mfrac><mi>&gamma;</mi><msqrt><mn>3</mn></msqrt></mfrac><mo>)</mo></mrow><mi>n</mi></msup><mo>&rsqb;</mo><mo>&lsqb;</mo><mn>1</mn><mo>+</mo><mi>C</mi><mi>ln</mi><mrow><mo>(</mo><mfrac><mover><mi>&gamma;</mi><mo>&CenterDot;</mo></mover><msub><mover><mi>&gamma;</mi><mo>&CenterDot;</mo></mover><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mfrac><mrow><mi>T</mi><mo>-</mo><msub><mi>T</mi><mn>0</mn></msub></mrow><mrow><msub><mi>T</mi><mi>m</mi></msub><mo>-</mo><msub><mi>T</mi><mn>0</mn></msub></mrow></mfrac><mo>)</mo></mrow><mi>m</mi></msup><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001117173850000016.GIF" wi="1518" he="95" /></maths>式(5)中,A为工件材料屈服强度,B为工件材料硬化模量,γ为应变,<img file="FDA0001117173850000017.GIF" wi="403" he="81" />n为加工硬化指数,C为应变率敏感性因数,<img file="FDA0001117173850000018.GIF" wi="29" he="47" />为应变率,<img file="FDA0001117173850000019.GIF" wi="449" he="82" />其中,Δy为剪切带厚度,剪切速度<img file="FDA00011171738500000110.GIF" wi="290" he="79" /><img file="FDA00011171738500000111.GIF" wi="43" he="47" />为参考应变率,T为剪切带温度,T<sub>0</sub>为参考温度,T<sub>m</sub>为材料溶化温度,m为热软化系数;然后,由上述得出的热功率P求得面热流密度<img file="FDA00011171738500000112.GIF" wi="265" he="101" />其中,l<sub>contact</sub>为铣刀与切屑的接触长度,<img file="FDA00011171738500000113.GIF" wi="464" he="83" />步骤二:利用ABAQUS软件对铣刀三维模型进行热传递有限元分析,包括:2‑1)根据所用铣刀建立铣刀三维模型及网格划分;2‑2)设定铣刀三维模型材料属性,包括铣刀材料的密度、比热容、热传导系数、杨氏模量、泊松比;2‑3)载荷及边界条件设定:以步骤一的面热流密度乘以热分配系数即g×A%作为载荷施加到铣刀三维模型刀刃上的刀屑接触区域,其中,A%的初设值为50%;对整个铣刀施加边界条件,边界条件为空气自然对流;2‑4)热传递有限元仿真分析后得到铣刀的温度场;步骤三:利用OMEGA软件得到刀具测温实验中刀杆上部多个测温点的温度:对刀具进行测温实验:在铣刀柄(1)上安装多个温度采集模块(4),在铣刀(7)的上部设有多个测温点(5),所述温度采集模块(4)具有热电偶插头(3),所述热电偶插头(3)与测温点(5)之间均分别连接有热电偶(11),所述温度采集模块(4)与所述计算机之间设有匹配的插接件;利用上述铣刀(7)对工件(14)进行铣削加工,铣削过程结束后,所述温度采集模块(4)将采集到的铣削过程铣刀的温度传输给计算机;所述计算机利用OMEGA软件得到刀杆上所述多个测温点的温度;步骤四:确定热分配系数A%的终值:根据步骤三中刀杆上多个测温点的位置将步骤二仿真结果中对应点的温度提取出来与步骤三得到的刀杆上所述多个测温点的温度进行对比,如果仿真结果大于测温实验结果则按照3%减小热分配系数A%,反之增大热分配系数A%,重复依次执行步骤2‑3)、2‑4)和步骤四,当满足仿真结果与测温实验结果相差±15%时,从而得出热分配系数A%。
地址 300072 天津市南开区卫津路92号