外圆切入式磨削工件圆度误差实时计算方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106181771A(43)申请公布日2016.12.07(21)申请号201610724244.3(22)申请日2016.08.25(71)申请人上海理工大学地址200093上海市杨浦区军工路516号(72)发明人姜晨程金义吴涛郝宇王鹏陆葳坪(74)专利代理机构上海申汇专利代理有限公司31001代理人吴宝根(51)Int.Cl.B24B49/02(2006.01)权利要求书2页说明书5页附图2页(54)发明名称外圆切入式磨削工件圆度误差实时计算方法(57)摘要本发明涉及一种外圆切入式磨削工件圆度误差实时计算方法，首先，采集磨削全过程中的声发射信号，并对该信号进行滤波处理；接着通过傅里叶变换提取所获得信号的谐波频率及对应的振幅；然后通过工件最终轮廓是砂轮对工件每圈叠加加工效果的总和，并结合谐波频率及振幅建立工件轮廓水平展开后的轮廓轨迹方程；最后应用经过变形处理的最小二乘法来计算工件圆度误差。该计算方法对外圆切入磨削过程中工件圆度误差进行实时计算，克服了现有外圆磨削对圆度仪的依赖，提高磨削加工效率，节约生产成本。CN106181771ACN106181771A权利要求书1/2页1.一种外圆切入式磨削工件圆度误差实时计算方法，其特征在于，包括以下步骤：1)磨削全过程的声发射信号采集为了准确采集声发射信号，将声发射传感器安装在外圆磨床尾架上，通过数据采集系统，由计算机采集获得，并对所采集的声发射信号进行滤波处理，采用切比雪夫II低通滤波器滤除干扰信号，磨削全过程声发射信号采集完成后保存好数据；2)提取声发射信号谐波频率及对应的振幅根据步骤1)获得的声发射信号，通过式n＝fT，式中n为磨削工件一圈时信号采集点个数，f为声发射传感器的采集频率，T为工件的转动周期，将磨削全过程的声发射信号按每圈进行分组，然后对每圈声发射信号进行傅里叶变换，提取每圈声发射信号的谐波频率及对应的振幅；3)建立工件轮廓水平展开后的轮廓轨迹方程由工件轮廓按轮廓上某一点截开水平展开后的轮廓形貌类似于谐波信号，从而建立模型：砂轮对工件每圈的磨削效果是多个正弦波形的叠加；工件在连续两圈的磨削过程中正弦波形存在一个滞后相位其表达式为：上式中：fw为工件的转动频率，fs为砂轮的转动频率，δ为fs/fw数值上的小数部分；工件被磨削第j圈时的磨削效果理论公式为：上式中：m为工件第j圈磨削时通过声发射信号分析得到的谐波频率个数，j为磨削的圈数，当j＝1时，第一圈磨削时的组成波的正弦波初相位都为0，fij为磨削第j圈时通过声发射信号提取的第i个谐波频率，aij为磨削第j圈时第i个谐波频率所对应的振幅，t为工件被磨削一周过程中的时间，T为工件的转动周期；由工件最终的轮廓是各圈累计加工效果叠加的总和，并考虑到磨削过程中存在噪声信号的干扰，可得工件水平展开后轮廓轨迹方程为：上式中：n表示工件被磨削的圈数，randn(t)为噪声信号函数。根据步骤2)提取的每圈声发射信号的谐波频率及对应的振幅带入式(3)可计算出工件磨削完成时其轮廓按水平展开后的轮廓计算表达式；4)工件圆度误差计算采用最小二乘圆法计算工件圆度误差，O为工件基圆圆心，O1为最小二乘圆圆心，点A为最小二乘圆上任意一点，以O为直角坐标系原点，O与截开点所在直线为X轴，OO1与X轴的夹角为ε，OA与X轴的夹角为θi；由于|OO1|远小于工件基圆半径，由几何关系近似得：|OA|＝ri≈r+acos(θi)+bsin(θi)(4)上式中：a为最小二乘圆圆心横坐标，b为最小二乘圆圆心纵坐标，r为最小二乘圆半径，θi为OA与X轴的夹角；2CN106181771A权利要求书2/2页为了计算水平展开后轮廓的圆度误差，将圆度误差计算模型水平展开，以工件基圆为基准，对基圆上沿截开点展开后水平拉直知工件基圆轨迹是一条水平线，建立直角坐标系，以工件基圆半径轨迹所在的直线为t轴，截开点为直角坐标系原点；由水平展开后的几何关系得最小二乘圆轨迹方程为：hi＝ri-r0＝r+acos(θi)+bsin(θi)-r0(5)上式中，r0为工件基圆半径。θi＝2πfwt(0≤t≤T)(6)联立式(5)、(6)可得：hi＝r+acos(2πfwt)+bsin(2πfwt)-r0(0≤t≤T)(7)根据最小二乘圆圆度误差的计算原理则有：求令f(a，b，r)取得最小值的最小二乘圆圆心坐标(a，b)及最小二乘圆半径r，由式(7)可知最小二乘圆轨迹与基圆半径轨迹的关系，从而最小二乘圆圆心轨迹的表达式为：∈i＝acos(2πfwt)+bsin(2πfwt)-r0(9)在(a，b)及r已知的情况下由式(3)、式(9)可得工件轮廓轨迹上各点到最小二乘圆圆心轨迹的距离的表达式为：h＝fmax(t)-fmi