具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，本发明提供了一种用于飞行器模态试验的关键参数识别方法，包括以下步骤：

航天飞行器的质量分布情况、部段直径、蒙皮厚度、部段长度、对接面刚度等参数会对全箭的动特性预示精度不同，其偏差对结果的影响也有很大的区别。将这些参数统称为结构参数。结构参数也是直接表达在结构动力学模型中的参数。结构参数与动特性参数性质上有很大区别：结构参数是结构局部特性的表现，数量很多；而动特性参数则是结构的整体动特性的表现，数量有限。由动特性参数确定结构参数是采用少量参数来估测大量参数，可靠性不高，而由结构参数来估测动特性参数有影响，但是理论上来说也是更加的可靠的，因此本实施例将飞行棋模态试验的分析重点调整至结构参数上，可保证飞行器模态试验能够更准确的描述飞行器的动特性。

步骤S1、将预示准确度和对动特性影响的敏感程度设定为分类双指标，并基于分类双指标对飞行器模态试验的结构参数进行分类；

步骤S1中，分类双指标设定结构参数分类的方法包括：

设定预示准确度阈值，将结构参数的预示准确度与预示准确度阈值比较，将结构参数基于预示准确度阈值进行二分类为预示准确度高的类别和预示准确度低的类别；

设定敏感程度阈值，将结构参数的预示准确度与敏感程度阈值比较，将结构参数基于敏感程度阈值进行二分类为对动特性影响的敏感程度高的类别和对动特性影响的敏感程度低的类别；

将结构参数基于预示准确度阈值进行二分类与基于敏感程度阈值进行二分类相结合实现将结构参数进行四分类得到A类别、B类别、C类别和D类别。

所述结构参数的四分类方法包括：

选取多个结构参数作为模型训练参数，并将所述模型训练参数的预示准确度和对动特性影响的敏感程度分别与所述预示准确度阈值和敏感程度阈值进行比较，其中，

若所述模型训练参数的预示准确度大于所述预示准确度阈值且对动特性影响的敏感程度大于敏感程度阈值，则所述模型训练参数标定为C类别；

若所述模型训练参数的预示准确度大于所述预示准确度阈值且对动特性影响的敏感程度小于敏感程度阈值，则所述模型训练参数标定为D类别；

若所述模型训练参数的预示准确度小于所述预示准确度阈值且对动特性影响的敏感程度大于敏感程度阈值，则所述模型训练参数标定为A类别；

若所述模型训练参数的预示准确度小于所述预示准确度阈值且对动特性影响的敏感程度小于敏感程度阈值，则所述模型训练参数标定为B类别；

选取多个已知类别的模型训练参数分别作为模型训练样本，其中，模型训练参数的预示准确度和对动特性影响的敏感程度作为模型训练样本的样本特征，模型训练参数的类别作为模型训练样本的样本标签；

将所述模型训练样本运用到贝叶斯分类器进行四分类训练得到参数分类模型，所述参数分类模型用于结构参数自动识别类别，其中，所述模型训练样本的样本特征作为贝叶斯分类器的输入，所述模型训练样本的样本标签作为贝叶斯分类器的输出；

将待分类的结构参数的预示准确度和对动特性影响的敏感程度输入至参数分类模型，输出待分类的结构参数的类别。

以上可将结构参数的分类有人为的阈值分类转变为模型自动分类，提高分类效率和精度，特别是在需要分类的结构参数数量大的情况下，优势更为明显。

如图2所示，结构参数四分类的结果包括A类别、B类别、C类别和D类别，其中，

A类别表征为结构参数具有预示准确度低且对动特性影响的敏感程度高的属性；

B类别表征为结构参数具有预示准确度低且对动特性影响的敏感程度低的属性；

C类别表征为结构参数具有预示准确度高且对动特性影响的敏感程度高的属性；

D类别表征为结构参数具有预示准确度高且对动特性影响的敏感程度低的属性。

具体的，典型的D类别的参数是壁厚较厚，结构简单部段的刚度参数，比如说固体发动机推力室的刚度特性。其刚度较大，对动特性影响不大，并且在振动过程中基本符合平截面假定，实际表现出的刚度与理论解十分接近。

典型的B类别的参数有大量对接螺栓连接的部段对接面刚度。该部段对接刚度很大，结构振动在此处响应很小，但是，考虑到螺栓预紧力、端框出的复杂局部结构，端框对接面接触情况等，此处的刚度预示十分困难。

C类别的参数举例为结构的质量，以及长度、直径等几何参数，这些参数对动特性影响十分敏感，但是却能够十分直接的获得或预示。

A类别的参数有蒙皮较薄部段的弯曲刚度，点式连接的部段对接面刚度等，这些参数即难以获取又会对动特性结果有较大影响。

步骤S2、基于步骤S1中结构参数的分类结果进行针对性的模态试验，并识别出模态试验的关键参数。

模态试验需要获得动特性预示结果，并保证动特性预示结果在预示偏差范围内，半实物模态试验的目的也是如此。在完成了对结构参数的四分类后可将结构参数分成A、B、C和D四个类别，即可建立考虑偏差的结构动力学模型，对落在A类别内的结构参数确定为关键参数考虑进行实物试验，将其从A区转至C区，以保证最终结构动力学模型的动特性预示结果偏差满足要求。

步骤S2，基于结构参数的分类结果进行针对性的模态试验的方法包括：

利用模态试验模型修正方法在结构参数确定出属于A类别的关键参数；

建立考虑偏差的结构动力学模型，利用结构动力学模型对属于A类别的关键参数进行实物试验，实现将关键参数从A类别转化至C类别，使得关键参数的预示准确度由低提至高，并使得结构动力学模型输出的动特性预示结果在预示偏差范围内。

实物试验的目的是提高对动特性影响的敏感程度高且预示准确度低的关键参数的预示准确度，以提高结构动力学模型的动特性预示精度。在结构参数内选择关键参数在很大程度上需要对物理结构有深刻的认识，所选关键参数应该是在物理上存在疑问并且与结构动力学模型特征相关的，实际上，对同一可疑点，通常有许多可供选择的修正参数。一般来说，应该在结构参数中选择对动特性敏感程度高，但是，反过来则是不对的。也就是说，动特性对某个结构参数敏感，并不是选择此参数作为实物试验的关键参数的充分条件。总之，在结构参数中选择实物试验的关键参数时应选择那些在结构动力学模型中已识别的不确定因素对应的结构参数，并且该结构参数对动特性影响的敏感程度高。

利用模态试验模型修正方法确定出关键参数的方法包括：

模态试验模型修正方法，可选择的方法有特征值方程平衡方法，子结构能量函数识别方法，最有子空间法，灵敏度抽样法等，均在本发明的保护范围内，本实施例为描述方便以特征值方程平衡方法进行关键参数的确定过程描述。

基于结构参数的偏差，建立的特征方程，特征方程为：

；

式中，为结构刚度矩阵，为矩阵增量，为质量矩阵，矩阵增量，矩阵增量和要保证方程的平衡，为模态振型矩阵，为模态质量阵；

将增量效应综合在一起，得到局部化的矩阵L，所述矩阵L公式为：

；

式中,模态振型矩阵包含无偏差状态对应的模态振型；为结构刚度矩阵， 为矩阵增量，为质量矩阵，矩阵增量，矩阵增量和要保证方程的平衡，为模态 振型矩阵，为模态质量阵；

通过检查局部化向量q，选出主要的建模误差，所述建模误差公式为：

；

式中，气权重数反应第h阶模态测试数据和分析数据之间的一致性程度，L为局 部化的矩阵，为L矩阵中的第i行第h列的数值，N为矩阵总数，为第i个建模误差，n为建 模误差总数；

设定误差阈值，将高于误差阈值的对应的结构参数重新利用结构动力学模型进 行动特性预示结果的偏差评估。

重新利用结构动力学模型进行动特性预示结果的偏差评估的方法包括：

若动特性预示结果在预示偏差范围内，则将高于误差阈值的对应的结构参数确 定为属于A类别的关键参数；

若动特性预示结果不在预示偏差范围内，则降低误差阈值增加高于误差阈值的 对应的结构参数的数量。

实物试验用于提高对动特性敏感程度高且预示准确度低的关键参数的预示准确度，以提高结构动力学模型的动特性预示精度。

高值相关的结构参数进行实物试验可有效提高结构动力学模型的动特性预示精 度，因此可设定一个误差阈值，用于将高于误差阈值的值定义为高高值，从而可获得高 值对应的结构参数重新利用结构动力学模型进行动特性预示结果的偏差评估，确定出高 值对应的结构参数是否可作为关键参数，若高值对应的结构参数无法作为关键参数，可对 误差阈值进行相对降低，使得高于误差阈值的值数量增加，即高值数量增加，则相应的 高值对应的结构参数数量也增加，再重新利用结构动力学模型进行动特性预示结果的偏 差评估，直至关键参数得以确定。

所述关键参数具有与A类别中参数相同的属性，所述关键参数为飞行器关键结构特性中的表现，以及所述飞行器模态试验的重点关注点，并且所述关键参数在由A类别转至C类别后通过结构动力学模型输出的动特性预示结果的精度得以提高，以保证利用关键参数进行的模态试验能够更准确的描述飞行器的动特性。

结构参数为飞行器结构特性的表现，以及直接表达在结构动力学模型中的参数，将模态试验关注点转化至利用结构特性评估飞行器动特性。

本发明利用由预示准确度和对动特性影响的敏感程度设定的分类双指标进行结构参数四分类，并通过在结构参数识别关键参数，实现对模态试验的精细化设计，将模态试验的关注点由整体动特性转向结构局部特征，可实现在结构有局部调整时，快速且准确的掌握结构动力学模型的参数，可通过局部模型的更换来获得新状态的动特性，无需重新进行模态试验，从而扩充了模态试验的灵活性。

以上实施例仅为本申请的示例性实施例，不用于限制本申请，本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内，对本申请做出各种修改或等同替换，这种修改或等同替换也应视为落在本申请的保护范围内。