CN113496097A - 一种基于sph的飞机油箱燃油晃动仿真分析方法 - Google Patents

一种基于sph的飞机油箱燃油晃动仿真分析方法 Download PDF

Info

Publication number
CN113496097A
CN113496097A CN202110484575.5A CN202110484575A CN113496097A CN 113496097 A CN113496097 A CN 113496097A CN 202110484575 A CN202110484575 A CN 202110484575A CN 113496097 A CN113496097 A CN 113496097A
Authority
CN
China
Prior art keywords
fuel
fuel tank
sph
airplane
aircraft
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
CN202110484575.5A
Other languages
English (en)
Other versions
CN113496097B (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 Jiaotong University
Original Assignee
Xian Jiaotong University
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 Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202110484575.5A priority Critical patent/CN113496097B/zh
Publication of CN113496097A publication Critical patent/CN113496097A/zh
Application granted granted Critical
Publication of CN113496097B publication Critical patent/CN113496097B/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/20Design optimisation, verification or simulation
    • G06F30/25Design optimisation, verification or simulation using particle-based methods
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60Testing or inspecting aircraft components or systems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Of Engines (AREA)

Abstract

本发明公开了一种基于SPH的飞机油箱燃油晃动仿真分析方法,旨在解决现有仿真分析方法无法准确仿真飞机在大过载和反飞工况下的飞机油箱燃油晃动的问题。为此目的,本发明中的方法包括:基于SPH飞机油箱燃油动态仿真;对飞机油箱燃油信息进行动态测量解算;根据动态测量解算结果进行误差分析。本发明计算效率高,可仿真特殊飞行姿态下的飞机油箱燃油晃动,并能够精确的动态结算飞机油箱燃油重心。

Description

一种基于SPH的飞机油箱燃油晃动仿真分析方法
技术领域
本发明属于飞机燃油系统领域,具体涉及一种基于SPH的飞机油箱燃油晃动仿真分析方法。
背景技术
燃油系统是飞机整套系统中必不可少的组成部分,其主要任务是贮存飞机所用燃油并为发动机持续输送燃料,以确保发动机在任何飞行状态下都能正常运转,因此燃油系统必须具有可靠性、可控性和连续性供油的能力。
目前,国内外通常使用燃油质量特性数据库进行燃油质量和重心的计算,而燃油质量特性数据库只能描述飞机油箱内燃油在一定液面高度与油面角下的油量。在飞机飞行过程中,飞机的飞行姿态和过载不断发生改变,导致燃油箱内的燃油发生晃动,油面角、燃油体积、燃油质量、燃油重心、油面高度等一系列燃油质量特性数据都不断变化,进而使发动机耗油、燃油箱转输油、全机重心等方面的计算产生误差,影响飞机的安全航行。需要一种高效可靠的方法对燃油的晃动进行仿真和分析,从而为设计人员对燃油质量特性进行测量解算的方面提供重要参考。
发明内容
本发明提供的一种基于SPH的飞机油箱燃油晃动仿真分析方法,能够精确仿真飞机在大过载工况下,燃油箱中燃油剧烈的晃动造成飞机重心偏移的场景,为飞机油量测量和重心计算提供数据。根据仿真出来的飞机油箱传感器,实时计算出传感器浸油高度,为测量燃油油量提供数据支持,同时通过仿真燃油的重心变化,揭示燃油箱液面和重心的变化规律。
为实现上述目的,本发明采用如下技术方案来实现的:
一种基于SPH的飞机油箱燃油晃动仿真分析方法,包括以下步骤:
步骤1)基于SPH对飞机油箱燃油进行动态仿真;
步骤2)对飞机油箱燃油信息进行动态测量解算;
步骤3)根据动态测量解算结果进行误差分析。
本发明进一步的改进在于,所述步骤1)基于SPH对飞机油箱燃油晃动进行动态仿真的方法包括:
步骤1.1),根据飞机油箱模型构建仿真平台的飞机油箱模型,将CATIA中的模型通过变换矩阵转换为仿真平台中的油箱模型;
步骤1.2),初始化燃油粒子:获取燃油粒子的质量、速度和位置信息;
步骤1.3),对飞机燃油施加动态载荷和飞行姿态信息;
步骤1.4),构建飞机燃油粒子晃动仿真计算模型;
步骤1.5),通过燃油粒子运动状态迭代分析,实时更新燃油粒子信息。
本发明进一步的改进在于,步骤1.2)中,设置燃油粒子的质量1.0、速度1.0、半径0.02,燃油粒子静止距离0.01,表面张力设置为0.1。
本发明进一步的改进在于,所述步骤1.4)包括:
1.4.1)计算渲染使用CUDA并行计算引擎与OpenGL渲染引擎;
1.4.2)约束方法使用粒子的密度约束和表面张力约束;
1.4.3)边界条件采取虚粒子法、压强修正法、速度修正法。
本发明进一步的改进在于,所述步骤2)对飞机油箱燃油信息进行动态测量解算的方法包括:
步骤2.1),获取飞机电容式燃油测量传感器位置,将传感器的实际空间坐标通过空间矩阵变换为仿真空间中的坐标,在飞机仿真空间中添加传感器;
步骤2.2),获取液位传感器最高燃油粒子的空间坐标和最低燃油粒子的空间坐标,最高处和最低处周围粒子密度大于给定阈值,以此排除燃油粒子飞溅的影响;
步骤2.3),通过转换矩阵将步骤2.2)中的两个空间坐标转换为实际空间的坐标,求出两点间的空间距离,即对应传感器感的浸油高度;
步骤2.4),通过俯仰角、横滚角和对应传感器的浸油高度在燃油质量数据库中插值解算出相应的油量;
步骤2.5),对每一个油箱中所有传感器测量的油量结果求均值,确定最终测试油量;
步骤2.6),计算所有粒子的形心,作为飞机燃油重心进行输出。
本发明进一步的改进在于,所述步骤3)根据动态测量解算结果进行误差分析的方法包括:
步骤3.1),对仿真得到的传感器高度分别进行±W误差计算,执行上述步骤2.4)和步骤 2.5),以此高度确定相应的测量油量,从而实时得到燃油油量的上界和下界;
步骤3.2),对上述得测得的油量进行滤波平滑操作。
本发明至少具有以下有益的技术效果:
(1)采用CUDA的并行计算,运行效率更高,采用openGL引擎渲染,效果更逼真;
(2)基于SPH的飞机燃油仿真,能够精确仿真飞机在大过载工况下,精确仿真出飞机油箱燃油动态变化的过程;
(3)根据仿真出来的传感器,实时计算出传感器的浸油高度,以此来求得油箱油量,同时能够仿真出燃油重心变化;
(4)不同于商业软件,基于此方法的软件可以做到自主可控。
附图说明
图1:飞机左翼油箱在仿真空间中的三角网格渲染模型图;
图2:(a)是飞机左翼油箱粒子态模型,(b)是燃油液态模型图;
图3:飞机飞行剖面图;
图4:仿真空间中的传感器渲染图;
图5:传感器的感应浸油高度示意图,传感器有的是斜置的;
图6:三个方向上重心变化表数据;其中图6(a)是x坐标上的重心变化图,图6(b)是y坐标上的重心变化图,图6(c)是z坐标上的重心变化图;
图7:平滑过的测试数据评估对比表数据,其中图7(a)是实测油量误差图,图7(b)是实测油量和理论油量对照图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明基于SPH对飞机油箱燃油晃动进行动态仿真,具体实施方式主要分为三个部分:
一、基于SPH对飞机油箱燃油晃动进行动态仿真:
1)首先根据飞机油箱模型构建仿真平台的飞机油箱模型,
CATIA中的模型通过变换矩阵转换为仿真平台中的油箱模型。先将CATIA中模型的坐标(px,py,pz)通过
Figure RE-GDA0003255237680000041
矩阵运算,得到点:
Figure RE-GDA0003255237680000042
其中:lx是CATIA模型中所有点的x坐标最小值,ly是CATIA模型中所有点的y坐标最小值、lz是CATIA模型中所有点的z坐标最小值,同理,得到空间中最小坐标点v0(lx,ly,lz) 与空间中最大坐标点v1(ux,uy,uz))。
通过归一化矩阵
Figure RE-GDA0003255237680000051
将导入的油箱模型进行缩放(s是缩放倍数)。最终得到仿真空间中的粒子坐标:
Figure RE-GDA0003255237680000052
依据这些坐标点在空间中渲染出飞机油箱的三角网格,如图1。
2)初始化燃油粒子:设置燃油粒子的质量1.0、速度1.0、半径0.02,燃油粒子静止距离 0.01(燃油粒子可被挤压变形),表面张力设置为0.1,并导入粒子在仿真空间中的位置,并设置其他燃油粒子参数属性。在本实施例中,a是飞机左翼油箱粒子模型,b是左翼燃油液态渲染模型(见图2);
3)对飞机燃油施加动态载荷和飞行姿态信息:设置x、y、z三个方向上的载荷ax,ay,az,并且设置飞机飞行过程中的俯仰角α和横滚角β,本实施例中,我们给出了飞行过程的剖面图(见图3);
4)构建飞机燃油粒子晃动仿真计算模型:对燃油粒子进行受力分析,分析燃油粒子密度、压强差、和速度差引起的受力,得到总的受力情况,并引入边界条件对燃油粒子进行边界修正。通过受力分析运动状态,实时渲染飞机燃油的变化过程,如图3,是飞机飞行剖面图。
二、对飞机油箱燃油信息进行动态测量解算:
1)根据飞机电容式燃油测量传感器位置,将传感器的实际空间坐标通过空间矩阵变换为仿真空间中的坐标,在飞机仿真空间中添加传感器,如图4;
2)获取液位传感器最高燃油粒子的空间坐标和最低燃油粒子的空间坐标,最高处和最低处周围粒子密度大于给定阈值,以此排除燃油粒子飞溅的影响,如图5,并将相应的坐标转化为CATIA模型中的模型的坐标,转化过程与步骤一中1)的过程相逆:
Figure RE-GDA0003255237680000061
此时,求出这两点的空间距离,即对应传感器的浸油高度;
3)通过俯仰角、横滚角和对应传感器的浸油高度在燃油质量数据库中插值解算出油箱i 的第j根传感器测得的油量vij;对每一个油箱中所有传感器测量的油量结果求均值:
Figure RE-GDA0003255237680000062
Figure RE-GDA0003255237680000063
确定最终总的油量是:
Figure RE-GDA0003255237680000064
4)计算所有粒子的形心(飞机燃油粒子作为质点的特性),作为飞机燃油重心进行输出:
Figure RE-GDA0003255237680000065
如图7所示,三张图表分别显式三个坐标轴方向上的坐标变化。
三、根据动态测量解算结果进行误差分析。
对仿真得到的传感器高度分别进行±L误差计算,L设置为5,如果超过满油时传感器高度,按满油计;小于0值时,计为0。以此高度确定相应的测量油量,从而实时得到燃油油量的上界和下界;此外,对上述得测得的油量进行滤波平滑操作,每十个数据进行一次平滑,结果如图7所示。
以上所揭露的仅为本发明一种基于SPH的飞机油箱燃油晃动仿真分析方法的较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明申请专利范围所作的等同变化,仍属本发明所涵盖范围。

Claims (6)

1.一种基于SPH的飞机油箱燃油晃动仿真分析方法,其特征在于,包括以下步骤:
步骤1)基于SPH对飞机油箱燃油进行动态仿真;
步骤2)对飞机油箱燃油信息进行动态测量解算;
步骤3)根据动态测量解算结果进行误差分析。
2.根据权利要求1所述的一种基于SPH的飞机油箱燃油晃动仿真分析方法,其特征在于,所述步骤1)基于SPH对飞机油箱燃油晃动进行动态仿真的方法包括:
步骤1.1),根据飞机油箱模型构建仿真平台的飞机油箱模型,将CATIA中的模型通过变换矩阵转换为仿真平台中的油箱模型;
步骤1.2),初始化燃油粒子:获取燃油粒子的质量、速度和位置信息;
步骤1.3),对飞机燃油施加动态载荷和飞行姿态信息;
步骤1.4),构建飞机燃油粒子晃动仿真计算模型;
步骤1.5),通过燃油粒子运动状态迭代分析,实时更新燃油粒子信息。
3.根据权利要求2所述的一种基于SPH的飞机油箱燃油晃动仿真分析方法,其特征在于,步骤1.2)中,设置燃油粒子的质量1.0、速度1.0、半径0.02,燃油粒子静止距离0.01,表面张力设置为0.1。
4.根据权利要求2所述的一种基于SPH的飞机油箱燃油晃动仿真分析方法,其特征在于,所述步骤1.4)包括:
1.4.1)计算渲染使用CUDA并行计算引擎与OpenGL渲染引擎;
1.4.2)约束方法使用粒子的密度约束和表面张力约束;
1.4.3)边界条件采取虚粒子法、压强修正法、速度修正法。
5.根据权利要求4所述的一种基于SPH的飞机油箱燃油晃动仿真分析方法,其特征在于,所述步骤2)对飞机油箱燃油信息进行动态测量解算的方法包括:
步骤2.1),获取飞机电容式燃油测量传感器位置,将传感器的实际空间坐标通过空间矩阵变换为仿真空间中的坐标,在飞机仿真空间中添加传感器;
步骤2.2),获取液位传感器最高燃油粒子的空间坐标和最低燃油粒子的空间坐标,最高处和最低处周围粒子密度大于给定阈值,以此排除燃油粒子飞溅的影响;
步骤2.3),通过转换矩阵将步骤2.2)中的两个空间坐标转换为实际空间的坐标,求出两点间的空间距离,即对应传感器感的浸油高度;
步骤2.4),通过俯仰角、横滚角和对应传感器的浸油高度在燃油质量数据库中插值解算出相应的油量;
步骤2.5),对每一个油箱中所有传感器测量的油量结果求均值,确定最终测试油量;
步骤2.6),计算所有粒子的形心,作为飞机燃油重心进行输出。
6.根据权利要求5所述的一种基于SPH的飞机油箱燃油晃动仿真分析方法,其特征在于,所述步骤3)根据动态测量解算结果进行误差分析的方法包括:
步骤3.1),对仿真得到的传感器高度分别进行±W误差计算,执行上述步骤2.4)和步骤2.5),以此高度确定相应的测量油量,从而实时得到燃油油量的上界和下界;
步骤3.2),对上述得测得的油量进行滤波平滑操作。
CN202110484575.5A 2021-04-30 2021-04-30 一种基于sph的飞机油箱燃油晃动仿真分析方法 Active CN113496097B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110484575.5A CN113496097B (zh) 2021-04-30 2021-04-30 一种基于sph的飞机油箱燃油晃动仿真分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110484575.5A CN113496097B (zh) 2021-04-30 2021-04-30 一种基于sph的飞机油箱燃油晃动仿真分析方法

Publications (2)

Publication Number Publication Date
CN113496097A true CN113496097A (zh) 2021-10-12
CN113496097B CN113496097B (zh) 2024-04-02

Family

ID=77998098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110484575.5A Active CN113496097B (zh) 2021-04-30 2021-04-30 一种基于sph的飞机油箱燃油晃动仿真分析方法

Country Status (1)

Country Link
CN (1) CN113496097B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114611209A (zh) * 2022-03-09 2022-06-10 南通大学 一种基于回溯算法的飞行器质心平衡供油策略分析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2191142C1 (ru) * 2001-10-10 2002-10-20 Открытое акционерное общество "Техприбор" Бортовая топливоизмерительная система с компенсацией по характеристическим параметрам топлива
RU2191356C1 (ru) * 2001-10-10 2002-10-20 Открытое акционерное общество "Техприбор" Бортовая топливоизмерительная система с компенсацией по динамической диэлектрической проницаемости топлива
CN105468851A (zh) * 2015-11-26 2016-04-06 中国航空工业集团公司沈阳飞机设计研究所 一种确定飞机动态重量特性的方法
CN111680365A (zh) * 2020-06-16 2020-09-18 郑州大学 一种联合CATIA与Amesim的飞机燃油箱油面变化特性仿真方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2191142C1 (ru) * 2001-10-10 2002-10-20 Открытое акционерное общество "Техприбор" Бортовая топливоизмерительная система с компенсацией по характеристическим параметрам топлива
RU2191356C1 (ru) * 2001-10-10 2002-10-20 Открытое акционерное общество "Техприбор" Бортовая топливоизмерительная система с компенсацией по динамической диэлектрической проницаемости топлива
CN105468851A (zh) * 2015-11-26 2016-04-06 中国航空工业集团公司沈阳飞机设计研究所 一种确定飞机动态重量特性的方法
CN111680365A (zh) * 2020-06-16 2020-09-18 郑州大学 一种联合CATIA与Amesim的飞机燃油箱油面变化特性仿真方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
唐浩等: "导弹油箱燃油晃动仿真分析", 《无线互联科技》, no. 3, pages 1 - 4 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114611209A (zh) * 2022-03-09 2022-06-10 南通大学 一种基于回溯算法的飞行器质心平衡供油策略分析方法
CN114611209B (zh) * 2022-03-09 2023-12-01 南通大学 一种基于回溯算法的飞行器质心平衡供油策略分析方法

Also Published As

Publication number Publication date
CN113496097B (zh) 2024-04-02

Similar Documents

Publication Publication Date Title
CN109444350B (zh) 一种基于无人机的大气污染物监测传感器的布局方法
CN105157590B (zh) 一种基于三维激光扫描技术的建构物健康监测系统
CN112100835B (zh) 一种适用于复杂流动的高效高精度翼型绕流数值模拟方法
CN110188759A (zh) 一种在数字图像相关法中应变场子区动态选择方法
CN110633790B (zh) 基于卷积神经网络的飞机油箱剩余油量测量方法和系统
CN107908913B (zh) 基于并行计算机的地球动力数模方法
CN112414668B (zh) 一种风洞试验数据静气弹修正方法、装置、设备及介质
CN108536912A (zh) 一种输电塔结构力学分析及其App制作的方法
CN102779356B (zh) 曲面网格化系统及方法
CN111339615A (zh) 一种基于两步插值的飞机油箱重心计算与补偿方法
CN113496097B (zh) 一种基于sph的飞机油箱燃油晃动仿真分析方法
CN110162823A (zh) 考虑气动面曲面效应和法向运动的非定常气动力计算方法
CN107204040A (zh) 多点地质统计学建模方法及装置、计算机存储介质
CN109657297A (zh) 一种气动载荷与结构有限元模型双向传递方法
Jude et al. An Octree-based, Cartesian CFD Solver for Helios on CPU and GPU Architectures.
CN109446573A (zh) 一种构建多维路面仿真模型的方法
CN103617291B (zh) 储层成因单元界面等效表征方法
CN102445147A (zh) 影像量测机台测头模拟系统及方法
CN110688805A (zh) 一种汽车外流场流动数值的分析方法
CN108595792B (zh) 一种飞机油箱燃油测量传感器布局优化方法
CN115272594A (zh) 一种基于geotools的等值面生成方法
CN106909730B (zh) 基于同伦映射算法的建筑物三维模型仿真方法及系统
CN115563894A (zh) 一种复杂分层及拼接数字岩心的渗透率预测方法
CN115422649A (zh) 吊挂结构的腐蚀疲劳寿命预测方法及相关设备
CN107644139A (zh) 一种从cad模型到cae模型的属性映射方法

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