CN114218684A - 一种民机载荷谱中动态载荷放大因子确定方法 - Google Patents

一种民机载荷谱中动态载荷放大因子确定方法 Download PDF

Info

Publication number
CN114218684A
CN114218684A CN202111657925.XA CN202111657925A CN114218684A CN 114218684 A CN114218684 A CN 114218684A CN 202111657925 A CN202111657925 A CN 202111657925A CN 114218684 A CN114218684 A CN 114218684A
Authority
CN
China
Prior art keywords
load
damage
dynamic
cycle
flight
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
CN202111657925.XA
Other languages
English (en)
Other versions
CN114218684B (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.)
AVIC First Aircraft Institute
Original Assignee
AVIC First Aircraft Institute
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 AVIC First Aircraft Institute filed Critical AVIC First Aircraft Institute
Priority to CN202111657925.XA priority Critical patent/CN114218684B/zh
Publication of CN114218684A publication Critical patent/CN114218684A/zh
Application granted granted Critical
Publication of CN114218684B publication Critical patent/CN114218684B/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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本申请属于航空全机载荷谱编制技术领域,特别涉及一种民机载荷谱中动态载荷放大因子确定方法。包括:确定需要考虑结构动态响应的飞行使用任务段;获取各个飞行使用任务段在1g状态下的平飞载荷,并将1g状态下的平飞载荷作为基准载荷;基于阵风功率谱密度方法,得到柔性结构在各个飞行使用任务段下的载荷时域历程;根据等损伤折算方法,对载荷循环进行损伤计算,得到载荷循环对应的柔性结构损伤值;根据柔性结构损伤值以及载荷循环中最大载荷循环造成的损伤值,得到结构动态响应造成的损伤比值,即为载荷频次放大的次数,根据等损伤折算方法,得到新生成的载荷循环,从而计算得到动态载荷放大因子。

Description

一种民机载荷谱中动态载荷放大因子确定方法
技术领域
本申请属于航空全机载荷谱编制技术领域,特别涉及一种民机载荷谱中动态载荷放大因子确定方法。
背景技术
民机结构柔性较大,在地面滑行和空中飞行中结构会产生明显的振动,从而引起明显的交变载荷(或载荷相对静态有一定的增大),合理、准确的确定飞机动态载荷放大因子,对于结构耐久性设计、分析具有重要作用。《民机结构耐久性与损伤容限设计手册(上册)》中3.5节“使用载荷及载荷谱的详细说明”的要求,其中,对襟翼放下离场、初期爬升、后期爬升、巡航、初期下降、后期下降和襟翼放下进场7种与突风有关的使用情况,明确规定应该用一个动态系数乘静态载荷或静态应力来考虑动态效应。
飞机在地面着陆撞击、滑行以及空中飞行过程中引起的结构动态响应,对结构的耐久性/损伤容限特性影响比较严重,是结构耐久性载荷设计所关注的一个重要内容,为此在进行结构耐久性/损伤容限分析时,需要考虑结构动态响应的影响。一般情况下,在大飞机结构耐久性/损伤容限分析时,对结构动态响应影响的处理方法有两种,一种为不改变载荷量级,只放大载荷的作用频次;另一种为不改变载荷作用频次,只对载荷幅值进行放大。
部分参考资料对结构的动态载荷放大因子的定义如下:
DMF=Δ载荷动态幅值/Δ载荷静态幅值
动态载荷放大因子(DMF:Dynamic Magnify Factor)就是弹性飞机动响应的最大载荷幅值与刚性飞机响应的最大载荷幅值之比。这里的载荷可以是剪力、弯矩、扭矩、某特定部位的载荷系数等。其中基准载荷是指飞机稳态飞行状态下的载荷,如1g平飞状态。动载荷计算将给出相对于基准载荷的增量载荷时间历程。按照刚性、柔性模型进行动态载荷放大因子计算时,往往得出的动态载荷放大因子较大,在几个机型动态载荷放大因子确定过程中都得到了验证,该种分析方法存在一定的瑕疵,得到的数值难以直接应用到飞机的设计、分析中,且部分机型没有相似机型的经验数值,因此提出一种新的动态载荷放大因子确定方法显得尤为重要。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供了一种民机载荷谱中动态载荷放大因子确定方法,以解决现有技术存在的至少一个问题。
本申请的技术方案是:
一种民机载荷谱中动态载荷放大因子确定方法,包括:
步骤一、通过对飞机典型飞行使用任务剖面进行分析,确定需要考虑结构动态响应的飞行使用任务段;
步骤二、基于飞机气动、惯性载荷计算模型,获取各个飞行使用任务段在1g状态下的平飞载荷,并将1g状态下的平飞载荷作为基准载荷;
步骤三、以基准载荷为基础,在弹性计算模型中,施加了3m/s阵风激励,基于阵风功率谱密度方法,得到柔性结构在各个飞行使用任务段下的载荷时域历程;
步骤四、根据柔性结构的载荷时域历程,对载荷时域历程进行双参数雨流计数,得到动态响应下的载荷循环大小和次数,根据等损伤折算方法,对载荷循环进行损伤计算,得到载荷循环对应的柔性结构损伤值;
步骤五、根据柔性结构损伤值以及载荷循环中最大载荷循环造成的损伤值,得到结构动态响应造成的损伤比值,即为载荷频次放大的次数,根据等损伤折算方法,将最大载荷循环的载荷幅值乘以DMF,确保新的单个载荷循环下与动态响应下的载荷历程造成的损伤相同,得到新生成的载荷循环,从而计算得到动态载荷放大因子。
在本申请的至少一个实施例中,步骤一中,需要考虑结构动态响应的飞行使用任务段包括:襟翼放下离场、初期爬升、后期爬升、巡航、初期下降、后期下降以及襟翼放下进场。
在本申请的至少一个实施例中,步骤四中,所述按照等损伤折算方法,对载荷循环进行损伤计算,得到载荷循环对应的柔性结构损伤值包括:
载荷循环对应的柔性结构损伤值为:
Figure BDA0003448868950000031
Figure BDA0003448868950000032
其中,
Figure BDA0003448868950000033
为第雨流法计及的循环幅值,Ni为每级载荷循环数,m为材料参数。
发明至少存在以下有益技术效果:
本申请的民机载荷谱中动态载荷放大因子确定方法,根据飞机使用任务特点出发,结合飞机气动、惯性载荷计算模型和弹性模型,从而得到动态载荷放大因子,本申请理论依据正确,分析步骤清晰、简单,工程概念明确。
附图说明
图1是本申请一个实施方式的飞机全机坐标系一个角度视图;
图2是本申请一个实施方式的飞机全机坐标系另一个角度视图;
图3是本申请一个实施方式的机翼5肋切面弯矩载荷时域历程示意图;
图4是本申请一个实施方式的机翼5肋切面载荷循环示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
在本申请的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请保护范围的限制。
下面结合附图1至图4对本申请做进一步详细说明。
本申请提供了一种民机载荷谱中动态载荷放大因子确定方法,包括以下步骤:
步骤一、通过对飞机典型飞行使用任务剖面进行分析,确定需要考虑结构动态响应的飞行使用任务段;
步骤二、基于飞机气动、惯性载荷计算模型,获取各个飞行使用任务段在1g状态下的平飞载荷,并将1g状态下的平飞载荷作为基准载荷;
步骤三、以基准载荷为基础,在弹性计算模型中,施加了3m/s阵风激励,基于阵风功率谱密度方法,得到柔性结构在各个飞行使用任务段下的载荷时域历程;
步骤四、根据柔性结构的载荷时域历程,对载荷时域历程进行双参数雨流计数,得到动态响应下的载荷循环大小和次数,根据等损伤折算方法,对载荷循环进行损伤计算,得到载荷循环对应的柔性结构损伤值;
步骤五、根据柔性结构损伤值以及载荷循环中最大载荷循环造成的损伤值,得到结构动态响应造成的损伤比值,即为载荷频次放大的次数,根据等损伤折算方法,将最大载荷循环的载荷幅值乘以DMF,确保新的单个载荷循环下与动态响应下的载荷历程造成的损伤相同,得到新生成的载荷循环,从而计算得到动态载荷放大因子。
在本申请的一个实施方式中,已知某民机巡航任务段的飞行参数。
首先,步骤一、通过对飞机典型飞行使用任务剖面进行分析,确定民机需要考虑动态响应的飞行任务段包括:襟翼放下离场、初期爬升、后期爬升、巡航、初期下降、后期下降以及襟翼放下进场7种;
以机翼结构巡航状态为例进行分析计算;
步骤二中,基于飞机气动、惯性载荷计算模型,得到巡航1g状态下机翼各个肋站位的载荷:
MX5=100000N·m
Fz5=10000N
Tx5=20000N·m
其中,MX5为机翼5肋站位处的垂向弯矩,Fz5为机翼5肋站位处的垂向剪力,Tx5为机翼5肋站位处的航向扭矩。
步骤三中,以基准载荷为基础,在弹性计算模型中,施加了3m/s阵风激励,基于阵风功率谱密度方法,得到机翼5肋站位在巡航任务段下的载荷时域历程曲线,由于机翼设计、分析的关键力素为弯矩,输出机翼弯矩的载荷时域历程,图3给出了动态响应的载荷时域历程,其中3m/s最大阵风响应载荷为:
MX5′=110000N·m
Fz5=11000N
Tx5=21000N·m
步骤四中,(1)依据得到的柔性结构的载荷时域历程,对载荷时域历程进行双参数雨流计数,得到动态响应下的载荷循环大小和次数,其中动态响应载荷的循环有:
C1(最大载荷循环,1g状态-最大载荷状态-1g状态,即MX5=100000N·m—MX5′=110000N·m—MX5=100000N·m);
C2(100500,109500);
C3(101000,109000);
C4(101500,103286.1);
(2)按照等损伤折算方法:根据奥丁aoding公式,对载荷循环进行损伤计算:
Figure BDA0003448868950000051
Figure BDA0003448868950000052
其中,
Figure BDA0003448868950000061
为第雨流法计及的循环幅值,Ni为每级载荷循环数,m为材料参数,对于铝合金取4.0。
得到C1、C2、C3、C4循环对应的损伤值分别为D1、D2、D3、D4(对应的次数分别为N1、N2为1、N3为1、N4为1,一般N1为1),则动态响应过程中的总损伤Dtotal值为:
Dtotal=D1×N1+D2×N2+D3×N3+D4×N4
如果循环数多于4个,依此类推计算。
最后,步骤五中,(1)根据得到的机翼5肋站位的损伤值大小,与载荷循环中最大载荷循环造成的损伤比值,即为结构动态响应造成的附件的损伤比值,该比值即为载荷频次放大的次数:
Nmagnify=Dtotal/(D1×N1)
其中,Nmagnify为机翼5肋站位处垂向弯矩对应的载荷放大频次;
(2)通过等损伤折算方法,将最大载荷循环的载荷幅值乘以DMF,确保新的单个载荷循环下与动态响应下的载荷历程造成的损伤相同:
Dc1′=Dtotal/N1
其中,Dc1′为新生成的载荷循环:1g状态-最大载荷状态-1g状态,即MX5=100000N·m—MX5″=115000N·m—MX5=100000N·m
按照动态放大因子的定义,可以得到该站位的动态载荷放大因子为:
DMFmagnify=Δelasticrigid=(MX5″-MX5)/(MX5′-MX5)=1.5
其中,Δelastic为载荷动态幅值增量,Δrigid为载荷静态幅值增量。
本申请的民机载荷谱中动态载荷放大因子确定方法,根据飞机使用任务特点出发,结合飞机气动、惯性载荷计算模型和弹性模型,从而得到动态载荷放大因子,本申请理论依据正确,分析步骤清晰、简单,工程概念明确。本申请解决了民机载荷谱中动态载荷放大因子精准计算的难题。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。

Claims (3)

1.一种民机载荷谱中动态载荷放大因子确定方法,其特征在于,包括:
步骤一、通过对飞机典型飞行使用任务剖面进行分析,确定需要考虑结构动态响应的飞行使用任务段;
步骤二、基于飞机气动、惯性载荷计算模型,获取各个飞行使用任务段在1g状态下的平飞载荷,并将1g状态下的平飞载荷作为基准载荷;
步骤三、以基准载荷为基础,在弹性计算模型中,施加了3m/s阵风激励,基于阵风功率谱密度方法,得到柔性结构在各个飞行使用任务段下的载荷时域历程;
步骤四、根据柔性结构的载荷时域历程,对载荷时域历程进行双参数雨流计数,得到动态响应下的载荷循环大小和次数,根据等损伤折算方法,对载荷循环进行损伤计算,得到载荷循环对应的柔性结构损伤值;
步骤五、根据柔性结构损伤值以及载荷循环中最大载荷循环造成的损伤值,得到结构动态响应造成的损伤比值,即为载荷频次放大的次数,根据等损伤折算方法,将最大载荷循环的载荷幅值乘以DMF,确保新的单个载荷循环下与动态响应下的载荷历程造成的损伤相同,得到新生成的载荷循环,从而计算得到动态载荷放大因子。
2.根据权利要求1所述的民机载荷谱中动态载荷放大因子确定方法,其特征在于,步骤一中,需要考虑结构动态响应的飞行使用任务段包括:襟翼放下离场、初期爬升、后期爬升、巡航、初期下降、后期下降以及襟翼放下进场。
3.根据权利要求1所述的民机载荷谱中动态载荷放大因子确定方法,其特征在于,步骤四中,所述按照等损伤折算方法,对载荷循环进行损伤计算,得到载荷循环对应的柔性结构损伤值包括:
载荷循环对应的柔性结构损伤值为:
Figure FDA0003448868940000011
Figure FDA0003448868940000021
其中,
Figure FDA0003448868940000022
为第雨流法计及的循环幅值,Ni为每级载荷循环数,m为材料参数。
CN202111657925.XA 2021-12-30 2021-12-30 一种民机载荷谱中动态载荷放大因子确定方法 Active CN114218684B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111657925.XA CN114218684B (zh) 2021-12-30 2021-12-30 一种民机载荷谱中动态载荷放大因子确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111657925.XA CN114218684B (zh) 2021-12-30 2021-12-30 一种民机载荷谱中动态载荷放大因子确定方法

Publications (2)

Publication Number Publication Date
CN114218684A true CN114218684A (zh) 2022-03-22
CN114218684B CN114218684B (zh) 2022-09-20

Family

ID=80707186

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111657925.XA Active CN114218684B (zh) 2021-12-30 2021-12-30 一种民机载荷谱中动态载荷放大因子确定方法

Country Status (1)

Country Link
CN (1) CN114218684B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114878123A (zh) * 2022-07-12 2022-08-09 中国飞机强度研究所 飞机疲劳强度测试中基于加速模型的组合载荷谱编制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110909501A (zh) * 2019-11-20 2020-03-24 中国核动力研究设计院 一种系统动力分析中载荷放大因子的计算方法
CN112528478A (zh) * 2020-12-04 2021-03-19 中国航空工业集团公司成都飞机设计研究所 一种无人机阵风载荷谱快速编制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110909501A (zh) * 2019-11-20 2020-03-24 中国核动力研究设计院 一种系统动力分析中载荷放大因子的计算方法
CN112528478A (zh) * 2020-12-04 2021-03-19 中国航空工业集团公司成都飞机设计研究所 一种无人机阵风载荷谱快速编制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
吕程等: "某型民机后机身疲劳试验若干加速方法的研究与应用", 《工程与试验》 *
宋迎东等: "发动机机动飞行类综合载荷谱研究", 《航空动力学报》 *
李唐等: "基于Fokker F27机群载荷谱损伤分散性计算分析", 《北京航空航天大学学报》 *
王长江等: "使用情况差异对民机载荷谱地空地损伤分散性的影响", 《机械强度》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114878123A (zh) * 2022-07-12 2022-08-09 中国飞机强度研究所 飞机疲劳强度测试中基于加速模型的组合载荷谱编制方法
CN114878123B (zh) * 2022-07-12 2022-09-16 中国飞机强度研究所 飞机疲劳强度测试中基于加速模型的组合载荷谱编制方法

Also Published As

Publication number Publication date
CN114218684B (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
Raveh CFD-based models of aerodynamic gust response
CN106529093A (zh) 一种针对大展弦比机翼的气动/结构/静气弹耦合优化方法
CN112528478B (zh) 一种无人机阵风载荷谱快速编制方法
Murua Flexible aircraft dynamics with a geometrically-nonlinear description of the unsteady aerodynamics
CN111274648B (zh) 一种民用飞机前缘襟翼的分布式飞行载荷设计方法
CN114218684B (zh) 一种民机载荷谱中动态载荷放大因子确定方法
Miller et al. Evaluation of the hinge moment and normal force aerodynamic loads from a seamless adaptive compliant trailing edge flap in flight
Precup et al. The design, construction, and tests of a concept aeroelastic wind tunnel model of a high-lift variable camber continuous trailing edge flap (HL-VCCTEF) wing configuration
CN110455477A (zh) 一种固体火箭舱段结构振动载荷谱的获取方法
Koning et al. Using RotCFD to Predict Isolated XV-15 Rotor Performance
EP2924598A1 (en) A method for determining a structural response of a flow body to an atmospheric disturbance
Abhishek et al. Prediction and fundamental understanding of stall loads in UH-60A pull-up maneuver
Klimmek et al. Aircraft loads-an important task from pre-design to loads flight testing
Benedict et al. Inclusion of aeroelastic twist into the CFD analysis of the twin-engine NASA common research model
Jacobson et al. Carefree hybrid methodology for rotor hover performance analysis
Chen et al. ZONA6 versus the Doublet-Lattice method for unsteady aerodynamics on lifting surfaces
Gov et al. Geometrically nonlinear model for gust response of very flexible wings using segmental modes
CN113449375B (zh) 一种复合材料桨叶疲劳寿命的半解析计算方法
Chen et al. A hybrid Navier Stokes-free wake method for modeling tandem rotors
CN117787120A (zh) 一种用于飞行仿真的突风载荷时域计算方法
CN113536456B (zh) 一种基于动力学响应辨识的大柔性结构降阶模型建模方法
Moreno et al. A rigid horizontal tail wind tunnel test for high transonic mach and high frequency unsteady pressure acquisition
Neves et al. Stall in ground effect using the unsteady vortex lattice method with Kirchhoff-based correction
Camberos et al. An efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment
Delrieux From design to flight testing: overview of rotorcraft acoustic research at onera for industrial applications.

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