CN104933259A - 一种飞机高升力系统载荷计算方法 - Google Patents

一种飞机高升力系统载荷计算方法 Download PDF

Info

Publication number
CN104933259A
CN104933259A CN201510367111.0A CN201510367111A CN104933259A CN 104933259 A CN104933259 A CN 104933259A CN 201510367111 A CN201510367111 A CN 201510367111A CN 104933259 A CN104933259 A CN 104933259A
Authority
CN
China
Prior art keywords
centerdot
load
hitch
lift system
point
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
CN201510367111.0A
Other languages
English (en)
Other versions
CN104933259B (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 CN201510367111.0A priority Critical patent/CN104933259B/zh
Publication of CN104933259A publication Critical patent/CN104933259A/zh
Application granted granted Critical
Publication of CN104933259B publication Critical patent/CN104933259B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Complex Calculations (AREA)

Abstract

本发明涉及航空系统设计领域,具体涉及一种飞机高升力系统载荷计算方法,以解决现有计算方法效率低,无法满足设计要求问题。高升力系统中的增升装置与悬挂机构多连接点连接,悬挂机构与高升力系统中的多个运动机构点接触,依据传力路线上各环节的受力与约束情况,可以把传力路线分解为有机联系的增升装置-悬挂机构、悬挂机构-运动机构子系统,再基于各子系统设计理念和力学特征,分步解决复杂系统的传力计算问题,计算速度快、计算结果准确可靠,可以完全避免有限元软件所需的复杂而又冗长的过程。

Description

一种飞机高升力系统载荷计算方法
技术领域
本发明涉及航空系统设计领域,具体涉及一种飞机高升力系统载荷计算方法。
背景技术
现有技术中,采用商用有限元软件NASTRAN或MARC进行飞机高升力系统载荷计算,但是,对于设计载荷情况多而复杂的高升力系统,以往的方法计算效率低,甚至不能完成高升力系统的载荷计算工作,制约着我国飞机系统设计。
发明内容
本发明的目的是提供一种飞机高升力系统载荷计算方法,以解决现有计算方法效率低,无法满足设计要求问题。
本发明的技术方案是:
一种飞机高升力系统载荷计算方法,所述高升力系统中的增升装置与悬挂机构多连接点连接,所述悬挂机构与所述高升力系统中的多个运动机构点接触,包括如下步骤:
步骤一、将所述高升力系统的传力路线分解为增升装置-悬挂机构、悬挂机构-运动机构;
步骤二、为所述高升力系统施加预定载荷,使所述增升装置受力,并将该力通过所述多连接点传递至所述悬挂机构;
步骤三、获取所述多连接点的刚度,根据所述预定载荷和所述多连接点的刚度得到各连接点的载荷,从而得到所述悬挂机构的总载荷;
步骤四、根据所述悬挂机构与所述多个运动机构的力学平衡方程,得到所述悬挂机构与每一个所述运动机构的接触点载荷。
可选地,在所述步骤三中:
首先根据最少势能原理公式(1):
δ Π = ∂ Π ∂ a 1 δa 1 + ∂ Π ∂ a 2 δa 2 + ... + ∂ Π ∂ a n δa n = 0 - - - ( 1 ) ;
得到公式(2):
∂ Π ∂ a ≡ K a - P = 0 - - - ( 2 ) ,
其中,P为所述预定载荷,K为各所述连接点的刚度矩阵;
再取位移泛函(3)为:
u ≈ u ~ = Σ i = 1 n N i a i = N a - - - ( 3 ) ,
其中,N为形函数,u为位移泛函,i为连接点个数;
再根据公式(4)得到各连接点的载荷Π:
Π = 1 2 u T K u - - - ( 4 ) ;
其中,uT为u的转置矩阵。
可选地,在所述步骤四中,所述悬挂机构与所述运动机构的接触点为三个,对应的载荷分别为滑轨上前滚轮载荷Fq、为滑轨上后前滚轮载荷Fh、啮合齿廓之间的正压力(缝翼)Ft,所述悬挂机构与所述运动机构的力学平衡方程为:
F q · s i n α + μ 1 F q · s i n ( π 2 + α ) + F h · s i n β + μ 2 F h · s i n ( π 2 + β ) + F t · s i n θ - μ c o s θ · F t + F Y = 0.0 F q · c o s α + μ 1 F q · c o s ( π 2 + α ) + F h · c o s β + μ 2 F h · c o s ( π 2 + β ) + F t c o s θ - μ s i n θ · F t + F Z = 0.0 R h g × ( 0.07 2.55 F q + 0.07 1.80 F h ) + F t · L - μ · F t · L 1 + M x = 0.0
其中,FY、fz、mx为所述悬挂机构的总载荷在X、Y、Z三个方向上的分载荷,α是Fq的方向角,β为Fh的方向角,θ为Ft的方向角,Rhg为滑轨外径,L为Ft到滑轨轴线距离,L1为齿轮与齿条间摩擦力到滑轨轴线距离,μ为齿轮摩擦系统。
本发明的有益效果:
本发明的高升力系统载荷计算方法,依据传力路线上各环节的受力与约束情况,可以把传力路线分解为有机联系的增升装置-悬挂机构、悬挂机构-运动机构子系统,再基于各子系统设计理念和力学特征,分步解决复杂系统的传力计算问题,计算速度快、计算结果准确可靠,可以完全避免有限元软件所需的复杂而又冗长的过程。
附图说明
图1为本发明所述实施例的流程图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。
本发明提供了一种飞机高升力系统载荷计算方法;其中,高升力系统中的增升装置与悬挂机构通过多连接点连接,悬挂机构与高升力系统中的多个运动机构点接触。
本发明的方法包括如下步骤:
步骤S101、将所述高升力系统的传力路线分解为增升装置-悬挂机构、悬挂机构-运动机构。
增升装置-高升力系统是复杂流场下的一个复杂传力系统,通常是气动外载-翼面(缝翼或前中后襟翼)-悬挂机构-非线性环节(运动机构)和驱动装置构成了一条可分级研究的传力路线。依据传力路线上各环节的受力与约束情况,可以把传力路线分解为有机联系的至少两个子系统,即增升装置-悬挂机构、悬挂机构-运动机构子系统。
后续再基于各子系统设计理念和力学特征,应用力学变分原理或二次规划进行力学分析,从而分步解决复杂系统的传力计算问题,计算速度快、计算结果准确可靠,可以完全避免有限元软件所需的复杂而又冗长的过程。。
步骤S102、为所述高升力系统施加预定载荷,使所述增升装置受力,并将该力通过所述多连接点传递至所述悬挂机构。
步骤S103、获取所述多连接点的刚度,根据所述预定载荷和所述多连接点的刚度得到各连接点的载荷,从而得到所述悬挂机构的总载荷。
由热力学第一定律,如令物体的动能为Ek,应变能为U,则在微小的δt时间间隔内,物体从一种状态过渡到另一种状态时,总能量的变化为:
δEk+δU=δW+δQ,
其中,δW为体力Fbi与面力pi所完成的功,δQ为物体由其周围介质所吸收(或向外散发)的热量。
假定弹性变形过程是绝热的,则对于静力平衡问题有δQ=0,δEK=0,δU=δW。
由虚位移原理:在外力作用下处于平衡状态的变形体,当给予物体微小虚位移时,外力的总虚功等于物体的总虚应变能,即δU=δW。
由变分原理,微分方程的边值问题可化为一泛函的极值问题:
由里兹法,取位移泛函为:
u ≈ u ~ = Σ i = 1 n N i a i = N a ,
其中,N为形函数,u为位移泛函,i为连接点个数。
泛函的变分为零相当于将泛函对所包含的待定参数进行全微分,并令所得的方程等于零,即:
δ Π = ∂ Π ∂ a 1 δa 1 + ∂ Π ∂ a 2 δa 2 + ... + ∂ Π ∂ a n δa n = 0 ,
由于δa1,δa2,…是任意的,满足上式时必然有…都等于零。因此可以得到一组方程为:
∂ Π ∂ a = ∂ Π ∂ a 1 ∂ Π ∂ a 2 . . . ∂ Π ∂ a n = 0 ,
这是与待定参数a的个数相等的方程组,用以求解a;对于二次泛函,上式退化为一组线性方程为:
∂ Π ∂ a ≡ K a - P = 0 ,
其中,P为所述预定载荷,K为各所述连接点的刚度矩阵。
再根据如下公式得到各连接点的载荷Π:
Π = 1 2 u T K u ,
其中,uT为u的转置矩阵。
对于基于二次规划的传力计算方法,有以下原理:
分析表明,悬挂机构的变形能为正定二次型,由等式约束二次规划的一般形式:
min f ( x ) = 1 2 x T H x + c T x ,
s.t.Ax=b,
H∈Rn×m为对称阵。特别当H正定时,目标函数为凸函数,线性约束下可行域又是凸集,二次规划为凸二次规划,此时,二次规划具有如下良好的性质:
K-T条件不仅是最优解的必要条件而且是充分条件;局部最优解就是全局最优解。
基于以上原理及分析(对于程序接口、方程组和二次规划的解法和策略等不再赘述),针对该传力计算问题,对于增生装置-悬挂机构子系统传力计算,基于变分法与二次规划原理实现复杂翼面上气动分布载荷到挂点集中载荷的计算。通过该方法,可以计算得到增升装置各悬挂机构(如:主襟翼各挂点、缝翼各滑轨挂点)载荷,每个运动机构各挂点(连接点)载荷合成即得到其总载荷。
步骤S104、根据所述悬挂机构与所述多个运动机构的力学平衡方程,得到所述悬挂机构与每一个所述运动机构的接触点载荷。
悬挂机构-非线性环节(运动机构传力问题)子系统传力计算,在已知总载情况下,通过简化分析、隔离体分析,建立滑轨或滑轮架平衡方程即可完成非线性环节的传力计算。构建可考虑变摩擦情况的高升力系统运动机构全参数方程组:
F q · s i n α + μ 1 F q · s i n ( π 2 + α ) + F h · s i n β + μ 2 F h · s i n ( π 2 + β ) + F t · s i n θ - μ c o s θ · F t + F Y = 0.0 F q · c o s α + μ 1 F q · c o s ( π 2 + α ) + F h · c o s β + μ 2 F h · c o s ( π 2 + β ) + F t c o s θ - μ s i n θ · F t + F Z = 0.0 R h g × ( 0.07 2.55 F q + 0.07 1.80 F h ) + F t · L - μ · F t · L 1 + M x = 0.0
其中,FY、fz、mx为所述悬挂机构的总载荷在X、Y、Z三个方向上的分载荷,α是Fq的方向角,β为Fh的方向角,θ为Ft的方向角,Rhg为滑轨外径,L为Ft到滑轨轴线距离,L1为齿轮与齿条间摩擦力到滑轨轴线距离,μ为齿轮摩擦系统。通过以上计算即可完成:悬挂机构传力、运动机构各接触体载荷、高升力系统载荷计算。
实施实例
某飞机高升力系统设计载荷计算。
在采用商用有限元软件MSC/nastran进行大运飞机高升力系统设计载荷计算时,出现了计算速度慢、结果奇异等问题,不能满足设计要求。这些问题是由于nastran软件的非线性模块的缺陷造成的,其他的非线性商用有限元软件,无不需要通过有限元建模、载荷分配、非线性分析等过程来完成。对于飞机设计多载荷情况、复杂系统构型、多缝道、多舵面结构等来说,这些过程太过复杂冗长。
经过深入研究,本专利给出了基于变分法与二次规划的高升力系统载荷计算方法,计算速度快,可以完全避免有限元软件所需的复杂而又冗长的过程,完成了大运飞机的高升力系统设计载荷计算。
飞行试验与地面验证试验表明:本专利给出的计算方法准确可靠。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (3)

1.一种飞机高升力系统载荷计算方法,所述高升力系统中的增升装置与悬挂机构多连接点连接,所述悬挂机构与所述高升力系统中的多个运动机构点接触,其特征在于,包括如下步骤:
步骤一、将所述高升力系统的传力路线分解为增升装置-悬挂机构、悬挂机构-运动机构;
步骤二、为所述高升力系统施加预定载荷,使所述增升装置受力,并将该力通过所述多连接点传递至所述悬挂机构;
步骤三、获取所述多连接点的刚度,根据所述预定载荷和所述多连接点的刚度得到各连接点的载荷,从而得到所述悬挂机构的总载荷;
步骤四、根据所述悬挂机构与所述多个运动机构的力学平衡方程,得到所述悬挂机构与每一个所述运动机构的接触点载荷。
2.根据权利要求1所述的飞机高升力系统载荷计算方法,其特征在于,在所述步骤三中:
首先根据最少势能原理公式(1):
δ Π = ∂ Π ∂ a 1 δa 1 + ∂ Π ∂ a 2 δa 2 + ... + ∂ Π ∂ a n δa n = 0 - - - ( 1 ) ;
得到公式(2):
∂ Π ∂ a ≡ K a - P = 0 - - - ( 2 ) ,
其中,P为所述预定载荷,K为各所述连接点的刚度矩阵;
再取位移泛函(3)为:
u ≈ u ~ = Σ i = 1 n N i a i = N a - - - ( 3 ) ,
其中,N为形函数,u为位移泛函,i为连接点个数;
再根据公式(4)得到各连接点的载荷Π:
Π = 1 2 u T K u - - - ( 4 ) ;
其中,uT为u的转置矩阵。
3.根据权利要求1或2所述的飞机高升力系统载荷计算方法,其特征在于,在所述步骤四中,所述悬挂机构与所述运动机构的接触点为三个,分别为滑轨上前滚轮载荷Fq、为滑轨上后前滚轮载荷Fh、啮合齿廓之间的正压力(缝翼)Ft,所述悬挂机构与所述运动机构的力学平衡方程为:
F q · s i n α + μ 1 F q · s i n ( π 2 + α ) + F h · s i n β + μ 2 F h · s i n ( π 2 + β ) + F t · s i n θ - μ c o s θ · F t + F Y = 0.0 F q · c o s α + μ 1 F q · c o s ( π 2 + α ) + F h · c o s β + μ 2 F h · c o s ( π 2 + β ) + F t c o s θ - μ s i n θ · F t + F Z = 0.0 R h g × ( 0.07 2.55 F q + 0.07 1.80 F h ) + F t · L - μ · F t · L 1 + M x = 0.0
其中,FY、fz、mx为所述悬挂机构的总载荷在X、Y、Z三个方向上的分载荷,α是Fq的方向角,β为Fh的方向角,θ为Ft的方向角,Rhg为滑轨外径,L为Ft到滑轨轴线距离,L1为齿轮与齿条间摩擦力到滑轨轴线距离,μ为齿轮摩擦系统。
CN201510367111.0A 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法 Active CN104933259B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510367111.0A CN104933259B (zh) 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510367111.0A CN104933259B (zh) 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法

Publications (2)

Publication Number Publication Date
CN104933259A true CN104933259A (zh) 2015-09-23
CN104933259B CN104933259B (zh) 2018-04-13

Family

ID=54120425

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510367111.0A Active CN104933259B (zh) 2015-06-29 2015-06-29 一种飞机高升力系统载荷计算方法

Country Status (1)

Country Link
CN (1) CN104933259B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777743A (zh) * 2016-12-28 2017-05-31 石家庄铁道大学 液压悬挂结构的有限元建模方法
CN107515109A (zh) * 2017-08-04 2017-12-26 中国航空工业集团公司西安飞机设计研究所 一种试验件总载测量方法
CN108009334A (zh) * 2017-11-22 2018-05-08 中国航空工业集团公司西安飞机设计研究所 一种飞机货运系统吊轨接头载荷计算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1882474A (zh) * 2003-11-12 2006-12-20 德国空中客车公司 飞机高升力系统的驱动系统中的载荷限定方法
US20070051183A1 (en) * 2004-09-16 2007-03-08 Bernd Schievelbusch Device for detecting synchronous errors of high-lift surfaces on aircraft
CN103502095A (zh) * 2011-04-28 2014-01-08 空中客车德国运营有限责任公司 用于飞行器的高升力系统和用于影响飞行器的高升力特性的方法
CN203756912U (zh) * 2013-12-04 2014-08-06 中国飞机强度研究所 一种载荷分解装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1882474A (zh) * 2003-11-12 2006-12-20 德国空中客车公司 飞机高升力系统的驱动系统中的载荷限定方法
US20070051183A1 (en) * 2004-09-16 2007-03-08 Bernd Schievelbusch Device for detecting synchronous errors of high-lift surfaces on aircraft
CN103502095A (zh) * 2011-04-28 2014-01-08 空中客车德国运营有限责任公司 用于飞行器的高升力系统和用于影响飞行器的高升力特性的方法
CN203756912U (zh) * 2013-12-04 2014-08-06 中国飞机强度研究所 一种载荷分解装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘彦生 等: "民机高升力系统操作试验的符合性验证和设计方法研究", 《航空工程进展》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106777743A (zh) * 2016-12-28 2017-05-31 石家庄铁道大学 液压悬挂结构的有限元建模方法
CN107515109A (zh) * 2017-08-04 2017-12-26 中国航空工业集团公司西安飞机设计研究所 一种试验件总载测量方法
CN108009334A (zh) * 2017-11-22 2018-05-08 中国航空工业集团公司西安飞机设计研究所 一种飞机货运系统吊轨接头载荷计算方法
CN108009334B (zh) * 2017-11-22 2021-03-30 中国航空工业集团公司西安飞机设计研究所 一种飞机货运系统吊轨接头载荷计算方法

Also Published As

Publication number Publication date
CN104933259B (zh) 2018-04-13

Similar Documents

Publication Publication Date Title
Ameduri et al. Morphing wings review: Aims, challenges, and current open issues of a technology
Nicassio et al. Low energy actuation technique of bistable composites for aircraft morphing
Ntouanoglou et al. 4D printing prospects for the aerospace industry: a critical review
Levy et al. Summary of data from the fifth AIAA CFD drag prediction workshop
CN103077259A (zh) 高超声速导弹多场耦合动力学一体化仿真分析方法
Reimer et al. CFD-based gust load analysis for a free-flying flexible passenger aircraft in comparison to a DLM-based approach
CN113723027A (zh) 一种针对弹性飞机的静气动弹性计算方法
CN104933259A (zh) 一种飞机高升力系统载荷计算方法
Wang et al. Optimum component reassignment for balanced systems with multi-state components operating in a shock environment
Pecora et al. KRISTINA: Kinematic rib-based structural system for innovative adaptive trailing edge
Kamal et al. A review of aerodynamics influence on various car model geometry through CFD techniques
Yin et al. Investigation of gear walk suppression while maintaining braking performance in a main landing gear
He et al. Aerodynamic optimisation of a high-speed train head shape using an advanced hybrid surrogate-based nonlinear model representation method
Majid et al. Status and challenges on design and implementation of camber morphing mechanisms
CN102129501B (zh) 一种钢板弹簧动力学仿真方法
Jonsson et al. Development of flutter constraints for high-fidelity aerostructural optimization
Chang et al. Reliability analysis and experimental verification of landing-gear steering mechanism considering environmental temperature
Yang et al. Static aeroelastic modeling and rapid analysis of wings in transonic flow
Zixuan et al. Predictive reentry guidance with aerodynamic parameter online correction
Morton et al. Static and dynamic aeroelastic simulations using Kestrel-A CREATE aircraft simulation tool
CN107357282A (zh) 一种基于多维隐马尔可夫模型的飞行控制系统评估方法
Morton et al. Rigid, maneuvering, and aeroelastic results for Kestrel-a CREATE simulation tool
Peters et al. A Data-Driven Reduced Order Model for Rotor Optimization
Delprete et al. Model-Based Design of Aircraft Landing Gear System
Huang et al. An improved adaptive sampling and experiment design method for aerodynamic 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