CN111077846B - 一种薄壁叶片余量去除顺序规划方法 - Google Patents

一种薄壁叶片余量去除顺序规划方法 Download PDF

Info

Publication number
CN111077846B
CN111077846B CN201911292011.0A CN201911292011A CN111077846B CN 111077846 B CN111077846 B CN 111077846B CN 201911292011 A CN201911292011 A CN 201911292011A CN 111077846 B CN111077846 B CN 111077846B
Authority
CN
China
Prior art keywords
blade
deformation
machining
track
thin
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.)
Active
Application number
CN201911292011.0A
Other languages
English (en)
Other versions
CN111077846A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201911292011.0A priority Critical patent/CN111077846B/zh
Publication of CN111077846A publication Critical patent/CN111077846A/zh
Application granted granted Critical
Publication of CN111077846B publication Critical patent/CN111077846B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/404Numerical 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 control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)

Abstract

本发明涉及一种薄壁叶片余量去除顺序规划方法,属于航空发动机叶片的高品质精密制造技术领域,特别涉及到薄壁叶片余量去除与轨迹规划方法,用于控制薄壁叶片在数控加工中产生的弹性变形误差,以此提高叶片的加工精度。该方法针对加工过程中叶片曲面在不同位置的刚度变化,利用弹性变形原理与法向变形场,给出不同余量去除顺序的叶片变形程度的评定方法,并结合加工轨迹的连续性要求,给出变形最小的加工轨迹。采用该方法生成的加工轨迹进行加工后,叶片的加工误差有了明显的下降,提高了加工精度。

Description

一种薄壁叶片余量去除顺序规划方法
技术领域
本发明属于航空发动机叶片高品质精密制造技术领域,特别涉及到薄壁叶片余量去除与轨迹规划方法,用于控制薄壁叶片在数控加工中产生的弹性变形误差,以此提高叶片的加工精度。
背景技术
叶片类零件是航空发动机的关键零件,具有外形复杂、易变形、加工难度大等特点。尤其是叶片型面的自由曲面与整体弱刚度的相互耦合,使得叶片在铣削过程中很难保证加工精度。
专利CN201410360136.3公开了一种薄壁件阶梯对称铣削加工方法,通过对薄壁类结构件采用阶梯状的非对称式铣削加工,保证未切削区域对切削区域的支撑,达到提高工件加工过程中的刚度,减小加工变形。然而该方法并不适用于自由曲面加工。专利CN201410526635.5公开了一种利用球头刀铣削自由曲面的加工方法,以张量的方式度量曲面上任意切触点的切削宽度,以此提高加工效率。但该方法并没有考虑叶片类自由曲面的弱刚性问题,没有给出自由曲面类薄壁零件的切削轨迹规划方法。
因此,发明一种薄壁叶片余量去除顺序规划方法显得非常必要。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种薄壁叶片余量去除顺序规划方法,通过叶片建模仿真技术获取叶片整体的刚度矩阵并计算目标曲面的变形场,同时将法向变形场作为动态加工过程中的变形计算基准,以此评定不同加工轨迹变形程度,最后结合加工轨迹的连续性要求给出变形最小的加工轨迹。
技术方案
一种薄壁叶片余量去除顺序规划方法,其特征在于步骤如下:
步骤1:将叶根装夹的叶片进行有限元划分,得到叶片的刚度矩阵KG
步骤2:依据叶片整体的刚度矩阵,计算待加工曲面的变形场;
步骤2.1:将叶片待加工曲面的有限元节点作为离散点{P1,P2,…,Pm},并在离散点对应的法向N=[nx,ny,nz]T施加单位力进行节点的位移测试,计算公式如下:
KG(δ(P1),δ(P2),…,δ(Pm))=(F(P1),F(P2),…,F(Pm))
其中,δ(Pm)表示为所有的节点在单位力F(Pm)的作用下的位移向量;
步骤2.2:选取离散点{P1,P2,…,Pm}对应的变形量{δ12,…,δm},依据叶片曲面在弧长参数域的(u,v)方向,采用双线性插值计算待加工曲面的法向变形场Gn
步骤2.3:考虑到变形误差被认定为法线方向的位移,变形误差的法向变形场表示为Gn,n
Gn,n=NTGn
步骤3:评定给定的加工轨迹变形程度,并规划最优的加工轨迹;
步骤3.1:将离散点{P1,P2,…,Pm}按照给定的加工轨迹组成与加工时刻s(s=0,1,…,m)相关的序列,表示为:
CP=[Pi,Pj,…,Pk]
其中,CP为加工轨迹序列;同时CP(s)=Pi,用s=0表示未加工;
步骤3.2:利用变形误差的法向变形场Gn,n计算不同加工时刻s下的局部变形量q(s),公式如下,
Figure BDA0002319361710000021
步骤3.3:根据局部变形量计算全过程变形量Q,公式如下:
Figure BDA0002319361710000022
步骤4:对于不同的加工轨迹,相应的序列CP对应的全过程变形量Q最小的加工轨迹是加工变形最小的轨迹,且叶片类薄壁零件的最优加工轨迹就是从叶尖向叶根行进给的平行切削轨迹。
有益效果
本发明提出的一种薄壁叶片余量去除顺序规划方法,该方法针对加工过程中叶片曲面在不同位置的刚度变化,利用弹性变形原理与法向变形场,给出不同余量去除顺序的叶片变形程度的评定方法,并结合加工轨迹的连续性要求,给出变形最小的加工轨迹。采用该方法生成的加工轨迹进行加工后,叶片的加工误差有了明显的下降,提高了加工精度。
采用该方法得到的薄壁叶片余量去除顺序可有有效的减小加工过程中的弹性变形误差,相比于平行于叶尖与叶根方向的平行切削轨迹,从叶尖向叶根行进给的平行切削轨迹的平均误差与最大误差分别减小了55.8%与71.7%,提高了加工精度。
附图说明
图1从叶尖向叶根行进给的平行切削轨迹
图2平行于叶尖与叶根方向的平行切削轨迹
具体实施方式
现结合实施例、附图对本发明作进一步描述:
步骤1.本实施例中,将50×40×1mm叶片在叶根装夹,并进行有限元划分,得到叶片的刚度矩阵KG
步骤2.依据叶片整体的刚度矩阵,计算待加工曲面的变形场。
步骤2.1将叶片待加工曲面的有限元节点作为离散点{P1,P2,…,Pm},并在离散点对应的法向N=[nx,ny,nz]T施加单位力进行节点的位移测试,计算公式如下,
KG(δ(P1),δ(P2),…,δ(Pm))=(F(P1),F(P2),…,F(Pm))
其中,δ(Pm)表示为所有的节点在单位力F(Pm)的作用下的位移向量。
步骤2.2选取离散点{P1,P2,…,Pm}对应的变形量{δ12,…,δm},依据叶片曲面在弧长参数域的(u,v)方向,采用双线性插值计算待加工曲面的法向变形场Gn
步骤2.3考虑到变形误差被认定为法线方向的位移,变形误差的法向变形场表示为Gn,n
Gn,n=NTGn
步骤3.评定给定的加工轨迹变形程度,并规划最优的加工轨迹。
步骤3.1.将离散点{P1,P2,…,Pm}按照给定的加工轨迹组成与加工时刻s(s=0,1,…,m)相关的序列,表示为,
CP=[Pi,Pj,…,Pk]
其中,CP为加工轨迹序列。同时CP(s)=Pi,用s=0表示未加工。
步骤3.2.利用变形误差的法向变形场Gn,n计算不同加工时刻s下的局部变形量q(s),公式如下,
Figure BDA0002319361710000041
步骤3.3.根据局部变形量计算全过程变形量Q,公式如下,
Figure BDA0002319361710000042
步骤4.对于不同的加工轨迹,相应的序列CP对应的全过程变形量Q最小的加工轨迹是加工变形最小的轨迹。分别计算从叶尖向叶根行进给的平行切削轨迹与平行于叶尖与叶根方向的平行切削轨迹的全过程变形量Q1与Q2,分别是Q1=10.75与Q2=25.47。为了使全过程变形量Q最小,同时考虑到轨迹的连续性要求,叶片类薄壁零件的最优加工轨迹就是全过程变形量Q1=10.75的从叶尖向叶根行进给的平行切削轨迹。
本实施例中,在三轴数控机床上对叶片采用两种切削轨迹进行加工,并通过检测得到,相比于平行于叶尖与叶根方向的平行切削轨迹,从叶尖向叶根行进给的平行切削轨迹的平均误差与最大误差分别减小了55.8%与71.7%。

Claims (1)

1.一种薄壁叶片余量去除顺序规划方法,其特征在于步骤如下:
步骤1:将叶根装夹的叶片进行有限元划分,得到叶片的刚度矩阵KG
步骤2:依据叶片整体的刚度矩阵,计算待加工曲面的变形场;
步骤2.1:将叶片待加工曲面的有限元节点作为离散点{P1,P2,…,Pm},并在离散点对应的法向N=[nx,ny,nz]T施加单位力进行节点的位移测试,计算公式如下:
KG(δ(P1),δ(P2),…,δ(Pm))=(F(P1),F(P2),…,F(Pm))
其中,δ(Pm)表示为所有的节点在单位力F(Pm)的作用下的位移向量;
步骤2.2:选取离散点{P1,P2,…,Pm}对应的变形量{δ1,δ2,…,δm},依据叶片曲面在弧长参数域的(u,v)方向,采用双线性插值计算待加工曲面的法向变形场Gn
步骤2.3:考虑到变形误差被认定为法线方向的位移,变形误差的法向变形场表示为Gn,n
Gn,n=NTGn
步骤3:评定给定的加工轨迹变形程度,并规划最优的加工轨迹;
步骤3.1:将离散点{P1,P2,…,Pm}按照给定的加工轨迹组成与加工时刻s相关的序列,s=0,1,…,m,表示为:
CP=[Pi,Pj,…,Pk]
其中,CP为加工轨迹序列;同时CP(s)=Pi,用s=0表示未加工;
步骤3.2:利用变形误差的法向变形场Gn,n计算不同加工时刻s下的局部变形量q(s),公式如下,
Figure FDA0003229841790000011
步骤3.3:根据局部变形量计算全过程变形量Q,公式如下:
Figure FDA0003229841790000012
步骤4:对于不同的加工轨迹,相应的序列CP对应的全过程变形量Q最小的加工轨迹是加工变形最小的轨迹,且叶片类薄壁零件的最优加工轨迹就是从叶尖向叶根行进给的平行切削轨迹。
CN201911292011.0A 2019-12-16 2019-12-16 一种薄壁叶片余量去除顺序规划方法 Active CN111077846B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911292011.0A CN111077846B (zh) 2019-12-16 2019-12-16 一种薄壁叶片余量去除顺序规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911292011.0A CN111077846B (zh) 2019-12-16 2019-12-16 一种薄壁叶片余量去除顺序规划方法

Publications (2)

Publication Number Publication Date
CN111077846A CN111077846A (zh) 2020-04-28
CN111077846B true CN111077846B (zh) 2021-10-15

Family

ID=70314663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911292011.0A Active CN111077846B (zh) 2019-12-16 2019-12-16 一种薄壁叶片余量去除顺序规划方法

Country Status (1)

Country Link
CN (1) CN111077846B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112192150B (zh) * 2020-09-21 2022-11-04 贵州新安航空机械有限责任公司 一种开式叶轮加工工艺
CN112861070A (zh) * 2021-01-29 2021-05-28 苏州千机智能技术有限公司 整体叶环叶片上多余材料区域的实际加工曲面重构方法及系统
CN114799357B (zh) * 2022-05-26 2022-11-29 成都飞机工业(集团)有限责任公司 一种腹板加工方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8545142B2 (en) * 2008-03-06 2013-10-01 University Of North Carolina At Charlotte Deformation machining systems and methods
CN103128312B (zh) * 2011-12-05 2016-06-01 沈阳黎明航空发动机(集团)有限责任公司 一种薄壁盘零件控制变形加工工艺
CN103084639A (zh) * 2013-01-21 2013-05-08 西北工业大学 一种基于非均匀余量的增强薄壁叶片工艺刚度的方法
CN103990840B (zh) * 2014-04-22 2016-05-04 西北工业大学 叶片双向变余量铣削方法
CN104227103A (zh) * 2014-07-24 2014-12-24 西安航空学院 薄壁件阶梯对称铣削加工方法

Also Published As

Publication number Publication date
CN111077846A (zh) 2020-04-28

Similar Documents

Publication Publication Date Title
CN111077846B (zh) 一种薄壁叶片余量去除顺序规划方法
CN111159825B (zh) 一种薄壁叶片切削轨迹参数优化方法
TWI414376B (zh) 一種五軸曲面側銑加工系統及其刀具路徑規劃方法
US8103375B2 (en) Fillet machining system
CN107728577B (zh) 基于薄壁曲面加工变形的瞬时切削量规划方法
CN103163836B (zh) 直纹曲面刀具轨迹规划方法
CN103212760B (zh) 一种等截面叶型带冠整体叶轮叶间流道加工方法
Ayesta et al. Optimum electrode path generation for EDM manufacturing of aerospace components
CN102581705B (zh) 一种铣削加工工件误差的分度补偿方法
CN112255966A (zh) 一种窄长类自由曲面零件加工轨迹自适应生成方法
CN110989490B (zh) 一种基于轮廓误差的工件最优安装位置的获取方法
CN108581384A (zh) 一种基于UG和Vericut的整体式叶轮四轴车铣加工方法
CN113239461A (zh) 一种非对称结构复杂薄壁零件变形控制方法
CN101249571A (zh) 多曲面岛屿五轴螺旋加工方法
CN109590559B (zh) 一种五轴电火花加工闭式叶轮轨迹规划方法
CN114888702A (zh) 一种压气机叶片数控抛光方法
CN103419540A (zh) 基于路径单元的变形曲面自适应投影加工方法
CN113996867A (zh) 基于变余量补偿的焊接类叶盘夹持残留区域光顺加工方法
CN105127696B (zh) 一种汽轮机调节级叶片加工方法
CN111063020B (zh) 基于PowerMill软件的工件三维轮廓加工方法
CN109605120B (zh) 一种提高航空叶片加工良品率的方法
CN113835397B (zh) 基于b样条曲线和路径积分的线性数控加工路径平滑方法
CN112296766B (zh) 一种整体叶盘叶片型面的数控柔性抛光方法
CN109765843B (zh) 基于镜像法和三次样条插值的薄壁件加工误差补偿方法
CN111026035B (zh) 一种基于曲率变化求解旋风铣削叶片刀位点的方法

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