CN105468870B - 一种小型螺旋翼的升力的计算方法 - Google Patents

一种小型螺旋翼的升力的计算方法 Download PDF

Info

Publication number
CN105468870B
CN105468870B CN201510980906.9A CN201510980906A CN105468870B CN 105468870 B CN105468870 B CN 105468870B CN 201510980906 A CN201510980906 A CN 201510980906A CN 105468870 B CN105468870 B CN 105468870B
Authority
CN
China
Prior art keywords
lift
screw wing
wing
angular velocity
air
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.)
Active
Application number
CN201510980906.9A
Other languages
English (en)
Other versions
CN105468870A (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.)
China Aerospace Times Electronics Corp
Original Assignee
China Aerospace Times Electronics Corp
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 China Aerospace Times Electronics Corp filed Critical China Aerospace Times Electronics Corp
Priority to CN201510980906.9A priority Critical patent/CN105468870B/zh
Publication of CN105468870A publication Critical patent/CN105468870A/zh
Application granted granted Critical
Publication of CN105468870B publication Critical patent/CN105468870B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明提供了一种小型螺旋翼的升力的计算方法，步骤如下：(1)假设空气角速度为零，测试螺旋翼稳定前的角速度和升力，计算螺旋翼升力系数C；(2)测试稳态下任一时刻螺旋翼的角速度Ωs和升力Fs，确定空气角速度与螺旋翼角速度的比例关系式中的常数α的值；(3)利用龙格‑库塔法计算空气角速度ωn(tn)；(4)利用步骤(3)计算的空气角速度ωn(tn)以及实测的螺旋翼角速度Ω(tn)计算螺旋翼的升力其中ρ为空气密度，N为螺旋翼的叶片数。与一般的螺旋翼升力计算方法相比，该方法能够实时计算动态过程中具有不同叶片数的螺旋翼的升力。

Description

一种小型螺旋翼的升力的计算方法
技术领域
本发明涉及一种螺旋翼升力和转速关系确定方法，属于飞行器控制技术领域。
背景技术
近年来得益于微机电(MEMS)和计算机技术的发展，出现了微小型的姿态和位置测量传感器，使得小型多旋翼无人机摆脱了无法自动控制的困境。旋翼飞行器具有垂直起降、定点悬停和原地转向，机动灵活等其他飞行器所不具有的优点，在一些地理环境复杂，空间狭小，和需要持续关注某一地区的情况下，旋翼飞行器无疑是最合适的选择。多旋翼无人机越来越受到军方、科研人员和大众消费者的重视。螺旋翼是旋翼无人机的动力装置，研究其工作的动态特性有利于飞行器的控制器设计。
旋翼是通过自身的旋转，推动气体向旋翼下方流动，从而产生向上的反作用力。其作用力的大小与转速成正比例关系。目前，对螺旋翼的建模多采用转速的二次方与升力成正比的关系式，该模型是描述螺旋翼在稳定工作状态时角速度和升力的关系，不能反映螺旋翼的具有角加速度时的角速度与升力变化关系。旋翼无人机的控制多采用负反馈的方法，每一次控制输出一次控制信号并且都经过一个控制周期，所以螺旋翼并不能持续地处于稳态工作状态，因此以前的螺旋翼升力模型并不准确。
发明内容
本发明解决的技术问题是：克服现有技术的不足，提供一种小型螺旋翼的升力的计算方法。
本发明的技术方案是：一种小型螺旋翼的升力的计算方法，步骤如下：
(1)假设空气角速度为零，测试螺旋翼稳定前的角速度和升力，计算螺旋翼升力系数C；
(2)测试稳态下任一时刻螺旋翼的角速度Ωs和升力Fs，确定空气角速度与螺旋翼角速度的比例关系式中的常数α的值；
(3)根据步骤(1)、(2)的结果，利用龙格-库塔法计算空气角速度ωn(tn)；具体公式为：
其中惯性环节系数延时系数比例系数Kn＝α·Ω(tn)，α为常数，Ω(tn)为测试的tn时刻的螺旋翼角速度；h 是时间步长，系数ci、ai、bij的取值如下：
(4)利用步骤(3)计算的空气角速度ωn(tn)以及实测的螺旋翼角速度Ω(tn)计算螺旋翼的升力其中ρ为空气密度，N为螺旋翼的叶片数。
所述步骤(1)中螺旋翼从开始转动到转了1/N圈的时间内，空气角速度为零。
所述步骤(2)中的
本发明与现有技术相比的有益效果是：
本发明建立了螺旋翼的动态模型，目的是用于解决螺旋翼在非稳态运动时的升力变化，并且考虑到了螺旋翼具有多个叶片的情况使本模型的应用范围更广。由于螺旋翼转速改变时与升力到达稳态时会经过一个过度过程。过渡过程的越短螺旋翼的升力越快的到达稳态，使得进行一次控制的周期越短。所以本模型能够为控制器带宽的选取提供依据。
附图说明
图1本方法的流程图。
图2仿真曲线和本方法得到的曲线对比。
具体实施方式
下面结合附图1对本发明做详细说明，具体步骤如下：
1、测试螺旋翼在空气角速度为零的情况下螺旋翼的某一角速度Ωo和其升力Fo。螺旋翼的升力系数为C的计算公式为：
螺旋翼从始转动到转了1/N圈的时间内，空气角速度为零。
2、测试螺旋翼在稳态情况下某个转速Ωs时的升力Fs。α为比例系数Kn中的常数，其具体公式为：
3、计算惯性环节系数Tn的值和延迟系数τn的值，其具体公式为：
4、tn时刻的升力Fn等于该时刻螺旋翼的角速度Ω(tn)与空气角速度ωn差值的平方，乘以螺旋翼升力系数C、空气密度ρ和叶片数量的一半。tn时刻螺旋翼升力的计算方法如下：
(1)利用龙格-库塔法计算空气角速度ωn(tn)；具体公式为：
其中惯性环节系数延时系数比例系数Kn＝α·Ω(tn)，α为常数，Ω(tn)为测试的tn时刻的螺旋翼角速度；h 是时间步长，系数ci、ai、bij的取值如下：
(2)tn时刻螺旋翼的升力Fn的具体公式为：
本发明中涉及的测试时的螺旋翼角速度和升力可以利用仿真软件仿真得到，也可以利用目前常用的测试设备进行测试。附图2中灰色曲线(箭头指向的曲线)为仿真的结果，黑色曲线为本方法的结果。传统的升力计算方法只是计算螺旋翼稳态时的升力，忽略了到达稳态之前的过度过程，本方补充了传统方法在这一方面的欠缺并且适用于具有N个叶片的螺旋翼。
本发明未详细说明部分属于本领域技术人员公知常识。

Claims (3)

1.一种小型螺旋翼的升力的计算方法，其特征在于步骤如下：
(1)假设空气角速度为零，测试螺旋翼稳定前的角速度Ωo和升力Fo，计算螺旋翼升力系数C；
(2)测试稳态下任一时刻螺旋翼的角速度Ωs和升力Fs，确定空气角速度与螺旋翼角速度的比例关系式中的常数α的值；
(3)根据步骤(1)、(2)的结果，利用龙格-库塔法计算空气角速度ωn(tn)；具体公式为：
其中惯性环节系数延时系数比例系数Kn＝α·Ω(tn)，α为常数，Ω(tn)为测试的tn时刻的螺旋翼角速度；h是时间步长，系数ci、ai、bij的取值如下：
(4)利用步骤(3)计算的空气角速度ωn(tn)以及实测的螺旋翼角速度Ω(tn)计算螺旋翼的升力其中ρ为空气密度，N为螺旋翼的叶片数。
2.根据权利要求1所述的一种小型螺旋翼升力的计算方法，其特征在于：所述步骤(1)中螺旋翼从开始转动到转了1/N圈的时间内，空气角速度为零。
3.根据权利要求1所述的一种小型螺旋翼升力的计算方法，其特征在于：所述步骤(2)中的
CN201510980906.9A 2015-12-23 2015-12-23 一种小型螺旋翼的升力的计算方法 Active CN105468870B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510980906.9A CN105468870B (zh) 2015-12-23 2015-12-23 一种小型螺旋翼的升力的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510980906.9A CN105468870B (zh) 2015-12-23 2015-12-23 一种小型螺旋翼的升力的计算方法

Publications (2)

Publication Number Publication Date
CN105468870A CN105468870A (zh) 2016-04-06
CN105468870B true CN105468870B (zh) 2018-11-23

Family

ID=55606565

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510980906.9A Active CN105468870B (zh) 2015-12-23 2015-12-23 一种小型螺旋翼的升力的计算方法

Country Status (1)

Country Link
CN (1) CN105468870B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109292110A (zh) * 2018-08-23 2019-02-01 淮阴工学院 一种翼伞纵向气动系数估算系统及其估算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103853893A (zh) * 2014-03-26 2014-06-11 中国民用航空总局第二研究所 一种从点云数据检测飞机姿态的参数模型匹配方法
CN104401480A (zh) * 2014-11-06 2015-03-11 南京航空航天大学 涵道式倾转飞行器
CN104504255A (zh) * 2014-12-11 2015-04-08 北京航天控制仪器研究所 一种螺旋翼升力和阻力力矩的确定方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6854344B2 (en) * 2002-12-19 2005-02-15 The Johns Hopkins University Sample retrieval device for aerosol collection
US8774982B2 (en) * 2010-08-26 2014-07-08 Leptron Industrial Robotic Helicopters, Inc. Helicopter with multi-rotors and wireless capability

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103853893A (zh) * 2014-03-26 2014-06-11 中国民用航空总局第二研究所 一种从点云数据检测飞机姿态的参数模型匹配方法
CN104401480A (zh) * 2014-11-06 2015-03-11 南京航空航天大学 涵道式倾转飞行器
CN104504255A (zh) * 2014-12-11 2015-04-08 北京航天控制仪器研究所 一种螺旋翼升力和阻力力矩的确定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《旋翼流场流动机理研究》;陈文轩;《直升机技术》;20111231(第2期);第1-14页 *

Also Published As

Publication number Publication date
CN105468870A (zh) 2016-04-06

Similar Documents

Publication Publication Date Title
Lu et al. Real-time simulation system for UAV based on Matlab/Simulink
CN109614633A (zh) 一种复合式旋翼飞行器非线性建模方法及配平方法
CN109446582B (zh) 一种考虑地球自转的高精度降阶平稳滑翔动力学建模方法
CN105242679A (zh) 一种四旋翼飞行器的控制系统设计方法
Wei et al. Frequency-domain system identification and simulation of a quadrotor controller
Gonzalez-Rocha et al. Measuring atmospheric winds from quadrotor motion
Scanavino Design and testing methodologies for UAVs under extreme environmental conditions
Zhao et al. Design and implementation of an innovative airborne electric propulsion measure system of fixed-wing UAV
Zhao Development of a dynamic model of a ducted fan VTOL UAV
Ostler et al. Performance flight testing of small, electric powered unmanned aerial vehicles
Kumar et al. Novel wireless sensing platform for experimental mapping and validation of ship air wake
CN105468870B (zh) 一种小型螺旋翼的升力的计算方法
CN211685678U (zh) 一种多旋翼无人机实时尾迹的仿真分析系统
Zhao et al. Geometry shape selection of NACA airfoils for Mars rotorcraft
Saderla et al. Parameter estimation of UAV from flight data using neural network
CN112441253A (zh) 一种多旋翼无人机实时尾迹的仿真分析系统及方法
CN105259760A (zh) 一种小型四轴飞行器控制器设计方法
Geldenhuys Aerodynamic development of a contra-rotating shrouded rotor system for a UAV
Yeo et al. Experimental characterization of lift on a rigid flapping wing
Boyet ESWIRP: European strategic wind tunnels improved research potential program overview
Prasad et al. Finite state inflow models for a coaxial rotor in hover
Pérez et al. Free-Vortex Wake and CFD Simulation of a Small Rotor for a Quadcopter at Hover
Woodbury et al. Flight test results of Observer/Kalman Filter Identifi [|# 12#|] cation of the Pegasus unmanned vehicle
Thu et al. Enhanced Performance Prediction of Hydrogen Fuel Cell Powered eVTOL UAV
KAYA Estimation of Aerodynamic Loads of a Propeller through Improved Blade Element and Momentum Theory and Propeller Design Optimization

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