CN113221479A - 一种考虑降雨天气的无人机动力学建模方法 - Google Patents
一种考虑降雨天气的无人机动力学建模方法 Download PDFInfo
- Publication number
- CN113221479A CN113221479A CN202110498319.1A CN202110498319A CN113221479A CN 113221479 A CN113221479 A CN 113221479A CN 202110498319 A CN202110498319 A CN 202110498319A CN 113221479 A CN113221479 A CN 113221479A
- Authority
- CN
- China
- Prior art keywords
- rainfall
- force
- unmanned aerial
- aerial vehicle
- respect
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Geometry (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Fluid Mechanics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种考虑降雨天气的无人机动力学建模方法,通过考虑降雨参数、计算雨滴撞击过程中所产生的力和力矩、考虑水膜附着的无人机气动特性,对无人机进行飞行动力学建模;以小扰动线化模型为基础,考虑降雨引起的附加力和力矩对于纵向飞行性能的影响,建立降雨影响的无人机飞行动力学模型。本发明的方法能够有效得到降雨天气下无人机飞行性能情况,从而为相关设计提供指导思想,提高无人机的全天候飞行能力。
Description
技术领域
本发明属于无人机技术领域,具体涉及一种考虑降雨天气的无人机动力学建模方法,可以用于无人机在降雨情况下飞行性能和飞行控制方法的研究。
背景技术
无人机在降雨条件下飞行时会对其飞行性能产生巨大的影响,包括各个舵面、机翼、尾翼的气动性能,从而恶化飞行的稳定性操纵性,增加飞行阻力减小升力,这将会严重影响长航时无人机全天候飞行的能力。因此为了能够应对这一环境,必须认识和了解这种天候条件下对于无人机飞行能力的影响,改善无人机在这种天候条件下的适应性,从而使得无人机的全天候飞行能力得到提升。
发明内容
本发明针对提升无人机全天候飞行能力,提出了一种考虑降雨天气的无人机动力学建模方法,通过计算降雨参数、计算雨滴产生力和力矩、考虑水膜造成气动影响,建立无人机飞行动力学模型。
具体的技术方案为:
一种考虑降雨天气的无人机动力学建模方法,包括以下步骤:
(1)以单位时间内的降雨量评价降雨的严重程度,计算空气中液态水含量等参数;
(2)假设所有雨滴以同一直径和速度下落,根据空气中液态水含量等参数,通过理论方法计算雨滴撞击飞机时产生的力与力矩;
(3)采用计算流体力学方法得到水膜附着以后引起的升力系数变化量、升力线斜率变化量以及阻力系数变化量;
(4)将雨滴对于飞机所造成的撞击力和力矩表示在无人机初始的平衡方程中,并进行线化处理,得到降雨时的纵向运动小扰动线化方程;建立降雨影响的无人机飞行动力学模型。
进一步的详细步骤为:
(1)以单位时间内的降雨量评价降雨的严重程度,计算空气中液态水含量等参数,空气中液态水含量LWC的计算公式如下所示:
其中,ρw为水的密度;R雨滴半径,单位:mm;N0=0.08cm-4,n=41,m=-0.21;
(2)假设所有雨滴以同一直径和速度下落,根据空气中液态水含量等参数,通过理论方法计算雨滴撞击飞机时产生的力与力矩;雨滴撞击飞机产生附加力大小可以表示为:雨滴撞击飞机产生附加力大小可以表示为:
俯仰力矩表达式如下:
(3)采用计算流体力学方法结合工程方法得到水膜附着以后引起的升力系数变化量、升力线斜率变化量以及阻力系数变化量:
(4)将雨滴对于飞机所造成的撞击力和力矩表示在无人机初始的平衡方程中,并进行线化处理,得到降雨时的纵向运动小扰动线化方程;建立降雨影响的无人机飞行动力学模型。
纵向方程组的状态变量为x=[Δu w q Δθ]T;控制变量为c=[ΔδT Δδe]T;
考虑降雨影响后的矩阵系数为:
其中,
其中,Vk和α分别为无人机运动的速度和迎角。Xu是X轴方向的力关于前飞速度的变化率,m是机体质量,Xw是X轴方向的力关于w的变化率;XV是机体坐标系下X方向的速度,是降雨附加力在该方向下分力;Zu是Z轴方向的力关于前飞速度的变化率,是Z轴方向的力关于的变化率;ZV是机体坐标系下Z方向的速度,是降雨附加力在该方向下分力;Mu是X轴方向的力矩关于前飞速度的变化率,是降雨附加力矩关于速度的导数;是X方向的力关于油门开度的导数;是X方向的力关于升降舵偏角的导数;是Z方向的力关于升降舵偏角的导数;是俯仰力矩关于油门开度的导数;是俯仰力矩关于升降舵偏角的导数;Iy是y方向转动惯量。
本发明的方法通过考虑降雨天气对于飞行动力学模型的影响,能够有效得到降雨天气下无人机飞行性能情况,从而为相关设计提供指导思想,提高无人机的全天候飞行能力。
附图说明
图1为本发明流程图;
图2为实施例在不同降雨量下,正向偏转升降舵单位角度后,飞机俯仰角响应曲线;
图3为实施例在不同降雨量下,正向偏转升降舵单位角度后,沿着机体坐标系x轴的飞行速度u响应曲线;
图4为实施例在不同降雨量下,正向偏转升降舵单位角度后,飞机飞行速度V响应曲线。
具体实施方式
下面结合附图和实施例对本发明的考虑降雨天气的飞行动力学建模方法进一步详细说明。
具体步骤如图1所述,本实施例针对某型无人机进行考虑降雨影响的飞行动力学建模。
采用的降雨情况分别为50mm/h降雨率、75mm/h降雨率以及100mm/h降雨率。
其中,ρw为水的密度,R雨滴直径(单位:mm),n、m、N0为实验测定常数。
N0=0.08cm-4,n=41,m=-0.21;
第二步,计算雨滴撞击飞机产生的力与力矩,
第三步,采用工程估算方法得到水膜附着以后引起的升力系数变化量、升力线斜率变化量以及阻力系数变化量;
第四步,以小扰动线化模型为基础,考虑降雨引起的附加力和力矩对于纵向飞行性能的影响:
纵向方程组的状态变量为x=[Δu w q Δθ]T;控制变量为c=[ΔδT Δδe]T;
考虑降雨影响后的矩阵系数为:
经仿真分析,正向偏转升降舵单位角度后,得到无人机俯仰角、沿着机体坐标系x轴的飞行速度u以及飞机飞行速度V的响应曲线。
如图2、图3和图4所示,结果表明当降雨率小于75mm/h时,雨滴所带来的影响是可以克服的,大于75mm/h,则会使无人机出现低头现象,造成无人机不可避免的出现操纵问题。
Claims (4)
1.一种考虑降雨天气的无人机动力学建模方法,其特征在于,包括以下步骤:
(1)以单位时间内的降雨量评价降雨的严重程度,计算空气中液态水含量参数;
(2)假设所有雨滴以同一直径和速度下落,根据空气中液态水含量参数,通过理论方法计算雨滴撞击飞机时产生的力与力矩;
(3)采用计算流体力学方法得到水膜附着以后引起的升力系数变化量、升力线斜率变化量以及阻力系数变化量;
(4)将雨滴对于飞机所造成的撞击力和力矩表示在无人机初始的平衡方程中,并进行线化处理,得到降雨时的纵向运动小扰动线化方程;建立降雨影响的无人机飞行动力学模型。
3.根据权利要求1所述的一种考虑降雨天气的无人机动力学建模方法,其特征在于,步骤(2)中雨滴撞击飞机产生力大小表示为:
[Vxb Vyb Vzb]T表示机体坐标系下的飞行速度;[Vrxb Vryb Vrzb]T表示机体坐标系下的水滴落速;
俯仰力矩表达式如下:
纵向方程组的状态变量为x=[Δu w q Δθ]T;控制变量为c=[ΔδT Δδe]T;其中,Δu是速度变化量,w是z方向速度,q是俯仰角速度,Δθ是俯仰角变化量,ΔδT是油门开度变化量,Δδe是升降舵变化量;
考虑降雨影响后的矩阵系数为:
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110498319.1A CN113221479B (zh) | 2021-05-08 | 2021-05-08 | 一种考虑降雨天气的无人机动力学建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110498319.1A CN113221479B (zh) | 2021-05-08 | 2021-05-08 | 一种考虑降雨天气的无人机动力学建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113221479A true CN113221479A (zh) | 2021-08-06 |
CN113221479B CN113221479B (zh) | 2022-02-01 |
Family
ID=77091775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110498319.1A Active CN113221479B (zh) | 2021-05-08 | 2021-05-08 | 一种考虑降雨天气的无人机动力学建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113221479B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113761811A (zh) * | 2021-09-09 | 2021-12-07 | 北京航空航天大学 | 一种考虑复杂风场的无人机柔性飞行动力学分析方法 |
CN113867399A (zh) * | 2021-10-28 | 2021-12-31 | 北京航空航天大学 | 一种考虑降雨影响的无人机稳定飞行控制方法 |
CN113987894A (zh) * | 2021-12-30 | 2022-01-28 | 北京航空航天大学 | 一种基于降雨天气条件的无人机动力学建模方法 |
CN114065670A (zh) * | 2021-11-30 | 2022-02-18 | 北京航空航天大学 | 一种考虑降雨影响的无人机气动导数快速辨识方法 |
CN114564045A (zh) * | 2022-04-28 | 2022-05-31 | 北京航空航天大学 | 一种考虑降雨和阵风条件的无人机飞行控制律设计方法 |
CN114942595A (zh) * | 2022-07-25 | 2022-08-26 | 西安爱生技术集团有限公司 | 一种考虑降雨影响的无人机阵风响应建模和分析方法 |
CN115828420A (zh) * | 2022-11-21 | 2023-03-21 | 南京航空航天大学 | 一种吸气式高超声速无人机动态可靠性分析方法 |
CN116300668A (zh) * | 2023-05-26 | 2023-06-23 | 北京航空航天大学杭州创新研究院 | 一种针对降雨干扰的四旋翼无人机分层抗干扰控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100222945A1 (en) * | 2009-02-27 | 2010-09-02 | Hans-Gerd Giesseler | Method and device for determining aerodynamic characteristics of an aircraft |
CN107092718A (zh) * | 2017-03-17 | 2017-08-25 | 中国人民解放军陆军航空兵学院 | 一种飞行器飞行中遭遇降雨时的数值模拟方法 |
US20170356925A1 (en) * | 2015-01-15 | 2017-12-14 | Laurent BERDOULAT | Method for correcting the calculation of a flight characteristic of an aeroplane by taking vertical wind into account, method for calculating the drag coefficient |
CN110908278A (zh) * | 2019-11-12 | 2020-03-24 | 北京航空航天大学 | 一种折叠翼飞行器的动力学建模与稳定控制方法 |
-
2021
- 2021-05-08 CN CN202110498319.1A patent/CN113221479B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100222945A1 (en) * | 2009-02-27 | 2010-09-02 | Hans-Gerd Giesseler | Method and device for determining aerodynamic characteristics of an aircraft |
US20170356925A1 (en) * | 2015-01-15 | 2017-12-14 | Laurent BERDOULAT | Method for correcting the calculation of a flight characteristic of an aeroplane by taking vertical wind into account, method for calculating the drag coefficient |
CN107092718A (zh) * | 2017-03-17 | 2017-08-25 | 中国人民解放军陆军航空兵学院 | 一种飞行器飞行中遭遇降雨时的数值模拟方法 |
CN110908278A (zh) * | 2019-11-12 | 2020-03-24 | 北京航空航天大学 | 一种折叠翼飞行器的动力学建模与稳定控制方法 |
Non-Patent Citations (2)
Title |
---|
ZHENLONG WU 等: "Heavy rain effects on aircraft longitudinal stability and control determined from numerical simulation data", 《SAGE》 * |
黄成涛 和 王立新: "降雨对飞机纵向飞行品质的量化影响", 《飞行力学》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113761811A (zh) * | 2021-09-09 | 2021-12-07 | 北京航空航天大学 | 一种考虑复杂风场的无人机柔性飞行动力学分析方法 |
CN113867399A (zh) * | 2021-10-28 | 2021-12-31 | 北京航空航天大学 | 一种考虑降雨影响的无人机稳定飞行控制方法 |
CN113867399B (zh) * | 2021-10-28 | 2024-06-07 | 北京航空航天大学 | 一种考虑降雨影响的无人机稳定飞行控制方法 |
CN114065670B (zh) * | 2021-11-30 | 2024-04-16 | 北京航空航天大学 | 一种考虑降雨影响的无人机气动导数快速辨识方法 |
CN114065670A (zh) * | 2021-11-30 | 2022-02-18 | 北京航空航天大学 | 一种考虑降雨影响的无人机气动导数快速辨识方法 |
CN113987894A (zh) * | 2021-12-30 | 2022-01-28 | 北京航空航天大学 | 一种基于降雨天气条件的无人机动力学建模方法 |
CN113987894B (zh) * | 2021-12-30 | 2022-03-11 | 北京航空航天大学 | 一种基于降雨天气条件的无人机动力学建模方法 |
CN114564045A (zh) * | 2022-04-28 | 2022-05-31 | 北京航空航天大学 | 一种考虑降雨和阵风条件的无人机飞行控制律设计方法 |
CN114942595A (zh) * | 2022-07-25 | 2022-08-26 | 西安爱生技术集团有限公司 | 一种考虑降雨影响的无人机阵风响应建模和分析方法 |
CN114942595B (zh) * | 2022-07-25 | 2022-11-18 | 西安爱生技术集团有限公司 | 一种考虑降雨影响的无人机阵风响应建模和分析方法 |
CN115828420B (zh) * | 2022-11-21 | 2024-03-29 | 南京航空航天大学 | 一种吸气式高超声速无人机动态可靠性分析方法 |
CN115828420A (zh) * | 2022-11-21 | 2023-03-21 | 南京航空航天大学 | 一种吸气式高超声速无人机动态可靠性分析方法 |
CN116300668A (zh) * | 2023-05-26 | 2023-06-23 | 北京航空航天大学杭州创新研究院 | 一种针对降雨干扰的四旋翼无人机分层抗干扰控制方法 |
CN116300668B (zh) * | 2023-05-26 | 2023-08-08 | 北京航空航天大学杭州创新研究院 | 一种针对降雨干扰的四旋翼无人机分层抗干扰控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113221479B (zh) | 2022-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113221479B (zh) | 一种考虑降雨天气的无人机动力学建模方法 | |
CN109871628B (zh) | 一种用于评估水陆两栖飞机适航符合性的仿真计算系统和方法 | |
CN113987894B (zh) | 一种基于降雨天气条件的无人机动力学建模方法 | |
CN105138003B (zh) | 多操纵面无人机直接升力控制方法 | |
CN109703768B (zh) | 一种基于姿态/轨迹复合控制的软式空中加油对接方法 | |
CN109426146A (zh) | 高超声速飞行器的高阶非奇异Terminal滑模控制方法 | |
CN109703769B (zh) | 一种基于预瞄策略的空中加油对接控制方法 | |
CN106970531B (zh) | 倾转翼垂直起降无人机模态转换控制策略确定方法 | |
CN112960101A (zh) | 一种极简超声速飞翼布局飞行器 | |
CN114564045A (zh) | 一种考虑降雨和阵风条件的无人机飞行控制律设计方法 | |
CN110568765A (zh) | 面向攻角跟踪的高超声速飞行器非对称输出受限控制方法 | |
CN111290278B (zh) | 一种基于预测滑模的高超声速飞行器鲁棒姿态控制方法 | |
CN109466778B (zh) | 一种基于姿态参数求导的弹射座椅俯仰横滚姿态控制方法 | |
Tao et al. | Modeling and control of a powered parafoil in wind and rain environments | |
CN111045440A (zh) | 一种高超声速飞行器俯冲段快速滚转控制方法 | |
CN114578861A (zh) | 一种利用阵风环境的无人机飞行控制策略设计方法 | |
CN113191099A (zh) | 一种考虑结冰影响的无人机动力学建模方法 | |
CN115933733A (zh) | 一种固定翼无人机纵向高度速度解耦非线性控制方法 | |
CN115079574A (zh) | 一种挠性高超声速飞行器分布式故障补偿方法 | |
CN113805602B (zh) | 一种考虑阵风影响的无人机飞行高度控制方法 | |
CN113051662B (zh) | 一种基于cfd和datcom的折叠翼尖变体飞行器的气动建模及其性能评估方法 | |
CN113341710B (zh) | 一种飞行器敏捷转弯复合控制方法和应用 | |
CN109634110B (zh) | 一种基于非线性干扰观测器的尾座式无人机过渡段鲁棒控制方法 | |
CN207607638U (zh) | 基于主动流动控制技术的飞翼布局飞行器航向控制装置 | |
CN110989397B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |