-
摩擦力矩检测设备常用PID控制,文中系统采用Z-N-PID算法实现对激光雷达轴系摩擦力矩的测量。Z-N-PID算法其控制参数由人为事先计算确定,因而对实际工况下输入信号特性的随机变化难于自适应[9-10],导致设备重复测量精度>0.0025 N·m,极限误差>0.05 N·m,不满足设备使用需求。而GA-BP神经网络可以通过学习和存贮大量的输入-输出模式映射关系,获得更加精准的网络预测,但因激光雷达轴系具有不同尺寸的俯仰轴系和方位轴系,且在不同转速下工作,所以,在激光雷达轴系摩擦力矩测试设备初期设计使用中,不能获取不同激光雷达轴系所对应的全部PID控制参数,无法对BP神经网络进行完整的训练,导致设备在测试轴系时数据波动较大,且重复测量精度较低。故提出了一种在线学习自更新的控制算法,以适应不同工况下的摩擦力矩测试。
所提出的在线学习自更新控制算法,首先通过GA算法初始化BP神经网络权值,然后分别建立激光雷达轴系摩擦力矩设备测试主轴实际转速、理想转速、转速误差及转速误差变化率的数据云,最后将数据云数据传递至GA-BP算法中,以获取最佳控制参数。自更新控制算法流程图如图1所示。
-
在所提出的自更新控制算法的数据云建立过程中,为避免数据冗余,影响运算速度,使用距离信息与局部密度作为删除数据的标准。
(1) 密度信息
在测试电机实际转速、转速误差、转速误差变化率和理想转速的数据云形成时,使用高斯核函数[11]计算其密度,定义第
$ i $ 个数据的局部密度为:$$ \gamma _k^i{\text{ = }}\sum\limits_1^n {\exp \left( - \frac{{{{\left\| {\left. {x_k^i - x_k^j} \right\|} \right.}^2}}}{{2{\sigma ^2}}}\right)} $$ (1) 式中:
$ {\sigma ^2} $ 为截断距离;$ \left\| {\left. {x_k^i - x_k^j} \right\|} \right. $ 为$ x_k^i $ 到$x_k^j$ 之间的欧式距离。设置采样频率为100 Hz,并通过
$ {x_i} $ 记录样本数据,将之前的局部密度按时间顺序即可得到密度集合$ {\gamma _k} $ :$$ \left\{ {\left. {{\gamma _i}} \right\}} \right._{i = 1}^n = \left\{ {\left. {{\gamma _1},{\gamma _2},{\gamma _3} \cdot \cdot \cdot ,{\gamma _n}} \right\}} \right. $$ (2) (2) 距离信息
取数据
$ {X_i} $ 与比其本身大的局部数据云之间的平均距离作为距离信息:$$ \begin{split} \\ {D_{{x_i}}} = \frac{1}{N}\sum\limits_{{x_i} \in {S_i}}^N {\left\| {\left. {x_{_k}^i - x_k^j} \right\|} \right.} \end{split} $$ (3) 式中:
$ {S_i} $ 表示密度大于$ {X_i} $ 的前$ N $ 个数据云的集合。通过公式(3)对局部密度中除最大值以外所有局部数据云计算距离,最大值处用1补齐,将之前的距离信息按时间排序即可得到距离信息集合:
$$ \left\{ {\left. {{D_i}} \right\}} \right._{i = 1}^n = \left\{ {\left. {{D_1},{D_2},{D_3}, \cdot \cdot \cdot ,{D_n}} \right\}} \right. $$ (4) (3) 数据剔除
在数据云形成过程中,主要剔除两类数据。其一,测试主轴在启动或变速会产生瞬时突变的转速,产生误差较大的数据。其二,局部密度之间的距离可能随数据的累计而变低。为避免数据冗余、提高系统鲁棒性、加快神经网络运算速度,需对后续新数据样本进行判断,其评判标准[12]为:
$$ \alpha {\text{ = }}{{{D_{{x_i}}}} \mathord{\left/ {\vphantom {{{D_{{x_i}}}} {\gamma _{_k}^i}}} \right. } {\gamma _{_k}^i}} $$ (5) 若
$ \alpha $ 接近于0或远大于1,则说明该点局部密度与其他密度值较为接近或偏差较大,对于此类数据应当剔除。 -
为保证BP网络输出精度和速度,选取测试主轴实际转速、转速系统误差、转速误差变化率和理想转速形成的数据云作为输入[13]。BP神经网络采用一层隐含层的网络形式,设置学习速率为0.5。由输入、输出参数的数量设定输入层和输出层节点数分别为4和3,隐含层节点数为5。输出层选取利于前向传输的Sigmoid函数作为激活函数,隐藏层选取以0为中心的tanh函数作为激活函数。
使用GA算法优化BP神经网络初始权值时,首先需对个体
$ i $ 进行编码,编码的长度为:$$ H = b \times c + c \times l + c + l $$ (6) 式中:
$ b $ 表示输入层节点数;$ c $ 表示隐含层节点数;$ l $ 表示输出层节点数。经计算可知,编码长度为43。
而后确定适应度函数,其公式如下:
$$ f = \frac{1}{{\alpha \displaystyle\sum\limits_{i = 1}^l {{{\left( {{r_i} - {y_i}} \right)}^2}} }} $$ (7) 式中:
$ f $ 表示适应度函数;$ {r_i} $ 表示第$ i $ 个节点系统理想输出;$ {y_i} $ 表示第$ i $ 个节点系统预测输出;$ \alpha $ 表示系数。此外,初始化设定种群规模N=100,最大迭代次数Tmax=50,变异率F=0.01,交叉率Cr=0.6。
通过GA算法获得适应度最优的个体,解码后即可获取BP神经网络的连接权值和阈值。
-
为分析文中所提出的自更新控制算法的系统响应时间与抗干扰性能,搭建Z-N-PID控制系统模型和自更新控制算法控制系统模型,对激光雷达轴系摩擦力矩检测设备的测试主轴进行仿真实验。
在离散模型的情况下,增量式数字PID控制算法表示为:
$$ \begin{split} \Delta u\left( k \right) = &{K_p}\left[ {e\left( k \right) - e\left( {k - 1} \right)} \right] + {K_i}e\left( k \right) + {K_d}\\ &\left[ {e\left( k \right) - 2e\left( {k - 1} \right) + e\left( {k - 2} \right)} \right] \end{split} $$ (8) 根据Ziegler-Nichols(Z-N)[14]方法,整定公式如下:
$$ {K_p}{\text{ = }}0.6{K_m},{K_d} = \frac{{{K_p}\pi }}{{4{\omega _m}}},K_i = \frac{{{K_p}{\omega _m}}}{\pi } $$ (9) 式中:
$ {K_m} $ 为系统开始震荡时的$ {K_p} $ 值;$ {\omega _m} $ 为震荡频率。计算初始值为:Kp=1.005、Ki=2.328、Kd=0.119。设置实验初始条件为单位阶跃响应,在10 s处添加0.2单位的瞬时干扰,分别获取两种算法系统响应曲线。系统响应曲线对比图如图2所示,系统响应曲线对比数据见表1。
表 1 系统响应曲线对比数据表
Table 1. Comparison data sheet of system response curve
Item Initial stage Interference phase Overshoot Adjust time/s Overshoot Adjust time/s Z-N-PID control algorithm 7.40% 5.968 17.70% 4.588 Self-updating control algorithm 0 7.197 10.40% 8.77 对表1数据进行分析可知,使用自更新控制算法的控制系统相比于使用Z-N-PID算法的控制系统拥有更低的超调量,但调整时间相对较慢。对于该检测系统而言,较小的超调量意味着系统产生的瞬时冲击较小,即对测量数据的影响也相应地变小,此外也不易导致设备的损坏,延长使用寿命。相对较慢的调整时间并不影响设备最终测量精度,故该自更新控制算法可以更好地满足该设备的需求。
-
为分析自更新控制算法在测量激光雷达轴系动摩擦力矩时,所采集的数据波动大小和重复测量精度,分别使用自更新控制算法的控制系统与Z-N-PID算法的控制系统,通过不含温箱部分的基于传递测量法的激光雷达轴系摩擦力矩检测设备,进行动摩擦力矩测试。激光雷达轴系摩擦力矩实验平台如图3所示。该设备可测的转速范围为0~50 rpm;能够在整周或规定角度范围内检测轴系摩擦力矩;能够检测轴承部件和轴系在负载下的摩擦力矩。
为减少润滑剂、杂质以及其他外界因素的影响[15],在实验前对被测轴系进行清洁。设置实验间隔为10 min,两侧为SKF7207深沟球轴承,测量时间为30 s,采样频率为500 Hz,室温26 ℃。根据参考文献[16],取稳定后(第16~26 s)平均值作为一次测量的结果,进行了100组以上实验,在10、30、50 rpm下各选取12次。对上述两种算法获取的实验数据分别去除最大值与最小值后,首先任取一组实验数据样本分析其结果,所取出的一组动摩擦力矩测量结果如图4所示。
图 4 测量结果。(a) 10 rpm动摩擦力矩测量结果;(b) 30 rpm动摩擦力矩测量结果;(c) 50 rpm动摩擦力矩测量结果
Figure 4. Measurement result. (a) 10 rpm dynamic friction torque measurement results; (b) 30 rpm dynamic friction torque measurement results; (c) 50 rpm dynamic friction torque measurement results
由图4可知,在激光雷达轴系动摩擦力矩检测的过程中,由于润滑剂、杂质等因素的影响,在开始阶段所采集的动摩擦力矩值会产生一定的过冲量,而后逐步趋于稳定。
取图4所展示的数据样本中第16~26 s的数据值,进行数据波动分析,其分析结果如表2所示。
表 2 采集点数据波动比较
Table 2. Comparison of data fluctuations at collection points
Rotating speed Algorithm Average
value/N·mLimit
error/N·mStandard deviation 10 rpm Z-N-PID control
algorithm0.005796 0.045 0.009981 Self-updating control algorithm 0.006355 0.042 0.009482 30 rpm Self-updating control algorithm 0.025631 0.049 0.011012 Self-updating control algorithm 0.020903 0.032 0.007471 50 rpm Self-updating control algorithm 0.041192 0.033 0.007800 Self-updating control algorithm 0.041027 0.026 0.004631 由表2可知,使用自更新控制算法的控制系统在转速为10 rpm的实验中出现最大极限误差,为0.042 N·m,该组数据标准差为0.009482。用Z-N-PID算法的控制系统在转速为30 rpm的实验中出现最大极限误差,为0.049 N·m,该组数据标准差为0.011012。经计算可知,在10、30、50 rpm的转速下,使用自更新控制算法的控制系统相比于使用Z-N-PID算法的控制系统的标准差分别降低了5.00%、32.16%、40.63%。
以上数据表明:在激光雷达轴系动摩擦力矩测量的稳态过程中,使用自更新控制算法的控制系统相比于使用Z-N-PID算法的控制系统所采集的数据拥有更好的收敛性,数据波动较小。
而后,对动摩擦力矩测量起始阶段产生的过冲量进行分析。分别取两种方法在抽取的10组样本中起始阶段最大值的均值作为平均过冲量,动摩擦力矩测量平均过冲量对比如图5所示。
从图5可以看出,系统的过冲量随转速的增加而增加。在10、30、50 rpm的转速下,使用自更新控制算法的控制系统相比于Z-N-PID算法对系统的过冲量均有较强的抑制力。经计算可知,在上述转速下,使用自更新控制算法的控制系统的平均过冲值为0.054667 N·m,使用Z-N-PID算法的控制系统的平均过冲值为0.062667 N·m,自更新控制算法相比Z-N-PID算法对过冲量的抑制力平均提高12.77%。
最后,对所抽取的10组样本摩擦力矩值进行分析,其结果如图6所示。
由图6可以知,使用自更新控制算法的控制系统,测试结果的波动随转速的增加而降低。经计算可知,使用自更新控制算法的控制系统在转速为10 rpm的实验中出现最大极限误差,为0.001948 N·m,在转速为50 rpm的实验中出现最小极限误差,为0.001052 N·m。使用Z-N-PID算法的控制系统在转速为30 rpm的实验中出现最大极限误差,为0.006855 N·m,在转速为50 rpm的实验中出现最小极限误差,为0.001865 N·m。此外,在10、30、50 rpm的转速下,使用自更新控制算法的控制系统相比于使用Z-N-PID算法的控制系统的重复测量误差分别降低了24.20%、71.66%、43.59%。
以上数据表明,在测量激光雷达轴系动摩擦力矩时,使用自更新控制算法的控制系统相比于使用Z-N-PID算法的控制系统可保持较高的重复测量精度,性能较为稳定,鲁棒性较好。
Self-updating control algorithm and verification of lidar shafting friction torque detection equipment
-
摘要: 对测量激光雷达轴系动摩擦力矩数据波动较大、重复测量精度低等问题开展研究,将基于测试主轴控制系统构建的数据云与GA-BP算法进行融合,提出了激光雷达轴系摩擦力矩检测设备自更新控制算法。以测试主轴的实际转速、理想转速、转速误差和转速误差变化率构建数据云,使用密度及距离信息实现数据的添加和删除,通过GA-BP算法实现在线控制参数的整定。以激光雷达轴系摩擦力矩测量设备测试主轴和被测轴系为研究对象,通过仿真实验证明该方法与使用Z-N-PID算法的控制系统相比,提高了系统抗干扰性。通过激光雷达轴系动摩擦力矩检测设备进行摩擦力矩检测,实验结果表明,所提出的自更新控制算法相比于Z-N-PID算法平均过冲量降低了12.77%,稳定后数据标准差降低了5.00%~40.63%,重复测量误差降低了24.20%~71.66%。Abstract: We studied the problems of large fluctuations in the dynamic friction torque data of lidar shafting and low precision of repeated measurements. The data cloud constructed based on the test spindle control system was fused with the GA-BP algorithm, and a self-updating control algorithm for lidar shafting friction torque detection equipment was proposed. The data cloud was constructed based on the actual speed, ideal speed, speed error and speed error change rate of the test spindle. The density and distance information were used to add and delete data, and the online control parameters were adjusted by the GA-BP algorithm. Taking the lidar shafting friction torque detection equipment to test the main shaft and the measured shaft system as the research object, the simulation experiment proves that this method improves the system’s anti-interference performance compared with the control system using the Z-N-PID algorithm. The friction torque is detected by lidar shafting friction torque detection equipment. The experimental results show that the proposed self-updating control algorithm reduces the average overshoot by 12.77% compared with the Z-N-PID algorithm, the data standard deviation after stability is reduced by 5.00%-40.63%, and the repeated measurement error is reduced by 24.20%-71.66%.
-
Key words:
- lidar shafting /
- friction torque /
- data cloud /
- self-updating control algorithm
-
表 1 系统响应曲线对比数据表
Table 1. Comparison data sheet of system response curve
Item Initial stage Interference phase Overshoot Adjust time/s Overshoot Adjust time/s Z-N-PID control algorithm 7.40% 5.968 17.70% 4.588 Self-updating control algorithm 0 7.197 10.40% 8.77 表 2 采集点数据波动比较
Table 2. Comparison of data fluctuations at collection points
Rotating speed Algorithm Average
value/N·mLimit
error/N·mStandard deviation 10 rpm Z-N-PID control
algorithm0.005796 0.045 0.009981 Self-updating control algorithm 0.006355 0.042 0.009482 30 rpm Self-updating control algorithm 0.025631 0.049 0.011012 Self-updating control algorithm 0.020903 0.032 0.007471 50 rpm Self-updating control algorithm 0.041192 0.033 0.007800 Self-updating control algorithm 0.041027 0.026 0.004631 -
[1] Rudloff L, Arghir M, Bonneau O, et al. Experimental analyses of a first generation foil bearing: startup torque and dynamic coefficients [J]. Journal of Engineering for Gas Turbines and Power-Transactions of the Asme, 2011, 133(9): 092501. [2] Hampson M R, Roberts E W, Watters R B, et al. Speed dependence of bearing torque-anomalous effects with grease lubrication[C]//15th European Space Mechanisms and Tribology Symposium, Sep 15–27, Noordwijk, The Netherlands, 2013. [3] Sim K, Lee Y, Song J W, et al. Effect of cooling flow on thermal performance of a gas foil bearing floating on a hot rotor [J]. Journal of Mechanical Science and Technology, 2018, 32(5): 1939-1954. doi: 10.1007/s12206-018-0401-8 [4] Fan J, Zhong J, Zhao J, et al. BP neural network tuned PID controller for position tracking of a pneumatic artificial muscle [J]. Technol Health Care, 2015, 23(S2): S231-S238. [5] 胡近朱, 高强, 侯远龙, 等. 某交流伺服系统的神经网络自抗扰控制[J]. 电光与控制, 2021, 28(2): 91-96. Hu Jinzhu, Gao Qiang, Hou Yuanlong, et al. Improved ADRC of AC servo system based on neural network [J]. Electronics Optics & Control, 2021, 28(2): 91-96. (in Chinese) [6] Li Xiaolu, Zeng Jingjing, Wang Hao, et al. Design and real-time imaging technology of three-dimensional scanning LiDAR [J]. Infrared and Laser Engineering, 2019, 48(5): 0503004. (in Chinese) doi: 10.3788/IRLA201948.0503004 [7] Guo Jingjing, Fei Xiaoyan, Ge Peng, et al. High-resolution three-dimensional imaging based on all-fiber photon-counting lidar system [J]. Infrared and Laser Engineering, 2021, 50(7): 20210162. (in Chinese) doi: 10.3788/IRLA20210162 [8] Gao Y, Ghasr M T, Zoughi R. Effects of and compensation for translational position error in microwave synthetic aperture radar imaging systems [J]. IEEE Transactions on Instrumentation and Measurement, 2020, 69(4): 1205-1212. doi: 10.1109/TIM.2019.2910340 [9] 高歆杨, 柯芳, 邹伟, 等. 基于模糊控制策略的快反镜伺服控制研究[J]. 兵工学报, 2020, 41(8): 1-10. Gao Xinyang, Ke Fang, Zou Wei, et al. Servo control of fast steering mirror based on fuzzy control strategy [J]. Acta Armamentarii, 2020, 41(8): 1-10. (in Chinese) [10] 吴辉, 司晨, 姜湖海, 等. 光电跟踪平台的观测自适应控制器[J]. 光学精密工程, 2021, 29(6): 1311-1320. doi: 10.37188/OPE.2020.0700 Wu Hui, Si Chen, Jiang Huhai, et al. Observational adaptive controller for photoelectric tracking platform [J]. Optics and Precision Engineering, 2021, 29(6): 1311-1320. (in Chinese) doi: 10.37188/OPE.2020.0700 [11] Yang J, Wu Z, Peng K, et al. Parameter selection of gaussian kernel SVM based on local density of training set [J]. Inverse Problems in Science and Engineering, 2021, 29(4): 536-548. doi: 10.1080/17415977.2020.1797716 [12] Zhou H, Liu H, Zhang Y, et al. An outlier detection algorithm based on an integrated outlier factor [J]. Intelligent Data Analysis, 2019, 23(5): 975-990. doi: 10.3233/IDA-184227 [13] 王涛, 翟绪恒, 孟丽岩. 在线自适应神经网络算法及参数鲁棒性分析[J]. 振动与冲击, 2019, 38(8): 210-217. Wang Tao, Zhai Xuheng, Meng Liyan. An online adaptive neural network algorithm and its parameters robustness analysis [J]. Journal of Vibration and Shock, 2019, 38(8): 210-217. (in Chinese) [14] Xu Huixiang, Kong Guoli. Development of pressure control system for laser infrared multipass cell using Ziegler-Nichols-PID algorithm [J]. Infrared and Laser Engineering, 2020, 49(9): 20190551. (in Chinese) doi: 10.3788/IRLA20190551 [15] Cyriac F, Lugt P M, Bosman R. Yield stress and low-tem-perature start-up torque of lubricating greases [J]. Tribology Letters, 2016, 63(1): 1-10. doi: 10.1007/s11249-016-0692-9 [16] Du Zhaofei. Researh on friction torque measurement of rolling bearing under radial load[D]. Wuhan: Huazhong University of Science and Technology, 2017. (in Chinese)