CN106940746A - 基于薄壁件的铣削颤振稳定性预测的并行时域方法 - Google Patents
基于薄壁件的铣削颤振稳定性预测的并行时域方法 Download PDFInfo
- Publication number
- CN106940746A CN106940746A CN201710150868.3A CN201710150868A CN106940746A CN 106940746 A CN106940746 A CN 106940746A CN 201710150868 A CN201710150868 A CN 201710150868A CN 106940746 A CN106940746 A CN 106940746A
- Authority
- CN
- China
- Prior art keywords
- thin
- wall part
- time domain
- milling
- cutting
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
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)
- Numerical Control (AREA)
Abstract
本发明公开了一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,具有如下步骤:获取刀具和薄壁件的模态参数;建立薄壁件铣削过程在时域上的运动微分方程并将其在一个刀齿周期内离散化;建立运动微分方程在一个刀齿周期内离散化的数学模型;得出不同铣刀主轴转速下的薄壁件铣削的颤振稳定性临界轴向切深;绘制颤振稳定性叶瓣图和三维颤振稳定性叶瓣图。本发明考虑了铣削过程在一个刀齿周期内的时变性,在整个加工过程中,薄壁件的各阶固有频率、刚度、阻尼比的时变性,刀具与薄壁件的模态,因此,本发明能更加准确、真实地预测薄壁件铣削的颤振;由于结合了并行计算理论,本发明的计算时间随着线程数的增加而成倍缩短,因此,本发明具有高效性。
Description
技术领域
本发明属于薄壁件铣削加工稳定性预测技术领域,具体涉及基于薄壁件的铣削颤振稳定性预测的并行时域方法。
背景技术
在影响铣削加工表面质量、拖慢铣削加工速度的众多因素中,颤振是最重要的一个因素,与普通工件不同,薄壁件刚度极差,更容易发生颤振,因此,颤振稳定性预测理论方法对于薄壁件铣削加工有重要意义。
铣削过程在一个刀齿周期内是不断变化着的,这种在刀齿周期内的时变性需要被考虑,在整个加工过程中,随着薄壁工件的材料去除与形状变化,工件的各阶固有频率、刚度、阻尼比也在变化,这种工件模态的时变性也要被考虑,薄壁件的刚度通常比刀具要小,因此需要同时考虑工件与刀具的模态。
目前,国内在研究薄壁件铣削稳定性方面对于铣削过程在刀齿周期内的时变性考虑的不够充分,这极大地影响了薄壁件的加工质量与加工精度。在研究中没有建立刀齿周期内的时变性模型与相对传递函数模型,也没有考虑工件传递函数随着薄壁工件的材料去除与形状变化。因此。目前的预测方法并不能准确地预测薄壁件铣削加工颤振稳定性。
发明内容
根据上述提出的技术问题,而提供一种基于薄壁件的铣削颤振稳定性预测的并行时域方法。本发明采用的技术手段如下:
一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,具有如下步骤:
S1、把加工过程划分成几个阶段,通过模态实验获取在各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数;
S2、建立薄壁件铣削过程在时域上的运动微分方程并将其在一个刀齿周期内离散化;
S3、通过各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数,建立薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型;
S4、根据并行计算理论,求解薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型,得出不同铣刀主轴转速下的薄壁件铣削的颤振稳定性临界轴向切深;
S5、根据一个加工阶段的主轴转速与颤振稳定性临界轴向切深绘制颤振稳定性叶瓣图,根据所有加工阶段的颤振稳定性叶瓣图绘制三维颤振稳定性叶瓣图。
以主轴转速为横轴,以颤振稳定性临界轴向切深为纵轴,可以绘制出某一个加工阶段的颤振稳定性叶瓣图,再以加工阶段编号为斜轴,可以绘制出表征整个加工过程的时变性的三维颤振稳定性叶瓣图。
所述步骤S1包括以下步骤:
S11、把加工过程划分成几个阶段;
S12、对各个阶段的不同薄壁件形状进行模态实验,
把刀具安装到刀柄上。把测力仪安装到机床平台上,因为后续实验要安装测力仪,所以为了保持模态数据的一致性,测量模态时也要装上测力仪。把虎钳安装在测力仪上,用虎钳夹紧薄壁件。用3M双面胶把加速度传感器粘贴到待测物体(刀具或薄壁件)上(沿X方向或Y方向),连接激振力锤、采集卡、传感器与计算机。开启采集软件CutPRO。用激振力锤沿Y方向敲击待测物体,保存采集结果数据。把传感器粘到另一个方向,再次用激振力锤沿Y方向敲击待测物体,保存采集结果数据。结束;
S13、根据模态实验数据计算各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数。
所述步骤S2包括以下步骤:
S21、建立薄壁件铣削过程在时域上的运动微分方程,
单模态系统在时域上的运动微分方程可以由式(1)表示:
S22、将薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化,把一个刀齿周期T分成m小段,那么第i小段时间就可以近似用式(2)表示:
也可以将式(2)表示为式(3):
其中:
在求解时,[q(t-T)]会被[qi-m]与[qi-m+1]以插值的方式近似求出,考虑到实际计算时离散化的时间间隔会非常小,所以把代入由拉格朗日一次插值形成的插值函数即可满足精度。于是式(3)就变成了式(4):
所述步骤S3包括以下步骤:
将t=ti+1带入式(4),得到式(5):
其中:
构造(2m+4)维向量zi:
构造(2m+4)×(2m+4)阶矩阵Di:
于是有:
简记为:
zi+1=[Di]×zi (8)
那么对整个一个刀齿周期:
zi+m=[Φ]×zi=[Dm-1]×[Dm-2]×[Dm-3]×[Dm-4]…[D1]×[D0]×zi (9)
根据Floquet定理,只有当[Φ]的所有特征值的模都小于1时,系统才是稳定的。
接下来,根据刀具的x,y方向的各阶模态ωni,m(xx),ki,m(xx),ξi,m(xx),ωni,m(yy),ki,m(yy),ξi,m(yy),与薄壁件的y方向的各阶模态ωni,w(yy),ki,w(yy),ξi,w(yy),
于是就要把式(8)扩展成下式:
式中和中的k均代表系统的某一阶模态。
然后式(9)就变成了下式:
于是就有:
式(11)为薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型,[Φ]为传递矩阵;
此时Floquet定理仍然是适用的,因为系统稳定就代表每一阶模态都稳定,即对于任意的k,的所有特征值的模都小于1。而式(12)的所有特征值的集合显然就是每个 的所有特征值的并集,所以仍然有以下结论:只有当[Φ]的所有特征值的模都小于1时,系统才是稳定的。所述并行计算理论为把扫描主轴转速的过程分片,每一个分片线程并行独立计算,得到转速——颤振稳定性临界轴向切深切深数据元组,最后同步线程,将每一片的转速——临界轴向切深切深数据元组组合起来,形成全部的转速——颤振稳定性临界轴向切深切深数据元组
所述每一个分片线程的计算步骤如下:
S41、初始化主轴转速,令主轴转速n等于初始转速值nbegin;
S42、判断总体终止条件,如果主轴转速n等于终止转速值nend,计算结束,否则,执行步骤S43;
S43、初始化扫描轴向切深,令扫描轴向切深ai等于初始扫描切深值abegin;
S44、判断轴向切深扫描终止条件,如果扫描轴向切深ai等于终止扫描切深值aend,执行步骤S45,否则,执行步骤S46;
S45、增加主轴转速,把主轴转速n加上一个转速步长Δn,然后,执行步骤S42;
Δn的大小根据具体情况而定,Δn越小,叶瓣图就越细腻,比如转速范围是10到20000rpm,那么Δn=100rpm就很细腻了。如果转速范围是10到2000rpm,那么Δn=10rpm就很细腻了。
S46、针对当前的扫描轴向切深ai,计算传递矩阵;
S47、计算传递矩阵的所有特征值λ1,λ2,λ3,…,选择其中模最大λi的模的值为上界值λm;
S48、如果λm等于1,执行步骤S49,否则,执行步骤S410;
S49、把步骤S46中当前的扫描轴向切深ai记录为主轴转速n下的颤振稳定性临界轴向切深,执行步骤S410;
S410、增加扫描轴向切深,把扫描轴向切深ai加上一个切深步长Δai,然后,执行步骤S44。
本发明具有以下优点:
本发明提出基于薄壁件的铣削颤振稳定性预测的并行时域方法,考虑了铣削过程在一个刀齿周期内的时变性,也考虑了在整个加工过程中,薄壁件的各阶固有频率、刚度、阻尼比的时变性,并且同时考虑刀具与薄壁件的模态,因此,本发明能更加准确、真实地预测薄壁件铣削的颤振;表1比较了不同线程数量下基于薄壁件的铣削颤振稳定性预测的并行时域方法的计算时间,可以看出由于结合了并行计算理论,本发明的计算时间随着线程数的增加而成倍缩短,因此,本发明具有高效性。
表1
计算时间(s) | |
单线程 | 48.05 |
2线程 | 27.36 |
3线程 | 23.49 |
基于上述理由本发明可在薄壁件铣削加工稳定性预测等技术领域广泛推广。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明的具体实施方式中基于薄壁件的铣削颤振稳定性预测的并行时域方法的流程图。
图2是本发明的具体实施方式中加工过程中各个阶段的薄壁件形状。
图3是本发明的具体实施方式中基于薄壁件的铣削颤振稳定性预测的并行时域方法的动力学模型的空间结构示意图。
图4是本发明的具体实施方式中基于薄壁件的铣削颤振稳定性预测的并行时域方法的动力学模型的俯视图。
图5是本发明的具体实施方式中每一个分片线程的计算步骤的流程图。
图6是本发明的具体实施方式中预测整个铣削加工过程的颤振稳定性叶瓣图。
具体实施方式
如图1-图6所示,一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,
具有如下步骤:
S1、把加工过程划分成几个阶段,通过模态实验获取在各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数:
S11、把加工过程划分成3个阶段,0阶段,1阶段,2阶段,每个阶段的薄壁件的形状如图2所示;
S12、对各个阶段的不同薄壁件形状进行模态实验,
把刀具安装到刀柄上。把测力仪安装到机床平台上,因为后续实验要安装测力仪,所以为了保持模态数据的一致性,测量模态时也要装上测力仪。把虎钳安装在测力仪上,用虎钳夹紧薄壁件。用3M双面胶把加速度传感器粘贴到待测物体(刀具或薄壁件)上(沿X方向或Y方向),连接激振力锤、采集卡、传感器与计算机。开启采集软件CutPRO。用激振力锤沿Y方向敲击待测物体,保存采集结果数据。把传感器粘到另一个方向,再次用激振力锤沿Y方向敲击待测物体,保存采集结果数据。结束;
S13、根据模态实验数据计算各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数,
刀具的模态参数如表2所示,
表2
在0阶段薄壁件的模态参数如表3示,
表3
一阶固有频率(Hz) | 343.5427 |
一阶刚度(N/m) | 2.1061e6 |
一阶阻尼比 | 4.8029e-2 |
在1阶段薄壁件的模态参数如表4示,
表4
一阶固有频率(Hz) | 357.2227 |
一阶刚度(N/m) | 6.83111e5 |
一阶阻尼比 | 4.8989e-2 |
在2阶段薄壁件的模态参数如表5示,
表5
一阶固有频率(Hz) | 408.725 |
一阶刚度(N/m) | 3.87417e6 |
一阶阻尼比 | 6.117e-2 |
二阶固有频率(Hz) | 504 |
二阶刚度(N/m) | 1.02448e7 |
二阶阻尼比 | 2.87698e-2 |
三阶固有频率(Hz) | 2292.1365 |
三阶刚度(N/m) | 4.6962e7 |
三阶阻尼比 | 1.46152e-2 |
四阶固有频率(Hz) | 4623.6416 |
四阶刚度(N/m) | 4.45815e8 |
四阶阻尼比 | 1.7302e-3 |
S2、建立薄壁件铣削过程在时域上的运动微分方程并将其在一个刀齿周期内离散化:
S21、建立薄壁件铣削过程在时域上的运动微分方程,
单模态系统在时域上的运动微分方程可以由式(1)表示:
S22、将薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化,把一个刀齿周期T分成m小段,那么第i小段时间就可以近似用式(2)表示:
也可以将式(2)表示为式(3):
其中:
在求解时,[q(t-T)]会被[qi-m]与[qi-m+1]以插值的方式近似求出,考虑到实际计算时离散化的时间间隔会非常小,所以把代入由拉格朗日一次插值形成的插值函数即可满足精度。于是式(3)就变成了式(4):
S3、通过各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数,建立薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型:
将t=ti+1带入式(4),得到式(5):
其中:
构造(2m+4)维向量zi:
构造(2m+4)×(2m+4)阶矩阵Di:
于是有:
简记为:
zi+1=[Di]×Zi (8)
那么对整个一个刀齿周期:
Zi+m=[Φ]×zi=[Dm-1]×[Dm-2]×[Dm-3]×[Dm-4]…[D1]×[D0]×zi(9)
根据Floquet定理,只有当[Φ]的所有特征值的模都小于1时,系统才是稳定的。
接下来,根据刀具的x,y方向的各阶模态ωni,m(xx),ki,m(xx),ξi,m(xx),ωni,m(yy),ki,m(yy),ξi,m(yy),与薄壁件的y方向的各阶模态ωni,w(yy),ki,w(yy),ξi,w(yy),
于是就要把式(8)扩展成下式:
式中和中的k均代表系统的某一阶模态。
然后式(9)就变成了下式:
于是就有:
式(11)为薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型,[Φ]为传递矩阵;
此时Floquet定理仍然是适用的,因为系统稳定就代表每一阶模态都稳定,即对于任意的k,的所有特征值的模都小于1。而式(12)的所有特征值的集合显然就是每个 的所有特征值的并集,所以仍然有以下结论:只有当|Φ]的所有特征值的模都小于1时,系统才是稳定的;
S4、根据并行计算理论,求解薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型,得出不同铣刀主轴转速下的薄壁件铣削的颤振稳定性临界轴向切深,
所述并行计算理论为把扫描主轴转速的过程分片,每一个分片线程并行独立计算,得到转速——颤振稳定性临界轴向切深切深数据元组,最后同步线程,将每一片的转速——临界轴向切深切深数据元组组合起来,形成全部的转速——颤振稳定性临界轴向切深切深数据元组。
所述每一个分片线程的计算步骤如下:
S41、初始化主轴转速,令主轴转速n等于初始转速值nbegin;
S42、判断总体终止条件,如果主轴转速n等于终止转速值nend,计算结束,否则,执行步骤S43;
S43、初始化扫描轴向切深,令扫描轴向切深ai等于初始扫描切深值abegin;
S44、判断轴向切深扫描终止条件,如果扫描轴向切深ai等于终止扫描切深值aend,执行步骤S45,否则,执行步骤S46;
S45、增加主轴转速,把主轴转速n加上一个转速步长Δn,然后,执行步骤S42;
S46、针对当前的扫描轴向切深ai,计算传递矩阵;
S47、计算传递矩阵的所有特征值λ1,λ2,λ3,…,选择其中模最大λi的模的值为上界值λm;
S48、如果λm等于1,执行步骤S49,否则,执行步骤S410;
S49、把步骤S46中当前的扫描轴向切深ai记录为主轴转速n下的颤振稳定性临界轴向切深,执行步骤S410;
S410、增加扫描轴向切深,把扫描轴向切深ai加上一个切深步长Δai,然后,执行步骤S44。;
S5、根据一个加工阶段的主轴转速与颤振稳定性临界轴向切深绘制颤振稳定性叶瓣图,根据所有加工阶段的颤振稳定性叶瓣图绘制三维颤振稳定性叶瓣图。
以主轴转速为横轴,以颤振稳定性临界轴向切深为纵轴,可以绘制出某一个加工阶段的颤振稳定性叶瓣图,再以加工阶段编号为斜轴,可以绘制出表征整个加工过程的时变性的三维颤振稳定性叶瓣图。
在三维颤振稳定性叶瓣图中,曲面之下的点(某一加工阶段,某一铣刀主轴转速,某一铣刀轴向切深)所代表的工况就不会造成颤振,曲面之上的点所代表的工况就会造成颤。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (5)
1.一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,其特征在于具有如下步骤:
S1、把加工过程划分成几个阶段,通过模态实验获取在各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数;
S2、建立薄壁件铣削过程在时域上的运动微分方程并将其在一个刀齿周期内离散化;
S3、通过各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数,建立薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型;
S4、根据并行计算理论,求解薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化的数学模型,得出不同铣刀主轴转速下的薄壁件铣削的颤振稳定性临界轴向切深;
S5、根据一个加工阶段的主轴转速与颤振稳定性临界轴向切深绘制颤振稳定性叶瓣图,根据所有加工阶段的颤振稳定性叶瓣图绘制三维颤振稳定性叶瓣图。
2.根据权利要求1所述的一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,其特征在于:所述步骤S1包括以下步骤:
S11、把加工过程划分成几个阶段;
S12、对各个阶段的不同薄壁件形状进行模态实验;
S13、根据模态实验数据计算各个阶段的开始时刻的刀具的模态参数与薄壁件的模态参数。
3.根据权利要求1所述的一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,其特征在于:所述步骤S2包括以下步骤:
S21、建立薄壁件铣削过程在时域上的运动微分方程;
S22、将薄壁件铣削过程在时域上的运动微分方程在一个刀齿周期内离散化。
4.根据权利要求1所述的一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,其特征在于:所述并行计算理论为把扫描主轴转速的过程分片,每一个分片线程并行独立计算,得到转速——颤振稳定性临界轴向切深切深数据元组,最后同步线程,将每一片的转速——临界轴向切深切深数据元组组合起来,形成全部的转速——颤振稳定性临界轴向切深切深数据元组。
5.根据权利要求4所述的一种基于薄壁件的铣削颤振稳定性预测的并行时域方法,其特征在于:所述每一个分片线程的计算步骤如下:
S41、初始化主轴转速,令主轴转速n等于初始转速值nbegin;
S42、判断总体终止条件,如果主轴转速n等于终止转速值nend,计算结束,否则,执行步骤S43;
S43、初始化扫描轴向切深,令扫描轴向切深ai等于初始扫描切深值abegin;
S44、判断轴向切深扫描终止条件,如果扫描轴向切深ai等于终止扫描切深值aend,执行步骤S45,否则,执行步骤S46;
S45、增加主轴转速,把主轴转速n加上一个转速步长Δn,然后,执行步骤S42;
S46、针对当前的扫描轴向切深ai,计算传递矩阵;
S47、计算传递矩阵的所有特征值λ1,λ2,λ3,…,选择其中模最大λi的模的值为上界值λm;
S48、如果λm等于1,执行步骤S49,否则,执行步骤S410;
S49、把步骤S46中当前的扫描轴向切深ai记录为主轴转速n下的颤振稳定性临界轴向切深,执行步骤S410;
S410、增加扫描轴向切深,把扫描轴向切深ai加上一个切深步长Δai,然后,执行步骤S44。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710150868.3A CN106940746A (zh) | 2017-03-14 | 2017-03-14 | 基于薄壁件的铣削颤振稳定性预测的并行时域方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710150868.3A CN106940746A (zh) | 2017-03-14 | 2017-03-14 | 基于薄壁件的铣削颤振稳定性预测的并行时域方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106940746A true CN106940746A (zh) | 2017-07-11 |
Family
ID=59469874
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710150868.3A Pending CN106940746A (zh) | 2017-03-14 | 2017-03-14 | 基于薄壁件的铣削颤振稳定性预测的并行时域方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106940746A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107807526A (zh) * | 2017-10-31 | 2018-03-16 | 上海交通大学 | 一种基于稳定性仿真智能抑制加工颤振的方法 |
CN108519758A (zh) * | 2018-03-21 | 2018-09-11 | 南京工业大学 | 一种避免干涉和颤振的多轴数控加工球头刀轴矢量规划方法 |
CN108846242A (zh) * | 2018-07-09 | 2018-11-20 | 西北工业大学 | 基于预拉应力施加的薄壁件铣削颤振抑制方法 |
CN108984891A (zh) * | 2018-07-09 | 2018-12-11 | 西北工业大学 | 基于预应力施加的薄壁件铣削稳定性改善方法 |
CN109093130A (zh) * | 2018-08-08 | 2018-12-28 | 西北工业大学 | 基于分区域加工的薄壁环形零件粗加工工艺优化方法 |
CN112380726A (zh) * | 2020-11-26 | 2021-02-19 | 华中科技大学无锡研究院 | 基于模态耦合颤振的机器人铣削临界稳定切深预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657607A (zh) * | 2015-02-10 | 2015-05-27 | 北京理工大学 | 一种薄壁件支撑装置及铣削稳定性预测方法 |
CN105447286A (zh) * | 2016-01-28 | 2016-03-30 | 东北大学 | 一种预测叶轮颤振的方法 |
CN105608288A (zh) * | 2016-01-20 | 2016-05-25 | 东北大学 | 一种基于过阻尼效应铣削颤振稳定性预测方法 |
CN106156477A (zh) * | 2015-04-28 | 2016-11-23 | 河南理工大学 | 薄壁件动态铣削稳定性叶瓣图高精度预测方法 |
-
2017
- 2017-03-14 CN CN201710150868.3A patent/CN106940746A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657607A (zh) * | 2015-02-10 | 2015-05-27 | 北京理工大学 | 一种薄壁件支撑装置及铣削稳定性预测方法 |
CN106156477A (zh) * | 2015-04-28 | 2016-11-23 | 河南理工大学 | 薄壁件动态铣削稳定性叶瓣图高精度预测方法 |
CN105608288A (zh) * | 2016-01-20 | 2016-05-25 | 东北大学 | 一种基于过阻尼效应铣削颤振稳定性预测方法 |
CN105447286A (zh) * | 2016-01-28 | 2016-03-30 | 东北大学 | 一种预测叶轮颤振的方法 |
Non-Patent Citations (3)
Title |
---|
DING Y 等: "A full-discretization method for prediction of milling stability", 《INTERNATIONAL JOURNAL OF MACHINE TOOLS & MANUFACTURE》 * |
张雪薇 等: "薄壁零件铣削三维颤振稳定性建模与分析", 《东北大学学报( 自然科学版)》 * |
徐帅: "面向钛合金薄壁件铣削颤振稳定性及加工误差研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107807526A (zh) * | 2017-10-31 | 2018-03-16 | 上海交通大学 | 一种基于稳定性仿真智能抑制加工颤振的方法 |
CN108519758A (zh) * | 2018-03-21 | 2018-09-11 | 南京工业大学 | 一种避免干涉和颤振的多轴数控加工球头刀轴矢量规划方法 |
CN108519758B (zh) * | 2018-03-21 | 2020-12-22 | 南京工业大学 | 一种避免干涉和颤振的多轴数控加工球头刀轴矢量规划方法 |
CN108846242A (zh) * | 2018-07-09 | 2018-11-20 | 西北工业大学 | 基于预拉应力施加的薄壁件铣削颤振抑制方法 |
CN108984891A (zh) * | 2018-07-09 | 2018-12-11 | 西北工业大学 | 基于预应力施加的薄壁件铣削稳定性改善方法 |
CN108984891B (zh) * | 2018-07-09 | 2019-04-16 | 西北工业大学 | 基于预应力施加的薄壁件铣削稳定性改善方法 |
CN109093130A (zh) * | 2018-08-08 | 2018-12-28 | 西北工业大学 | 基于分区域加工的薄壁环形零件粗加工工艺优化方法 |
CN112380726A (zh) * | 2020-11-26 | 2021-02-19 | 华中科技大学无锡研究院 | 基于模态耦合颤振的机器人铣削临界稳定切深预测方法 |
CN112380726B (zh) * | 2020-11-26 | 2023-06-30 | 华中科技大学无锡研究院 | 基于模态耦合颤振的机器人铣削临界稳定切深预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106940746A (zh) | 基于薄壁件的铣削颤振稳定性预测的并行时域方法 | |
CN106802630A (zh) | 基于薄壁件的铣削颤振稳定性预测的并行频域方法 | |
CN102873381B (zh) | 一种基于动力学模型的高速铣削工艺参数优化方法 | |
CN106965032B (zh) | 薄壁件铣削颤振抑制方法 | |
CN104484506B (zh) | 一种基于可靠性叶瓣图的车削加工颤振预测方法 | |
TWI472402B (zh) | Tool flutter monitoring method | |
CN106141815A (zh) | 一种基于ar模型的高速铣削颤振在线辨识方法 | |
CN109093130A (zh) | 基于分区域加工的薄壁环形零件粗加工工艺优化方法 | |
CN102554326A (zh) | 一种基于叶轮叶片动态刚度的铣削精加工方法 | |
CN106217130A (zh) | 复杂曲面加工过程中铣刀状态在线监测与预警方法 | |
CN103823787B (zh) | 基于微分求积法的多车刀并行车削稳定性判定方法 | |
CN106126930A (zh) | 一种基于二分法的机床加工稳定性边界快速求解方法 | |
CN107423489A (zh) | 薄壁件铣削过程稳定性快速预测方法 | |
Mustapha et al. | A hybrid analytical model for the transverse vibration response of a micro-end mill | |
CN112016203B (zh) | 基于分段Hermite插值多项式和整体离散策略预测铣削稳定性的方法 | |
Ahmed Al-Dolaimy | Effect of cutting parameters on surface roughness in turning operations | |
DE102020130135A1 (de) | Ratterbewertungssystem | |
CN109048466B (zh) | 一种基于多频变转速的铣削颤振抑制方法 | |
CN109940462B (zh) | 铣刀切削振动变化特性的检测与高斯过程模型构建方法 | |
CN109933940A (zh) | 基于滚刀主轴振动响应模型的滚齿工艺参数优化方法 | |
CN106363463A (zh) | 基于占能比的铣削加工颤振在线监测方法 | |
CN109357826A (zh) | 一种车削刀具-工件非线性振动系统耦合特征的研究方法 | |
Yu | General influence coefficient algorithm in balancing of rotating machinery | |
Das et al. | A study on effect of process parameters on vibration of cutting tool in turning operation | |
CN109902377B (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170711 |
|
RJ01 | Rejection of invention patent application after publication |