固体火箭发动机喷管燃气舵铰链力矩参变分析

以固体火箭发动机喷管燃气舵绕流流场为研究对象,采用数值模拟技术,通过离散求解流动方程及相关湍流模型方程,对包含舵基、舵片和发动机喷管在内的复杂流场进行了变参数分析,获得并揭示了舵基、舵片附近的流动现象和流动特点.验证了所引入的数值模拟方法的有效性,分析结果对舵基、舵片的结构安装以及控制机构的设计具有借鉴意义.     

基于数值模拟的燃气舵轴向位置优化分析

为了优化燃气舵的安装位置,应用fluent软件对燃气舵三维湍流流场结构进行了分析,在此基础上就燃气舵安装位置对燃气舵流场的影响进行了数值模拟.通过数值模拟得到燃气舵在不同轴向位置、耳片长度和舵偏角下绕流流场的复杂波系结构,绘制了燃气舵在不同舵偏角下升力的变化曲线,给出了燃气舵可以提供的升力参量随舵偏角度变化的规律.通过分析,确定了燃气舵的轴向安装位置,数值模拟和实验结果相符.     

三维舵面绕流的N-S方程计算研究

基于粘性非结构混合网格,采用格心格式的有限体积法,应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟。操纵面铰链力矩的计算结果与风洞试验结果的对比分析表明,计算方法对粘性流动具有良好的适应性。对舵面铰链力矩随迎角、马赫数、雷诺数变化的规律也进行了分析研究。     

二维舵面绕流的N-S方程数值计算研究

本文基于结构化网格,应用Baldwin-Lomax湍流模型进行绕舵面流动的二维N-S方程数值模拟,将部分计算结果与风洞试验结果进行了对比分析研究,分析了计算方法的可靠性和计算结果的合理性.为开展三维舵面绕流的N-S方程数值模拟和铰链力矩数值计算研究奠定了基础.     

三维舵面绕流的N-S方程数值计算研究

通过自动生成复杂构型的非结构混合网格。采用格心格式的有限体积法以适应不同的网格单元，应用Spalart-Allmaras湍流模型进行绕三维舵面流动的N-S方程数值模拟，将部分计算铰链力矩结果与风洞试验结果进行了对比分析。说明了计算方法对粘性流动良好的适应性,分析了迎角、马赫数、雷诺数对舵面铰链力矩的影响以及计算方法的可靠性和计算结果的合理性。     

基于Kriging模型的高超声速舵面优化设计

采用RANS方程进行了高超声速舵面的气动特性计算与平面形状优化设计研究。为了提高优化效率,采用拉丁超立方抽样试验设计方法得到样本点,并建立Kriging模型替代费时的流动数值模拟。以最大化舵面操纵力矩和最小化平面面积为目标,以阻力不增为约束条件,进行了舵面平面形状优化设计。优化结果表明,优化舵面的平面面积减少了4.78%,力矩绝对值在设计点增加了24.69%,且在较大的攻角范围均有提高。将初始舵面和优化舵面分别装配到方形截面导弹上,验证计算表明,优化舵面的操纵效率更高。     

新型推力矢量燃气舵抛离技术研究

本文针对燃气舵推力矢量存在的问题进行了抛离技术研究.在研究中采用新型技术实现燃气舵的抛离方案:即利用发动机尾流燃气的压力作为锁止力,利用弹性元件的动力作为驱动,实现燃气舵的抛离.通过调整弹性系数和预紧力,可以实现对抛离时机的自动控制.     

燃气舵外围流场计算

对燃气舵沿气流方向的横截面进行其外围的二维流场计算。计算中将燃气舵外围流场划分成三个区域,I区为钝头体头部绕流脱体激波后的亚跨声速区,Ⅱ区为脱体激波后的超声速区,Ⅲ区为膨胀波后的流动区。对各个流动区,分别用条带法、特征线法和普朗特尔-迈耶公式计算,获得各区的流动参数,为燃气舵的气动力、气动加热、烧蚀、传热、应力等计算提供流场数据。     

碳/酚醛燃气舵热化学烧蚀过程数值研究

为了研究碳/酚醛燃气舵热化学烧蚀问题,在Fluent平台上采用UDF二次开发方法,对碳/酚醛燃气舵二维非定常流固热耦合过程进行了数值仿真研究。对几何建模、边界条件、热解气体逸出过程以及热化学烧蚀导致的边界退移等问题进行了详尽的描述,并选取了合适的计算模型。对不同舵偏角下燃气舵温度分布、热解气体逸出通量以及热化学烧蚀等问题进行了分析研究。结果表明:燃气舵表面化学反应以及热解气体的逸出过程能够有效降低燃气舵壁面以及内部温度,下降温度最高值约为325K;燃气舵前缘一直是热化学烧蚀最严重区域,达到0.5mm,迎风面烧蚀量次之,而背风面几乎不存在热化学烧蚀现象;迎风面热化学烧蚀会随着舵偏角的增大而变得更为严重。     

双发飞翼布局保形尾喷管燃气舵概念设计

针对飞翼布局操纵效能低、稳定性差等不足,基于大展弦比双发动机布局飞翼无人机,开展了保形尾喷管燃气舵概念设计,分别设计了一种燃气方向舵与燃气升降舵。结合可靠的计算流体力学(CFD)计算方法,获取了全机纵向和横航向气动特性,并与常规气动舵设计进行了对比。研究表明:燃气方向舵比传统阻力方向舵具有更好的偏航操纵效率,且基本不影响全机升阻特性和俯仰力矩特性;燃气升降舵与内翼段开设升降舵时俯仰操纵效率接近,且在小迎角下具有显著的俯仰操纵效率,带来的附加阻力较小。     

直升机旋翼轴弯矩特征及应用

针对球柔性桨毂构型的直升机,分析了各支臂mΩ挥舞力动载、挥舞力静载不平衡量和摆振力动载、离心力动载、离心力静载不平衡量产生的旋翼轴弯矩的幅值、频率、相位和分布规律。对某直升机制定了飞行载荷测量方案,获取了测量数据,分析了旋翼轴弯矩的特征,建立了有效性判据,判别了无效数据,编制了实测载荷谱。某直升机的分析结果表明：在旋翼坐标系下旋翼轴动弯矩远大于静弯矩,旋翼各支臂挥舞力静载不平衡量和和离心力静载不平衡量较小;旋翼轴弯矩以1Ω载荷为主,主要由旋翼各支臂1Ω的挥舞力、摆振力和离心力产生;合成弯矩呈大小相等的旋转弯矩特征,沿旋翼轴不同周向位置的弯矩幅值相近,相位差与周向的夹角相同。     

两车速度对超车过程轿车气动特性的影响

利用动网格技术对超车过程中两种相对速度下车辆周围非定常流场进行了数值模拟,得到超车车辆和被超车辆受到的侧向力、侧倾力矩和横摆力矩在超车过程中的变化情况,并根据两种速度情况下的流场特性,分析其对气动力特性的影响.计算结果表明:超车速度越大,被超车辆最大侧向力、侧倾力矩和横摆力矩越大;超车车辆最大侧向力、侧倾力矩越大,横摆力矩越小.      

考虑转捩的跨声速气冷涡轮叶片气热耦合计算

为了研究转捩对气热耦合计算的影响,在B-L代数模型与SST(shear stress transport)k-ω二方程模型的基础上,增加了两类基于间歇因子的转捩模型:代数AGS(Abu-Gharmam Shaw)模型与一方程间歇因子输运方程.选取NASA-MARKⅡ跨声速叶片为算例,分别采用全湍流模型与加入转捩的模型进行气热耦合计算.数值计算结果与试验对比表明由于能够预测附面层中的转捩过程,采用转捩模型的耦合计算得到的结果与试验吻合最好,由于在叶片壁面附近的网格较粗,采用间歇因子输运方程的转捩模型计算的结果要逊于采用代数转捩模型的结果.      

两车速度对超车过程轿车气动特性影响的研究

利用动网格技术对超车过程中两种相对速度下车辆周围非定常流场进行了数值模拟,得到超车车辆和被超车辆受到的侧向力、侧倾力矩和横摆力矩在超车过程中的变化情况,并根据两种速度情况下的流场特性,分析其对气动力特性的影响。计算结果表明：超车速度越大,被超车辆最大侧向力、侧倾力矩和横摆力矩越大;超车车辆最大侧向力、侧倾力矩越大,横摆力矩越小。     

负载力对高压导叶控制偏差影响试验研究

针对某发动机交付试车检查高压导叶控制规律时,在转速上升和下降过程中出现导叶控制规律偏差大的情况,通过对高压导叶控制原理和试车数据分析,得出高压导叶控制系统负载力为影响偏差的主要因素。利用地面试验设备,采用对顶缸加载方式,开展负载力对导叶控制偏差影响试验,复现了发动机试车过程中控制规律偏差问题。试验结果表明:负载力越大,导叶控制偏差越大。     

热斑在1-1/2级涡轮内的非定常迁移数值模拟

对典型的1-1/2级高压涡轮叶栅,通过求解二维非定常N-S方程研究了燃烧室出口热斑对涡轮一级动叶和二级导叶流场和温度场的非定常影响。计算表明热斑的存在对涡轮级型面时均压力的影响微小,但是会显著地增大一级动叶和二级导叶型面压力随时间波动的幅度。计算印证了热斑会导致涡轮一级动叶压力面严重过热的结论,并显示在二级导叶上会出现与一级动叶类似的压力面过热现象。      

变几何涡轮叶栅叶端小翼的气动性能

针对大子午扩张变几何涡轮在可调静叶转动时旋转轴端严重恶化端区流场的问题,提出在可调静叶的机匣端部应用小翼结构的方法以克服这一问题并减少叶端间隙泄漏流动。应用数值方法和标准k-ω两方程湍流模型,并结合低速风洞试验,首先研究了可调静叶栅小翼端部流场及损失分布,并考虑了可调静叶转动的影响,随后给出了叶端凹槽状小翼结构,并评估了其气动性能以及对间隙变化的敏感性。研究结果表明：在可调静叶栅中应用叶端小翼不但可以避免可调静叶转动时旋转轴端恶化端区流场,还降低了叶端间隙泄漏驱动力,从而使得可调静叶在所有转角下都具有较好的端区流动性能,并且叶端小翼结合凹槽结构可以进一步减少间隙泄漏,总体上可调静叶栅总压损失系数降低了8.9%。     

基于非线性试验气动力的飞机静气动弹性响应分析

分别使用线性和非线性静气动弹性分析方法,对某飞机纵向静气动弹性响应特性随飞行动压、攻角、平尾偏度和纵向过载变化的趋势进行了分析,并将分析结果和飞行试验进行了比较。线性方法的气动力计算使用亚音速偶极子格网法。非线性方法由风洞试验提供刚体气动力,并使用线性气动力影响系数矩阵对其进行弹性化处理。通过对两种方法的计算结果进行对比分析可以看出:①使用刚体试验气动力的非线性分析方法能够获得和飞行试验比较一致的结果;②在线性方法所提供的结果中,翼面的弯曲变形及剪力、弯矩和扭矩等部分参数能为初步设计提供大致参考,但气动力系数的弹性增量、翼面的扭转变形、翼面的剪力、弯矩和扭矩的弹性部分可能不能为型号设计提供正确的指导;③不论刚体气动力为线性还是非线性,气动力系数随动压变化的趋势均呈线性。     

Aerodynamic force and moment measurement of 10° half-angle cone in JF12 shock tunnel

An aerodynamic force and moment measurement was conducted in JF12 long-testduration detonation-driven shock tunnel of Institute of Mechanics,Chinese Academy of Sciences.The test duration of JF12 is 100–130 ms.The nominal Mach number is 7.0 and the exit diameter of the contoured nozzle is 2.5 m.The total enthalpy is 2.5 MJ/kg which duplicates the hypersonic flight conditions of Mach number 7.0 at 35 km altitude.The test model is the standard aerodynamic force model of 10° half-angle sharp cone.The length of the test model is 1500 mm and the weight is 57 kg.The aerodynamic forces were measured with a six-component strain balance.The angles of attack were set to be à5°,0°,5°,10° and 14°,respectively.The experimental results show that in the 100–130 ms test duration,the signals of strain balance have 3–4 complete vibration cycles.So,the aerodynamic forces and moments can be obtained directly by averaging the signals of balance without acceleration compensation.The force measurement error of repeatability of JF12 is less than 2%.The aerodynamic force coefficients of JF12 are in good agreement with those of conventional hypersonic wind tunnels.For this test model at Mach number 7.0 and total enthalpy of 2.5 MJ/kg,the real-gas effects on aerodynamic force characteristics are not very evident.     

平板上有叶片绕流的源项模型

将源项模型应用于叶片涡流发生器的数值模拟中, 以叶片在当地流场中产生的升力作为源项一个分量原型, 并引入了法向力作为源项另一个分量原型.该源项由叶片的形状、大小、位置以及来流条件等参数决定.在数值模拟中, 叶片位于平板上, 其高度约为当地边界层厚度的1/5.在相同计算条件下, 对"升力-法向力"源项模型计算结果和全网格计算结果及实验数据进行比较.源项模型可以得到与全网格数值模拟基本一致的诱导速度场和涡量场, 而网格量减少了65%.