CN105468814B - 一种考虑全机气动干扰的直升机飞行特性计算方法 - Google Patents
一种考虑全机气动干扰的直升机飞行特性计算方法 Download PDFInfo
- Publication number
- CN105468814B CN105468814B CN201510716964.0A CN201510716964A CN105468814B CN 105468814 B CN105468814 B CN 105468814B CN 201510716964 A CN201510716964 A CN 201510716964A CN 105468814 B CN105468814 B CN 105468814B
- Authority
- CN
- China
- Prior art keywords
- helicopter
- rotor
- component
- aerodynamic interference
- trim
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 43
- 230000008878 coupling Effects 0.000 claims abstract description 41
- 238000010168 coupling process Methods 0.000 claims abstract description 41
- 238000005859 coupling reaction Methods 0.000 claims abstract description 41
- 230000004044 response Effects 0.000 claims abstract description 32
- 239000002245 particle Substances 0.000 claims abstract description 31
- 238000000205 computational method Methods 0.000 claims abstract description 26
- 238000013461 design Methods 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims description 36
- 230000008859 change Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 230000006698 induction Effects 0.000 claims description 8
- 230000006641 stabilisation Effects 0.000 claims description 8
- 238000011105 stabilization Methods 0.000 claims description 8
- 238000002360 preparation method Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000004321 preservation Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 abstract description 12
- 230000001133 acceleration Effects 0.000 abstract description 4
- 230000005405 multipole Effects 0.000 abstract description 3
- 230000006872 improvement Effects 0.000 description 5
- 241001672694 Citrus reticulata Species 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Traffic Control Systems (AREA)
- Navigation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开的一种考虑全机气动干扰的直升机飞行特性计算方法,包括耦合气动干扰的直升机配平计算方法、耦合气动干扰的直升机操纵响应计算方法和耦合气动干扰的直升机稳定性计算方法。采用涡面元法和粘性涡粒子法建立的直升机全机气动干扰计算模型,并利用多极展开和涡元重构两种加速技术提高计算效率,完成直升机全机干扰流场计算,同时解决了直升机全机干扰流场计算效率和精度问题。将该全机气动干扰计算模型与直升机飞行动力学模型耦合,建立耦合计算方法,形成一套彻底摆脱气动干扰风洞试验的高效高精度飞行特性计算方法,从而达到提高直升机飞行特性计算精度,缩短气动布局设计周期的目的。本计算方法适用于直升机可用包线内的所有飞行状态。
Description
技术领域
本发明涉及直升机总体设计技术领域,具体涉及一种考虑全机气动干扰的直升机飞行特性计算方法,可以应用于直升机概念和初步设计阶段的直升机气动布局设计过程中飞行性能和飞行品质计算分析。
背景技术
直升机全机气动布局设计是总体设计中的主要内容之一,它包括直升机各部件的气动外形、参数和相对位置的确定。直升机是否具有良好的飞行特性(包括配平特性、稳定性、操纵性及飞行性能等)在很大程度上决定于气动布局,因此气动布局设计在直升机总体设计中处于核心地位,而气动布局设计过程中飞行特性的准确计算是关键。由于直升机气动布局紧凑,机身、尾翼等处于旋翼尾流当中,气动干扰问题严重,对直升机的飞行特性具有重要影响,因此为获得最佳气动布局,在飞行特性计算中必须考虑各部件之间的气动干扰。
现有考虑气动干扰的直升机飞行特性计算方法主要分为两类,第一类为目前主要采用的方法。第一类方法先通过组合部件风洞试验确定气动干扰因子,然后在计算机身、尾翼等部件气动力时利用试验获得的干扰因子计入干扰速度的影响。该类方法的主要缺点在于:
(1)依赖于风洞试验,效率低、经济性差。直升机气动干扰与许多飞行状态参数密切相关,必须进行大量不同状态的试验;其次,在迭代设计过程中直升机气动布局可能不断变化,需进行各种变布局的风洞试验,而风洞试验本身周期长,费用高,因此这类计算方法效率低、经济性差。
(2)依赖于经验公式,难以包括所有飞行状态,人为因素较大。直升机气动干扰与许多飞行状态参数密切相关,如飞行重量、高度、空速、姿态、旋翼桨盘角度等等,而风洞试验状态并不能覆盖所有飞行状态,因此通常的方法是利用试验数据根据人为经验建立气动干扰因子与飞行状态参数的函数关系,进而再应用于飞行特性计算中。
第二类方法是采用旋翼固定或自由尾迹模型计算旋翼尾迹对各部件的气动干扰因子,进而在机身、尾翼等部件气动力计算时考虑旋翼的干扰作用。这类方法的主要缺点是气动干扰计算精度差,主要原因在于:
(1)旋翼尾迹模型基于势流理论,忽略粘性耗散,人为设定参数对结果影响大,如涡核半径等。
(2)只考虑旋翼尾流对其它部件的单向干扰作用,未考虑全机气动干扰,即所有气动部件及尾流的相互干扰作用。
另外,现有的直升机计算流体力学(CFD)数值计算方法尽管可较准确的计算各部件之间的气动干扰,但是由于直升机全机气动外形复杂、运动部件较多、涡流特征显著,CFD计算方法依赖于大规模、精细的、分块的复杂网格系统和各种尾涡捕捉方法,计算量巨大,即使采用大规模并行计算方法,其计算周期仍以周甚至月为度量单位,难以满足日常工程设计中飞行特性计算的需求。
发明内容
本发明的目的是为了解决直升机气动布局工程设计中现有的飞行特性计算方法的不足,提出了一种不依赖于气动干扰风洞试验数据,且计算效率和精度较高的直升机飞行特性计算方法。
本发明的方法采用基于涡面元法和粘性涡粒子法的直升机全机气动干扰计算模型(为简化描述,本发明后续部分将该模型简称为直升机全机气动干扰计算模型),进行直升机全机干扰流场计算。直升机全机气动干扰计算模型中包括了多极展开和涡元重构两种加速技术提高计算效率,同时解决了直升机全机干扰流场计算效率和精度问题。将直升机全机气动干扰计算模型与直升机飞行动力学模型耦合,建立耦合计算方法,形成了一套彻底摆脱气动干扰风洞试验的高效高精度飞行特性计算方法,从而达到提高直升机飞行特性计算精度,缩短气动布局设计周期的目的。
本发明的一种考虑全机气动干扰的直升机飞行特性计算方法,通过如下步骤实现:
11)数据准备:对直升机各孤立部件三维几何模型物面进行网格划分(物面网格为三角形或四边形平面网格),导出各孤立部件网格单元节点坐标信息数据表,制作直升机设计参数以及孤立部件气动特性数据表,设定直升机稳定飞行状态参数;
12)数据输入:将直升机设计参数数据表和孤立部件气动特性数据表输入直升机飞行动力学模型,将各孤立部件网格单元节点坐标信息数据表输入直升机全机气动干扰计算模型;
13)耦合气动干扰的直升机配平计算:首先进行直升机飞行动力学模型配平计算输出配平值,将配平值传递给直升机全机气动干扰计算模型,然后进行直升机全机气动干扰计算模型配平计算,并将直升机全机气动干扰计算模型配平计算输出的旋翼和尾桨桨盘诱导速度分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型,重复上述计算过程,直到直升机飞行动力学模型配平计算输出的配平值收敛,得到直升机稳定飞行状态配平值;
14)耦合气动干扰的直升机操纵响应计算:根据步骤13)输出的配平值对直升机飞行动力学模型和直升机全机气动干扰计算模型进行初始化,设定操纵响应计算时间、时间步长以及直升机操纵输入变化量时间历程,首先进行1个时间步的直升机飞行动力学模型时间积分,将时间积分得到的直升机运动状态量信息传递给直升机全机气动干扰计算模型,然后进行直升机全机气动干扰计算模型时间积分,将全机气动干扰计算模型时间积分后输出的旋翼和尾桨桨盘诱导速度分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型,重复上述操纵响应计算过程,直到达到设定的计算时间,得到直升机的操纵响应;
15)耦合气动干扰的直升机稳定性计算:根据步骤13)输出的配平值对直升机飞行动力学模型和直升机全机气动干扰计算模型进行初始化,设定直升机运动状态量和操纵输入量的小扰动量,采用数值差分求导的方法,对直升机飞行动力学模型进行线化,获得直升机气动导数矩阵和操纵导数矩阵,对直升机气动导数矩阵进行特征值和特征向量计算,得到直升机各运动模态的稳定根。
作为上述技术方案的进一步改进,所述步骤13)耦合气动干扰的直升机配平计算具体过程为:
43)配平初值计算:将直升机稳定飞行状态参数输入直升机飞行动力学模型,进行稳定飞行配平计算,输出配平初值;
44)全机气动干扰计算模型配平计算:将步骤43)输出的配平初值或步骤45)输出的配平值输入直升机全机气动干扰计算模型,进行全机气动干扰计算模型配平计算,输出旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量;
45)飞行动力学模型配平计算:将步骤44)输出旋翼和尾桨桨盘诱导入流分布和各部件受到其它部件干扰速度量输入直升机飞行动力学模型,进行稳定飞行配平计算,输出配平值;
46)耦合配平计算收敛判断:若步骤45)仅执行过一次,计算步骤43)输出的配平初值与步骤45)输出的配平值之间的相对残差值,若步骤45)执行次数大于1,计算前后两次步骤45)输出配平值之间的相对残差值,当相对残差值大于某小量(通常设定为10-3量级),配平计算未收敛,重复所述步骤44)和45),当相对残差值小于等于某小量,配平计算已收敛,执行步骤47);
47)耦合配平计算结果输出:输出最后一次执行步骤45)配平值和最后一次执行步骤44)过程中各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息以及旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量,耦合气动干扰的直升机配平计算结束。
作为上述技术方案的进一步改进,所述步骤14)耦合气动干扰的直升机操纵响应计算具体过程为:
52)耦合配平计算:执行耦合气动干扰的直升机配平计算方法所有步骤,即步骤41)~步骤46),直到耦合气动干扰的直升机配平计算结束;
53)初始化一:将步骤52)输出的配平值转化为直升机飞行动力学模型中各状态量和操纵输入量的初始值,将步骤52)输出的旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量输入直升机飞行动力学模型,设定当前时刻t=0,完成直升机飞行动力学模型初始化;
54)初始化二:将步骤52)输出的配平值、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息输入直升机全机气动干扰计算模型,完成直升机全机气动干扰计算模型初始化;
55)直升机飞行动力学模型时间积分:根据直升机操纵输入变化量时间历程数据表,读入t时刻直升机操纵输入变化量,采用DE/STEP常微分方程积分算法,对直升机飞行动力学模型进行1个时间步长积分,得到t+△t时刻直升机运动状态量信息,并将此运动状态量信息传递给直升机全机气动干扰计算模型;
56)直升机全机气动干扰计算模型时间积分:对直升机全机气动干扰计算模型进行1个时间步长积分,得到t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度,并将t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型;
57)时间积分判断和输出:令t=t+△t,输出步骤55)得到的直升机运动状态量信息,若t小于T,重复步骤55)和56),否者,直升机操纵响应计算结束。
作为上述技术方案的进一步改进,所述步骤15)耦合气动干扰的直升机稳定性计算具体过程为:
62)耦合配平计算:执行耦合气动干扰的直升机配平计算方法所有步骤,即步骤41)~步骤46),直到耦合气动干扰的直升机配平计算结束;
63)初始化一:将步骤62)输出的配平值转化为直升机飞行动力学模型中各状态量和操纵输入量的初始值,将步骤62)输出的旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量输入直升机飞行动力学模型,设定当前时刻t=0,完成直升机飞行动力学模型初始化;
64)初始化二:将步骤62)输出的配平值、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息输入直升机全机气动干扰计算模型,完成直升机全机气动干扰计算模型初始化;
65)t时刻气动导数和操纵导数矩阵计算:采用数值求导方法利用直升机飞行动力学模型,计算t时刻气动导数和操纵导数矩阵,具体步骤为,先计算t时刻直升机运动状态量时间变化率,在t时刻直升机运动状态量和操纵输入量的基础上,根据设定的直升机运动状态量和操纵输入量的小扰动量,分别对每一个直升机运动状态量和操纵输入量进行扰动,计算t时刻直升机运动状态量扰动后的时间变化率,将扰动后和扰动前的每个一个运动状态量时间变化率与每一个直升机运动状态量和操纵输入量的小扰动量进行数值求导计算,获得t时刻直升机的气动导数矩阵和操纵导数矩阵,并记录保存;
66)直升机飞行动力学模型时间积分:保持步骤62)输出的配平值中的操纵输入量,采用DE/STEP常微分方程积分算法,对直升机飞行动力学模型进行1个时间步长积分,得到t+△t时刻直升机运动状态量信息,并将此运动信息传递给直升机全机气动干扰计算模型;
67)直升机全机气动干扰计算模型时间积分:对直升机全机气动干扰计算模型进行1个时间步长积分,得到t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度,并将t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型;
68)气动导数和操纵导数矩阵计算判断:令t=t+△t,若t小于旋翼旋转周期,重复步骤65)和67),若t等于旋翼旋转周期,转入步骤69);
69)稳定性计算:将执行N次步骤65)保存的N个气动导数矩阵和操纵导数矩阵进行代数平均,获得设定稳定飞行状态下直升机的气动和操纵导数矩阵,对气动导数矩阵进行特征值和特征向量计算,获得直升机各运动模态的稳定根,耦合气动干扰的直升机稳定性计算结束。
作为上述技术方案的进一步改进,所述孤立部件包括直升机旋翼、机身、平尾、垂尾以及尾桨。
作为上述技术方案的另一种改进,所述稳定飞行状态参数包括直升机的飞行高度、飞行速度、爬升角和偏航角速率。
作为上述技术方案的另一种改进,所述配平初值和配平值包括直升机体轴系下机体的3个线速度、3个姿态角、3个旋翼操纵输入量、1个尾桨操纵输入量、旋翼桨盘锥度角、旋翼桨盘后倒角和旋翼桨盘侧倒角,若考虑了旋翼桨叶摆振运动,则还包括旋翼集合型、前进型和后退型摆振角。
作为上述技术方案的另一种改进,所述直升机运动状态量包括直升机体轴系下机体的3个线速度、3个角速度、3个姿态角、旋翼参考桨叶方位角、旋翼各片桨叶挥舞角、旋翼各片桨叶挥舞角速度、旋翼各片桨叶变距角、旋翼各片桨叶变距角速率,若考虑了旋翼桨叶摆振运动,则还包括旋翼各片桨叶摆振角、旋翼各片桨叶摆振角速度。
作为上述技术方案的另一种改进,所述各部件受到其它部件干扰速度量包括机身气动中心受到旋翼、平尾、垂尾和尾桨涡面元及其尾迹诱导的气流速度量,平尾气动中心受到旋翼、机身、垂尾和尾桨涡面元及其尾迹诱导的气流速度量,垂尾气动中心受到旋翼、机身、平尾和尾桨涡面元及其尾迹诱导气流速度量。
作为上述技术方案的另一种改进,所述旋翼和尾桨桨盘诱导入流分布包括旋翼和尾桨桨盘平面受到旋翼、机身、平尾、垂尾和尾桨涡面元及其尾迹诱导的气流速度分布。
与现有方法相比,本发明公开的方法的效果和优点如下:
(1)彻底摆脱直升机飞行特性计算依赖气动干扰试验数据和气动干扰经验修正的现状。
(2)计算精度高:在低速和机动飞行状态,直升机配平值、异轴响应计算精度比已有理论计算方法提高10%以上。
(3)计算效率高:采用普通个人电脑,稳定平飞状态配平计算时间约1-3小时。
(4)本计算方法适用于直升机可用包线内的所有飞行状态。
附图说明
图1是直升机全机气动干扰计算模型配平计算步骤流程图;
图2是直升机全机气动干扰计算模型时间积分步骤流程图;
图3是耦合气动干扰的直升机配平计算方法步骤流程图;
图4是耦合气动干扰的直升机操纵响应计算方法步骤流程图;
图5是耦合气动干扰的直升机稳定性计算方法步骤流程图;
图6是某直升机各部件物面网格图;
图7.1‐7.4是某直升机配平计算结果图,其中,图7.1是旋翼总距计算结果图,图7.2是旋翼需用功率计算结果图,图7.3是纵向周期变距杆量计算结果图,图7.4是俯仰姿态角计算结果图;
图8.1‐8.4是某直升机操纵输入变化量时间历程及操纵响应计算结果图,其中,图8.1是操纵输入计算结果图,图8.2是滚转角速度响应计算结果图,图8.3是俯仰角速度响应计算结果图,图8.4是偏航角速度响应计算结果图。
具体实施方式
本发明的直升机飞行动力学模型,是在已有通用的常规直升机飞行动力学模型基
础上改造而来的。通用常规直升机飞行动力学模型一般可表示为:
其中,y=(yf,ym,ytr,ymi,ytri,ye)为直升机状态量向量,为由直升机旋翼横向周期变距角B1、旋翼纵向周期变距角A1、旋翼总距角θ0、尾桨总距角θtr以及它们的时间变化率所组成的操纵输入向量。直升机状态量向量由机体运动状态量向量yf、旋翼运动状态量向量ym、尾桨运动状态量向量ytr、旋翼桨盘入流状态量向量ymi、尾桨桨盘入流状态量向量ytri以及发动机状态量向量ye组成。其中yf中包括机体3个线速度分量、3个角速度分量、3个姿态角;ym中一般需包括旋翼参考桨叶方位角、旋翼转速、旋翼各片桨叶的挥舞角和角速度、旋翼各片桨叶的摆振角和角速度(若考虑摆振运动);ytr中一般需包括尾桨参考桨叶方位角、尾桨各片桨叶的挥舞角和角速度(若考虑挥舞运动);若采用动态入流模型,ymi和ytri应分别包括旋翼和尾桨桨盘各阶入流系数;ye与所选择的发动机模型相关。
本发明对通用常规直升机飞行动力学模型进行改造,以满足直升机全机气动干扰计算模型与直升机飞行动力学模型耦合求解的需要。改造后的直升机飞行动力学模型表示为:
其中t为时间,[υmi]和[υtri]分别表示t时刻旋翼和尾桨桨盘平面诱导入流速度分布,[υia]表示t时刻机身、平尾、垂尾等部件气动中心分别受到其它部件引起的干扰气流速度,而改造后的直升机飞行动力学模型中的直升机状态量向量y=(yf,ym,ytr,ye)中不再包括旋翼和尾桨桨盘入流状态量。
本发明中所述直升机飞行动力学模型均指改造后的直升机飞行动力学模型,所述直升机飞行动力学模型配平计算和时间积分计算,均指通用常规直升机飞行动力学模型已有的配平和时间积分算法。
本发明的直升机全机气动干扰计算模型基于涡面元法和粘性涡粒子法而建立,它由直升机各部件物面网格(各部件主要包括旋翼、机身、平尾、垂尾、尾桨等气动面)、各部件物面涡面元、各部件尾迹涡粒子构成。
如图1所示,直升机全机气动干扰计算模型配平计算包括以下步骤:
21)数据准备:对直升机各孤立部件三维几何模型物面进行网格划分(物面网格为三角形或四边形平面网格),导出各孤立部件网格单元节点坐标信息数据表,设定直升机稳定飞行配平值;
22)数据输入:将各孤立部件网格单元节点坐标信息数据表输入直升机全机气动干扰计算模型;
23)初始化:采用涡面元法中面元影响系数矩阵计算方法,计算直升机各孤立部件涡面元影响系数矩阵,设定当前时刻t=0,时间步长△t,要求△t为旋翼旋转周期的1/N,将直升机稳定飞行配平值转化旋翼旋转周期内以△t为步长N组离散的直升机运动状态量信息;
24)执行N步全机气动干扰气动模型时间积分,记录每步执行后各部件尾迹涡粒子位置和强度;
25)直升机全机气动干扰计算模型配平收敛判断:若步骤24)仅执行过1次,则执行步骤24),若步骤24)执行次数大于1,判断步骤24)前后两次执行所记录的各部件尾迹涡粒子位置和强度的相对残差值,当相对残差值大于某小量(通常设定为10-3量级),配平计算未收敛,重复所述步骤14),当相对残差值小于等于某小量,配平计算已收敛,执行步骤26);
26)输出结果:计算旋翼和尾桨桨盘诱导入流分布,计算各部件受到其它部件干扰速度量,并输出,结束。
如图2所示,所述直升机全机气动干扰计算模型时间积分由t时刻到t+△t时刻一步积分包括以下步骤:
31)输入t时刻直升机运动状态量信息;
32)计算各部件涡面元强度:根据t时刻直升机运动状态量信息和各部件尾迹对各部件涡面元控制点诱导速度,利用各部件影响系数矩阵,计算t时刻各部件涡面元强度;
33)各部件生成新的尾迹涡粒子:根据t时刻直升机运动状态量信息和t-△t时刻及t时刻各部件涡面元强度,确定各部件新生成尾迹涡粒子的位置和强度;
34)计算各部件尾迹涡粒子运动速度和梯度:基于多极展开加速技术计算t时刻各部件尾迹涡粒子运动速度和梯度;
35)各部件尾迹涡粒子位置和涡强时间积分;根据t时刻各部件尾迹涡粒子运动速度和梯度,采用数值积分算法,获得t+△t时刻各部件尾迹涡粒子的位置和强度;
36)涡粒子重构:采用涡粒子重构加速技术,对t+△t时刻各部件尾迹涡粒子进行合并和分裂;
37)计算各部件涡面元控制点诱导速度:根据t+△t时刻各部件尾迹涡粒子的位置和强度,计算t+△t时刻各部件尾迹对各部件涡面元控制点诱导速度,令t=t+△t,则t时刻到t+△t时刻直升机全机气动干扰计算模型时间积分结束。
如图3所示,耦合气动干扰的直升机配平计算方法包括以下步骤:
41)数据准备:对直升机各孤立部件三维几何模型物面进行网格划分(物面网格为三角形或四边形平面网格),导出各孤立部件网格单元节点坐标信息数据表,制作直升机设计参数以及孤立部件气动特性数据表,设定直升机稳定飞行状态参数;
42)数据输入:将直升机设计参数数据表和孤立部件气动特性数据表输入直升机飞行动力学模型,将各孤立部件网格单元节点坐标信息数据表输入直升机全机气动干扰计算模型;
43)配平初值计算:将直升机稳定飞行状态参数输入直升机飞行动力学模型,进行稳定飞行配平计算,输出配平初值;
44)全机气动干扰计算模型配平计算:将步骤43)输出的配平初值或步骤45)输出的配平值输入直升机全机气动干扰计算模型,进行全机气动干扰计算模型配平计算,输出旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量;
45)飞行动力学模型配平计算:将步骤44)输出旋翼和尾桨桨盘诱导入流分布和各部件受到其它部件干扰速度量输入直升机飞行动力学模型,进行稳定飞行配平计算,输出配平值;
46)耦合配平计算收敛判断:若步骤45)仅执行过一次,计算步骤43)输出的配平初值与步骤45)输出的配平值之间的相对残差值,若步骤45)执行次数大于1,计算前后两次步骤45)输出配平值之间的相对残差值,当相对残差值大于某小量(通常设定为10-3量级),配平计算未收敛,重复所述步骤44)和45),当相对残差值小于等于某小量,配平计算已收敛,执行步骤47);
47)耦合配平计算结果输出:输出最后一次执行步骤45)配平值和最后一次执行步骤44)过程中各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息以及旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量,耦合气动干扰的直升机配平计算结束。
如图4所示,耦合气动干扰的直升机操纵响应计算方法包括以下步骤:
51)数据准备:制作直升机操纵输入变化量时间历程数据表,设定操纵响应计算时间T和时间步长△t;
52)耦合配平计算:执行耦合气动干扰的直升机配平计算方法所有步骤,即步骤41)~步骤46),直到耦合气动干扰的直升机配平计算结束;
53)初始化一:将步骤52)输出的配平值转化为直升机飞行动力学模型中各状态量和操纵输入量的初始值,将步骤52)输出的旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量输入直升机飞行动力学模型,设定当前时刻t=0,完成直升机飞行动力学模型初始化;
54)初始化二:将步骤52)输出的配平值、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息输入直升机全机气动干扰计算模型,完成直升机全机气动干扰计算模型初始化;
55)直升机飞行动力学模型时间积分:根据直升机操纵输入变化量时间历程数据表,读入t时刻直升机操纵输入变化量,采用DE/STEP常微分方程积分算法,对直升机飞行动力学模型进行1个时间步长积分,得到t+△t时刻直升机运动状态量信息,并将此运动状态量信息传递给直升机全机气动干扰计算模型;
56)直升机全机气动干扰计算模型时间积分:对直升机全机气动干扰计算模型进行1个时间步长积分,得到t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度,并将t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型;
57)时间积分判断和输出:令t=t+△t,输出步骤55)得到的直升机运动状态量信息,若t小于T,重复步骤55)和56),否者,直升机操纵响应计算结束。
如图5所示,耦合气动干扰的直升机稳定性计算方法包括以下步骤:
61)数据准备:设定时间步长△t,要求△t为旋翼旋转周期的1/N,设定直升机运动状态量和操纵输入量的小扰动量;
62)耦合配平计算:执行耦合气动干扰的直升机配平计算方法所有步骤,即步骤41)~步骤46),直到耦合气动干扰的直升机配平计算结束;
63)初始化一:将步骤62)输出的配平值转化为直升机飞行动力学模型中各状态量和操纵输入量的初始值,将步骤62)输出的旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量输入直升机飞行动力学模型,设定当前时刻t=0,完成直升机飞行动力学模型初始化;
64)初始化二:将步骤62)输出的配平值、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息输入直升机全机气动干扰计算模型,完成直升机全机气动干扰计算模型初始化;
65)t时刻气动导数和操纵导数矩阵计算:采用数值求导方法利用直升机飞行动力学模型,计算t时刻气动导数和操纵导数矩阵,具体步骤为,先计算t时刻直升机运动状态量时间变化率,在t时刻直升机运动状态量和操纵输入量的基础上,根据设定的直升机运动状态量和操纵输入量的小扰动量,分别对每一个直升机运动状态量和操纵输入量进行扰动,计算t时刻直升机运动状态量扰动后的时间变化率,将扰动后和扰动前的每个一个运动状态量时间变化率与每一个直升机运动状态量和操纵输入量的小扰动量进行数值求导计算,获得t时刻直升机的气动导数矩阵和操纵导数矩阵,并记录保存;
66)直升机飞行动力学模型时间积分:保持步骤62)输出的配平值中的操纵输入量,采用DE/STEP常微分方程积分算法,对直升机飞行动力学模型进行1个时间步长积分,得到t+△t时刻直升机运动状态量信息,并将此运动信息传递给直升机全机气动干扰计算模型;
67)直升机全机气动干扰计算模型时间积分:对直升机全机气动干扰计算模型进行1个时间步长积分,得到t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度,并将t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型;
68)气动导数和操纵导数矩阵计算判断:令t=t+△t,若t小于旋翼旋转周期,重复步骤65)和67),若t等于旋翼旋转周期,转入步骤69);
69)稳定性计算:将执行N次步骤65)保存的N个气动导数矩阵和操纵导数矩阵进行代数平均,获得设定稳定飞行状态下直升机的气动和操纵导数矩阵,对气动导数矩阵进行特征值和特征向量计算,获得直升机各运动模态的稳定根,耦合气动干扰的直升机稳定性计算结束。
应用例:
本实例中,将本发明提出的一种考虑全机气动干扰的直升机飞行特性计算方法应用于实际的某直升机配平、操纵响应和稳定性计算中。
71)配平计算
如图6所示,应用的模拟直升机包括机身1、旋翼2、平尾3、垂尾4和尾桨5。照具体实施方式中耦合气动干扰的直升机配平计算方法实施步骤,首先采用CATIA或Gabit等商用软件中的网格划分功能对各孤立部件三维几何模型物面进行网格划分(如图6所示),导出各孤立部件网格单元节点坐标信息数据表,制作直升机设计参数表(如表1所示)和各孤立部件气动特性数据表(数据量太大,在本应用例中未给出),设定8组直升机稳定平直飞行状态参数(如表2所示)。然后,按照耦合气动干扰的直升机配平计算方法余下各步骤对8组稳定飞行状态参数分别进行计算,得到了对应稳定飞行状态的配平值,主要配平值如图7.1至7.4所示,图中同时给出了飞行试验结果和不考虑全机气动干扰的飞行动力学模型计算结果,由图中对比可知本本发明提出的方法,大幅提高了直升机配平计算精度,尤其是在低速飞行状态。
表1
表2
72)操纵响应计算
按照具体实施方式中耦合气动干扰的直升机操纵响应计算方法实施步骤,首先设定直升机前飞右压杆操纵输入变化量时间历程(如图8.1),设定操纵响应计算时间T=6s和时间步长△t=0.01s,选择表2中的第8组稳定飞行状态参数,按照耦合气动干扰的直升机操纵响应计算方法余下各步骤进行操纵响应计算,得到直升机的操纵响应(如图8.2-图8.4),图中同时给出了飞行试验结果和不考虑全机气动干扰的飞行动力学模型计算结果,由图中对比可知本发明提出的方法有效提高了机动飞行操纵响应计算精度,尤其是异轴响应的计算精度。
73)稳定性计算
按照具体实施方式中耦合气动干扰的直升机稳定性计算方法实施步骤,首先设定直升机运动状态量和操纵输入量的小扰动量,对于速度量纲设置扰动量为6.5×10-3,对于角度和角速度量纲扰动量设置为6.5×10-2,设定时间步长△t=0.00727,则N=32,选择表2中的第8组稳定飞行状态参数按照耦合气动干扰的直升机操纵响应计算方法余下各步骤进行稳定性计算,得到直升机各运动模态的稳定根,如表3所示。
表3
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进,这些改进也应视为本发明的保护范围。
Claims (8)
1.一种考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:包括考虑直升机全机气动干扰的配平、操纵响应和稳定性计算方法,所述一种考虑全机气动干扰的直升机飞行特性计算方法,通过如下步骤实现:
11)数据准备:对直升机各孤立部件三维几何模型物面进行网格划分,导出各孤立部件网格单元节点坐标信息数据表,制作直升机设计参数以及孤立部件气动特性数据表,设定直升机稳定飞行状态参数;
12)数据输入:将直升机设计参数数据表和孤立部件气动特性数据表输入直升机飞行动力学模型,将各孤立部件网格单元节点坐标信息数据表输入直升机全机气动干扰计算模型;
13)耦合气动干扰的直升机配平计算:首先进行直升机飞行动力学模型配平计算输出配平值,将配平值传递给直升机全机气动干扰计算模型,然后进行直升机全机气动干扰计算模型配平计算,并将直升机全机气动干扰计算模型配平计算输出的旋翼和尾桨桨盘诱导速度分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型,重复上述计算过程,直到直升机飞行动力学模型配平计算输出的配平值收敛,得到直升机稳定飞行状态配平值;
14)耦合气动干扰的直升机操纵响应计算:根据步骤13)输出的配平值对直升机飞行动力学模型和直升机全机气动干扰计算模型进行初始化,设定操纵响应计算时间、时间步长以及直升机操纵输入变化量时间历程,首先进行1个时间步的直升机飞行动力学模型时间积分,将时间积分得到的直升机运动状态量信息传递给直升机全机气动干扰计算模型,然后进行直升机全机气动干扰计算模型时间积分,将全机气动干扰计算模型时间积分后输出的旋翼和尾桨桨盘诱导速度分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型,重复上述操纵响应计算过程,直到达到设定的计算时间,得到直升机的操纵响应;
15)耦合气动干扰的直升机稳定性计算:根据步骤13)输出的配平值对直升机飞行动力学模型和直升机全机气动干扰计算模型进行初始化,设定直升机运动状态量和操纵输入量的小扰动量,采用数值差分求导的方法,对直升机飞行动力学模型进行线化,获得直升机气动导数矩阵和操纵导数矩阵,对直升机气动导数矩阵进行特征值和特征向量计算,得到直升机各运动模态的稳定根;所述步骤13)耦合气动干扰的直升机配平计算具体过程为:43)配平初值计算:将直升机稳定飞行状态参数输入直升机飞行动力学模型,进行稳定飞行配平计算,输出配平初值;
44)全机气动干扰计算模型配平计算:将步骤43)输出的配平初值输入直升机全机气动干扰计算模型,进行全机气动干扰计算模型配平计算,输出旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量;如果是从步骤46返回步骤44,则将步骤45)输出的配平值输入直升机全机气动干扰计算模型,进行全机气动干扰计算模型配平计算,输出旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量;
45)飞行动力学模型配平计算:将步骤44)输出旋翼和尾桨桨盘诱导入流分布和各部件受到其它部件干扰速度量输入直升机飞行动力学模型,进行稳定飞行配平计算,输出配平值;
46)耦合配平计算收敛判断:若步骤45)仅执行过一次,计算步骤43)输出的配平初值与步骤45)输出的配平值之间的相对残差值,若步骤45)执行次数大于1,计算前后两次步骤45)输出配平值之间的相对残差值,当相对残差值大于某小量,配平计算未收敛,重复所述步骤44)和45),当相对残差值小于等于某小量,配平计算已收敛,执行步骤47);
47)耦合配平计算结果输出:输出最后一次执行步骤45)配平值和最后一次执行步骤44)过程中各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息以及旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量,耦合气动干扰的直升机配平计算结束;
所述步骤14)耦合气动干扰的直升机操纵响应计算具体过程为:
52)耦合配平计算:执行耦合气动干扰的直升机配平计算方法所有步骤,即步骤43)~步骤46),直到耦合气动干扰的直升机配平计算结束;
53)初始化一:将步骤52)输出的配平值转化为直升机飞行动力学模型中各状态量和操纵输入量的初始值,将步骤52)输出的旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量输入直升机飞行动力学模型,设定当前时刻t=0,完成直升机飞行动力学模型初始化;
54)初始化二:将步骤52)输出的配平值、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息输入直升机全机气动干扰计算模型,完成直升机全机气动干扰计算模型初始化;
55)直升机飞行动力学模型时间积分:根据直升机操纵输入变化量时间历程数据表,读入t时刻直升机操纵输入变化量,采用DE/STEP常微分方程积分算法,对直升机飞行动力学模型进行1个时间步长积分,得到t+△t时刻直升机运动状态量信息,并将此运动状态量信息传递给直升机全机气动干扰计算模型;
56)直升机全机气动干扰计算模型时间积分:对直升机全机气动干扰计算模型进行1个时间步长积分,得到t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度,并将t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型;
57)时间积分判断和输出:令t=t+△t,输出步骤55)得到的直升机运动状态量信息,若t小于T,重复步骤55)和56),否者,直升机操纵响应计算结束。
2.根据权利要求1所述的考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:所述步骤15)耦合气动干扰的直升机稳定性计算具体过程为:
62)耦合配平计算:执行耦合气动干扰的直升机配平计算方法所有步骤,即步骤43)~步骤46),直到耦合气动干扰的直升机配平计算结束;
63)初始化一:将步骤62)输出的配平值转化为直升机飞行动力学模型中各状态量和操纵输入量的初始值,将步骤62)输出的旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量输入直升机飞行动力学模型,设定当前时刻t=0,完成直升机飞行动力学模型初始化;
64)初始化二:将步骤62)输出的配平值、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度数据信息输入直升机全机气动干扰计算模型,完成直升机全机气动干扰计算模型初始化;
65)t时刻气动导数和操纵导数矩阵计算:采用数值求导方法利用直升机飞行动力学模型,计算t时刻气动导数和操纵导数矩阵,具体步骤为,先计算t时刻直升机运动状态量时间变化率,在t时刻直升机运动状态量和操纵输入量的基础上,根据设定的直升机运动状态量和操纵输入量的小扰动量,分别对每一个直升机运动状态量和操纵输入量进行扰动,计算t时刻直升机运动状态量扰动后的时间变化率,将扰动后和扰动前的每个一个运动状态量时间变化率与每一个直升机运动状态量和操纵输入量的小扰动量进行数值求导计算,获得t时刻直升机的气动导数矩阵和操纵导数矩阵,并记录保存;
66)直升机飞行动力学模型时间积分:保持步骤62)输出的配平值中的操纵输入量,采用DE/STEP常微分方程积分算法,对直升机飞行动力学模型进行1个时间步长积分,得到t+△t时刻直升机运动状态量信息,并将此运动信息传递给直升机全机气动干扰计算模型;
67)直升机全机气动干扰计算模型时间积分:对直升机全机气动干扰计算模型进行1个时间步长积分,得到t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量、各部件涡面元强度、各部件尾迹结构中涡粒子位置和强度,并将t+△t时刻旋翼和尾桨桨盘诱导入流分布、各部件受到其它部件干扰速度量传递给直升机飞行动力学模型;
68)气动导数和操纵导数矩阵计算判断:令t=t+△t,若t小于旋翼旋转周期,重复步骤65)和67),若t等于旋翼旋转周期,转入步骤69);
69)稳定性计算:将执行N次步骤65)保存的N个气动导数矩阵和操纵导数矩阵进行代数平均,获得设定稳定飞行状态下直升机的气动和操纵导数矩阵,对气动导数矩阵进行特征值和特征向量计算,获得直升机各运动模态的稳定根,耦合气动干扰的直升机稳定性计算结束。
3.根据权利要求1所述的考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:所述孤立部件包括直升机旋翼、机身、平尾、垂尾以及尾桨;、所述稳定飞行状态参数包括直升机的飞行高度、飞行速度、爬升角和偏航角速率。
4.根据权利要求1所述的考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:所述配平初值和配平值包括直升机体轴系下机体的3个线速度、3个姿态角、3个旋翼操纵输入量、1个尾桨操纵输入量、旋翼桨盘锥度角、旋翼桨盘后倒角和旋翼桨盘侧倒角,若考虑了旋翼桨叶摆振运动,则还包括旋翼集合型、前进型和后退型摆振角。
5.根据权利要求1所述的考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:所述直升机运动状态量包括直升机体轴系下机体的3个线速度、3个角速度、3个姿态角、旋翼参考桨叶方位角、旋翼各片桨叶挥舞角、旋翼各片桨叶挥舞角速度、旋翼各片桨叶变距角、旋翼各片桨叶变距角速率,若考虑了旋翼桨叶摆振运动,则还包括旋翼各片桨叶摆振角、旋翼各片桨叶摆振角速度。
6.根据权利要求1所述的考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:所述各部件受到其它部件干扰速度量包括机身气动中心受到旋翼、平尾、垂尾和尾桨涡面元及其尾迹诱导的气流速度量,平尾气动中心受到旋翼、机身、垂尾和尾桨涡面元及其尾迹诱导的气流速度量,垂尾气动中心受到旋翼、机身、平尾和尾桨涡面元及其尾迹诱导气流速度量。
7.根据权利要求1所述的考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:所述旋翼和尾桨桨盘诱导入流分布包括旋翼和尾桨桨盘平面受到旋翼、机身、平尾、垂尾和尾桨涡面元及其尾迹诱导的气流速度分布。
8.根据权利要求1所述的考虑全机气动干扰的直升机飞行特性计算方法,其特征在于:步骤13)、14)、15)中,直升机飞行动力学模型和直升机全机气动干扰计算模型实现了耦合计算,即直升机飞行动力学模型的配平或时间积分计算输出结果将作为直升机全机气动干扰计算模型配平或时间积分计算的输入量,而反过来直升机全机气动干扰计算模型的配平或时间积分计算输出结果将作为直升机飞行动力学模型配平或时间积分计算的输入量,两者实现耦合配平计算或耦合时间积分计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510716964.0A CN105468814B (zh) | 2015-10-29 | 2015-10-29 | 一种考虑全机气动干扰的直升机飞行特性计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510716964.0A CN105468814B (zh) | 2015-10-29 | 2015-10-29 | 一种考虑全机气动干扰的直升机飞行特性计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105468814A CN105468814A (zh) | 2016-04-06 |
CN105468814B true CN105468814B (zh) | 2018-11-09 |
Family
ID=55606510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510716964.0A Expired - Fee Related CN105468814B (zh) | 2015-10-29 | 2015-10-29 | 一种考虑全机气动干扰的直升机飞行特性计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105468814B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106407731B (zh) * | 2016-11-30 | 2019-05-07 | 中国直升机设计研究所 | 一种气动干扰流场数值计算方法 |
CN108228921B (zh) * | 2016-12-13 | 2021-07-13 | 北京空天技术研究所 | 风洞试验装置及其设计方法 |
CN106650077B (zh) * | 2016-12-15 | 2021-07-16 | 中国航空工业集团公司西安飞机设计研究所 | 一种弹性飞机尾涡遭遇动响应分析方法 |
CN107909883A (zh) * | 2017-11-03 | 2018-04-13 | 中国民用航空飞行学院 | 罗宾逊r44直升机多任务飞行训练器 |
CN108090302B (zh) * | 2018-01-04 | 2021-11-30 | 北京中航智科技有限公司 | 一种直升机飞行力学模拟方法及系统 |
CN108873862B (zh) * | 2018-06-15 | 2021-06-29 | 上海航天控制技术研究所 | 一种针对飞行器的控制系统稳定性的综合评估方法 |
CN112182752B (zh) * | 2020-09-25 | 2022-11-18 | 中国直升机设计研究所 | 一种直升机飞行姿态预测方法 |
CN112668092B (zh) * | 2020-12-07 | 2023-10-31 | 沈阳飞机设计研究所扬州协同创新研究院有限公司 | 一种耦合气动干扰的飞行器混合配平分析方法 |
CN113884268B (zh) * | 2021-12-08 | 2022-02-22 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种全动平尾直升机机身纵向气动特性试验及分析方法 |
CN115655642B (zh) * | 2022-12-09 | 2023-04-07 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种有效评估直升机旋翼气动性能的试验方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6092990A (en) * | 1997-06-05 | 2000-07-25 | Mcdonnell Douglas Helicopter Company | Oscillating air jets for helicopter rotor aerodynamic control and BVI noise reduction |
CN104657529A (zh) * | 2013-11-25 | 2015-05-27 | 中国直升机设计研究所 | 一种从飞行性能角度进行直升机气动布局参数设计的方法 |
CN104881510A (zh) * | 2015-02-13 | 2015-09-02 | 南京航空航天大学 | 一种直升机旋翼/尾桨气动干扰数值仿真方法 |
CN104899365A (zh) * | 2015-05-27 | 2015-09-09 | 南京航空航天大学 | 一种可减小气动干扰不利影响的直升机气动布局优化方法 |
-
2015
- 2015-10-29 CN CN201510716964.0A patent/CN105468814B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6092990A (en) * | 1997-06-05 | 2000-07-25 | Mcdonnell Douglas Helicopter Company | Oscillating air jets for helicopter rotor aerodynamic control and BVI noise reduction |
CN104657529A (zh) * | 2013-11-25 | 2015-05-27 | 中国直升机设计研究所 | 一种从飞行性能角度进行直升机气动布局参数设计的方法 |
CN104881510A (zh) * | 2015-02-13 | 2015-09-02 | 南京航空航天大学 | 一种直升机旋翼/尾桨气动干扰数值仿真方法 |
CN104899365A (zh) * | 2015-05-27 | 2015-09-09 | 南京航空航天大学 | 一种可减小气动干扰不利影响的直升机气动布局优化方法 |
Non-Patent Citations (3)
Title |
---|
Dynamic Inflow for Practical Applications;Peter D et al.;《Journal of the American Helicopter Society》;19881231;第64-68页 * |
一种适用于直升机概念设计的飞行动力学分析方法;桂敬冉;《南京航空航天大学学报》;20150430;第47卷(第2期);第259-265页 * |
适用于直升机飞行力学分析的三维空间大气紊流模型;吉洪蕾;《航空学报》;20140725;第35卷(第7期);第1825-1835页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105468814A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105468814B (zh) | 一种考虑全机气动干扰的直升机飞行特性计算方法 | |
CN109614633B (zh) | 一种复合式旋翼飞行器非线性建模及线性化配平方法 | |
Ventura Diaz et al. | High-fidelity computational aerodynamics of multi-rotor unmanned aerial vehicles | |
CN113567083B (zh) | 一种全动平尾直升机多部件气动干扰特性试验模拟方法 | |
CN104881510B (zh) | 一种直升机旋翼/尾桨气动干扰数值仿真方法 | |
CN108090302B (zh) | 一种直升机飞行力学模拟方法及系统 | |
CN105676672B (zh) | 扑翼飞行器复合飞行策略仿真建模方法 | |
Öhrle et al. | Compound helicopter X 3 in high-speed flight: correlation of simulation and flight test | |
CN109969426A (zh) | 一种针对复合推力构型直升机的升力分配方法及系统 | |
CN112699622A (zh) | 一种固定翼无人机的紧密编队气动耦合效应建模方法 | |
CN112733276A (zh) | 倾转旋翼无人机的倾转机构动力学优化方法及系统 | |
CN115758940A (zh) | 一种适用于倾转旋翼飞行器着舰飞行特性分析方法 | |
Zhao et al. | Physics-based modeling of viscous ground effect for rotorcraft applications | |
Singh et al. | Dynamic stall modeling using viscous vortex particle method for coaxial rotors | |
Cook et al. | Modelling the flight dynamics of the hang glider | |
Mathur et al. | Wind Tunnel Testing and Aerodynamic Characterization of a QuadPlane Uncrewed Aircraft System | |
CN111881632A (zh) | 一种直升机风限图确定方法及系统 | |
Zhao et al. | Geometry shape selection of NACA airfoils for Mars rotorcraft | |
CN212220549U (zh) | 采用机载流量传感增强风中旋翼俯仰控制的飞行器 | |
CN112329155A (zh) | 交叉双旋翼无人直升机的仿真方法和装置 | |
Pérez et al. | Free-Vortex Wake and CFD Simulation of a Small Rotor for a Quadcopter at Hover | |
O'Donnell et al. | Aerodynamic parameter estimation from wind tunnel testing of a small UAS | |
PRUDDEN | Rotor aerodynamic interaction effects for multirotor unmanned aircraft systems in forward flight | |
Mathew et al. | Albatross and Falcon inspired Bionic UAV: An Aerodynamic Analysis | |
Perez et al. | Computational study of a small rotor at hover using CFD and UVLM |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181109 Termination date: 20201029 |
|
CF01 | Termination of patent right due to non-payment of annual fee |