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

主要研究俯仰激励下液体大幅晃动问题，将ALE(arbitrary Lagrange-Euler)运动学描述引入到Navier-Stokes方程中，推导了俯仰激励下液体大幅晃动数值模拟计算公式，并利用Galerkin加权余量法推导了有限元数值离散方程，采用ALE有限元方法对方形贮腔中的液体大幅晃动进行了数值模拟计算，对结果进行了比较分析，揭示了俯仰激励下液体大幅晃动问题的非线性现象。     

三维液体非线性晃动及其复杂现象

主要讨论三维液体非线性晃动问题 。将ALE（任意的拉格朗日-欧拉）运动学描述引人到Navier-Stokes方程的分步有限元计算格式中;在时间域上采用分步离散方法中的速度修正格式,利用Galerkin加权余量方法得到了系统的有限元离散方程;推导了考虑表面张力效应时有限元边界条件的弱积分形式;模拟了三维液体的非线性晃动问题,得到了一系列三维液体非线性晃动的复杂现象.进一步模拟了考虑表面张力效应以及在微重力环境下三维液体的非线性晃动,揭示了考虑表面张力效应以及在微重力环境下液体非线性晃动的重要特征.井将所得结论与现有的实验结果进行了比较.从而证实了该方法的有效性与正确性.     

三维液体非线性晃动动力学特性的数值模拟

主要讨论圆筒形贮腔中三维液体非线性晃动问题,将任意的拉格朗日-欧拉（即ArbitraryLagrangian-Eulerian,简称ALE）运动学描述引入到Navies-Stokes方程中,在时间域上采用一种速度和压力的分步计算格式进行时间离散;在空间域上利用Galerkin加权余量法对系统方程进行数值离散;得到了数值计算粘性不可压液体非线性晃动的ALE分步有限元法的计算格式,推导了三维液体自由液面上结点法向矢量的数值计算方法,模拟了圆筒形贮腔（包括带圆环形隔板的圆筒形贮腔）中三维液体的非线性晃动;并得到了一些重要的非线必不知所云 性,通过数值模拟结果与实验结果的比较,证明实了本文方法的可靠性与有效性。     

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

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

微重环境中三维液体非线性晃动的数值模拟

将ALE（任意的拉格朗日-欧拉）运动学描述关系引入到Navier-Stokes方程中,在时间域上采用分步离散方法中的速度修正格式,利用Galerkin加权余量方法推导了系统的有限元数值离散方程;推导了考虑表面张力效应时有限元边界件的弱积分形式。模拟了考虑表面张力情况下圆筒形贮腔中液体的非线性晃动,揭示了考虑表面张力效应时液体非线性晃动的重要特征。     

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

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

液体三维晃动特征问题的有限元数值计算方法

本文采用有限元方法数值求解任意刚性容器内液体三维晃动的固有频率和模态。通过建立液体晃动特征问题的泛函极值原理，编制了四面体等参单元有限元程序，计算了平放圆柱腔内三维液体晃动的特征频率，并将矩形容器、球腔、带“十”字隔板球形容器内的液体三维晃动计算结果与解析解、实验结果和二维有限元数值解进行了比较，程序的正确性得到了验证。     

任意三维贮箱内液体晃动的模态分析及其等效力学模型

本文采用有限元方法建立了数值求解任意三维贮箱内液体晃动固有频率和摸态、动力特征参数的通用程序，利用ANSYS有限元前处理模块．使得程序的通用性达到最大化。通过对液体受迫晃动的结果分析，发现了三维晃动中贮箱外激励与模态响应之间的复杂关系．并且提出了建立液体受迫晃动的等效力学模型的方法。本文以油罐车充液贮箱为例，数值研究了贮箱端部不同所引起的动力学差异，给出了晃动频率和晃动质量的变化曲线。     

低重环境下俯仰运动圆柱贮箱内液体晃动

考虑低重环境下由于表面张力的影响使得圆柱贮箱内液体呈现弯曲自由液面的情况,以俯仰激励下液体晃动的占优模态振型函数为晃动速度势的基函数,利用傅里叶-贝塞尔级数对贮箱受俯仰激励时的自由液面处的运动边界条件进行展开,得到能描述晃动系统本质的广义状态方程,并分别给出了固有频率、晃动波高、晃动力和晃动力矩等晃动特征的计算方法. 通过具体算例得到了俯仰激励下贮箱内液体晃动特征的动态响应,同时验证了文中方法的收敛性、可行性和正确性.      

充液系统液体-多体耦合动力响应分析

提出了充液系统的液体-多体耦合力学模型，基于ALE有限元法和多体系统动力学理论，发展了液体-多体耦合动力响应分析的一种有效方法. 对于液体子系统，将其运动分解为随同贮箱的大位移运动和相对贮箱的大幅晃动，引入贮箱固连参考系中的任意拉格朗日-欧拉(ALE)运动学描述，建立了贮箱固连非惯性参考系中液体的ALE有限元方程，对液体有限元方程的缩聚大大减少了液体子系统的计算规模. 为了计及液体阻尼的影响，引入了液体修正的Rayleigh阻尼，避免了伪阻尼力的出现. 对于多体子系统，应用多体系统动力学理论建立动力学方程. 在此基础上详细导出了液体-多体耦合动力学方程，并采用预估-多重校正算法(PMA)和时间步长控制算法进行迭代求解，既保证了迭代收敛，又提高了计算效率. 所给算例成功求解了液体运输车辆系统的液体-多体耦合动力响应，深入分析了有关参数对系统动力响应的影响，获得了一些结论.     

带环形隔板的圆柱储箱内液体晃动阻尼分析

根据现有的阻尼理论，在线性势流的假设下，分别对带刚性和弹性隔板的圆柱形储液箱内液体晃动做了 有限元仿真计算. 计算结果与实验测定结果相比较发现，在晃动频率上符合得很好，但阻尼 比相差较大. 同时, 通过ALE有限元数值仿真，对储箱内无隔板和有隔板情况下的 黏性流体晃动做了比较，初步推定隔板对液体晃动的``滞阻'的一个重要表现形式为涡耗散.     

Stable response of low-gravity liquid non-linear sloshing in a circle cylindrical tank

Under pitch excitation,the sloshing of liquid in circular cylindrical tank includes planar motion,rotary motion and rotary motion inside planar motion.The boundaries between stable motion and unstable motion depend on the radius of the tank,the liquid height,the gravitational intension,the surface tensor and the sloshing damping.In this article,the differential equations of nonlinear sloshing are built first. And by variational principle,the Lagrange function of liquid pressure is constructed in volume intergration form.Then the velocity potential function is expanded in series by wave height function at the free surface.The nonlinear equations with kinematics and dynamics free surface boundary conditions through variation are derived.At last,these equations are solved by multiple-scales method.The influence of Bond number on the global stable response of nonlinear liquid sloshing in circular cylinder tank is analyzed in detail.The result indicates that variation of amplitude frequency response characteristics of the system with Bond,jump,lag and other nonlinear phenomena of liquid sloshing are investigated.     

基于多尺度方法的平动圆柱贮箱航天器刚——液耦合动力学研究

以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.      

低重环境下俯仰运动圆柱贮箱中液体非线性晃动

在低重力环境下,用变分原理建立了液体晃动的压力体积分形式的Lagrange函数,并将速度势函数在自由液面处作波高函数的级数展开,从而导出自由液面运动学和动力学边界条件非线性方程组;最后用四阶Runge-Kutta法求解非线性方程组。计算结果表明,随俯仰激励频率的逐渐变化,由于面外主模态和次生模态同时失稳,致使整个系统各阶模态和波高函数由稳态运动过渡为不稳定运动。