CN105414616B - 螺旋铣孔过程中切削力预报和稳定性判别方法 - Google Patents
螺旋铣孔过程中切削力预报和稳定性判别方法 Download PDFInfo
- Publication number
- CN105414616B CN105414616B CN201510875935.9A CN201510875935A CN105414616B CN 105414616 B CN105414616 B CN 105414616B CN 201510875935 A CN201510875935 A CN 201510875935A CN 105414616 B CN105414616 B CN 105414616B
- Authority
- CN
- China
- Prior art keywords
- cutting
- milling
- cutter
- force
- tooth
- 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
Links
- 238000005520 cutting process Methods 0.000 title claims abstract description 132
- 238000003801 milling Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 claims abstract description 35
- 230000008569 process Effects 0.000 claims abstract description 15
- 230000010354 integration Effects 0.000 claims abstract description 12
- 238000002474 experimental method Methods 0.000 claims abstract description 11
- 230000008929 regeneration Effects 0.000 claims abstract description 11
- 238000011069 regeneration method Methods 0.000 claims abstract description 11
- 230000014509 gene expression Effects 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 5
- 238000010008 shearing Methods 0.000 claims description 4
- 238000012417 linear regression Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000007654 immersion Methods 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000003754 machining Methods 0.000 abstract description 12
- 230000001808 coupling effect Effects 0.000 abstract description 6
- 238000010276 construction Methods 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 14
- 238000005516 engineering process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23C—MILLING
- B23C3/00—Milling particular work; Special milling operations; Machines therefor
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23C—MILLING
- B23C2265/00—Details of general geometric configurations
- B23C2265/40—Spiral
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Numerical Control (AREA)
Abstract
本发明提供了一种螺旋铣孔过程中切削力预报和稳定性判定方法,包括如下步骤:依据待加工孔尺寸,选定刀具及加工参数,生成螺旋铣加工刀路;基于刀齿轨迹圆弧假设,依据相邻刀位的刀具几何位置关系,得到刀齿切入切出角及未变形切厚的解析计算表达式;采用三轴单齿铣削实验标定和三轴插铣实验分别标定得到侧刃与底刃切削力系数;将以上所求参数代入二元切削力模型计算得到侧刃与底刃切削力,并将两者求和实现螺旋铣孔过程的瞬态切削力的精确预报;建立了综合考虑铣刀结构模态耦合效应和动态切厚再生效应的三自由度铣削动力学方程,并基于数值积分方法判定加工过程的稳定性。本发明提高了螺旋铣孔过程中切削力预报和稳定性判定的准确性。
Description
技术领域
本发明涉及螺旋铣孔加工技术,具体地,涉及一种螺旋铣孔过程中切削力预报和稳定性判别方法。
背景技术
在航空航天领域,叠层构件装配中制孔数量多、难度大、要求高、工作繁重。传统“钻-扩-铰”的复杂工序流程,加工稳定性差、作业效率低。因此,一种新型的高效高精度制孔技术—螺旋铣孔应运而生。螺旋铣孔过程中刀具自转的同时围绕着孔中心轴线做公转运动并保持轴向进给。这种特殊的运动方式决定了螺旋铣孔的优势:
1、突破了传统钻孔技术中一把刀具加工同一直径孔的限制,实现了单一刀具加工出不同直径的孔和复杂的台阶孔,并可对孔径实施在线补偿。因此,采用螺旋铣孔工艺可以有效减少换刀次数、提高加工效率;
2、螺旋铣孔是断续切削,不仅排屑容易,而且有利于刀具散热,可降低刀具因热量积累而造成的高温磨损现象,可提高刀具的使用寿命。另外,整个铣孔过程可采用微量润滑甚至空冷方式来实现冷却,属绿色加工范畴;
3、螺旋铣孔的轴向力远远小于传统钻孔,使得制孔有良好的出口质量。另外,加工过程中排屑空间大,能够有效防止切屑对已加工表面的损害,有利于提高制孔质量。
然而,螺旋铣孔过程中刀具侧刃和底刃同时参与切削,其加工机理比较复杂。为了在保证加工精度及加工质量的前提下,最大限度地提高材料去除率,这就需要合理选择加工工艺参数。切削力是铣削加工过程中的一个非常重要的物理量,其大小直接影响加工状态和加工表面质量。切削过程的稳定性分析有助于选取合理的加工参数从而有效避免再生颤振的发生,实现加工过程平稳运行。因此,开展螺旋铣孔加工切削力学与动力学建模具有十分重要的意义。
文献“Z.Li,Q.Liu,X.Ming,X.Wang,Y.Dong,Cutting force prediction andanalytical solution of regenerative chatter stability for helical millingoperation,Int.J.Adv.Manuf.Technol.73(2014)433–442.”和文献“C.Liu,G.Wang,M.S.Dargusch,Mechanics and Dynamics of Helical Milling Operations,Vestnik-Journal Mech.Eng.60(2014)716–724.”同时建立了螺旋铣孔切削力预报模型和加工颤振预报模型。在切削力预报方面,两者考虑了底刃切削力,但忽略了底刃切削力系数与侧刃切削力系数的差异,另外,两者的未变形切厚计算均基于经典正弦函数积假设,这直接限制了其切削力预报精度。在加工颤振预报方面,前者直接忽略了z向的颤振,后者虽然考虑了z向的再生效应,但其将x,y向和z向分作两个子系统来考虑,忽略了x,y向和z向的耦合效应。因此,两者稳定性判定的精度皆难以得到保证。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种螺旋铣孔过程中切削力预报和稳定性判定方法。
根据本发明提供的螺旋铣孔过程中切削力预报和稳定性判定方法,包括如下步骤:
步骤1:根据已知的待加工孔直径DB、深度LB、刀具直径Dt、刀齿数N、主轴转速n、切向每齿进给fzt以及轴向螺距ap生成螺旋铣加工刀路;
步骤2:基于刀齿轨迹圆弧假设,跟据相距一个每齿进给的两个刀位处的刀具几何位置关系,得到刀齿切入角、切出角及未变形切厚的解析计算表达式;
步骤3,采用三轴单齿铣削实验标定得到切向r、径向t及轴向a的侧刃铣削力系数,包括剪切力系数KS,qc(q=r,t,a)和犁切力系数KS,qe(q=r,t,a);采用多次不同进给的三轴插铣实验,基于线性回归法拟合得到切向r、径向t及轴向a的底刃铣削力系数,包括剪切力系数KB,qc(q=r,t,a)和犁切力系数KB,qe(q=r,t,a);
步骤4,将步骤2中得到的刀齿的切入角、切出角和未变形切厚的计算值、步骤3中标定的侧刃切削力系数和底刃铣削力系数,代入二元切削力模型计算得到侧刃的切向铣削力、径向铣削力、轴向铣削力和底刃的切向铣削力、径向铣削力、轴向铣削力;然后将侧刃的切向铣削力、径向铣削力、轴向铣削力和底刃的切向铣削力、径向铣削力、轴向铣削力的坐标变换到工件坐标系下的X、Y、Z轴方向,两者求和即得螺旋铣孔过程中瞬态切削力;
步骤5,建立螺旋铣孔过程的三自由度时滞动力学方程,并进行状态空间变换得到状态空间方程;
步骤6,根据第二类Volterra积分方程经典的数值积分算法,利用梯形公式,得到相应离散动态映射,进而构造了单个刀齿切削周期上的状态转移矩阵,根据Floquet理论判定该铣削系统的稳定性。
优选地,所述步骤2中的刀齿切入切出角及未变形切厚计算方法具体如下:
首先,根据步骤1确定螺旋刀路直径Dh和轴向每齿进给率fza:
Dh=DB-Dt (1)
对于刀具侧刃,刀齿的切入角和切出角依据相邻两刀位刀具几何关系可确定如下:
其中:
依据刀齿轨迹圆弧近似假设,刀具侧刃的瞬时未变形切厚可确定如下:
其中,hS,j(t)为第j个刀齿在t时刻的未变形切厚值,为第j个刀齿在t时刻的周向浸入角。
对于刀具底刃,其在螺旋铣孔过程中一直参与切削,则切宽即为刀具半径,其瞬时未变形切厚hB即为轴向每齿进给量fza:
hB=fza (6)
优选地,所述步骤4中计算切削力的具体公式如下:
由于轴向切削深度很小,因此忽略刀具螺旋角的影响;基于二元机械力学模型,作用在刀齿j上的切向、径向和轴向的侧刃力FS,q,j(t)表示为:
式中,是窗函数,用于判断当前刀齿是否参数切削:
aptan(t)为侧刃切削深度,其随刀具回转角度变化,表示为:
其中,为刀齿浸入角;
通过坐标变换,并对每个刀齿上切削力求和,可得t时刻作用在刀具侧刃上的切削力:
其中,
其中,θ0为刀位点初始位置角;
基于二元机械力学模型,作用在刀齿j上的切向、径向和轴向的底刃力FB,q,j(t)表示为:
其中,为底刃第j个刀齿在t时刻未变形切厚值;
通过坐标变换,并对每个刀齿上切削力求和,可得t时刻作用在刀具底刃上的切削力:
其中:
将侧刃力与底刃力求和即得t时刻作用在整个刀具上的总切削力:
F(t)=FS(t)+FB(t) (16)
优选地,所述步骤5中,所述的三自由度动力学方程构建如下:
其中:
其中,M、C、K分别表示刀具的模态质量、阻尼、刚度矩阵;q(t)为刀具模态坐标,且振型系数在刀尖点出归一,即q(t)=[x(t),y(t),z(t)]T,x(t)为刀具x方向位移,y(t)为刀具y方向位移,z(t)为刀具z方向位移;
FD(t)为切厚再生效应引起的动态力,表示如下:
其中,T为时滞量且等于刀齿切削周期,系数矩阵和的表达式如下:
式(17)中的静态力项F(t)不影响其稳定性,故略去该项;同时,另p(t)=Mq(t)+Cq(t)/2;记则式(17)可以转化为如下状态空间形式:
其中:
由于B(t)[x(t)-x(t-T)]是齐次方程的非齐次项,式(31)的解可以表示为:
其中,t0为初始时刻。
优选地,所述步骤6中包含如下子步骤:
步骤6.1,等距离散刀齿通过周期T为m个小时间区段,即T=mτ,其中m为正整数;响应的每个采用时间点为ti=t0+(i-1)τ,其中i=1,...,m+1;τ为时间区段;
根据第二类Volterra积分方程经典的数值积分算法,利用梯形公式,式(34)可表示为:
步骤6.2,变换式(35)可得到下面的离散动态映射:
其中:
其中Bi=B(ti),i=1,...,m+1,I为单位矩阵;
步骤6.3,单个刀齿切削周期上的状态转移矩阵Φ为:
根据Floquet理论,若Φ的所有特征值的模均小于1,则系统是稳定的;若Φ中任一特征值的模大于1,则系统是不稳定的。
与现有技术相比,本发明具有如下的有益效果:
1、本发明基于刀齿轨迹圆弧假设,依据相邻刀齿的刀具几何位置关系得到了刀齿切入切出角及未变形切厚的精确解析计算表达式,建立了同时包含侧刃切削和底刃切削的铣削力模型,提高了切削力预报结果的准确;
2、本发明建立了综合考虑铣刀结构模态耦合效应和动态切厚再生效应的三自由度铣削动力学方程,并基于数值积分方法分析了加工过程的稳定性。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明的流程示意图;
图2为本发明中螺旋铣孔加工过程示意图;
图3为本发明中螺旋铣孔加工过程几何参数提取示意图;
图4为本发明中刀具绕待加工孔中心回转一周过程中的x,y,z方向的瞬态切削力;
图5为本发明中螺旋铣孔加工过程的稳定性图谱。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
本实施例提供了一种螺旋铣孔切削力预报及稳定性判定方法,其流程图如图1所示。首先,依据待加工孔尺寸,选定刀具及加工参数,生成螺旋铣加工刀路;基于刀齿轨迹圆弧假设,依据相邻刀位的刀具几何位置关系,得到刀齿切入切出角及未变形切厚的解析计算表达式;采用三轴单齿铣削实验标定和三轴插铣实验分别标定得到侧刃与底刃切削力系数;然后,将以上所求参数代入二元切削力模型计算得到侧刃与底刃切削力,并将两者求和实现螺旋铣孔过程的瞬态切削力的精确预报;在此基础上,建立了综合考虑铣刀结构模态耦合效应和动态切厚再生效应的三自由度铣削动力学方程,并基于数值积分方法判定加工过程的稳定性。
以下以螺旋铣直径为16mm的孔为例进行具体说明。
步骤1,螺旋铣孔刀路生成:已知待加工孔直径DB=16mm、深度LB=20mm,毛坯材料为AL7075,选取刀具参数:刀具直径Dt=10mm、刀齿数N=2,设定螺旋铣加工参数:主轴转速n=3000rpm、切向每齿进给fzt=0.15mm、轴向螺距ap=1mm、铣削方式为逆铣。利用商用CAM软件生成螺旋铣加工刀路,如图2所示。
步骤2,刀齿切入切出角及未变形切厚计算:依据步骤1的相关参数可以确定螺旋刀路直径Dh和轴向每齿进给率fza:
Dh=DB-Dt (1)
对于刀具侧刃,基于刀齿轨迹圆弧近似假设,依据相距一个每齿进给的两个刀位处的刀具几何位置关系(如图3所示),确定刀齿切入角和切出角如下:
其中:
确定瞬时未变形切厚如下:
对于刀具底刃,其在螺旋铣孔过程中一直参与切削,则切宽即为刀具半径,其瞬时未变形切厚即为轴向每齿进给量:
hB=fza (6)
步骤3,切削力系数标定:采用三轴单齿铣削实验标定得到切向、径向及轴向的侧刃切削力系数KS,rc=168、KS,tc=796、KS,ac=222和犁切力系数KS,re=30.8、KS,te=27.7、KS,ae=1.5;采用一系列不同进给下的三轴插铣实验,基于线性回归法拟合得到切向、径向及轴向的底刃切削力系数KB,rc=123、KB,tc=805、KB,ac=261和犁切力系数KB,re=0.13、KB,te=2.23、KB,ae=1.94。
步骤4,瞬态切削力预报:考虑到螺旋铣孔过程中,轴向切削深度很小,因此忽略刀具螺旋角的影响。将步骤2中得到的刀齿切入切出角信息和未变形切厚计算值、步骤3中标定的铣削力系数,代入二元切削力模型,计算得到作用在刀齿j上的切向、径向和轴向的侧刃力:
式中,是窗函数,用于判断当前刀齿是否参数切削:
aptan(t)为侧刃切削深度,其随刀具回转角度变化,表示为:
通过坐标变换,并对每个刀齿上切削力求和,可得t时刻作用在刀具侧刃上的切削力:
其中:
基于二元机械力学模型,作用在刀齿j上的切向、径向和轴向的底刃力表示为:
通过坐标变换,并对每个刀齿上切削力求和,可得t时刻作用在刀具底刃上的切削力:
其中:
将侧刃力与底刃力求和即得t时刻作用在整个刀具上的总切削力(如图4所示):
F(t)=FS(t)+FB(t) (16)
步骤5,螺旋铣孔过程动力学建模:建立综合考虑铣刀结构模态耦合效应和动态切厚再生效应的三自由度铣削动力学方程:
其中:
M,C,K表示刀具的模态质量、阻尼和刚度矩阵。参考大多数文献的假设,忽略x、y、z向的结构模态耦合效应,即只保留M,C,K矩阵的主项。通过模态试验可得x方向的固有频率f=921.0Hz,阻尼比ξ=0.047,刚度k=1.13N·μm-1,y方向的固有频率f=979.6Hz,阻尼比ξ=0.031,刚度k=1.11N·μm-1,z方向的固有频率f=800.0Hz,阻尼比ξ=0.05,刚度k=12N·μm-1,由此即可计算出M,C,K矩阵的主项。q(t)为刀具模态坐标,且振型系数在刀尖点出归一,即q(t)=[x(t),y(t),z(t)]T。FD(t)为切厚再生效应引起的动态力,表示如下:
其中T为时滞量且等于刀齿切削周期,系数矩阵和的表达式如下:
式(17)中的静态力项F(t)不影响其稳定性,故略去该项。同时,使用类似于广泛运用于哈密顿系统中的变换式,即,另p(t)=Mq(t)+Cq(t)/2。记则式(17)可以转化为如下状态空间形式:
其中:
考虑到B(t)[x(t)-x(t-T)]是齐次方程的非齐次项,式(31)的解可以表示为:
步骤6,加工系统稳定性判定:首先,等距离散刀齿通过周期T为m个小时间区段,即T=mτ(其中m为正整数)。响应的每个采用时间点为ti=t0+(i-1)τ,其中i=1,...,m+1。
根据第二类Volterra积分方程经典的数值积分算法,利用梯形公式,式(34)可表示为:
变换式(35)可得到下面的离散动态映射:
其中:
其中Bi=B(ti),i=1,...,m+1。
然后,单个刀齿切削周期上的状态转移矩阵可表示为:
最后,根据Floquet理论,若Φ的所有特征值的模均小于1,则系统是稳定的;若Φ中任一特征值的模大于1,则系统是不稳定的。据此,绘制出螺旋铣孔过程的稳定性图谱如图5所示。
本发明提供了一种螺旋铣孔过程中切削力预报和稳定性判定方法,包括如下步骤:依据待加工孔尺寸,选定刀具及加工参数,生成螺旋铣加工刀路;基于刀齿轨迹圆弧假设,依据相邻刀位的刀具几何位置关系,得到刀齿切入切出角及未变形切厚的解析计算表达式;采用三轴单齿铣削实验标定和三轴插铣实验分别标定得到侧刃与底刃切削力系数;将以上所求参数代入二元切削力模型计算得到侧刃与底刃切削力,并将两者求和实现螺旋铣孔过程的瞬态切削力的精确预报;建立了综合考虑铣刀结构模态耦合效应和动态切厚再生效应的三自由度铣削动力学方程,并基于数值积分方法判定加工过程的稳定性。本发明提高了螺旋铣孔过程中切削力预报和稳定性判定的准确性。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (5)
1.一种螺旋铣孔过程中切削力预报和稳定性判定方法,其特征在于,包括如下步骤:
步骤1:根据已知的待加工孔直径DB、深度LB、刀具直径Dt、刀齿数N、主轴转速n、切向每齿进给fzt以及轴向螺距ap生成螺旋铣加工刀路;
步骤2:基于刀齿轨迹圆弧假设,根据相邻刀位的刀具几何位置关系,得到刀齿切入角、切出角及未变形切厚的解析计算表达式;
步骤3,采用三轴单齿铣削实验标定得到切向r、径向t及轴向a的侧刃铣削力系数,包括剪切力系数KS,qc(q=r,t,a)和犁切力系数KS,qe(q=r,t,a);采用多次不同进给的三轴插铣实验,基于线性回归法拟合得到切向r、径向t及轴向a的底刃铣削力系数,包括剪切力系数KB,qc(q=r,t,a)和犁切力系数KB,qe(q=r,t,a);
步骤4,将步骤2中得到的刀齿的切入角、切出角和未变形切厚的计算值、步骤3中标定的侧刃切削力系数和底刃铣削力系数,代入二元切削力模型计算得到侧刃的切向铣削力、径向铣削力、轴向铣削力和底刃的切向铣削力、径向铣削力、轴向铣削力;然后将侧刃的切向铣削力、径向铣削力、轴向铣削力和底刃的切向铣削力、径向铣削力、轴向铣削力的坐标变换到工件坐标系下的X、Y、Z轴方向,两者求和即得螺旋铣孔过程中瞬态切削力;
步骤5,建立螺旋铣孔过程的三自由度时滞动力学方程,并进行状态空间变换得到状态空间方程;
步骤6,根据第二类Volterra积分方程经典的数值积分算法,利用梯形公式,得到相应离散动态映射,进而构造了单个刀齿切削周期上的状态转移矩阵,根据Floquet理论判定该铣削系统的稳定性。
2.根据权利要求1所述的螺旋铣孔过程中切削力预报和稳定性判定方法,其特征在于,所述步骤2中的刀齿切入切出角及未变形切厚计算方法具体如下:
首先,根据步骤1确定螺旋刀路直径Dh和轴向每齿进给率fza:
Dh=DB-Dt (1)
对于刀具侧刃,刀齿的切入角和切出角依据相邻两刀位刀具几何关系可确定如下:
其中:
依据刀齿轨迹圆弧近似假设,刀具侧刃的瞬时未变形切厚可确定如下:
其中,hS,j(t)为第j个刀齿在t时刻的未变形切厚值,为第j个刀齿在t时刻的周向浸入角;
对于刀具底刃,其在螺旋铣孔过程中一直参与切削,则切宽即为刀具半径,其瞬时未变形切厚hB即为轴向每齿进给量fza:
hB=fza (6) 。
3.根据权利要求2所述的螺旋铣孔过程中切削力预报和稳定性判定方法,其特征在于,所述步骤4中计算切削力的具体公式如下:
由于轴向切削深度很小,因此忽略刀具螺旋角的影响;基于二元机械力学模型,作用在刀齿j上的切向、径向和轴向的侧刃力FS,q,j(t)表示为:
式中,是窗函数,用于判断当前刀齿是否参数切削:
aptan(t)为侧刃切削深度,其随刀具回转角度变化,表示为:
其中,为刀齿浸入角;
通过坐标变换,并对每个刀齿上切削力求和,可得t时刻作用在刀具侧刃上的切削力:
其中,
其中,θ0为刀位点初始位置角;
基于二元机械力学模型,作用在刀齿j上的切向、径向和轴向的底刃力FB,q,j(t)表示为:
其中,为底刃第j个刀齿在t时刻未变形切厚值;
通过坐标变换,并对每个刀齿上切削力求和,可得t时刻作用在刀具底刃上的切削力:
其中:
将侧刃力与底刃力求和即得t时刻作用在整个刀具上的总切削力:
F(t)=FS(t)+FB(t) (16) 。
4.根据权利要求3所述的螺旋铣孔过程中切削力预报和稳定性判定方法,其特征在于,所述步骤5中,所述的三自由度动力学方程构建如下:
其中:
其中,M、C、K分别表示刀具的模态质量、阻尼、刚度矩阵;q(t)为刀具模态坐标,且振型系数在刀尖点处归一,即q(t)=[x(t),y(t),z(t)]T,x(t)为刀具x方向位移,y(t)为刀具y方向位移,z(t)为刀具z方向位移;
FD(t)为切厚再生效应引起的动态力,表示如下:
其中,T为时滞量且等于刀齿切削周期,系数矩阵和的表达式如下:
式(17)中的静态力项F(t)不影响其稳定性,故略去该项;同时,另p(t)=Mq(t)+Cq(t)/2;记则式(17)可以转化为如下状态空间形式:
其中:
由于B(t)[x(t)-x(t-T)]是齐次方程的非齐次项,式(31)的解可以表示为:
其中,t0为初始时刻。
5.根据权利要求4所述的螺旋铣孔过程中切削力预报和稳定性判定方法,其特征在于,所述步骤6中包含如下子步骤:
步骤6.1,等距离散刀齿通过周期T为m个小时间区段,即T=mτ,其中m为正整数;响应的每个采用时间点为ti=t0+(i-1)τ,其中i=1,...,m+1;τ为时间区段;
根据第二类Volterra积分方程经典的数值积分算法,利用梯形公式,式(34)可表示为:
步骤6.2,变换式(35)可得到下面的离散动态映射:
其中:
其中Bi=B(ti),i=1,...,m+1,I为单位矩阵;
步骤6.3,单个刀齿切削周期上的状态转移矩阵Φ为:
根据Floquet理论,若Φ的所有特征值的模均小于1,则系统是稳定的;若Φ中任 一特征值的模大于1,则系统是不稳定的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510875935.9A CN105414616B (zh) | 2015-12-02 | 2015-12-02 | 螺旋铣孔过程中切削力预报和稳定性判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510875935.9A CN105414616B (zh) | 2015-12-02 | 2015-12-02 | 螺旋铣孔过程中切削力预报和稳定性判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105414616A CN105414616A (zh) | 2016-03-23 |
CN105414616B true CN105414616B (zh) | 2017-09-15 |
Family
ID=55493396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510875935.9A Active CN105414616B (zh) | 2015-12-02 | 2015-12-02 | 螺旋铣孔过程中切削力预报和稳定性判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105414616B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106502202B (zh) * | 2017-01-06 | 2018-10-16 | 大连理工大学 | 一种球头铣刀与导向叶片接触区域的半解析建模方法 |
CN106956026A (zh) * | 2017-03-15 | 2017-07-18 | 广东长盈精密技术有限公司 | 螺旋铣圆方法 |
CN109708875B (zh) * | 2019-01-24 | 2020-11-10 | 北华大学 | 一种基于变分模态分解Volterra模型奇异值熵的转子故障诊断方法 |
CN109909806B (zh) * | 2019-03-22 | 2020-06-19 | 南京理工大学 | 一种提升机器人铣削稳定域的方法 |
CN111390253B (zh) * | 2020-04-26 | 2021-02-09 | 兰州理工大学 | 一种叠层材料制孔加工轨迹规划方法 |
CN112091292B (zh) * | 2020-09-14 | 2022-09-02 | 中国航发贵州黎阳航空动力有限公司 | 一种余量孔扩孔加工方法 |
CN112733298B (zh) * | 2021-01-12 | 2022-04-29 | 天津大学 | 基于螺旋铣孔的混联机器人不同位姿下加工性能评价方法 |
CN114043312B (zh) * | 2021-11-26 | 2022-10-11 | 广东海洋大学 | 一种螺旋铣孔状态异常的在线检测识别方法及系统 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4434137B2 (ja) * | 2005-12-20 | 2010-03-17 | 株式会社日立製作所 | 傾斜面の切削加工における切削自励振動の安定限界軸方向切込み量の算出方法 |
CN101653841B (zh) * | 2009-08-27 | 2010-12-29 | 西北工业大学 | 铣削过程稳定域判定方法 |
CN101853324B (zh) * | 2010-06-24 | 2012-07-04 | 西北工业大学 | 圆周铣过程中铣削力建模方法 |
CN102059381B (zh) * | 2010-11-04 | 2012-08-01 | 西北工业大学 | 周铣过程中加工误差预测方法 |
CN103559550B (zh) * | 2013-09-09 | 2016-10-05 | 西北工业大学 | 多模态耦合下的铣削稳定域预测方法 |
CN103955169B (zh) * | 2014-04-11 | 2017-01-18 | 上海交通大学 | 五轴数控侧铣加工铣削力预报方法 |
CN104657607B (zh) * | 2015-02-10 | 2017-12-15 | 北京理工大学 | 一种薄壁件铣削稳定性预测方法 |
CN104778333B (zh) * | 2015-04-27 | 2017-09-12 | 华中科技大学 | 一种螺旋立铣刀正交车铣加工三维稳定性建模方法 |
-
2015
- 2015-12-02 CN CN201510875935.9A patent/CN105414616B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105414616A (zh) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105414616B (zh) | 螺旋铣孔过程中切削力预报和稳定性判别方法 | |
CN103646141B (zh) | 一种平底螺旋立铣刀正交车铣轴类零件的切削力建模方法 | |
Engin et al. | Mechanics and dynamics of general milling cutters.: Part II: inserted cutters | |
Shi et al. | An improved cutting power-based model for evaluating total energy consumption in general end milling process | |
CN106156477B (zh) | 薄壁件动态铣削稳定性叶瓣图高精度预测方法 | |
Han et al. | Precise prediction of forces in milling circular corners | |
JP2011206894A (ja) | 工具軌跡生成装置、工具軌跡算出方法および工具軌跡生成プログラム | |
de Oliveira et al. | Evaluating the influences of the cutting parameters on the surface roughness and form errors in 4-axis milling of thin-walled free-form parts of AISI H13 steel | |
CN108563848B (zh) | 一种平底螺旋立铣刀的铣削力建模方法 | |
CN113962105B (zh) | 一种无颤振精加工铣削过程的高效参数优化方法 | |
CN108804749B (zh) | 薄壁件铣削过程阻尼建模方法 | |
Taner et al. | Generalized cutting force model in multi-axis milling using a new engagement boundary determination approach | |
Parsian et al. | A mechanistic approach to model cutting forces in drilling with indexable inserts | |
CN114186175A (zh) | 振动作用下高能效铣刀主切削力能耗动态特性的解算方法 | |
CN111045384B (zh) | 数控加工方法、加工装置与数控机床 | |
CN113486475A (zh) | 一种圆柱齿轮滚齿加工切削力的预测方法 | |
CN111176209B (zh) | 型腔螺旋铣削加工进给率与转速离线规划方法 | |
Suzuki et al. | Force prediction in cutting operations with self-propelled rotary tools considering bearing friction | |
Zhang et al. | Study of tool trajectory in blisk channel ECM with spiral feeding | |
CN110727246B (zh) | 一种基于刀位文件的刀具与工件瞬时接触轮廓提取方法 | |
CN111299668B (zh) | 一种不等齿距铣刀的齿间角确定方法 | |
CN111112755B (zh) | 一种变传动比齿条插齿力的预测方法 | |
CN114599468B (zh) | 具有螺旋角过渡的铣削刀具 | |
Chen et al. | FEM-based prediction and control of milling deformation for a thin-wall web of Ti-6Al-4V alloy | |
Rekers et al. | Cutting force reduction in the milling of aluminum alloys with serrated cutting tool edges |
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 |