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

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

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

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

纵向参数激励下平动刚-液耦合系统稳定性

纵向参数激励下液体晃动稳定性是航天器动力学中一个广受关注的问题,然而在以往的研究中没有考虑液体晃动与航天器运动之间的耦合作用对系统参数振动稳定性的影响. 建立了用液体晃动等效单摆模型描述的纵向激励下平动刚-液耦合系统的Mathieu方程,采用摄动法确定了耦合系统1/2亚谐波振动和谐波振动的激励幅-频稳定性边界. 研究发现,液体晃动与主刚体横向运动的耦合作用扩大了参数振动不稳定区, 并使其向高频移动,影响的程度随等效晃动质量的减小而减小;液体晃动模态阻尼对1/2亚谐波振动不稳定区的缩小作用远弱于对谐波振动不稳定区的缩小作用. 对耦合系统第1阶液体晃动模态1/2亚谐波振动响应的研究表明:当纵向激励参数在不稳定区内时,可能引起主刚体的纵横耦合振动现象.      

燃料消耗下充液航天器等效动力学建模与分析

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.      

EQUIVALENT DYNAMICS MODELING AND ANALYSIS OF LIQUID-FILLED SPACECRAFT WITH FUEL CONSUMPTION 1)

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.     

考虑几何非线性和热效应的刚-柔耦合动力学

温度增高和温度梯度会引起梁的纵向、横向变形位移,在一定程度上对刚-柔耦合规律产生影响.该文考虑热应变,从平面梁的非线性的应变与位移关系式出发,建立了刚体运动、弹性变形和温度相互耦合的有限元离散的热传导方程和动力学方程.研究热流作用下的中心刚体-简支梁系统的刚-柔耦合动力学性质,揭示了几何非线性项和热应变对弹性变形和刚体运动影响.     

微重力下圆柱形贮箱内液体晃动的分岔现象

为了得到微重力下液体晃动的特性,晃动模态用具有水平静液面的液体晃动模态近似. 液体晃动取前5 个重要模态,用Lagrange 方程推导了横向力作用下液体晃动和航天器结构的无量纲的耦合动力学方程并进行数值模拟. 模拟了液体晃动模态随外力振幅和频率变化产生的分岔现象,并分析了系统参数,如Bond数、接触角、接触角迟滞、充液高度、频率比、质量比以及外力的周期形式和方向等对晃动模态的分岔的影响.     

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

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

中心刚体-柔性梁耦合系统离散模型的研究

采用数值仿真对由中心刚体、柔性梁组成的刚-柔耦合系统的动力学离散模型进行了研究.考虑到刚柔-耦合系统的控制方程没有精确解析解,只能寻求数值解,最广泛使用的离散方法是有限元,但其广义坐标数目过于庞大,因此本文探讨了采用经典结构动力学中不同边界的模态函数离散动边界下刚柔耦合动力学方程的可行性及各自的优劣,得到刚柔耦合系统的模态缩减规律.     

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

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

Dynamic response of liquid-filled rectangular tank with elastic appendages under pitching excitation

Nonlinear dynamics of liquid-filled rectangular tank with elastic appendages are studied. Based on the assumption of ideal fluid, the coupling dynamic equations of rigid tank, elastic appendages and liquid fuel are derived using H-O principle. In the case of pitch excitation, the modified potential function and wave height function are introduced to describe the moving boundary of fluid, then Galerkin's method is used to discretize the dynamic equations into ordinary differential equations. The natural frequencies of the coupling system are formulated in liquid depth, the length of the tank, etc. The formulae are confirmed by numerical simulations, which also show that the effects of liquid and elastic appendages on the attitude angular of rigid.     

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.     

A coupled BE–FE study for evaluation of seismically isolated cylindrical liquid storage tanks considering fluid–structure interaction

The effects of base-isolation on the seismic response of cylindrical vertical flexible liquid storage tanks subjected to horizontal seismic ground motion are presented in this paper. The whole system consists of two main parts: the base-isolation component, and the fluid–structure interaction subsystem. Dynamic analysis of liquid storage tank is achieved through the use of finite shell elements for the structure and internal boundary elements for the liquid region. The boundary element equations are employed to obtain an equivalent liquid mass matrix which is then coupled with the shell structure mass matrix, resulting in the coupled equations of motion. Finally, the coupled equations of motion are connected with the base-isolation system to observe the whole system behavior. A bilinear hysteretic element is used to illustrate the base-isolation system. The analysis is performed in the time domain in which, hydrodynamic interaction is taken into account. It is shown that the seismic response of the isolated tank could be significantly reduced compared with the fixed-foundation tanks. Parametric studies are carried out to study the effects of different system parameters on the effectiveness of the base-isolation system. These parameters are the tank geometry aspect ratio (height to radius), the flexibility of the isolation system, the liquid surface displacement variations, and the tank wall flexibility. It is observed that the seismic isolation is more effective in slender tanks in comparison with broad tanks. Furthermore, it is shown that the isolation efficiency is partially more significant in rigid tanks. It is also noticeable that flexible isolators considerably reduce the seismic response in comparison with stiff isolators. Despite the foregoing advantages, liquid surface displacement increases due to seismic isolation, especially in slender tanks.     

非结构动网格在三维可动边界问题中的应用

研究用于非结构动网格的弹簧近似方法，采用顶点弹簧描述，导出并讨论了弹簧倔强系数的取值。通过引入边界修正和扭转效应修正，对标准弹簧近似方法进行了改进，转动翼型算例的结果表明，改进后的方法大大提高了网格变形能力和网格质量，应用该动网格方法耦合求解基于（Arbitrary Lagrangian-Eulerian,ALE）描述的三维Euler方程，模拟了作俯抑振动的矩形机翼绕流，计算结果与实验数据及文献计算结果十分一致，作为多个自由刚与流体耦合运动问题的简单例证，耦合刚体动力学方程，模拟了激波与双立方体的相互作用，得到了非定常流场结构，研究表明，基于弹簧近似的非结构动网格与有限体积流式流场解算器相结合，是模拟包含运动边界的非定常流动问题的有效方法。     

Nonlinear multimodal model for TLD of irregular tank geometry and small fluid depth

Tuned liquid dampers (TLDs) utilize sloshing fluid to absorb and dissipate structural vibrational energy. TLDs of irregular or complex tank geometry may be required in practice to avoid tank interference with fixed structural or mechanical components. The literature offers few analytical models to predict the response of this type of TLD, particularly when the fluid depth is small. In this paper, a multimodal model is developed utilizing a Boussinesq-type modal theory which is valid for small TLD fluid depths. The Bateman–Luke variational principle is employed to develop a system of coupled nonlinear ordinary differential equations which describe the fluid response when the tank is subjected to base excitation. Energy dissipation is incorporated into the model from the inclusion of damping screens. The fluid model is used to describe the response of a 2D structure–TLD system when the structure is subjected to external loading and the TLD tank geometry is irregular.Shake table experiments are conducted on a rectangular and chamfered tank subjected to unidirectional base excitation. Comparisons of the experimental and predicted sloshing forces and energy dissipation per cycle indicate that the model is able to predict the fluid response at fluid depth ratios greater than h/L=0.10. Next, structure–TLD system tests are conducted and it is found that the model can predict the structural and TLD responses. The simulated and experimental results show that the TLD tank transfers energy between orthogonal structural sway modes.     

Quasi-Euler-Poinsot motion of a rigid body

The phase-plane method of nonlinear oscillation is used to discuss the influence of the small dissipation upon the Euler-Poinsot motion of a rigid body about a fixed point. The equations of phase coordinates are applied instead of Eulerian equations, and the global characteristics of the motion of rigid body are analysed according to the distribution and the type of the singular points. A Chaplygin's sphere on a rough plane, a rigid body in viscous medium and one with a cavity filled with viscous fluid are discussed as examples. It is shown that the motions of rigid bodies dissipated by various physical factors have a common qualitative character. The rigid body tends to make a permanent rotation about the principal axis of the largest moment of inertia. The transitive process can change from oscillatory to aperiodic with the decrease in dissipation.     

充液矩形贮箱弹性箱底板应力与变形分析

基于相容拉格朗日-欧拉法,通过对流场与弹性固体间流固耦合作用的分析,得到了矩形贮箱弹性底板流固耦合系统的自由振动微分方程。将伯努利方程与外加激励条件、速度势函数耦合到自由振动方程中,采用迦辽金积分法,给出了矩形贮箱在流体作用下的应力与变形的解析解。讨论了弹性底板的抗弯刚度、结构尺寸、底板材料参数及流体深度等因素对底板应力与变形的影响。研究结果表明:在液体晃动非线性激励作用下,贮箱底板的应力和变形随着液体深度、板长的增大而增大,随着板厚的减小而增大,且成非线性变化关系;底板的变形及应力与底板的材料常数相关,其中板厚的变化对其变形和应力影响要比板长及液体深度的影响显著得多。本文结果可为工程实际中矩形贮箱的设计提供参考。     

A boundary element method for small-amplitude viscous fluid sloshing in couple with structural vibration

A boundary element method is presented for the coupled motion analysis of structural vibration with small-amplitude fluid sloshing in two-dimensional space. The linearized Navier-Stokes equations are considered in frequency domain and transformed into boundary integral equations. An appropriate fundamental solution for the Helmholtz equation with pure imaginary constant is found. The condition of zero-stress is imposed on the free surface, and non-slip condition of fluid particles is 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. Some typical numerical examples are presented.     

Coupled vibration of a partially fluid-filled cylindrical container with a cylindrical internal body

In the present paper a method is proposed to investigate the effects of a rigid internal body on the coupled vibration of a partially fluid-filled cylindrical container. The internal body is a thin-walled and open-ended cylindrical shell. The internal body is concentrically and partially submerged inside a container. The radial and axial distances between the internal body and the container are filled with fluid. Along the contact surface between the container and the fluid, the compatibility requirement for the fluid–structure interactions is applied and the Rayleigh–Ritz method is used to calculate the natural frequencies and modes of a partially fluid-filled cylindrical container. The fluid domain is continuous, simply connected, and non-convex. The fluid is assumed to be incompressible and inviscid. The velocity potential for fluid motion is formulated in terms of eigenfunction expansions for two distinct fluid regions. The resulting equations are solved by using the Galerkin method. The results from the proposed method are in good agreement with experimental and numerical solutions available in the literature for the partially water-filled cylindrical container without internal body. A finite element analysis is also used to check the validity of the present method for the partially water-filled cylindrical container with internal body. The effects of the fluid level, internal body radius, and internal body length on the natural frequencies of the coupled system are also investigated.