一种基于载荷与状态结合的飞行器结构强度数字孪生方法与流程

1.本发明涉及飞行器结构健康监测与载荷识别技术领域，特别涉及在飞行器易损伤位置下的结构健康监测与载荷识别，同时根据识别出的系统状态参数与载荷参数，构建飞行器数字孪生样机模型，实现对飞行器结构的动态监测与虚拟演化。


背景技术：

2.现代大型复杂结构在全寿命周期内受到疲劳、意外撞击等不安全因素影响，结构健康监测是保障其安全可靠运行的重要手段。然而，现有结构健康监测方法存在着针对复合材料多失效模式下安全状态表征欠完备、不确定性条件下损伤参数反演精度低和传感器测点不确定性优化可靠性低等问题。
3.现代飞行器在复杂多变的服役环境中不可避免地会受到一些疲劳载荷、偶然的冲击载荷等多种不安全因素的影响，使得飞行器结构受到损伤甚至破坏，从而导致灾难性的事故。欧美国家从1970年以来在结构健康监测与安全状态评估这个领域已经有了较大的发展，并且其理论研究成果已经在服役的一些先进飞行器上得到了广泛的应用。同时，数字孪生技术也广泛的应用于智能制造、城市管理等诸多问题中。结合数字孪生技术，可对孪生对象实现动态监测，历史溯源追查，任务分配与动态管理，未来行为预测等多种应用，具有良好的发展前景。


技术实现要素：

4.本发明提供一种针对飞行器结构载荷与状态的动态孪生方法，建立飞行器结构的数字孪生样机模型，实现对飞行器结构所受载荷与结构健康状态的动态监测。本发明充分考虑实际工程问题中普遍存在的复合材料失效与结构损伤模式，构建能够合理表征典型机翼结构在静力与疲劳载荷作用下的结构损伤模型，结合结构载荷反演技术，建立了基于结构载荷与状态的数字孪生模型。
5.本发明采用的技术方案为：一种基于结构载荷与状态的数字孪生模型建模方法，该方法实现步骤如下：
6.第一步：根据实际物理模型，由结构约束条件，确定结构未知外载荷作用点，根据模型几何参数确立结构载荷识别的网格模型与结构状态识别的网格模型，形成结构数字孪生模型计算网格。
7.第二步：基于结构载荷识别网格模型，利用tiam算法，对考虑不确定性的结构数字模型构造分布式动态载荷的时域识别模型。首先，应用正交多项式拟合结构分布式动态载荷，建立结构在时域内的逆向模型，进而得到结构动响应与正交多项式系数之间的联系，利用结构上有限测点的动响应信息及结构的动态特性实现对正交多项式系数的求取，从而构造完成基于正交多项式展开的确定性结构分布式动态载荷时域识别模型；
8.第三步：根据结构状态识别网格模型识别结构损伤位置信息。结构损伤的位置可以通过损伤前后模态应变能的变化率来确定。通过计算结构损伤前后的模态应变能，计算
出损伤定位指标，即模态应变能变化率(modal strain energy change ratio，msecr)：
[0009][0010]
其中，
[0011][0012]
其中n代表选用的模态阶数，m代表单元数，
h
f
i
代表未损伤结构的第i阶模态频率，
d
f
i
代表损伤结构的第i阶模态频率，公式中的权重系数
h
f
i2
/
d
f
i2
会减小模态频率变化小的模态的影响。具有较大指标的单元被视为损伤单元。
[0013]
第四步：在确定候选损伤单元损伤定位指标的基础上，通过优化方法确定损伤单元中每层的损伤程度和损伤类型，采用优化方法识别损伤程度。通过模式搜索法(hook and jeeves，hj)来求解优化问题。通过搜索谷线或脊线方向来加速优化。重复这两个步骤直到满足精度要求。该方法的整个过程可分为三部分，即未损伤结构的模型修正、损伤定位和损伤量化。
[0014]
第五步：通过构建的结构数字模型网格，对结构上位置载荷作用点上识别的载荷与结构内部各节点损伤状态进行输出，形成物理结构的动态数字孪生模型。同时，考虑结构外部载荷识别更新频率与内部损伤状态更新频率的差异性，建立动态更新模型对不同周期下得到的外部载荷，结构内部损伤状态以及已知传感器获取的动态数据进行同步更新。
[0015]
第六步：在动态数字孪生模型基础上，建立静载作用下含温度影响的飞行器典型结构渐进损伤演化模型，建立针对复合材料的典型的渐进损伤模型的应力分析，损伤判据，以及损伤后材料的性能退化方案。建立静载作用下含温度影响的复合材料渐进损伤仿真方法，利用有限元方法对渐进损伤模型进行损伤演化迭代计算。
[0016]
本发明与现有技术相比的优点在于：
[0017]
本发明提供了一种基于载荷与状态结合的飞行器结构强度数字孪生方法，并涵盖结构载荷识别与结构状态监测的结构动态监测方法，弥补和完善了传统结构健康监测对结构损伤单一变量量化的局限性。所建立的数字孪生样机系统可以更加直观的对飞行器中典型结构的受载情况与结构损伤状态等参数进行同步监测，对结构的损伤破坏可进行虚拟演化与寿命预测。
附图说明
[0018]
图1是本发明针对飞行器结构强度数字孪生的流程图；
[0019]
图2是本发明针对飞行器数字孪生模型损伤演化的流程图。
具体实施方式
[0020]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅为本发明的一部分实施例，而不是全部的实施例，基于本发明中的实施例，本领域的普通技术人员在不付出创造性劳动的前提下所获得的所有其他实施例，都属于本发明的保护范围。
[0021]
如图1所示，本发明提出了一种基于载荷与状态结合的飞行器结构强度数字孪生方法，包括以下步骤：
[0022]
(1)根据实体结构模型，由结构约束条件，确定结构未知外载荷作用点，同时建立结构的数字模型，根据模型几何参数确立结构载荷识别的网格模型与结构状态识别的网格模型。
[0023]
(2)基于结构载荷识别的网格模型，利用tiam算法对考虑不确定性的结构构造分布式动态载荷的识别模型进行分布式动态载荷识别。首先，应用正交多项式拟合结构分布式动态载荷，建立结构在时域内的逆向模型，进而得到结构动响应与正交多项式系数之间的联系，利用结构上有限测点的动响应信息及结构的动态特性实现对正交多项式系数的求取，从而构造完成基于正交多项式展开的确定性结构分布式动态载荷时域识别模型；其中，分布式动态载荷关于正交多项式的表达式为：
[0024]
f(x,t)＝[p0(x) p1(x) p2(x)
ꢀ…ꢀ
p
j
(x)][a0(t) a1(t) a2(t)
ꢀ…ꢀ
a
j
(t)]
t
[0025]
式中，f(x,t)为结构第x个节点在t时刻的动态载荷；[p0(x) p1(x) p2(x)
ꢀ…ꢀ
p
j
(x)]为第x个节点处的j阶正交多项式向量；[a0(t) a1(t) a2(t)
ꢀ…ꢀ
a
j
(t)]
t
表示在t时刻的正交多项式系数向量，记为a(t)。
[0026]
之后，将结构的不确定性参数用区间定量化描述，对正交多项式系数向量在不确定性参数区间中心值处进行一阶泰勒级数展开，进而将不确定性结构的分布式动态载荷识别问题转换为两类确定性问题，即在不确定性参数中心值处正交多项式的系数向量的求解，和正交多项式系数向量关于每一个不确定性参数在其中心值处的灵敏度的计算，从而得到任意时刻正交多项式的系数向量的上、下界值；其中，正交多项式的系数向量的上、下界值按如下的方式求解：假设结构含有m个不确定性参数，记结构的m维不确定性参数的区间向量为b
i
＝[b
l
,b
u
],b
i
∈b
i
＝[b
il
,b
iu
],i＝1,2,
…
,m，其中b
l
为区间变量b
i
的下界，b
u
为区间变量b
i
的上界；b
il
为区间变量b
i
第m个分量b
i
的下界，b
iu
为区间变量b
i
第m个分量b
i
的上界；进而对区间变量b
i
进行如下的变换后得到：
[0027][0028]
其中，不确定性参数的区间半径记为
△
b＝0.5
×
(b
u
‑
b
l
),
△
b
i
＝0.5
×
(b
iu
‑
b
il
),i＝1,2,
…
m，不确定性参数的区间中心值记为b
c
＝0.5
×
(b
u
b
l
),标准区间变量为λ＝[
‑
1，1]，t时刻正交多项式系数向量的上界值a
u
(t)与下界值a
l
(t)可直接通过下列显式求得：
[0029][0030][0031]
式中，a(t,b
c
)为不确定性参数中心值处正交多项式系数向量，为多项式系数向量关于第i个不确定性参数在中心值处的灵敏度，可以利用差分代替微分的方法近似表示，即δb
i
为关于不确定性参数b
i
的摄动；
[0032]
进一步，利用求解的正交多项式系数向量的界值，得到考虑不确定性的结构分布式动态载荷f
i
(x,t,b)的上界值f
u
(x,t)和下界值f
l
(x,t)：
[0033][0034][0035]
式中，p
k
(x)为第x个节点处的第k项正交多项式，a
ku
(t)和a
kl
(t)分别为t时刻第k项正交多项式系数的上界值和下界值；从而完成利用泰勒级数展开法tiam算法对考虑不确定性的结构分布式动态载荷的识别模型。
[0036]
(3)根据结构状态识别网格模型识别结构损伤位置信息。结构损伤的位置可以通过损伤前后模态应变能的变化率来确定。每个单元的模态应变能(modal strain energy，mse)可以按照如下公式计算：
[0037][0038]
其中mse
ij
表示第i阶模态第j个单元的单元应变能，k
j
代表第j个单元的刚度矩阵，代表第i阶模态振型。
[0039]
通过计算结构损伤前后的模态应变能，计算出损伤定位指标，即模态应变能变化率(modal strain energy change ratio，msecr)：
[0040][0041]
其中，
[0042][0043]
其中n代表选用的模态阶数，m代表单元数，
h
f
i
代表未损伤结构的第i阶模态频率，
d
f
i
代表损伤结构的第i阶模态频率，公式中的权重系数
h
f
i2
/
d
f
i2
会减小模态频率变化小的模态的影响。具有较大指标的单元被视为损伤单元。
[0044]
(4)在确定候选损伤单元损伤定位指标的基础上，通过优化方法确定损伤单元中每层的损伤程度和损伤类型，这里采用优化方法识别损伤程度：
[0045]
find
[0046]
min
[0047]
s.t.
[0048]
其中代表第k层的三个损伤参数，m代表待辨识损伤单元个数，l代表层合板的层数，和代表参考模型和实际损伤模型第i阶模态振型，f
r
和f
d
代表参考模型和实际损伤模型的频率向量，μ1和μ2分别代表模态振型和频率对应的加权系数。
[0049]
采用模式搜索法(hook and jeeves，hj)方法来求解优化问题。通过搜索谷线或脊
线方向来加速优化。该方法包括两种类型的移动方式：一种是探索性移动，即对每个坐标方向的探索，以找到一个快速的方向。另一种是模式移动，它是通过一系列探索移动来沿快速下降方向进行加速运动。重复探索性移动和模式移动两种方式直到计算得到满足精度要求的损伤程度结果。方法的流程图如图1所示。该方法的整个过程可分为三部分，即未损伤结构的模型修正、损伤定位和损伤量化。
[0050]
(5)通过构建的结构数字模型，对结构上位置载荷作用点上识别的载荷与结构内部各节点损伤状态进行输出，形成结构的动态数字孪生模型。同时，考虑结构外部载荷识别更新频率与内部损伤状态更新频率的差异性，建立动态更新模型对不同周期下得到的外部载荷，结构内部损伤状态以及已知传感器获取的动态数据进行同步更新。
[0051]
(6)在动态数字孪生模型基础上，建立静载作用下含温度影响的飞行器典型结构渐进损伤演化模型，如图2所示。针对复合材料的典型的渐进损伤模型应主要包含三个部分：应力分析，损伤判据，以及损伤后材料的性能退化方案。在应力分析中，建立组份材料在损伤前、后的本构模型，用以确定材料内部的应力场分布。在损伤判断过程中，引入合适的损伤判据用以判断损伤时刻与损伤类型。在损伤材料的性能退化方案中，建立相应的损伤演化方案，以描述材料在初始损伤发生后的性能衰减规律。在具体的求解中，这三个部分统一在一个完整的循环迭代过程当中，在每一个步长中不断循环求解，直到材料发生最终的整体失效。建立静载作用下含温度影响的复合材料渐进损伤仿真方法，利用有限元方法建立渐进损伤模型迭代计算。
[0052]
综上所述，本发明提出了一种针对复合材料层合板结构在疲劳载荷作用下、在给定的约束条件下，充分考虑结构多失效模式及其相关性的系统可靠性分析方法。首先，根据层合板结构材料以及载荷等情况的具体特征，结合二维hashin准则计算每个失效元在不同失效模式下的失效概率，认为最大失效概率层发生破坏；其次，根据建立的单向层合板退化模型，在疲劳加载过程中对材料属性进行退化，对退化后的模型进行重分析，计算各失效元失效概率，确定主失效路径；最后，对循环加载结束的结构进行剩余强度校核，在此基础上建立疲劳系统可靠性指标，为结构之后的分析和优化提供指导和参考。
[0053]
以上仅是本发明的具体步骤，对本发明的保护范围不构成任何限制；其可扩展应用结构多失效模式的可靠性求解领域，凡采用等同变换或者等效替换而形成的技术方案，均落在本发明权利保护范围之内。
[0054]
本发明未详细阐述部分属于本领域技术人员的公知技术。