CN109492237A - 一种气动系数获得方法 - Google Patents
一种气动系数获得方法 Download PDFInfo
- Publication number
- CN109492237A CN109492237A CN201710817189.7A CN201710817189A CN109492237A CN 109492237 A CN109492237 A CN 109492237A CN 201710817189 A CN201710817189 A CN 201710817189A CN 109492237 A CN109492237 A CN 109492237A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- aerodynamic
- angle
- aerodynamic coefficient
- flight
- 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.)
- Pending
Links
- 238000002360 preparation method Methods 0.000 title claims abstract description 11
- 230000001133 acceleration Effects 0.000 claims abstract description 14
- 238000005553 drilling Methods 0.000 claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 3
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000005096 rolling process Methods 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000004088 simulation Methods 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
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- 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)
- Aviation & Aerospace Engineering (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Mathematical Analysis (AREA)
- Transportation (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明涉及一种气动系数获得方法,包括获得飞行器气动系数;解算得到气动参数;计算飞行攻角α、飞行侧滑角β;插值得到插值直角坐标气动系数;获得极坐标体系气动系数;得到飞行器的加速度和角加速度;得到飞行器下一时刻的飞行速度和飞行姿态角;解算飞行器下一时刻的气动系数的步骤。本发明的方法,节省了CFD计算、风洞试验等耗费的研制费用和时间成本。
Description
技术领域
本发明属于航空航天领域,具体涉及一种气动系数获得方法。
背景技术
根据飞行器体轴系与气流空间几何角度关系,飞行器使用的气动系数分为两种体系:直角坐标体系和极坐标体系。直角坐标体系气动系数随攻角α和侧滑角β两维变化,数据点在空间上呈现矩阵分布,如图1所示。极坐标体系气动系数随总攻角αΣ和滚转角Φ两维变化,数据点在空间上呈现圆形分布,靠近圆心处的数据较密,远离圆心处的数据较稀疏,如图2所示。由图1和图2的数据点分布图可知,直角坐标体系气动系数和极坐标体系气动系数,除了少量的数据点重合外,绝大多数的数据点是不重合的,分布规律差异较大。
直角坐标体系气动系数和极坐标体系气动系数之间存在差异,导致飞行仿真和控制回路设计等在使用两种坐标体系的气动系数方式上存在较大差异,只能使用相同坐标体系下的气动系数,两种坐标体系下的气动系数不能通用,要改变飞行仿真和控制回路设计的坐标体系,必须通过风洞试验、CFD计算等方式获得相同坐标体系下的气动系数。
发明内容
本发明的目的是:解决直角坐标体系气动系数不能在极坐标体系下使用的问题,找出一种气动系数获得方法。
考虑到现有技术的上述问题,根据本发明公开的一个方面,本发明采用以下技术方案:
一种气动系数获得方法,包括以下步骤:
第一步,通过风洞试验获得飞行器直角坐标体系下的气动系数;
第二步,获得飞行器飞行参数,并根据飞行参数解算得到气动参数;其中,所述的飞行参数为飞行速度和飞行姿态角,所述的气动参数为总攻角α∑和滚转角Φ;
第三步,根据总攻角α∑和滚转角Φ计算直角坐标体系下的飞行攻角α、飞行侧滑角β:
tanα=tanα∑cosΦ
sinβ=sinα∑sinΦ;
第四步,根据所述的飞行攻角α、飞行侧滑角β,插值得到插值直角坐标气动系数,所述的插值直角坐标气动系数包括轴向力系数法向力系数横向力系数滚转力矩系数偏航力矩系数和俯仰力矩系数
第五步,根据所述的插值直角坐标气动系数和转换公式获得极坐标体系气动系数,所述的极坐标体系气动系数为轴向力系数法向力系数 Cξ(αΣ,Φ)、横向力系数Cη(αΣ,Φ)、滚转力矩系数偏航力矩系数mξ(αΣ,Φ)、俯仰力矩系数mη(αΣ,Φ);其中的转换公式如下:
和
第六步,使用第五步得到的极坐标体系气动系数进行动力学解算,得到飞行器的加速度和角加速度;
第七步,使用第六步得到的飞行器的加速度和角加速度,进行运动学解算,得到飞行器下一时刻的飞行速度和飞行姿态角;
第八步,根据下一时刻的飞行速度和飞行姿态角解算得到飞行器下一时刻的气动系数,并反复执行步骤一至八获得各个时刻的气动系数,供仿真使用。
有益效果:
本发明通过把直角坐标体系下的气动系数转换成适于极标体系使用的气动系数,节省了CFD计算、风洞试验等耗费的研制费用和时间成本。
附图说明
图1气动系数直角坐标形式;
图2气动系数极坐标形式。
具体实施方式
下面结合实施例对本发明作进一步地详细说明,但本发明的实施方式不限于此。
本发明的一种气动系数获得方法,包括以下步骤:
第一步,通过风洞试验获得飞行器直角坐标体系下的气动系数;
第二步,获得飞行器飞行参数,并根据飞行参数解算得到气动参数;其中,所述的飞行参数为飞行速度和飞行姿态角,所述的气动参数为总攻角α∑和滚转角Φ;
第三步,根据总攻角α∑和滚转角Φ计算直角坐标体系下的飞行攻角α、飞行侧滑角β:
tanα=tanα∑cosΦ
sinβ=sinα∑sinΦ;
第四步,根据所述的飞行攻角α、飞行侧滑角β,插值得到插值直角坐标气动系数,所述的插值直角坐标气动系数包括轴向力系数法向力系数横向力系数滚转力矩系数偏航力矩系数和俯仰力矩系数
第五步,根据所述的插值直角坐标气动系数和转换公式获得极坐标体系气动系数,所述的极坐标体系气动系数为轴向力系数法向力系数 Cξ(αΣ,Φ)、横向力系数Cη(αΣ,Φ)、滚转力矩系数偏航力矩系数mξ(αΣ,Φ)、俯仰力矩系数mη(αΣ,Φ);其中的转换公式如下:
和
第六步,使用第五步得到的极坐标体系气动系数进行动力学解算,得到飞行器的加速度和角加速度;
第七步,使用第六步得到的飞行器的加速度和角加速度,进行运动学解算,得到飞行器下一时刻的飞行速度和飞行姿态角;
第八步,根据下一时刻的飞行速度和飞行姿态角解算得到飞行器下一时刻的气动系数,并反复执行步骤一至八获得各个时刻的气动系数,供仿真使用。
Claims (4)
1.一种气动系数获得方法,包括以下步骤:
第一步,通过风洞试验获得飞行器直角坐标体系下的气动系数;
第二步,获得飞行器飞行参数,并根据飞行参数解算得到气动参数;其中,所述的飞行参数为飞行速度和飞行姿态角,所述的气动参数为总攻角α∑和滚转角Φ;
第三步,根据总攻角α∑和滚转角Φ计算直角坐标体系下的飞行攻角α、飞行侧滑角β;
第四步,根据所述的飞行攻角α、飞行侧滑角β,插值得到插值直角坐标气动系数,所述的插值直角坐标气动系数包括轴向力系数法向力系数横向力系数滚转力矩系数偏航力矩系数和俯仰力矩系数
第五步,根据所述的插值直角坐标气动系数和转换公式获得极坐标体系气动系数;
第六步,使用第五步得到的极坐标体系气动系数进行动力学解算,得到飞行器的加速度和角加速度;
第七步,使用第六步得到的飞行器的加速度和角加速度,进行运动学解算,得到飞行器下一时刻的飞行速度和飞行姿态角;
第八步,根据下一时刻的飞行速度和飞行姿态角解算得到飞行器下一时刻的气动系数,并反复执行步骤一至八获得各个时刻的气动系数,供仿真使用。
2.如权利要求1所述的一种气动系数获得方法,其特征在于,第三步中采用如下方法计算直角坐标体系下的飞行攻角α、飞行侧滑角β:
tanα=tanα∑cosΦ
sinβ=sinα∑sinΦ。
3.如权利要求2所述的一种气动系数获得方法,其特征在于,第五步中所述的极坐标体系气动系数为轴向力系数法向力系数Cξ(αΣ,Φ)、横向力系数Cη(αΣ,Φ)、滚转力矩系数偏航力矩系数mξ(αΣ,Φ)、俯仰力矩系数mη(αΣ,Φ)。
4.如权利要求3所述的一种气动系数获得方法,其特征在于,第五步中的转换公式如下:
和
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710817189.7A CN109492237A (zh) | 2017-09-12 | 2017-09-12 | 一种气动系数获得方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710817189.7A CN109492237A (zh) | 2017-09-12 | 2017-09-12 | 一种气动系数获得方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109492237A true CN109492237A (zh) | 2019-03-19 |
Family
ID=65687729
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710817189.7A Pending CN109492237A (zh) | 2017-09-12 | 2017-09-12 | 一种气动系数获得方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109492237A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110717239A (zh) * | 2019-07-30 | 2020-01-21 | 蓝箭航天空间科技股份有限公司 | 一种箭(弹)体分布气动特性计算方法 |
CN111985169A (zh) * | 2020-08-26 | 2020-11-24 | 北京中科宇航技术有限公司 | 一种近轴对称飞行器气动特性建模方法 |
CN112595487A (zh) * | 2020-12-29 | 2021-04-02 | 中国航天空气动力技术研究院 | 一种旋转飞行器动态测力风洞试验数据的处理系统及方法 |
CN115127768A (zh) * | 2022-09-01 | 2022-09-30 | 中国航空工业集团公司哈尔滨空气动力研究所 | 基于舵面可控模型的飞行器纵向运动模拟风洞试验方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU744634A1 (ru) * | 1978-01-05 | 1980-06-30 | Ленинградский Ордена Ленина Электротехнический Институт Им. В.И. Ульянова (Ленина) | Устройство дл преобразовани пр моугольных координат в пол рные |
CN106228014A (zh) * | 2016-07-27 | 2016-12-14 | 江西洪都航空工业集团有限责任公司 | 一种导弹气动系数的获取方法 |
-
2017
- 2017-09-12 CN CN201710817189.7A patent/CN109492237A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU744634A1 (ru) * | 1978-01-05 | 1980-06-30 | Ленинградский Ордена Ленина Электротехнический Институт Им. В.И. Ульянова (Ленина) | Устройство дл преобразовани пр моугольных координат в пол рные |
CN106228014A (zh) * | 2016-07-27 | 2016-12-14 | 江西洪都航空工业集团有限责任公司 | 一种导弹气动系数的获取方法 |
Non-Patent Citations (4)
Title |
---|
XILIN YANG等: ""A 3D Collision Avoidance Strategy for UAVs in a Non-Cooperative Environment"", 《J INTELL ROBOT SYST》 * |
宋振铎等: "击顶折转轻型反坦克导弹弹道建模及仿真", 《兵工学报》 * |
杨黔龙等: "引入侧向加速度反馈的侧滑转弯控制方法", 《探测与控制学报》 * |
郝辉等: ""导弹常用空间直角坐标系间转换方法"", 《四川兵工学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110717239A (zh) * | 2019-07-30 | 2020-01-21 | 蓝箭航天空间科技股份有限公司 | 一种箭(弹)体分布气动特性计算方法 |
CN110717239B (zh) * | 2019-07-30 | 2020-07-03 | 蓝箭航天空间科技股份有限公司 | 一种箭体或弹体分布气动特性计算方法 |
CN111985169A (zh) * | 2020-08-26 | 2020-11-24 | 北京中科宇航技术有限公司 | 一种近轴对称飞行器气动特性建模方法 |
CN111985169B (zh) * | 2020-08-26 | 2024-04-02 | 北京中科宇航技术有限公司 | 一种近轴对称飞行器气动特性建模方法 |
CN112595487A (zh) * | 2020-12-29 | 2021-04-02 | 中国航天空气动力技术研究院 | 一种旋转飞行器动态测力风洞试验数据的处理系统及方法 |
CN115127768A (zh) * | 2022-09-01 | 2022-09-30 | 中国航空工业集团公司哈尔滨空气动力研究所 | 基于舵面可控模型的飞行器纵向运动模拟风洞试验方法 |
CN115127768B (zh) * | 2022-09-01 | 2022-11-08 | 中国航空工业集团公司哈尔滨空气动力研究所 | 基于舵面可控模型的飞行器纵向运动模拟风洞试验方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106228014B (zh) | 一种导弹气动系数的获取方法 | |
CN109492237A (zh) | 一种气动系数获得方法 | |
CN109297666B (zh) | 一种基于两套运动机构的级间分离风洞试验装置和试验方法 | |
Zhang et al. | Modeling and attitude control of bi-copter | |
CN109635494B (zh) | 一种飞行试验与地面仿真气动力数据综合建模方法 | |
CN109612676B (zh) | 基于飞行试验数据的气动参数反算方法 | |
Pshikhopov et al. | Mathematical model of robot on base of airship | |
CN111695193B (zh) | 一种全局相关三维气动力数学模型的建模方法及系统 | |
CN114608794A (zh) | 一种模型风洞虚拟飞行试验气动力系数测量方法 | |
CN110398339B (zh) | 一种满足多级全自由飞风洞试验相似律的条件确定方法 | |
CN112307683B (zh) | 火箭侧向喷流干扰确定方法、终端及存储介质 | |
CN113848963B (zh) | 一种飞行控制系统的控制律参数设计方法 | |
CN106096091B (zh) | 一种飞机运动模拟方法 | |
CN109141802B (zh) | 用于捕获轨迹试验中的外挂物控制律的模拟方法 | |
CN108225323B (zh) | 基于偏差影响方向组合确定落区边界的方法、介质和设备 | |
CN109190248B (zh) | 一种用于滑翔飞行器的滑翔射程解析方法及解析系统 | |
CN116643578B (zh) | 一种微小型尾座式无人机多模态统一控制方法 | |
Shen et al. | A 6DOF mathematical model of parachute in Mars EDL | |
CN111159957B (zh) | 一种钝体外形过渡流区气动特性的快速预测方法 | |
CN105799949B (zh) | 一种亚轨道卫星的压心设计方法、姿态控制方法和系统 | |
CN110414110B (zh) | 一种用于飞行失速状态下的飞机受力仿真方法 | |
CN108414185B (zh) | 对称飞行器风洞试验数据零点误差处理方法 | |
CN111985169A (zh) | 一种近轴对称飞行器气动特性建模方法 | |
CN110928321A (zh) | 一种四旋翼无人机姿态的鲁棒控制方法 | |
CN114564047A (zh) | 一种考虑气象条件的无人机等速飞行控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190319 |