发明名称 一种汽车变速器智能化在线质量检测方法
摘要 一种汽车变速器智能化在线质量检测方法属于汽车噪声和振动控制领域。针对汽车变速器在线交检试验的重要性及复杂性,本文在信号采集方面实现了振动电压信号与转速脉冲信号的同步采集。由于在汽车变速器的交检试验过程采集到的振动电压信号为非稳态信号,用常规的信号处理方法会出现“频率混叠”现象,本发明将阶次方法与倒频谱方法相结合进行变速器非稳态信号的处理进而提取特征信号,既解决了“频率混叠”问题又消除了功率谱的边频带,更便于观察与提取特征信号。在质量检测方面使用了支持向量机方法,同时利用标准极值曲线方法不断丰富与更新支持向量机模型,使该模型包含更多的信息,达到了检测速度快且准确的目的。
申请公布号 CN102539151B 申请公布日期 2014.05.14
申请号 CN201210016832.3 申请日期 2012.01.18
申请人 北京工业大学 发明人 付胜;徐斌;蔡莎莎
分类号 G01M13/02(2006.01)I 主分类号 G01M13/02(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 刘萍
主权项 1.一种汽车变速器智能化在线质量检测方法,其特征在于,包括以下步骤:(1)将汽车变速器安装在测试试验台上,将编码器安装在试验台的输入轴上,将振动传感器水平安装在汽车变速器上;当试验台带动汽车变速器旋转时,开始采集转速信号;(2)设编码器每转发出w个脉冲,即输入轴每转过360/w度角,编码器就发出一个脉冲,计数器的值就增加1;每隔Δt的时间读取一次计数器上的脉冲数,则假设在时刻t<sub>i</sub>读取计数器上的脉冲数为m<sub>i</sub>,t<sub>i+1</sub>时刻读取计数器上的脉冲数为m<sub>i+1</sub>,则当前的转速值<img file="FDA0000462062020000011.GIF" wi="315" he="115" />(3)假设系统测试转速范围为n<sub>0</sub>~n<sub>t</sub>,n<sub>0</sub>为转速下限,为最高转速的30%,n<sub>t</sub>为转速的上限,为最高转速的80%;当转速n<sub>0</sub>≤n<sub>i</sub>≤n<sub>t</sub>时,将计数器清零,开始采集振动电压信号,同时重新采集转速脉冲信号;设时域平均法中的平均次数为k,每次分析振动电压信号点数为x<sub>0</sub>,分析转数为r;则当脉冲数达到n<sub>总</sub>=w·k·r时,将采集的振动电压信号与转速脉冲存储于文件中,文件名为当前时间,然后继续进行同步采集振动电压信号与转速脉冲信号;设振动电压信号的采集频率为f<sub>1</sub>,脉冲信号的读取频率为f<sub>2</sub>;(4)信号处理与信号采集同时进行,信号采集模块不停的检测是否存在新的数据文件存在,若存在则将相同文件名的振动电压信号S(t<sub>i</sub>){S(t<sub>1</sub>),S(t<sub>2</sub>),S(t<sub>3</sub>),……S(t<sub>N1</sub>)}与转速脉冲信号M(p){M(1),M(2),M(3),……M(N2)}提取出来进行处理,其中S(t<sub>i</sub>)为每一个时刻<img file="FDA0000462062020000012.GIF" wi="141" he="139" />对应的电压值,M(p)为每一时刻<img file="FDA0000462062020000013.GIF" wi="152" he="130" />对应的脉冲个数,i为电压值序列号,p脉冲数序列号,N1为振动电压信号文件中电压值的个数,N2为脉冲信号文件中脉冲读取的次数;(5)由于输入轴每转过1度角计数器的值就加1,所以脉冲个数M(p)也是轴转过的角度;将M(p)与对应的时间序列t<sub>p</sub>生成曲线;对信号分析时,采用r转为数据源,即分析转数,每次分析电压个数为x<sub>0</sub>,则每转取的电压个数为<img file="FDA0000462062020000014.GIF" wi="96" he="116" />由于需要等角度间隔取数据,所以每转过<img file="FDA0000462062020000015.GIF" wi="177" he="116" />角度时取一个数据点;通过插值得到每一个数据点所对应的时刻t(g),即等角度时间序列t(g);(6)将振动电压信号S(t<sub>i</sub>)与对应的时间t<sub>i</sub>生成曲线,利用等角度时间序列t(g)将振动电压信号进行拉格朗日插值得到每个等角度时刻t(g)对应的振动电压信号R(θ),其中θ为等间隔转过的角度;(7)对等角度的振动电压信号R(θ)进行时域平均处理,加强周期信号,r为分析转数,k为平均次数;则通过插值得到的等角度的振动电压信号R(θ)为r·k转的电压信号,电压信号个数为x<sub>0</sub>·k;时域平均法即以分析x<sub>0</sub>个振动电压信号为周期截取R(θ)信号,共截取k段,然后将各段对应点相加再进行平均得到输出信号,即:<maths num="0001"><![CDATA[<math><mrow><mi>Q</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>R</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>+</mo><mi>R</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>+</mo><mi>h</mi><mo>)</mo></mrow><mo>+</mo><mi>R</mi><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>x</mi></mrow><mn>0</mn></msub><mo>+</mo><mi>h</mi><mo>)</mo></mrow><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><mi>R</mi><mrow><mo>(</mo><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msub><mi>x</mi><mn>0</mn></msub><mo>+</mo><mi>h</mi><mo>)</mo></mrow></mrow><mi>k</mi></mfrac><mo>,</mo></mrow></math>]]></maths>其中Q(h)为时域平均法后输出的电压值,h为输出的电压值序列号,h=1,2,3…x<sub>0</sub>;(8)通过步骤(7)的Q(h)求取两个时域特征值:均方根幅值、峰值如下所示:均方根幅值(X<sub>1</sub>):<img file="FDA0000462062020000022.GIF" wi="522" he="169" />h=1,2,3…x<sub>0</sub>,x<sub>0</sub>为分析点数;峰值(X<sub>2</sub>):Peak=max|Q(h)|,h=1,2,3…x<sub>0</sub>;(9)将步骤(7)的Q(h)进行倒频谱分析得到谱图C(q<sub>a</sub>),横坐标阶次频率<img file="FDA0000462062020000023.GIF" wi="157" he="117" />纵坐标为电压幅值C(q<sub>a</sub>),其中a=1,2,3…x<sub>0</sub>;r为分析转数;(10)由于变速器类型不同,内部的齿轮轴承参数就不同,进而计算出的特征阶次也不同;所以在数据库中变速器的型号对应特定的特征阶次,根据变速器型号查询数据库得出变速器轴、齿轮、轴承的特定阶次向量:X=[X<sub>i1</sub>,X<sub>i2</sub>,..X<sub>in</sub>,X<sub>j1</sub>,X<sub>j2</sub>,..X<sub>jm</sub>,X<sub>iL1</sub>,X<sub>iL2</sub>,X<sub>iL3</sub>,X<sub>iL4</sub>,X<sub>iR1</sub>,X<sub>iR2</sub>,X<sub>iR3</sub>,X<sub>iR4</sub>],其中下标i为第i根轴特征阶次,共有n根轴;下标j为第j对啮合齿轮特征阶次,共有m对啮合齿轮;下标L1-L4分别为左或前轴承的外圈,内圈,滚动体自转,保持架特征阶次;下标R1-R4分别为右或后轴承的外圈,内圈,滚动体自转,保持架特征阶次;(11)支持向量机模型中选择算法复杂度较小的多项式核函数:K=(X<sup>r</sup>·y+1)<sup>12</sup>,其中X<sup>r</sup>为特征值向量X<sup>r</sup>=[X<sub>1</sub>,X<sub>2</sub>,X<sub>i1</sub>,X<sub>i2</sub>,..X<sub>in</sub>,X<sub>j1</sub>,X<sub>j2</sub>,..X<sub>jm</sub>,X<sub>iL1</sub>,X<sub>iL2</sub>,X<sub>iL3</sub>,X<sub>iL4</sub>,X<sub>iR1</sub>,X<sub>iR2</sub>,X<sub>iR3</sub>,X<sub>iR4</sub>],y为输出向量y=[123…N],N为变速器缺陷类型数目;根据训练样本确定支持向量机模型中缺陷类型,并将时域特征值与频域特征值合并作为特征值向量X<sup>r</sup>=[X<sub>1</sub>,X<sub>2</sub>,X<sub>i1</sub>,X<sub>i2</sub>,..X<sub>in</sub>,X<sub>j1</sub>,X<sub>j2</sub>,..X<sub>jm</sub>,X<sub>iL1</sub>,X<sub>iL2</sub>,X<sub>iL3</sub>,X<sub>iL4</sub>,X<sub>iR1</sub>,X<sub>iR2</sub>,X<sub>iR3</sub>,X<sub>iR4</sub>]输入到支持向量机模型中进行检测,判断当前变速器是否存在缺陷及缺陷类型位置;(12)若步骤(11)的判断结果为存在缺陷,则输出变速器的缺陷类型及位置;若步骤(11)判断没有缺陷,则利用极值曲线方法进行二次判断,将步骤(9)中得到的谱线与数据库中的极值曲线相比较,判断是否超出极值曲线,若没有超出,则判断不存在缺陷;若存在超出部分,则查询数据库得到超出部分所对应的缺陷类型及位置,并且将该缺陷类型及位置增加到上述支持向量机模型中;(13)系统运行过程中,信号采集与信号处理处在两个线程中同时进行,只要变速器转速在测试范围内,信号采集就不断的进行,每个周期结束后,都将信号以时间为文件名存储于硬盘中;而信号处理模块则不停的检测硬盘中是否存在未处理的信号,若存在则重复步骤(4)至步骤(12)。
地址 100124 北京市朝阳区平乐园100号