[发明专利]一种飞行器燃油系统重心位置的数字化实时确定方法

说明书

技术领域

本发明涉及飞行器燃油系统重心确定技术领域，具体涉及一种飞行器燃油系统重心位置的数字化实时确定方法。

背景技术

飞行器飞行过程中，飞行姿态的变化会引起飞行器燃油系统的晃动，燃油晃动改变了燃油的重心，进而引起全机重心分布的变化。快速、实时获取飞行器燃油重心位置的变化，确定全机重心分布对操控飞行器飞行稳定性至关重要。常规采用计算流体动力学求解飞行器燃油晃动的方法，在求解燃油控制方程时需要划分计算网格，而实际燃油晃动时其自由液面极其不规则，甚至会发生破碎、飞溅，引起网格畸变，最终导致计算终止。

发明内容

为了克服上述现有技术的缺点，本发明的目的是提供一种飞行器燃油系统重心位置的数字化实时确定方法，提高了燃油系统重心位置的实时解算效率和精度。

为了达到上述目的，本发明采取的技术方案为：

一种飞行器燃油系统重心位置的数字化实时确定方法，包括以下步骤：

1)根据飞行器燃油系统贮油箱的实际结构尺寸，即总长、总宽、总高构建贮油箱结构的数字化模型，此模型限定了燃油空间运动界限；

2)获取贮油箱内的燃油油量信息和初始时刻燃油液位信息；

3)将贮油箱内的连续性燃油粒子化近似与离散成一系列具有独立质量的相互作用的燃油粒子，这些燃油粒子即为后续的解算对象，燃油粒子所占据的空间即为后续的粒子化燃油数字计算域；

4)施加粒子化燃油数字计算域的边界条件；

5)基于飞行器飞行姿态，确定燃油系统运动参数；

6)实时解算粒子化燃油系统数字计算模型，即燃油粒子控制方程；

7)实时解算燃油系统重心位置坐标(x_t,y_t,z_t)；

8)最终获取燃油系统重心位置坐标(x_t,y_t,z_t)随时间的变化曲线。

所述的步骤3)中连续性燃油粒子化近似与离散，包含以下步骤：

3.1)采用光滑核函数对描述燃油特性的连续函数，即N-S方程进行核近似；

具体为采用下式对描述燃油特性的连续函数f(x)进行核近似：

其中，W为任一光滑核函数，h为决定光滑核函数支持域尺寸的光滑长度，x为核函数中心点燃油粒子的空间坐标，x’为燃油粒子x的支持域内其他任一燃油粒子的空间坐标，Ω为燃油粒子x的支持域，即粒子化燃油数字计算域；

3.2)应用粒子近似对核近似方程进行近似估算，方法是采用核函数支持域内所有粒子相关值的叠加求和取代核近似方程中相对应的积分项；

具体是将粒子j处的无穷小体元dx'用粒子j的体积ΔV_j代替，实现方程的粒子化近似，则函数f(x)经过粒子化近似后的表达式为：

其中，N为离散后燃油粒子的总数，

若粒子j的密度为ρ_j，则粒子j的质量m_j＝ρ_j·ΔV_j，

函数f(x)进一步修改为则粒子i处描述燃油特性的连续函数f(x)的粒子化近似表达式为：

所述的步骤4)中施加粒子化燃油数字计算域的边界条件，具体包括燃油粒子空间运动界限边界的施加和燃油自由液面边界的判定，具体为：

针对燃油粒子空间运动界限边界的施加，具体方法为：一种是在运动界限处设置一组虚粒子，借助虚粒子对邻近运动界限的真实粒子产生的排斥力来阻止真实粒子穿透运动界限；另一种是在运动界限外部设置镜像粒子，镜像粒子与内部真实粒子关于运动界限对称，镜像粒子与真实粒子速度相反，通过施加压力梯度防止真实粒子穿越运动界限；

针对燃油自由液面边界的判定，具体方法为：一种是因自由液面处的粒子密度是通过周围粒子密度的加权平均确定的，若某一粒子的密度小于实际粒子密度，则认定该粒子位于自由液面，并将该粒子密度强制等于实际粒子密度；另一种方法是若某一粒子支持域内的粒子数量小于内部粒子相同尺度支持域内的粒子数量，则认定该粒子位于自由液面。

所述的步骤7)中实时解算燃油系统重心位置坐标(x_t,y_t,z_t)，任一t时刻燃油系统的重心空间位置的三维坐标(x_t,y_t,z_t)分别由以下公式确定：