一种变体机翼蒙皮变形重构方法、计算机存储介质及设备技术

本发明专利技术属于变体机翼变形测量领域，涉及一种变体机翼蒙皮变形重构方法、计算机存储介质及设备。该方法包括：根据变体机翼的主被动变形，将蒙皮划分为固定区蒙皮与变形区蒙皮。固定区蒙皮采用一般方法重构，结果作为变形区蒙皮重构的边界条件；将变形区的内部支撑化简为边界条件并放置单元节点，利用单元节点的自由度来表征骨架驱动引起的变形；根据仿真应变云图在部分单元内表面布置双层应变传感器；将传感器测值进行弯曲、拉压解耦，构造误差泛函，取变分得到单元的节点自由度，叠加固定区蒙皮的重构结果，插值得到蒙皮的全场变形情况。本发明专利技术能够解决现有技术中蒙皮重构方法难以适用于变体机翼蒙皮重构、重构结果误差较大的技术问题。问题。问题。

【技术实现步骤摘要】
一种变体机翼蒙皮变形重构方法、计算机存储介质及设备


[0001]本专利技术属于变体机翼变形测量领域，涉及一种基于内表面有限应变点测量的变体机翼蒙皮变形重构方法、计算机存储介质及设备，更具体地，涉及一种利用变体机翼蒙皮内表面离散的应变数据分区重构出固定区蒙皮变形与变形区蒙皮变形并叠加得到蒙皮全场变形的重构方法。

技术介绍

[0002]变体机翼是一种通过主动改变自身气动外形、气动特性，以适应不同的飞行情况的智能机翼，它使得环境适应性飞行器成为了一种可能，是未来先进飞行器的重要特征和发展方向之一。区别于传统机翼的铰链式非连续变形结构，目前变体机翼变形主要以翼型变弯度、变厚度和变扭转等中小尺度的变形为主，通过驱动机翼蒙皮的连续变形实现以较小的驱动能耗就能显著改变飞行器整体气动性能和环境适应性的目的。连续变化的蒙皮面相较于传统的铰链式分体变形，能够显著地改变飞行器的升阻比和提高燃油功耗。因此，与现有的对固定翼整体结构进行简化后的变形重构不同，蒙皮实时全场变形的高效高精感知是变体机翼能够有效控制变形量，达到精准控制飞行器气动特性的关键。
[0003]变体机翼蒙皮的结构、变形具有以下特点：整体的蒙皮结构可分为固定区和变形区，如图1所示，固定区蒙皮的变形与固定翼蒙皮一致，为受外载荷引起的被动变形；变形区蒙皮的变形主要为受驱动结构引起的主动变形，变形时需要满足表面切方向的弹性变形来达到所需要的变形量，而在表面法向具有更大的刚性以确保实际飞行过程中不会因外界变化的气压而影响表面流线型，因此变形区蒙皮材料一般为复合材料，且材料在切向与法向需要表现出不同的材料特性；为了实现连续变形，机翼内驱动装置密集，如图2所示，蒙皮的内表面与大量骨架接触，具有复杂的边界条件。现有变形感知方法一般会将机翼整体简化为悬臂板等结构，在分析静载变形时这样的简化可以显著减少计算量，但当考虑气动特性且重构对象为变体机翼时，变形区蒙皮的外形是飞行器调整启动性能以适应环境的关键，不能将外形简单简化。
[0004]变体机翼变形区蒙皮的变形特点要求重构方法能适应复杂材料特性与变形情况，且变体机翼内置的复杂驱动装置使得蒙皮内表面空间狭窄，无法满足在外表面以及大面积布置传感器的空间要求，只能采用内表面粘贴式的有限点接触式测量方式。因此，亟需开发新的基于有限点的接触式变形测量方法，将变体机翼固定区与变形区分区重构并叠加，实现蒙皮全场变形重构。

技术实现思路

[0005]针对现有技术的以上缺陷或改进需求，本专利技术提供了一种变体机翼蒙皮变形重构方法、计算机存储介质及设备，其目的在于，解决现有技术中蒙皮重构方法难以适用于变体机翼蒙皮重构、重构结果误差较大的技术问题。
[0006]为实现上述目的，按照本专利技术的一个方面，提供了一种变体机翼蒙皮变形重构方
法，包括如下步骤：
[0007]步骤1：根据变体机翼的主被动变形，将整体机翼蒙皮划分为固定区蒙皮与变形区蒙皮，固定区蒙皮的小变形通过化简为悬臂板进行重构，并将固定区蒙皮的重构结果作为变形区蒙皮大变形重构时的固定边界条件；
[0008]步骤2：根据变形区蒙皮内部的驱动骨架结构特点，将内部骨架结构化简为边界条件并放置单元节点，使单元节点与驱动骨架重合，利用单元节点的自由度来表征骨架驱动引起的变形；
[0009]步骤3：对变形区蒙皮进行变形仿真，以获取该结构在真实运行过程中的整体应变分布云图；
[0010]步骤4：根据应变云图的分布情况在应变最大值超过阈值的变形区蒙皮单元内表面布置双层应变传感器测量；
[0011]步骤5：在变体机翼进行自主变形时，变形区蒙皮的变形形式包含拉压变形和弯曲变形的耦合，因此将双层传感器测得的应变数据进行弯曲、拉压应变的解耦，并利用解耦得到的实测应变与节点自由度导出的理论应变构造误差泛函，取变分，得到的驻点为单元局部坐标系内的节点自由度，再通过插值方式得到变形区蒙皮的变形情况，在边界位置叠加固定区蒙皮的变形，得到变体机翼蒙皮整体结构的全场变形情况；
[0012]其中，单元局部坐标系内的节点自由度u
e
求解如下：
[0013][0014]上式中，K
e
是一个仅与单元中传感器布点位置和单元形函数有关的函数，f
e
是一个与单元中传感器布点位置、单元形函数和实际测量值有关的函数；B
m
、B
b
和B
s
分别为理论膜应变矩阵、理论曲率矩阵、理论横向剪切应变矩阵；h为蒙皮中性面和双层应变传感器中的上应变片的距离；ω为虚拟布置系数，λ为横向剪切应变惩罚系数；e
ε
、κ
ε
分别为蒙皮中性面处的测量膜应变与测量曲率。
[0015]进一步地，
[0016]步骤2中，利用单元节点的自由度来表征骨架驱动引起的变形的方法如下：
[0017]根据Mindlin中厚板理论构建四节点逆壳单元，每个节点分别具有沿x y z轴方向的三个平动自由度u v w和转动自由度θ
x θ
y θ
z
，则每个节点的位移向量表示为：
[0018][0019]整个单元的位移向量表示为：
[0020][0021]上式中，上标e表示在单元内分析；下表i表示节点编号，i＝1，2，3，4。
[0022]进一步地，
[0023]将单元中性面上任意一点(x，y)的膜变形通过形函数节点插值的形式表示为：
[0024][0025]上式中，N
k
，L
i
，M
i
为单元插值形函数，其具体形式为：
[0026][0027][0028]其中，
[0029][0030]s，t为某一节点在局部坐标系中的等参坐标，假设单元的长宽尺寸为2a
×
2b则：
[0031][0032]则单元内任意一点(x，y，z)的理论位移表示为：
[0033][0034]进一步地，
[0035]根据一阶剪切变形理论，用式(7)的理论位移表示应变，得到单元的理论应变如下：
[0036][0037][0038][0039]上式中ε、κ、γ分别表示单元的线应变、曲率和切应变，下标x0、y0、x
y
表示x方向、y方向和xy方向，分别表示对x方向和y方向求偏导数，则e(u
e
)、κ(u
e
)、g(u
e
)分别表示单元中性面的理论膜应变向量、理论曲率向量和理论横向剪切应变向量。
[0040]进一步地，
[0041][0042]其中，k＝i＝1，2，3，4。
[0043]进一步地，
[0044]步骤5中，基于机翼变形区蒙皮内表面单元内双层应变传感器获得内表面相应位置的实测应变值：
[0045][0046]上式中上标
+
‑
分别表示双层应变传感器的上应变片和下应变片，其中ε0°
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种变体机翼蒙皮变形重构方法，其特征在于，包括如下步骤：步骤1：根据变体机翼的主被动变形，将整体机翼蒙皮划分为固定区蒙皮与变形区蒙皮，固定区蒙皮的小变形通过化简为悬臂板进行重构，并将固定区蒙皮的重构结果作为变形区蒙皮大变形重构时的固定边界条件；步骤2：根据变形区蒙皮内部的驱动骨架结构特点，将内部支撑结构化简为边界条件并放置单元节点，使单元节点与驱动骨架重合，利用单元节点的自由度来表征骨架驱动引起的变形；步骤3：对变体机翼进行变形仿真，以获取变形区蒙皮结构在真实运行过程中的整体应变分布云图；步骤4：根据应变云图的分布情况在应变最大值超过阈值的变形区蒙皮单元内表面布置双层应变传感器测量；步骤5：在变体机翼进行自主变形时，变形区蒙皮的变形形式包含拉压变形和弯曲变形的耦合，因此将双层传感器测得的应变数据进行拉压应变、弯曲的解耦，并利用解耦得到的实测应变与节点自由度导出的理论应变构造误差泛函，取变分，得到的驻点为单元局部坐标系内的节点自由度，再通过插值方式得到变形区蒙皮的变形情况，在边界位置叠加固定区蒙皮的变形，得到变体机翼蒙皮整体结构的全场变形情况；其中，单元局部坐标系内的节点自由度u
e
求解如下：上式中，K
e
是一个仅与单元中传感器布点位置和单元形函数有关的函数，f
e
是一个与单元中传感器布点位置、单元形函数和实际测量值有关的函数；B
m
、B
b
和B
s
分别为理论膜应变矩阵、理论曲率矩阵、理论横向剪切应变矩阵；h为蒙皮中性面和双层应变传感器中的上应变片的距离；ω为虚拟布置系数，λ为横向剪切应变惩罚系数；e
ε
、κ
ε
分别为蒙皮中性面处的测量膜应变与测量曲率。2.如权利要求1所述的变体机翼蒙皮变形重构方法，其特征在于，步骤2中，利用单元节点的自由度来表征骨架驱动引起的变形的方法如下：根据Mindlin中厚板理论构建四节点逆壳单元，每个节点分别具有沿xyz轴方向的三个平动自由度uvw和转动自由度θ
x
θ
y
θ
z
，则每个节点的位移向量表示为：整个单元的位移向量表示为：上式中，上标e表示在单元内分析；下表i表示节点编号，i＝1，2，3，4。3.如权利要求2所述的变体机翼蒙皮变形重构方法，其特征在于，将单元中性面上任意一点(x，y)的膜变形通过形函数节点插值的形式表示为：
上式中，N
k
，L
i
，M
i
为单元插值形函数，其具体形式为：为单元插值形函数，其具体形式为：
其中，j＝1,2,3,4且j≠i；s，t为某一节点在局部坐标系中的等参坐标，假设单元的长宽尺寸为2a
×
2b则：则单元内任意一点(x，y，z)的理论位移表示为：4.如权利要求3所述的变体机翼蒙皮变形重构方法，其特征在于，根据一阶剪切变形理论，用式(7)的理论位移表示应变，得到单元的理论应变如下：根据一阶剪切变形理论，用式(7...

【专利技术属性】
技术研发人员：冀晶晶，朱家烽，黄永安，张帆，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：