主权项 |
一种外圆切入式磨削工件圆度误差实时计算方法,其特征在于,包括以下步骤:1)磨削全过程的声发射信号采集为了准确采集声发射信号,将声发射传感器安装在外圆磨床尾架上,通过数据采集系统,由计算机采集获得,并对所采集的声发射信号进行滤波处理,采用切比雪夫II低通滤波器滤除干扰信号,磨削全过程声发射信号采集完成后保存好数据;2)提取声发射信号谐波频率及对应的振幅根据步骤1)获得的声发射信号,通过式n=fT,式中n为磨削工件一圈时信号采集点个数,f为声发射传感器的采集频率,T为工件的转动周期,将磨削全过程的声发射信号按每圈进行分组,然后对每圈声发射信号进行傅里叶变换,提取每圈声发射信号的谐波频率及对应的振幅;3)建立工件轮廓水平展开后的轮廓轨迹方程由工件轮廓按轮廓上某一点截开水平展开后的轮廓形貌类似于谐波信号,从而建立模型:砂轮对工件每圈的磨削效果是多个正弦波形的叠加;工件在连续两圈的磨削过程中正弦波形存在一个滞后相位<img file="FDA0001091476110000011.GIF" wi="59" he="49" />其表达式为:<img file="FDA0001091476110000012.GIF" wi="1170" he="115" />上式中:f<sub>w</sub>为工件的转动频率,f<sub>s</sub>为砂轮的转动频率,δ为f<sub>s</sub>/f<sub>w</sub>数值上的小数部分;工件被磨削第j圈时的磨削效果理论公式为:<img file="FDA0001091476110000013.GIF" wi="1689" he="87" />上式中:m为工件第j圈磨削时通过声发射信号分析得到的谐波频率个数,j为磨削的圈数,当j=1时,第一圈磨削时的组成波的正弦波初相位都为0,f<sub>ij</sub>为磨削第j圈时通过声发射信号提取的第i个谐波频率,a<sub>ij</sub>为磨削第j圈时第i个谐波频率所对应的振幅,t为工件被磨削一周过程中的时间,T为工件的转动周期;由工件最终的轮廓是各圈累计加工效果叠加的总和,并考虑到磨削过程中存在噪声信号的干扰,可得工件水平展开后轮廓轨迹方程为:<img file="FDA0001091476110000021.GIF" wi="1646" he="177" />上式中:n表示工件被磨削的圈数,randn(t)为噪声信号函数。根据步骤2)提取的每圈声发射信号的谐波频率及对应的振幅带入式(3)可计算出工件磨削完成时其轮廓按水平展开后的轮廓计算表达式;4)工件圆度误差计算采用最小二乘圆法计算工件圆度误差,O为工件基圆圆心,O<sub>1</sub>为最小二乘圆圆心,点A为最小二乘圆上任意一点,以O为直角坐标系原点,O与截开点所在直线为X轴,OO<sub>1</sub>与X轴的夹角为ε,OA与X轴的夹角为θ<sub>i</sub>;由于|OO<sub>1</sub>|远小于工件基圆半径,由几何关系近似得:|OA|=r<sub>i</sub>≈r+acos(θ<sub>i</sub>)+bsin(θ<sub>i</sub>) (4)上式中:a为最小二乘圆圆心横坐标,b为最小二乘圆圆心纵坐标,r为最小二乘圆半径,θ<sub>i</sub>为OA与X轴的夹角;为了计算水平展开后轮廓的圆度误差,将圆度误差计算模型水平展开,以工件基圆为基准,对基圆上沿截开点展开后水平拉直知工件基圆轨迹是一条水平线,建立直角坐标系,以工件基圆半径轨迹所在的直线为t轴,截开点为直角坐标系原点;由水平展开后的几何关系得最小二乘圆轨迹方程为:h<sub>i</sub>=r<sub>i</sub>‑r<sub>0</sub>=r+acos(θ<sub>i</sub>)+bsin(θ<sub>i</sub>)‑r<sub>0</sub> (5)上式中,r<sub>0</sub>为工件基圆半径。θ<sub>i</sub>=2πf<sub>w</sub>t (0≤t≤T) (6)联立式(5)、(6)可得:h<sub>i</sub>=r+acos(2πf<sub>w</sub>t)+bsin(2πf<sub>w</sub>t)‑r<sub>0</sub> (0≤t≤T) (7)根据最小二乘圆圆度误差的计算原理则有:<maths num="0001"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>∫</mo><mn>0</mn><mi>T</mi></msubsup><msup><mrow><mo>(</mo><mi>Ω</mi><mo>-</mo><msub><mi>h</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001091476110000031.GIF" wi="1174" he="102" /></maths>求令f(a,b,r)取得最小值的最小二乘圆圆心坐标(a,b)及最小二乘圆半径r,由式(7)可知最小二乘圆轨迹与基圆半径轨迹的关系,从而最小二乘圆圆心轨迹的表达式为:∈<sub>i</sub>=acos(2πf<sub>w</sub>t)+bsin(2πf<sub>w</sub>t)‑r<sub>0</sub> (9)在(a,b)及r已知的情况下由式(3)、式(9)可得工件轮廓轨迹上各点到最小二乘圆圆心轨迹的距离的表达式为:<img file="FDA0001091476110000032.GIF" wi="1738" he="82" />h=f<sub>max</sub>(t)‑f<sub>min</sub>(t) (11)求出式(10)的最大值和最小值并带入式(11)求出h,h值即为工件的圆度误差。 |