求旋转椭球形贮箱内液体小幅晃动基频的一种新方法-动力学与控制学报.PDF

第９卷第３期２０１１年９月 动力 学 与控 制 学报 Ｖｏｌ．９Ｎｏ．３ １６７２６５５３／２０１１／０９ ／２４９４ Ｓｅｐ．２０１１ ⑶ ＪＯＵＲＮＡＬＯＦＤＹＮＡＭＩＣＳＡＮＤＣＯＮＴＲＯＬ 求旋转椭球形贮箱内液体小幅晃动基频的一种新方法 杨旦旦　岳宝增 （北京理工大学宇航学院力学系，北京　１０００８１） 摘要　针对旋转椭球形贮箱内液体小幅晃动，建立了原点位于与箱于箱内静液面接触线处相切的圆锥的顶 点的球坐标系，用高斯超几何级数解析表达速度势和波高的模态函数，通过伽辽金方法把变分方程转变为 一个标准的特征值问题形式的频率方程，求解了不同尺寸比例的椭球形贮箱在不同的充液比和不同的Ｂｏｎｄ 数情况下液体小幅晃动的基频，把所求结果和已有的理论和实验结果进行对照． 关键词　变分原理，　高斯超几何级数，　晃动基频 化为只需在箱壁与液体接触线处满足的简单边界 引言 条件，这样就有利于利用伽辽金方法求解液体晃动 航天器上大量液体的存在及其工况的多样复 频率．严格意义上说，这是一种近似的方法，但由于 ［１］ 杂性，导致液体的低频晃动可能与航天器的结构模 即使对很小的接触角存在的误差都很小 ，因此不 态或控制系统的特征频率相交藕，使航天器产生共 碍于问题的求解．所得结果和其他方法计算所得结 振等动力学不稳定现象，因此微重力下液体晃动的 果及实验结果进行照，验证了本文的方法的正确 频率一直是航天器动力学与控制问题研究的热点． 性． 微重力情况下，考虑表面张力的效应，静止液面是 １　液体小幅晃动的模型 向上弯曲的，如图１所示． １．１　计算模型 使用如图２所示的球坐标系统，静液面，受扰 液面，和箱壁用下式表示： 图１　微重力下旋转椭球形贮箱内静液面形状 Ｆｉｇ．１　Ｍｅｎｉｓｃｕｓｉｎｓｐｈｅｒｏｉｄｔａｎｋｓｉｎｍｉｃｒｏｇｒａｖｉｔｙ 本文基于微重力下液体在贮箱中的静止形状， 用高斯超几何级数作为液体势函数和波高的模态 函数研究了椭球形贮箱内液体的晃动基频．为避免 柱坐标系出现的静止液面轴向坐标是径向坐标的 图２　旋转椭球形贮箱和坐标系统 二值函数和满足箱壁处液体位移与箱壁相切的相 Ｆｉｇ．２　Ｓｐｈｅｒｏｉｄｔａｎｋａｎｄｃｏｏｒｄｉｎａｔｅｓｙｓｔｅｍｓ 容性条件，采用了坐标原点位于与箱于箱内静液面 Ｍ：Ｒ＝Ｒ（）； Ｍ θ 接触线处相切的圆锥的顶点的球坐标系，通过对拉 Ｆ：Ｒ＝Ｒ（，，ｔ）＝Ｒ（）＋（，，ｔ）； Ｆ θφ Ｍ θ ζθφ 普拉斯方程的推导，得到液体势函数和波高函数的 Ｗ：Ｒ＝Ｒ（） Ｗ θ 解析模态表达，同时把箱壁上须满足的边界条件转