CN115438602A - 一种移动风场环境中弹性飞机风场气动载荷的确定方法 - Google Patents

一种移动风场环境中弹性飞机风场气动载荷的确定方法 Download PDF

Info

Publication number
CN115438602A
CN115438602A CN202211242691.7A CN202211242691A CN115438602A CN 115438602 A CN115438602 A CN 115438602A CN 202211242691 A CN202211242691 A CN 202211242691A CN 115438602 A CN115438602 A CN 115438602A
Authority
CN
China
Prior art keywords
aerodynamic
determining
aircraft
wind field
elastic
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
CN202211242691.7A
Other languages
English (en)
Other versions
CN115438602B (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 CN202211242691.7A priority Critical patent/CN115438602B/zh
Publication of CN115438602A publication Critical patent/CN115438602A/zh
Application granted granted Critical
Publication of CN115438602B publication Critical patent/CN115438602B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • 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
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请属于弹性飞机动载荷设计技术领域，特别涉及一种移动风场环境中弹性飞机风场气动载荷的确定方法。该方法主要包括步骤S1、基于飞行马赫数、飞机动力学模态数据及非定常气动力模型确定弹性飞机的气动力影响系数矩阵；步骤S2、在任一仿真时刻，确定风速、飞机速压以及弹性飞机气动网格的下洗矩阵；步骤S3、根据气动力影响系数矩阵及弹性飞机气动网格的下洗矩阵，计算频域中的风场气动力系数；步骤S4、基于飞机速压及频域中的风场气动力系数，确定时域中的风场气动力系数，返回步骤S2，更新仿真时刻。本申请能够更好地模拟移动风场从飞机尾部作用过程中的弹性飞机气动载荷变化，满足了工程中的设计需求。

Description

一种移动风场环境中弹性飞机风场气动载荷的确定方法
技术领域
本申请属于弹性飞机动载荷设计技术领域，特别涉及一种移动风场环境中弹性飞机风场气动载荷的确定方法。
背景技术
现代大型飞机结构的柔度不断降低，弹性结构模态振动对飞机的动响应影响有时候是不可忽略的。移动风场与传统风场不一样，移动风场以音速传播、且携带了更强的扰动能量，移动风场有可能激起弹性模态的振动，进而影响飞机的整体运动特性。
目前，已经有关于移动风场环境中刚性飞机气动载荷的确定方法，如申请号为202110680212.9的中国发明专利中，该专利公开了一种用于移动风场环境中飞行器动响应的确定方法。但是对于弹性飞机，普通阵风对飞机的传统非定常激气动力计算方法不能考虑风场的传播速度和水平风速，且不能处理风场从飞机后方追击而来的过程。而移动风场同时具有水平风速、垂直风速的空间特点，且以超音速/音速传播，因此需要建立移动风场环境中的弹性飞机气动载荷的确定方法，并建立移动风场环境中的弹性飞机动响应的确定方法。
发明内容
为了解决上述技术问题至少之一，本申请提供了一种移动风场环境中弹性飞机风场气动载荷的确定方法，可以考虑移动风场包围飞机过程对气动载荷的影响，以满足工程中的设计需求，用于移动风场从飞机尾部追击、作用过程中弹性飞机所受风场气动载荷的确定。
本申请移动风场环境中弹性飞机风场气动载荷的确定方法，主要包括：
步骤S1、基于飞行马赫数、飞机动力学模态数据及非定常气动力模型确定弹性飞机的气动力影响系数矩阵；
步骤S2、在任一仿真时刻，确定风速、飞机速压以及弹性飞机气动网格的下洗矩阵；
步骤S3、根据所述气动力影响系数矩阵及所述弹性飞机气动网格的下洗矩阵，计算频域中的风场气动力系数；
步骤S4、基于所述飞机速压及所述频域中的风场气动力系数，确定时域中的风场气动力系数，返回步骤S2，更新仿真时刻。
优选的是，步骤S1进一步包括：
步骤S11、确定移动风场相对于弹性飞机的初始作用位置；
步骤S12、确定当前海拔高度上的音速、大气密度及飞机的真空速；
步骤S13、计算移动风场从飞机后方包围飞机的包围速度；
步骤S14、对弹性飞机的有限元模型进行动力学求解，获得模态数据；
步骤S15、基于弹性飞机的非定常气动力模型，计算任意第i个气动网格下洗控制点到移动风场初始位置的时间差；
步骤S16、基于气动弹性力学中的亚音速偶极子格网法计算弹性飞机的气动力影响系数矩阵。
优选的是，步骤S2进一步包括：
步骤S21、确定当前的仿真时刻参数；
步骤S22、基于移动风场初始的水平风速及垂向风速，插值出当前的仿真时刻的水平风速与垂向风速；
步骤S23、计算飞机的速压；
步骤S24、对于任意第i个气动网格下洗控制点，确定当前的仿真时刻各气动网格下洗控制点所受到的垂向风速；
步骤S25、确定各气动网格下洗控制点处诱导的下洗；
步骤S26、将下洗公式中的分子分母拆开，并将分母部分与所述速压组合，以修正飞机速压；
步骤S27、将所述分子部分进行Fourier变换；
步骤S28、将所有气动网格在频域上受到的垂向移动速度整理为向量形式，形成弹性飞机气动网格的下洗矩阵。
优选的是，步骤S3进一步包括：
步骤S31、根据所述气动力影响系数矩阵及所述弹性飞机气动网格的下洗矩阵，采用亚音速偶极子格网法计算所有气动网格上的压力系数；
步骤S32、确定频域中的移动风场直接引起的广义气动力系数向量。
优选的是，步骤S32进一步包括：
基于所述飞机动力学模态数据确定所有气动网格下洗控制点处的模态矩阵；
基于所述非定常气动力模型确定所有气动网格的面积矩阵；
基于所述模态矩阵、面积矩阵及所述压力系数确定所述广义气动力系数向量。
优选的是，步骤S4进一步包括：
步骤S41、利用Fourier逆变换将广义气动力系数向量转换为移动风场对飞机的时域广义激励力系数向量；
步骤S42、基于修正后的飞机速压及所述时域广义激励力系数向量，计算移动风场对飞机的时域广义气动载荷向量，即所述时域中的风场气动力系数。
本申请通过引入每个网格的时间滞后项，计算了移动风场从尾部作用于飞机过程中对每个网格的不同下洗影响，是一种考虑了移动风场包围弹性飞机过程对气动载荷影响的计算方法，根据上述方法可以更好地模拟移动风场从飞机尾部作用过程中的弹性飞机气动载荷变化，满足了工程中的设计需求。
附图说明
图1是本申请移动风场环境中弹性飞机风场气动载荷的确定方法的一实施方式的流程图。
图2是本申请一优选实施例的1-cos型移动风场引起的机翼气动弯矩示意图。
图3是本申请一优选实施例的1-cos型移动风场引起的平尾气动弯矩示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚，下面将结合本申请实施方式中的附图，对本申请实施方式中的技术方案进行更加详细的描述。在附图中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施方式是本申请一部分实施方式，而不是全部的实施方式。下面通过参考附图描述的实施方式是示例性的，旨在用于解释本申请，而不能理解为对本申请的限制。基于本申请中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本申请保护的范围。下面结合附图对本申请的实施方式进行详细说明。
本申请移动风场环境中弹性飞机风场气动载荷的确定方法，如图1所示，主要包括：
步骤S1、基于飞行马赫数、飞机动力学模态数据及非定常气动力模型确定弹性飞机的气动力影响系数矩阵；
步骤S2、在任一仿真时刻，确定风速、飞机速压以及弹性飞机气动网格的下洗矩阵；
步骤S3、根据所述气动力影响系数矩阵及所述弹性飞机气动网格的下洗矩阵，计算频域中的风场气动力系数；
步骤S4、基于所述飞机速压及所述频域中的风场气动力系数，确定时域中的风场气动力系数；
步骤S5、基于所述时域中的风场气动力系数，计算飞机动响应。
在一些可选实施方式中，步骤S1主要包括：
S11，给出移动风场的水平风速Ux、垂向风速Uz，二者随时间变化；给出移动风场相对于飞机的初始作用位置x0；
S12，给出飞行海拔高度H、飞行速度Ma，计算当前海拔高度上的音速VS和大气密度ρ、飞机的真空速V；
S13，计算移动风场从飞机后方的包围速度V△，如公式(1)所示：
V△＝VS-V (1)；
S14，对飞机的有限元模型进行动力学求解，获得模态数据；
S15，建立飞机的非定常气动力模型，并计算任意第i个气动网格下洗控制点xi到移动风场初始位置x0的时间差△ti为：
Figure BDA0003885292120000041
S16，根据步骤S12中的飞行马赫数Ma、步骤S14中的模态数据、步骤S15中的非定常气动力模型，采用传统气动弹性力学中的亚音速偶极子格网法计算弹性飞机的气动力影响系数矩阵D。
在一些可选实施方式中，如图2所示，步骤S2主要包括：
S21，对于任一仿真时刻t；
S22，根据步骤S11中的水平风速Ux、垂向风速Uz，插值出时刻t的水平风速Ux(t)、垂向风速Uz(t)；
S23，根据步骤S12中的大气密度ρ、飞机的真空速V、步骤S22中的水平风速Ux(t)、垂向风速Uz(t)，计算飞机的速压
Figure BDA0003885292120000042
为：
Figure BDA0003885292120000051
S24，对于任意第i个气动网格下洗控制点，根据步骤S11中的垂向风速Uz和步骤S15中的时间差△ti，插值出仿真时刻t时该下洗控制点所受到的垂向风速Uzi(t)为：
Uzi(t)＝Uz(t+△ti) (4)；
S25，仿真时刻t时，移动风场在任意第i个气动网格下洗控制点处诱导的下洗wgi(t)计算如公式(5)所示：
Figure BDA0003885292120000052
S26，将公式(5)的分子与分母拆开，并且将分母与公式(3)组合(这里主要是对公式(3)中的速压进行修正，在后续步骤S42中用到)，分别整理为公式(6)、公式(7)：
Figure BDA0003885292120000053
Figure BDA0003885292120000054
S27，将公式(6)进行Fourier变换，如公式(8)所示：
Figure BDA0003885292120000055
其中，i为虚数；ω为振动圆频率；L为参考长度；k为减缩频率，其计算如公式(9)所示：
Figure BDA0003885292120000056
S28，将所有气动网格在频域上受到的垂向移动速度整理为向量形式，如公式(10)所示：
Figure BDA0003885292120000057
其中，m为气动网格的总数目。
在一些可选实施方式中，步骤S3主要包括：
S31，根据步骤S16中的气动力影响系数矩阵D、步骤S28中的垂向移动风速向量
Figure BDA0003885292120000058
采用亚音速偶极子格网法计算所有气动网格上的压力系数△Cp为：
Figure BDA0003885292120000059
S32，根据步骤S14中的模态数据可以计算出所有气动网格下洗控制点处的模态矩阵ΦH，根据步骤S15中的飞机非定常气动力模型可以计算出所有气动网格的面积矩阵S，根据ΦH、S和步骤S31中的压力系数△Cp可以计算频域中的移动风场直接引起的广义气动力系数向量Gw(iω)为：
Figure BDA0003885292120000061
其中，
Figure BDA0003885292120000062
为广义非定常气动力系数矩阵。
在一些可选实施方式中，步骤S4主要包括：
S41，根据步骤S32中的广义气动力系数向量Gw(iω)，利用Fourier逆变换技术，得到移动风场对飞机的时域广义激励力系数向量Gw(t)为：
Figure BDA0003885292120000063
S42，根据步骤S26中修正过的飞机速压
Figure BDA0003885292120000064
步骤S41中的时域广义激励力系数向量Gw(t)，计算移动风场对飞机的时域广义气动载荷向量Fw(t)为：
Figure BDA0003885292120000065
如图2及图3所示，分别给出了尾追1-cos型移动风场引起的机翼气动弯矩和平尾气动弯矩的仿真结果，可以看出，本申请在步骤S15、步骤S25、步骤S28中，通过引入每个网格的时间滞后项，计算了移动风场从尾部作用于飞机过程中对每个网格的不同下洗影响，是一种考虑了移动风场包围弹性飞机过程对气动载荷影响的计算方法，根据上述方法可以更好地模拟移动风场从飞机尾部作用过程中的弹性飞机气动载荷变化，满足了工程中的设计需求。
以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。

Claims (6)

1.一种移动风场环境中弹性飞机风场气动载荷的确定方法，其特征在于，包括：
步骤S1、基于飞行马赫数、飞机动力学模态数据及非定常气动力模型确定弹性飞机的气动力影响系数矩阵；
步骤S2、在任一仿真时刻，确定风速、飞机速压以及弹性飞机气动网格的下洗矩阵；
步骤S3、根据所述气动力影响系数矩阵及所述弹性飞机气动网格的下洗矩阵，计算频域中的风场气动力系数；
步骤S4、基于所述飞机速压及所述频域中的风场气动力系数，确定时域中的风场气动力系数，返回步骤S2，更新仿真时刻。
2.如权利要求1所述的移动风场环境中弹性飞机风场气动载荷的确定方法，其特征在于，步骤S1进一步包括：
步骤S11、确定移动风场相对于弹性飞机的初始作用位置；
步骤S12、确定当前海拔高度上的音速、大气密度及飞机的真空速；
步骤S13、计算移动风场从飞机后方包围飞机的包围速度；
步骤S14、对弹性飞机的有限元模型进行动力学求解，获得模态数据；
步骤S15、基于弹性飞机的非定常气动力模型，计算任意第i个气动网格下洗控制点到移动风场初始位置的时间差；
步骤S16、基于气动弹性力学中的亚音速偶极子格网法计算弹性飞机的气动力影响系数矩阵。
3.如权利要求1所述的移动风场环境中弹性飞机风场气动载荷的确定方法，其特征在于，步骤S2进一步包括：
步骤S21、确定当前的仿真时刻参数；
步骤S22、基于移动风场初始的水平风速及垂向风速，插值出当前的仿真时刻的水平风速与垂向风速；
步骤S23、计算飞机的速压；
步骤S24、对于任意第i个气动网格下洗控制点，确定当前的仿真时刻各气动网格下洗控制点所受到的垂向风速；
步骤S25、确定各气动网格下洗控制点处诱导的下洗；
步骤S26、将下洗公式中的分子分母拆开，并将分母部分与所述速压组合，以修正飞机速压；
步骤S27、将所述分子部分进行Fourier变换；
步骤S28、将所有气动网格在频域上受到的垂向移动速度整理为向量形式，形成弹性飞机气动网格的下洗矩阵。
4.如权利要求1所述的移动风场环境中弹性飞机风场气动载荷的确定方法，其特征在于，步骤S3进一步包括：
步骤S31、根据所述气动力影响系数矩阵及所述弹性飞机气动网格的下洗矩阵，采用亚音速偶极子格网法计算所有气动网格上的压力系数；
步骤S32、确定频域中的移动风场直接引起的广义气动力系数向量。
5.如权利要求4所述的移动风场环境中弹性飞机风场气动载荷的确定方法，其特征在于，步骤S32进一步包括：
基于所述飞机动力学模态数据确定所有气动网格下洗控制点处的模态矩阵；
基于所述非定常气动力模型确定所有气动网格的面积矩阵；
基于所述模态矩阵、面积矩阵及所述压力系数确定所述广义气动力系数向量。
6.如权利要求1所述的移动风场环境中弹性飞机风场气动载荷的确定方法，其特征在于，步骤S4进一步包括：
步骤S41、利用Fourier逆变换将广义气动力系数向量转换为移动风场对飞机的时域广义激励力系数向量；
步骤S42、基于修正后的飞机速压及所述时域广义激励力系数向量，计算移动风场对飞机的时域广义气动载荷向量，即所述时域中的风场气动力系数。
CN202211242691.7A 2022-10-11 2022-10-11 一种移动风场环境中弹性飞机风场气动载荷的确定方法 Active CN115438602B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211242691.7A CN115438602B (zh) 2022-10-11 2022-10-11 一种移动风场环境中弹性飞机风场气动载荷的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211242691.7A CN115438602B (zh) 2022-10-11 2022-10-11 一种移动风场环境中弹性飞机风场气动载荷的确定方法

Publications (2)

Publication Number Publication Date
CN115438602A true CN115438602A (zh) 2022-12-06
CN115438602B CN115438602B (zh) 2023-08-04

Family

ID=84250388

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211242691.7A Active CN115438602B (zh) 2022-10-11 2022-10-11 一种移动风场环境中弹性飞机风场气动载荷的确定方法

Country Status (1)

Country Link
CN (1) CN115438602B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096088A (zh) * 2016-05-31 2016-11-09 中国航空工业集团公司西安飞机设计研究所 一种螺旋桨飞机螺旋颤振分析方法
CN111324991A (zh) * 2019-12-10 2020-06-23 中国飞机强度研究所 一种地面颤振试验中气动力模型的重构方法
CN113761646A (zh) * 2021-06-18 2021-12-07 中国航空工业集团公司西安飞机设计研究所 一种用于移动风场环境中飞行器动响应的确定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096088A (zh) * 2016-05-31 2016-11-09 中国航空工业集团公司西安飞机设计研究所 一种螺旋桨飞机螺旋颤振分析方法
CN111324991A (zh) * 2019-12-10 2020-06-23 中国飞机强度研究所 一种地面颤振试验中气动力模型的重构方法
CN113761646A (zh) * 2021-06-18 2021-12-07 中国航空工业集团公司西安飞机设计研究所 一种用于移动风场环境中飞行器动响应的确定方法

Also Published As

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

Similar Documents

Publication Publication Date Title
CN105629725B (zh) 一种后缘舵滑翔飞行器的弹性运动建模方法
CN113392599B (zh) 一种弹性飞行器动响应的确定方法
CN113868771B (zh) 一种考虑结构和气动非线性的飞行动力学建模方法
CN108595756B (zh) 大包线飞行干扰估计的方法及装置
Dimitriadis et al. A modal frequency-domain generalised force matrix for the unsteady Vortex Lattice method
Pfifer et al. System Identification of a small flexible aircraft-invited
Zhao Flutter suppression of a high aspect-ratio wing with multiple control surfaces
Zeng et al. Experimental model-based aeroelastic control for flutter suppression and gust-load alleviation
Yang et al. Transonic flutter suppression for a three-dimensional elastic wing via active disturbance rejection control
Livne et al. Design, construction, and tests of an aeroelastic wind tunnel model of a variable camber continuous trailing edge flap (VCCTEF) concept wing
Qian et al. Active flutter suppression of a multiple-actuated-wing wind tunnel model
CN112580241A (zh) 一种基于结构降阶模型的非线性气动弹性动响应分析方法
CN111474852B (zh) 一种压电驱动变形机翼的离散滑模控制方法
CN113221237B (zh) 一种基于降阶建模的大迎角颤振分析方法
Mattaboni et al. Active flutter suppression for a three-surface transport aircraft by recurrent neural networks
CN115438602A (zh) 一种移动风场环境中弹性飞机风场气动载荷的确定方法
CN115438603B (zh) 一种移动风场环境中弹性飞机风场动响应的确定方法
Liu et al. Continuous dynamic simulation for morphing wing aeroelasticity
CN115840992A (zh) 一种弹性飞行器飞行仿真方法、系统、计算机存储介质及终端
Haley et al. Generalized predictive control for active flutter suppression
Tong et al. Dynamic response analysis under atmospheric disturbances for helicopters based on elastic blades
CN114707370A (zh) 一种适用于弹性飞机的飞行仿真方法
CN115964795A (zh) 一种基于干扰观测器的变体飞行器变形控制方法
CN108614916A (zh) 一种快速补偿大型轻薄有源相控阵天线阵面变形的方法
Scott et al. An aeroelastic analysis of a thin flexible membrane

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