CN111143987A - 一种飞机高升力系统动力学建模方法 - Google Patents

一种飞机高升力系统动力学建模方法 Download PDF

Info

Publication number
CN111143987A
CN111143987A CN201911347722.3A CN201911347722A CN111143987A CN 111143987 A CN111143987 A CN 111143987A CN 201911347722 A CN201911347722 A CN 201911347722A CN 111143987 A CN111143987 A CN 111143987A
Authority
CN
China
Prior art keywords
lift system
flap
high lift
modeling method
speed reducer
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
CN201911347722.3A
Other languages
English (en)
Other versions
CN111143987B (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.)
Xian Aircraft Design and Research Institute of AVIC
Original Assignee
Xian Aircraft Design and Research Institute of AVIC
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 Xian Aircraft Design and Research Institute of AVIC filed Critical Xian Aircraft Design and Research Institute of AVIC
Priority to CN201911347722.3A priority Critical patent/CN111143987B/zh
Publication of CN111143987A publication Critical patent/CN111143987A/zh
Application granted granted Critical
Publication of CN111143987B publication Critical patent/CN111143987B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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

Abstract

本发明公开了一种飞机高升力系统动力学建模方法,包括步骤1:计算各个作动器之间的所有扭力管的刚度keq;步骤2:建立传动线系中的减速器的数学模型;步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;步骤4:计算每个复合机构及襟翼的转动惯量;步骤5:通过拉格朗日法等动力学分析方法,建立高升力系统动力学模型,本发明解决了飞机高升力系统研制中的动力学问题,对于铁鸟试验中的故障能给出精确的故障机理,在高升力系统设计中,解决了减速器保护门限设计、系统故障载荷计算等关键问题,提高了研制效率,极大地降低了高升力系统的成本,解决了研制瓶颈问题。

Description

一种飞机高升力系统动力学建模方法
技术领域
本发明属于飞机高升力系统动力学设计与验证技术领域,尤其涉及一种飞机高升力系统动力学建模方法。
背景技术
大型喷气式运输机、客机需要高效的高升力系统与增升装置,高升力系统设计是现代飞机的设计内容。由于增升装置的质量往往较大,动力学效应在设计中是不可忽略的,地面试验与飞行中经常会动力学因素引起的故障或问题。我国研制的的ARJ21-700飞机、C919飞机的高升力系统一般采购国外知名供应商货架产品,高升力系统动力学设计问题国外供应商一般都以技术秘密不予告知,这样的飞机高升力系统往往比飞机设计实际需求的高升力系统要重很多,作动部件、驱动部件等载荷保护门限高,高升力系统经济性不好,价格非常昂贵。这样的技术局面,需要我国自主研发的大型喷气式飞机项目牵引,需要在精密复杂高升力系统设计研制中逐步发现技术问题并研究解决。
发明内容
本发明的目的:提出一种飞机高升力系统动力学建模方法,给出大型喷气式飞机高升力系统在铁鸟试验中暴露出来的系统复杂动力学现象的故障机理,对于高升力系统的动力学特征在设计研究、试验、系统使用中的影响等提供技术手段。
本发明的技术方案:
一种飞机高升力系统动力学建模方法,包括以下步骤:
步骤1:计算各个作动器之间的所有扭力管的刚度keq
步骤2:建立传动线系中的减速器的数学模型;
步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;
步骤4:计算每个复合机构的转动惯量J1及襟翼的转动惯量J2
步骤5:通过拉格朗日法等动力学分析方法,建立高升力系统动力学模型。
步骤1所述的计算各个作动器之间的所有扭力管的刚度keq,计算公式为:
Figure BDA0002333857010000021
其中
Figure BDA0002333857010000022
G是扭力管材料的剪切模量,Ii是第i根扭力管横截面的惯性矩,Li是第i根扭力管长度。
步骤2所述的传动线系中的减速器的数学模型为:
Figure BDA0002333857010000023
式中,θin(t)是减速器输入端角位移,n是减速器减速比、θout(t)是减速器输出端角位移,
Figure BDA0002333857010000024
是减速器输入端角加速度,
Figure BDA0002333857010000025
是减速器输出端角加速度。
步骤3所述的组合数学模型为
Figure BDA0002333857010000026
式中θtube(t)是复合机构输入端角位移、θflap(t)是襟翼的角位移、
Figure BDA0002333857010000027
是复合机构输入端角加速度、
Figure BDA0002333857010000028
襟翼的角加速度。
步骤3所述的复合机构包括减速器、运动副及运动机构。
步骤4所述的计算每个复合机构的转动惯量J1,计算公式为:J1=mh×r2,其中mh是运动机构滑轮架质量,r是运动机构滑轮架回转半径。
步骤4所述的计算襟翼的转动惯量J2,计算公式为:
Figure BDA0002333857010000029
其中m为襟翼质量,a是襟翼重心到襟翼一端的运动机构运动平面距离,b是襟翼重心到襟翼另一端的运动机构运动平面距离。
步骤5所述的高升力系统动力学模型为:
Figure BDA00023338570100000210
式中,[m]是转动惯量矩阵,
Figure BDA00023338570100000211
是广义加速度矩阵,[c]是阻尼矩阵,
Figure BDA00023338570100000212
是广义速度矩阵,[k]是刚度矩阵,{q(t)}是广义位移矩阵,{Q(t)}是广义力矩阵。
所述的复合机构输入端角位移θtube(t)的计算公式为:
Figure BDA00023338570100000213
式中,ni是复合机构中减速器的减速比,a、b、c为采用最小二乘法拟合得到的系数。
本发明的有益效果:本发明提供了一种飞机高升力系统动力学建模方法,与传统的有限元方法比,模型自由度非常少,通过数学软件求解及其方便。能解决飞机高升力系统研制中的动力学问题,对于铁鸟试验中的故障能给出精确的故障机理,在高升力系统设计中,解决了减速器保护门限设计、系统故障载荷计算等关键问题,提高了研制效率,极大地降低了高升力系统的成本,解决了研制瓶颈问题。
具体实施方式
一种飞机高升力系统动力学建模方法,包括以下步骤:
步骤1:计算各个作动器之间的所有扭力管的刚度keq,计算公式为:
Figure BDA0002333857010000031
其中
Figure BDA0002333857010000032
G是扭力管材料的剪切模量,Ii是第i根扭力管横截面的惯性矩,Li是第i根扭力管长度。
步骤2:建立传动线系中的减速器的数学模型;
传动线系中的减速器的数学模型为:
Figure BDA0002333857010000033
式中,θin(t)是减速器输入端角位移,n是减速器减速比、θout(t)是减速器输出端角位移,
Figure BDA0002333857010000034
是减速器输入端角加速度,
Figure BDA0002333857010000035
是减速器输出端角加速度。
步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;复合机构包括减速器、运动副及运动机构,组合数学模型为
Figure BDA0002333857010000036
式中θtube(t)是复合机构输入端角位移、θflap(t)是襟翼的角位移、
Figure BDA0002333857010000037
是复合机构输入端角加速度、
Figure BDA0002333857010000038
襟翼的角加速度。
复合机构输入端角位移θtube(t)的计算公式为:
Figure BDA0002333857010000039
式中,ni是复合机构中减速器的减速比,a、b、c为采用最小二乘法拟合得到的系数。
步骤4:计算每个复合机构及襟翼的转动惯量;
每个复合机构的转动惯量,计算公式为:J1=mh×r2,其中mh是运动机构滑轮架质量,r是运动机构滑轮架回转半径。
襟翼的转动惯量,计算公式为:
Figure BDA0002333857010000041
其中m为襟翼质量,a是襟翼重心到襟翼一端的运动机构运动平面距离,b是襟翼重心到襟翼另一端的运动机构运动平面距离。
步骤5:通过拉格朗日法等动力学分析方法,建立高升力系统动力学模型。
高升力系统动力学模型为:
Figure BDA0002333857010000042
式中,[m]是转动惯量矩阵,
Figure BDA0002333857010000043
是广义加速度矩阵,[c]是阻尼矩阵,
Figure BDA0002333857010000044
是广义速度矩阵,[k]是刚度矩阵,{q(t)}是广义位移矩阵,{Q(t)}是广义力矩阵。
本发明提供了一种飞机高升力系统动力学建模方法,与传统的有限元方法比,模型自由度非常少,通过数学软件求解及其方便。能解决飞机高升力系统研制中的动力学问题,对于铁鸟试验中的故障能给出精确的故障机理,在高升力系统设计中,解决了减速器保护门限设计、系统故障载荷计算等关键问题,提高了研制效率,极大地降低了高升力系统的成本,解决了研制瓶颈问题。

Claims (9)

1.一种飞机高升力系统动力学建模方法,其特征在于:包括以下步骤:
步骤1:计算各个作动器之间的所有扭力管的刚度keq
步骤2:建立传动线系中的减速器的数学模型;
步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;
步骤4:计算每个复合机构的转动惯量J1及襟翼的转动惯量J2
步骤5:通过拉格朗日法等动力学分析方法,建立高升力系统动力学模型。
2.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤1所述的计算各个作动器之间的所有扭力管的刚度keq,计算公式为:
Figure FDA0002333855000000011
其中
Figure FDA0002333855000000012
G是扭力管材料的剪切模量,Ii是第i根扭力管横截面的惯性矩,Li是第i根扭力管长度。
3.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤2所述的传动线系中的减速器的数学模型为:
Figure FDA0002333855000000013
式中,θin(t)是减速器输入端角位移,n是减速器减速比、θout(t)是减速器输出端角位移,
Figure FDA0002333855000000014
是减速器输入端角加速度,
Figure FDA0002333855000000015
是减速器输出端角加速度。
4.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤3所述的组合数学模型为
Figure FDA0002333855000000016
式中θtube(t)是复合机构输入端角位移、θflap(t)是襟翼的角位移、
Figure FDA0002333855000000017
是复合机构输入端角加速度、
Figure FDA0002333855000000018
襟翼的角加速度。
5.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤3所述的复合机构包括减速器、运动副及运动机构。
6.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤4所述的计算每个复合机构的转动惯量J1,计算公式为:J1=mh×r2,其中mh是运动机构滑轮架质量,r是运动机构滑轮架回转半径。
7.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤4所述的计算襟翼的转动惯量J2,计算公式为:
Figure FDA0002333855000000021
其中m为襟翼质量,a是襟翼重心到襟翼一端的运动机构运动平面距离,b是襟翼重心到襟翼另一端的运动机构运动平面距离。
8.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤5所述的高升力系统动力学模型为:
Figure FDA0002333855000000022
式中,[m]是转动惯量矩阵,
Figure FDA0002333855000000023
是广义加速度矩阵,[c]是阻尼矩阵,
Figure FDA0002333855000000024
是广义速度矩阵,[k]是刚度矩阵,{q(t)}是广义位移矩阵,{Q(t)}是广义力矩阵。
9.根据权利要求4所述的一种飞机高升力系统动力学建模方法,其特征在于:所述的复合机构输入端角位移θtube(t)的计算公式为:
Figure FDA0002333855000000025
式中,ni是复合机构中减速器的减速比,a、b、c为采用最小二乘法拟合得到的系数。
CN201911347722.3A 2019-12-24 2019-12-24 一种飞机高升力系统动力学建模方法 Active CN111143987B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911347722.3A CN111143987B (zh) 2019-12-24 2019-12-24 一种飞机高升力系统动力学建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911347722.3A CN111143987B (zh) 2019-12-24 2019-12-24 一种飞机高升力系统动力学建模方法

Publications (2)

Publication Number Publication Date
CN111143987A true CN111143987A (zh) 2020-05-12
CN111143987B CN111143987B (zh) 2023-08-04

Family

ID=70519635

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911347722.3A Active CN111143987B (zh) 2019-12-24 2019-12-24 一种飞机高升力系统动力学建模方法

Country Status (1)

Country Link
CN (1) CN111143987B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111639428A (zh) * 2020-05-28 2020-09-08 上海电气风电集团股份有限公司 带襟翼的风力发电机的动力学仿真方法、设备及介质

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0258926A1 (en) * 1986-08-18 1988-03-09 Strijense Kunststof Technieken B.V. Wind turbine rotor with two rotor blades
CA2235307A1 (en) * 1995-10-24 1997-05-01 Hans-Jurgen Bothe Hybrid aircraft
US6755205B1 (en) * 2002-09-12 2004-06-29 Woodward Governor Company Method to stabilize a nozzle flapper valve
US20070041823A1 (en) * 2005-08-22 2007-02-22 Miller Donald C Fluid energy converter
WO2010118919A1 (de) * 2009-04-15 2010-10-21 Airbus Operations Gmbh System und verfahren zur bestimmung von lokalen beschleunigungen, dynamischen lastverteilungen und aerodynamischen daten bei einem luftfahrzeug
CA2830352A1 (en) * 2011-03-30 2012-10-04 The Society Of Japanese Aerospace Companies High-lift device of flight vehicle
CN203111496U (zh) * 2013-01-29 2013-08-07 中国航空工业集团公司西安飞机设计研究所 一种结合同步及异步技术的高升力控制系统
US20140308123A1 (en) * 2013-04-11 2014-10-16 Sikorsky Aircraft Corporation Bearing for rotor blade
JP2015149060A (ja) * 2014-01-10 2015-08-20 株式会社神戸製鋼所 周期的な運動を行う機構を備える機器のモデルの作成方法、及びそのモデルを用いた機器の故障診断方法
US20150307191A1 (en) * 2014-04-28 2015-10-29 Daedalus Flight Systems, LLC Flapping wing aerial vehicles
CN106096088A (zh) * 2016-05-31 2016-11-09 中国航空工业集团公司西安飞机设计研究所 一种螺旋桨飞机螺旋颤振分析方法
CN106197909A (zh) * 2016-06-28 2016-12-07 北京空间飞行器总体设计部 一种谐波减速器传动刚度测试方法
CN106741863A (zh) * 2016-11-17 2017-05-31 中国商用飞机有限责任公司 飞机的高升力系统
CN106845048A (zh) * 2017-04-24 2017-06-13 北京航空航天大学 一种计入摩擦和齿侧间隙的内啮合齿轮轴减速器非线性动力学建模方法
CN107092750A (zh) * 2017-04-24 2017-08-25 北京航空航天大学 一种内啮合齿轮轴减速器轮齿损伤故障的非线性动力学建模方法
CN108275525A (zh) * 2018-01-17 2018-07-13 中北大学 基于动力学分析的高速电梯运行参数监测方法
CN108883824A (zh) * 2016-03-23 2018-11-23 冯春魁 飞行器的数据的获取、处理及飞行状况监控的方法及系统
CA3017202A1 (en) * 2017-11-27 2019-05-27 Goodrich Actuation Systems Limited Improved system for detecting a mechanical fault in a rotating shaft
CN110083994A (zh) * 2019-05-31 2019-08-02 吉林大学 一种行星混联式混合动力系统扭振优化方法
US20190300208A1 (en) * 2018-03-27 2019-10-03 Massachusetts Institute Of Technology Methods and apparatus for in-situ measurements of atmospheric density

Patent Citations (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0258926A1 (en) * 1986-08-18 1988-03-09 Strijense Kunststof Technieken B.V. Wind turbine rotor with two rotor blades
CA2235307A1 (en) * 1995-10-24 1997-05-01 Hans-Jurgen Bothe Hybrid aircraft
US6755205B1 (en) * 2002-09-12 2004-06-29 Woodward Governor Company Method to stabilize a nozzle flapper valve
US20070041823A1 (en) * 2005-08-22 2007-02-22 Miller Donald C Fluid energy converter
WO2010118919A1 (de) * 2009-04-15 2010-10-21 Airbus Operations Gmbh System und verfahren zur bestimmung von lokalen beschleunigungen, dynamischen lastverteilungen und aerodynamischen daten bei einem luftfahrzeug
CA2830352A1 (en) * 2011-03-30 2012-10-04 The Society Of Japanese Aerospace Companies High-lift device of flight vehicle
WO2012132420A1 (ja) * 2011-03-30 2012-10-04 社団法人日本航空宇宙工業会 飛行体の高揚力装置
CN203111496U (zh) * 2013-01-29 2013-08-07 中国航空工业集团公司西安飞机设计研究所 一种结合同步及异步技术的高升力控制系统
US20140308123A1 (en) * 2013-04-11 2014-10-16 Sikorsky Aircraft Corporation Bearing for rotor blade
JP2015149060A (ja) * 2014-01-10 2015-08-20 株式会社神戸製鋼所 周期的な運動を行う機構を備える機器のモデルの作成方法、及びそのモデルを用いた機器の故障診断方法
US20150307191A1 (en) * 2014-04-28 2015-10-29 Daedalus Flight Systems, LLC Flapping wing aerial vehicles
CN108883824A (zh) * 2016-03-23 2018-11-23 冯春魁 飞行器的数据的获取、处理及飞行状况监控的方法及系统
CN106096088A (zh) * 2016-05-31 2016-11-09 中国航空工业集团公司西安飞机设计研究所 一种螺旋桨飞机螺旋颤振分析方法
CN106197909A (zh) * 2016-06-28 2016-12-07 北京空间飞行器总体设计部 一种谐波减速器传动刚度测试方法
CN106741863A (zh) * 2016-11-17 2017-05-31 中国商用飞机有限责任公司 飞机的高升力系统
CN106845048A (zh) * 2017-04-24 2017-06-13 北京航空航天大学 一种计入摩擦和齿侧间隙的内啮合齿轮轴减速器非线性动力学建模方法
CN107092750A (zh) * 2017-04-24 2017-08-25 北京航空航天大学 一种内啮合齿轮轴减速器轮齿损伤故障的非线性动力学建模方法
CA3017202A1 (en) * 2017-11-27 2019-05-27 Goodrich Actuation Systems Limited Improved system for detecting a mechanical fault in a rotating shaft
CN108275525A (zh) * 2018-01-17 2018-07-13 中北大学 基于动力学分析的高速电梯运行参数监测方法
US20190300208A1 (en) * 2018-03-27 2019-10-03 Massachusetts Institute Of Technology Methods and apparatus for in-situ measurements of atmospheric density
CN110083994A (zh) * 2019-05-31 2019-08-02 吉林大学 一种行星混联式混合动力系统扭振优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
朱凯: "某型机内襟翼运动机构方案设计", pages 031 - 46 *
李玉川: "襟翼作动器控制系统设计及安全性评估", pages 031 - 870 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111639428A (zh) * 2020-05-28 2020-09-08 上海电气风电集团股份有限公司 带襟翼的风力发电机的动力学仿真方法、设备及介质
CN111639428B (zh) * 2020-05-28 2023-12-29 上海电气风电集团股份有限公司 带襟翼的风力发电机的动力学仿真方法、设备及介质

Also Published As

Publication number Publication date
CN111143987B (zh) 2023-08-04

Similar Documents

Publication Publication Date Title
CN107544261B (zh) 不确定环境接触下的可重构机器人分散学习最优控制方法
CN101508112B (zh) 三自由度搬运工业机器人多目标优化设计参数的获取方法
CN109991856B (zh) 一种机器人驾驶车辆集成协调控制方法
CN108663940A (zh) 基于集总复合估计的飞行器神经网络学习控制方法
CN105260566B (zh) 一种直升机操纵解耦设计方法
CN205383949U (zh) 一种空气舵与燃气舵的联动结构
CN110889169B (zh) 基于多体系统传递矩阵法的舵面系统非线性颤振模型建模方法
CN106372332B (zh) 一种基于典型工况的轨道交通齿轮箱强度疲劳的评价方法
CN102117362B (zh) 滑流影响下的轻型飞机水平尾翼设计载荷确定方法
CN107292005B (zh) 一种基于模态解耦技术的轨道交通变流器振动优化方法
CN104091485B (zh) 一种双电机驱动的负载模拟器
CN111143987A (zh) 一种飞机高升力系统动力学建模方法
CN111474852B (zh) 一种压电驱动变形机翼的离散滑模控制方法
CN109515748B (zh) 一种飞机襟翼非线性力加载方法
CN108628335A (zh) 一种飞翼无人机的航向控制方法
CN109241553B (zh) 一种机械臂关节实时动力学建模方法
CN109592064B (zh) 飞机与机械操纵系统变形差异对机动操纵影响设计方法
CN106777438B (zh) 直升机典型旋转部件的空间多体运动仿真分析方法
Wu et al. Control of hybrid machines with 2-DOF for trajectory tracking problems
JP6660038B1 (ja) 軸トルク制御装置
Ruan et al. Influence of clearance and structural coupling parameters on shimmy stability of landing gear
CN102673801B (zh) 井字型反扭矩舵结构布局方法
CN106840723A (zh) 液压摇臂加载装置及其加载方法
Grodzki et al. Modelling of UAV's composite structures and prediction of safety factor
CN116540554B (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