CN112571150B - 一种用于监测薄板齿轮的薄板加工状态的非线性方法 - Google Patents

一种用于监测薄板齿轮的薄板加工状态的非线性方法 Download PDF

Info

Publication number
CN112571150B
CN112571150B CN202011431396.7A CN202011431396A CN112571150B CN 112571150 B CN112571150 B CN 112571150B CN 202011431396 A CN202011431396 A CN 202011431396A CN 112571150 B CN112571150 B CN 112571150B
Authority
CN
China
Prior art keywords
thin plate
acceleration
machining
gear
vibration
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
CN202011431396.7A
Other languages
English (en)
Other versions
CN112571150A (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.)
Central South University
Original Assignee
Central South 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 Central South University filed Critical Central South University
Priority to CN202011431396.7A priority Critical patent/CN112571150B/zh
Publication of CN112571150A publication Critical patent/CN112571150A/zh
Application granted granted Critical
Publication of CN112571150B publication Critical patent/CN112571150B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B23MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23QDETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
    • B23Q17/00Arrangements for observing, indicating or measuring on machine tools

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Electrical Machinery Utilizing Piezoelectricity, Electrostriction Or Magnetostriction (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种用于监测薄板齿轮的薄板加工状态的非线性方法,包括如下步骤:S1,获得加工中的工件振动加速度随时间函数的曲线;S2,获得薄板齿轮‑刀具系统的固有频率,阻尼比;S3,建立薄板齿轮‑刀具系统的动力学响应方程;S4,获得多组薄腹板齿轮‑刀具系统的响应特征;S5,获得多组振动能量;S6,形成单位向量内的能量分布曲线;S7,若监测到实际加工中的振动能量不在以上的振动能量分布曲线内,则为异常工况。本发明通过超声振动辅助加工中超声振动激励工件‑刀具系统后的频响特征,建立了超声加工中薄板状态的非线性监测方法,采集一组正常工况下的数据即可建立超声加工薄板齿轮的加工状态模型,可提高薄板齿轮加工稳定性及加工效率。

Description

一种用于监测薄板齿轮的薄板加工状态的非线性方法
技术领域
本发明涉及机加工领域,特别是涉及一种用于监测薄板齿轮的薄板加工状态的非线性方法。
背景技术
超声振动辅助加工能够显著降低切削力与切削温度,提高加工中的稳定性,刀具的寿命及加工效率,其广泛用于各类高强度高硬度材料的加工。但是在超声振动辅助加工中,较传统切削加工,多了超声振动源及在其激励下的非线性响应。其导致了更加复杂的力-热耦合机制,使得薄板加工状态复杂,难以预测。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。为此,本发明提出一种用于监测薄板齿轮的薄板加工状态的非线性方法,能够。
根据本发明的实施例的一种用于监测薄板齿轮的薄板加工状态的非线性方法,包括如下步骤:S1,开展超声振动辅助切削加工薄板试验,获得加工中的工件振动加速度随时间函数的曲线;S2,建立薄板加工动力学方程,获得薄板齿轮-刀具系统的固有频率ωn,阻尼比ξ;S3,依据S2获得的固有频率和阻尼比,建立薄板齿轮-刀具系统的动力学响应方程;S4,采集多个不同时刻下的位移,速度,加速度,获得多组薄板齿轮-刀具系统的响应特征;S5,对多组位移、速度和加速度,在时域内对单位切削力系数进行积分获得振动能量,将多组的三个振动能量分别进行叠加,设为ζi,i=1,…,11;S6,将S5获得的振动能量归一化处理,将振动能量的最大值归为单位向量,形成单位向量内的能量分布曲线;S7,若监测到实际加工中的振动能量不在以上的振动能量分布曲线内,则为异常工况。
根据本发明实施例的一种用于监测薄板齿轮的薄板加工状态的非线性方法,至少具有如下技术效果:通过超声振动辅助加工中超声振动激励工件-刀具系统后的频响特征,建立了超声加工中薄板状态的非线性监测方法,可建立超声加工薄板齿轮的加工状态模型,可提高薄板齿轮加工稳定性及加工效率。
在本发明一些实施例中,S2中的加工动力学方程如下:
Figure GDA0003396293120000021
在本发明一些实施例中,其中
Figure GDA0003396293120000022
Figure GDA0003396293120000023
x1(t)=xw(t)+xc(t)+av
在本发明一些实施例中,其中sv、vv、av的关系式如下:
sv=hvsin(2πfvt+θ)
vv=2πhvfvcos(2πfvt+θ)
av=-4π2hv(fv)2sin(2πfvt+θ)。
在本发明一些实施例中,S3中动力学响应方程如下:
Figure GDA0003396293120000024
在本发明一些实施例中,步骤S5中,位移、速度、加速度在时域内对单位切削力系数进行积分获得振动能量设定为ηy,i
Figure GDA0003396293120000031
Figure GDA0003396293120000032
Figure GDA0003396293120000033
Figure GDA0003396293120000034
Figure GDA0003396293120000035
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
具体实施方式
本发明实施例的一种用于监测薄板齿轮的薄板加工状态的非线性方法,包括步骤S1、S2、S3、S4、S5、S6、S7。
S1,开展超声振动辅助切削加工薄板试验,获得加工中的工件振动加速度随时间函数的曲线。
S2,建立薄板加工动力学方程,获得薄板齿轮-刀具系统的固有频率ωn,阻尼比ξ。具体过程如下:
加工动力学方程为:
Figure GDA0003396293120000036
其中ωn为薄板-刀具系统的固有频率,ap为切削深度,切削参数设定,已知,ξ为系统的阻尼比,mt为薄板-刀具系统的模态质量,通过建立有限元模型,获得模态质量,x1(t)为薄板-刀具系统的当前相对位移;τ是相邻两刀齿之间的迟滞时间,x1(t-τ)为相邻两刀之间的位移延迟。k(t)为超声振动辅助加工切削力系数,其可以通过超声振动辅助车削切削力试验,通过线性回归获得。
Figure GDA0003396293120000041
为薄板-刀具间的相对速度,
Figure GDA0003396293120000042
薄板-刀具间的相对加速度。定义切削速度方向为正向。依据采集到的薄板齿轮的加速度
Figure GDA0003396293120000043
速度
Figure GDA0003396293120000044
位移xw(t),刀具的加速度
Figure GDA0003396293120000045
速度
Figure GDA0003396293120000046
位移xc(t),通过叠加法即可获得薄板齿轮-刀具的相对加速度
Figure GDA0003396293120000047
速度
Figure GDA0003396293120000048
位移x1(t)。其中:
Figure GDA0003396293120000049
Figure GDA00033962931200000410
x1(t)=xw(t)+xc(t)+av
依据采集的多组工件与刀具间的相对位移、速度和加速度,至少2组。相对位移、速度和加速度可通过步骤S1获得的时间-加速度曲线对时间积分获得,或通过传感器检测获得。
求解方程(1),获得薄板齿轮-刀具系统的固有频率ωn,阻尼比ξ。
超声振动辅助加工中,刀具的超声振动位移、速度、加速度分别为:
sv=hvsin(2πfvt+θ);
vv=2πhvfvcos(2πfvt+θ);
av=-4π2hv(fv)2sin(2πfvt+θ);
其中hv为超声振动装置的振幅,与超声装置相关,已知的;π为圆周率;θ为初始相位角,设为0;t为时间,即设定切削开始的时刻为0,从初始时刻到关注时刻间的时间。fv为超声振动装置的频率,已知。
S3,依据S2获得的固有频率和阻尼比,建立薄板齿轮-刀具系统的动力学响应方程。薄板齿轮-刀具系统的动力学响应方程为:
Figure GDA0003396293120000051
其中初始值当t=0,
Figure GDA0003396293120000052
求解以上方程,获得薄板齿轮-刀具系统的位移y(t),速度
Figure GDA0003396293120000053
加速度
Figure GDA0003396293120000054
即获得了在当前固有频率下的薄板齿轮-刀具系统的响应特征。
S4,重复步骤S2和S3,采集多个不同时刻下的位移,速度,加速度,获得多组薄腹板齿轮-刀具系统的响应特征;位移,速度,加速度可通过位移传感器,速度传感器,加速度传感器进行采集,重复S2与S3的多次,至少11次,采集不同时刻下的位移,速度,加速度,则可获得11组固有频率下的薄腹板齿轮-刀具系统的响应特征。其中频率需要覆盖机床许用转速内的激励频率。位移,速度,加速度分别设为yi(t),
Figure GDA0003396293120000055
S5,对多组位移、速度和加速度,在时域内对单位切削力系数进行积分获得振动能量。将多组的三个振动能量分别进行叠加,设为ζi,i=1,…,11。
具体的,将S4中的11个频率下的位移,速度,加速度,分别在时域内对单位切削力系数进行积分获得对应的振动能量,设定为ηy,i
Figure GDA0003396293120000056
Figure GDA0003396293120000057
Figure GDA0003396293120000061
Figure GDA0003396293120000062
其中t1为初始积分时刻,t2为积分终止时刻。
将同个频率下的三个能量进行叠加,设为某个频率下的能量和ζi,i=1,…,11,其中
Figure GDA0003396293120000063
总能量ζ为
Figure GDA0003396293120000064
S6,将以上振动能量归一化处理,将振动能量的最大值归为单位向量。则形成单位向量内的能量分布曲线(函数);横坐标为频率,纵坐标为能量比值λ,频率振动能量/总能量,即
Figure GDA0003396293120000065
S7,若监测到实际加工中的振动能量不在以上的振动能量分布曲线内,则为异常工况。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。

Claims (2)

1.一种用于监测薄板齿轮的薄板加工状态的非线性方法,其特征在于,包括如下步骤:
S1,开展超声振动辅助切削加工薄板试验,获得加工中的工件振动加速度随时间函数的曲线;
S2,建立薄板加工动力学方程,获得薄板齿轮-刀具系统的固有频率ωn,阻尼比ξ;
S3,依据S2获得的固有频率和阻尼比,建立薄板齿轮-刀具系统的动力学响应方程;
S4,采集多个不同时刻下的位移,速度,加速度,获得多组薄板齿轮-刀具系统的响应特征;
S5,对多组位移、速度和加速度,在时域内对单位切削力系数进行积分获得振动能量,将多组的三个振动能量分别进行叠加,设为ζi,i=1,…,11;
S6,将S5获得的振动能量归一化处理,将振动能量的最大值归为单位向量,形成单位向量内的能量分布曲线;
S7,若监测到实际加工中的振动能量不在以上的振动能量分布曲线内,则为异常工况;
其中S2中的加工动力学方程如下:
Figure FDA0003396293110000011
Figure FDA0003396293110000012
Figure FDA0003396293110000013
x1(t)=xw(t)+xc(t)+av
其中sv、vv、av的关系式如下:
sv=hvsin(2πfvt+θ)
vv=2πhvfvcos(2πfvt+θ)
av=-4π2hv(fv)2sin(2πfvt+θ);
S3中动力学响应方程如下:
Figure FDA0003396293110000021
ωn为薄板-刀具系统的固有频率,ap为切削深度,ξ为系统的阻尼比,mt为薄板-刀具系统的模态质量,x1(t)为薄板-刀具系统的当前相对位移;τ是相邻两刀齿之间的迟滞时间,x1(t-τ)为相邻两刀之间的位移延迟;k(t)为超声振动辅助加工切削力系数,
Figure FDA0003396293110000022
为薄板-刀具间的相对速度,
Figure FDA0003396293110000023
薄板-刀具间的相对加速度;薄板齿轮的加速度
Figure FDA0003396293110000024
速度
Figure FDA0003396293110000025
位移xw(t);刀具的加速度
Figure FDA0003396293110000026
速度
Figure FDA0003396293110000027
位移xc(t);sv、vv和av分别为刀具的超声振动位移、速度、加速度;hv为超声振动装置的振幅;π为圆周率;θ为初始相位角;t为时间;fv为超声振动装置的频率;薄板齿轮-刀具系统的位移y(t),速度
Figure FDA0003396293110000028
加速度
Figure FDA0003396293110000029
2.根据权利要求1所述的用于监测薄板齿轮的薄板加工状态的非线性方法,其特征在于,步骤S5中,位移、速度、加速度在时域内对单位切削力系数进行积分获得振动能量设定为ηy,i
Figure FDA00033962931100000210
i=1,…,11;
Figure FDA00033962931100000211
Figure FDA0003396293110000031
Figure FDA0003396293110000032
Figure FDA0003396293110000033
i=1,…,11;
t1为初始积分时刻,t2为积分终止时刻;
ζi为某个频率下的能量和。
CN202011431396.7A 2020-12-09 2020-12-09 一种用于监测薄板齿轮的薄板加工状态的非线性方法 Active CN112571150B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011431396.7A CN112571150B (zh) 2020-12-09 2020-12-09 一种用于监测薄板齿轮的薄板加工状态的非线性方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011431396.7A CN112571150B (zh) 2020-12-09 2020-12-09 一种用于监测薄板齿轮的薄板加工状态的非线性方法

Publications (2)

Publication Number Publication Date
CN112571150A CN112571150A (zh) 2021-03-30
CN112571150B true CN112571150B (zh) 2022-02-01

Family

ID=75130715

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011431396.7A Active CN112571150B (zh) 2020-12-09 2020-12-09 一种用于监测薄板齿轮的薄板加工状态的非线性方法

Country Status (1)

Country Link
CN (1) CN112571150B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1346053A (en) * 1970-01-30 1974-02-06 Street R Indexing mechanism
CN104285139A (zh) * 2012-04-24 2015-01-14 Skf公司 轴承监测方法以及系统
CN106903547A (zh) * 2017-04-20 2017-06-30 中南大学 适用于高速铣削加工的电主轴主动振动控制装置及系统
CN107451382A (zh) * 2017-09-19 2017-12-08 湖南工学院 高速切削加工工件表面形貌的控制方法
CN108319962A (zh) * 2018-01-29 2018-07-24 安徽大学 一种基于卷积神经网络的刀具磨损监测方法
CN208117409U (zh) * 2018-04-26 2018-11-20 中南大学 用于监测切削加工及切屑形成中刀具状态的装置
CN109894923A (zh) * 2019-03-07 2019-06-18 中南大学 一种超声振动辅助加工中刀具稳定性控制方法
CN109922909A (zh) * 2016-09-09 2019-06-21 株式会社NejiLaw 切削头、切削工具及切削系统
CN109967789A (zh) * 2019-03-11 2019-07-05 杭州电子科技大学 一种超声拉刀、高效精密超声拉削装置及其拉削工艺
CN110712066A (zh) * 2019-10-22 2020-01-21 湖南工学院 适用于深孔内螺纹加工中刀具状态的监测方法
CN111843615A (zh) * 2020-06-29 2020-10-30 中南大学 一种超声振动辅助加工中材料的断裂韧性的快速识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109960865B9 (zh) * 2019-03-20 2023-05-05 湘潭大学 薄板无网格Galerkin结构动力响应分析的GPU加速方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1346053A (en) * 1970-01-30 1974-02-06 Street R Indexing mechanism
CN104285139A (zh) * 2012-04-24 2015-01-14 Skf公司 轴承监测方法以及系统
CN109922909A (zh) * 2016-09-09 2019-06-21 株式会社NejiLaw 切削头、切削工具及切削系统
CN106903547A (zh) * 2017-04-20 2017-06-30 中南大学 适用于高速铣削加工的电主轴主动振动控制装置及系统
CN107451382A (zh) * 2017-09-19 2017-12-08 湖南工学院 高速切削加工工件表面形貌的控制方法
CN108319962A (zh) * 2018-01-29 2018-07-24 安徽大学 一种基于卷积神经网络的刀具磨损监测方法
CN208117409U (zh) * 2018-04-26 2018-11-20 中南大学 用于监测切削加工及切屑形成中刀具状态的装置
CN109894923A (zh) * 2019-03-07 2019-06-18 中南大学 一种超声振动辅助加工中刀具稳定性控制方法
CN109967789A (zh) * 2019-03-11 2019-07-05 杭州电子科技大学 一种超声拉刀、高效精密超声拉削装置及其拉削工艺
CN110712066A (zh) * 2019-10-22 2020-01-21 湖南工学院 适用于深孔内螺纹加工中刀具状态的监测方法
CN111843615A (zh) * 2020-06-29 2020-10-30 中南大学 一种超声振动辅助加工中材料的断裂韧性的快速识别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
切削加工状态监测技术的现状与发展;薛鸿健等;《机械与电子》;19900930;第36-39页 *
高精密球加工过程中的状态监测;安平;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20140630;全文 *

Also Published As

Publication number Publication date
CN112571150A (zh) 2021-03-30

Similar Documents

Publication Publication Date Title
Xiang et al. Ultrasonic longitudinal-torsional vibration-assisted cutting of Nomex® honeycomb-core composites
CN111644900B (zh) 一种基于主轴振动特征融合的刀具破损实时监测方法
Wan et al. Study of static and dynamic ploughing mechanisms by establishing generalized model with static milling forces
Haiyan et al. A mechanistic model for cutting force in helical milling of carbon fiber-reinforced polymers
Wan et al. Investigation on milling chatter identification at early stage with variance ratio and Hilbert–Huang transform
Choi et al. On-line chatter detection using wavelet-based parameter estimation
Pratap et al. Modeling cutting force in micro-milling of ti-6al-4 v titanium alloy
Zhang et al. Modeling of the effect of tool edge radius on surface generation in elliptical vibration cutting
Attanasio et al. Finite element simulation of high speed micro milling in the presence of tool run-out with experimental validations
CN104778333A (zh) 一种螺旋立铣刀正交车铣加工三维稳定性建模方法
CN105414616B (zh) 螺旋铣孔过程中切削力预报和稳定性判别方法
Ding et al. Dynamic surface generation modeling of two-dimensional vibration-assisted micro-end-milling
WO2020143203A1 (zh) 一种深孔镗削加工颤振的在线监测与抑制方法
Liu et al. Hybrid model-and signal-based chatter detection in the milling process
CN104794337A (zh) 一种基于边界条件判断的正交车铣加工端面刃切削力建模方法
CN112571150B (zh) 一种用于监测薄板齿轮的薄板加工状态的非线性方法
CN105160128A (zh) 一种曲线端铣加工过程切削力的预测方法
Maroju et al. Investigations on feasibility of low-frequency vibration-assisted turning
CN109376440A (zh) 一种高效铣削弱刚度零件的加工控制方法
Li et al. Chatter detection and identification in high-efficient robotic milling CFRP composites using acoustic emission technique
Zhai et al. Influence of cutting parameters on force coefficients and stability in plunge milling
Shi et al. A chatter mitigation technique in milling based on H∞-ADDPMS and piezoelectric stack actuators
CN107679335A (zh) 考虑刀具振动下动态切屑厚度的实时切削力系数计算方法
Susanto et al. Vibration analysis in milling of thin-walled workpieces using Hilbert-Huang transform
Chen et al. Quantitative impacts of regenerative vibration and abrasive wheel eccentricity on surface grinding dynamic performance

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