CN113971375A - 一种阵风环境下的无人机气动力预测方法 - Google Patents

一种阵风环境下的无人机气动力预测方法 Download PDF

Info

Publication number
CN113971375A
CN113971375A CN202111336188.3A CN202111336188A CN113971375A CN 113971375 A CN113971375 A CN 113971375A CN 202111336188 A CN202111336188 A CN 202111336188A CN 113971375 A CN113971375 A CN 113971375A
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
gust
speed
heave
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
Application number
CN202111336188.3A
Other languages
English (en)
Other versions
CN113971375B (zh
Inventor
李道春
阚梓
赵仕伟
申童
姚卓尔
向锦武
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202111336188.3A priority Critical patent/CN113971375B/zh
Publication of CN113971375A publication Critical patent/CN113971375A/zh
Application granted granted Critical
Publication of CN113971375B publication Critical patent/CN113971375B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10Noise analysis or noise optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Computer Hardware Design (AREA)
  • General Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Geometry (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Fluid Mechanics (AREA)
  • Algebra (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种阵风环境下的无人机气动力预测方法,包括:以无人机的随机运动作为输入,以计算流体力学方法计算的无人机随机运动的气动力响应作为输出,通过深度神经网络,建立阵风环境无人机气动力预测模型。本发明的优点是:能够有效得到阵风环境下无人机气动力情况,从而为相关设计提供指导思想,提高无人机的全天候飞行能力。

Description

一种阵风环境下的无人机气动力预测方法
技术领域
本发明涉及无人机技术领域,特别涉及一种阵风环境下的无人机气动力预测方法。
背景技术
阵风是无人机使用过程中经常遇到的自然现象,阵风使无人机增加额外的附加力矩和附加升力,可能会显著改变结构的动力特性和气动力特性,从而恶化飞行的稳定性操纵性,对无人机的飞行性能和结构安全造成不利影响,将会严重影响长航时无人机全天候飞行的能力。因此阵风响应分析是无人机设计过程中的一个重要问题。传统的方法在计算无人机阵风响应气动力时,多采用估算的方法,但是这与实际情况可能存在较大的差距,尤其是非定常气动力计算,估算的方法往往对升力估算精确度不高。但是如果采用CFD的方法计算气动力,往往又消耗大量时间和计算资源。所以发展和应用一种既可以保证非定常气动力计算精度,又能对非定常气动力进行快速计算的方法非常必要。
发明内容
本发明针对现有技术的缺陷,提供了一种阵风环境下的无人机气动力预测方法。
为了实现以上发明目的,本发明采取的技术方案如下:
一种阵风环境下的无人机气动力预测方法,包括以下步骤:
第一步,根据阵风作用于无人机的速度,采用相对速度法,相应的转换为无人机的升沉运动,模拟阵风环境;设定阵风垂向速度为v(t),无人机垂向获得相对速度b(t)=-v(t)。
第二步,无人机的升沉运动速度b(t)服从高斯分布,而且其功率谱密度均匀分布,因此与高斯白噪声相同,对该高斯白噪声型的无人机随机运动信号进行滤波,滤波范围根据阵风的频率确定,得到滤波后的无人机升沉幅度,如下所示:
Hrand=filter(wag(P))
其中,Hrand表示无人机升沉速度,wag(P)为P dBW强度的高斯白噪声,filter为滤波函数;
采用同样的方法,生成测试集的无人机升沉随机运动速度Hrand 2,其中高斯白噪声的强度小于P;
结合计算流体力学,采用RANS方法,选择Spalart-Allmaras湍流模型进行流场求解,无人机随机升沉运动速度Hrand使用动网格技术进行实现,即通过设置机翼为刚体运动,指定机翼在垂直方向上的速度,该速度为随机升沉运动速度Hrand的关于时间的导数。
第三步,无人机升沉随机运动速度Hrand作为训练样本集的输入,相应的气动力响应作为训练样本集的输出,带入长短期记忆神经网络模型:
it=σ(Wi·yt-1,ut+bi)
Figure BDA0003350579270000021
ft=σ(Wf·yt-1,ut+bf)
Figure BDA0003350579270000022
ot=σ(Wo·yt-1,ut+bo)
yt=ot*tanh(Ct)
Figure BDA0003350579270000023
Figure BDA0003350579270000031
其中,it,Ot和ft分别代表输入门、输出门遗忘门,Ct代表t时刻的单元状态,yt和ut分别代表t时刻的输入和输出;W是权重矩阵,b代表偏置,下标t表示t时刻各模型中的参数,O表示输入门模型中参数,f表示遗忘门中的参数,C表示各单元中的参数;σ和tanh代表激活函数。
经迭代后,可得阵风环境无人机气动力预测的神经网络模型,测试集的无人机升沉随机运动Hrand 2带入神经网络模型,验证训练模型的参数最优。
第四步,以无人机的相对沉浮运动速度b(t)作为输入,通过步骤三训练的气动力神经网络模型,快速得到阵风影响的无人机气动力。
与现有技术相比,本发明的优点在于:
通过考虑阵风环境对于无人机气动力的影响,能够有效得到阵风环境下的无人机飞行气动性能情况,从而为相关设计提供指导思想,提高无人机的全天候飞行能力。
附图说明
图1是本发明实施例流程图;
图2为本发明实施例生成的无人机升沉随机运动;
图3为本发明实施例无人机升沉随机运动的升力系数响应曲线;
图4为本发明实施例通过神经网络预测得到的升力系数曲线。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下根据附图并列举实施例,对本发明做进一步详细说明。
如图1所示,本实施例针对某型无人机进行迎角阶跃型阵风环境下的气动力预测。
无人机以水平速度V前飞,突然受到Ug的微下冲气流的影响,即受到迎角阶跃型阵风。
第一步,由于受到迎角阶跃型阵风,即无人机的运动为以Ug的速度向上运动,升沉运动速度为b(t)=Ug;其中,V=34m/s,Ug=5m/s,攻角阶跃约为8度。
第二步,无人机的随机升沉运动速度为滤波后的高斯白噪声数据,如图2所示;采用同样的方法,生成测试集的升沉随机运动。结合计算流体力学方法,用动网格技术实现无人机随机升沉运动速度Hrand,计算气动力响应,如图3所示。
第三步,无人机升沉随机运动速度作为训练样本集的输入,相应的气动力响应作为训练样本集的输出,经迭代后,可得阵风环境无人机气动力预测的神经网络模型。测试集的无人机升沉随机运动速度带入神经网络模型,验证训练模型的参数最优。
第四步,以无人机以Ug的速度向上运动作为输入,通过步骤第三步训练的气动力神经网络模型,快速得到阵风影响的无人机气动力,如图4所示。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (1)

1.一种阵风环境下的无人机气动力预测方法,其特征在于,包括以下步骤:
第一步,根据阵风作用于无人机的速度,采用相对速度法,相应的转换为无人机的升沉运动,模拟阵风环境;设定阵风垂向速度为v(t),无人机垂向获得相对速度b(t)=-v(t);
第二步,无人机的升沉运动速度b(t)服从高斯分布,而且其功率谱密度均匀分布,因此与高斯白噪声相同,对该高斯白噪声型的无人机随机运动信号进行滤波,滤波范围根据阵风的频率确定,得到滤波后的无人机升沉幅度,如下所示:
Hrand=filter(wag(P))
其中,Hrand表示无人机升沉速度,wag(P)为P dBW强度的高斯白噪声,filter为滤波函数;
采用同样的方法,生成测试集的无人机升沉随机运动速度Hrand2,其中高斯白噪声的强度小于P;
结合计算流体力学,采用RANS方法,选择Spalart-Allmaras湍流模型进行流场求解,无人机随机升沉运动速度Hrand使用动网格技术进行实现,即通过设置机翼为刚体运动,指定机翼在垂直方向上的速度,该速度为随机升沉运动速度Hrand的关于时间的导数;
第三步,无人机升沉随机运动速度Hrand作为训练样本集的输入,相应的气动力响应作为训练样本集的输出,带入长短期记忆神经网络模型:
it=σ(Wi·yt-1,ut+bi);
Figure FDA0003350579260000011
ft=σ(Wf·yt-1,ut+bf);
Figure FDA0003350579260000012
ot=σ(Wo·yt-1,ut+bo);
yt=ot*tanh(Ct);
Figure FDA0003350579260000021
Figure FDA0003350579260000022
其中,it,Ot和ft分别代表输入门、输出门遗忘门,Ct代表t时刻的单元状态,yt和ut分别代表t时刻的输入和输出;W是权重矩阵,b代表偏置,下标t表示t时刻各模型中的参数,O表示输入门模型中参数,f表示遗忘门中的参数,C表示各单元中的参数;σ和tanh代表激活函数;
经迭代后,可得阵风环境无人机气动力预测的神经网络模型,测试集的无人机升沉随机运动Hrand2带入神经网络模型,验证训练模型的参数最优;
第四步,以无人机的相对沉浮运动速度b(t)作为输入,通过步骤三训练的气动力神经网络模型,快速得到阵风影响的无人机气动力。
CN202111336188.3A 2021-11-12 2021-11-12 一种阵风环境下的无人机气动力预测方法 Active CN113971375B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111336188.3A CN113971375B (zh) 2021-11-12 2021-11-12 一种阵风环境下的无人机气动力预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111336188.3A CN113971375B (zh) 2021-11-12 2021-11-12 一种阵风环境下的无人机气动力预测方法

Publications (2)

Publication Number Publication Date
CN113971375A true CN113971375A (zh) 2022-01-25
CN113971375B CN113971375B (zh) 2024-04-09

Family

ID=79589683

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111336188.3A Active CN113971375B (zh) 2021-11-12 2021-11-12 一种阵风环境下的无人机气动力预测方法

Country Status (1)

Country Link
CN (1) CN113971375B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115841062A (zh) * 2023-02-20 2023-03-24 北京航空航天大学 考虑气动结构耦合的大展弦比无人机气动导数预测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB813297A (en) * 1956-05-07 1959-05-13 Franklin Institute A system for alleviating the effects of gusts on aircraft
US20090171634A1 (en) * 2005-12-06 2009-07-02 Airbus Deutschland Gmbh Method for Reconstructing Gusts and Structural Loads at Aircraft, in Particular Passenger Aircraft
US20120059638A1 (en) * 2009-05-14 2012-03-08 Alston Wind, S.L.U. Method And System For Predicting The Occurrence Of A Wind Gust At A Wind Turbine
CN108516101A (zh) * 2018-04-28 2018-09-11 北京航空航天大学 一种主动和被动相结合的固定翼飞机阵风减缓的控制方法
CN112528478A (zh) * 2020-12-04 2021-03-19 中国航空工业集团公司成都飞机设计研究所 一种无人机阵风载荷谱快速编制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB813297A (en) * 1956-05-07 1959-05-13 Franklin Institute A system for alleviating the effects of gusts on aircraft
US20090171634A1 (en) * 2005-12-06 2009-07-02 Airbus Deutschland Gmbh Method for Reconstructing Gusts and Structural Loads at Aircraft, in Particular Passenger Aircraft
US20120059638A1 (en) * 2009-05-14 2012-03-08 Alston Wind, S.L.U. Method And System For Predicting The Occurrence Of A Wind Gust At A Wind Turbine
CN108516101A (zh) * 2018-04-28 2018-09-11 北京航空航天大学 一种主动和被动相结合的固定翼飞机阵风减缓的控制方法
CN112528478A (zh) * 2020-12-04 2021-03-19 中国航空工业集团公司成都飞机设计研究所 一种无人机阵风载荷谱快速编制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
朱龙俊等: "针对阵风干扰的低空无人机定高控制系统设计", 计算机应用与软件, vol. 37, no. 08, 12 August 2020 (2020-08-12), pages 79 - 84 *
赵炜等: "耦合滑流太阳能无人机阵风响应特性研究", 飞行力学, vol. 38, no. 03, 14 March 2020 (2020-03-14), pages 11 - 17 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115841062A (zh) * 2023-02-20 2023-03-24 北京航空航天大学 考虑气动结构耦合的大展弦比无人机气动导数预测方法

Also Published As

Publication number Publication date
CN113971375B (zh) 2024-04-09

Similar Documents

Publication Publication Date Title
CN109871628B (zh) 一种用于评估水陆两栖飞机适航符合性的仿真计算系统和方法
Lee et al. Wind tunnel testing of a helicopter fuselage and rotor in a ship airwake
US20170046968A1 (en) Flight simulation modeling of aircraft dynamic stall aerodynamics
Izraelevitz et al. A novel degree of freedom in flapping wings shows promise for a dual aerial/aquatic vehicle propulsor
CN103895831A (zh) 一种船舶鳍/翼鳍减横摇抗饱和控制装置及其控制方法
Mirzaei et al. Toward design and fabrication of wind-driven vehicles: procedure to optimize the threshold of driving forces
CN108873923B (zh) 应急指挥控制固定翼无人机舰面紧急起飞控制方法
CN113971375A (zh) 一种阵风环境下的无人机气动力预测方法
CN112380623A (zh) 一种水陆两栖飞机飞行品质适航符合性的评估方法
CN105487386A (zh) 一种在载荷布放强扰下的uuv自适应模糊滑模控制方法
CN113687597B (zh) 一种基于二阶波浪漂移力的无人船路径跟踪方法
Haas et al. Large-eddy simulation of airborne wind energy farms
Wharington Heuristic control of dynamic soaring
Viola et al. A CFD-based wing sail optimisation method coupled to a VPP
Wu et al. Analysis and verification of a miniature dolphin-like underwater glider
CN115685764B (zh) 变翼展飞行器任务自适应的抗干扰跟踪控制方法及系统
CN114003052B (zh) 一种基于动态补偿系统的固定翼无人机纵向运动鲁棒自适应控制方法
CN109398594B (zh) 一种水翼船爬浪控制方法
Lee et al. Wind tunnel testing to determine unsteady loads on a helicopter fuselage in a ship airwake
Mok et al. Performance of a windsurfing sail under steady condition
CN115841062B (zh) 考虑气动结构耦合的大展弦比无人机气动导数预测方法
Kubrynski Design of a flapped laminar airfoil for high performance sailplane
Doherty et al. Use of global drag rise boundaries to investigate ill-posed transonic airfoil optimization
Messinger et al. Applications of System Identification with Sparse Bayesian Regression Discovery of Unmodeled Dynamics of an Airship
CN113917938B (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