基于风洞试验数据与CFD计算的无人机控制矩阵的修正方法与流程

本发明涉及一种基于风洞试验数据与CFD计算的无人机控制矩阵的修正方法，同时兼顾风洞实验数据及CFD计算结果进行参考，能够在工程上快速实现常规布局和非常规布局飞行控制矩阵的计算，主要在飞行器气动设计过程中使用，属于航空飞行器气动设计技术领域。

背景技术：

飞行器设计过程中，控制率的设计通常需要气动协作部门提供，一旦飞行控制矩阵提供误差较大，会导致设计周期的延长和设计的精度的下降。而飞行控制矩阵，主要是由气动部门提供的，传统的经验公式的利用局限于常规的气动布局，对于非常规布局及近似常规布局的外形计算上存在一定的误差，且随着未来气动布局(基于气动设计优化得到的外形)的发展误差会越来越大。因此有必要在满足工程设计精度和时间周期要求的基础上，对传统经验公式进行必要的修正。

技术实现要素：

本发明的技术解决问题是：克服现有技术的不足，提供了一种基于风洞试验数据与CFD计算的无人机控制矩阵的修正方法，实现了无人机控制矩阵工程方法应用，可靠快速地实现了风洞试验数据与CFD计算数据与传统经验公式，最大程度满足飞行控制矩阵工程精度和周期要求。

本发明的技术解决方案是：基于风洞试验数据与CFD计算的无人机控制矩阵的修正方法，步骤如下：

(1)根据无人机的外形参数构建数学模型，并进行CFD计算，得到无人机的气动参数；

(2)判断数学模型是否属于常规气动布局，若是进入步骤(3)，若不是进入步骤(4)；

(3)根据上述确定的数学模型，利用传统的经验公式计算控制矩阵，并用(1)中气动参数中的C_L,C_D以及C_m结合传统的经验公式，对控制矩阵中的相应元素进行修正，得到新的控制矩阵，利用该控制矩阵调整无人机PID控制参数；所述C_L为飞行器升力系数，C_D飞行器阻力系数，C_m为飞行器俯仰力矩系数；

(4)根据上述确定的数学模型，选择利用传统的经验公式或datcom计算控制矩阵；对基于步骤(1)中数学模型加工得到的模型进行风洞试验，利用风洞试验数据中的C_Lα,C_mα以及对控制矩阵中的相应元素进行修正，得到新的控制矩阵，利用该控制矩阵调整无人机PID控制参数；

所述C_Lα飞行器全机升力系数对迎角导数，C_mα飞行器全机俯仰力矩系数对迎角导数，因迎角变化引起的飞机升力系数导数，因迎角变化引起的飞机俯仰力矩系数导数。

当能够获取无人机的外场实验数据时，利用外场实验数据对(3)或(4)得到的控制矩阵进行修正，得到最终控制矩阵。

在步骤(1)中对建立的数学模型进行简化，得到翼身融合体加尾翼的简化数学模型。

对简化的数学模型进行cfd计算，判断升阻比、力矩特性、升力系数、阻力系数是否满足设计指标，若满足，则对简化的数学模型继续执行步骤(2)；否则，重新对数学模型进行简化。

所述的简化处理为将无人机中机身长度方向小于等效弦长的10％的凸起部件进行光滑处理，或者将V型尾翼修改成T型尾翼或者平尾加双立尾或者平尾加单立尾；或者将尾撑杆去除；或者将进气道唇口与机身融合处理；或者将尾喷口与机身融合处理；当翼梢小翼面积小于机翼面积10％时，将翼梢小翼直接去除，或者当翼梢小翼面积大于等于机翼面积10％时，将翼梢小翼去除后，在机翼前缘线和后缘线不变的情况下，增加展长，机翼面积增加量为翼梢小翼面积；或者将外挂物或螺旋桨忽略处理。

气道唇口与机身融合处理为融合表面的引导线的二阶导数连续。

尾喷口与机身融合处理为尾喷口与机身后缘为封闭结构，避免拉维尔喷管效应。

采用等效投影面积法将V型尾翼修改成T型尾翼或者平尾加双立尾或者平尾加单立尾。

等效投影面积法为将V型尾翼投影到水平面和铅直面，得到投影后的两个面积，使修改后段翼型在水平面和铅直面的投影面积与其相等。

当采用V型尾翼修改为T型尾翼时，无人机的使用环境必须为中高空环境。

本发明与现有技术相比的有益效果是：

(1)通过基于风洞试验数据与CFD计算的无人机控制矩阵的修正方法，有效地将气动设计与控制参数设计(PID参数)相结合，避免了以往飞行器设计研制过程中气动设计与控制参数设计联系脱节现象。针对飞行器气动布局不断推陈出新，新布局不断涌现，原有单纯的传统方法估算精度越来越低的现状，基于风洞试验数据与CFD计算的无人机控制矩阵的修正方法，有效地利用新手段借助于风洞试验，CFD数据，以及datcom多维度的对精度进行极大地提升。

(2)本文的简化手段相比于其他简化手段，有效地保证了CFD计算过程中的参数的精度，以及利用datcom软件计算过程中的精度。

(3)由于V型尾翼设计便于控制系统进行设计，未来气动布局的设计中出现的频率会显著增加。本发明根据使用环境(飞行高度)，发动机位置，以及尾翼与主翼相对位置，以及是否飞翼布局，来确定是否修改成T型尾翼或者平尾加双立尾或者平尾加单立尾的。

(4)本发明创造性的将已有风洞数据库中的数据引入非常规布局的控制矩阵计算中，将原有方法精度得到提高。

(5)根据外场的飞行试验，以及自由飞实验数据，获得的气动参数对飞行控制矩阵进行进一步修正，进而获得控制矩阵。未来可用于对该型号飞行器的改进型，无人机PID控制参数设计使用。

附图说明

图1为本发明流程图。

具体实施方式

下面结合附图对本发明的具体实施方式进行进一步的详细描述。

飞行器气动设计过程中飞行控制矩阵，通常需要精确提供用于飞行器PID设计，为保证设计周期和工程设计的进度，通常传统上采用的手段是简单的经验公式估算。在使用过程中，本发明创造性的将风洞数据和CFD数据引入其中，作为修正非常规布局和常规气动外形的依据。

如图1所示，一种基于风洞试验数据与CFD计算的无人机控制矩阵的修正方法，实施步骤如下：

(1)根据无人机的外形参数构建数学模型，并进行CFD计算，得到无人机的气动参数，例如C_L,C_D以及C_m；所述C_L为飞行器升力系数，C_D飞行器阻力系数，C_m为飞行器俯仰力矩系数；

在保证飞行器设计指标参数的要求下，为了便于对无人机控制矩阵进行求解，为PID控制设计提出控制矩阵。本发明对建立的数学模型进行简化，得到翼身融合体加尾翼的简化数学模型，进而对简化的数学模型进行CFD计算，判断升阻比、力矩特性、升力系数、阻力系数是否满足设计指标，若满足，则对简化的数学模型继续执行步骤(2)；否则，重新对数学模型进行简化。

上述提及的简化可以采用下述方式：

A将无人机中机身长度方向小于等效弦长的10％的凸起部件进行光滑处理，

B将V型尾翼修改成T型尾翼或者平尾加双立尾或者平尾加单立尾；

C将尾撑杆去除；

D将进气道唇口与机身融合处理；

E将尾喷口与机身融合处理；

F当翼梢小翼面积小于机翼面积10％时，将翼梢小翼直接去除，或者当翼梢小翼面积大于等于机翼面积10％时，将翼梢小翼去除后，在机翼前缘线和后缘线不变的情况下，增加展长，机翼面积增加量为翼梢小翼面积；

G将外挂物或螺旋桨忽略处理。

其中，气道唇口与机身融合处理为融合表面的引导线的二阶导数连续。尾喷口与机身融合处理为尾喷口与机身后缘为封闭结构，避免拉维尔喷管效应。

本发明采用等效投影面积法将V型尾翼修改成T型尾翼或者平尾加双立尾或者平尾加单立尾。等效投影面积法为将V型尾翼投影到水平面和铅直面，得到投影后的两个面积，使修改后段翼型在水平面和铅直面的投影面积与其相等。对简化以后得到的尾翼，在保证等效投影面积法下，相对位置本发明中不做特殊要求，按照常规要求完成简化即可。

当采用V型尾翼修改为T型尾翼时，无人机的使用环境必须为中高空环境。

(2)判断数学模型是否属于常规气动布局，若是进入步骤(3)，若不是进入步骤(4)；上述常规气动布局为主翼在前水平尾翼在后，有一个或两个垂尾，例如，ARJ21。

(3)根据上述确定的数学模型，利用传统的经验公式(具体参见飞机设计手册)计算控制矩阵，并用(1)中气动参数中的C_L,C_D以及C_m对控制矩阵中的相应元素进行修正，对C_L,C_D以及C_m的参数项进行替代，得到新的控制矩阵，利用该控制矩阵调整无人机PID控制参数，所述C_L为飞行器升力系数，C_D飞行器阻力系数，C_m为飞行器俯仰力矩系数；

(4)根据上述确定的数学模型，选择利用传统的经验公式或datcom计算控制矩阵；对基于简化数学模型加工得到的模型进行风洞试验，利用风洞试验数据中的C_Lα,C_mα以及对控制矩阵中的相应元素进行修正，对C_Lα,C_mα以及的参数项进行替代，得到新的控制矩阵，利用该控制矩阵作为原始输入参数应用于无人机PID控制参数的设计；所述C_Lα飞行器全机升力系数对迎角导数，C_mα飞行器全机俯仰力矩系数对迎角导数，因迎角变化引起的飞机升力系数导数，因迎角变化引起的飞机俯仰力矩系数导数。

当能够获取无人机的外场实验数据时，利用外场实验数据对(3)或(4)得到的控制矩阵进行修正，得到最终控制矩阵。未来可用于对该型号飞行器的改进型，无人机PID控制参数设计使用。

本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。