基于混合神经网络模型的大变形柔性体动态应力补偿方法
【技术领域】
[0001] 本发明属于降落伞测量技术和计算机机器学习领域，具体涉及一种基于混合神经 网络模型（HybridNeuralNetworkModel,HNNM)的大变形柔性体动态应力补偿方法。
【背景技术】
[0002] 目前，在降落伞测量技术、气动力学和结构力学基础上建立起来的 CFD(ComputationalFluidDynamics)-MSD(MassSpringDamper)親合模型，利用计算流体 力学和结构动力学的基础知识，进行数值求解，模拟大变形柔性体工作过程的结构和流场 变化情况。这种模型可以较好地反映大变形柔性体所受的气动弹性力本质。但大变形柔性 体透气量、弹性系数及阻尼系数的不确定性是其缺陷所在，特别是耦合模型是建立在半理 论半经验的基础上，大变形柔性体应力、应变初始参数的设定主要是基于某些假设，模型最 终的计算结果缺乏真实实验数据来进行量化验证。
[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神经网络收敛速度慢，且可能收敛至局部最小值等问题。
【附图说明】
[0012] 图1本发明的具体操作流程图。
[0013] 图2混合神经网络模型。
[0014] 图3实施例中反向传播神经网络模型交叉验证示意图。
[0015] 图4实施例中4种不同风速条件下混合神经网络模型动态应力补偿曲线，其中，图 (a)为风速为20m/s，图（b)为风速为30m/s，图（c)为风速为40m/s，图（d)为风速为50m/ So
【具体实施方式】
[0016]结合图1，本发明的一种基于混合神经网络模型的大变形柔性体动态应力补偿方 法，包括如下步骤：
[0017] 步骤1、采集大变形柔性体静态试验数据，采集的内容包括：采集应力传感器输出 电压、环境温度、大变形柔性体纵向应变、大变形柔性体纵向应力；所述大变形柔性体静态 试验数据具体为X= (Xl，x2,e，〇)T，其中，Xl为应力传感器输出电压，x2为环境温度，e 为大变形柔性体纵向应变，其计算公式为e=AL/L，该式中，AL大变形柔性体纵向伸长 量，L为大变形柔性体长度，〇为大变形柔性体纵向应力，其计算公式为〇 =FAw*t)，该 式中，F为大变形柔性体纵向载荷，w为大变形柔性体宽度、t为大变形柔性体厚度；样本数 据均作归一化处理。
[0018]步骤2、构建混合神经网络模型f(Xl，x2)，所述混合神经网络模型包括知识基模型 fiOOtOg和反向传播神经网络模型f3(Xl，x2)，混合神经网络模型的期望输出为大变形 柔性体纵向应力，函数关系式表示为：
[0019]f(x1；x2) =fj(xj)f2 (x2) +f3 (x1；x2)
[0020] 其中符号的具体含义及计算公式如下：
[0021] f\(Xl)描述大变形柔性体力学模型的力学特征，函数关系式表示为：
，其中，Wpw2、w3、w4均为待确定权值参数；
[0022] f2(x2)描述大变形柔性体力学模型的环境温度特征，函数关系式表示为：
其中，E为大变形柔性体弹性模 量，其计算公式为E= 〇 /e，Einit、Efinal分别为静态试验温度取值范围内大变形柔性体弹 性模量的初始值、最终值，dx为时间步长，w5、w6为待确定权值参数;W(wuw2, . . .w6)T为 知识基模型的权值向量；
[0023]f3(Xl，x2)为描述大变形柔性体力学模型的纵向应力与知识基模型的误差， 形式为反向传播神经网络，网络结构为2X3X1，即一个输入层、一个隐层和一个输 出层，输入层神经元个数为2,输出层神经元个数为1，均采用线性恒等函数作为神 经元激励函数；隐层神经元个数为n，以一组逐次增高的幂函数矿1作为神经元激励 函数，i表示隐层的第i个神经元;W2=(w7,w8,. . .，wn+6)T为反向传播神经网络模 型隐层神经元与输出层神经元之间的权值向量；反向传播神经网络的函数关系式表 示为：
，式中，X:=(xx2)，
为激励函数。
[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°C~ 125°C)的实验环境下，使用高精度万能试验机输出大小不等的拉力F，沿大变形柔性体纵 向方向作用于大变形柔性体两端，采集应力传感器输出电压、环境温度、大变形柔性体纵向 伸长量。本发明所用的应力传感器已经在专利《一种测量柔性织物应力作用的应变传感器》 中公开，专利号为201320332161. 1。
[0037] 计算大变形柔性体纵向应变，其计算公式为e=AL/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中的混合 神经网络模型进行训练，具体包括以下步骤：
[0