CN104615876A - 一种复合材料切削热分配系数的计算方法 - Google Patents

一种复合材料切削热分配系数的计算方法 Download PDF

Info

Publication number
CN104615876A
CN104615876A CN201510042392.2A CN201510042392A CN104615876A CN 104615876 A CN104615876 A CN 104615876A CN 201510042392 A CN201510042392 A CN 201510042392A CN 104615876 A CN104615876 A CN 104615876A
Authority
CN
China
Prior art keywords
heat
centerdot
cutter
formula
cutting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510042392.2A
Other languages
English (en)
Other versions
CN104615876B (zh
Inventor
贾振元
殷俊伟
王福吉
陈晨
朱浩杰
宿友亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201510042392.2A priority Critical patent/CN104615876B/zh
Publication of CN104615876A publication Critical patent/CN104615876A/zh
Application granted granted Critical
Publication of CN104615876B publication Critical patent/CN104615876B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Numerical Control (AREA)
  • Machine Tool Sensing Apparatuses (AREA)

Abstract

本发明一种复合材料切削热分配系数的计算方法属于复合材料切削加工技术领域,涉及一种将碳纤维复合材料纤维方向考虑在内的切削热分配系数的计算方法。该方法根据复合材料不同纤维方向、不同工艺参数及刀具几何参数下所产生的已加工表面回弹量,并结合傅里叶传导定律及其传热反问题的求解方法,得到传入刀具的热量,再通过对加工过程消耗总能量的计算,最终获得加工过程刀具和工件之间的热分配系数。本发明推导的公式中,将纤维方向和刀具角度全部考虑在内,得到刀具前角、刀具后角、切深、纤维方向对切削热分配系数的影响关系;在加工不同纤维方向的复合材料过程中,准确计算出刀具与工件的切削热分配系数,使计算出的温度场分布更接近实际。

Description

一种复合材料切削热分配系数的计算方法
技术领域
本发明属于复合材料切削加工技术领域,涉及一种将碳纤维复合材料纤维方向考虑在内的切削热分配系数的计算方法。
背景技术
复合材料中的树脂基体是一种对温度极其敏感的热固性高分子材料,在高温的影响下极易发生软化、失去强度甚至碳化分解,这就容易使材料中的增强纤维失去支撑和保护,在受到切削力的作用时加剧损伤的产生。而切削热分配系数是决定切削热传入工件热量大小的关键,若能发明一种加工过程中的切削热分配系数计算方法,对研究工件温度分布、减小切削热影响、提高加工效率具有重要的应用价值。由于复合材料具有各向异性、层叠特征,导致复合材料切削热传导过程与金属存在差异,导致金属加工过程中切削热分配系数计算方法不适用于复合材料。
针对这一问题,Olga Klinkova、Sylvain Drapier等人发表的《Characterization of friction properties at theworkmaterial/cutting tool interface during the machining ofrandomly structured carbon fibers reinforced polymer withcarbide tools under dry conditions》一文在《TribologyInternational》2011年第44期第2050到2058页中用材料的热导率、密度、比热容来计算碳纤维复合材料磨擦试验过程中切削热流入磨头的比例系数;Chung-Shin,Je-Ee Ho,Chang-Hu Chan andBaw-Chiang Hwang发表的《Prediction of Cutting Temperature inTurning Carbon Fiber Reinforced Plastics Composites with WornTools》一文在《Journal of Applied Sciences》2011年第22期第3698到3671页中利用传热反问题的计算方法对碳纤维复合材料车削过程中切削热传入切屑的比例进行研究。但是,上述研究中不仅计算量大、计算复杂,而且没有考虑纤维方向对热分配系数的影响,导致计算结果误差较大。
发明内容
本发明针对现有方法无法考虑纤维方向对切削热分配系数影响的缺陷,发明一种复合材料切削热分配系数的计算方法。该方法根据复合材料不同纤维方向、不同工艺参数及刀具几何参数下所产生的已加工表面回弹量,并结合傅里叶传导定律及其传热反问题的求解方法,得到传入刀具的热量。由于复合材料加工过程塑性变形很小,主要的热源来自刀具与复合材料之间的摩擦,通过对加工过程消耗总能量的计算,最终获得加工过程刀具和工件之间的热分配系数。使用本发明推导的公式将纤维方向和刀具角度全部考虑在内,计算方法简单准确。
本发明采用的技术方案是一种复合材料切削热分配系数的计算方法,该方法根据复合材料不同纤维方向、不同工艺参数及刀具几何参数下所产生的已加工表面回弹量,并结合傅里叶传导定律及其传热反问题的求解方法,得到传入刀具的热量,再通过对加工过程消耗总能量的计算,最终获得加工过程刀具和工件之间的热分配系数;计算方法采用如下步骤:
步骤1:利用傅里叶热传导定律推导出切削热分配给刀具部分的比例系数表达式为,
λ 1 = k · A · ΔT L · W - - - ( 1 )
式中,λ1为切削热分配给刀具部分比例系数,W为切削过程产生的总热量,A为刀具与工件接触面积,k为刀具的热导率,ΔT为刀具测温点与刀刃之间的温度差,L为刀具上测温点与刀刃之间的距离。
切削热分配给工件部分的比例系数表达式为,
λ 2 = 1 - k · A · ΔT L · W - - - ( 2 )
式中,λ2为切削热分配给工件部分的比例系数。
步骤2:根据测力仪测量所得的平均水平切削力与设定的切削速度的乘积计算加工过程中单位时间内产生的总热量,表达式为,
W=Ft·V   (3)
式中,Ft表示平均水平切削力,V表示切削速度。
步骤3:根据不同纤维方向复合材料的已加工表面回弹量不同计算刀具与工件的接触面积。考虑复合材料纤维方向、刀尖圆弧与工件接触的长度、回弹的已加工表面与刀具后刀面的接触长度、刀具前角和刀具后角得到刀具与工件接触面积的表达式,
A = [ 90 + α + γ 180 π R e + 1 sin α · ( 3 E y 4 E * ( θ ) R e 1 2 ) 2 3 ] · l t - - - ( 4 )
式中,θ为纤维方向,定义为从切削方向顺时针旋转到纤维所处位置而经过的角度,E*(θ)为刀具和复合材料的等效法向弹性模量,Fy平均法向接触力,Re为刀尖圆弧半径,lt为工件厚度,α为刀具后角,γ为刀具前角。
步骤4:根据传热反问题求解方法得到刀刃处与测温点处温度差的关系表达式,
ΔT=1.3Ttest-Ttest=0.3Ttest   (5)
式中,Ttest为测温点处的温度值。
步骤5:将式(3)、式(4)、式(5)带入公式(1)中,得到切削热分配给刀具的比例系数关系式为,
λ 1 = 0.3 k T test l t L · F t · V · [ 90 + α + γ 180 π R e + 1 sin α · ( 3 E y 4 E * ( θ ) R e 1 2 ) 2 3 ] - - - ( 6 )
切削热分配给工件的比例系数关系式为,
λ 2 = 1 - 0.3 k T test l t L · F t · V · [ 90 + α + γ 180 π R e + 1 sin α · ( 3 E y 4 E * ( θ ) R e 1 2 ) 2 3 ] - - - ( 7 )
完成复合材料切削热分配系数的计算。
本发明有益效果是:由于复合材料加工过程中切削热分配系数与纤维方向和刀具角度有关,使用本发明推导的公式将纤维方向和刀具角度全部考虑在内,并且得到刀具前角、刀具后角、切深、纤维方向对切削热分配系数的影响关系。不仅计算简单,而且计算结果准确,可以使计算出的温度分布更接近实际情况。
附图说明
图1为本发明实施例的刀具与工件接触示意图。图中γ表示刀具前角,α表示刀具后角,θ表示纤维方向,Re为刀尖圆弧半径,lt为工件厚度,Δh为已加工表面回弹量。
图2为刀具与工件接触面积随纤维方向变化趋势。图中横坐标表示纤维方向,单位为度;纵坐标表示刀具与工件的接触面积,单位为mm2
图3为刀具的切削热分配系数随纤维方向变化趋势。图中横坐标表示纤维方向,单位为度;纵坐标表示刀具的切削热分配系数。
图4为工件的切削热分配系数随纤维方向变化趋势。图中横坐标表示纤维方向,单位为度;纵坐标表示工件的切削热分配系数。
具体实施方式
下面结合附图和技术方案详细说明本发明的具体实施,
本实施例中纤维方向θ定义为从切削方向顺时针旋转到纤维所处位置而经过的角度。本实施例具体针对纤维方向θ=0°、45°、90°、135°四种典型纤维方向,在刀具前角γ=25°,刀具后角α=17°,切削速度为500mm/min、切削深度为150μm的情况下,计算加工过程中切削热分配系数。
切削热分配系数的计算方法包含以下步骤:
步骤1:建立切削热分配系数与各影响因素之间的关系表达式。利用傅里叶热传导定律推导出加工过程中切削热分配给刀具的比例系数表达式,见上述公式(1),切削热分配给工件的比例系数表达式,见上述公式(2)。
步骤2:计算单位时间内加工过程中产生的总热量W。通过测量得到的加工过程平均水平切削力与切削速度的乘积计算加工过程消耗的总能量。由于加工过程消耗的总能量除部分转化为切屑动能、断裂能以外,其余都转化为热量,可近似地认为切削时所消耗的能量全部转化为热量。即加工过程消耗的总能量全部转化为热量。
本实施例中,利用测力仪测得从切削方向顺时针旋转到纤维所处位置而经过的纤维方向θ角度分别为0°、45°、90°、135°四种典型纤维方向对应的平均水平切削力Ft大小依次为46.905N、108.188N、251.448N、184.048N。将这些值分别带入公式(3)可以得到本实施例中纤维方向θ=0°、45°、90°、135°四种典型纤维方向,在加工过程中单位时间内产生的总热量W分别为0.391J、0.902J、2.095J、1.534J。
步骤3:计算刀具与工件的接触面积A。由于刀刃存在一定的钝圆,工件已加工表面会产生回弹,并且切削不同纤维方向复合材料工件时对应的已加工表面回弹量不同,导致刀具后刀面与工件已加工表面的接触长度不同,进而导致刀具与工件接触面积不同,如附图1所示。
将刀具前角γ=25°,刀具后角α=17°,刀尖圆弧半径54μm,工件厚度3mm,纤维方向θ=0°、45°、90°、135°四种典型纤维方向下测量得到的法相接触力0.817N,14.829N,11.21N,20.24,四种纤维方向对应的刀具和复合材料的等效法向弹性模量值9.59GPa,15.4GPa,136GPa,15.4GPa,一起带入公式(4),分别得到本实施例中纤维方向θ=0°、45°、90°、135°四种典型纤维方向对应的刀具与工件接触面积分别为0.421mm2、0.596mm2、0.42mm2、0.646mm2。刀具与工件接触面积随纤维方向变化趋势见附图2,从附图2可以看出,不同纤维方向下刀具与工件的接触面积明显不同。
步骤4:计算刀具测温点与刀刃之间的温度差ΔT。本实施例分别测量切削纤维方向θ=0°、45°、90°、135°四种典型纤维方向时,与刀具切削刃距离L为2.5mm处的温度值Ttest,测得四种纤维方向下对应的Ttest值分别为16.7℃,17.7℃,20.5℃,19.5℃。将测得的Ttest值分别代入公式(5),得到本实施例中切削纤维方向θ=0°、45°、90°、135°四种典型纤维方向工件时,刀具测温点与刀刃之间的温度差分别为5.01℃,5.31℃,6.15℃,5.85℃。
步骤5:将上述各步骤计算得到的数据全部带入公式(6)中,得到本实施例中纤维方向θ=0°、45°、90°、135°四种典型纤维方向对应的刀具的切削热分配系数分别为0.0908、0.059、0.0405、0.0808;工件的切削热分配系数分别为0.9092、0.941、0.9595、0.9192。刀具的切削热分配系数随纤维方向变化趋势见附图3,从附图3中可以看出,刀具的切削热分配系数随纤维方向呈现先减小后增大的变化趋势,且在切削90°纤维方向时刀具的切削热分配系数最小。工件的切削热分配系数随纤维方向变化趋势见附图4,从附图4可以看出,工件的切削热分配系数随纤维方向呈现先增大后减小的变化趋势,且在切削90°纤维方向时工件的切削热分配系数最大。
通过本实施例的计算可以得出:在复合材料切削过程中,工件的切削热分配系数远大于刀具的切削热分配系数,即加工过程中大量的切削热会传入工件,造成树脂基体软化,加剧损伤产生。此外,切削热分配系数随纤维方向变化十分明显。应用本发明的计算方法不仅计算过程简单,还可以准确计算出加工不同纤维方向复合材料过程中刀具与工件的切削热分配系数,使计算出的温度场分布更接近实际。

Claims (1)

1.一种复合材料切削热分配系数的计算方法,其特征是,该方法根据复合材料不同纤维方向、不同工艺参数及刀具几何参数下所产生的已加工表面回弹量,结合傅里叶传导定律及其传热反问题的求解方法,得到传入刀具的热量,再通过对加工过程消耗总能量的计算,最终获得加工过程刀具和工件之间的热分配系数;计算方法采用如下步骤:
步骤1:利用傅里叶热传导定律推导出切削热分配给刀具部分的比例系数表达式为,
λ 1 = k · A · ΔT L · W - - - ( 1 )
式中,λ1为切削热分配给刀具部分比例系数,W为切削过程产生的总热量,A为刀具与工件接触面积,k为刀具的热导率,ΔT为刀具测温点与刀刃之间的温度差,L为刀具上测温点与刀刃之间的距离;
切削热分配给工件部分的比例系数表达式为,
λ 2 = 1 - k · A · ΔT L · W - - - ( 2 )
式中,λ2为切削热分配给工件部分的比例系数;
步骤2:根据测力仪测量所得的平均水平切削力与设定的切削速度的乘积计算加工过程中单位时间内产生的总热量,表达式为,
W=Ft·V    (3)
式中,Ft表示平均水平切削力,V表示切削速度;
步骤3:根据不同纤维方向复合材料的已加工表面回弹量不同计算刀具与工件的接触面积考虑复合材料纤维方向、刀尖圆弧与工件接触的长度、回弹的已加工表面与刀具后刀面的接触长度、刀具前角和刀具后角得到刀具与工件接触面积的表达式:
A = [ 90 + α + γ 180 π R e + 1 sin α · ( 3 F y 4 E * ( θ ) R e 1 2 ) 2 3 ] · l t - - - ( 4 )
式中,θ为纤维方向,定义为从切削方向顺时针旋转到纤维所处位置而经过的角度,E*(θ)为刀具和复合材料的等效法向弹性模量,Fy平均法向接触力,Re为刀尖圆弧半径,lt为工件厚度,α为刀具后角,γ为刀具前角:
步骤4:根据传热反问题求解方法得到刀刃处与测温点处温度差的关系表达式为:
ΔT=1.3Ttest-Ttest=0.3Ttest    (5)
式中,Ttest为测温点处的温度值:
步骤5:将式(3)、式(4)、式(5)带入公式(1)中,得到切削热分配给刀具的比例系数关系式为:
λ 1 = 0.3 k T test l t L · F t · V · [ 90 + α + γ 180 π R e + 1 sin α · ( 3 F y 4 E * ( θ ) R e 1 2 ) 2 3 ] - - - ( 6 )
切削热分配给工件的比例系数关系式为:
λ 2 = 1 - 0.3 k T test l t L · F t · V · [ 90 + α + γ 180 π R e + 1 sin α · ( 3 F y 4 E * ( θ ) R e 1 2 ) 2 3 ] - - - ( 7 )
完成复合材料切削热分配系数的计算。
CN201510042392.2A 2015-01-28 2015-01-28 一种复合材料切削热分配系数的计算方法 Expired - Fee Related CN104615876B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510042392.2A CN104615876B (zh) 2015-01-28 2015-01-28 一种复合材料切削热分配系数的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510042392.2A CN104615876B (zh) 2015-01-28 2015-01-28 一种复合材料切削热分配系数的计算方法

Publications (2)

Publication Number Publication Date
CN104615876A true CN104615876A (zh) 2015-05-13
CN104615876B CN104615876B (zh) 2017-06-23

Family

ID=53150317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510042392.2A Expired - Fee Related CN104615876B (zh) 2015-01-28 2015-01-28 一种复合材料切削热分配系数的计算方法

Country Status (1)

Country Link
CN (1) CN104615876B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104834786A (zh) * 2015-05-15 2015-08-12 大连理工大学 碳纤维复合材料去除过程的细观仿真建模方法
CN105138842A (zh) * 2015-08-26 2015-12-09 大连理工大学 一种碳纤维复合材料钻削损伤的表征方法
CN105352631A (zh) * 2015-10-27 2016-02-24 沈阳航空航天大学 一种加工过程中摩擦热占整体切削热比例的测定方法
CN109408926A (zh) * 2018-10-12 2019-03-01 大连理工大学 求解复杂结构多维瞬态非线性热传导反问题的方法
CN110907494A (zh) * 2019-12-12 2020-03-24 河南科技大学 一种用于检测摩擦副热分配系数的检测系统及检测方法
CN111730408A (zh) * 2020-07-17 2020-10-02 重庆工商大学 一种考虑切屑换热的干切削加工强化传热方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103324781A (zh) * 2013-05-27 2013-09-25 华中科技大学 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法
CN103639895A (zh) * 2013-11-15 2014-03-19 湖南大学 一种基于传热反问题的磨削温度测量方法
CN104268343A (zh) * 2014-09-28 2015-01-07 北京理工大学 一种用于端铣切削的切削力预测及温度预测的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103324781A (zh) * 2013-05-27 2013-09-25 华中科技大学 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法
CN103639895A (zh) * 2013-11-15 2014-03-19 湖南大学 一种基于传热反问题的磨削温度测量方法
CN104268343A (zh) * 2014-09-28 2015-01-07 北京理工大学 一种用于端铣切削的切削力预测及温度预测的方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CUI Y X.ET: "Dynamic calibration of the cutting temperature sensor of NiCr/NiSi thin-film thermocouple", 《CHINESE JOURNAL OF MECHANICAL ENGINEERING》 *
刘旺玉等: "高速切削热分配的有限元模拟", 《工具技术》 *
崔云先等: "测温刀具温度场分布及其动态性能仿真预测", 《大连交通大学学报》 *
杜美憬: "高速切削刀具内的传热分析", 《中国优秀硕士学位论文全文数据库工程科技I辑》 *
殷俊伟: "NiCr_NiSi薄膜热电偶测温刀具及其传热数学模型研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104834786A (zh) * 2015-05-15 2015-08-12 大连理工大学 碳纤维复合材料去除过程的细观仿真建模方法
CN104834786B (zh) * 2015-05-15 2018-07-13 大连理工大学 碳纤维复合材料去除过程的细观仿真建模方法
CN105138842A (zh) * 2015-08-26 2015-12-09 大连理工大学 一种碳纤维复合材料钻削损伤的表征方法
CN105138842B (zh) * 2015-08-26 2017-12-05 大连理工大学 一种碳纤维复合材料钻削损伤的表征方法
CN105352631A (zh) * 2015-10-27 2016-02-24 沈阳航空航天大学 一种加工过程中摩擦热占整体切削热比例的测定方法
CN105352631B (zh) * 2015-10-27 2017-11-10 沈阳航空航天大学 一种加工过程中摩擦热占整体切削热比例的测定方法
CN109408926A (zh) * 2018-10-12 2019-03-01 大连理工大学 求解复杂结构多维瞬态非线性热传导反问题的方法
CN109408926B (zh) * 2018-10-12 2023-04-07 大连理工大学 求解复杂结构多维瞬态非线性热传导反问题的方法
CN110907494A (zh) * 2019-12-12 2020-03-24 河南科技大学 一种用于检测摩擦副热分配系数的检测系统及检测方法
CN110907494B (zh) * 2019-12-12 2022-02-15 河南科技大学 一种用于检测摩擦副热分配系数的检测系统及检测方法
CN111730408A (zh) * 2020-07-17 2020-10-02 重庆工商大学 一种考虑切屑换热的干切削加工强化传热方法
CN111730408B (zh) * 2020-07-17 2021-10-15 重庆工商大学 一种考虑切屑换热的干切削加工强化传热方法

Also Published As

Publication number Publication date
CN104615876B (zh) 2017-06-23

Similar Documents

Publication Publication Date Title
CN104615876A (zh) 一种复合材料切削热分配系数的计算方法
Xu et al. Research on chip formation parameters of aluminum alloy 6061-T6 based on high-speed orthogonal cutting model
CN101412196A (zh) 基于刀具角度与切削用量参数变化的车削力预测方法
Li et al. Quadratic curve heat flux distribution model in the grinding zone
CN105643024A (zh) 一种车削大螺距螺纹轴向分层切削方法、刀具磨损测试方法及其力热载荷计算方法
Guo et al. Machinability of wood fiber/polyethylene composite during orthogonal cutting
CN101853324A (zh) 圆周铣过程中铣削力建模方法
CN108115357A (zh) 一种大型水轮机叶片的在线加工方法
CN103105332B (zh) 用正交切削法测定材料断裂韧性的方法
Horman et al. Temperature influence on wear characteristics and blunting of the tool in continuous wood cutting process
Wu et al. 3D FEM simulation of chip breakage in turning AISI1045 with complicate-grooved insert
CN105843171A (zh) 基于计算机软件的单点金刚石车削加工中表面粗糙度的预测方法
CN105787179A (zh) 用于碳纤维增强单向层合板二维切削温度的建模方法
Wei et al. Research on milling forces during high-speed milling of wood-plastic composites
CN104128646B (zh) 工件材料最小切削厚度的确定方法
Wei et al. Study on surface roughness model of 3D ultrasonic vibration–assisted turning driven by a single actuator
CN105181508A (zh) 难加工材料去除量和刀具磨损之间的匹配模型
Li et al. Modeling of instantaneous cutting force for large pitch screw with vibration consideration of the machine tool
Chang et al. Experimental and FEM study of residual stresses during ultra-precision turning of aluminum 2024-T3
Xu et al. Experimental research on the critical conditions and critical equation of chip splitting when turning a C45E4 disc workpiece symmetrically with a high-speed steel double-edged turning tool
Yang et al. Investigation on the temperature field under the action of the blunt tool edge for precision cutting of titanium alloys
Han et al. Experiment study on deflection of aluminum alloy thin-wall workpiece in milling process
Shi et al. Evaluation of a regression prediction model for surface roughness of wood–polyethylene composite (WPC)
Xu et al. Effect of Softening Treatment on Cutting Force during Slicing the Veneers of Common Fast-growing Wood.
Xiang et al. Study on milling force of high volume fraction SiCp/Al composites with ultrasonic longitudinal and torsional vibration high speed milling

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170623

CF01 Termination of patent right due to non-payment of annual fee