非线性晃动问题的ALE边界元方法

利用ＡＬＥ（任意的Ｌａｇｒａｎｇｅ－Ｅｕｌｅｒ）边界元方法数值求解了具有自由液面的非线性晃动问题，即受外力激励下流体的非线性振动问题。把ＡＬＥ有限元方法的思想应用到边界元方法中，得到了ＡＬＥ边界元方法。对于自由液面的非线性动力边界条件，应用Ｇａｌｅｒｋｉｎ加权方法进行了有限元数值离散。为了增加求解精度，对动力边界条件提出了增加误差修正项的数值求解方法。对时间变量采用Ｎｅｗｍａｒｋ方法进行离散。推导了系统非线性方程的预测－多次校正法迭代格式。进行了算例分析与比较，得到了令人比较满意的结果。     

分隔贮箱内液体晃动参数的识别方法

分隔贮箱内液体晃动问题可以用一个两自由度非线性阻尼系统予以描述。本文采用等效线性化方法处理非线性阻尼参数，给出了一种利用稳态正弦基础激励晃动试验技术确定分隔贮箱内液体晃动等效力学模型的参数识别方法。方法通过引入晃动模态的频响函数剩余影响的拟合项实现两阶晃动模态的分离识别。数值仿真计算表明，这一方法给出的参数识别结果可以满足工程应用的精度要求     

贮箱内推进剂微重力下晃动的有限元模拟及试验

采用有限体积(VOF)法 任意拉格朗日．欧拉(ALE)法，建立分析液体火箭贮箱内推进剂大幅晃动的数学模型，给出了计算的基本方程。对常重力和低重力下不同形状贮箱内推进剂的非线性晃动进行了数值模拟，并用落塔法进行了微重力下贮箱内液体晃动试验。结果表明，数值模拟法与落体试验的结果相近，在工程上有一定的实用价值。     

非轴对称贮箱液体晃动的等效力学模型

针对非轴对称贮箱液体晃动问题,从一般意义上推导了作三维运动的任意三维贮箱内液体晃动的等效力学模型,采用有限元方法建立计算等效力学模型参数的数值算法,利用矩阵相似变换的性质分离液体零频以消除液体刚度矩阵的奇异性,从而将等效力学模型参数的计算归结为广义特征值问题和线性代数方程组的求解,具有较高的计算效率。通过非轴对称贮箱液体晃动的算例说明,不同的等效晃动方向可能不正交,晃动平面可能不过液体质心,因此非轴对称贮箱液体晃动的等效力学模型比轴对称贮箱液体晃动的更加复杂。最后以三棱柱形充液贮箱为例,采用商用CFD（计算流体动力学）分析软件FLOW-3D验证了在贮箱作典型三维运动时等效力学模型的有效性。
     

低重环境下圆柱贮箱内液体晃动特性研究的一种解析法

考虑低重环境下圆柱贮箱内由于表面张力影响而呈现弯曲自由液面的情况,利用Fourier Bessel级数对贮箱受横向激励时的自由液面处的边界条件进行展开,得到液体晃动系统的广义状态方程,给出了固有频率、晃动波高、晃动力和晃动力矩等晃动特性的计算公式。通过数值算例与文献结果对比,在验证文中方法正确性的同时,具体研究了各晃动特性随充液深度、外激励频率和Bond数等参数变化的规律。     

火箭推进剂液体晃动关机响应的数值仿真

本文采用VOF方法对火箭发动机关机阶段所引起的液体晃动进行了数值仿真，以氧化剂贮箱为例，得到了贮箱无阻尼环和带阻尼环时的晃动仿真图和液体质心变化曲线，检验了VOF方法处理碎波的能力。     

卫星液体晃动气浮台仿真试验系统的运动分析

本文通过边界元数值法求解了部分充液自旋球腔内的液体晃动问题。以流体运动的基本方程和系统运动的Ｅｕｌｅｒ动力学方程为基础，考虑了贮箱偏置、涡旋、重力及Ｃｏｒｉｏｌｉｓ力等因素对流体晃动和系统运动状态的影响，求解出液体的速度场，并在此基础上估算液体的能量耗散率和系统的章动时间常数     

在低微重力下卫星液体晃动力学应用研究

研究充液航天器任意旋转对称偏置贮箱在低微重力情况下的小幅液体晃动问题,通过引入惯性张力强度因子,采用混合加权剩余法,建立低重、微重以及失重情况统一求解的有限元数值模型,并编制具有工程实用性和通用性的计算机分析软件,针对几种贮箱的液体晃动问题进行了计算,应用表明该软件可为复杂航天器的液体晃动设计提供有效的工具。     

用多维模态理论分析航天器贮箱液体有限幅晃动力

首次将多维模态理论应用到求解航天工程中常见的圆柱贮箱液体非线性晃动问题中。针对贮箱作水平横向运动，根据Narirrmnov-Moiseev的三阶渐近假设关系，通过选取主导模态以及确定它们的阶次关系，将一般形式的无穷维模态系统降为5维渐近模态系统。在模态系统基础之上，推导出了液体作用于贮箱壁的力的简洁表达式。数值仿真结果表明，相对于线性晃动理论结果，有限幅晃动液体将对贮箱壁产生更大的横向作用力，同时对贮箱底部还作用有明显的纵向力。可将此公式应用于工程实际问题中部分充液贮箱液体有限幅晃动力的估算。     

圆柱贮箱内液体非线性稳态晃动实验及动力学特性分析

针对圆柱贮箱内液体非线性稳态晃动力和力矩动态信号的采集与分析问题,设计了一种能同时实现充液贮箱整体固定安装和液体动态晃动力、力矩精确测量的实验装置,给出了相应的实验数据采集、分析和处理方法。通过开展大量实验,将实验数据与CFD仿真结果进行细致对比后发现:在旋转晃动状态下,CFD仿真结果存在明显的计算偏差,不便被独立的用于完整揭示贮箱内液体晃动的真实特性。最后,归纳了贮箱内液体的非线性稳态晃动行为随横向外激励频率的变化规律,并对比分析了驻波晃动、拍振晃动、多类复杂的旋转晃动等多种不同形式的非线性稳态晃动行为对运载工具控制系统的任务执行效率和稳定性的影响。     

低重力环境下复杂结构航天器动特性研究的实现途径

文章从试验技术与理论计算结合的角度,对解决复杂结构航天器低重力环境下大幅度液体晃动非线性动态特性问题和模拟月球重力环境下月球探测器动力学性能的研究,进行了建设性的论述。     

微重力下球形贮箱内液体晃动有限元模拟

以卫星球形贮箱为研究对象,用任意拉格朗日-欧拉(ALE)法建立了数学模型,给出了计算的基本方程,并基于MSC-Dytran有限元程序,对常重力和微重力环境中贮箱内推进剂的晃动进行了数值模拟。结果表明,微重力环境中液体更易发生大幅晃动,甚至产生飞溅,工程中常用的等效力学模型已不再适用,而有限元数值模拟方法可较好地解决该问题,为卫星及大充液复杂结构航天器控制系统的设计和分析提供理论参考。     

低重环境下旋转轴对称贮箱内液体晃动研究

针对低重环境下旋转轴对称贮箱自由晃动和受侧向、纵向激励时的箱内液体晃动,用变分原理推导了基于低重环境下静液面形状的液体晃动的积分形式的动力学控制方程,由拉普拉斯方程得到用高斯超几何级数解析表达的速度势和波高,进而得到液体晃动的模态运动方程。分析了球形贮箱和Cassini贮箱内液体自由晃动基频,贮箱受侧向激励时液体晃动波高振幅、晃动对贮箱产生的晃动力和力矩、晃动的等效力学模型,以及贮箱受纵向激励时的晃动波高振幅、可能出现的参激共振问题。     

球形贮箱中三维液体大幅晃动数值模拟

提出一种新颖有效的网格移动算法,使用任意拉格朗日-欧拉（ALE）有限元方法对球形贮箱中三维液体大幅晃动问题进行数值模拟。对接触线网格及自由液面网格指定其运动轨道,对湿壁面网格采用代数更新算法,对域内网格采用Laplace方程对网格速度进行平滑。同时采用基于特征线的时间分裂（CBS）格式来处理对流效应引起的数值振荡。数值计算结果与已有试验数据进行了对比,两者吻合良好,证明了网格移动算法与计算格式的可行性与准确性。     

液体大幅晃动情形的航天器刚液耦合动响应仿真分析

基于任意拉格朗日-欧拉有限元(ALE)方法处理液体大幅晃动问题,并系统地建立了航天器刚-液耦合动力学的数值模型.通过液体受惯性力、刚体受流体作用力与力矩来实现系统的耦合,采用交错算法对液体晃动模块与刚体动力学模块进行迭代求解.将刚-液耦合系统的数值计算结果与物理试验及液体刚化计算结果进行了比较,验证了本文方法的适用性和...     

航天器贮箱大幅液体晃动三维质心面等效模型研究

针对航天器贮箱内的大幅液体晃动建立了三维质心面等效模型,此模型将贮箱内液体等效为位于液体质心处的质点,该质心点可以在质心约束面内任意运动,同时利用质心点动力学方程及接触碰撞模型近似计算液体晃动对贮箱的作用力。质心面等效模型的数值仿真结果与基于VOF方法的流体动力学数值计算结果基本吻合,从而验证了质心面等效模型用于航天器贮箱大幅液体晃动分析的有效性。     

航天器贮箱内液体晃动阻尼研究(一):理论分析

在理想流假设及线性化的自由表面边界条件的基础上,对刚性容器中液体小幅晃动进行了流体相似律的研究。得到了晃动的固有频率及各部分阻尼随特征长度、特征运动粘性系数以及特征重力加速的变化关系。此外还分析了壁面阻尼,受污染的自由液面阻尼以及液体的内部阻尼这几部分晃动阻尼的相对关系,以及它们随各个特征参量的变化情况。分析发现各部分的晃动阻尼比取决于一个无量纲参数,内部阻尼所占的比例随着运动粘性系数的增大、特征尺寸的减小、特征重力加速度的减少而增大,因此在某些情况下内部阻尼不可忽略。     

Equivalent mechanical model for liquid sloshing during draining

For a spacecraft draining liquid fuel during a continuous thruster maneuver, a brief equivalent method is proposed to model the time-varying properties of liquid sloshing for dynamics and control design. The sloshing liquid during draining is equivalent to a set of mechanical model with variable parameters. The model parameters for sample filling ratios are determined by an efficient finite element method according to equivalent principles, while the parameters for other filling ratios are obtained by piecewise linear interpolation. Using the proposed model, forces and torques acted on a Cassini shaped tank by the inside liquid during draining are investigated undergoing several typical motions. Verifications and comparative studies are done with Computational Fluid Dynamics simulations, which confirm the accuracy of the brief model while the sloshing amplitude is small and the flow rate is low.