CN113761646B - 一种用于移动风场环境中飞行器动响应的确定方法 - Google Patents
一种用于移动风场环境中飞行器动响应的确定方法 Download PDFInfo
- Publication number
- CN113761646B CN113761646B CN202110680212.9A CN202110680212A CN113761646B CN 113761646 B CN113761646 B CN 113761646B CN 202110680212 A CN202110680212 A CN 202110680212A CN 113761646 B CN113761646 B CN 113761646B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- grid
- determining
- aerodynamic
- wind field
- 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
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)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本申请属于飞行器仿真技术领域,特别涉及一种用于移动风场环境中飞行器动响应的确定方法。该方法包括气动载荷计算参数初始化;更新飞行器表面各个网格当地的“等效时刻”;插值得到各个网格所受到的波阵风速度;计算各个网格处的迎角、侧滑角、动压;插值得到各个网格形心处的压力系数;对所有网格进行积分,得到飞行器的外部气动载荷。本发明提供了一种改进的移动风场环境中飞行器外部气动载荷的确定方法,考虑了移动风场包围飞行器过程对气动载荷的影响,简化了气动载荷的确定流程,提高了气动载荷的确定精度。
Description
技术领域
本申请属于飞行器仿真技术领域,特别涉及一种用于移动风场环境 中飞行器动响应的确定方法。
背景技术
动态响应一般是指控制***在典型输入信号的作用下,其输出量从 初始状态到最终状态的响应。
移动风场是一种特殊的突风,与普通风场模型相比有以下3点差异: (1)移动风场具有高速传播的特点,一般以超音速或音速扩张,而普通 风场一般不考虑自身的传播速度;(2)移动风场携带的扰动能量更大, 这表现在它的风场速度较大;(3)移动风场的速度是空间的,可以同时 沿着飞机的航向和法向扰动。当移动风场以音速(或略高于音速)从后方包围飞机时,相对速度较小,尾翼会先受到移动风场的作用,然后才 是机翼受到扰动。
目前,移动风场环境中飞行器动响应没有较好的方法。移动风场环 境中飞行器外部气动载荷的确定方法普遍采用面元法,气动载荷计算技 术及流程复杂,计算精度较低,尤其是在气动非线性区偏差将更大。
发明内容
为了解决上述技术问题,本申请提供了一种用于移动风场环境中飞 行器动响应的确定方法,包括:
步骤S1、对划分为网格的飞行器进行所有网格的初始气动力计算, 确定飞行器外部的气动力及气动力矩;
步骤S2、基于所述气动力及气动力矩进行动响应求解,获得在体轴 系上投影的飞行器各速度分量;
步骤S3、对给定的移动风场数据进行插值,获得飞行器各网格处受 到的与所述体轴系各轴对应的各方向的风场速度;
步骤S4、根据飞行器各速度分量及各方向的风场速度,确定飞行器 各网格处的等效速度分量;
步骤S5、根据所述等效速度分量确定飞行器各网格处的等效气动迎 角及等效侧滑角;
步骤S6、基于所述等效气动迎角及等效侧滑角,对压力分布数据进 行插值,获得飞行器各网格处的压力系数;
步骤S7、基于所述压力系数确定飞行器各网格处的气动力及气动力 矩,返回步骤S1,采用新的气动力及气动力矩替换初始的气动力及气动 力矩,进行下一时刻的飞行器动响应计算,直至仿真结束。
优选的是,步骤S2中,采用四阶龙格-库塔法对飞行器进行动响应求 解。
优选的是,步骤S3进一步包括:
步骤S11、确定仿真时间t;
步骤S12、基于移动风场包围飞行器的速度VΔ确定飞行器各网格在移 动风场中的等效时刻ti:ti=t-Δxi/VΔ,其中,Δxi为第i个网格形心与飞行 器尾翼后边界的水平距离;
步骤S13、根据移动风场数据插值获得各网格处所受到的水平移动风 场速度Ux(ti)、侧向移动风场速度Uy(ti)、垂向移动风场速度Uz(ti)。
优选的是,步骤S4进一步包括:
步骤S41、获取大地坐标系到飞机体轴系的变换矩阵Lbg;
步骤S42、基于所述变换矩阵确定飞行器各网格处的等效速度分量 ui(t)、vi(t)、wi(t):
其中,u(t)、v(t)、w(t)为飞行器各速度分量,Ux(ti)、Uy(ti)、Uz(ti)为各 方向的风场速度。
优选的是,所述变换矩阵根据步骤S2动响应求解过程中的气动特性 数据计算。
优选的是,步骤S6中,所述压力分布数据包括飞行器的原始基础压 力分布数据及升降舵偏转增量压力分布数据,所述飞行器的原始基础压 力分布数据及升降舵偏转增量压力分布数据均通过风洞试验或CFD仿真 获得。
优选的是,在步骤S6中进行插值过程前,进一步包括:
获取飞行器的飞行马赫数;
进行飞行器1g配平,获得飞行器配平迎角及配平升降舵偏度;
获取进行飞行器表面网格划分时的各网格的形心无因次化位置。
优选的是,步骤S7进一步包括:
步骤S71、根据飞行器各网格处的压力系数、等效速压、物面法向矢 量及网格面积确定飞行器各网格处的气动力;
步骤S72、根据飞行器各网格处的气动力及形心位置矢量确定飞行器 各网格处的气动力矩;
步骤S73、根据飞行器各网格处的气动力及气动力矩确定飞行器外部 的气动力及气动力矩;
其中,物面法向矢量、网格面积及形心位置矢量均在网格划分时确 定。
优选的是,所述飞行器表面划分的网格为三角形网格。
本申请考虑移动风场包围飞行器过程对气动载荷的影响,采用风洞 试验及CFD仿真得到的压力分布数据,简化了气动载荷的确定流程,提高 了飞行器气动载荷确定精度及动响应仿真精度。
附图说明
图1是本申请用于移动风场环境中飞行器动响应的确定方法的一优 选实施例的流程图。
图2是本申请图1所示实施例的俯仰角加速度动响应图。
图3是本申请图1所示实施例的飞行器法向载荷系数动响应图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本 申请实施方式中的附图,对本申请实施方式中的技术方案进行更加详细 的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件 或具有相同或类似功能的元件。所描述的实施方式是本申请一部分实施 方式,而不是全部的实施方式。下面通过参考附图描述的实施方式是示 例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本 申请中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下 所获得的所有其他实施方式,都属于本申请保护的范围。下面结合附图 对本申请的实施方式进行详细说明。
本申请提供了一种用于移动风场环境中飞行器动响应的确定方法, 主要包括:步骤S1、对划分为网格的飞行器进行所有网格的初始气动力 计算,确定飞行器外部的气动力及气动力矩;步骤S2、基于所述气动力 及气动力矩进行动响应求解,获得在体轴系上投影的飞行器各速度分量; 步骤S3、对给定的移动风场数据进行插值,获得飞行器各网格处受到的 与所述体轴系各轴对应的各方向的风场速度;步骤S4、根据飞行器各速 度分量及各方向的风场速度,确定飞行器各网格处的等效速度分量;步 骤S5、根据所述等效速度分量确定飞行器各网格处的等效气动迎角及等 效侧滑角;步骤S6、基于所述等效气动迎角及等效侧滑角,对压力分布 数据进行插值,获得飞行器各网格处的压力系数;步骤S7、基于所述压 力系数确定飞行器各网格处的气动力及气动力矩,返回步骤S1,采用新 的气动力及气动力矩替换初始的气动力及气动力矩,进行下一时刻的飞 行器动响应计算,直至仿真结束。
图1给出了一个具体的实施流程,参考图1,本申请的工作原理是: 将整个飞行器按三角网格进行划分,计算得到所有网格的形心位置矢量、 物面法向矢量、面积、形心无因次化位置;按照移动风场时间轴线推进, 更新飞行器表面各个网格当地的“等效时刻”,按照“等效时刻”插值 获得各网格处的移动风场速度,可以详细地考虑移动风场对不同位置网 格的影响,进而考虑网格进出移动风场的动态效应;针对不同网格,分 别计算网格当地的迎角、侧滑角、速压等参数,进而插值获得网格处的 压力系数,这对提高飞行器外部气动载荷的计算精度有积极作用;飞行 器的原始压力分布数据通过风洞试验或CFD仿真获得,数据精度高,计算 结果精度高,这对提高移动风场环境中的飞行器外部气动载荷的计算精度有积极作用;将所有网格的气动载荷累加起来得到飞行器的气动载荷, 代入飞行器动力学方程进行动响应求解。
以下结合算例进行说明。
飞行器表面三角网格的划分可以采用CATIA等商业软件完成,网格 的数量要尽量多、以提高飞行器外部气动载荷的计算精度;通过CATIA 等商业软件将飞行器表面划分为三角形网格,需要导出网格节点坐标、 网格组成节点编号等信息,根据右手法则计算网格的上述5个参量。在一 个具体实施例中,将某飞行器表面划分三角形网格,全机共得到26892个 网格。
1.2、计算参数初始化:输入飞行器飞行高度H、飞行器飞行速度VA(对 应飞行马赫数M)、飞行器质量m、飞行器气动特性数据。例如飞行器 飞行高度H为1500m,飞行器飞行速度VA为200.7m/s(对应飞行马赫数M 为0.6),飞行器质量为80t。
1.3、移动风场包围飞行器速度计算:根据步骤1.2中的飞行器飞行高 度H计算出当前高度上的音速VS、大气密度ρ,移动风场包围飞行器的速 度VΔ计算公式如(1)所示:
VΔ=VS-VA (1)
依据飞行器飞行高度,采用标准大气计算方法得到当前高度上的音 速VS、大气密度ρ;当移动风场从后方包围飞行器时,包围速度需按照上 式进行计算。例如,飞行高度H为1500m处的大气密度ρ为1.0581kg/m3; 音速VS为334.5m/s,根据公式计算移动风场包围飞行器的速度VΔ为 133.8m/s。
1.4、飞行器1g配平计算:采用步骤1.2中的飞行器飞行高度H、飞行 器质量m、飞行器气动特性数据,对飞行器进行1g配平,得到飞行器配平 迎角αtrim、配平升降舵偏度δetrim。本实施例得到飞行器配平迎角αtrim为 1.2°、配平升降舵偏度δetrim为-2.3°。
1.5、计算飞行器全机气动力:对飞行器所有网格的气动力进行计算, 得到飞行器外部的气动力FS(t)、气动力矩MS(t);FS(t)在X轴、Y轴、Z轴 上的投影分量分别为FxS(t)、FyS(t)、FzS(t),MS(t)在X轴、Y轴、Z轴上的 投影分量分别为MxS(t)、MyS(t)、MzS(t)。
1.6、飞行器动响应求解:根据步骤1.5中的FxS(t)、FyS(t)、FzS(t)、MxS(t)、 MyS(t)、MzS(t),采用四阶龙格-库塔法对飞行器进行动响应求解,计算公 式如(2)和(3)所示:
其中,g为重力加速度;u、v和w分别为飞行器速度在体轴系X轴、 Y轴、Z轴上的投影;p、q和r分别为飞行器角速度在体轴系X轴、Y轴、 Z轴上的投影;Ixx、Iyy和Izz分别为飞行器相对体轴系X轴、Y轴、Z轴的 惯性矩;Izx为飞行器的惯性积。
1.7、仿真时间更新:按照移动风场的时间轴线,更新动响应仿真时 间t。
1.8、网格“等效时刻”的更新:任意第i个网格在移动风场中的“等 效时刻”ti计算如公式(4)所示;
ti=t-Δxi/VΔ (4)
其中,Δxi为第i个网格形心与飞行器尾翼后边界的水平距离。
需要说明的是,网格的位置差异引起的时间滞后影响,越是靠近飞 行器头部方向的网格、越晚进入移动风场,相对尾翼上的网格具有时间 滞后。
1.9、网格处的移动风场速度插值:任意第i个网格处所受到的水平移 动风场速度Ux(ti)、侧向移动风场速度Uy(ti)、垂向移动风场速度Uz(ti)根据 移动风场数据插值获得。
需要说明的是,如果网格未进入风场、或者已经退出风场,则移动 风场的速度为零。
1.10、网格当地等效速度更新:任意第i个网格处的等效速度分量ui(t)、 vi(t)、wi(t)计算如公式(5)所示:
其中,Ux(ti)、Uy(ti)、Uz(ti)分别表示移动风场的水平分量、侧向分量、 垂向分量在时刻处ti的风速;Lbg为大地坐标系到飞机体轴系的变换矩阵, Lbg计算如公式(6)所示:
1.11、网格当地等效空速确定:根据步骤1.10中的ui(t)、vi(t)、wi(t), 任意第i个网格处的等效空速Vi计算公式如(7)所示:
1.12、网格当地等效气动迎角αi确定:根据步骤1.10中的ui(t)、vi(t)、 wi(t),根据步骤1.11中的Vi,任意第i个网格处的等效气动迎角αi、等效侧 滑角βi计算公式如(8)所示:
1.13、网格当地等效速压确定:根据步骤1.3中的大气密度ρ、步骤1.11 中的Vi,任意第i个网格处的等效速压Qi计算公式如(9)所示:
1.14、网格处的压力系数确定:根据步骤1.1中的网格形心无因次化 位置步骤1.2中的飞行马赫数M、步骤1.4中的δetrim、步骤1.12 中的等效气动迎角αi、等效侧滑角βi,任意第i个网格处压力系数Cpi(t)通 过对飞行器的原始基础压力分布数据Cp0(通过风洞试验或CFD仿真获得) 和升降舵偏转增量压力分布数据ΔCp0(通过风洞试验或CFD仿真获得) 插值得到,确定公式如(10)所示:
1.15、网格气动载荷的确定:根据步骤1.1中的形心位置矢量ri、物面 法向矢量ni、面积si,根据步骤1.13中的等效速压Qi、步骤1.14中的压力 系数Cpi(t),任意第i个网格的气动力fi和气动力矩mi的计算公式如(11)所 示:
1.16、飞行器外部气动力的确定:根据步骤1.13中的气动力fi和气动 力矩mi,通过对全部网格积分得到飞行器外部的气动力FS(t)、气动力矩 MS(t),具体的计算公式分别如(12)、(13)所示:
其中,N为飞行器外表面的所有网格数目;FS(t)在X轴、Y轴、Z轴上 的投影分量分别为FxS(t)、FyS(t)、FzS(t),MS(t)在X轴、Y轴、Z轴上的投 影分量分别为MxS(t)、MyS(t)、MzS(t)。
1.17,回到步骤1.5,进行迭代计算,直至仿真结束,输出仿真结果。
图2是本发明一个实施例的俯仰角加速度动响应图,横轴为时间,纵 轴为俯仰角加速度,带方框符号的曲线表示考虑移动风场包围飞行器过 程的俯仰角加速度动响应,带三角符号的曲线表示忽略移动风场包围飞 行器过程的俯仰角加速度动响应,虚线曲线表示“1-cos”移动风场激励 信号。
图3是本发明一个实施例的飞行器法向载荷系数动响应图,横轴为时 间,纵轴为飞行器法向载荷系数,带方框符号的曲线表示考虑移动风场 包围飞行器过程的法向载荷系数动响应,带三角符号的曲线表示忽略移 动风场包围飞行器过程的法向载荷系数动响应,虚线曲线表示“1-cos” 移动风场激励信号。
从图2、图3中可以看出:对于考虑移动风场包围飞行器过程的计算 情况,飞行器的俯仰角加速度动响应要更为剧烈一些,且法向载荷系数 出现了明显的双峰值现象。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不 局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内, 可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此, 本申请的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种用于移动风场环境中飞行器动响应的确定方法,其特征在于,包括:
步骤S1、对划分为网格的飞行器进行所有网格的初始气动力计算,确定飞行器外部的气动力及气动力矩;
步骤S2、基于所述气动力及气动力矩进行动响应求解,获得在体轴系上投影的飞行器各速度分量;
步骤S3、对给定的移动风场数据进行插值,获得飞行器各网格处受到的与所述体轴系各轴对应的各方向的风场速度;
步骤S4、根据飞行器各速度分量及各方向的风场速度,确定飞行器各网格处的等效速度分量;
步骤S5、根据所述等效速度分量确定飞行器各网格处的等效气动迎角及等效侧滑角;
步骤S6、基于所述等效气动迎角及等效侧滑角,对压力分布数据进行插值,获得飞行器各网格处的压力系数;
步骤S7、基于所述压力系数确定飞行器各网格处的气动力及气动力矩,返回步骤S1,采用新的气动力及气动力矩替换初始的气动力及气动力矩,进行下一时刻的飞行器动响应计算,直至仿真结束;
其中,步骤S3进一步包括:
步骤S11、确定仿真时间t;
步骤S12、基于移动风场包围飞行器的速度VΔ确定飞行器各网格在移动风场中的等效时刻ti:ti=t-Δxi/VΔ,其中,Δxi为第i个网格形心与飞行器尾翼后边界的水平距离;
步骤S13、根据移动风场数据插值获得各网格处所受到的水平移动风场速度Ux(ti)、侧向移动风场速度Uy(ti)、垂向移动风场速度Uz(ti)。
2.如权利要求1所述的用于移动风场环境中飞行器动响应的确定方法,其特征在于,步骤S2中,采用四阶龙格-库塔法对飞行器进行动响应求解。
4.如权利要求3所述的用于移动风场环境中飞行器动响应的确定方法,其特征在于,所述变换矩阵根据步骤S2动响应求解过程中的气动特性数据计算。
5.如权利要求1所述的用于移动风场环境中飞行器动响应的确定方法,其特征在于,步骤S6中,所述压力分布数据包括飞行器的原始基础压力分布数据及升降舵偏转增量压力分布数据,所述飞行器的原始基础压力分布数据及升降舵偏转增量压力分布数据均通过风洞试验或CFD仿真获得。
6.如权利要求5所述的用于移动风场环境中飞行器动响应的确定方法,其特征在于,在步骤S6中进行插值过程前,进一步包括:
获取飞行器的飞行马赫数;
进行飞行器1g配平,获得飞行器配平迎角及配平升降舵偏度;
获取进行飞行器表面网格划分时的各网格的形心无因次化位置。
7.如权利要求1所述的用于移动风场环境中飞行器动响应的确定方法,其特征在于,步骤S7进一步包括:
步骤S71、根据飞行器各网格处的压力系数、等效速压、物面法向矢量及网格面积确定飞行器各网格处的气动力;
步骤S72、根据飞行器各网格处的气动力及形心位置矢量确定飞行器各网格处的气动力矩;
步骤S73、根据飞行器各网格处的气动力及气动力矩确定飞行器外部的气动力及气动力矩;
其中,物面法向矢量、网格面积及形心位置矢量均在网格划分时确定。
8.如权利要求1所述的用于移动风场环境中飞行器动响应的确定方法,其特征在于,所述飞行器表面划分的网格为三角形网格。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110680212.9A CN113761646B (zh) | 2021-06-18 | 2021-06-18 | 一种用于移动风场环境中飞行器动响应的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110680212.9A CN113761646B (zh) | 2021-06-18 | 2021-06-18 | 一种用于移动风场环境中飞行器动响应的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113761646A CN113761646A (zh) | 2021-12-07 |
CN113761646B true CN113761646B (zh) | 2022-09-06 |
Family
ID=78787465
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110680212.9A Active CN113761646B (zh) | 2021-06-18 | 2021-06-18 | 一种用于移动风场环境中飞行器动响应的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113761646B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115438603B (zh) * | 2022-10-11 | 2023-08-04 | 中国航空工业集团公司西安飞机设计研究所 | 一种移动风场环境中弹性飞机风场动响应的确定方法 |
CN115438602B (zh) * | 2022-10-11 | 2023-08-04 | 中国航空工业集团公司西安飞机设计研究所 | 一种移动风场环境中弹性飞机风场气动载荷的确定方法 |
CN116522429B (zh) * | 2023-07-05 | 2023-09-12 | 西安羚控电子科技有限公司 | 飞行器任务载荷显示特效生成方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102460319A (zh) * | 2009-04-15 | 2012-05-16 | 空中客车运营有限公司 | 用于确定飞行器的局部加速度、动力学载荷分布和空气动力学数据的***和方法 |
CN103577649A (zh) * | 2013-11-13 | 2014-02-12 | 中国航空工业集团公司西安飞机设计研究所 | 运输类飞机货物空投时货舱地板载荷的确定方法 |
FR2995426A1 (fr) * | 2012-09-11 | 2014-03-14 | Airbus Operations Sas | Methode de simulation des charges aerodynamiques instationnaires sur une structure externe d'avion. |
CN105426617A (zh) * | 2015-11-25 | 2016-03-23 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机外挂气动载荷计算方法 |
-
2021
- 2021-06-18 CN CN202110680212.9A patent/CN113761646B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102460319A (zh) * | 2009-04-15 | 2012-05-16 | 空中客车运营有限公司 | 用于确定飞行器的局部加速度、动力学载荷分布和空气动力学数据的***和方法 |
FR2995426A1 (fr) * | 2012-09-11 | 2014-03-14 | Airbus Operations Sas | Methode de simulation des charges aerodynamiques instationnaires sur une structure externe d'avion. |
CN103577649A (zh) * | 2013-11-13 | 2014-02-12 | 中国航空工业集团公司西安飞机设计研究所 | 运输类飞机货物空投时货舱地板载荷的确定方法 |
CN105426617A (zh) * | 2015-11-25 | 2016-03-23 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机外挂气动载荷计算方法 |
Non-Patent Citations (1)
Title |
---|
一种用于获得飞机部件气动载荷的网格向量法;韩鹏等;《空气动力学学报》;20180815(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113761646A (zh) | 2021-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113761646B (zh) | 一种用于移动风场环境中飞行器动响应的确定方法 | |
CN102073755B (zh) | 近空间高超声速飞行器运动控制仿真方法 | |
RU2755843C2 (ru) | Система и способ оценки воздушной скорости летательного аппарата на основании модели накопления данных о погоде. | |
Lu et al. | Real-time simulation system for UAV based on Matlab/Simulink | |
CN113408215B (zh) | 一种用于移动风场环境中飞行器气动载荷的确定方法 | |
CN113848963B (zh) | 一种飞行控制***的控制律参数设计方法 | |
CN108363843A (zh) | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 | |
CN111125935B (zh) | 一种面向临近空间飞行器的仿真*** | |
CN108153330A (zh) | 基于可行域约束的无人飞行器三维航迹自适应跟踪方法 | |
CN113392599B (zh) | 一种弹性飞行器动响应的确定方法 | |
CN106096091B (zh) | 一种飞机运动模拟方法 | |
WO2024140288A1 (zh) | 多旋翼风阻系数标定方法、设备及计算机可读存储介质 | |
CN111508282A (zh) | 低空无人机农田作业飞行障碍物冲突检测方法 | |
CN114020019A (zh) | 飞行器的制导方法与装置 | |
CN104931970B (zh) | 一种面向机载气象雷达仿真的三维云场生成方法 | |
WO2024140280A1 (zh) | 一种多旋翼空速解算方法、设备及计算机可读存储介质 | |
CN113625730B (zh) | 一种基于超扭滑模的四旋翼自适应容错控制方法 | |
CN106570242A (zh) | 低空大动压整体式整流罩高速分离流固耦合仿真方法 | |
EP3734544B1 (en) | Systems and methods for video display | |
CN116560249A (zh) | 一种机动飞行的高自由度简化建模与轨迹仿真方法 | |
CN115686036A (zh) | 一种基于预设性能的变外形飞行器多维复合控制方法 | |
CN115048724A (zh) | 一种基于b型样条的变体空天飞行器气动系数在线辨识方法 | |
CN112762960A (zh) | 一种飞行器所处风场的在线计算方法 | |
CN116894302B (zh) | 一种基于精算迎风面积的leo航天器大气阻力算法 | |
Fruncillo et al. | Effects of Comprehensive Added Masses Modeling on Airship Equations of Motion and Dynamic Stability |
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 |