CN109746762A - 一种深孔镗削加工颤振的在线监测与抑制方法 - Google Patents
一种深孔镗削加工颤振的在线监测与抑制方法 Download PDFInfo
- Publication number
- CN109746762A CN109746762A CN201910013374.XA CN201910013374A CN109746762A CN 109746762 A CN109746762 A CN 109746762A CN 201910013374 A CN201910013374 A CN 201910013374A CN 109746762 A CN109746762 A CN 109746762A
- Authority
- CN
- China
- Prior art keywords
- boring
- flutter
- formula
- processing
- deep hole
- 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
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23B—TURNING; BORING
- B23B47/00—Constructional features of components specially designed for boring or drilling machines; Accessories therefor
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, 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
- B23Q11/00—Accessories fitted to machine tools for keeping tools or parts of the machine in good working condition or for cooling work; Safety devices specially combined with or arranged in, or specially adapted for use in connection with, machine tools
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, 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/00—Arrangements for observing, indicating or measuring on machine tools
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Machine Tool Sensing Apparatuses (AREA)
Abstract
本发明公开了一种深孔镗削加工颤振的在线监测方法,颤振的监测和抑制方法:先构建深孔镗削加工的在线监测***;构建浮动镗刀镗削加工动力学模型;建立驱动电机电流信号与镗削加工动态特性的关系;再针对已知***外部输入的情况,来完全确定***在未来各个时刻的状态,通过对状态变量的描述和求解建立外部输入输出变量和内部状态变量之间的关系;然后将颤振信号与正常加工信号进行对比,实现对早期颤振的在线监测。本发明中的在线监测和抑制方法对深孔镗削过程中的颤振进行实时监测和抑制,可以提高工件的加工精度,降低工件的不合格率。
Description
技术领域
本发明涉及深孔镗削加工技术领域,更具体的说是涉及一种深孔镗削加工颤振的在线监测与抑制方法。
背景技术
随着科技的发展,深孔的运用越来越广泛,例如:工业中齿轮轴、曲轴和喷油器提供润滑油的孔;武器工业的枪管炮管;航空工业中发动机的冷却孔;医疗工业中的空心植入物或者外科医疗器械等等。深孔零件在机械加工中属于较高精度要求的零件,精镗常作为深孔零件最后一道工序用来保证孔的精度,但由于镗杆的悬伸量较大,往往会导致镗削***的结构刚度很低,并且在加工的过程中很容易发生颤振。然而颤振会极大地影响精镗的精度,严重时甚至会使零件报废。因此,在线监测颤振成为了提高深孔加工效率和加工精度的关键因素。
目前,深孔镗削加工中主要应用的监测方法分为直接法和间接法。直接法使用外部传感器如麦克风或声发射传感器来实现对颤振的在线监测,直接法装置简单,但是直接法中麦克风和声发生传感器往往会受到附近噪声的影响,当刀具进入零件加工时,在一定频率下的声音信号会被放大,从而导致虚假的颤振报警,降低了监测方法的准确性。间接法通过对安装在主轴、支架的传感器信号(即力、扭矩、振动)进行评估来检测颤振,大多数安装在主轴和支架上的传感器虽然能够保证良好的精度,但是,在安装后会降低主轴的动态刚度。
因此,研究出一种既可以避免降低主轴的动态刚度又能准确的对颤振现象进行在线监测和抑制的方法是本领域技术人员亟需解决的问题。
发明内容
有鉴于此,本发明提供了一种可以避免降低主轴的动态刚度,且监测准确性高的深孔镗削加工颤振的在线监测与抑制方法。
为了实现上述目的,本发明采用如下技术方案:一种深孔镗削加工颤振的在线监测方法,颤振的在线监测方法如下:
步骤(1),构建深孔镗削加工颤振的在线监测***,其包括:深孔镗床、电流放大器、数据采集卡、工控机和控制器;所述电流放大器、数据采集卡、工控机和控制器依次排布且电性连接;
所述深孔镗床上固定连接夹具部、电机及镗杆,所述电机与所述镗杆固定连接,工件置于所述夹具部内,所述镗杆的端部设置浮动镗刀,所述浮动镗刀对所述工件进行深孔镗削加工;所述夹具部及电机分别与所述电流放大器电性连接;所述夹具部及电机分别与所述控制器电性连接;
步骤(2),构建浮动镗刀镗削加工动力学模型:包括以下步骤:
(2.1),所述浮动镗刀为对称结构,选取其中一半刀具作为研究对象;
(2.2),依据动力学公式:对刀具进行动力学分析;其中,M为***质量、C为阻尼系数、k为刚度矩阵,β为刀具切削方向与竖直平面的夹角,F(t)为颤振状态下切削力的大小,x(t)为浮动镗刀位移量,为浮动镗刀镗削速度,为浮动镗刀镗削加速度;
(2.3),在再生型颤振的影响下计算得到浮动镗刀实际切削量为:y(t)=y0-[x(t)-x(t-T)];其中,y0为理论切削量,单位mm,T为镗刀旋转的周期,单位s,x(t)为浮动镗刀位移量,x(t-T)为镗削一个周期后浮动镗刀位移量;
(2.4),假设所述颤振状态下切削力的大小为:F(t)=ksby(t);其中,b为镗削宽度系数,ks为镗削刚度系数;
(2.5),将所述步骤(2.3)和(2.4)中的公式代入步骤(2.2)中的动力学公式中,对y(t)进行拉普拉斯变换得到传递函数1+(1-e-sT)ksbΦ(s)=0;其中,Φ(s)是对x(t)进行拉普拉斯变换得到传递函数,其中,ζ为镗削振动***的阻尼比,p为镗削***的固有频率,拉普拉斯变换是将一个有参数实数t(t≥0)的函数转换为一个参数为复数s的函数;
(2.6),依据所述步骤(2.5),当s=ji时,再生型颤振处于临界状态,利用Nyquist稳定性判据,将s=ji带入到公式1+(1-e-sT)ksbΦ(s)=0中计算得到再生型颤振临界状态下的公式:其中,j为虚部常数,i为虚数单位;
(2.7),对所述步骤(2.6)中的临界状态公式进行求解得到:其中
(2.8),根据所述步骤(2.6)中的临界状态公式以及镗削加工周期T的计算公式,计算得到镗削加工周期T:
(2.9),根据所述步骤(2.6)中的临界状态公式以及步骤(2.8)中的加工周期T,计算得到临界状态下的转速和极限切削厚度:
步骤(3),建立驱动电机电流信号与镗削加工动态特性的关系,包括如下步骤:
(3.1),计算主轴电机转矩和有效电流的比例关系:Mm=Ki=Mf+Mc+Cω+Jε;其中,Mm为主轴电机转矩,k为转矩常数,i为主轴电机输出有效电流,Mf为摩擦切割转矩,Mc为机床传动***等效到电机上的转矩,ω为角速度,J为转动惯量,ε为角加速度,C为***等效阻尼;
(3.2),假设镗削***净切削功率为p,***载荷磨损系数为k,Mfω=kp;
(3.3),镗削***随着时间的变化,机床动态特性也会发生改变,假设随时间变化的机床动态特性为:[x1(t),x2(t),...,xn(t)];
(3.4),根据所述步骤(2.4)和步骤(3.3)计算得到切削力与机床动态特性的关系式:F(t)=ksby(t)=f[x1(t),x2(t),...,xn(t)];f表示F(t)是关于x1(t),x2(t),...,xn(t)的一个函数;
(3.5),根据所述步骤(3.4)中的公式计算得到镗削***的净切削功率:p=F(t)v(t)=f[x1(t),x2(t),...,xn(t)]v(t);其中,v(t)为加工过程中的瞬时速度;
(3.6),将所述步骤(3.2)-(3.5)的公式带入到步骤(3.1)的公式中,并进行整理得到:Kiω=kf[x1(t),x2(t),...,xn(t)]v(t)+Mcω+Cω2+Jεω;由上式可知,***的动态特性发生改变会直接影响驱动电机电流信号的改变;
步骤(4),提供一种基于状态空间法估计电流放大器和***结构模态的干扰,其中,状态空间建立具体包括如下步骤:
(4.1),所述步骤(3.6)中整理得到的公式:Kiω=kf[x1(t),x2(t),...,xn(t)]v(t)+Mcω+cω2+Jεω;为驱动电机电流信号与***动态特性相关的关系式;
(4.2),***加工连续时域中的空间模型由状态矩阵A、输入矩阵B、输出矩阵C、直接传递矩阵D组成;
(4.3),假设镗削加工***是一个线性定常***,该***的状态方程和输出方程的经验公式为:x(t)=Ax(t)+Bu(t)+w(t),y(t)=Cx(t)+Du(t)+v(t),其中,w(t)是***干扰噪声,v(t)是测量噪声,u(t)是***中的输入向量;电流放大器和***结构模态的干扰导致的扰动转矩为τ(t),在实际计算中为简化计算不考虑噪声的干扰,且***中直接传递矩阵D也忽略不计;
(4.4),在不考虑噪声的情况下镗削加工***的传递函数为:其中I是单位矩阵;
(4.5),设在线监测***中状态矩阵为A1、输入矩阵为B1、输出矩阵为C1,***中反馈增益系数为U,***检测的周期为T,根据所述步骤(4.3)中的经验公式,计算得到监测***的状态函数为x1(t)=A1x1(t)+B1u1(t)=A1x1(t-T)+U[y(t)-y1(t)];输出函数为:y1(t)=C1x1(t);
(4.6),实际监测***与理论检测***的误差为:e(t)=C1x1(t)-Cx(t),经过一个周期后的监测***误差为:e(t+T)=(A1-UC1)e(t),由误差公式表明滤波器误差传递与输入矩阵B1无关;
(4.7),根据所述步骤(4.5)中的公式,得到监测***的传递函数为:
(4.8),根据所述步骤(4.4)、(4.7)中的公式得到测量转矩与扰动转矩之间的传递函数为:
步骤(5),通过步骤(4)去除电流放大器与结构模态的干扰后,将电流放大器采集到的实际加工电流信号与理论加工电流信号进行对比,实现对早期颤振的在线监测。
采用上述技术方案的有益效果是,本发明中通过使用电流传感器采集驱动电机电流信号,采用流形学习算法提取颤振特征向量实现深孔镗削的在线监测,及时发现颤振现象并采取相应的措施进行抑制,提高镗削加工的精度。
优选的,所述步骤5中首先通过采用流形学习算法,对所述步骤(4)中采集到的驱动电机电流信号进行降维处理,再将采集到的颤振信号与正常加工信号进行对比,然后提取颤振信号的特征向量,观察加工中电流信号的变化,进而实现对早期颤振的在线监测。
采用上述技术方案的有益效果是,选用上述方法对颤振现象进行在线监测,可以实时监测颤振的现象,并及时的对颤振现象进行抑制,提高镗削加工的精度,更好的满足深孔镗削加工精度的要求。
再生颤振是指在振动很大的场合,多数是由于再生效应的反馈机制所引起的自激振动。
一种深孔镗削加工颤振的抑制方法,包括以下步骤:
步骤(1),采用所述深孔镗削加工颤振的在线监测方法获得主轴转速、切削厚度的关系式,其中,影响主轴转速、切削厚度的因素有:***阻尼C、***刚度k、周期T;
步骤(2),通过改变所述步骤(1)中的参数值,进而实现对颤振的抑制。
采用上述技术方案的有益效果是,通过对影响颤振因素中的参数进行相应的调整,实现对颤振现象的抑制;通过对颤振现象的监测和抑制,使镗削加工的精度得到提升,降低工件的不合格率。
本发明的有益效果:
(1)本发明基于状态空间法,建立测量转矩与扰动转矩之间的传递函数,可以减少电流放大器和***结构模态干扰的影响,提高驱动电机电流信号的观测准确性,进而通过电流传感器采集驱动电机电流信号,并采用流形学习算法提取颤振特征向量实现深孔镗削的在线监测;
(2)本发明通过综合考虑再生型颤振机理、浮动镗刀结构,构建深孔镗削动力学模型,建立主轴转速与极限切深之间的关系,并适当调整相应的镗削参数,进而实现对深孔镗削颤振的抑制;
(3)本发明中通过对镗削加工过程中颤振现象的在线监测和抑制,进而使镗削加工的精度更加准确,更好的满足工件精度的要求,降低工件的不合格率,提高工作效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的镗削颤振在线监测示意图;
图2附图为本发明提供的镗削在线检测结构图;
图3附图为本发明提供的浮动镗刀动力学模型;
图4附图为本发明提供的镗削加工***反馈示意图;
图5附图为本发明提供的镗削加工***传递框图;
图6附图为本发明提供的监测***传递***框图。
其中,图中,
1-夹具;2-工件;3-镗杆;4-电机;5-底座;6-电流放大器;7-数据采集卡;8-工控机;9-控制器。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种深孔镗削加工颤振的在线监测方法,颤振的在线监测方法如下:
步骤(1),构建深孔镗削加工颤振的在线监测***,其包括:深孔镗床、电流放大器6、数据采集卡7、工控机8和控制器9;电流放大器6、数据采集卡7、工控机8和控制器9依次排布且电性连接;
深孔镗床上固定连接夹具部1、电机4及镗杆3,电机4与镗杆3固定连接,工件2置于夹具部1内,镗杆3的端部设置浮动镗刀,浮动镗刀对工件2进行深孔镗削加工;夹具部1及电机4分别与电流放大器6电性连接;夹具部1及电机4分别与控制器9电性连接;电机上还连接电流传感器,电流传感器与电流放大器电性连接;
在电机4处安装电流传感器,采集电机4的输出电流信号,电流传感器与电流放大器6电性连接,将电流信号放大至可以读取的范围;电流信号放大器6与数据采集卡7连接,将放大后的数据进行保存和初步处理;然后在工控机上对电流信号进行判断,之后再通过控制器对镗削加工进行控制。
步骤(2),构建浮动镗刀镗削加工动力学模型:包括以下步骤:
(2.1),浮动镗刀为对称结构,选取其中一半刀具作为研究对象;
(2.2),依据动力学公式:对刀具进行动力学分析;其中,M为***质量、C为阻尼系数、k为刚度矩阵,β为刀具切削方向与竖直平面的夹角,F(t)为颤振状态下切削力的大小,x(t)为浮动镗刀位移量,为浮动镗刀镗削速度,为浮动镗刀镗削加速度;
(2.3),在再生型颤振的影响下计算得到浮动镗刀实际切削量为:y(t)=y0-[x(t)-x(t-T)];其中,y0为理论切削量,单位mm,理论切削量是假设得到的,T为镗刀旋转的周期,单位s,x(t)为浮动镗刀位移量,x(t-T)为镗削一个周期后浮动镗刀位移量;
(2.4),假设颤振状态下切削力的大小为:F(t)=ksby(t);其中,b为镗削宽度系数;ks为镗削刚度系数;
(2.5),对步骤(2.3)和(2.4)中的公式代入步骤(2.2)中的动力学公式中,对y(t)进行拉普拉斯变换得到传递函数1+(1-e-sT)ksbΦ(s)=0;其中,Φ(s)是对x(t)进行拉普拉斯变换得到传递函数,其中,ζ为镗削振动***的阻尼比,p为镗削***的固有频率,
(2.6),依据所述步骤(2.5),当s=ji时,再生型颤振处于临界状态,利用Nyquist稳定性判据,将s=ji带入到公式1+(1-e-sT)ksbΦ(s)=0中计算得到再生型颤振临界状态下的公式:其中,j为虚部常数,i为虚数单位;
(2.7),对步骤(2.6)中的临界状态公式进行求解得到:其中
(2.8),根据步骤(2.6)中的临界状态公式以及镗削加工周期T的计算公式,计算得到镗削加工周期T:
(2.9),根据步骤(2.6)中的临界状态公式以及步骤(2.8)中的加工周期T,计算得到临界状态下的转速和极限切削厚度:
步骤(3),建立驱动电机电流信号与镗削加工动态特性的关系,如图3所示,其中,D表示实际镗削表面;E表示理想镗削表面;F表示上一转镗削表面;G表示上一转理想镗削表面,包括如下步骤:
(3.1),计算主轴电机转矩和有效电流的比例关系:Mm=Ki=Mf+Mc+Cω+Jε;其中,Mm为主轴电机转矩,k为转矩常数,i为主轴电机输出有效电流,Mf为摩擦切割转矩,Mc为机床传动***等效到电机上的转矩,ω为角速度,J为转动惯量,ε为角加速度,C为***等效阻尼;
(3.2),假设镗削***净切削功率为p,***载荷磨损系数为k,Mfω=kp;
(3.3),镗削***随着时间的变化,机床动态特性也会发生改变,假设随时间变化的机床动态特性为:[x1(t),x2(t),...,xn(t)];
(3.4),根据步骤(2.4)和步骤(3.3)计算得到切削力与机床动态特性的关系式:F(t)=ksby(t)=f[x1(t),x2(t),...,xn(t)];
(3.5),根据步骤(3.4)中的公式计算得到镗削***的净切削功率:p=F(t)v(t)=f[x1(t),x2(t),...,xn(t)]v(t);
(3.6),将步骤(3.2)-(3.5)的公式带入到步骤(3.1)的公式中,并进行整理得到:Kiω=kf[x1(t),x2(t),...,xn(t)]v(t)+Mcω+Cω2+Jεω;由上式可知,***的动态特性发生改变会直接影响驱动电机电流信号的改变;
步骤(4),提供一种基于状态空间法估计电流放大器和***结构模态的干扰,进而提高驱动电机电流信号的观测的准确性,具体针对已知***外部输入的情况,来完全确定***在未来各个时刻的状态,通过对状态变量的描述和求解建立外部输入输出变量和内部状态变量之间的关系;镗削加工***反馈如图4所示,其中,状态空间建立具体包括如下步骤:
(4.1),步骤(3.6)中整理得到的公式:Kiω=kf[x1(t),x2(t),...,xn(t)]v(t)+Mcω+Cω2+Jεω;为驱动电机电流信号与***动态特性相关的关系式;
(4.2),***加工连续时域中的空间模型由状态矩阵A、输入矩阵B、输出矩阵C、直接传递矩阵D组成;
(4.3),假设镗削加工***是一个线性定常***,该***的状态方程和输出方程的经验公式为:x(t)=Ax(t)+Bu(t)+w(t),y(t)=Cx(t)+Du(t)+v(t),其中,w(t)是***干扰噪声,v(t)是测量噪声,u(t)是***中的输入向量;电流放大器和***结构模态的干扰导致的扰动转矩为τ(t),在实际计算中为简化计算不考虑噪声的干扰,且***中直接传递矩阵D也忽略不计,镗削加工***传递框图如图5所示,其中,H表示扰动转矩;I表示额定电流;
(4.4),在不考虑噪声的情况下镗削加工***的传递函数为:其中I是单位矩阵;
(4.5),设在线监测***中状态矩阵为A1、输入矩阵为B1、输出矩阵为C1,***中反馈增益系数为U,***检测的周期为T,根据步骤(4.3)中的经验公式,计算得到监测***的状态函数为x1(t)=A1x1(t)+B1u1(t)=A1x1(t-T)+U[y(t)-y1(t)];输出函数为:y1(t)=C1x1(t);
(4.6),实际监测***与理论检测***的误差为:e(t)=C1x1(t)-Cx(t),经过一个周期后的监测***误差为:e(t+T)=(A1-UC1)e(t),由误差公式表明滤波器误差传递与输入矩阵B1无关;
(4.7)建立监测***信号传递***框图,如图6所示,其中,I表示额定电流;J表示测量转矩;
(4.8),根据步骤(4.5)中的公式,得到监测***的传递函数为:
(4.9),根据步骤(4.4)、(4.7)中的公式得到测量转矩与扰动转矩之间的传递函数为:
步骤(5),首先通过采用流形学习算法,通过步骤(4)去除电流放大器与结构模态的干扰后,再对步骤(4)中采集到的驱动电机电流信号进行降维处理,再将电流放大器采集到的实际加工电流与理论加工电流进行对比,然后提取颤振信号的特征向量,观察加工中电流信号的变化,进而实现对早期颤振的在线监测。其中,理论电流信号在加工过程中是不变的定值。
一种深孔镗削加工颤振的抑制方法,包括以下步骤:
步骤(1),采用深孔镗削加工颤振的在线监测方法获得主轴转速、切削厚度的关系式,其中,影响主轴转速、切削厚度的因素有:***阻尼C、***刚度k、周期T;
步骤(2),通过改变步骤(1)中的参数值,进而实现对颤振的抑制。
进一步地,通过改变稳定性叶瓣图(SLD)选择合适的工艺参数a、τ,可以避免颤振问题;利用特殊刀具几何形状或主轴变速技术,通过改变加工周期可以减少再生效应;通过设计特殊的切削刃几何形状可以增加过程阻尼;***的刚度可以通过重新设计***、采用特殊夹具、使用加强筋、采用高性能材料等不同的方法提高;通过被动(使用阻尼器或高内阻尼材料)或主动技术(主动结构颤振抑制,主动刀具,主动主轴***和主动夹具)来提高***的阻尼,例如在镗杆上安装压电驱动器或采用磁/电流变液提高***阻尼;周期需要在加工过程中改变主轴转速,进而得到改变。通过改变上述参数可以实现对颤振的抑制。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (3)
1.一种深孔镗削加工颤振的在线监测方法,其特征在于,颤振的在线监测方法如下:
步骤(1),构建深孔镗削加工颤振的在线监测***,其包括:深孔镗床、电流放大器、数据采集卡、工控机和控制器;所述电流放大器、数据采集卡、工控机和控制器依次电性连接;
所述深孔镗床上固定连接夹具部、电机及镗杆,所述电机与所述镗杆固定连接,工件置于所述夹具部内,所述镗杆的端部设置浮动镗刀,所述浮动镗刀对所述工件进行深孔镗削加工;所述夹具部及电机分别与所述电流放大器电性连接;所述夹具部及电机分别与所述控制器电性连接;
步骤(2),构建浮动镗刀镗削加工动力学模型:包括以下步骤:
(2.1),所述浮动镗刀为对称结构,选取其中一半刀具作为研究对象;
(2.2),依据动力学公式:对刀具进行动力学分析;其中,M为***质量、C为阻尼系数、k为刚度矩阵,β为刀具切削方向与竖直平面的夹角,F(t)为颤振状态下切削力的大小,x(t)为浮动镗刀位移量,为浮动镗刀镗削速度,为浮动镗刀镗削加速度;
(2.3),在再生型颤振的影响下计算得到浮动镗刀实际切削量为:y(t)=y0-[x(t)-x(t-T)];其中,y0为理论切削量,单位mm,T为镗刀旋转的周期,单位s,x(t)为浮动镗刀位移量,x(t-T)为镗削一个周期后浮动镗刀位移量;
(2.4),假设所述颤振状态下切削力的大小为:F(t)=ksby(t);其中,b为镗削宽度系数,ks为镗削刚度系数;
(2.5),将所述步骤(2.3)和(2.4)中的公式代入步骤(2.2)中的动力学公式中,对y(t)进行拉普拉斯变换得到传递函数1+(1-e-sT)ksbΦ(s)=0;其中,Φ(s)是对x(t)进行拉普拉斯变换得到传递函数,其中,ζ为镗削振动***的阻尼比,p为镗削***的固有频率,
(2.6),依据所述步骤(2.5),当s=ji时,再生型颤振处于临界状态,利用Nyquist稳定性判据,将s=ji带入到公式1+(1-e-sT)ksbΦ(s)=0中计算得到再生型颤振临界状态下的公式:(-j2+2ζpji+p2)M+[1-cos(jT)+i sin(jT)]ksb cosβ=0;其中,j为虚部常数,i为虚数单位;
(2.7),对所述步骤(2.6)中的临界状态公式进行求解得到:其中
(2.8),根据所述步骤(2.7)中的临界状态公式推导得到镗削加工周期T:
a为任意常数;
(2.9),根据所述步骤(2.6)中的临界状态公式以及步骤(2.8)中的加工周期T,计算得到临界状态下的转速和极限切削厚度:
步骤(3),建立驱动电机电流信号与镗削加工动态特性的关系,包括如下步骤:
(3.1),计算主轴电机转矩和有效电流的比例关系:Mm=Ki=Mf+Mc+Cω+Jε;其中,Mm为主轴电机转矩,k为转矩常数,i为主轴电机输出有效电流,Mf为摩擦切割转矩,Mc为机床传动***等效到电机上的转矩,ω为角速度,J为转动惯量,ε为角加速度,C为***等效阻尼;
(3.2),假设镗削***净切削功率为p,***载荷磨损系数为k,Mfω=kp;
(3.3),镗削***随着时间的变化,机床动态特性也会发生改变,假设随时间变化的机床动态特性为:[x1(t),x2(t),...,xn(t)];
(3.4),根据所述步骤(2.4)和步骤(3.3)计算得到切削力与机床动态特性的关系式:F(t)=ksby(t)=f[x1(t),x2(t),...,xn(t)];
(3.5),根据所述步骤(3.4)中的公式计算得到镗削***的净切削功率:p=F(t)v(t)=f[x1(t),x2(t),...,xn(t)]v(t);其中,v(t)为加工过程中的瞬时速度;
(3.6),将所述步骤(3.2)-(3.5)的公式带入到步骤(3.1)的公式中,并进行整理得到:Kiω=kf[x1(t),x2(t),...,xn(t)]v(t)+Mcω+Cω2+Jεω;由上式可知,***的动态特性发生改变会直接影响驱动电机电流信号的改变;
步骤(4),提供一种基于状态空间法估计电流放大器和***结构模态的干扰,其中,状态空间建立具体包括如下步骤:
(4.1),所述步骤(3.6)中整理得到的公式:Kiω=kf[x1(t),x2(t),...,xn(t)]v(t)+Mcω+Cω2+Jεω;为驱动电机电流信号与***动态特性相关的关系式;
(4.2),***加工连续时域中的空间模型由状态矩阵A、输入矩阵B、输出矩阵C、直接传递矩阵D组成;
(4.3),假设镗削加工***是一个线性定常***,该***的状态方程和输出方程的经验公式为:x(t)=Ax(t)+Bu(t)+w(t),y(t)=Cx(t)+Du(t)+v(t),其中,w(t)是***干扰噪声,v(t)是测量噪声,u(t)是***中的输入向量;电流放大器和***结构模态的干扰导致的扰动转矩为τ(t),在实际计算中为简化计算不考虑噪声的干扰,且***中直接传递矩阵D也忽略不计;
(4.4),在不考虑噪声的情况下镗削加工***的传递函数为:其中I是单位矩阵;
(4.5),设在线监测***中状态矩阵为A1、输入矩阵为B1、输出矩阵为C1,***中反馈增益系数为U,***检测的周期为T,根据所述步骤(4.3)中的经验公式,计算得到监测***的状态函数为x1(t)=A1x1(t)+B1u1(t)=A1x1(t-T)+U[y(t)-y1(t)];输出函数为:y1(t)=C1x1(t);
(4.6),实际监测***与理论检测***的误差为:e(t)=C1x1(t)-Cx(t),经过一个周期后的监测***误差为:e(t+T)=(A1-UC1)e(t),由误差公式表明滤波器误差传递与输入矩阵B1无关;
(4.7),根据所述步骤(4.5)中的公式,得到监测***的传递函数为:
(4.8),根据所述步骤(4.4)、(4.7)中的公式得到测量转矩与扰动转矩之间的传递函数为:
步骤(5),通过步骤(4)去除电流放大器与结构模态的干扰后,将电流放大器采集到的实际加工电流与理论加工电流进行对比,实现对早期颤振的在线监测。
2.根据权利要求1所述的一种深孔镗削加工颤振的在线监测方法,其特征在于,所述步骤5中首先通过采用流形学习算法,对所述步骤(4)中采集到的驱动电机电流信号进行降维处理,再将电流放大器采集到的实际加工电流与理论加工电流进行对比,然后提取颤振信号的特征向量,观察加工中电流信号的变化,进而实现对早期颤振的在线监测。
3.一种深孔镗削加工颤振的抑制方法,其特征在于,包括以下步骤:
步骤(1),采用权利要求1至2任一项所述深孔镗削加工颤振的在线监测方法获得主轴转速、切削厚度的关系式,其中,影响主轴转速、切削厚度的因素有:***阻尼C、***刚度k、周期T;
步骤(2),通过改变所述步骤(1)中的参数值,进而实现对颤振的抑制。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910013374.XA CN109746762B (zh) | 2019-01-07 | 2019-01-07 | 一种深孔镗削加工颤振的在线监测与抑制方法 |
PCT/CN2019/095846 WO2020143203A1 (zh) | 2019-01-07 | 2019-07-12 | 一种深孔镗削加工颤振的在线监测与抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910013374.XA CN109746762B (zh) | 2019-01-07 | 2019-01-07 | 一种深孔镗削加工颤振的在线监测与抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109746762A true CN109746762A (zh) | 2019-05-14 |
CN109746762B CN109746762B (zh) | 2021-02-02 |
Family
ID=66404569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910013374.XA Expired - Fee Related CN109746762B (zh) | 2019-01-07 | 2019-01-07 | 一种深孔镗削加工颤振的在线监测与抑制方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109746762B (zh) |
WO (1) | WO2020143203A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110346130A (zh) * | 2019-07-19 | 2019-10-18 | 北京理工大学 | 一种基于经验模态分解和时频多特征的镗削颤振检测方法 |
CN110434676A (zh) * | 2019-07-29 | 2019-11-12 | 北京理工大学 | 一种多传感器时频特征融合的镗削颤振监测方法 |
CN110712066A (zh) * | 2019-10-22 | 2020-01-21 | 湖南工学院 | 适用于深孔内螺纹加工中刀具状态的监测方法 |
WO2020143203A1 (zh) * | 2019-01-07 | 2020-07-16 | 北京理工大学 | 一种深孔镗削加工颤振的在线监测与抑制方法 |
CN112180833A (zh) * | 2020-09-07 | 2021-01-05 | 清华大学 | 数控机床及其自适应加工的控制方法和装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5855824A (ja) * | 1981-09-30 | 1983-04-02 | Toshiba Corp | 振動監視装置 |
US7341410B2 (en) * | 2003-03-10 | 2008-03-11 | Foster-Miller, Inc. | Dynamical instrument for machining |
CN101722438A (zh) * | 2008-10-28 | 2010-06-09 | 大隈株式会社 | 机床的振动抑制方法和装置 |
CN102490086A (zh) * | 2011-10-28 | 2012-06-13 | 浙江大学 | 镗杆工作状态实时监测*** |
CN104076734A (zh) * | 2014-06-26 | 2014-10-01 | 天津职业技术师范大学 | 一种铣削加工颤振在线寻优方法 |
CN105678043A (zh) * | 2014-11-17 | 2016-06-15 | 天津职业技术师范大学 | 一种考虑刚度时变的大切除率铣削颤振监测方法 |
CN107330375A (zh) * | 2017-06-05 | 2017-11-07 | 天津大学 | 基于经验模态分解和样本熵特征的机床颤振在线监测方法 |
CN108490789A (zh) * | 2018-05-08 | 2018-09-04 | 西安交通大学 | 一种电主轴铣削颤振鲁棒主动控制方法及其主动控制*** |
CN108846242A (zh) * | 2018-07-09 | 2018-11-20 | 西北工业大学 | 基于预拉应力施加的薄壁件铣削颤振抑制方法 |
CN109048466A (zh) * | 2018-09-03 | 2018-12-21 | 西安交通大学 | 一种基于多频变转速的铣削颤振抑制方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3421562B2 (ja) * | 1997-12-26 | 2003-06-30 | オークマ株式会社 | 刃具の振れ検出方法 |
JP2011154588A (ja) * | 2010-01-28 | 2011-08-11 | Brother Industries Ltd | 工作機械 |
CN106802630B (zh) * | 2017-03-14 | 2019-01-04 | 东北大学 | 基于薄壁件的铣削颤振稳定性预测的并行频域方法 |
CN107457609B (zh) * | 2017-08-25 | 2019-12-03 | 西安交通大学 | 基于刚度变化的铣削颤振抑制方法及铣削颤振优化*** |
CN109746762B (zh) * | 2019-01-07 | 2021-02-02 | 北京理工大学 | 一种深孔镗削加工颤振的在线监测与抑制方法 |
-
2019
- 2019-01-07 CN CN201910013374.XA patent/CN109746762B/zh not_active Expired - Fee Related
- 2019-07-12 WO PCT/CN2019/095846 patent/WO2020143203A1/zh active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5855824A (ja) * | 1981-09-30 | 1983-04-02 | Toshiba Corp | 振動監視装置 |
US7341410B2 (en) * | 2003-03-10 | 2008-03-11 | Foster-Miller, Inc. | Dynamical instrument for machining |
CN101722438A (zh) * | 2008-10-28 | 2010-06-09 | 大隈株式会社 | 机床的振动抑制方法和装置 |
CN102490086A (zh) * | 2011-10-28 | 2012-06-13 | 浙江大学 | 镗杆工作状态实时监测*** |
CN104076734A (zh) * | 2014-06-26 | 2014-10-01 | 天津职业技术师范大学 | 一种铣削加工颤振在线寻优方法 |
CN105678043A (zh) * | 2014-11-17 | 2016-06-15 | 天津职业技术师范大学 | 一种考虑刚度时变的大切除率铣削颤振监测方法 |
CN107330375A (zh) * | 2017-06-05 | 2017-11-07 | 天津大学 | 基于经验模态分解和样本熵特征的机床颤振在线监测方法 |
CN108490789A (zh) * | 2018-05-08 | 2018-09-04 | 西安交通大学 | 一种电主轴铣削颤振鲁棒主动控制方法及其主动控制*** |
CN108846242A (zh) * | 2018-07-09 | 2018-11-20 | 西北工业大学 | 基于预拉应力施加的薄壁件铣削颤振抑制方法 |
CN109048466A (zh) * | 2018-09-03 | 2018-12-21 | 西安交通大学 | 一种基于多频变转速的铣削颤振抑制方法 |
Non-Patent Citations (1)
Title |
---|
张之敬: "基于再生理论的微小型正交车铣颤振", 《清华大学学报(自然科学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020143203A1 (zh) * | 2019-01-07 | 2020-07-16 | 北京理工大学 | 一种深孔镗削加工颤振的在线监测与抑制方法 |
CN110346130A (zh) * | 2019-07-19 | 2019-10-18 | 北京理工大学 | 一种基于经验模态分解和时频多特征的镗削颤振检测方法 |
CN110434676A (zh) * | 2019-07-29 | 2019-11-12 | 北京理工大学 | 一种多传感器时频特征融合的镗削颤振监测方法 |
CN110712066A (zh) * | 2019-10-22 | 2020-01-21 | 湖南工学院 | 适用于深孔内螺纹加工中刀具状态的监测方法 |
CN112180833A (zh) * | 2020-09-07 | 2021-01-05 | 清华大学 | 数控机床及其自适应加工的控制方法和装置 |
CN112180833B (zh) * | 2020-09-07 | 2021-09-07 | 清华大学 | 数控机床及其自适应加工的控制方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
WO2020143203A1 (zh) | 2020-07-16 |
CN109746762B (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109746762A (zh) | 一种深孔镗削加工颤振的在线监测与抑制方法 | |
Zhu et al. | Recent progress of chatter prediction, detection and suppression in milling | |
Quintana et al. | Chatter in machining processes: A review | |
Siddhpura et al. | A review of chatter vibration research in turning | |
Tao et al. | Timely chatter identification for robotic drilling using a local maximum synchrosqueezing-based method | |
Liu et al. | Early chatter detection in gear grinding process using servo feed motor current | |
Wu et al. | Chatter suppression in fast tool servo-assisted turning by spindle speed variation | |
CN111618658B (zh) | 一种面向无颤振高效铣削的主轴转速自适应调整方法 | |
CN105965320A (zh) | 一种高速铣削电主轴颤振智能检测与主动抑制装置 | |
Gagnol et al. | Stability-based spindle design optimization | |
CN102195545B (zh) | 电动机的控制装置 | |
CN102853979A (zh) | 电主轴半主动振动控制试验台架、***及电主轴控制方法 | |
CN106644475A (zh) | 一种微型气体轴承静动态特性及稳定性测试方法 | |
CN108062452A (zh) | 一种弧形齿面蜗杆减速机动态性能评价及优化方法 | |
Bayly et al. | Analysis of tool oscillation and hole roundness error in a quasi-static model of reaming | |
Guo et al. | Comprehensive effect of multi-parameters on vibration in high-speed precision milling | |
Zhai et al. | Influence of cutting parameters on force coefficients and stability in plunge milling | |
CN111230590A (zh) | 拉削加工稳定性监测方法及*** | |
CN112977394B (zh) | 一种抑制发动机扭矩脉动的方法和混合动力*** | |
Bouchareb et al. | Effect of the interaction between depth of cut and height-to-width ratio of a workpiece on vibration amplitude during face milling of C45 steel | |
CN109048466B (zh) | 一种基于多频变转速的铣削颤振抑制方法 | |
Li et al. | The effect of tool rotation on regenerative chatter in line boring | |
Daghini et al. | Design, implementation and analysis of composite material dampers for turning operations | |
CN109702554A (zh) | 一种自适应振动塑形铣削颤振抑制方法 | |
Ding et al. | Adaptive removal of time-varying harmonics for chatter detection in thin-walled turning |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210202 Termination date: 20220107 |