CN112329245A - 一种基于最大剪应力原理的切削力建模系数标定方法 - Google Patents

一种基于最大剪应力原理的切削力建模系数标定方法 Download PDF

Info

Publication number
CN112329245A
CN112329245A CN202011241933.1A CN202011241933A CN112329245A CN 112329245 A CN112329245 A CN 112329245A CN 202011241933 A CN202011241933 A CN 202011241933A CN 112329245 A CN112329245 A CN 112329245A
Authority
CN
China
Prior art keywords
force
cutting
angle
normal
coefficient
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.)
Pending
Application number
CN202011241933.1A
Other languages
English (en)
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.)
Beijing Research Institute of Precise Mechatronic Controls
Original Assignee
Beijing Research Institute of Precise Mechatronic Controls
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 Beijing Research Institute of Precise Mechatronic Controls filed Critical Beijing Research Institute of Precise Mechatronic Controls
Priority to CN202011241933.1A priority Critical patent/CN112329245A/zh
Publication of CN112329245A publication Critical patent/CN112329245A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Turning (AREA)

Abstract

一种基于最大剪应力原理的切削力建模系数标定方法,包括:建立坐标系,将切削力分解为切向力、轴向力和径向力,在法平面坐标系下采用最大剪应力原理,提出一种切削力系数标定方法,本发明能够解析切削力模型系数和切削几何参数之间的约束关系,揭示切削过程中切削力的变化规律。

Description

一种基于最大剪应力原理的切削力建模系数标定方法
技术领域
本发明涉及一种基于最大剪应力原理的切削力建模系数标定方法,属于机械制造技术领域。
背景技术
切削过程是多物理场耦合下的材料变形过程,因此,切削力建模过程中往往根据需求进行一定的简化,以提高计算效率。Merchant最初提出的剪切模型中,不考虑金属侧向流动,且认为仅在剪切平面上产生切削变形。Lee假设材料是塑性固体,采用滑移场线理论建立正交切削计算模型。Shaw重新考虑塑性带上的应力分布情况,提出一种改进后的滑移场线计算模型。Moufki提出一种考虑热力耦合效应的切削力计算模型。然而,上述切削力计算模型均不适用于斜角切削过程。
因而,研究材料斜角切削过程中的切削力变化规律,提出切削力模型系数标定方法,为生产加工提供恰当的切削参数,对促进制造技术的发展是非常重要的。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出了一种基于最大剪应力原理的切削力建模系数标定方法,建立包含关键几何参数的力学计算模型,构型简洁,提高了力学模型的适应性。
本发明的技术方案是:
一种基于最大剪应力原理的切削力建模系数标定方法,包括步骤如下:
1)使用刀具切削圆柱体金属件,采集获得多组切削力试验数据;
2)根据任意一组切削力试验数据,将切削力在前刀面上分解为切向力FT、轴向力FA和径向力FR;然后,根据刀具刃倾角λs,将切向力FT、轴向力FA和径向力FR转换为法向切向力FTn、法向轴向力FRn和法向径向力FAn
3)根据切削力试验数据获得切削合力F,从而确定切削合力F在法平面上的投影与切削平面之间夹角ξn、切削合力F与法平面之间的夹角ξλ、法向摩擦角βn;根据切削合力F,确定剪力Fs,从而确定剪力Fs与法平面之间的夹角
Figure BDA0002768694450000021
剪力Fs在法平面上的投影与切削平面的夹角
Figure BDA0002768694450000022
4)根据刀具刃倾角λs,根据步骤3)所述剪力Fs和剪力Fs在法平面上的投影与切削平面的夹角
Figure BDA0002768694450000023
确定剪应力τs
5)根据步骤3)所述法向摩擦角βn、剪力Fs在法平面上的投影与切削平面的夹角
Figure BDA0002768694450000024
根据步骤4)所述剪应力τs,确定切向力系数KR、径向力系数KA和轴向力系数KT
步骤2)所述将切向力FT、轴向力FA和径向力FR转换为法向切向力FTn、法向轴向力FRn和法向径向力FAn的方法,具体为:
(FTn,FRn,FAn)T=T(λs)(FT,FR,FA)T
Figure BDA0002768694450000025
确定步骤3)所述ξn的方法,具体为:
Figure BDA0002768694450000026
确定步骤3)所述ξλ的方法,具体为:
sinξλ=sinβasinη
tan(ξnn)=tanβacosη
其中,η为流屑角。
确定步骤3)所述法向摩擦角βn的方法,具体为:
βn=ξnn
其中,γn为法向前角。
确定步骤3)所述剪力Fs在法平面投影与切削平面之间夹角
Figure BDA0002768694450000027
的方法,具体为:
Figure BDA0002768694450000031
Figure BDA0002768694450000032
确定步骤3)所述剪力Fs的方法,具体为:
Figure BDA0002768694450000033
确定步骤4)所述剪应力τs的方法,具体为:
Figure BDA0002768694450000034
Figure BDA0002768694450000035
其中,f为进给量,ap为径向切削深度,As为剪切面积。
步骤5)所述确定切向力系数KR、径向力系数KA和轴向力系数KT的方法,具体为:
Figure BDA0002768694450000036
其中,η为流屑角,γn为刀具的前角。
本发明与现有技术相比的有益效果是:
1)本发明针对圆柱形工件切削加工过程,采用机械力学式的建模方法,建立包含关键几何参数的力学计算模型,构型简洁,提高了力学模型的适应性。
2)本发明建立了基于最大剪应力原理的斜角切削模型,解析切削力模型系数之间的约束关系,揭示了车削力演变过程规律。
3)本发明提出了一种切削力系数标定方法,仅需数次试验即可完成车削力模型系数标定过程,降低了试验成本,提高了建模效率。
附图说明
图1为本发明方法流程图。
具体实施方式
丝杠切削加工过程中,由于工件轮廓曲线曲率变化,切削力容易出现大幅度变化,对刀具和工件造成切削载荷冲击,加剧刀片磨损、破损,甚至引起刀具折断,导致工件加工质量不高,合格率偏低。建立切削力模型,研究切削力变化规律,对实现高效优质生产具有重要意义。
本发明一种基于最大剪应力原理的切削力建模系数标定方法,包括步骤如下:
1)使用刀具以不同的切削参数切削圆柱体金属件,采集获得多组切削力试验数据;大于或得等于20组;
2)根据任意一组切削力试验数据,将切削力在前刀面上分解为切向力FT、轴向力FA和径向力FR;然后,根据刀具刃倾角λs,将切向力FT、轴向力FA和径向力FR转换为法向切向力FTn、法向轴向力FRn和法向径向力FAn,具体为:
(FTn,FRn,FAn)T=T(λs)(FT,FR,FA)T
Figure BDA0002768694450000041
3)根据切削力试验数据获得切削合力F,从而确定切削合力F在法平面上的投影与切削平面之间夹角ξn、切削合力F与法平面之间的夹角ξλ、法向摩擦角βn;根据切削合力F,确定剪力Fs,从而确定剪力Fs与法平面之间的夹角
Figure BDA0002768694450000043
剪力Fs在法平面上的投影与切削平面的夹角
Figure BDA0002768694450000044
Figure BDA0002768694450000042
确定步骤3)所述ξn的方法,具体为:
Figure BDA0002768694450000051
确定步骤3)所述ξλ的方法,具体为:
sinξλ=sinβasinη
tan(ξnn)=tanβacosη
其中,βa为摩擦角,η为流屑角(已知,由切削力试验数据中获得)。
确定步骤3)所述法向摩擦角βn的方法,具体为:
βn=ξnn
其中,γn为法向前角(已知,由切削力试验数据中获得)。
确定步骤3)所述剪力Fs在法平面投影与切削平面之间夹角
Figure BDA0002768694450000059
的方法,具体为:
Figure BDA0002768694450000052
Figure BDA0002768694450000053
其中,
Figure BDA0002768694450000054
为剪力Fs与法平面之间夹角。
4)根据刀具刃倾角λs,根据步骤3)所述剪力Fs和剪力Fs在法平面上的投影与切削平面的夹角
Figure BDA0002768694450000055
确定剪应力τs,具体为:
Figure BDA0002768694450000056
Figure BDA0002768694450000057
其中,f为进给量,ap为径向切削深度,As为剪切面积。
5)根据步骤3)所述法向摩擦角βn、剪力Fs在法平面上的投影与切削平面的夹角
Figure BDA0002768694450000058
根据步骤4)所述剪应力τs,确定切向力系数KR、径向力系数KA和轴向力系数KT,具体为:
Figure BDA0002768694450000061
其中,η为流屑角,γn为刀具的前角。
下面结合附图和具体实施方式对本发明做进一步详细的描述。
在圆柱形工件端面建立坐标系OXYZ,其中O点位于端面圆心,平面XOY与端面重合。切削力在前刀面上的分力可表示为切向力FT、轴向力FA和径向力FR,经矩阵变换,切削力在坐标系OXYZ中可表示为:
Figure BDA0002768694450000062
车削过程中,切向力FT、轴向力FA和径向力FR大小与每次进刀前后的切削层面积SA相关,可表示为:
Figure BDA0002768694450000063
式中,γ0为刀具前角,λs为刃倾角,KR为径向力系数,KA为轴向力系数,KT为切向力系数。
斜角切削方式下,KR、KA和KT可表示为:
Figure BDA0002768694450000071
其中,βn为摩擦角,η为流屑角,
Figure BDA00027686944500000710
为剪切角,τs为剪应力,γn为前角;
根据最大剪应力原理,作用力和剪切面之间的夹角应为π/4,故有:
Figure BDA0002768694450000072
Figure BDA0002768694450000073
式中,Fs为剪力,F为合力,ξλ为F与法平面之间的夹角,
Figure BDA0002768694450000074
为Fs与法平面之间夹角,ξn为F在法平面上的投影与切削平面之间夹角,
Figure BDA0002768694450000075
为Fs在法平面投影与切削平面之间夹角。
从式(4)和(5)可得:
Figure BDA0002768694450000076
Figure BDA0002768694450000077
剪应力τs可表示为:
Figure BDA0002768694450000078
式中,As为剪切面积,可表示为:
Figure BDA0002768694450000079
式中,f为进给量,ap为切削深度,λs为刃倾角。
FTn、FRn和FAn分别为法向切向力、法向径向力和法向轴向力,可表示为:
(FTn,FRn,FAn)T=T(λs)(FT,FR,FA)T (10)
式中,T(λs)为变换矩阵,可表示为
Figure BDA0002768694450000081
另外
Figure BDA0002768694450000082
切削加工过程中,作用在刀具前刀面上的载荷有正向压力Fn 和摩擦力Ff,可以得到:
Ffsinη=Fsinξλ (13)
Fcosη=Fntan(ξnn) (14)
进一步有:
sinξλ=sinβasinη (15)
tan(ξnn)=tanβacosη (16)
其中,η为流屑角;βn为摩擦角
ξnn=βn (17)
式中,βa为摩擦角,βn为法向摩擦角,η为流屑角,γn为法向前角,可表示为:
γn=arctan(tanγo)cosλs (18)
式中,γ0为刀具前角。
切削力系数标定步骤如图1所示:
步骤1取第i次采集的切削力数据,获得[FR,FA,FT],代入式(10)转换成法向切削力数据[FRn,FAn,FTn]。
步骤2将切削力数据代入式(12),(15)~(17),解出系数ξn、ξλ、βa和βn。
步骤3将系数ξλ和ξn代入式(6)和(7),解出系数
Figure BDA0002768694450000083
Figure BDA0002768694450000084
步骤4由(FRC,FAC,FTC)T求得合力F,代入式(4)求得剪力Fs
步骤5法向剪切角
Figure BDA0002768694450000091
代入式(8)和式(9),可分别解出As和τs
步骤6重复步骤1~5,不少于20组系数βn、
Figure BDA0002768694450000092
和τs
步骤7取出存储系数βn、
Figure BDA0002768694450000093
和τs,分别将其表示成切削层面积SA的函数。
步骤8使用上述量,利用公式(3)确定将切向力系数KR、径向力系数KA和轴向力系数KT表示成切削层面积SA的函数。
本发明说明书中未作详细描述的内容属本领域专业技术人员的公知技术。

Claims (9)

1.一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,包括步骤如下:
1)使用刀具切削圆柱体金属件,采集获得多组切削力试验数据;
2)根据任意一组切削力试验数据,将切削力在前刀面上分解为切向力FT、轴向力FA和径向力FR;然后,根据刀具刃倾角λs,将切向力FT、轴向力FA和径向力FR转换为法向切向力FTn、法向轴向力FRn和法向径向力FAn
3)根据切削力试验数据获得切削合力F,从而确定切削合力F在法平面上的投影与切削平面之间夹角ξn、切削合力F与法平面之间的夹角ξλ、法向摩擦角βn;根据切削合力F,确定剪力Fs,从而确定剪力Fs与法平面之间的夹角
Figure FDA0002768694440000011
剪力Fs在法平面上的投影与切削平面的夹角
Figure FDA0002768694440000012
4)根据刀具刃倾角λs,根据步骤3)所述剪力Fs和剪力Fs在法平面上的投影与切削平面的夹角
Figure FDA0002768694440000013
确定剪应力τs
5)根据步骤3)所述法向摩擦角βn、剪力Fs在法平面上的投影与切削平面的夹角
Figure FDA0002768694440000014
根据步骤4)所述剪应力τs,确定切向力系数KR、径向力系数KA和轴向力系数KT
2.根据权利要求1所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,步骤2)所述将切向力FT、轴向力FA和径向力FR转换为法向切向力FTn、法向轴向力FRn和法向径向力FAn的方法,具体为:
(FTn,FRn,FAn)T=T(λs)(FT,FR,FA)T
Figure FDA0002768694440000015
3.根据权利要求1所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,确定步骤3)所述ξn的方法,具体为:
Figure FDA0002768694440000021
4.根据权利要求2所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,确定步骤3)所述ξλ的方法,具体为:
sinξλ=sinβasinη
tan(ξnn)=tanβacosη
其中,η为流屑角。
5.根据权利要求2所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,确定步骤3)所述法向摩擦角βn的方法,具体为:
βn=ξnn
其中,γn为法向前角。
6.根据权利要求1所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,确定步骤3)所述剪力Fs在法平面投影与切削平面之间夹角
Figure FDA0002768694440000022
的方法,具体为:
Figure FDA0002768694440000023
Figure FDA0002768694440000024
7.根据权利要求3~6任一所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,确定步骤3)所述剪力Fs的方法,具体为:
Figure FDA0002768694440000025
8.根据权利要求7所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,确定步骤4)所述剪应力τs的方法,具体为:
Figure FDA0002768694440000026
Figure FDA0002768694440000027
其中,f为进给量,ap为径向切削深度,As为剪切面积。
9.根据权利要求7所述的一种基于最大剪应力原理的切削力建模系数标定方法,其特征在于,步骤5)所述确定切向力系数KR、径向力系数KA和轴向力系数KT的方法,具体为:
Figure FDA0002768694440000031
其中,η为流屑角,γn为刀具的前角。
CN202011241933.1A 2020-11-09 2020-11-09 一种基于最大剪应力原理的切削力建模系数标定方法 Pending CN112329245A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011241933.1A CN112329245A (zh) 2020-11-09 2020-11-09 一种基于最大剪应力原理的切削力建模系数标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011241933.1A CN112329245A (zh) 2020-11-09 2020-11-09 一种基于最大剪应力原理的切削力建模系数标定方法

Publications (1)

Publication Number Publication Date
CN112329245A true CN112329245A (zh) 2021-02-05

Family

ID=74316554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011241933.1A Pending CN112329245A (zh) 2020-11-09 2020-11-09 一种基于最大剪应力原理的切削力建模系数标定方法

Country Status (1)

Country Link
CN (1) CN112329245A (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107085412A (zh) * 2017-03-30 2017-08-22 北京精密机电控制设备研究所 一种控制曲线车削加工切削力突变的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107085412A (zh) * 2017-03-30 2017-08-22 北京精密机电控制设备研究所 一种控制曲线车削加工切削力突变的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
罗智文: "基于斜角切削的难加工材料曲线加工切削力建模与参数优化", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》, no. 6, pages 022 - 81 *

Similar Documents

Publication Publication Date Title
CN109783968B (zh) 基于多工步的金属切削加工过程的三维有限元仿真方法
CN106897484B (zh) 一种高温合金高效低应力磨削工艺参数的获得方法
Kong et al. Broaching performance of superalloy GH4169 based on FEM
CN102880756A (zh) 薄壁叶片精密铣削加工变形补偿方法
CN109940459B (zh) 一种高效铣刀损伤的多尺度识别方法
CN101722399B (zh) 一种斜面球阀的制造方法
CN107505850B (zh) 一种刀具换刀判断方法
CN107451382B (zh) 高速切削加工工件表面形貌的控制方法
CN101491844B (zh) 圆周铣削过程铣削力系数与刀具径向偏心标定方法
CN101653921A (zh) 凸轮轴数控磨削轮廓误差补偿方法
CN104493636A (zh) 一种用于提高铣削稳定性的金属深冷加工方法
CN109754332B (zh) 基于切削力的机床铣削加工过程的能耗模型建模方法
CN106372368B (zh) 一种提高大型薄壁回转结构件车加工精度的方法
CN109857061A (zh) 一种基于热力耦合模型的工件表面残余应力调控方法
CN113239461B (zh) 一种非对称结构复杂薄壁零件变形控制方法
Tomczak et al. Skrew rolling of balls in multiple helical impressions
CN109299581B (zh) 一种结合曲面插值的端面铣刀铣削力预测方法
CN110681696A (zh) 高含硫钢轧制工艺
CN112329245A (zh) 一种基于最大剪应力原理的切削力建模系数标定方法
Ozel et al. Modeling of high speed machining processes for predicting tool forces, stresses and temperatures using FEM simulations
CN105631072B (zh) 一种铝合金粗加工刀具设计方法
KR100902863B1 (ko) 임펠러의 황삭가공을 위한 공구조합 선정방법
CN105945372A (zh) 一种美制内花键制作方法
CN115472246A (zh) 一种高温合金斜轧穿孔过程离层缺陷的预测方法及模型
Gowd et al. Empirical modelling of hard turning process of Inconel using response surface methodology

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination