CN111143987B - 一种飞机高升力系统动力学建模方法 - Google Patents
一种飞机高升力系统动力学建模方法 Download PDFInfo
- Publication number
- CN111143987B CN111143987B CN201911347722.3A CN201911347722A CN111143987B CN 111143987 B CN111143987 B CN 111143987B CN 201911347722 A CN201911347722 A CN 201911347722A CN 111143987 B CN111143987 B CN 111143987B
- Authority
- CN
- China
- Prior art keywords
- lift system
- flap
- speed reducer
- dynamics
- matrix
- 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
Links
Classifications
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种飞机高升力系统动力学建模方法,包括步骤1:计算各个作动器之间的所有扭力管的刚度keq;步骤2:建立传动线系中的减速器的数学模型;步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;步骤4:计算每个复合机构及襟翼的转动惯量;步骤5:通过拉格朗日法等动力学分析方法,建立高升力系统动力学模型,本发明解决了飞机高升力系统研制中的动力学问题,对于铁鸟试验中的故障能给出精确的故障机理,在高升力系统设计中,解决了减速器保护门限设计、系统故障载荷计算等关键问题,提高了研制效率,极大地降低了高升力系统的成本,解决了研制瓶颈问题。
Description
技术领域
本发明属于飞机高升力系统动力学设计与验证技术领域,尤其涉及一种飞机高升力系统动力学建模方法。
背景技术
大型喷气式运输机、客机需要高效的高升力系统与增升装置,高升力系统设计是现代飞机的设计内容。由于增升装置的质量往往较大,动力学效应在设计中是不可忽略的,地面试验与飞行中经常会动力学因素引起的故障或问题。我国研制的的ARJ21-700飞机、C919飞机的高升力系统一般采购国外知名供应商货架产品,高升力系统动力学设计问题国外供应商一般都以技术秘密不予告知,这样的飞机高升力系统往往比飞机设计实际需求的高升力系统要重很多,作动部件、驱动部件等载荷保护门限高,高升力系统经济性不好,价格非常昂贵。这样的技术局面,需要我国自主研发的大型喷气式飞机项目牵引,需要在精密复杂高升力系统设计研制中逐步发现技术问题并研究解决。
发明内容
本发明的目的:提出一种飞机高升力系统动力学建模方法,给出大型喷气式飞机高升力系统在铁鸟试验中暴露出来的系统复杂动力学现象的故障机理,对于高升力系统的动力学特征在设计研究、试验、系统使用中的影响等提供技术手段。
本发明的技术方案:
一种飞机高升力系统动力学建模方法,包括以下步骤:
步骤1:计算各个作动器之间的所有扭力管的刚度keq;
步骤2:建立传动线系中的减速器的数学模型;
步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;
步骤4:计算每个复合机构的转动惯量J1及襟翼的转动惯量J2;
步骤5:通过拉格朗日法等动力学分析方法,建立高升力系统动力学模型。
步骤1所述的计算各个作动器之间的所有扭力管的刚度keq,计算公式为:其中/>G是扭力管材料的剪切模量,Ii是第i根扭力管横截面的惯性矩,Li是第i根扭力管长度。
步骤2所述的传动线系中的减速器的数学模型为:式中,θin(t)是减速器输入端角位移,n是减速器减速比、θout(t)是减速器输出端角位移,/>是减速器输入端角加速度,/>是减速器输出端角加速度。
步骤3所述的组合数学模型为式中θtube(t)是复合机构输入端角位移、θflap(t)是襟翼的角位移、/>是复合机构输入端角加速度、/>襟翼的角加速度。
步骤3所述的复合机构包括减速器、运动副及运动机构。
步骤4所述的计算每个复合机构的转动惯量J1,计算公式为:J1=mh×r2,其中mh是运动机构滑轮架质量,r是运动机构滑轮架回转半径。
步骤4所述的计算襟翼的转动惯量J2,计算公式为:其中m为襟翼质量,a是襟翼重心到襟翼一端的运动机构运动平面距离,b是襟翼重心到襟翼另一端的运动机构运动平面距离。
步骤5所述的高升力系统动力学模型为:
式中,[m]是转动惯量矩阵,是广义加速度矩阵,[c]是阻尼矩阵,/>是广义速度矩阵,[k]是刚度矩阵,{q(t)}是广义位移矩阵,{Q(t)}是广义力矩阵。
所述的复合机构输入端角位移θtube(t)的计算公式为:
式中,ni是复合机构中减速器的减速比,a、b、c为采用最小二乘法拟合得到的系数。
本发明的有益效果:本发明提供了一种飞机高升力系统动力学建模方法,与传统的有限元方法比,模型自由度非常少,通过数学软件求解及其方便。能解决飞机高升力系统研制中的动力学问题,对于铁鸟试验中的故障能给出精确的故障机理,在高升力系统设计中,解决了减速器保护门限设计、系统故障载荷计算等关键问题,提高了研制效率,极大地降低了高升力系统的成本,解决了研制瓶颈问题。
具体实施方式
一种飞机高升力系统动力学建模方法,包括以下步骤:
步骤1:计算各个作动器之间的所有扭力管的刚度keq,计算公式为:其中/>G是扭力管材料的剪切模量,Ii是第i根扭力管横截面的惯性矩,Li是第i根扭力管长度。
步骤2:建立传动线系中的减速器的数学模型;
传动线系中的减速器的数学模型为:式中,θin(t)是减速器输入端角位移,n是减速器减速比、θout(t)是减速器输出端角位移,/>是减速器输入端角加速度,/>是减速器输出端角加速度。
步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;复合机构包括减速器、运动副及运动机构,组合数学模型为式中θtube(t)是复合机构输入端角位移、θflap(t)是襟翼的角位移、/>是复合机构输入端角加速度、/>襟翼的角加速度。
复合机构输入端角位移θtube(t)的计算公式为:
式中,ni是复合机构中减速器的减速比,a、b、c为采用最小二乘法拟合得到的系数。
步骤4:计算每个复合机构及襟翼的转动惯量;
每个复合机构的转动惯量,计算公式为:J1=mh×r2,其中mh是运动机构滑轮架质量,r是运动机构滑轮架回转半径。
襟翼的转动惯量,计算公式为:其中m为襟翼质量,a是襟翼重心到襟翼一端的运动机构运动平面距离,b是襟翼重心到襟翼另一端的运动机构运动平面距离。
步骤5:通过拉格朗日法等动力学分析方法,建立高升力系统动力学模型。
高升力系统动力学模型为:
式中,[m]是转动惯量矩阵,是广义加速度矩阵,[c]是阻尼矩阵,/>是广义速度矩阵,[k]是刚度矩阵,{q(t)}是广义位移矩阵,{Q(t)}是广义力矩阵。
本发明提供了一种飞机高升力系统动力学建模方法,与传统的有限元方法比,模型自由度非常少,通过数学软件求解及其方便。能解决飞机高升力系统研制中的动力学问题,对于铁鸟试验中的故障能给出精确的故障机理,在高升力系统设计中,解决了减速器保护门限设计、系统故障载荷计算等关键问题,提高了研制效率,极大地降低了高升力系统的成本,解决了研制瓶颈问题。
Claims (7)
1.一种飞机高升力系统动力学建模方法,其特征在于:包括以下步骤:
步骤1:计算各个作动器之间的所有扭力管的刚度keq;
步骤2:建立传动线系中的减速器的数学模型;传动线系中的减速器的数学模型为:式中,θin(t)是减速器输入端角位移,n是减速器减速比、θout(t)是减速器输出端角位移,/>是减速器输入端角加速度,/>是减速器输出端角加速度;
步骤3:建立传动线系与增升装置之间的复合机构的组合数学模型;组合数学模型为式中θtube(t)是复合机构输入端角位移、θflap(t)是襟翼的角位移、是复合机构输入端角加速度、/>襟翼的角加速度;
步骤4:计算每个复合机构的转动惯量J1及襟翼的转动惯量J2;
步骤5:通过拉格朗日法动力学分析方法,建立高升力系统动力学模型。
2.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤1所述的计算各个作动器之间的所有扭力管的刚度keq,计算公式为:其中/>G是扭力管材料的剪切模量,Ii是第i根扭力管横截面的惯性矩,Li是第i根扭力管长度。
3.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤3所述的复合机构包括减速器、运动副及运动机构。
4.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤4所述的计算每个复合机构的转动惯量J1,计算公式为:J1=mh×r2,其中mh是运动机构滑轮架质量,r是运动机构滑轮架回转半径。
5.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤4所述的计算襟翼的转动惯量J2,计算公式为:其中m为襟翼质量,a是襟翼重心到襟翼一端的运动机构运动平面距离,b是襟翼重心到襟翼另一端的运动机构运动平面距离。
6.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:步骤5所述的高升力系统动力学模型为:
式中,[m]是转动惯量矩阵,是广义加速度矩阵,[c]是阻尼矩阵,/>是广义速度矩阵,[k]是刚度矩阵,{q(t)}是广义位移矩阵,{Q(t)}是广义力矩阵。
7.根据权利要求1所述的一种飞机高升力系统动力学建模方法,其特征在于:所述的复合机构输入端角位移θtube(t)的计算公式为:
式中,ni是复合机构中减速器的减速比,a、b、c为采用最小二乘法拟合得到的系数。
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 CN111143987A (zh) | 2020-05-12 |
CN111143987B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111639428B (zh) * | 2020-05-28 | 2023-12-29 | 上海电气风电集团股份有限公司 | 带襟翼的风力发电机的动力学仿真方法、设备及介质 |
Citations (16)
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 |
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 | 中国航空工业集团公司西安飞机设计研究所 | 一种结合同步及异步技术的高升力控制系统 |
JP2015149060A (ja) * | 2014-01-10 | 2015-08-20 | 株式会社神戸製鋼所 | 周期的な運動を行う機構を備える機器のモデルの作成方法、及びそのモデルを用いた機器の故障診断方法 |
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 | 吉林大学 | 一种行星混联式混合动力系统扭振优化方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
MX2008002475A (es) * | 2005-08-22 | 2008-09-10 | Viryd Technologies Inc | Convertidor de energia fluida. |
US20140308123A1 (en) * | 2013-04-11 | 2014-10-16 | Sikorsky Aircraft Corporation | Bearing for rotor blade |
US10017248B2 (en) * | 2014-04-28 | 2018-07-10 | University Of Maryland, College Park | Flapping wing aerial vehicles |
US11242162B2 (en) * | 2018-03-27 | 2022-02-08 | Massachusetts Institute Of Technology | Methods and apparatus for in-situ measurements of atmospheric density |
-
2019
- 2019-12-24 CN CN201911347722.3A patent/CN111143987B/zh active Active
Patent Citations (17)
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 |
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 | 中国航空工业集团公司西安飞机设计研究所 | 一种结合同步及异步技术的高升力控制系统 |
JP2015149060A (ja) * | 2014-01-10 | 2015-08-20 | 株式会社神戸製鋼所 | 周期的な運動を行う機構を備える機器のモデルの作成方法、及びそのモデルを用いた機器の故障診断方法 |
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 | 中北大学 | 基于动力学分析的高速电梯运行参数监测方法 |
CN110083994A (zh) * | 2019-05-31 | 2019-08-02 | 吉林大学 | 一种行星混联式混合动力系统扭振优化方法 |
Non-Patent Citations (1)
Title |
---|
李玉川.襟翼作动器控制系统设计及安全性评估.《中国优秀硕士学位论文全文数据库(电子期刊)》.2017,C031-870. * |
Also Published As
Publication number | Publication date |
---|---|
CN111143987A (zh) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2407954B1 (fr) | Dispositif actionneur hexapode | |
JP5877164B2 (ja) | 航空機用トラクタ | |
CN102358430B (zh) | 舰载机捕捉拦阻装置 | |
CN111143987B (zh) | 一种飞机高升力系统动力学建模方法 | |
CN105260566B (zh) | 一种直升机操纵解耦设计方法 | |
CN110561425A (zh) | 绳索驱动的柔性机器人力位混合控制方法及系统 | |
CN102117362B (zh) | 滑流影响下的轻型飞机水平尾翼设计载荷确定方法 | |
CN109711008A (zh) | 一种飞机重心包线计算方法 | |
CN107697272A (zh) | 一种无人机舵面驱动机构 | |
CN108162935A (zh) | 一种轨道车辆自适应阻力的制动控制方法 | |
CN109094811A (zh) | 多点定位锁紧、快速释放机构 | |
CN106965953B (zh) | 一种起落架收放气动载荷模拟装置及其加载方法 | |
CN106840723A (zh) | 液压摇臂加载装置及其加载方法 | |
CN111284722A (zh) | 一种无人机空中回收系统 | |
CN104627828B (zh) | 带载下放装置的脱钩控制方法、设备、系统及工程机械 | |
CN113467252B (zh) | 基于负载广义位移的双无人机吊运系统非线性自适应控制方法 | |
CN103818815B (zh) | 发动机吊装结构 | |
CN113071704B (zh) | 模拟机翼变形的试验方法和系统 | |
CN109029898A (zh) | 风洞试验列车滑动加速装置及加速系统 | |
CN114476123A (zh) | 一种用于软管式空中加油对接模拟试验装置的设计方法 | |
CN105352703B (zh) | 基于锥轮导向的弯刀滑块变角机构及变角方法 | |
CN109814585B (zh) | 近似线性化控制的空间绳系组合体小角度摆动抑制方法 | |
CN107145646A (zh) | 一种计算飞机尾撬支反力的方法 | |
CN110006615A (zh) | 基于惯性原理驱动的高速车辆风洞模型的运行装置 | |
CN105700537A (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 |