[发明专利]一种制动盘固有频率在线检测方法

摘要

本发明公开一种制动盘固有频率在线检测方法。采用截断奇异值分解降噪技术提高对噪声的抑制能力。对算法的算例仿真结果表明，基于奇异值截断的LSCF算法具有优异的抗噪声干扰能力，识别结果准确，稳定性高。在某型号制动盘生产线上的试验应用结果表明系统能准确判断制动盘产品质量，且具有稳定、高效的优点。

主权项

1.一种制动盘固有频率在线检测方法，其特征在于具有高的识别精度和稳定性，包括以下步骤：1)使用力锤敲击制动盘；2)通过PULSE采集制动盘的力锤激励力信号和粘贴在制动盘上传感器的振动加速度响应信号，根据采集得到的信号得到各响应点和激励点之间的频率响应函数；3)对频率响应函数进行截断奇异值分解降噪：3‑1)将测得的频率响应函数转化为时域的单位脉冲响应函数；3‑2)将转化后的单位脉冲响应函数看作一个含有噪声干扰的原始信号h(tk)；其中：tk为第k个时间点，Nf为频率响应函数的线数；根据Takens相空间重构理论，将所述的原始信号映射到m×n维相空间内，m＜n，得到重构的相空间轨道矩阵DH：相空间轨道矩阵DH为Hankel矩阵，且m+n‑1＝Nf；相空间轨道矩阵DH表示为DH＝D+W，其中D为信号的轨道矩阵，W为噪声的轨道矩阵；3‑3)对相空间轨道矩阵DH进行奇异值分解：svd(DH)＝[U，S，V]   (2)U为m×m阶酉矩阵，V为n×n阶酉矩阵，S是半正定m×n阶对角矩阵，DH的奇异值从大到小排列在S的对角线上，根据Frobenious范数意义下矩阵最佳逼近定理，截断其前p个奇异值而其他奇异值置为零，利用奇异值分解的逆过程得到D′H矩阵，即D′H＝U×Sp×VH，Sp为截断后部分奇异值置零的m×n阶矩阵；D′H矩阵即为截断阶次为p的情况下对相空间轨道矩阵DH的最佳逼近；3‑4)截断阶次p可由信号的奇异熵增量确定，截断后信号奇异熵的表达式为：式中，ΔEi为奇异熵在阶次i处的增量，其值可通过下式计算：式中，λi为从大到小排列的第i个较大奇异值；通过控制ΔEi的大小来确定截断阶次，在保留信号有效信息的同时去除噪声；3‑5)将m×n维矩阵D′_H向右下补零成为N_f×N_f维矩阵则得到降噪后信号h′(t_k)，4)将降噪后的时域信号转化至频域后采用LSCF算法进行参数识别4‑1)频率响应函数的同分母模型可表示为：其中，H(ωk)是频率响应函数，k为序列号；k＝1，2，...，Nf；ωk为第k个频率点的频率，ωk＝kΔω，Δω为频率分辨率，br和ar分别是分子和分母多项式系数，Zk是多项式基函数，2N是多项式的阶次；频率响应函数的模态阶数为N，r为多项式的阶次；采用离散时间域模型，设离散时间间隔为Δt，则Δt＝2π/(NfΔω)   (7)Zk可表示为：式中，j为复数符号；4‑2)将实测频率响应函数替代频率响应函数H(ω_k)，基于公式(6)构造误差函数e(ω_k)如下：4‑3)构造相应的最小二乘方程组如下：Jθ＝[X Y]θ＝0   (10)式中J＝[X Y]且：θb＝[‑b0 ‑b1 … ‑b2N]T   (14)θa＝[a0 a1 … a2N]T   (15)式中，aq为分母多项式系数，bq为分子多项式系数，q＝0、1……2N；将等式(10)两端左乘JH并正规化方程如下：式中，Re表示取复数的实部，R＝Re(XHX)∈R(2N+1)×(2N+1)，S＝Re(XHY)∈R(2N+1)×(2N+1)，T＝Re(YHY)∈R(2N+1)×(2N+1)；由式(16)第一行可使分子多项式系数θb由分母多项式系数θa表示为：θb＝‑R‑1Sθa   (17)将θb的表达式代入式(16)第二行可得缩减的正规化方程如下：(T‑STR‑1S)θa＝0   (18)取a2N＝1代入式(18)，即可求解出分母多项式系数θa；4‑4)在θa确定以后，即可求解频率响应函数H(ωk)的分母多项式对应的方程4‑5)对矩阵A进行特征值分解如下：式中，Φ为特征向量矩阵，Λ为特征值对角矩阵；分母多项式对应方程的根位于特征值矩阵Λ的对角线上，且Λ_i是第i个特征值，λ_i即为系统的第i个极点；从λi的实部中即可得到实测频率响应函数的各阶固有频率；从λi的虚部中即可得到实测频率响应函数的各阶阻尼比。