CN107967384A - 一种基于二维仿真分析的铣削残余应力精确获取方法 - Google Patents

一种基于二维仿真分析的铣削残余应力精确获取方法 Download PDF

Info

Publication number
CN107967384A
CN107967384A CN201711189158.8A CN201711189158A CN107967384A CN 107967384 A CN107967384 A CN 107967384A CN 201711189158 A CN201711189158 A CN 201711189158A CN 107967384 A CN107967384 A CN 107967384A
Authority
CN
China
Prior art keywords
residual stress
milling
simulation
model
triangle
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
CN201711189158.8A
Other languages
English (en)
Other versions
CN107967384B (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201711189158.8A priority Critical patent/CN107967384B/zh
Publication of CN107967384A publication Critical patent/CN107967384A/zh
Application granted granted Critical
Publication of CN107967384B publication Critical patent/CN107967384B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Abstract

本发明提供了一种基于二维仿真分析的铣削残余应力精确获取方法,包括:步骤(1)根据实际情况确定所需仿真参数;步骤(2)建立二维铣削残余应力仿真模型;步骤(3)初次计算铣削仿真残余应力;步骤(4)多次迭代计算获取最终铣削残余应力。本发明能够降低铣削残余应力仿真所需硬件要求,提高仿真效率,减少仿真计算时间,且与常规二维仿真相比,能获得更为精确的仿真结果。

Description

一种基于二维仿真分析的铣削残余应力精确获取方法
技术领域
本发明涉及金属切削加工过程中残余应力获取的技术领域,具体涉及一种基于二维仿真分析的铣削残余应力精确获取方法,适用于各类材料的残余应力铣削仿真试验。
背景技术
金属在切削加工过程中不可避免地引入加工残余应力。加工残余应力是衡量工件表面质量的一项重要指标,严重影响着零件的静力强度、疲劳强度及抗应力腐蚀能力,进而影响零件的使用寿命。同时,残余应力还影响着零件的形状精度和尺寸稳定性,在航空薄壁件的加工中,铣削加工残余应力的存在及其在后续使用过程中的不稳定性会造成应力松弛和再分布,从而引起构件的加工变形。鉴于残余应力的性质、大小及其分布对零件使用性能的重要影响,对加工后残余应力的分布规律进行定量分析具有重要意义。然而现有的研究主要集中在车削、磨削加工残余应力,对铣削加工残余应力研究很少,迫切需要一种铣削残余应力的精确获取方法。
发明内容
本发明主要解决的技术问题为:主流的铣削仿真过程中,三维切削仿真方法计算量大,耗时长,需要浪费极大的资源配置,并不符合现有阶段的科研设计节奏;二维仿真数值提取的区域处于切屑区域,仿真结果与铣削工件表面完整性相差较大,不能满足现有科研设计需求的问题。本发明考虑铣削工艺特点,结合二维铣削仿真特性,提出一种基于二维仿真分析的铣削残余应力精确获取方法,实现了计算精度和计算代价的双重考量。
本发明的目的是针对现有技术的不足,提供一种基于二维仿真分析的铣削残余应力精确获取方法,该方法步骤如下:
步骤(1)、根据实际情况确定所需仿真参数,需要确定的仿真参数为加工刀具的前角α、前角β、刀具半径r,被加工工件材料的本构模型、密度、热导率、杨氏模量、热容、泊松比,切削的铣削深度a、铣削宽度b;
步骤(2)、建立二维铣削残余应力仿真模型,需建立仿真模型分为刀具模型和工件模型,其中刀具模型为三角形,其中一个锐角角度为90-α-β,工件模型是三角形和矩形的组合,上部分为直角三角形,直角边长度分别为铣削宽度b(mm)和c(mm),下部分为矩形,高度为1(mm),宽度为c(mm),其中
步骤(3)、初次计算铣削仿真残余应力,需按照铣削角度要求摆放刀具,刀具刀尖沿工件三角形与矩形的交界线运动,切削深度深度由b(mm)逐渐减为0(mm);
步骤(4)、多次迭代计算获取最终铣削残余应力,需根据需求确定迭代次数N,每次迭代时残余应力从工件矩形区域中选取,第n次迭代选取的点距三角形与矩形区域交界的直角边距离三角形区域另一直角边将提取的残余应力作为初始应力代入下次计算中,直至N次迭代完成。
本发明的原理:使用多次迭代的方法,将某次仿真计算的残余应力作为下一次仿真的初始应力带入计算,以此提高二维铣削仿真的计算精度。
本发明与现有技术相比的有益效果是:首先,相较于其他铣削残余应力数值计算方法,本发明采用有限元仿真进行计算,使用难度较低,不需要很强的理论知识,且能够适应更加复杂的加工条件;其次,相较于其他二维铣削残余应力计算方法,本发明提供方法拥有更高的计算精度;最后,相较于其他三位铣削残余应力,能够极大的提升计算速度,降低资源使用率。本发明可操作性强、精确度高,只需经过简单的分析计算和数据迭代,即可得到接近真实工况的铣削残余应力,为铣削机理研究和表面完整性控制提供可靠数据。
附图说明
图1为本发明一种基于二维仿真分析的铣削残余应力精确获取方法流程图;
图2为本发明的模型示意图。
具体实施方式
下面结合附图以及具体实施例进一步说明本发明。
本发明提供一种基于二维仿真分析的铣削残余应力精确获取方法,其方法流程如下:
步骤(1)、根据实际情况确定所需仿真参数;
所述步骤(1)中需要确定的仿真参数为加工刀具的前角α、前角β、刀具半径r,被加工工件材料的本构模型、密度、热导率、杨氏模量、热容、泊松比,切削的铣削深度a、铣削宽度b。
步骤(2)、建立二维铣削残余应力仿真模型;
所述步骤(2)中的仿真模型分为刀具模型和工件模型,其中刀具模型为三角形,其中一个锐角角度为90-α-β,工件模型是三角形和矩形的组合,上部分为直角三角形,直角边长度分别为b(mm)和c(mm),下部分为矩形,高度为1(mm),宽度为c(mm),其中,
步骤(3)、初次计算铣削仿真残余应力;
所述步骤(3)中按照铣削角度要求摆放刀具,刀具刀尖沿工件三角形与矩形的交界线运动,切削深度深度由b(mm)逐渐减为0(mm)。
步骤(4)、多次迭代计算获取最终铣削残余应力。
所述步骤(4)中需根据需求确定迭代次数N,每次迭代时残余应力从工件矩形区域中选取,第n次迭代选取的点距三角形与矩形区域交界的直角边距离三角形区域另一直角边将提取的残余应力作为初始应力代入下次计算中,直至N次迭代完成。
实施例1
下面举例具体说明本发明方法,但本发明的保护范围不限于下述实例:
步骤(1):刀具前角5°,刀具后角16°,刀具半径为2mm;被加工材料为钛合金,本构模型为:
密度为4.43g/cm3、热导率为7.03W·m-1·℃-1、杨氏模量为110Gpa、热容为505J/kg·℃、泊松比0.34;铣削宽度为0.2mm,铣削深度为0.4mm。
步骤(2):建立刀具模型,其中一个锐角角度为69°;建立三角形和矩形的组合的工件模型,上部分为直角三角形,直角边长度分别为0.2(mm)和1.3(mm),下部分为矩形,高度为1(mm),宽度为1.3(mm)。
步骤(3):按照金属切削原理中前、后角定义摆放刀具,且刀具刀尖沿工件三角形与矩形的交界线运动,切削深度深度由0.2(mm)逐渐减为0(mm)。
步骤(4):假定迭代次数为4,每次迭代时残余应力从工件矩形区域中选取,第n次迭代选取的点距三角形与矩形区域交界的直角边距离三角形区域另一直角边将提取的残余应力作为初始应力代入下次计算中,直至4次迭代完成,数据如表1所示。
表1迭代数据
最终获得铣削残余应力为259.966Mpa。

Claims (1)

1.一种基于二维仿真分析的铣削残余应力精确获取方法,其特征在于,该方法步骤如下:
步骤(1)、根据实际情况确定所需仿真参数,需要确定的仿真参数为加工刀具的前角α、前角β、刀具半径r,被加工工件材料的本构模型、密度、热导率、杨氏模量、热容、泊松比,切削的铣削深度a、铣削宽度b;
步骤(2)、建立二维铣削残余应力仿真模型,需建立仿真模型分为刀具模型和工件模型,其中刀具模型为三角形,其中一个锐角角度为90-α-β,工件模型是三角形和矩形的组合,上部分为直角三角形,直角边长度分别为铣削宽度b(mm)和c(mm),下部分为矩形,高度为1(mm),宽度为c(mm),其中
步骤(3)、初次计算铣削仿真残余应力,需按照铣削角度要求摆放刀具,刀具刀尖沿工件三角形与矩形的交界线运动,切削深度深度由b(mm)逐渐减为0(mm);
步骤(4)、多次迭代计算获取最终铣削残余应力,需根据需求确定迭代次数N,每次迭代时残余应力从工件矩形区域中选取,第n次迭代选取的点距三角形与矩形区域交界的直角边距离三角形区域另一直角边将提取的残余应力作为初始应力代入下次计算中,直至N次迭代完成。
CN201711189158.8A 2017-11-24 2017-11-24 一种基于二维仿真分析的铣削残余应力精确获取方法 Active CN107967384B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711189158.8A CN107967384B (zh) 2017-11-24 2017-11-24 一种基于二维仿真分析的铣削残余应力精确获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711189158.8A CN107967384B (zh) 2017-11-24 2017-11-24 一种基于二维仿真分析的铣削残余应力精确获取方法

Publications (2)

Publication Number Publication Date
CN107967384A true CN107967384A (zh) 2018-04-27
CN107967384B CN107967384B (zh) 2020-05-26

Family

ID=61999773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711189158.8A Active CN107967384B (zh) 2017-11-24 2017-11-24 一种基于二维仿真分析的铣削残余应力精确获取方法

Country Status (1)

Country Link
CN (1) CN107967384B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112304478A (zh) * 2020-07-30 2021-02-02 北京航空航天大学 一种基于蠕变轮廓法的残余应力测试方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101403909A (zh) * 2008-11-12 2009-04-08 山东理工大学 三角网格细分曲面nc刀轨生成方法
CN102622489A (zh) * 2012-03-26 2012-08-01 上海交通大学 基于acis平台的五轴侧铣加工切削力预测方法
US20170096902A1 (en) * 2015-10-01 2017-04-06 Rolls-Royce Plc Method of applying a thermal barrier coating to a metallic article and a thermal barrier coated metallic article

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101403909A (zh) * 2008-11-12 2009-04-08 山东理工大学 三角网格细分曲面nc刀轨生成方法
CN102622489A (zh) * 2012-03-26 2012-08-01 上海交通大学 基于acis平台的五轴侧铣加工切削力预测方法
US20170096902A1 (en) * 2015-10-01 2017-04-06 Rolls-Royce Plc Method of applying a thermal barrier coating to a metallic article and a thermal barrier coated metallic article

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
曾红 等: "铣削加工曲面残余应力有限元分析", 《机床与液压》 *
董兆伟 等: "铣削过程中残余应力仿真分析", 《北京航空航天大学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112304478A (zh) * 2020-07-30 2021-02-02 北京航空航天大学 一种基于蠕变轮廓法的残余应力测试方法
CN112304478B (zh) * 2020-07-30 2021-09-28 北京航空航天大学 一种基于蠕变轮廓法的残余应力测试方法

Also Published As

Publication number Publication date
CN107967384B (zh) 2020-05-26

Similar Documents

Publication Publication Date Title
CN104289748B (zh) 一种大型薄壁蒙皮自适应等壁厚铣削系统及其加工方法
CN101412196B (zh) 基于刀具角度与切削用量参数变化的车削力预测方法
CN105160059A (zh) 一种基于bp和ga的叶片加工切削用量优化选择方法
CN103084639A (zh) 一种基于非均匀余量的增强薄壁叶片工艺刚度的方法
CN109783968B (zh) 基于多工步的金属切削加工过程的三维有限元仿真方法
CN104759942A (zh) 一种薄壁零件的铣削变形在线测量与补偿加工方法
Grigoriev et al. Creep-feed grinding: an overview of kinematics, parameters and effects on process efficiency/Globoko brusenje: pregled kinematike, parametrov in vplivov na procesno ucinkovitost
CN106125666A (zh) 以切削力波动为约束的曲面加工刀具轨迹规划方法
CN105269402A (zh) 一种基于铣削加工的钛合金材料表面粗糙度预测方法
CN104318022A (zh) 一种预测工件表面粗糙度并提高切削效率的方法
CN105069249A (zh) 一种隔框类整体结构件优化切削路径的方法
CN106682281A (zh) 基于最大切削力的铣削瞬时切削力预测方法
CN110321607A (zh) 一种分段式化铣样板制备方法
CN107180131B (zh) 用于确定多点激光冲击强化薄壁件变形曲率半径的方法
CN104227103A (zh) 薄壁件阶梯对称铣削加工方法
CN107967384A (zh) 一种基于二维仿真分析的铣削残余应力精确获取方法
DONG et al. Study on machining deformation of aircraft monolithic component by FEM and experiment
CN103177166B (zh) 一种基于多项式拟合的模锻件毛坯设计方法
CN104014633B (zh) 一种基于有限元分析法的去除棒料芯部缺陷的模拟方法及基于该模拟方法的冲孔方法
CN107944087A (zh) 薄壁件铣削变形的快速实现方法及系统
CN103624761A (zh) 锅炉筒体管孔划线样板及其制作方法
Zhou et al. Influence of cutting and clamping forces on machining distortion of diesel engine connecting rod
CN106777641B (zh) 一种带倒直角的加强板快速建模方法
CN109605120A (zh) 一种提高航空叶片加工良品率的方法
CN206305833U (zh) 一种用于加工片状产品的cnc夹治具

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