摘要
在全局摩阻测量中,薄油膜技术可以很好地表征表面摩阻的分布情况。用特定波长的紫外线照射添加了荧光显色分子的不同厚度的油膜,油膜将发出不同的亮度。利用该原理通过检测受激发的荧光油膜灰度值可解算出相应油膜的厚度。本次采用BP神经网络及极限学习机(Extreme learning machine, ELM)神经网络搭建模型完成了荧光油膜厚度与灰度关系的预测,运用Hopfield神经网络完成了相应参数的辨识。实验表明,ELM神经网络模型、BP神经网络模型及插值法模型的预测误差分别为5.150%、5.485%和5.935%。通过Hopfield神经网络辨识,光源功率、光距和曝光系数等影响因素的参数误差率控制在1%左右,达到实际工程运用的要求。与传统插值法相比,通过神经网络可获得更高的精度,为荧光油膜灰度与厚度研究提供了一种可行的方法。
利用荧光油膜在激光照射下的显色反应,由荧光油膜灰度表征出油膜厚度这一关键因素,进而通过得到油膜厚度解算该处的摩阻分布。早在1977年,由Tanne
在标定过程中,一般选择插值法或数据拟合的方式来处理数据使其平滑。传统的插值法或者数据拟合的方式在一般厚度(如毫米级)的条件下比较实用,但是在面对高精度(如微米级)时,则需要进行算法改进来达到精度要求,其改进的过程较为繁琐,也存在着一定的困难,而两者都有各自的缺点,比如数据拟合方式,由于该方法不一定通过所有数据点,在精度方面也存在一定的局限性;传统插值法在靠近数据两端点插值时,会存在“龙格现象
本次采用神经网络对荧光油膜建立模型,针对神经网络对单输入、单输出模型有较好预测的特点,与传统插值法结果进行比对,验证了通过神经网络建立预测模型,有精度更高、预测更准确的优势。
本次采用一元全区插值法进行荧光油膜灰度与厚度关系的模型建立,其数学模型表示为
(1) |
式中:从属于个结点,为不同于的其他结点。通过7次拉格朗日插值计算求得插值点处的函数近似值,且插值点在该个结点中(=8)选取满足
(2) |
当插值点位于区间外时,则属特殊情况,将取该端的个结点进行插值计算。
BP神经网络为多层前馈网络,遵循误差反向传播训练机制,围绕梯度下降法中心思想展开,采用梯度搜索方
BP网络结构如图

图1 BP神经网络逼近的结构
Fig.1 Structure of BP neural network approximation

图2 用于逼近的BP网络
Fig.2 BP network for approximation
前向传播是网络的输出运算过程,隐层神经元的输入过程实际为所有输入的加权函数求和过程,令隐层神经个数为,则有
(3) |
式中:为网络层间连接权值;为限层结点;为输入层结点。
取,则网络的输出与预想输出的误差为
(4) |
误差性能指标函数为
(5) |
逆向传播通过学习算法,对各层间的权值进行调整。输出层与隐层的连接权值计算为
(6) |
式中为激发函数。时刻网络的权值为
(7) |
隐层及输入层连接权值计算为
(8) |
(9) |
时刻网络权值为
(10) |
当上一次的权值对本次权值产生了影响,就需加入约束因子,此时的权值有
(11) |
(12) |
式中:为学习速率,为约束因子,,。雅可比矩阵(输入输出灵敏度信息)算法表达式为
(13) |
极限学习机(Extreme learning machine, ELM)神经网络较BP神经网络具有更高的精度和运算速
(14) |
式中:为激活函数,、分别为输入、输出权重,为偏置。该网络目标函数为
(15) |
则可得矩阵形式
(16) |
式中:为隐层输出,为输出权重,为期望值。训练过程中得到
(17) |
可得损失目标函数
(18) |
假定用于参数辨识的约束函数与Hopfield网络约束函数相等,对其状态方程矩阵进行分解变换可得网络的权值矩阵和神经元外输入矩阵,将其代入网络状态方程式进行运算,收敛状态可得到较稳定的辨识参数结果。令一个系统的状态方程为
(19) |
式中:、为待辨识的参数矩阵,取,,令为状态矢量,为控制输入,则其参数辨识过程就是对向量的求解过程。设可调系统为
(20) |
式中:,,取,则可得到状态偏差方程
(21) |
定义基于状态偏差变化率的系统参数辨识函数为
(22) |
(23) |
(24) |
Hopfield网络约束函数趋于收敛的过程,就是估计矩阵和接近实际矩阵和的过
(25) |
式中:,;为调整系数;为网络收敛因子。
令,并取,则状态方程可以转换为
(26) |
Hopfield网络的约束函数为
(27) |
由于,取Hopfield网络的输出等于待辨识参数,则
(28) |
(29) |
取参数辨识的约束函数与Hopfield网络约束函数相
(30) |
通过,可以换算为
(31) |
(32) |
通过
(33) |
(34) |
将式(
(35) |
荧光油膜的标定过程是一个数据实时采集的过程,本次通过搭建载玻片斜面的方法进行油膜灰度与厚度相关信息的采

图3 荧光油膜灰度图
Fig.3 Grayscale of fluorescent oil film

图4 油膜网格图
Fig.4 Oil film grid
由盖破片和载玻片构成的斜面所采集的数据并不是直接为荧光油膜灰度与厚度的关系数据,而是灰度与横坐标长度的关系,还需进行三角变换处理,如
(36) |

图5 斜面三角变换
Fig.5 Transformation of slope trigonometric
序号 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|
油膜灰度值 | 32 | 34 | 36 | 40 | 41 | 44 | 47 | 52 | 54 | 57 |
油膜厚度值 | 0.000 00 | 0.005 53 | 0.006 64 | 0.007 14 | 0.010 51 | 0.014 29 | 0.021 43 | 0.026 55 | 0.028 57 | 0.033 19 |
序号 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 |
油膜灰度值 | 59 | 62 | 65 | 67 | 71 | 72 | 74 | 77 | 80 | 82 |
油膜厚度值 | 0.038 16 | 0.042 68 | 0.049 23 | 0.050 01 | 0.062 50 | 0.063 05 | 0.071 44 | 0.078 58 | 0.096 79 | 0.100 01 |
序号 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 | 30 |
油膜灰度值 | 84 | 86 | 88 | 90 | 92 | 93 | 95 | 97 | 99 | 100 |
油膜厚度值 | 0.114 3 | 0.128 59 | 0.148 78 | 0.157 16 | 0.171 46 | 0.185 73 | 0.200 02 | 0.274 89 | 0.285 74 | 0.300 03 |
序号 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|
油膜灰度值 | 33 | 35 | 43 | 58 | 68 | 79 | 85 | 87 | 96 | 101 |
油膜厚度值 | 0.003 87 | 0.006 08 | 0.012 72 | 0.035 95 | 0.057 15 | 0.085 72 | 0.121 44 | 0.142 87 | 0.228 60 | 0.347 35 |
将标定实验得到荧光油膜灰度与厚度初始数据(

图6 训练过程图
Fig.6 Training process diagram

图7 神经网络训练
Fig.7 Neural network training

图8 误差曲线
Fig.8 Error curves
将训练好的模型通过

图9 BP神经网络预测结果
Fig.9 Prediction results of BP neural network
ELM网络较BP网络具有运算速度快的优势,且精度较高,训练模型的预测结果如

图10 ELM神经网络预测结果
Fig.10 Prediction results of ELM neural network
绘制BP神经网络、ELM网络及传统插值法预测曲线,如

图11 数据预测对比
Fig.11 Comparison of data prediction
测试灰度 | 真实 厚度/mm | BP预测 厚度/mm | ELM预测 厚度/mm | 插值预测 厚度/mm | BP误差/% | ELM误差/% | 插值误差/% |
---|---|---|---|---|---|---|---|
33 | 0.003 87 | 0.003 91 | 0.003 62 | 0.003 39 | 1.033 | 6.460 | 12.40 |
35 | 0.006 08 | 0.004 83 | 0.006 63 | 0.006 75 | 20.56 | 9.046 | 11.02 |
43 | 0.012 72 | 0.011 90 | 0.012 12 | 0.013 81 | 6.447 | 4.171 | 8.569 |
58 | 0.035 95 | 0.034 52 | 0.035 54 | 0.035 67 | 3.978 | 1.140 | 0.779 |
68 | 0.057 15 | 0.054 71 | 0.053 29 | 0.053 59 | 4.270 | 6.750 | 6.229 |
79 | 0.085 72 | 0.088 16 | 0.090 43 | 0.092 02 | 2.846 | 5.495 | 7.350 |
85 | 0.121 44 | 0.117 21 | 0.118 99 | 0.120 90 | 3.483 | 2.070 | 0.445 |
87 | 0.142 87 | 0.140 96 | 0.131 86 | 0.139 26 | 1.337 | 7.706 | 2.527 |
96 | 0.228 60 | 0.240 36 | 0.233 49 | 0.233 67 | 2.783 | 2.139 | 2.218 |
101 | 0.347 35 | 0.319 17 | 0.324 68 | 0.320 23 | 8.113 | 6.527 | 7.808 |
平均误差/% | — | — | — | — | 5.485 | 5.150 | 5.935 |
由
Hopfield神经网络进行模型的参数辨识,需要根据经验公式或者手册给出初定的参数
(37) |
式中:为荧光亮度(灰度值);为初始光强;为荧光油膜的厚度;为荧光效率等参数。因为影响油膜灰度值的因素不止油膜本身的厚度因素,还有周围环境因素以及油膜的自身特性,如浓度,所以对该线性模型做变换,可以得到
(38) |
因影响因素可能为非线性系统和常数,所以、系数定义为影响因素的权重系数,根据经验式,可以得到假定影响因素参数,如
影响因素 | 符号 | 类型 | 假定权重初始值 |
---|---|---|---|
功率 | W | 非线性 | 0.13 |
光距 | R | 非线性 | 0.21 |
曝光系数 | T | 常数 | 33.3 |
通过Hopfield神经网络参数辨识,可得以下辨识结果,如图

图12 紫外光源功率参数辨识
Fig.12 Power parameter identification of UV light source

图13 紫外光源距离参数辨识
Fig.13 Distance parameter identification of UV light source

图14 曝光系数参数辨识
Fig.14 Identification of exposure coefficient parameters
通过Hopfield神经网络所辨识的结果参数如
影响因素 | 初始权重 | 辨识后的权重 | 误差百分比/% |
---|---|---|---|
功率 | 0.13 | 0.128 9 | 0.85 |
光距 | 0.21 | 0.207 8 | 1.05 |
曝光系数 | 33.3 | 33.129 6 | 0.51 |
通过Hopfield神经网络进行参数辨识,能够很好地找出功率、光距和曝光系数等参数对荧光油膜灰度成像的影响程度,且辨识结果误差范围在1%左右,达到实际工程应用的标准。
经过实验分析论证,ELM神经网络模型、BP神经网络模型及插值法模型的预测误差分别为5.150%、5.485%和5.935%。通过神经网络搭建模型较传统插值法精度有所提高;而ELM神经网络具备比BP神经网络运算速度快、精度更高的优势。利用Hopfield神经网络进行参数辨识,其辨识得到参数的误差均控制在1%左右,后续可以继续调整Hopfield神经网络的参数,使误差率达到期望的效果。与传统插值法进行对比,利用神经网络搭建模型可以获得更高的精度,为后续荧光油膜灰度与厚度研究提供了一种可行的方法。
References
Tanner L. A skin friction meter, using the viscosity balance principle,suitable for use with flatrocurved metal surfaces[J]. Phys E:Scientific Instr, 1977, 10: 278-284. [百度学术]
Monson D J, Matree G G. Boundary-layer transition and global skin friction measurement with an oil-fringe imaging technique[J]. Society of Automotive Engineers, 1993, 102: 1829-1843. [百度学术]
Brown J L, Naughton J W. The thin oil film equation[M]. Washington DC: NASA_AMES Research Center, 1999. [百度学术]
Liu Tianshu, Montefort J. Global luminescent oil-friction meter[J]. AIAA Journal, 2008,46(2):476-485. [百度学术]
李鹏,明晓.风力机叶片的全局表面摩擦力测量的荧光油膜法[J].南京航空航天大学学报,2011,43(5):581-585. [百度学术]
Li Peng, Ming Xiao. Fluorescent oil film method for global surface friction measurement of wind turbine blades[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2011,43(5): 581-585. [百度学术]
李鹏.全局表面摩擦应力直接测量技术研究[D].南京:南京航空航天大学,2012. [百度学术]
Li Peng. Research on direct measurement technology of global surface friction stress[D]. Nanjing: Nanjing University of Aeronautics and Astronautics,2012. [百度学术]
张雨浓,李名鸣,陈锦浩,等. 龙格现象难题破解之系数与阶次双确定方法[J]. 计算机工程与应用,2013,49(3): 44-49. [百度学术]
Zhang Yunong, Li Mingming, Chen Jinhao, et al. Double determination of coefficient and order for solving the Runge problem[J]. Computer Science and Technology, 2013,49(3): 44-49. [百度学术]
刘冰心, 王宁, 张冬. 基于BP神经网络的智能电网配电系统改进算法的研究[J]. 现代电子技术, 2012, 35(21): 143-144,148. [百度学术]
Liu Bingxin, WANG Ning, ZHANG Dong. Research on improved algorithm of smart grid distribution system based on BP neural network[J]. Modern Electronic Technology, 2012, 35(21): 143-144,148. [百度学术]
霍冠良,宁志华.基于BP神经网络的含褶皱复合材料强度预测[J].南京航空航天大学学报,2020,52(3): 460-467. [百度学术]
Huo Guanliang, Ning Zhihua. BP neural network for strength prediction of folded composites[J]. Journal of Nanjing University of Aeronautics & Astronautics,2020,52(3): 460-467. [百度学术]
Zhai J, Shao Q, Wang X. Architecture selection of ELM networks based on sensitivity of hidden nodes[J]. Neural Processing Letters, 2016, 44(2): 471-489. [百度学术]
卓新建, 章祥荪, 徐宗本. 关于Hopfield-型神经网络收敛性和求能量函数极小解的可靠性的研究[J]. 系统科学与数学, 2001(2): 152-162. [百度学术]
Zhuo Xinjian. ZHANG Xiangsun, XU Zongben. Research on convergence of Hopfield-type neural networks and reliability of minimum solutions of energy functions[J]. Systems Science and Mathematics, 2001(2): 152-162. [百度学术]
章毅, 周明天, 王平安. 关于神经网络的能量函数 [J]. 计算机研究与发展, 1999, 36(7): 794-799. [百度学术]
Zhang Yi, ZHOU Mingtian, WANG Pingan. On the energy function of neural networks[J]. Computer Research and Development, 1999, 36(7): 794-799. [百度学术]