火箭橇系统的摩擦力分析与计算

根据火箭橇动态试验非线性和非定常特点,以及火箭橇与轨道的摩擦特性,对某型火箭橇系统进行了动力学分析、建模和流场数值模拟,界定了不同运动状态的动力学特点,得到了系统的气动力特性。研究结果表明:在60～90m/s速度条件下,该型火箭橇阻力系数约为0.58,升力系数约为0.003,气动力主要表现为气动阻力,升力与火箭橇自重相比相对较小,约占自重的0.5%～1.2%。同时结果表明,火箭橇运动速度是影响摩擦因数的主要因素,摩擦因数随运动速度的增大而减小。在数值模拟和试验数据基础上,确定了幂函数形式的摩擦因数计算公式,公式以运动速度为底数,系数为2.554,指数为-0.756。并在多种工况下,对火箭橇运动参数进行预测和验证,与试验数据符合良好。      

封闭水域内固体火箭发动机尾流结构数值研究

为了研究固体火箭发动机在封闭水域中工作时尾流场特性参数变化规律,建立三维空间模型,忽略重浮力等体积力及燃气与水之间的换热传质,采用Realizable k-ε湍流模型与VOF(Volume of fluid)多相流模型,针对试验水箱内固体发动机在不同水深条件下的工作过程进行数值模拟,获得尾流流场燃气射流形貌,以及水箱壁...     

亚/超声速楔状流层流边界层速度与温度相似解及拟合解

利用相似变换获得了楔状流层流边界层无量纲流函数的3阶非线性常微分方程,用Runge-Kutta法求解微分方程获得了不同楔形角楔状流层流边界层无量纲速度随相似变量的变化曲线;推导了亚声速和超声速楔状流层流边界层无量纲温度关于相似变量的2阶线性齐次和非齐次微分方程,获得了温度分布的通解,恒壁温条件下亚声速楔状流和绝热壁面条件下超声速楔状流层流边界层无量纲温度解析解及指数函数形式的拟合解.以楔形角为0为例利用相似变换研究了超声速条件下气体压缩性及黏度随温度变化等因素对层流边界层速度与温度的影响,得出不可压缩常物性与可压缩变物性条件下无量纲速度相对误差绝对值小于9.8%的结论.研究表明:Pr越大贴近壁面处无量纲温度变化越剧烈;超声速条件下壁温低于绝热壁温时黏性耗散作用可以使层流边界层气体温度从壁面到主流间出现先升高后降低的变化.      

单轨火箭橇在轨动力特性数值分析

为系统性分析火箭橇在轨动力特性,采用三维Eluer-Bernouli梁单元对火箭橇系统进行离散,通过重构生成长程不平顺轨建立考虑轨道不平顺及滑靴磨损的靴轨非线性接触力模型,由Newmark-β结合Newton-Raphson局部迭代求解非线性动力学方程获得火箭橇在轨动力特性数值解,并通过试验验证,计算结果表明：马赫数为2的速度下火箭橇质心竖向过载峰值约为2 700g;火箭橇在轨占空比与航向速度成正比,马赫数为2的速度下竖向火箭橇在轨占空比为16%;高速段滑靴磨损量占全弹道磨损量的84%,靴轨单次接触磨损量与航向速度3次方成正比。     

地面效应对火箭橇发动机尾喷管流场特性的影响研究

摘要：基于火箭撬地面试验的需要，设计了大载荷高加速的短轨火箭橇动力系统，并对地面效应对火箭发动机尾喷管流场的影响进行了理论和试验研究。文中采用FLUENT计算软件对不同条件下的火箭发动机尾喷管外流场进行了三维数值仿真，分析了单枚火箭以及双枚火箭工作时，地面效应对发动机尾喷管流场特性的影响。针对本文研究的大载荷高加速的短轨火箭撬动力系统，地面效应的影响和双尾喷管之间的相互作用都是在1.5ms时开始显现，在2ms~3ms时比较明显，在5ms~20ms开始减小，最后形成稳定的尾部流场。对于喷管出口直径275mm、轴线距离地面不小于472mm和双喷管轴线距离不小于680mm的火箭橇动力系统，在火箭橇运动起来20ms之后，地面效应影响和双火箭发动机之间相互作用变得非常小。     

基于动网格的火箭橇刹车系统仿真研究

为了解决火箭橇刹车系统工作时的非稳态特性在流场分析中难以采用相对 运动理论及迎风阻力的算法问题,采用Fluent 动网格技术对缆绳水刹车以及戽斗水刹车 进行了非稳态工况仿真计算研究,得到了缆绳水刹车的刹车力特性曲线以及戽斗水刹车 的刹车力随时间变化曲线,缆绳刹车仿真结果与理论解析结果进行了对比吻合较好。戽 斗水刹车进行了仿真和试验值对比, 结果表明水戽斗在750m/s 速度下仿实误差为 3.92%,在590m/s 时仿实误差为6.01%,该计算为火箭橇刹车工程设计提供了算法参考。     

垂直交叉双圆柱绕流强迫振动数值模拟

采用虚拟体法数值模拟了平面单圆柱绕流强迫振动,以及空间垂直交叉双圆柱绕流下游圆柱强迫振动流场。验证了单圆柱强迫振动中的锁定状态以及相位突变现象,从而证实了该数值方法模拟振动流场的可靠性。研究了在雷诺数Re=150、间距为5倍圆柱直径、下游圆柱按正弦曲线振动时,对尾流场的影响。下游圆柱两端尾流场在振动的作用下,涡的横向间距增大,而中心尾流场由于受上游圆柱尾流的影响而保持原先的状态。下游圆柱在锁定区的振动使尾流场变得稳定,表现在流向二次涡结构的减少,以及尾流场的速度场频谱趋向单一化。     

基于马蹄涡模型的飞机尾流介电常数分布的快速建模方法

为便于研究飞机尾流电磁散射特性,提出了一种飞机尾流介电常数分布的快速建模方法。根据马蹄涡模型得到飞机尾流的诱导速度场,并利用Euler方程和等熵流方法求出了尾流的热力学参数分布,进而得到介电常数分布。同时,提供了一种将针对一种飞机计算出的结果外推到另一种飞机上的方法,避免了重复计算。该方法具有计算量小、涉及参数少的特点,是一种快速简便地模拟飞机尾流介电常数分布的方法。     

尾涡中涡强与轴向速度的相互关系

大型飞机的翼尖涡会对飞入它尾流中的其它飞机造成巨大的危害。文章希望找到尾涡核中的轴向速度超过自由流速度的条件，从而深入了解这种尾涡的特性。在风洞实验中，我们采用的机翼模型的型号是NACA0015，它的半翼展展弦比为0．80。我们利用三传感器热线探针来测量尾涡核中的轴向速度，并将之与测得的涡的环量进行了比较，结果发现，轴向速度与一个无量纲环量参数呈线性关系。当这一环量参数较小时，轴向速度小于自由流速度，而当环量增大时，轴向速度也将逐渐增大，甚至能够超出自由流速度70％。我们还发现，轴向速度的大小对实验中所采用的两种端帽结构——平型和圆型十分敏感。另外，如果利用涡的半径来修正这一环量参数，那么两种端帽结构下的数据将会重合为一条直线。     

自由射流中飞机进气道前方亚声速流场数值仿真研究

为了预测航空发动机高空模拟自由射流试验中飞机进气道-发动机组合体前方的亚声速流场特性,优化试验舱的气动设计,采用 CFD 方法在亚声速自由射流和真实大气飞行条件下对某战斗机进气道的外流场进行数值模拟。分析了进气道对前方气流压缩作用与飞行马赫数关系,对比和分析了分别在自由射流与真实大气中进气道前和进气道入口处的马赫数分布,确定了3种马赫数下进气道在自由射流中的最佳安装位置。比较发现,亚声速自由射流对真实高空大气飞行进行模拟,可以得到马赫数相似的流场。     

三维钝体结构风压分布与绕流特性的大涡模拟研究

针对风场作用下三维钝体结构的绕流特性与空间风压分布特性,采用大涡模型对其进行了数值模拟研究,以平均风速和脉动风速作为入口边界条件,在稳态分析的基础上对9种不同来流方向的工况进行了瞬态动力分析,再现了钝体绕流的冲撞、分离、重附着等现象,揭示了钝体绕流主、次涡循环区的形成、发展机理,提出了空间风压场分布的一般规律。研究表明:钝体绕流对流场风速剖面影响的强弱程度依次为尾流区、屋顶区、来流区,且对中心区域流场风速剖面的影响要大于两侧流场,并由中心区域向两侧逐渐减小;来流区钝体绕流特性与空间风压分布特性较屋顶区和尾流区稳定,尾流区涡旋的非对称生成、脱落对钝体左右两侧流场的冲撞、分离、重附着等具有重要的影响;来流方向与钝体呈45°或135°夹角时,对尾流风场的影响范围最小,但随着来流方向逐渐垂直于钝体,尾流影响区域也逐渐增大。     

斜流压气机串列转子流场特性分析

以斜流压气机串列转子为研究对象,运用CFD软件进行了数值模拟,获得了该压气机在100%和80%设计转速下S1流面流场、子午流场、阻塞工况及近失速工况的流场特性,为斜流压气机串列转子的设计和性能分析提供参考。研究结果表明:该斜流压气机在设计转速(69900 r/min)下超声速特性明显,特性曲线较陡峭;当转速小于80%设计转速时亚声速特性明显,特性曲线较平缓。随着转速的减小,压气机的稳定裕度逐渐增大。该斜流压气机串列转子叶根和叶中截面的损失主要来源于叶型尾缘的掺混损失和叶型吸力面的小范围激波损失;而叶尖截面的损失主要来源于叶型吸力面的大范围激波损失、激波与附面层的相互影响的损失和激波与叶尖泄漏流相互作用的损失。该斜流压气机进入近失速工况后,前排转子的激波强度进一步增大,并且叶间存在大范围低速区,造成流动损失进一步增大。     

可变形翼型超声速非定常气动力的研究

由小扰动线化方程出发,利用基于格林公式的边界元方法求解二维超声速流动速度场,进而对变形翼型气动特性进行理论分析和数值计算。计算结果显示虚拟质量力导致的非定常附加升力对厚度变化远比弦长或弯度变化敏感;弦长或弯度变化诱导的非定常效应影响很小,此时非定常气动升力和准定常气动升力近似相等,和不可压缩或亚声速时相类似,变形翼型的非定常气动升力近似等于准定常计算结果叠加上非定常附加升力,与具体的变形历史过程无关。同时还分析了超声速时翼型往复变形时的气动特性,和亚声速情况完全不同,随着马赫数升高,完成往复变形所需要外界输入的功逐渐减小,趋近一个常值。另外翼型完成往复变形所需做的功将不依赖于来流攻角,且与翼型参数的变化加速度a存在简单的数学关系:A ∝a～(1/2)。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

高速压气机叶栅纳秒脉冲等离子体流动控制仿真研究

为研究纳秒脉冲等离子体气动激励在高亚声速来流条件下抑制压气机叶栅流动的分离机制,建立了基于唯象学的模拟纳秒脉冲介质阻挡等离子体气动激励特性的热源模型,在微秒量级时间尺度上分析研究了纳秒脉冲等离子体气动激励对叶栅通道流动结构的影响机制,并初步探究了纳秒脉冲等离子体气动激励的流动控制规律。研究结果表明:基于唯象学的热源模型能够较好地模拟纳秒脉冲等离子体气动激励诱导产生冲击波的气动特性;纳秒脉冲等离子体气动激励诱导产生的冲击波在高亚声速来流条件下能够对叶栅通道流动结构产生较大影响,其影响规律与激励特征和流场特性有关;高亚声速来流条件下,在叶栅通道中施加纳秒脉冲等离子体气动激励能够降低通道出口总压损失,改变流场结构。     

火箭橇水刹车戽斗入水运动的阻力预测

火箭橇水刹车戽斗作为火箭橇高速试验的减速装置,附加质量较大,运行速度快,在利用水的冲击作用减速时,承受很大的载荷。为了确保火箭橇水刹车戽斗的安全,同时达到预期的减速效果,对火箭橇水刹车戽斗的运动情况进行数值仿真分析。借助Abaqus中的流固耦合技术,通过CEL方法对火箭橇水刹车戽斗入水运动情况进行分析,得到火箭橇水刹车戽斗在多种不同运动情况下的阻力曲线,并根据流体阻力理论计算公式,拟合出该火箭橇水戽斗入水的流体阻力计算公式。利用该公式可对该火箭橇水刹车戽斗在不同速度、不同初始水位高度下所受的阻力进行预测。     

飞行参数对射流矢量喷管内流场影响的数值模拟

采用三维雷诺平均N-S方程和k-ε湍流模型对固体火箭发动机燃气二次喷射矢量喷管内外流场进行数值模拟,研究飞行参数对矢量喷管内主/次流耦合流场的影响。结果表明,自由来流造成喷管出口附近压力变化,通过喷管内壁面附面层作用影响喷管内部流动特性。飞行马赫数为超声速,或者环境压强较低时,侧向力大小与飞行参数无关。在低空亚声速流中,侧向力随飞行马赫数增加迅速减小,随环境压强降低而减小。     

跨声速涡轮导叶吸力面换热特性数值研究

为了了解涡轮导叶吸力面在跨声速条件下的换热特性,采用数值模拟的方法,分析了出口马赫数对平面叶栅内流场与换热特性的影响,以及气膜出流对吸力面气膜冷却特性的影响。结果表明,跨声速条件下,斜激波导致的逆压梯度导致了吸力面层流边界层分离和转捩;亚声速条件下,吹风比从0.5增大至1.5时,转捩位置前移了约0.1倍弦长;跨声速条件的转捩位置随吹风比增大未发生变化,但是边界层分离现象被抑制,分离泡的尺寸明显变小。在吸力面小吹风比更容易获得更高的冷却效率;边界层的分离导致冷效率分布不同于亚声速条件,在分离区冷却效率迅速降低,在吹风比0.75时降低约50%。     

ARJ21飞机遭遇不同前机尾流的响应和安全性研究

为研究在实际飞行中国产ARJ21飞机遭遇前机尾流时受到的气动力、力矩及飞行安全问题，基于经典尾涡模型和空气动力学响应模型，对ARJ21飞机遭遇尾流的响应和安全性进行了分析。首先，通过尾涡耗散模型和速度模型，建立了不同前机尾涡速度场；然后，根据ARJ21飞机空气动力学响应模型进行求解，得到其升力和滚转力矩；最后，将计算结果与过载增量指标和滚转力矩系数指标进行比较，分析了国际民航组织(ICAO)规定的中型机尾流间隔下ARJ21飞机的飞行安全性。研究结果表明：当ARJ21飞机以ICAO规定的中型机尾流间隔遭遇尾流时，过载增量均小于0.15(即处于无颠簸状态),最大滚转力矩系数均小于0.048;当滚转力矩系数在0.048～0.070时，相比于ICAO规定的中型机尾流间隔，ARJ21飞机跟随中型机的尾流间隔可以缩减2.2～2.7 km,跟随重型机的尾流间隔可以缩减0.7～1.8 km,跟随超重型机的尾流间隔可以缩减2.4～3.4 km。     

空中加油吊舱尾流场计算分析

为了得到空中加油吊舱尾后流场的分布特性,采用数值模拟方法对某型吊舱尾流场特性进行了仿真计算,分析了飞行高度和速度对加油吊舱尾流场分布的影响,并对计算结果进行了验证。研究结果表明,高度的变化对流场中速度矢量和涡量分布的影响较小,速度的变化对静压速度矢量和涡量的分布影响较大。