技术领域
[0001] 本
发明属于
降落伞测量技术和计算机
机器学习领域，具体涉及一种基于混合神经网络模型(Hybrid Neural Network Model,HNNM)的大变形柔性体动态
应力补偿方法。
背景技术
[0002] 目前，在降落伞测量技术、
气动力学和结构力学
基础上建立起来的CFD(Computational Fluid Dynamics)-MSD(Mass Spring Damper)耦合模型，利用计算
流体力学和结构动力学的基础知识，进行数值求解，模拟大变形柔性体工作过程的结构和流场变化情况。这种模型可以较好地反映大变形柔性体所受的气动弹性力本质。但大变形柔性体透气量、弹性系数及阻尼系数的不确定性是其
缺陷所在，特别是耦合模型是建立在半理论半经验的基础上，大变形柔性体应力、应变初始参数的设定主要是基于某些假设，模型最终的计算结果缺乏真实实验数据来进行量化验证。
[0003] 针对大变形柔性体工作过程中存在的非线性问题，必须对其应力、应变测量
精度进行非线性补偿或校正，一般采用
硬件补偿联合
软件补偿的方式：一方面从硬件上对应力
传感器输出
信号予以补偿，以保证数据可靠、准确；另一方面要采用软件手段对影响应力传感器输出的其他可变因素进行相应的补偿。常用的软件补偿方法有最小二乘法、分段线性化法和基于神经网络的自适应方法。在各种神经网络方法中，反向传播(Back Propagation，BP)神经网络是近年来受到
模式识别、信号理论与控制、故障诊断等领域研究人员广泛关注的一种
前馈神经网络，具有逼近任意非线性函数的能力，在非线性系统建模中被广泛应用。但BP神经网络存在收敛速度慢，容易收敛于局部极小值，鲁棒性差，学习率、动量项系数和初始权值等参数难以调整等缺陷。
[0004] 混合神经网络模型主要用于建立机理模型的数学描述，是先验知识建模与神经网络建模的综合应用，使用混合神经网络模型描述大变形柔性体受力模型，能够较好地反映大变形柔性体的非线性动力学特性，逼近非线性应力、应变之间的非线性映射关系。但是，
现有技术中混合神经网络模型主要在小变形、大位移及大转动应用场景中对非线性应力、应变进行补偿，尚无应用于大变形柔性体动态应力补偿的相关描述。
发明内容
[0005] 本发明的目的是在于提供一种基于混合神经网络模型的大变形柔性体动态应力补偿方法，可以有效应用于大变形柔性体的动态应力补偿，尤其是降落伞、
气垫船等设备工作过程中动态应力补偿。
[0006] 实现本发明目的的技术解决方案为：一种基于混合神经网络模型的大变形柔性体动态应力补偿方法，包括如下步骤：
[0007] 步骤1、采集大变形柔性体静态试验数据，采集的内容包括：采集应力传感器
输出电压、环境
温度、大变形柔性体纵向应变、大变形柔性体纵向应力；
[0008] 步骤2、构建混合神经网络模型；
[0009] 步骤3、采用步骤1中的静态试验数据，对混合神经网络模型进行训练；
[0010] 步骤4、使用步骤3完成训练的混合神经网络模型，对动态试验或实际应用场景中的大变形柔性体的动态应力进行补偿。
[0011] 本发明与现有技术相比，其显著优点为：1)本发明的动态应力补偿方法，是一种基于混合神经网络模型的大变形柔性体动态应力补偿方法，可有效应用于大变形柔性体的动态应力补偿，尤其是降落伞、气垫船等设备工作过程中的动态应力补偿，同时为大变形柔性体流固耦合计算提供真实的实验数据验证支持。2)本发明利用知识基模型的引导特性，将先验知识引入大变形柔性体的动力学建模中，可有效减小混合神经网络规模，缩短建模周期；3)本发明利用BP神经网络良好的非线性函数逼近能力，可以有效补偿复杂工作环境下传感器输入输出非线性特性及外部因素导致的非线性误差，显著提高应力传感器的测量精度；4)本发明使用LM(Levenberg-Marquardt)
算法来优化BP神经网络，可以较好地解决BP神经网络收敛速度慢，且可能收敛至局部最小值等问题。
附图说明
[0013] 图2混合神经网络模型。
[0014] 图3
实施例中反向传播神经网络模型交叉验证示意图。
[0015] 图4实施例中4种不同
风速条件下混合神经网络模型动态应力补偿曲线，其中，图(a)为风速为20m/s，图(b)为风速为30m/s，图(c)为风速为40m/s，图(d)为风速为50m/s。
具体实施方式
[0016] 结合图1，本发明的一种基于混合神经网络模型的大变形柔性体动态应力补偿方法，包括如下步骤：
[0017] 步骤1、采集大变形柔性体静态试验数据，采集的内容包括：采集应力传感器输出电压、
环境温度、大变形柔性体纵向应变、大变形柔性体纵向应力；所述大变形柔性体静态试验数据具体为X＝(x1,x2,ε,σ)T，其中，x1为应力传感器输出电压，x2为环境温度，ε为大变形柔性体纵向应变，其计算公式为ε＝ΔL/L，该式中，ΔL大变形柔性体纵向伸长量，L为大变形柔性体长度，σ为大变形柔性体纵向应力，其计算公式为σ＝F/(w*t)，该式中，F为大变形柔性体纵向
载荷，w为大变形柔性体宽度、t为大变形柔性体厚度；样本数据均作归一化处理。
[0018] 步骤2、构建混合神经网络模型f(x1,x2)，所述混合神经网络模型包括知识基模型f1(x1)f2(x2)和反向传播神经网络模型f3(x1,x2)，混合神经网络模型的期望输出为大变形柔性体纵向应力，函数关系式表示为：
[0019] f(x1,x2)＝f1(x1)f2(x2)+f3(x1,x2)
[0020] 其中符号的具体含义及计算公式如下：
[0021] f1(x1)描述大变形柔性体力学模型的力学特征，函数关系式表示为：其中，w1、w2、w3、w4均为待确定权值参数；
[0022] f2(x2)描述大变形柔性体力学模型的环境温度特征，函数关系式表示为：其中，E为大变形柔性体
弹性模量，
其计算公式为E＝σ/ε，Einit、Efinal分别为静态试验温度取值范围内大变形柔性体弹性模量的初始值、最终值，dx为时间步长，w5、w6为待确定权值参数；W1＝(w1,w2,...w6)T为知识基模型的权值向量；
[0023] f3(x1,x2)为描述大变形柔性体力学模型的纵向应力与知识基模型的误差，形式为反向传播神经网络，网络结构为2×3×1，即一个
输入层、一个隐层和一个
输出层，输入层神经元个数为2，输出层神经元个数为1，均采用线性恒等函数作为神经元激励函数；隐层神经元个数为n，以一组逐次增高的幂函数xi-1作为神经元激励函数，i表示隐层的第i个神经元；W2＝(w7,w8,...,wn+6)T为反向传播神经网络模型隐层神经元与输出层神经元之间的权值向量；反向传播神经网络的函数关系式表示为：
式中，X1＝(x1,x2)， 为激励函数。
[0024] 步骤3、采用步骤1中的静态试验数据，对混合神经网络模型进行训练；具体包括以下步骤：
[0025] 步骤3-1、对知识基模型进行训练，具体为：
[0026] 将样本数据X中的应力传感器输出电压、环境温度作为知识基模型输入；将样本数据X中的大变形柔性体纵向应力作为知识基模型期望输出；采用非线性回归方法进行训练，使用拟
牛顿法中的BFGS算法进行
迭代，求解待确定权值向量W1；
[0027] 步骤3-2、对反向传播神经网络模型进行训练，具体为：
[0028] 将样本数据X中的应力传感器输出电压、环境温度作为反向传播神经网络模型输入；将对应样本数据X中大变形柔性体纵向应力与完成训练的知识基模型输出之差作为反向传播神经网络模型期望输出；采用反向传播算法进行训练，使用LM算法进行迭代，求解待确定权值向量W2。
[0029] 步骤4、使用步骤3完成训练的混合神经网络模型，对动态试验或实际应用场景中的大变形柔性体的动态应力进行补偿。具体为：
[0030] 步骤4-1、采集动态试验或实际应用场景中的大变形柔性体动态样本数据，作归一化处理；
[0031] 步骤4-2、将动态样本数据中的应力传感器输出电压、环境温度作为完成训练的知识基模型输入；将动态样本数据中的应力传感器输出电压、环境温度作为完成训练的反向传播神经网络模型输入；
[0032] 步骤4-3、将知识基模型输出与反向传播神经网络模型输出之和作为混合神经网络模型的输出，作反归一化处理，即为大变形柔性体的动态补偿应力。
[0033] 下面结合实施例对本发明做进一步详细的描述：
[0034] 本发明的思路是构建对动态应力进行非线性补偿的混合神经网络模型，采用该模型对动态实验或实际应用场景中的动态应力进行非线性补偿，操作流程如图1所示，具体分为以下几个步骤：
[0035] 步骤1、采集大变形柔性体静态试验数据。
[0036] 在静态实验环境中，除了在大变形柔性体的纵向方向表面安装应力传感器(载荷50kg)之外，还安装了TI公司的LM95172-Q1数字温度传感器。在不同环境温度(-40℃～125℃)的实验环境下，使用高精度万能试验机输出大小不等的拉力F，沿大变形柔性体纵向方向作用于大变形柔性体两端，采集应力传感器输出电压、环境温度、大变形柔性体纵向伸长量。本发明所用的应力传感器已经在
专利《一种测量柔性织物应力作用的应变传感器》中公开，专利号为201320332161.1。
[0037] 计算大变形柔性体纵向应变，其计算公式为ε＝ΔL/L，式中L取9.5mm；大变形柔性体纵向应力，其计算公式为σ＝F/(w*t)，式中w取6.5mm，t取0.9mm；静态试验温度取值范围内大变形柔性体弹性模量的初始值Einit＝0.98987，最终值Efinal＝-3.65353，时间步长dx＝0.55289。
[0038] 静态试验数据样本均作归一化处理。
[0039] 步骤2、构建混合神经网络模型，如图2所示。
[0040] 步骤3、采用步骤1中的静态试验数据样本(样本数m＝300)，对步骤2中的混合神经网络模型进行训练，具体包括以下步骤：
[0041] 步骤3-1、将样本数据集分成两个集合：训练集(样本数为m1，m1＝％70*m)和交叉验证集(样本数为m2，m2＝％30*m)；
[0042] 步骤3-2、对知识基模型进行训练，具体为：
[0043] 步骤3-2-1、设定训练最大迭代次数为500，学习率为0.02，收敛精度为0.5；设定知识基模型的其他参数：
[0044] 步骤3-2-2、基于训练集样本，使用非线性回归方法进行训练。采用拟牛顿法中的BFGS(Broyden-Fletcher-Goldfarb-Shanno)算法进行迭代：调节权值向量W1，使代价函数J1(W1)满足收敛精度，即 从而求解知识基模型的权值向量W1。代价函数为式中，λ1表示权重衰减参数。
[0045] 步骤3-2-4、将训练集样本带入未训练的知识基模型进行训练，经过109轮迭代后收敛，统计结果如下表所示。
[0046]
[0047] 此时J1(W1)＜0.5，即基于训练集的代价函数满足收敛精度。计算知识基模型权值向量的估计值为W1＝(0.04805,1.49642,-1.76484,1.09689,2.06986,1.30183)。
[0048] 步骤3-2-5、将交叉验证集样本带入完成训练的知识基模型进行验证，统计结果如下表所示。
[0049]
[0050] 此时J1(W1)＜0.5，即基于交叉验证样本集的代价函数满足收敛精度，拟合结果如图3所示。因此，f1(x1)f2(x2)即为满足收敛精度的知识基模型。
[0051] 步骤3-3、对反向传播神经网络模型进行训练，具体为：
[0052] 步骤3-3-1、设定训练最大迭代次数为1000，学习率为0.02，收敛精度为0.1；
[0053] 步骤3-3-2、基于训练集样本，使用LM算法对反向传播神经网络进行训练。通过迭代方式调节权值向量W2，使代价函数J2(W2)满足收敛精度，从而求解反向传播神经网络的权值向量W2。代价函数为 式中，λ2表示权重衰减参数；Y(i)表示大变形柔性体纵向应力σ(i)与知识基函数 之差，即
[0054] 步骤3-3-3、将训练集样本带入未训练的反向传播神经网络模型进行训练，隐层神经元
节点数为7时，经过69轮迭代后收敛，此时训练结果最好，J2(W2)＝0.06933，W2＝(0.65679，2.9914，1.01206，0.94953,1.05043，-1.64913，1.44707)。
[0055] 步骤3-3-4、将交叉验证集样本带入完成训练的反向传播神经网络模型进行验证，计算出代价函数J2(W2)＝0.08642。此时J2(W2)＜0.1，即代价函数满足收敛精度。因此，f3(x1,x2)即为满足收敛精度的反向传播神经网络模型。
[0056] 步骤4、使用步骤3完成训练的混合神经网络模型，对动态试验或实际应用场景中的大变形柔性体的动态应力进行补偿。
[0057] 为了比较不同类型神经网络算法的训练速度及收敛精度，将混合神经网络与标准BP神经网络(梯度下降)、改进BP神经网络(LM)及RBF(Radial Basis Function)算法进行对比。实验中，样本数为450，学习率为0.02，最大迭代次数为1000轮，目标误差精度为0.1；参数方面，RBF的径向基函数分布为1.1。实验环境为：Intel Core i5-4570(3.6GHz)，8.0GB RAM，Windows 8.1(32bit)，MATLAB R2014a。结果如下表所示，在这4种算法中，虽然RBF神经网络的收敛速度最快，但混合神经网络收敛速度仅次于RBF神经网络，且收敛精度最高，对动态应力的补偿效果最好。
[0058]网络类型 隐层神经元数量 训练步数 收敛误差
标准BP 12 1515 0.17431
改进BP 7 436 0.09794
RBF 5 139 0.08263
混合神经网络 7 231 0.06438
[0059] 如图4(a)所示，环境温度为25℃,风速为20m/s，将混合神经网络补偿前后结果作比较，整个时间序列种变化趋势基本保持一致；在峰值过后出现有规律的“呼吸”现象，这符合大变形柔性体工作过程的实际情况。此外，补偿算法考虑环境温度等因素对柔性体弹性模量的影响，因此在大变形柔性体应力变化过程中，补偿前后结果出现一定的差异性。相比之下，补偿结果更加具备合理性。
[0060] 图4(b)、图4(c)、图4(d)分别表示在环境温度为25℃，风速为30m/s，40m/s，50m/s条件下的大变形柔性体工作过程的应力变化曲线，采用本发明都取得了较好的补偿结果，这得益于本发明方法利用采集应力传感器输出电压、环境温度与大变形柔性体纵向应力之间的非线性关系来进行补偿运算，有较强的泛化能力。
[0061] 由上可知，采用本发明的方法能够准确有效地补偿大变形柔性体工作过程的动态应力，能够满足应用要求。
