机床主轴的常见状况和使用.docVIP

数控机床要求较高运动精度，其故障信息的载体与一般机械有所不同。大量诊断实例表明，从加工零件的质量上分析机床的运行状态，是诊断机床的主要思路。主轴是数控机床的关键部件，其运动精度直接影响工件的加工质量。当主轴不平衡运行时，产生的往复惯性力使轴瓦与主轴颈在油膜壁最薄处发生撞击，引起机床主轴的运动误差。主要包括主轴的纯径向跳动D、纯角度摆动和纯轴向窜动L。 　　时间序列分析法是以模式识别理论为基础，将动态过程为随机的复杂系统抽象为简单物理模型以便于实际系统运行状态的监测与诊断。现有文献大部分都是基于一维时间序列进行讨论，但对于许多实际问题，仅考虑一个因素并不能够充分揭示系统内部规律。为进一步解释各相关物理量对系统的影响关系，提高模型拟合精度，可引入多维时间序列模型。本文将多维时间序列应用于工件表面粗糙度Ra监测，分析了D、、L与Ra之间的影响关系。 　　1时序模型的选择常用时序模型有自回归模型（AR）、滑动平均模型（MA）以及自回归滑动模型（ARMA ）。后两类模型都可以用高阶的AR模型来逼近，其逼近程度取决于所取的自回归模型的阶次。由于AR模型的参数估计为线性回归过程，其计算简单、速度快，实际物理系统也往往是全极点系统，因此， AR模型的应用广泛，特别适用于机械故障诊断中。 　　对于AR模型讨论的前提是假设所获得的时间序列是平稳、正态分布的序列。但当机床运转处于某种隐患并且情况在不断恶化时，通过传感器得到的数据将是不平稳时间序列，含有某种随时间稳定发展的趋势。对于非平稳时序经过若干阶差分变换可以得到平稳时间序列，然后对其进行ARMA建模。为此提出非平稳时间序列ARIMA模型。 　　1 1 AR IMA模型ARIMA模型是平稳时间序列在非平稳状态下的拓展，它与非平稳序列{X t }密切相关， AR IMA模型经过有限次差分（ Y t = （ 1- B ）d X t， B为后移算子）后，{ Y t }即为ARMA序列。 　　适合ARIMA模型建模的数据具有缓慢衰变的正样本自相关系数函数的特征；经过差分变换后，适合ARMA建模的数据应具有快速下降的样本自相关系数函数的特征。 　　1 2多维AR模型概念从一维AR模型的概念出发可以定义多维AR模型。多维AR模型可以拟合多数实际情况，因为它可以逼近多维ARMA模型。设有m维的AR（ p）模型：z t = A 1 z t- 1 + A 2 z t- 2 + + A p z t- p + u t（1）式中： z t是m 1维零均值平稳序列；A j（ j = 1， 2，， p）是m m矩阵，其中A p 0； { u t }是m维白噪声序列，在不同时刻彼此无关，即E u t = 0 E u t u T s = S，t= s 0，t s（2）式中：S是m m阶正定阵，且E u t z t- j = 0（ j = 1，2，）。则多变量平稳的条件： det p i= 0 A i B i = 0的根全在单位圆外。如果{ z t }满足式（1）及平稳性条件，则{ z t }就称为m维AR（p ）序列。 　　2多维AR模型的建立2 1数据采集通过传感器得到的振动、转角信号都是连续信号，而建立时序模型需要对连续信号进行离散采样。 　　因此首先要确定采样间隔。若采样间隔不合理，则会产生不同频率谐波的混迭。当取得过大，采样数据将丢失观测信号中原有的相关关系，以致模型发生降阶，降低信号分辨率；当取得过小时，又会将高频噪声作为有用信号计入，使模型发生升阶，造成计算量增大。采用香农采样定理：当采样频率f s高于谐波的最高频率f max的2倍时，即可避免频混，为此常取= 1/2 5f m ax进行采样。 　　2 2多维AR模型的参数估计AR模型的参数估计法常用直接估计法，其中包括YuleW alker估计法和最小二乘估计法。相关矩估计方法是时间序列模型参数估计最常用的方法，它是建立在Y W方程基础上的。虽然该方法的估计精度没有最小二乘法那样高，但其计算简便实用。特别是对于具有正态特性的时间序列，在样本数据足够多的情况下，其估计精度可以与最小二乘估计法相当。本文采用选择Y W估算法。这就得到A j（ j = 1， 2，，p ）的矩估计。从理论上可知YuleW alker矩估计量有很多优点，如正定性、渐近无偏性、弱相容性和渐近正态性。 　　3多维AR模型定阶的FPE准则FPE准则（最小最终预报误差准则）是由样本对模型定阶。它是以模型输出的一步预报误差方差来判定模型阶次：一步预报误差方差阵的行列式越小，就认为模型拟合得越理想，这时的模型阶次认为是最佳的阶次。 　　由样本长度为N的m维观测数据所拟合的AR（p）模型见式（1）。^ T i（8）需要说明，如果FPE p（z t）的值随p的增大而单调上升，则可判定模型为一阶自