基于udp通讯协议的贮箱液体晃动联合仿真分析.docxVIP

基于udp通讯协议的贮箱液体晃动联合仿真分析 摘要：首先介绍了基于UDP通讯协议实现的贮箱液体晃动联合仿真方法.该方法采用FLUENT计算运动贮箱的液体晃动力、力矩、质心结果，星体动力学程序在获取贮箱晃动结果后通过控制算法和策略计算星体运动的线速度和角速度.两者之间的双向数据交换，通过C/C++语言开发的通讯接口来实现.然后借助上述联合仿真方法，开展了贮箱晃动联合仿真分析，并与在轨遥感数据进行了对比，然后分析了俯仰0～25°机动工况的4贮箱液体晃动特性，最后对比分析了带孔和无孔叶片在星体动力学控制下的液体晃动特性.非沿轨迹工况下，2、3号贮箱的晃动力和力矩峰值大于1、4号贮箱，4个贮箱的Y向晃动力峰值在1.05～1.96 N，力矩值在0.05～0.36 N·m;带孔叶片对液体的搅动作用要弱于无孔叶片.无孔叶片搅动的气泡出现了明显的拉伸与离心运动，并导致更大的流体流动速度. 0 引言在微重力环境下，充液航天器内液体晃动及其对控制系统的影响是当前国内外航天技术研究的重要课题．贮箱液体晃动对航天器的结构强度和运动状态都具有重要影响由于受到实验成本、数据测量等客观条件限制，目前很难采用实验方法进行在轨航天器贮箱液体晃动动力学行为研究．当前普遍采用数值仿真或理论方法来研究上述问题．数值方法一般采用连续流模拟方法对特定贮箱运动状态下的液体晃动特性进行计算和分析在与星体动力学状态的耦合控制方面，目前一般都采用理论方法来对贮箱晃动特性进行建模，然后与航天器控制系统进行耦合．基于液体小幅晃动假设的等效晃动力学模型是早期的理论模型．随后由于需要对大幅液体晃动条件(交会对接、大角度姿态机动等)进行分析，人们又陆续发展了大幅晃动等效力学模型，如复合摆等效力学模型综上所述，目前在卫星闭环姿态控制数学仿真中，一般采用等效力学理论模型对液体晃动特性进行求解1 联合仿真方法和数学模型1.1 星体动力学与贮箱液体晃动联合仿真本次研究的贮箱液体晃动联合仿真平台是由CFD软件FLUENT与星体动力学控制程序构成．FLUENT计算运动贮箱的液体晃动力、力矩、质心结果，星体动力学程序在获取贮箱晃动结果后通过控制算法和策略计算星体运动的线速度和角速度．两者之间的双向数据交换通过C/C++语言开发的接口来实现．星体动力学模块对贮箱受力的反馈控制作用通过接口程序传递给各贮箱晃动(FLUENT计算)的线速度和角速度功能来实现．本次研究开展的联合仿真数据交换量较小，更注重数据交互的实时性，因此采用了UDP协议作为Server端星体动力学计算和Client端液体晃动仿真数据的手段．FLU-ENT采用UDF基于C编译相应的宏来实现与星体动力学程序的数据交换、贮箱液体晃动并行计算．星体动力学程序与FLUENT软件的联合仿真平台组成见图1.星体动力学程序是以动态链接库形式进行集成，需要通过C/C++建立Server端程序将星体动力学动态链接库导入;同时为星体动力学仿真定义各种初始条件，并按照指定的时间间隔导入贮箱晃动仿真数据．与其它通信协议相比，UDP协议还具有程序架构简单、理论上可支持任意数量Client挂载Server等的优点．图2给出了星体动力学联合仿真框架．星体动力学程序与贮箱晃动CFD算例之间基于UDP通讯框架如图3所示．1.2 贮箱液体晃动数学模型不可压、不相溶、等温的气液两相流NavierStokes方程可表示为式(1)为连续性方程，式(2)为动量守恒方程．式(2)考虑了压力、湍流、表面张力、重力对两相流动的影响．u为流体速度，ρ为密度，F和g分别为表面张力和重力，μ为粘度．流动粘度和密度可根据气液体积分数值进行加强平均获得．(2) VOF方法不同流体分子间作用力的差异将导致相界面上出现表面张力，利用体积分数方法(volume of fluid,VOF)对气液两相及相界面进行描述根据体积分数ε可将流场分为3个区域:气相区(ε=0)、液相区(ε=1)、气液界面区(0ε1)．气液界面的流体属性，包括流体粘度、密度采用体积分数的线性加权平均来进行计算．速度压力耦合采用SIMPLE算法，使用PRESTO!压力插值方案，动量方程离散采用二阶迎风格式方案．为提高相间界面的清晰度，流体体积分数求解采用second-order离散方法．2 贮箱液体晃动联合仿真分析2.1 PMD贮箱建模与参数设置卫星推进分系统采用板式全管理表面张力贮箱．星体4贮箱布置结构见图4.MMH和MON-1推进剂为对角布置．假定1、4贮箱推进剂为MMH,2和3贮箱推进剂为MON-1.该贮箱主要由贮箱壳体和板式推进剂管理装置(板式PMD)组成．板式PMD是该贮箱的核心部件．为了简化结构，降低CFD计算耗时，本次研究将导流板、底部蓄液叶片