CN113962028A - 一种考虑结冰影响的无人机气动导数辨识方法 - Google Patents

一种考虑结冰影响的无人机气动导数辨识方法 Download PDF

Info

Publication number
CN113962028A
CN113962028A CN202111336221.2A CN202111336221A CN113962028A CN 113962028 A CN113962028 A CN 113962028A CN 202111336221 A CN202111336221 A CN 202111336221A CN 113962028 A CN113962028 A CN 113962028A
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
icing
wing
parameters
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
CN202111336221.2A
Other languages
English (en)
Other versions
CN113962028B (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 CN202111336221.2A priority Critical patent/CN113962028B/zh
Publication of CN113962028A publication Critical patent/CN113962028A/zh
Application granted granted Critical
Publication of CN113962028B publication Critical patent/CN113962028B/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/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • 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/045Combinations of networks
    • 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)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Algebra (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Computational Mathematics (AREA)
  • Medical Informatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种考虑结冰影响的无人机气动导数辨识方法,步骤1、确定大气参数和飞行参数的上限和下限;随机选取参数值,建立参数的随机数据集,计算外流场和水滴撞击特性,得到机翼结冰后的形状;步骤2、以随机数据集作为输入,进行归一化处理,采用卷积神经网络,构建结冰机翼形状预测模型;步骤3、指定大气参数和飞行参数,通过结冰机翼形状预测模型,得到结冰机翼形状;步骤4、采用计算流体力学方法,将常用翼型作为测试用例,采用RANS方法对动态气动性能预测能力。结冰机翼的无人机非定常受迫振荡运动,计算无人机的动态气动性能,辨识无人机的气动导数。本发明的优点是:有效得到无人机飞行性能情况,从而为相关设计提供指导思想。

Description

一种考虑结冰影响的无人机气动导数辨识方法
技术领域
本发明涉及无人机技术领域,特别涉及一种考虑结冰影响的无人机气动导数辨识方法。
背景技术
机翼结冰问题是影响航空飞行器飞行安全的重要隐患之一。无人机积冰是由于环境条件和飞行条件共同的影响结果。当发生积冰时,若处理不妥当,轻则恶化无人机的飞行性能、降低飞行品质,重则会发生无人机坠毁的事故,造成严重的损失。当无人机发生积冰现象时,会使机体总体质量增加,同时无人机尤其是机翼表面的气动外形发生变化,将严重影响无人机执行长航时全天候飞行任务的能力。实现长航时无人机全天候飞行,积冰现象是十分容易出现而且必须克服的一种天候条件。
传统的方法在计算无人机结冰机翼的气动导数时,多采用估算的方法,与实际情况可能存在较大的差距。采用试验的方法得到结冰机翼的气动导数时,成本较高。
发明内容
本发明针对现有技术的缺陷,提供了一种考虑结冰影响的无人机气动导数辨识方法。能对非定常气动导数进行快速计算的方法非常必要,对于提高积冰条件下无人机的飞行性能,发展无人机全天候飞行能力,具有重要而深远的意义。
为了实现以上发明目的,本发明采取的技术方案如下:
一种考虑结冰影响的无人机气动导数辨识方法,包括以下步骤:
步骤1、以飞行速度V、飞行攻角α、温度T、水滴直径D、液态水含量LWC,结冰时长T为参数,根据无人机对应的飞行状态和大气状态,选择各个参数的上限和下限;在各参数的范围内,随机选取参数值,建立参数的随机数据集,通过计算流体力学方法,计算外流场和水滴撞击特性计算,最后得到随机数据集对应的机翼结冰后的形状;
步骤2、以随机数据集作为输入,对应的机翼结冰后的形状作为输出,其中,对输入数据进行归一化处理,输出的机翼形状进校灰度化处理,采用卷积神经网络,训练更新神经网络卷积层和全连接层的权重,偏置相关参数,构建结冰机翼形状预测模型;
步骤3、指定飞行速度V、飞行攻角α、温度T、水滴直径D、液态水含量LWC和结冰时长T参数,通过结冰机翼形状预测模型,得到结冰机翼形状;
步骤4、采用计算流体力学方法,结合动网格技术,将常用翼型作为测试用例,采用RANS方法对动态气动性能预测能力。结冰机翼的无人机非定常受迫振荡运动,计算无人机的动态气动性能,辨识无人机的气动导数。
进一步地,所述计算无人机的动态气动性能具体为:通过俯仰阻尼动导数计算方法,结冰机翼受迫振荡运动采用随机运动形式,计算得到结冰机翼的非定常俯仰力矩系数Cm_ice,其表达式如下:
Figure BDA0003350584150000021
其中,α为攻角,q为俯仰角速度,Cm0_ice为零攻角时结冰机翼的俯仰力矩系数,
Figure BDA0003350584150000022
Cmα_ice
Figure BDA0003350584150000023
分别为静导数和组合动导数。
与现有技术相比,本发明的优点在于:
通过考虑结冰对于无人机气动导数的影响,能够有效得到无人机飞行性能情况,从而为相关设计提供指导思想,提高无人机的全天候飞行能力。
附图说明
图1是本发明实施例考虑结冰影响的无人机气动导数辨识方法流程图;
图2是本发明实施例结冰机翼形状预测模型建模流程图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下根据附图并列举实施例,对本发明做进一步详细说明。
如图1所示,一种考虑结冰影响的无人机气动导数辨识方法,包括以下步骤:
步骤1、以飞行速度V、飞行攻角α、温度T、水滴直径D、液态水含量LWC,结冰时长T为参数,根据无人机对应的飞行状态和大气状态,选择各个参数的上限和下限;在各参数的范围内,随机选取参数值,建立参数的随机数据集,通过计算流体力学方法,计算外流场和水滴撞击特性计算,最后得到随机数据集对应的机翼结冰后的形状;
步骤2、以随机数据集作为输入,对应的机翼结冰后的形状作为输出,其中,对输入数据进行归一化处理,输出的机翼形状进校灰度化处理,采用卷积神经网络,训练更新神经网络相关参数,构建结冰机翼形状预测模型,如图2所示;
步骤3、指定大气参数(飞行速度V、飞行攻角)和飞行参数(温度T、水滴直径D、液态水含量LWC,结冰时长T),通过结冰机翼形状预测模型,得到结冰机翼形状;
步骤4、采用计算流体力学方法,结合动网格技术,将常用翼型作为测试用例,采用RANS方法对动态气动性能预测能力,用于后续无人机的动态气动性能计算。结冰机翼的无人机非定常受迫振荡运动,计算无人机的动态气动性能,辨识无人机的气动导数。以俯仰阻尼动导数计算方法为例,结冰机翼受迫振荡运动采用随机运动形式,计算得到结冰机翼的非定常俯仰力矩系数Cm_ice,其表达式如下:
Figure BDA0003350584150000041
其中,α为攻角,q为俯仰角速度,Cm0_ice为零攻角时结冰机翼的俯仰力矩系数,
Figure BDA0003350584150000042
Cmα_ice
Figure BDA0003350584150000043
分别为静导数和组合动导数。
计算得到气动导数对比如下表所示,静导数C相比不考虑降雨影响绝对值减小了24.3%,组合动导数
Figure BDA0003350584150000044
相比不考虑降雨影响绝对值减小了13.3%。
Figure BDA0003350584150000045
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的实施方法,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (2)

1.一种考虑结冰影响的无人机气动导数辨识方法,其特征在于,包括以下步骤:
步骤1、以飞行速度V、飞行攻角α、温度T、水滴直径D、液态水含量LWC,结冰时长T为参数,根据无人机对应的飞行状态和大气状态,选择各个参数的上限和下限;在各参数的范围内,随机选取参数值,建立参数的随机数据集,通过计算流体力学方法,计算外流场和水滴撞击特性计算,最后得到随机数据集对应的机翼结冰后的形状;
步骤2、以随机数据集作为输入,对应的机翼结冰后的形状作为输出,其中,对输入数据进行归一化处理,输出的机翼形状进校灰度化处理,采用卷积神经网络,训练更新神经网络卷积层和全连接层的权重和偏置参数,构建结冰机翼形状预测模型;
步骤3、指定飞行速度V、飞行攻角α、温度T、水滴直径D、液态水含量LWC和结冰时长T参数,通过结冰机翼形状预测模型,得到结冰机翼形状;
步骤4、采用计算流体力学方法,结合动网格技术,将常用翼型作为测试用例,采用RANS方法对动态气动性能预测能力;结冰机翼的无人机非定常受迫振荡运动,计算无人机的动态气动性能,辨识无人机的气动导数。
2.根据权利要求1所述的一种考虑结冰影响的无人机气动导数辨识方法,其特征在于:所述计算无人机的动态气动性能具体为:通过俯仰阻尼动导数计算方法,结冰机翼受迫振荡运动采用随机运动形式,计算得到结冰机翼的非定常俯仰力矩系数Cm_ice,其表达式如下:
Figure FDA0003350584140000011
其中,α为攻角,q为俯仰角速度,Cm0_ice为零攻角时结冰机翼的俯仰力矩系数,
Figure FDA0003350584140000021
Cmα_ice
Figure FDA0003350584140000022
分别为静导数和组合动导数。
CN202111336221.2A 2021-11-12 2021-11-12 一种考虑结冰影响的无人机气动导数辨识方法 Active CN113962028B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111336221.2A CN113962028B (zh) 2021-11-12 2021-11-12 一种考虑结冰影响的无人机气动导数辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111336221.2A CN113962028B (zh) 2021-11-12 2021-11-12 一种考虑结冰影响的无人机气动导数辨识方法

Publications (2)

Publication Number Publication Date
CN113962028A true CN113962028A (zh) 2022-01-21
CN113962028B CN113962028B (zh) 2024-04-30

Family

ID=79470264

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111336221.2A Active CN113962028B (zh) 2021-11-12 2021-11-12 一种考虑结冰影响的无人机气动导数辨识方法

Country Status (1)

Country Link
CN (1) CN113962028B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130181093A1 (en) * 2011-01-20 2013-07-18 Nobuo Inoue Anti-icing system, wing, aircraft, and anti-icing method
KR20130115563A (ko) * 2012-04-12 2013-10-22 한국항공우주산업 주식회사 에어포일 임계 결빙형상을 결정하는 방법
KR101486785B1 (ko) * 2014-01-27 2015-02-04 부경대학교 산학협력단 강수형태와 얼음상황의 진단을 위한 강수온도 감지기
CN109558650A (zh) * 2018-11-09 2019-04-02 中国直升机设计研究所 直升机旋翼结冰对旋翼性能影响的分析方法
CN111738481A (zh) * 2020-04-01 2020-10-02 南京航空航天大学 基于bp神经网络的飞机结冰气象参数mvd预测方法
CN111976996A (zh) * 2020-08-26 2020-11-24 四川大学 一种无人机机翼分区防冰方法
CN113191099A (zh) * 2021-04-28 2021-07-30 北京航空航天大学 一种考虑结冰影响的无人机动力学建模方法
CN113434978A (zh) * 2021-06-28 2021-09-24 成都安世亚太科技有限公司 一种航空飞行器结冰可靠性智能评估方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130181093A1 (en) * 2011-01-20 2013-07-18 Nobuo Inoue Anti-icing system, wing, aircraft, and anti-icing method
KR20130115563A (ko) * 2012-04-12 2013-10-22 한국항공우주산업 주식회사 에어포일 임계 결빙형상을 결정하는 방법
KR101486785B1 (ko) * 2014-01-27 2015-02-04 부경대학교 산학협력단 강수형태와 얼음상황의 진단을 위한 강수온도 감지기
CN109558650A (zh) * 2018-11-09 2019-04-02 中国直升机设计研究所 直升机旋翼结冰对旋翼性能影响的分析方法
CN111738481A (zh) * 2020-04-01 2020-10-02 南京航空航天大学 基于bp神经网络的飞机结冰气象参数mvd预测方法
CN111976996A (zh) * 2020-08-26 2020-11-24 四川大学 一种无人机机翼分区防冰方法
CN113191099A (zh) * 2021-04-28 2021-07-30 北京航空航天大学 一种考虑结冰影响的无人机动力学建模方法
CN113434978A (zh) * 2021-06-28 2021-09-24 成都安世亚太科技有限公司 一种航空飞行器结冰可靠性智能评估方法

Also Published As

Publication number Publication date
CN113962028B (zh) 2024-04-30

Similar Documents

Publication Publication Date Title
US11401044B2 (en) Method and assistance system for detecting a degradation of flight performance
CN114065399B (zh) 一种考虑复杂气象条件下的无人飞行器飞行性能计算方法
CN113191099A (zh) 一种考虑结冰影响的无人机动力学建模方法
CN105314094A (zh) 飞机操纵面的闭环控制
US20170121028A1 (en) Enhancing engine performance to improve fuel consumption based on atmospheric ice particles
CN114065398B (zh) 一种大展弦比柔性飞行器飞行性能计算方法
CN114065670A (zh) 一种考虑降雨影响的无人机气动导数快速辨识方法
CN110641726A (zh) 一种快速确定飞机机翼安装角的方法
CN112380792A (zh) 一种水陆两栖飞机飞行性能适航符合性的评估方法
CN114065407B (zh) 一种基于机场风速多元感知的尾流间隔动态缩减方法
CN112380623A (zh) 一种水陆两栖飞机飞行品质适航符合性的评估方法
CN114004021B (zh) 用于飞行管理系统性能管理的巡航燃油流量计算方法
CN111976996A (zh) 一种无人机机翼分区防冰方法
CN113962028A (zh) 一种考虑结冰影响的无人机气动导数辨识方法
CN109614644B (zh) 一种外吹式襟翼布局飞机动力增升效果评估方法
CN104615863B (zh) 一种带控制面的三自由度机翼的颤振边界预测方法
CN112731958B (zh) 一种基于速度保护的机轮承载信号使用方法
CN112989725A (zh) 一种用于飞机结冰环境模拟的仿真方法
CN113221237A (zh) 一种基于降阶建模的大迎角颤振分析方法
CN114912284A (zh) 基于第一法则的飞行管理系统爬升性能预测算法
CN112526548B (zh) 一种基于测风激光雷达的降雨识别方法及装置
CN110963046A (zh) 一种确定无人机翼面结冰条件的方法
CN109871562B (zh) 一种基于Catia二次开发计算飞机水上迫降静水面漂浮特性的方法
Clarke et al. Trajectory generation for autonomous soaring UAS
CN115200826B (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