CN115659748B 一种基于cfd技术的飞行器姿态控制律设计方法 （南京航空航天大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN115659748B(45)授权公告日2025.07.04

(21)申请号202211354379.7

(22)申请日2022.10.31

(65)同一申请的已公布的文献号申请公布号CN115659748A

(43)申请公布日2023.01.31

(73)专利权人南京航空航天大学

地址210007江苏省南京市秦淮区御道街

29号

(72)发明人邓双厚秦帆支豪林肖天航

(74)专利代理机构江苏圣典律师事务所32237专利代理师周宁

(51)Int.CI.

GO6F30/23(2020.01)

GO6F30/28(2020.01)

B64F5/00(2017.01)

GO6F113/08(2020.01)

GO6F119/14(2020.01)

(56)对比文件

CN106019930A,2016.10.12

CN108595755A,2018.09.28

审查员许哲

权利要求书4页说明书10页附图4页

(54)发明名称

一种基于CFD技术的飞行器姿态控制律设计方法

(57)摘要

CN115659748B本发明公开了一种基于CFD技术的飞行器姿态控制律设计方法，属于计算流体力学、飞行动力学和控制技术领域，缩短了飞行试验周期并降低成本。本发明借助CFD仿真技术对飞行器进行静、动态气动特性计算，获得静稳定性和动稳定性导数；根据静、动导数得到飞行器纵向或横侧向运动舵偏角和姿态角之间的传递函数，并设计控制器实现飞行器俯仰、滚转和偏航姿态的控制；运用CFD与控制的耦合仿真技术获得飞行器姿态角与舵偏角随时间变化曲线，验证控制参数是否满足设计要求。若不满足则进一步通过系统辨识得到修正后的飞行器系统传递函数和控制参数，并运用CFD

CN115659748B

开始

开始

根据工程对象模型的外形划分网格

运用CFD技术进行静、动态气动特性计算

配平状态下静稳定性和动稳定性导数

飞行器纵向或横侧向运动能偏角与姿态角的传递函数

设计姿态运动控制器，获得初始控制参数一

运用CFD/控制耦合仿真技术进行仿真验证

飞行器姿态角和舵偏角变化曲线

S

是备满足设计双求

获得满足设计要求的控制参数

结束

系统辨识获得修正后的传递函数和控制参数

CN115659748B权利要求书1/4页

2

1.一种基于CFD技术的飞行器姿态控制律设计方法，其特征在于，包括以下步骤：

S10、根据工程对象模型的几何外形对其进行网格划分；

S20、运用CFD技术求解模型在不同姿态角和舵偏角下的气动力和力矩系数，并计算初始配平状态下的静稳定性导数和动稳定性导数；

S30、根据步骤S20所得静、动稳定性导数，获得飞行器纵向或横侧向运动舵偏角和姿态角之间的传递函数；

S40、针对飞行器运动舵偏角和姿态角的传递函数，设计姿态运动控制器并设定控制目标，通过控制系统仿真获得满足设计要求的控制参数，作为初始控制参数；

S50、运用CFD与控制系统耦合仿真技术对控制参数进行验证，获得飞行器姿态角和舵偏角随时间变化的曲线，以及非定常的气动力、力矩和流场特征的变化情况；

S60、判断飞行器舵偏角和姿态角随时间变化曲线是否满足设计要求，是则获得达到设计要求的控制参数并结束，否则执行步骤S70;

S70、对飞行器姿态角随舵偏角变化的响应曲线进行系统辨识，获得修正后的传递函数并得到相应的控制参数，继续执行步骤S50,直到满足设计要求。

2.根据权利要求1所述的基于CFD技术的飞行器姿态控制律设计方法，其特征在于，步骤S20具体为：采用CFD技术数值求解定常粘性RANS方程，计算飞行器在不同姿态角和舵偏角下的气动力和力矩；根据要求的初始姿态角，采用多项式曲线拟合的方式获取配平舵偏角以及配平状态下的气动力、力矩和静稳定性导数；通过迭代求解非定常RANS方程得到飞行器做强迫振荡运动时的气动力随时间变化情况，采用傅里叶级数进行拟合，进而辨识出俯仰、偏航和滚转方向的动导数。

3.根据权利要求1或2所述的基于CFD技术的飞行器姿态控制律设计方法，其特征在于，计算飞行器在不同姿态角和舵偏角下的气动力和力矩具体为：

在连续介质假设下，基于质量、动量和能量守恒定律，建立直角坐标系下三维雷诺平均N-S方程如下：

(1)

其中Q为守恒变