fluen软件在火箭气动特性计算中的应用.docxVIP

fluen软件在火箭气动特性计算中的应用 气动特性计算过程fluen-系统二次开发 气特征计算是钻头设计过程中非常重要的环节。它的精度直接影响到后续工作的结果，如计算轨迹计算、全球参数验证、结构强度审查和可靠性分析。本文针对火箭设计过程中的气动特性计算,在计算流体力学商业软件FLUENT的基础上进行二次开发。整个过程主要由三部分组成:首先对火箭外形进行几何建模和网格离散;其次根据火箭构型不同设定边界条件,形成计算参数文件,后台启动FLUENT,自动读入网格文件和参数文件,开始气动特性计算;最后对得到的气动数据进行修正和处理得到用户需要结果。 1 fluent软件在vsc环境下的第二次开发 1.1 复杂面板参数设置 FLUENT软件包含多种计算模型,可应用于航空航天、传热、化学反应等很多领域。为了适应各个行业和模块的需求,FLUENT软件设计了复杂的面板参数设置,计算中需要用户输入的操作较多。但对特定的研究人员来说,主要关心解决某一类专门的问题,因此也不必要在繁杂的参数设置上花费过多的时间。 本文针对火箭气动计算问题,将计算时的常用参数整理成文件格式,通过FLUENT直接读入该参数文件对计算参数进行设置,自动判断火箭构型。不需要手动输入,实现了火箭气动计算的参数化和自动化。 1.2 自动计算过程 批处理可以不受干扰、按一定的顺序执行多任务,使用需要任务提交的公共计算资源,在处理单元数目很大的任务时减少内存使用。 FLUENT软件的文本用户界面(TUI)采用Scheme语言编写,不仅可以控制所有的面板参数设置,而且可以显示解算器中的隐藏属性。TUI和Scheme语言可以通过以下方法执行:1)在FLUENT的控制窗口中输入命令;2)读入包含一系列命令的文件。 本文中采用第二种方式通过周期轮询来实现流程化和自动化的过程,流程图如图1所示。针对不同构型,不同马赫数,不同攻角的计算状态生成适合该状态的参数文件,然后在后台启动FLUENT,读入参数模板文件和网格文件进行计算。多个计算状态按一定次序依次执行,实现自动计算功能。对计算结果的保存、提取和处理也是自动化的过程,大大减少了用户的工作量,从而提高了设计效率。 1.3 单元面力量和最大刚度 壁面上沿指定的力矢量的合力,等于每个面上压力与粘性力和该指定力矢量的点乘之和。即: Fa=a·Fp+a·Fv(1) 其中a为指定的力矢量,Fp为压力,Fv为粘性力。 为了减少舍入误差,可以采用参考压力来标准化单元压力值。则壁面上施加的净压力就可用每个单元面上的力矢量和求得。即: Fp=-∑i=1n∑i=1n(p-pref)An=-∑i=1n∑i=1npAn+pref∑i=1nAn(2)∑i=1nAn(2) 其中n为该壁面上包含的小面总数,A为小面面积,n为该小面上的单位法向力。 关于一个指定点A的总力矩等于每个面上压力和粘性力与该力作用点B到点A的位移矢量rAB的叉乘之和。即: MA=rAB×Fp+rAB×Fv(3) 其中A为指定力矩中心,B为力作用点,rAB为位移矢量。力矩方向遵循叉乘右手法则。 无量纲化的法向力系数CN、阻力系数CD定义如下: CN=NqA,CD=DqA(4)CΝ=ΝqA,CD=DqA(4) 其中动压q=ρv2/2。 N,D分别为法向力、阻力,ρ,v,A和l分别为参考密度,参考速度,参考面积和参考长度。 压心系数的定义如下: Xcp=CM/CN(5) 2 vier-st对型的气动特性和阻力的计算 本文计算了两种火箭构型的气动特性。构型A为弹体加四尾翼,构型B在构型A基础上增加了四个助推。来流马赫数分别为Ma=0.4,0.6,0.8,1.0,1.2,2.0,4.0,6.0;来流攻角α=4°。计算时取弹体最大截面积为参考面积,弹体长度为参考长度。 文中采用基于有限体积法的三维定常Navier- Stokes层流模型,压力选用二阶离散方法,压力速度耦合采用SIMPLE格式,密度和动量均选用二阶迎风格式,能量选用一阶迎风格式。时间推进使用亚松弛技术,压力、动量和能量亚松弛因子均设为0.1。 图2为构型A表面网格及Y-Z平面剖面网格,网格总数为141,580个网格点。图3-图5分别给出了构型A在α=4°时,法向力系数、压心系数和阻力系数随马赫数的变化。从图可以看出,气动特性计算值与实验值吻合良好。 图6为构型B Y-Z平面剖面网格,网格划分为514,740网格点。图7-图9分别给出构型B当α=4°时,法向力系数、压心系数和阻力系数随马赫数的变化。由图可以看出,构型B在亚跨音速时的阻力计算值与实验值有一定的偏差。我们对网格进行修改,并且使用不同的参数设置,结果发现法向力的计算结果相近,但计算阻力的结果有差异。采用与构型B相似的构型进行计算,也发现相似的现象。可以根据实验数据对这个稳定、有固定