CN107831731B - 一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法 - Google Patents
一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法 Download PDFInfo
- Publication number
- CN107831731B CN107831731B CN201711049865.7A CN201711049865A CN107831731B CN 107831731 B CN107831731 B CN 107831731B CN 201711049865 A CN201711049865 A CN 201711049865A CN 107831731 B CN107831731 B CN 107831731B
- Authority
- CN
- China
- Prior art keywords
- cutting
- layer
- radius
- tool path
- arc
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/19—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/35—Nc in input of data, input till input file format
- G05B2219/35349—Display part, programmed locus and tool path, traject, dynamic locus
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,具体包括步骤一:输入刀轨和工件信息,获取外拐角参数;步骤二:优化外拐角最内层切削层;步骤三:优化剩余切削层;步骤四:计算另一侧切削层轮廓半径;步骤五:确定所有切削刀轨的的位置和起终点;步骤六:加入非切削移动刀轨,生成完整的局部循环;步骤七:结合原始刀轨,输出优化后的刀轨。本发明采用局部循环分层切削的加工方式,使得切削力热均衡化,减小了载荷冲击;本发明优化刀轨由圆弧与直线组成,切入切出平滑过渡,有效减小振动,提高加工表面质量;本发明通过约束最大切宽,避免了力热载荷过大。
Description
技术领域
本发明涉及一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,属于数控加工技术领域。
背景技术
在飞机壁板、梁、框以及模具型腔加工中,广泛存在拐角特征的加工。在外拐角的加工过程中,由于切入切出角的变化,使得拐角处的实际切削宽度大于直边的设计切宽。拐角处刀轨突变,切屑厚度与径向切深增大,造成切削力的突然增大,振动加剧,对刀具造成冲击,易产生崩刃断刀现象,加工质量降低,零件精度不能满足要求。为避免切削用量增加,切削力增大所产生的影响,需要对拐角处实际切削力进行仿真并进行刀轨优化。
拐角加工中,切削力直接影响刀具磨损和已加工表面质量,也是引发振动和颤振的重要原因,基于切削力仿真可以在数控加工之前对刀轨、切削参数进行优化。Xiong Han在学术期刊《International Journal of Machine Tools&Manufacture》发表论文“Precise prediction of forces in milling circular corners”研究了拐角加工中,瞬时切屑厚度变化,并基于单位切削力系数的力学模型,进行切削力仿真计算。但是,该方法需要数值迭代求解,而且并没有基于切削力仿真对切削用量进行约束优化。Shaochun Sui在学术期刊《International Journal of Advanced Manufacturing Technology》发表论文“Tool path generation and optimization method for pocket flank milling ofaircraft structural parts based on the constraints of cutting force anddynamic characteristics of machine tools”提出了一种变螺旋曲线的型腔加工刀规优化方法,该方法考虑了切削力与机床动态特性的约束,但需在拐角处单独循环走刀,且仅考虑了内拐角特征的加工。
因此,针对型腔外拐角数控铣削加工,需要对外拐角切削力进行仿真预适应,对刀轨形状进行优化,提出一种新的外拐角连续刀轨生成方法。
发明内容
本发明的目的是为了解决上述问题,提出一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,对于拐角处的切削过程,本方法基于瞬时切削力模型,仿真拐角加工过程中的最大瞬时切削力,在拐角加工中,最大瞬时切削力随刀轨实际切削宽度的变化而变化。故利用拐角切削过程中的最大实际切宽来评价拐角加工过程中切削力的变化程度,并以最大实际切宽小于等于名义切宽为优化目标,来分层优化拐角处的局部刀轨,使得加工过程中切削力变化均匀,避免切削力的过度变大以及由切宽变大引起的颤振现象。
一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,具体的优化步骤如下;
步骤一:输入刀轨和工件信息,获取外拐角参数;
步骤二:优化外拐角最内层切削层;
步骤三:优化剩余切削层;
步骤四:计算另一侧切削层轮廓半径;
步骤五:确定所有切削刀轨的的位置和起终点;
步骤六:加入非切削移动刀轨,生成完整的局部循环;
步骤七:结合原始刀轨,输出优化后的刀轨。
本发明的优点在于:
(1)采用局部循环分层切削的加工方式,使得切削力热均衡化,减小了载荷冲击;
(2)优化刀轨由圆弧与直线组成,切入切出平滑过渡,有效减小振动,提高加工表面质量;
(3)通过约束最大切宽,避免了力热载荷过大;
附图说明
图1是本发明的一种切削力仿真预适应的型腔内拐角数控铣削加工刀轨优化方法流程图。
图2是内拐角特征加工示意图;
图3是本发明的外拐角刀轨最后一层加工示意图。
图4是本发明的外拐角刀轨剩余切削层加工示意图。
图5是本发明的外拐角切削层圆弧半径计算示意图。
图6是本发明的外拐角刀轨位置点计算示意图。
图2、3、4、5、6中R0为加工前拐角圆弧半径,Rc为加工后拐角圆弧半径,R为刀具半径,拐角夹角为2θ,B点、D点为拐角直线圆弧相接点,ae为名义切宽,aem为最大实际切宽。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,流程如图1 所示,具体的优化步骤如下。
步骤一:输入刀轨和工件信息,获取外拐角参数;
假设切削过程中切深和进给速度保持不变,根据刀轨和工件信息得到拐角参数,如图2 所示。拐角夹角为2θ,加工前圆弧半径R0,拐角圆弧半径Rc,刀具半径R,刀轨名义切宽ae, L为切入刀轨直线长度,即直线AB的长度。切宽分层阈值λ,λ≥1,Rc>R。
以内轮廓圆弧圆心为原点、以拐角角平分线为Y轴作直角坐标系,则拐角部分起点为:
终点为:
通过以上计算,可得到外拐角加工,输入数控系统的G代码如下:
步骤二:优化外拐角最内层切削层;
按照实际的刀轨和加工前轮廓进行切削时,刀尖在对角线上时实际切宽达到最大,实际最大切宽aem为:
若aem>λae,则需要对刀轨进行优化。为保证分层切削后最后一层切削时最大切削力与 DE,D′E′保持一致,故采用如附图3所示的分层方法,首先将切削层分为两层,即最后切削层与第一个剩余切削层。最后切削层切削前左侧轮廓与D′E′相切,右侧轮廓与AB相切,且左右轮廓圆弧圆心同在Y轴,两圆弧相切于Y轴上一点。其中拐角夹角为2θ,加工前圆弧半径R0,拐角圆弧半径Rc,刀具半径R,刀轨名义切宽ae。同时,约束最大切宽aem等于名义切宽ae。根据拐角半径Rc确定最后一层切削前D′E′侧轮廓半径R′c,解下面关于R′c的一元二次方程组
一元二次方程的判别式为
方程解为
若aem≤λae,则不需要优化,转步骤五。
步骤三:优化剩余切削层;
步骤二优化之后剩余的切削部分采用由外而内的分层方式,以保证第1层切削层与进入拐角前部分的铣削力变化均匀。如附图4所示。
第i个剩余切削层在D′E′侧由Ri-1,R′c组成月牙形切削层(i=1,2,…),假设计算第i 层剩余切削层最大实际切宽aei,由下式计算当前剩余切削层实际最大切宽
当aei>λae时,新增第i个切削层,构造外侧轮廓半径为Ri-1,内轮廓半径为Ri的月牙形切削层,根据当前层外侧轮廓半径Ri-1和切宽ae求得内轮廓圆弧半径Ri,关于Ri的一元二次方程组如下
解一元二次方程组可得当前层内轮廓圆弧半径Ri,更新剩余切削层最大切宽aei直至 aei≤λae,不再新增切削层,此时共新增切削层m个。
最终结合步骤二的R′c轮廓,步骤二、三共新增了m+1个切削层,最终给出新增的切削层轮廓半径R1,R2,…,Rn,其中Rn=R′c,n=m+1。
步骤四:计算另一侧切削层轮廓半径;
n+1个新增的半径为R1,…,Rn圆弧为与D′E′相切的圆弧半径,还需要求出与对应的与AB 相切的圆弧半径。如附图5所示
设第k层切削层D′E′侧圆弧半径为Rk,(k=1,2,…,n),AB侧圆弧半径为Rk0,且Rk0圆弧 AB边和Rk圆弧同时相切,可得
由此可得AB侧圆弧半径R10,…,Rn0。
步骤五:确定所有切削刀轨的的位置和起终点;
如附图6所示,对于第k层切削刀轨,由D′E′侧切削圆弧半径Rk,AB侧切削圆弧半径Rk0以及刀具半径R可得:
终点为
将R1,…,Rn代入可得新增的各切削层的切削刀轨。
步骤六:加入非切削移动刀轨,生成完整的局部循环;
第k层切削刀轨均比第1层靠内,故非切削移动的圆弧采用与第1层相同的轨迹进行,其余部分以直线移动;然后与各切削层切削部分起终点相连接。
步骤七:结合原始刀轨,输出如下优化后的刀轨;
输入数控系统的G代码如下:
…
…
实施例:
如图所示,对于拐角特征为半径为30mm的外拐角,拐角角度为45°,加工前为圆弧轮廓,半径为29mm,刀轨切宽为5mm,刀具直径为D20的立铣刀,直线刀轨长20的刀轨实例,根据上述算法进行切削力均衡刀轨优化。
步骤一:输入刀轨和工件信息,获取外拐角参数;
根据刀轨和工件信息得到拐角夹角为2θ=45°,加工前圆弧半径R0=29,拐角圆弧半径 Rc=30,刀具半径R=10,刀轨名义切宽ae=5,直线刀轨部分长度L=20,λ=1。
原刀轨信息如下:
G01X44.609Y-3.170F1000
G01X36.955Y15.307
G03X-36.955Y15.307I-36.955J15.307
G01X-44.609Y-3.170
步骤二:优化外拐角最内层切削层;
按照实际的刀轨和加工前轮廓进行切削时,刀具在对角线上时实际切宽达到最大,实际最大切宽aem
由于aem>ae,则需要对刀轨进行优化。为保证分层切削后最后一层切削时最大切削力与 DE,D′E′保持一致,故采用如附图2所示的分层方法,最后一层切削前轮廓与D′E′相切,根据拐角半径Rc、ae确定最后一层切削前D′E′侧轮廓半径R′c,解下面关于R′c的一元二次方程组
解方程组(6)得R′c=34.346
步骤三:优化剩余切削层;
剩余的切削部分采用由外而内的分层方式,以保证第1层切削层与进入拐角前部分的铣削力变化均匀。
1.第1层剩余切削层
第1个剩余切削层在D′E′侧由R0,R′c组成月牙形切削层,计算第1层剩余切削层最大实际切宽ae1
新增第1个切削层,构造外侧轮廓半径为R0,内轮廓半径为R1的月牙形切削层,有关于 R1的一元二次方程组如下
解一元二次方程组(7)得R1=32.708
2.第2层剩余切削层
第2个剩余切削层在D′E′侧由R1,R′c组成月牙形切削层,计算第2层剩余切削层最大实际切宽ae2
故对剩余切削层的分层完成,共新增了m=1层,结合R′c,步骤1,2共新增了n=2个切削层,新增的切削层轮廓半径如下
R1=32.708
R2=34.346
步骤四:削层轮廓半径;
4个新增的半径为R1,…,R4圆弧为与D′E′相切的圆弧半径,还需要求出与对应的与AB相切的圆弧半径。
第1层切削层AB侧圆弧半径
第2层切削层AB侧圆弧半径
步骤五:确定所有切削刀轨的的位置和起终点;
对于第k层切削刀轨,由D′E′侧切削圆弧半径Rk,AB侧切削圆弧半径Rk0以及刀具半径R可得:
将R1,…,Rn代入可得新增的各切削层的切削刀轨。
步骤六:加入非切削移动刀轨,生成完整的局部循环;
第k层切削刀轨均比第1层靠内,故非切削移动的圆弧采用与第1层相同的轨迹进行,其余部分以直线移动;然后与各切削层切削部分起终点相连接。
步骤七:结合原始刀轨,输出如下优化后的刀轨;
输入数控系统的G代码如下:
G01 X[(Rc-R)cosθ+Lsinθ]Y[(Rc-R)sinθ-Lcosθ]Ffeedrate
…
…
G01 X[(Rc-R)cosθ]Y[(Rc-R)sinθ]
G03 X[-(Rc-R)cosθ]Y[(Rc-R)sinθ]I[-(Rc-R)cosθ]J[-(Rc-R)sinθ]
G01 X[-(Rc-R)cosθ-Lsinθ]Y[(Rc-R)sinθ-Lcosθ] 。
Claims (5)
1.一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,具体的优化步骤如下;
步骤一:输入刀轨和工件信息,获取外拐角参数;
具体为:
假设切削过程中切深和进给速度保持不变,根据刀轨和工件信息得到拐角参数,设拐角夹角为2θ,加工前圆弧半径R0,拐角圆弧半径Rc,刀具半径R,刀轨名义切宽ae,L为切入刀轨直线长度;切宽分层阈值λ,λ≥1,Rc>R;
以外轮廓圆弧圆心为原点、以拐角角平分线为Y轴作直角坐标系,则拐角部分起点为:
终点为:
步骤二:优化外拐角最内层切削层;
步骤三:优化剩余切削层;
步骤四:计算另一侧切削层轮廓半径;
具体为:
n+1个新增的半径为R1,…,Rn圆弧为与D′E′相切的圆弧半径,还需要求出对应的与AB相切的圆弧半径;
设第k层切削层D′E′侧圆弧半径为Rk,k=1,2,…,n,AB侧圆弧半径为Rk0,且Rk0圆弧AB边和Rk圆弧同时相切,可得
由此可得AB侧圆弧半径R10,…,Rn0;
步骤五:确定所有切削刀轨的的位置和起终点;
步骤六:加入非切削移动刀轨,生成完整的局部循环;
步骤七:结合原始刀轨,输出优化后的刀轨。
2.根据权利要求1所述的一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,所述的步骤二具体为:
按照实际的刀轨和加工前轮廓进行切削时,刀具在对角线上时实际切宽达到最大,实际最大切宽aem为:
若aem>λae,对刀轨进行优化,首先将切削层分为两层,即最后切削层与第一个剩余切削层;最后切削层切削前左侧轮廓与切削前退刀侧内轮廓左侧直线段相切,右侧轮廓与切削前退刀侧外轮廓右侧直线段相切,且左右轮廓圆弧圆心同在Y轴,两圆弧相切于Y轴上一点;其中拐角夹角为2θ,加工前圆弧半径R0,拐角圆弧半径Rc,刀具半径R,刀轨名义切宽ae;同时,约束最大切宽aem等于名义切宽ae;根据拐角半径Rc、确定最后一层切削前退刀侧内轮廓圆弧半径R′c,解下面关于R′c的一元二次方程组
一元二次方程的判别式为
方程解为
若aem≤λae,则不需要优化,转步骤五。
3.根据权利要求1所述的一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法,所述的步骤三具体为:
步骤二优化之后剩余的切削部分采用由外而内的分层方式;
第i个剩余切削层在D′E′侧由Ri-1,R′c组成月牙形切削层,假设计算第i层剩余切削层最大实际切宽aei,由下式计算当前剩余切削层实际最大切宽
当aei>λae时,新增第i个切削层,构造外侧轮廓半径为Ri-1,内轮廓半径为Ri的月牙形切削层,根据当前层外侧轮廓半径Ri-1和切宽ae求得内轮廓圆弧半径Ri,关于Ri的一元二次方程组如下
解一元二次方程组可得当前层内轮廓圆弧半径Ri,更新剩余切削层最大切宽aei直至aei≤λae,不再新增切削层,此时共新增切削层m个;
最终结合步骤二的R′c轮廓,步骤二、三共新增了m+1个切削层,最终给出新增的切削层轮廓半径R1,R2,…,Rn,其中Rn=R′c,n=m+1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711049865.7A CN107831731B (zh) | 2017-10-31 | 2017-10-31 | 一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711049865.7A CN107831731B (zh) | 2017-10-31 | 2017-10-31 | 一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107831731A CN107831731A (zh) | 2018-03-23 |
CN107831731B true CN107831731B (zh) | 2020-02-18 |
Family
ID=61651232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711049865.7A Expired - Fee Related CN107831731B (zh) | 2017-10-31 | 2017-10-31 | 一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831731B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112464469B (zh) * | 2020-11-25 | 2022-05-17 | 北京航空航天大学 | 用于加工研究的特征聚合零件结构及加工工艺优化方法 |
CN117170306B (zh) * | 2023-10-31 | 2024-03-15 | 成都飞机工业(集团)有限责任公司 | 一种刀具状态监控方法、系统、电子设备及存储介质 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184276B (zh) * | 2011-03-21 | 2013-02-06 | 西安交通大学 | 一种切削加工中进给速度的优化方法 |
CN104007696B (zh) * | 2014-03-19 | 2017-01-11 | 沈阳飞机工业(集团)有限公司 | 飞机结构件型腔复杂转角残留区域刀轨生成方法 |
CN104759942B (zh) * | 2015-04-22 | 2018-06-26 | 华中科技大学 | 一种薄壁零件的铣削变形在线测量与补偿加工方法 |
CN105160128A (zh) * | 2015-09-25 | 2015-12-16 | 北京理工大学 | 一种曲线端铣加工过程切削力的预测方法 |
CN105234466B (zh) * | 2015-10-28 | 2018-03-16 | 成都飞机工业(集团)有限责任公司 | 一种槽特征侧铣加工刀轨生成方法 |
CN105701323A (zh) * | 2016-03-23 | 2016-06-22 | 华中科技大学 | 一种插铣加工切削力预测建模方法 |
CN105955195B (zh) * | 2016-05-16 | 2018-05-22 | 哈尔滨理工大学 | 一种基于铣削力预测的曲面加工轨迹生成方法 |
CN106216747A (zh) * | 2016-07-22 | 2016-12-14 | 河北师范大学 | 一种整体叶轮五轴联动数控刀轨路径加工方法 |
-
2017
- 2017-10-31 CN CN201711049865.7A patent/CN107831731B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107831731A (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103064344B (zh) | 一种基于nurbs曲线插补的速度平滑控制方法 | |
CN101976060B (zh) | 基于机床动力学和曲线特性的nurbs插补方法 | |
Ming et al. | Four-axis trochoidal toolpath planning for rough milling of aero-engine blisks | |
KR101928419B1 (ko) | 공구경로의 생성방법 및 공작기계 | |
Lai et al. | On the development of a parametric interpolator with confined chord error, feedrate, acceleration and jerk | |
CA2747118C (en) | Method of curvature controlled data smoothing | |
Sun et al. | A real-time and look-ahead interpolation methodology with dynamic B-spline transition scheme for CNC machining of short line segments | |
CN105234466B (zh) | 一种槽特征侧铣加工刀轨生成方法 | |
Zhu et al. | Simultaneous optimization of tool path and shape for five-axis flank milling | |
CN105700466A (zh) | 高速数控加工轨迹的曲率光顺方法 | |
CN107831731B (zh) | 一种切削力仿真预适应的型腔外拐角数控铣削加工刀轨优化方法 | |
CN107831730B (zh) | 一种切削力仿真预适应的型腔内拐角数控铣削加工刀轨优化方法 | |
CN107765648A (zh) | 一种cnc加工的进给速度规划方法及装置 | |
CN104281099A (zh) | 一种考虑加工特性的nurbs直接插补方法及装置 | |
Banerjee et al. | Process planning for Floor machining of 2½D pockets based on a morphed spiral tool path pattern | |
CN106843146B (zh) | 一种自适应变增益轮廓误差补偿方法 | |
CN105005265A (zh) | 一种基于曲线光滑变形的数控加工进给率规划方法 | |
CN106094737B (zh) | 一种指定加工误差条件下的数控加工速度优化控制方法 | |
CN109918807A (zh) | 一种优化曲率的局部刀轨光顺方法 | |
Zhaoyu et al. | A variable-depth multi-layer five-axis trochoidal milling method for machining deep freeform 3D slots | |
Wei et al. | Modeling and machining of integral impeller based on NURBS curve | |
CN103136426B (zh) | 航空叶片圆弧形前后缘工艺模型生成方法 | |
TW200821784A (en) | Processing path method of converting straight line to NURBS curve of computer numerical control | |
CN113467376B (zh) | 一种面向多加工场景的多轴轨迹压缩方法 | |
CN113835397B (zh) | 基于b样条曲线和路径积分的线性数控加工路径平滑方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200218 Termination date: 20201031 |
|
CF01 | Termination of patent right due to non-payment of annual fee |