导弹油箱燃油晃动仿真分析
时间：2020-12-08 20:05:25 来源：雅意学习网 本文已影响 人 
摘 要：随着航空航天技术的发展，研究液体晃动显得越来越重要。液体晃动对飞行器燃油供给、飞行状态以及机体结构变形产生较大影响。如何准确模拟液面变化、计算液体对壁面的冲击力是目前存在的难点。针对这些问题，简要介绍了SPH方法的基本原理，对导弹油箱内燃油的晃动进行了仿真计算，模拟出液面变化。计算了液体重心位置的变化曲线，同时得到油箱与燃油之间的相互作用力，并对结果做出了分析。对导弹油箱内燃油晃动研究提供了一定的理论和工程参考价值。
关键词：光滑粒子动力学；燃油晃动；数值仿真；接触力
贮箱内液体的晃动是指液体自由表面由于受到外加扰动或激励而产生的运动。液体的晃动问题在航空航天、船舶、石油化工以及核动力等领域较为常见，从上世纪50年代起，各国工程师对这一问题已经给予了广泛关注[1]。储箱内带自由液面的液体晃动不仅会对贮箱壁产生冲击载荷，引起结构的疲劳破坏。对于载液比重大的载液系统，如航天器、大型民机以及油罐车等，在运动过程中，短时间内液体的大幅晃动引发的液体重心变化还会对这些载液系统的运动轨迹或飞行姿态产生不利的影响。依据相应的要求，油箱的几何尺寸应避免在飞行中使燃油产生剧烈晃动，使燃油晃动造成的重心变化在可接受范围内；由于飞机姿态变化所带来的燃油重心移动也处于可以接受的范畴[2]。在航空领域，飞机在起飞、着陆及飞行过程中，因气动紊流、发动机转子的不平衡，武器的发射等原因，会引起油箱内燃油的晃动。研究晃动产生的冲击力、晃动力矩所引发的油箱结构强度问题、液体晃动重心的变化对全机稳定性的影响是十分必要的[3]。
液体晃动是一种特殊的流动现象，拥有自由表面流动的相关特点，同时还与约束液体运动的壁面有着较为复杂的相互作用，此类问题的难点就在于对液面大幅晃动的准确模拟以及对液体与壁面相互作用的处理如何准确模拟液面大幅晃动以及处理[4]。因此，对液体晃动问题的求解极具挑战性。
光滑粒子力学（Smoothed Particle Hydrodynamic， SPH）方法广泛用于求解流体动力学问题，可以很好的适应自由边界的变化，并且对于自由表面大变形流动问题的处理效果也比较好。对于液体的边界条件处理，本文考虑了自由液面以及液体与固体壁面两种情况。自由边界问题，文章采用Dirichlet方法处理，根据表面粒子密度来判断特定粒子是否位于自由边界上。而对于液体与固体壁面的相互作用的区域，则在固体壁面上布置一组虚粒子，用来对内部粒子产生排斥力，从而对粒子与固体壁面之间的相互作用进行模拟。
文章基于SPH方法数值模拟导弹油箱受到激励时液体的大幅晃动，模拟了燃油液面变化，计算了燃油晃动过程中重心坐标的变化以及油箱和燃油相互作用力的变化，并对结果作了分析。
1 理论简介
1.1 SPH方法基本原理
SPH的理论基础是粒子方法[5-7]，使用一系列离散点或者粒子来描述连续的系统，系统中的离散的粒子将承载所有的物理信息。
对于粒子i的任意一个函数值都可以应用光滑函数对其支持域内所有粒子相应的函数值进行加权平均后获得，支持域是一个半径为kh的圆，k是确定粒子影响半径的系数，如图1所示，其中S是影响区域的边界。
通过粒子近似式，可以将连续积分形式下的函数转换成在支持域的粒子的离散化求和，通过这种方式实现无网格化。
1.2 核函数
本文研究选择的核函数为Monghan和Lattanzio提出的三次样条函数[8]，也称B样条函数。
为了满足一致性要求，上式中的系数αd在一维、二维和三维空间中分别取为 和 。
1.3 控制方程
拉格朗日描述下的粘性流体控制方程就是Navier-Stokes方程。
用希腊字母上标α、β表示坐标系的方向，可以通过张量理论中的指标法用来表示方程的叠加，并且可以在变化的拉格朗日空间内引入对总时间的导数，将N-S方程作粒子近似，并进行空间离散，得到粒子近似的N-S方程。
连续方程：
动量守恒方程：
式中，σ为总应力张量，由各向同性压力p和粘性应力 τ两部分组成，v为速度，ρ为密度，x为距离，σ为总应力张量。
1.4 边界条件
在液体和固体壁面的相互作用的区域内，液体不可以穿透贮箱壁，同时也不能脱离贮箱壁存在。因此在固壁边界上布置一组虚粒子，对邻近粒子施加排斥力，这样可以防止内部粒子穿透固壁边界，如图2所示。排斥力的方向和两个粒子的中心连线相同，排斥力的大小可以根据公式（5）计算得到[4]。
式中的n1和n2分别取值为12和4；r0是排斥力的作用范围，其值为离散后粒子的初始间距；D为常数，大小与最大速度的平方相近。
处理自由表面边界，则采用Dirichlet边界条件[11]，根据粒子i的粒子密度来判断是否位于自由边界，若粒子 密度满足式（6），那么认定其在自由表面上，并使该粒子密度等于流体表面的实际密度。
式中，ρ0是自由表面流体的实际密度，参数β的大小在0.85-0.98之间，本文将其取值为0.95。
2 模型的建立
2.1 模型说明
贮箱结构取自某型导弹的油箱[9-10]，几何模型如图3所示，左图为贮箱结构图，油箱的长约为830mm，宽约为650mm，高度约为560mm。油箱内部液体为煤油，液体体积约为0.148m3。
为方便将油箱内液体离散成均匀的SPH粒子，采用体积等效的方式，将油箱中的燃油用体积与其相等的形状简单的几何体代替。计算得到油箱内燃油的体积约为0.148m3，为保证体积等效的几何体完全包含于油箱内部，设计如图4所示长方体形状的水柱，燃油划分为5460个SPH单元，油箱划分为壳单元，不考虑油箱的变形，系统采用右手坐标系。
