充液航天器液体晃动和液固耦合动力学的研究与应用

随着火箭运载能力、卫星工作寿命和深空探测器任务复杂度的不断提高, 液体推进剂占航天器总质量的比重也不断增加. 液体推进剂的晃动影响着航天器的运动稳定性和姿轨控系统的可靠性, 是航天器动力学中一个备受关注的问题. 充液航天器中晃动的液体是一个分布参数系统, 理论上是无穷维的, 而工程上希望建立的数学模型是简单、低维的, 因此对液体晃动等效力学模型的研究经久不衰. 另外, 液体推进剂对航天器的结构动特性有着重要的影响, 在建立充液航天器的结构动力学模型时需要考虑液体推进剂与贮箱等结构的耦合效应. 本文首先结合液体晃动动力学理论和航天工程实际, 从理论研究、数值研究和实验研究等三个方面综述了国内外在充液航天器液体晃动动力学领域的研究现状, 并以此为基础介绍了航天工程中液体晃动等效力学模型的应用进展情况; 然后, 以液体运载火箭为例概述了国内外在充液航天器液固耦合建模方面的成果,介绍了求解液固耦合问题的数值方法和应用软件; 最后, 根据航天器工程的发展需求, 对充液航天器液体晃动和液固耦合动力学的进一步研究方向提出了一些建议.      

Damping computation of liquid sloshing in containers aboard spacecraft

Under the non-rotating assumption, a method for the calculation of damping of fuel sloshing with small amplitude in containers aboard spacecraft is proposed in the present paper. And we have presented an eigen-value equation for sloshing damping and frequency computation. This equation may be solved by Ritz or Galerkin methods for a container of simple geometry or by finite element method for a container of arbitrary geometric shape even with rigid baffles. The simulated results show that the equivalent damping coefficients is directly proportional to fuel‘s viscosity, whereas it almost exhibits no influence on sloshing frequencies. The drawback of the proposed method lies in expensive computation cost. Thus far, it hasn‘t yet be applied to a container with elastic baffles.     

储液罐动力学与控制研究进展

阐述了储液罐动力学与控制的工程应用背景, 从3个方面回顾了储液罐动力学与控制的研究进展, 即: 储液罐类液体晃动动力学、液体晃动等效力学模型和储液罐多体系统动力学与控制. 其中对储液罐类液体晃动动力学的研究成果从解析方法和数值方法两方面进行了概述; 在储液罐多体系统动力学中概述了车载、船载储液罐系统动力学和充液航天器固--液--控耦合动力学近年来的研究成果.对今后需要进一步开展的研究方向进行了展望.      

矩形贮箱内液体非线性晃动动力学建模与分析

基于理想流体的假设,根据H-O原理建立了充液贮箱刚体平动与液体非线 性晃动的耦合动力学方程,通过引入改进的势函数描述刚体和液体之间的动边界. 利用伽辽 金方法对动力学方程进行了离散. 针对液体非线性晃动情况,与ALE有限元方法、边界元方 法的结果进行了比较,验证了方法的可行性. 对刚体平动和液体非线性晃动耦合的情 况,数值模拟了多种外力激励下系统的响应. 利用等效力学模型解释了耦合系统固有频率升 高的现象.     

ACTUATION OF SLOSHING MODULATED FORCE AND MOMENT ON LIQUID CONTAINER DRIVEN BY JITTER ACCELERATIONS ASSOCIATED WITH SIEW MOTION IN MICROGRAVITY

ＡＣＴＵＡＴＩＯＮＯＦＳＬＯＳＨＩＮＧＭＯＤＵＬＡＴＥＤＦＯＲＣＥＡＮＤＭＯＭＥＮＴＯＮＬＩＱＵＩＤＣＯＮＴＡＩＮＥＲＤＲＩＶＥＮＢＹＪＩＴＴＥＲＡＣＣＥＬＥＲＡＴＩＯＮＳＡＳＳＯＣＩＡＴＥＤＷＩＴＨＳＩＥＷＭＯＴＩＯＮＩＮＭＩＣＲＯＧＲＡＶＩＴＹＲ...     

颗粒阻尼在航天器减振中的应用与前景

振动抑制是航天器结构设计关心的主要问题之一。本文介绍了一种新型的颗粒阻尼技术,分析了颗粒阻尼器的背景设计难点、减振原理及其影响因素,给出了几个应用实例,对其国内外应用研究现状进行了总结。针对颗粒阻尼器的基本设计原则可以认为：颗粒阻尼器在航天器的蜂窝夹层结构中具有天然的安装优势,而其不受环境限制的特点有望为航天器结构极端工况下的部／组件减振提供一条可行的路径。     

储液容器内液体晃荡的非线性动力学分析

基于非线性波动理论模型,求解储液容器内液体晃动的固有频率、模态及动力学响应问题。流体使用u_s-u_p状态方程,利用ABAQUS软件的自适应网格技术,建立储液容器液体晃动数学模型,通过施加水平简谐激励得到液体晃动的固有频率和模态,并与解析解对比,验证了该方法的准确性与可行性。然后,分析了矩形储液容器在多种激励作用下液体非线性晃动响应特性。

     

圆筒形贮腔中微重力液体非线性晃动的数值模拟

本文讨论低重力液体在圆筒形贮腔中的非线性晃动问题。将ALE（任意的拉格朗日－欧拉）运动学描述关系引入到Navier-Stokes方程中,在时间域上采用分步离散方法中的速度修正格式,利用Galerkinlk加权余量方法推导了系统的有限元数值离散方程;推导了考虑表面张力效应时有限元边界条件的弱积分形式。推导了自由液面上法向矢量的计算公式。模拟了圆筒形贮腔中低重力液体的非线性晃动,并得到了自由液面、波高变化、压力响应等非线性动力特性。揭示了微重力液体非线性晃动的重要特征并将所得结论与现有的实验结果进行了比较。从而证实了本文方法的有效性与正确性。     

Free bending vibration of annular cylindrical tank partially filled with liquid in the consideration of surface wave

ＦＲＥＥＢＥＮＤＩＮＧＶＩＢＲＡＴＩＯＮＯＦＡＮＮＵＬＡＲＣＹＬＩＮＤＲＩＣＡＬＴＡＮＫＰＡＲＴＩＡＬＬＹＦＩＬＬＥＤＷＩＴＨＬＩＱＵＩＤＩＮＴＨＥＣＯＮＳＩＤＥＲＡＴＩＯＮＯＦＳＵＲＦＡＣＥＷＡＶＥＺｈｏｕＤｉｎｇ（周叮）（ＮａｎｊｉｎｇＵｎｉｖｅ...     

俯仰激励下三维液体大幅晃动问题研究

主要研究俯仰激励下三维液体大幅晃动问题,将任意Lagrange-Euler法 (arbitrary Lagrange-Euler, ALE)运动学描述引入到Navier-Stokes方程中,推导了俯仰激励下液体大幅晃 动数值模拟计算公式,并利用Galerkin加权余量法推导了有限元数值离散方程和分步有限 元计算格式,采用ALE分步有限元方法对圆筒形贮腔中的液体大幅晃动进行了数值模拟计 算. 得到了波高、晃动力及晃动力矩等晃动特性的时间变化历程,并对结果进行了分析,揭 示了俯仰激励下三维液体大幅晃动问题的非线性现象.     

基于SPH法的二维矩形液舱晃荡研究

液体晃荡是一种复杂的流体运动现象,自由液面的存在使得该现象具有很强的非线性和随机性。针对二维矩形液舱在不同振幅水平激励下的纵荡问题,应用SPH法对其进行了数值研究。首先计算了小振幅激励下的纵荡,计算结果分别与线性理论解、文献VOF法结果及文献SPH法结果作了比较分析,以验证所建数值模型的合理性;然后计算了液舱在大振幅水平激励下的纵荡,着重分析了不同振幅下液体晃荡的速度向量图、液面波动时程、压强波动时程、动量波动时程以及波动的频谱图,并将计算所得液面波动结果与小振幅激励下的液面波动结果作了比较。分析结果表明,在大振幅水平激励下,液面波动的波峰值较小振幅下的结果有较为明显的增大,而波谷值则无过大的变化,总体波动幅值比小振幅下的结果大;随着激励幅值的增大,液面波动幅值呈现明显增大的趋势,压强的整体波动幅值也呈增大趋势,动量波动的均值亦有明显增大;波动能量随着激励幅值的增大而增大并向第一阶频率区域集中。SPH法对处理液体大幅晃荡这种具有自由表面大变形的问题有十分优越的特性。     

Motion stability dynamics for spacecraft coupled with partially filled liquid container

This paper presents a canonical Hamiltonian model of liquid sloshing for the container coupled with spacecraft. Elliptical shape of rigid body is considered as spacecraft structure. Hamiltonian system is an important form of mechanical system. It mostly used to stabilize the potential shaping of dynamical system. Free surface movement of liquid inside the container is called sloshing. If there is uncontrolled resonance between the motion of tank and liquid-frequency inside the tank then such sloshing can be a reason of attitude disturbance or structural damage of spacecraft. Equivalent mechanical model of simple pendulum or mass attached with spring for sloshing is used by many researchers. Mass attached with spring is used as an equivalent model of sloshing to derive the mathematical equations in terms of Hamiltonian model. Analytical method of Lyapunov function with Casimir energy function is used to find the stability for spacecraft dynamics. Vertical axial rotation is taken as the major axial steady rotation for the moving rigid body.     

Improved method for implementing contact angle condition in simulation of liquid sloshing under microgravity

The finite element simulation of liquid sloshing under microgravity is focused on in this paper. For this class of flows, an important issue is to implement the contact angle boundary condition appropriately. A novel method of adding infinitely small free-surface mesh just adjacent to the contact line is proposed here, which coincides with the physical definition of the contact angle. This free-surface mesh has its orientation determined according to the contact angle, and the mean curvature at the contact line can be computed using this orientation. Hence, surface tension force can be computed and incorporated as an essential boundary condition into the pressure solve. This method is convenient to be applied in three-dimensional simulations, and the use of relatively coarse meshes near the contact line is allowed. In order to validate the proposed method, several numerical simulations of flows in two-dimensional circular and three-dimensional cylindrical and spherical containers are presented, including calculation of equilibrium positions of the free surfaces, unsteady flows of liquid reorientation, and large-amplitude nonlinear liquid sloshing under lateral excitation in microgravity. Parts of the numerical results are compared with the theoretical and published experimental results.     

Experimental findings of the suppression of rotary sloshing on the dynamic response of a liquid storage tank

In liquid storage tanks, rotary sloshing occurs when the frequency of the lateral harmonic load is near the lowest frequency of the tank–liquid system. Rotary sloshing is a type of sloshing that modifies the tank response, which may cause instabilities of the tank wall. However, the consequences of rotary sloshing for the development of strain in the tank wall have not been elucidated. This paper presents an experimental determination of the effects of rotary sloshing on the development of strain and acceleration at various locations of a storage tank. A low-density-polyethylene tank containing water was tested using a shake table. Nine excitations with frequencies near the first free-vibration frequency of the tank–water system were employed. To suppress rotary sloshing, a high-density foam floating lid was utilised as a barrier. Results reveal rotary sloshing boosts not only the development of both hoop and axial strain but also the acceleration in the horizontal direction perpendicular to the excitation. The lid reduced the maximum hoop and axial strain by 500% and 400%, respectively compared to that when rotary sloshing occurs. Moreover, the lid suppressed the nonplanar sloshing by erasing the first three free-vibration frequencies of the tank–water system without the lid.     

流固耦合不可压物质点法及其在晃动问题中的应用

作为一种混合拉格朗日欧拉法,物质点法在流固耦合问题中具有重要的应用前景。对于自由液面的流动问题,基于物质点法框架已建立了弱可压物质点法和完全不可压物质点法,但在处理流固耦合问题时遇到了困难。弱可压物质点法由于采用可压缩状态方程,导致求解时间步长过小,压力振荡严重,产生了非物理的飞溅现象;完全不可压物质点法基于投影算法和不可压条件,消除了弱可压物质点法的压力振荡,提高了时间步长,但难以处理移动边界问题。基于变分形式的投影算法提出了一种新型流固耦合不可压物质点法,得到了体积加权的压力泊松方程PPE(Pressure Poisson Equation),解决了完全不可压物质点法无法处理不规则边界和移动边界的问题。采用流固耦合不可压物质点法研究了运动刚体容器中的液体晃动问题,并与已有实验和数值结果进行对比,验证了算法的正确性和精度。     

Numerical modeling of sloshing motion in a tuned liquid damper outfitted with a submerged slat screen

Tuned liquid dampers (TLDs) are among the most economical and effective passive damping devices. They have been increasingly used to reduce dynamic response and protect structures from failure due to external dynamic excitations. Slat screens are one of the most effective devices used to increase the inherent damping of a TLD, and to reduce the non‐linearity of the free surface motion. A numerical algorithm has been developed to solve the complete non‐linear, moving boundary flow problem in a TLD outfitted with slat screens. The model has been developed to handle conditions leading from small to large interfacial deformations without imposing any linearization assumptions. The numerical algorithm is based on the finite‐difference method. The free surface has been determined using the volume‐of‐fluid method and the donor–acceptor algorithm. The effect of the slat screens has been modeled explicitly by using the partial‐cell treatment method. The present algorithm has been validated against experimental data. The results indicated that the present algorithm is capable of providing accurate details of the flow field inside the TLD and through the screens. These details are essential to improve our understanding of the important parameters governing the performance of a TLD, and hence, to enhance our ability to design better TLDs. Copyright © 2010 John Wiley & Sons, Ltd.     

A boundary element method and spectral analysis model for small‐amplitude viscous fluid sloshing in couple with structural vibrations

A boundary element method (BEM) is presented for the coupled motion analysis of structural vibrations with small‐amplitude fluid sloshing in two‐dimensional space. The linearized Navier–Stokes equations are considered in the frequency domain and transformed into a Laplace equation and a Helmholtz equation with pure imaginary constant. An appropriate fundamental solution for the Helmholtz equation is provided. The conditions of zero stress are imposed on the free surface, and non‐slip conditions of fluid particles are imposed on the walls of the container. For rigid motion models, the expressions for added mass and added damping to the structural motion equations are obtained. Numerical examples are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

结构非线性振动及其控制系统的等效力学模型研究

研究了结构非线性振动及其控制系统的等效力学模型建立问题,对结构进行等效线性化处理;对状态变量的观测,仅测量结构的前几个位移。将随机激励和观测噪声等维化处理后,导出了结构非线性振动及其控制系统的等效力学模型。通过等效力学模型的应用研究,得出了多高层结构在风与地震作用下的多维与一维等效力学模型、等效振型力学模型及其线性力学模型,给出了等效力学模型在多高层结构自适应抗风抗震控制中的应用。     

Active vibration control and suppression for integrated structures

Spacesmictures,aircraftstrUctures,satellitesandsoonarerequiredtobelightinweightduetotherequirementofoperation.TheyarealsolightlydampedbecauseofthelowinternaldampingofthematerialsusedintheirconstrUction.Thus,theywillgeneraiClargeamplitudevibration,whichmayreducetheprecisionofoperationandaffectthePerformanceofoperation.Itisessentialtousesuitablecontrolsystemtocontrolthevibrationofsimctures.Sincethesesmicturesaredistributedparametersystemshavinganinfinitesetofvibrationmodes,thecontrolsystemwith…     

可燃液体爆燃特性及其抑制实验

设计了一套可燃液体爆炸特性实验装置,利用该装置和立式爆轰管对RP-5油料、RP-3油料及工业酒精的爆炸特性、1301惰性气体对这3种燃料的抑制进行了研究。结果表明:RP-5油料、RP-3油料及工业酒精爆炸的体积分数范围分别为1.53%～7.73%、0.82%～7.17%及 3.38%～18.25%;酒精云雾爆轰的临界起爆能为2.11 MJ/m2、爆速和爆压分别为1 609 m/s 、1 480 kPa,爆轰波传播的胞格宽度为14.5 mm,长度为16.2 mm。1301惰性气体对RP-5油料、RP-3油料及工业酒精的最小惰化体积分数分别为6.75%、6.8%及 5.56%;二氧化碳和氮气对RP-3油料的最小惰化体积分数分别为45%和49%;1301惰性气体对油料爆炸抑制效果明显好于二氧化碳与氮气。