1. 引言
我们目前在机床装配过程中常见的测量校正仪器有:百分/千分杠杆表、大理石方尺、激光干涉仪、水平仪、芯棒等,这些仪器虽然测量的精度能得到保证,但其价格昂贵且耗时费力,检测结果没有得到量化,装配质量的评判仍需要依靠工作人员的经验,机床装配难免会出现偏差。虽然目前一些数控系统针对机床装配问题开发了一些自诊断功能,但还局限于电气系统方面的故障,对于机床加工过程中产生的机械故障问题还没有有效的解决方案,且由于振动信息耦合性和机床机械结构的封闭特性使得借助外置传感器的方法有一定的局限性,反而让机床装配过程变得更加复杂 [1]。
通过数控机床的机械结构和工作原理可知,其进给轴传动系统是利用伺服电机连接驱动控制丝杠、齿轮箱、联轴器等部件的运转,如图1所示,从而带动工作台、刀具或加工物料的移动 [2]。这些部件的运行状态会直接影响伺服电机转子的运行,从而间接影响了伺服电机定子电流信号的变化。最初,国内外学者对电机故障和其负载电流的关联性进行了建模和分析。韩国浦项大学的Jeong和Cho [3] 利用建模探索出电流信号和机床倾斜角之间的关系;G.W. Younkin [4] 也提出利用该关联性矫正机床水平相较于使用水平仪更加精准快捷。高宏力等 [5] 采用恒模算法对负载电流的频域特性进行提取,延长了机床使用寿命。王淞灵 [6] 提出基于尺寸变动度的装配特征精度指标评价体系,对机床装配质量进行量化。张会端 [7] 利用拉格朗日方程推导进给轴系统各个部件的运动方程,进而分析出其稳定性及动态特征。因此,电流信号很适合用于机床传动部件的故障诊断,对电流信号的变化规律进行数据分析和提取,就可以有效准确地监测诊断传动系统的故障,从而达到评判机床装配质量的目的。
本文将介绍一种基于进给轴信号采样的检测方法,利用HNC-8采样模块对进给轴位置、电流和速度数据进行采样,分析伺服电机在不同进给速度下电流的波动和频谱特征,筛选适合分析装配质量的速度区间,再结合曲线拟合分析出常见的进给轴故障。并在HNC-8数控系统上开发了相应的诊断模块,以实现对装配质量的量化,最后通过现场实验进行了验证。由于机床机械结构复杂,关于机床装配质量的分析是一个很庞大的系统工程,所以在本文仅介绍进给轴系统的装配质量检测方法,但其它部件的装配质量问题也可以参考该方法来诊断。
Figure 1. The composition of feed shaft system
图1. 进给轴系统基本组成
2. 基于频域分析的速度区间筛选实验
2.1. 实验原理
由于电机齿槽效应 [8] 的影响,当转矩频率与电机定子或转子的机械共振频率一致时,齿槽转矩产生的振动和噪声将放大,与此同时,当进给轴运动处于低速爬行状态时,其产生的特征频率与调制波频率极其接近,因而很难辨别因机械故障产生的特征频率,所以需要筛选出适合分析装配质量的速度区间,并对采样数据进行滤波,尽可能抑制或防止其他频谱成份的干扰。
2.2. 实验过程
在三台机床上,每台机床的X、Y、Z轴分别以10 mm/min、1000 mm/min、4000 mm/min、8000 mm/min四种速度做全行程正反向运动,各运行三次,记录速度电流相应的特征频率,其中三台机床的各进给轴电机型号和规格都相同,共计测得108组数据,在此列举一台机床实验数据进行说明。
2.3. 实验对象配置参数
表1为该实验中一台机床的X/Y/Z进给轴配置参数。
Table 1. Feed shaft configuration parameters of the experimental machine tool
表1. 实验机床进给轴配置参数
2.4. 实验结果
表2为该实验中一台机床的X/Y/Z进给轴在不同速度区间下运转产生的特征频率。
Table 2. Characteristic frequency of feed shaft being tested
表2. 测试进给轴特征频率
2.5. 实验分析
· 在10 mm/min的低速爬行状态下,X、Y、Z进给轴电流频率均以极对数频率0.05 Hz为主,但实际频谱未检测到;爬行速度受噪声频率影响太大,且结构频率太低不容易检测到。
· 在1000 mm/min的低中速区间下,X、Y轴电流特征频率:5.4 Hz,速度特征频率:5.4 Hz、12.7 Hz;Z轴电流频谱特征:16.6 Hz,速度频谱特征:16.6 Hz、25.1 Hz。X、Y、Z进给轴频谱成份均以结构频率为主,以及少量的4~6倍频,受一定的杂波影响。
· 在4000 mm/min的中速区间下,X、Y轴频谱特征以20.8 HZ为主,Z轴以16.67 HZ为主。X、Y、Z进给轴频谱成份均以结构频率为主,受杂波影响较小。
· 在8000 mm/min的高速区间下,X、Y轴电流频谱特征:41.5 Hz,速度频谱特征:41.5 Hz、100.1 Hz。X、Y、Z进给轴频谱成份均以结构频率为主,以及少量的2~3倍频,受一定的杂波影响。
2.6. 实验结论
· 进给轴在低速爬行时,电流波动极对数频谱较明显,速度波动杂波较严重,频谱特征不容易提取。
· 低速频谱以高阶成份为主,随着转频增大,逐渐转为低阶,以结构频率为主。
· 进给轴在中速区间的速度电流均很平稳,以结构频率为主,受杂波影响较小。适合进一步分析故障特征频率。
3. 装配质量量化方法
3.1. 时域波动指标的应用
时域分析指分析采样数据的波动特征,计算出一系列量化指标,如电流/速度的负载比、均方差、波动率等,来评判机床运行过程中是否平稳。
3.1.1. 计算方法
· 电流负载比,电流平均值与伺服电机额定电流的比值。判断伺服电机在恒速段下是否超负载运转,计算公式如下,其中,Y为电流负载比,
电流/速度平均值,In为额定电流。
· 电流/速度均方差(标准差),描述进给轴在不同位置下负载电流/实际速度的离散程度。该值越小,说明进给运动越平稳,反之亦然。计算公式如下,其中,
为电流/速度均方差(标准差), 表示电流/速度平均值,xi表示某一位置下的负载电流/实际速度值,n为数据量。
· 电流/速度波动率,表示进给轴在不同位置下负载电流/实际速度的变化幅度,即如果进给轴在某一段短距离的行程中,负载电流/实际速度产生了较大的变化,那么该值就会变大。计算公式如下,其中,Y为电流/速度波动率,
为电流/速度均方差(标准差),In为额定电流。
3.1.2. 波动指标横向对比实验
仍使用速度区间筛选实验中的三台机床,其编号分别为5、8、9号,每个机床X、Y、Z轴分别以相同的进给速度做正反向运动,走完全行程一个来回为一次,每次设置的进给速度逐渐递增或递减,各进行若干次,对比三台机床进给轴随速度变化下的波动特征。这里仅对三台机床进给轴的电流均方差进行对比,以下图2~4分别为三台机床X、Y、Z轴电流均方差对比,图中所示的横坐标为速度系数,与进给速度正相关;纵坐标为电流均方差值。
Figure 2. Mean square error of the current on X-axis
图2. X轴电流均方差
Figure 3. Mean square error of the current on Y-axis
图3. Y轴电流均方差
Figure 4. Mean square error of the current on Z-axis
图4. Z轴电流均方差
3.1.3. 分析结果
· X、Y、Z三轴的电流均方差随着进给速度的增长有下降的趋势,速度波动趋于稳定。
· Y轴的波动9号机床最小,X与Z轴的波动8号机床最小。从理论上来说,8号机床的X、Z轴装配质量相较于其它两台机床的好;9号机床的Y轴装配质量相较于其他两台机床的好。
3.2. 时域曲线的应用
3.2.1. HNC-8采样模块与数据可视化
HNC-8具有独立的采样模块,使用数字量内部采集把生成的电控信号转化为实时性较高的数据如进给速度、实际位置、负载电流、主轴转速等,然后将转化的数据存放在配置好的采样通道内,供界面分析或获取,支持多通道,可最大同时配置32个采样通道,每个采样通道都可以储存不同的数据,最大可储存10万个数据点。在此基础上绘制进给轴全行程的负载电流、进给速度时域曲线。其中,实际位置为横坐标,负载电流/进给速度为纵坐标。
采样开始之前可以根据配置和进给轴的正负软限位自动生成G代码,因为需要分析进给轴在整条丝杠不同位置下的负载电流和实际进给速度的动态变化,则设定采样过程中的进给轴运动轨迹为全行程,即从负软限位到正软限位,再从正软限位返回到负软限位。进给轴运动由负软限位到正软限位为正向,相反,由正软限位到负软限位为负向。
3.2.2. 案例分析
客户现场的一台机床在加工工件存在接刀痕现象,为查明原因,使该机床的X、Y、Z进给轴以中速5000 mm/min做全行程的正反向运动,对该进给轴的负载电流、实际速度数据进行采样,绘制时域曲线,发现X/Y轴的速度均方差与其它相同型号配置的机床的基本一致,而Z轴较大,其时域曲线如图5(a)所示。可以看到该机床的Z轴在调整前,全程的速度波动不均匀,从实际位置的−500 mm~−300 mm、0 mm~−150 mm都有明显的速度波动,波动幅度在±100 mm/min左右。
经多次实验发现,Z轴在每次换向的过程中,速度有较大波动。图5(a)中的曲线也证明了这一点,所以优先从换向问题的角度查找原因。常见原因有四项:① 接线端子接触不良;② 螺钉松动;③ 伺服参数增益过大;④ 反向间隙过大。这里使用排除法来缩小确定范围。对比其它相同型号机床的伺服参数设置,发现该机床与其它机床一致,排除②;分别检查接线和螺钉松紧,排除① ②,所以猜测是反向间隙误差引起。在分析了导致反向间隙误差的原因 [9] 之后,发现Z轴的丝杠预紧力不适当,导致轴向窜动量超出允许范围。在调整之后,Z轴的时域曲线如图5(b)所示,全程的速度波动均匀,波动幅度在±25 mm/min左右,消除了接刀痕现象。
4. 丝杠装配质量故障诊断
经3.2.2的案例之后发现,使用时域曲线和计算波动指标的方法只能知道哪个进给轴工作有异常,但不能知道异常产生的根本原因,故障诊断仍需要靠调试或装配人员的经验,例如丝杠的装配。为了能够进一步诊断装配质量故障,在时域曲线的应用基础上提出丝杠装配质量故障诊断方法。
(a) 调整前(b) 调整后
Figure 5. The comparisons of time-domain curves (a-b)
图5. 时域曲线对比(a-b)
4.1. 原理
关于丝杠装配,常见的问题有:① 丝杠与导轨不平行、② 丝杠预紧力过大、③ 转子与丝杠不同轴等等,如图6所示 [10] [11]。
由于电流信号可以很好地反映进给轴传动系统运行状态,那么其时域曲线可对应丝杠的整体结构。绘制的时域曲线应是全行程的,能够保证进给运动从丝杠的一端到另一端。使用最小二乘法实现二次多项式拟合生成如下拟合曲线方程:
其中c为平均值,如果机床运行状态趋于稳定,系数a,b趋近于0,拟合曲线近似于一条直线如图7(a)所示;如果机床运行状态不稳定,则系数a,b的绝对值较大,当系数a不为0时,拟合曲线近似于一条抛物线如图7(b)所示,说明丝杠预紧力可能过大或过小或者丝杠机械结构产生了弯曲;当系数a趋近于0,系数b不为0时,拟合曲线近似于一条斜线如图7(c)所示,说明丝杠与导轨有可能不平行。
(a) 丝杠与导航不平行 (b) 电机转子与丝杠不同轴
Figure 6. The screw assembly quality problems (a-b)
图6. 丝杠装配问题(a-b)
(a) 系数a,b趋近于0(b) 系数a不为0(c) 系数a趋近于0,系数b不为0
Figure 7. Curve-fitting (a-c)
图7. 曲线拟合(a-c)
4.2. 实验验证
在客户现场选取一台机床,对该机床的单个轴丝杠分别进行预紧力过大、丝杠导轨不平行实验,为了便于实验,这里对X轴进行测试。使X轴以4000 mm/min的速度做全行程的正反向运动,观察X轴时域曲线拟合方程三个系数的值。
4.2.1. 预紧力过大实验
分别对X轴丝杠预紧力调整至45˚、90˚和不调整(正常)三种情况进行对比,其曲线拟合实验数据如表3所示。
Table 3. The experimental data of excessive preload
表3. 预紧力过大实验数据
由表3可知,调整丝杠预紧力会显著影响曲线拟合方程系数a的变化,丝杠预紧力越大,曲线拟合方程系数a系数绝对值越大;在丝杠预紧力不调整(正常)情况下,其系数a,b的绝对值趋近于0;负载电流平均值为1.117 A与系数c接近;调整丝杠预紧力对系数b产生的影响极小,丝杠预紧力与系数b应呈无关性。
4.2.2. 丝杠导轨不平行实验
分别对X轴丝杠相对导轨朝上下左右四个方向进行偏移,相对偏移量为0.2 mm,对比其时域拟合曲线方程系数,其曲线拟合实验数据如表4所示。
Table 4. The experimental data of the lead screw is not parallel to the guide rail
表4. 丝杠导轨不平行实验数据
由表4可知,调整丝杠相对导轨偏移量会显著影响曲线拟合方程系数b的变化;在丝杠预紧力不调整(正常)情况下,其系数a,b的绝对值趋近于0;负载电流平均值为1.117 A与系数c接近;调整丝杠预紧力对系数a产生的影响极小,丝杠与导轨不平行度与系数a应呈无关性。
4.2.3. 实验结论
· 拟合方程系数a绝对值大小与丝杠预紧力调整不当问题有关。
· 拟合方程系数b绝对值大小与丝杠与导轨不平行问题有关。
5. 结语
基于进给轴信号采样的装配质量检测方法提供了更直观的数据分析,从时域、频域、曲线拟合等多个维度评判进给轴电机和丝杠的装配质量,相较于传统仪器矫正方法方便灵活且不易出错,无需人工手动操作,为机床厂家节省了成本,提高了生产效率。该方法在经过机床实验验证之后,也确保了其可靠性和实用性。
基金项目
山东省重点研发计划创新工程项目(项目编号:2019JZZY010437)。
NOTES
*通讯作者。