[发明专利]静气动弹性低动态试验数据提取方法、装置、介质及设备

说明书

本发明公开了静气动弹性低动态试验数据提取方法、装置、介质及设备，该方法包括：步骤1：获得机翼的多组应变数据；步骤2：获得多个首次聚类的类心；步骤3：分别计算每组应变数据与每个聚类类心之间的距离；步骤4：针对每组应变数据，将应变数据归属到距离的最小值对应的聚类中；步骤5：针对每个聚类，判断聚类中的每组应变数据是否为无效的应变数据，并删除；步骤6：重新计算每个聚类的类心；步骤7：针对每组无效的应变数据，重复执行步骤3‑6，当再次判断无效的应变数据无效时，将最后一次删除无效的应变数据后重新计算的聚类的类心作为对应风速下的机翼应变量。本发明得到不同风速下的机翼应变，对于柔性飞行器静气动弹性试验尤为重要。

技术领域

本发明属于数据处理技术领域，具体涉及一种静气动弹性低动态试验数据提取方法、装置、电子设备及介质。

背景技术

现代飞行器设计对静气动弹性特性要求越来越严苛。风洞静气动弹性试验技术是测量飞行器不同飞行速度下机翼载荷与气动弹性变形，得到飞行器配平数据的重要手段。随着国内大展弦比、长航时飞行器的研制，机翼刚度越来越小，静气动弹性效应愈发明显。

现有技术中，用于飞行器风洞静气动弹性试验的数据处理方法多针对常规飞行器，常规飞行器机翼刚度大变形小、风洞试验中机翼基本不出现微幅振荡。可以采用均值估计的简化手段得到不同风速下机翼应变、根据经验估计机翼载荷、控制风速。

然而高空长航时飞行器多要求采用大展弦比柔性机翼，对静气动弹性特性提出了更高的要求。此类飞行器在风洞静气动弹性试验中机翼变形大，由于扰动机翼容易出现微幅振动，且柔性机翼微幅振动响应无法与风速变化完全同步。因此，传统数据处理方法的精度难以保证，已经不再适用；此类柔性机翼对载荷敏感性强，风速失之毫厘就会出现差之千里的破坏，根据经验控制风速的方法也难以适用。

因此，特别需要一种能够排除柔性机翼微幅振动干扰，大大降低由于柔性机翼和风速变化非同步响应所带来的误差的静气动弹性低动态试验数据提取方法。

发明内容

本发明的目的是提出一种能够排除柔性机翼微幅振动干扰，大大降低由于柔性机翼和风速变化非同步响应所带来的误差的静气动弹性低动态试验数据提取方法。

第一方面，本发明提供一种静气动弹性低动态试验数据提取方法，包括：步骤1：在每个预设攻角下，在预设时段内分k次逐渐提高风速，获得机翼的多组应变数据；步骤2：对所述机翼的多组应变数据进行聚类，获得多个首次聚类的类心；步骤3：分别计算每组所述应变数据与每个聚类类心之间的距离；步骤4：针对每组所述应变数据，比较所述应变数据与每个聚类类心之间的距离，获取所述距离的最小值，将所述应变数据归属到所述距离的最小值对应的聚类中；步骤5：针对每个聚类，计算所述聚类中应变数据的平均值，基于所述平均值分别判断所述聚类中的每组应变数据是否为无效的应变数据，并将无效的应变数据从所述聚类中删除；步骤6：重新计算每个聚类的类心；步骤7：针对每组无效的应变数据，重复执行步骤3-6，当再次判断所述无效的应变数据无效时或迭代次数等于预设次数时，将最后一次删除无效的应变数据后重新计算的聚类的类心作为对应风速下的机翼应变量。

可选的，采用下述公式获得首次聚类的类心：

C_j＝1/2(ε_j+ε_j+1)

j＝1,2,…,k

其中，C_j为第j个聚类的类心，ε_j为机翼的第j组应变数据，k为风速变化次数，n为应变数据的组数，j+1n。

可选的，采用下述公式计算每组所述应变数据与每个聚类类心之间的距离：

其中，d(ε_l,C_j)为第l组应变数据与第j个聚类类心的距离，ε_l为机翼的第l组应变数据。