轴对称降落伞小迎角稳定下降时流场特性

根据降落伞的特点,通过伞衣零厚度假设、伞轴对称假设和流场定常假设,建立轴对称降落伞的流体力学计算模型。在0°~5°的小迎角范围内,求解RNG(Renormalization Group)k-ε湍流模型下的N-S方程组,获得与有关单位试验相吻合的数值模拟结果。分析发现0°迎角时轴对称降落伞截面流场中鞍点与鞍点直接连接的状态是不稳定的。随着降落伞迎角增大,奇点类型和奇点间的连接方式将发生改变,降落伞物面附近的涡环里将出现极限环,截面流场中将出现新的鞍点和结点。迎角继续增大,极限环的面积将会减小直至消失,结点将附着在物面上,成为半结点。结果表明小迎角范围内降落伞截面流场的拓扑结构符合拓扑规律。文章揭示了小迎角范围内轴对称降落伞稳定下降时流场特性的演变规律,为进一步研究降落伞流场的流动机理和流固耦合问题打下基础。     

前体尾流对降落伞工作性能的影响

为研究前体尾流对降落伞工作性能的非定常影响,基于Realizablek-ε湍流模型采用PISO算法开展了物伞系统的非定常绕流数值计算,获得了精细的流场旋涡结构。在此基础上,研究了不同拖曳比下物伞系统的尾涡演变规律、流场分布规律以及伞衣气动特性变化。结果表明:前体尾涡导致伞衣入口处的涡量大小和方向时刻变化,随拖曳比增加,涡量黏性耗散增强,进入伞衣的旋涡强度逐渐减弱,伞衣入口形成稳定的负涡量区,伞衣尾涡脱离周期随之延长;拖曳比对尾涡区后端（伞衣入口处）流场压力的影响远大于前端,随拖曳比增加,流动形式逐渐由闭式转变为开式,流场的速度分布和压力分布更为对称,伞衣入口形成稳定的正压区,内外压差增加;当拖曳比大于9时,前体尾流对降落伞阻力系数和表面压强系数的影响减小。      

降落伞开伞过程的多结点模型仿真

根据降落伞的结构和力学特征,在轴对称假设下创建了伞衣及回收物系统的多结点结构模型.通过考虑应力,重力和气动力的作用效果,建立了用于无气流攻角平面圆形降落伞充气模拟的多结点结构模型动力学方程组.对开伞过程中的流场变化引入准定常假设,利用simple算法数值模拟求解RNG(Renormalization Group)k-ε湍流模型下的雷诺平均N-S(Navier-Stokes)方程以获得选定时刻的伞衣表面压力分布.结合多结点模型动力学方程组的解算代码和计算流体力学程序,采用流固耦合的方法对选定的平面圆形降落伞模型的开伞过程进行了动态仿真,得到了开伞过程中降落伞外形和特性的变化.通过结果分析和比较,证明了多结点模型的可行性,发展出了一种用于降落伞流固耦合计算的新方法.      

半球形伞流固耦合工程算法

降落伞流固耦合模拟非常难以实现,而一些工程设计任务要求迅速高效地完成.为了适应该要求,针对稳定下降阶段的半球形伞,提出了一种流固耦合工程算法.以轴对称伞形为伞衣初始形状,通过流场模拟得到结构模拟所需的伞衣压力系数分布,通过结构模拟得到新的伞衣外形,如此循环,直到伞衣中幅线和伞绳外形趋于稳定为止.结果表明,伞衣外形能迅速收敛,该方法能满足工程设计要求.     

降落伞伞衣载荷的性能试验

在改变伞衣面积和透气量的情况下,在定常风洞中对小型平面圆形伞的开伞过程进行了动态试验,研究了开伞过程中伞衣形状和伞衣所受载荷之间的动态关系.伞衣载荷不仅受降落伞阻力和附加质量变化率产生的力影响,同时,不稳定性及伞衣(绳)的抖动也是产生伞衣载荷的一个重要因素.稳定型伞衣充满后的载荷平均值约为最大开伞动载的60％左右,伞衣投影直径变化引起的呼吸现象也会引起载荷波动.透气量越大,降落伞稳定性越好,但降落伞的充气性能将会变差.该实验研究可为降落伞的载荷分析提供一定的依据.      

一种计算开伞动载的实用方法

依据海因里奇开伞动力理论,在计算个表充满时间的连续方程中,考虑了伞顶孔,伞顶部分大透气量织物面积,以及伞衣中、下部分织物透气量对开伞动载的影响。引用了非线性的阻力特征文化曲线［２］,计算伞衣充气过程中阻力特征的变化。并考虑了伞衣充气过程初始瞬时的阻力特征,假设在伞衣初始充气时期的一段时间中,进气口直径为常值,数值上等于［４（ＣＡ）０／Ｃ（ｓｔ）π）］（１／２）。计算实例表明,计算结果与实测数据相比较,其一致性是比较满意的。     

降落伞初始充气阶段数值模拟

首先建立了初始充气模型 ;将前体、伞绳及伞衣离散为一系列弹性连接的质量节点 ,伞衣内部流场为一维不可压流 ;最后利用模型计算了初始充气阶段伞衣形状变化、伞衣投影面积变化及伞绳张力     

不对称静子尾迹流场激振力分析及计算方法

为降低转子叶片在前排静子尾迹流场激振情况下的动应力,提高转子叶片疲劳寿命,对不对称静子叶片分布进行流场分析,研究在不对称静子分布非谐情况下,转子所受到的尾迹流场的激振力频谱及幅值.研究表明不对称静子分布非谐设计下,转子所受到的激振力由静子均布情况下的单频率高幅值激振转变为多频率成分低幅值激振力,通过对各个频率成分幅值的计算分析,得到了不对称静子分布非谐设计可有效降低转子叶片受到的尾流激振力,提高疲劳寿命的结论.同时,对不对称静子分布的三维流场计算,提出了简化计算方法,可大大降低计算时间和对计算机硬件要求.      

角区三维分离流附着鞍点拓扑结构及其演化

实验研究角区层流边界层三维定常分离附着鞍点拓扑结构的存在性及其演化规律.二维粒子图像测速(PIV,Particle Image Velocimetry)激光片光源设置于给定平板和不同模型组成的角区对称面,采用微距镜头捕捉对称面奇点附近的局部流动拓扑结构.实验证实在一定流动参数条件下多种形状模型与平板形成的角区均存在有别于经典分离的附着鞍点拓扑结构;随着雷诺数增加或者随着模型钝度的变化,附着鞍点拓扑结构与经典分离鞍点拓扑结构之间存在着一定的演化规律;对附着鞍点拓扑结构与Lighthill经典三维分离模式的关系进行了分析,表明Lighthill经典三维分离模式也适用于分析附着鞍点拓扑结构.     

分离式飞机应急数据记录跟踪系统设计与试验

分离式飞机应急数据记录跟踪系统具备智能弹射与分离、拖曳式跟踪拍摄、缓降与应急漂浮和数据传输等功能,针对弹射和缓降等过程进行了系统设计和无人机试验验证。同时,针对伞-囊组合体的特点,分析了气囊尾流区中伞衣阻力系数的变化规律。结果表明:气囊半径和伞衣名义直径是影响伞衣阻力系数的主要因素;伞衣阻力系数随气囊半径增大而下降,随伞衣名义直径增大而上升;在气动力分析和数值模拟的基础上,确定了伞衣阻力系数的计算公式。无人机试验完成各项设计功能,系统整体方案合理可行,为后续工程应用提供了重要参考。      

圆柱体减阻技术及其机理初步研究

利用氢气泡显示技术观测空心圆柱体在开有不同数目的细缝以及细缝在不同的攻角情况下的流态,对此类空心圆柱体的减阻机理进行了初步探讨.细缝是否改变了圆柱体的尾迹流场是其能否减阻的主要判据,据此得到细缝的攻角 α满足|α|≤60°时可以降低圆柱体的阻力,且当|α|=20°时,效果最佳.此时位于下游的细缝将产生明显的二维射流,极大地影响尾迹的流态.而当细缝的攻角满足70°≤|α|≤90°时,细缝对阻力几乎没有影响.      

射流角度对平板横向射流的影响

采用数值和试验方法研究了射流角度对平板横向射流流动结构和工作特性的影响,将得出的规律应用于射流控制矢量喷管上.在小型风洞试验台上进行试验,用纹影方法来观察实验模型的流场结构,通过静压测点来测量实验模型的壁面压力.研究结果表明:数值与试验结果吻合较好;对平板横向射流,增大射流角度能增大射流上游的分离区,弓形激波位置更靠前,角度增加到一定大小,流场结构变化不再明显;对射流控制矢量喷管进行数值模拟得出,增大射流角度能有效提高喷管的推力矢量性能,在NPR为 4.6,SPR为0.7条件下,射流角度从90°增加到130°,推力矢量性能提高28.3%.     

含无质量系绳的卫星系统平面运动和常规动力学

研究绳系卫星系统在同一轨道平面内运动，由子星－母星之间相对位置向量描述，亦即由该向量的长度变化和方位角的化而决定，在引入距离速率控制算法之后，系统的动力学主要由方位角运动而定，当母星轨道为圆形时，系统运动状态具有极限点，而在椭圆轨道下则有极限环，采用了非线性动态系统的方法和技术来计算极限状态和分析它们的稳定特性，用本文中方法对TSS－1以及其他绳系卫星系统方案作了数值模拟。     

解Stokes问题的区域分解算法

对用非协调混合元离散的Ｓｔｏｋｅｓ问题，建立了延拓定理并分析了一个两子域不重叠区域分解算法，利用已有的对二阶椭圆型问题有效的预处理器，借助Ｂｒａｍｂｌｅ ＪＨ等的思想，将原鞍点问题变为一个在适当内积意义下的对称正定问题，并估计了它的条件数。     

基于凯恩方程的无人机伞降回收动力学建模与仿真

在无人机的伞降回收过程中,无人机与降落伞一直都处于实时的动平衡状态,两者在伞降回收过程中的耦合关系及其复杂,因此很难建立精准的无人机伞降回收动力学模型。针对该问题,将伞降回收系统划分为降落伞和无人机分别进行处理。针对时变对象降落伞,通过阻力面积随充气时间的变化关系建立其动力学模型。针对无人机,首先,基于多体动力学思路,将其划分为左右机翼和机身的多体系统,通过平板绕流系数优化其伞降过程中的大迎角动力学模型;然后,通过偏速度矩阵将各体的动力学模型引入伞降回收系统质心;最终,基于凯恩方程推导并建立了伞降回收系统六自由度模型,并引入海拔高度和风力对无人机伞降回收的影响。通过数值仿真与实验数据的对比,可以发现两者具有较好的一致性,该动力学模型能够为无人机的伞降回收提供指导。      

飞行器进入火星大气的流场预测

针对火星着陆探测器进入-下降-着陆过程的高超声速进入阶段, 利用三维并 行程序求解流体动力学Navier-Stokes方程与化学反应动力学模型, 分析火星 科学实验室进入火星大气时探测器周围的流场结构、化学非平衡效应影响和气 动特性变化规律. 结果表明, 对于完全气体模型, 来流的热力学性质参数选 取影响激波位置和强度. 在化学非平衡效应影响下, 探测器头部激波脱体距离 大幅减小, 驻点压力变化不大, 波后温度显著降低. CO₂在激波后大量分解, 消耗相当能量. 流线结构显示, 探测器尾迹流动中存在复杂的旋涡运动等流动 分离现象.      

角区三维漩涡流动的新分离结构

分别用数值计算和PIV(Particle Image Velocimetry)实验方法研究角区层流边界层三维定常分离的流动结构,证实了角区确实存在有别于传统分离现象的附着鞍点结构,对称面上游的流线并非从壁面向上抬起从壁面分离,而是经由空间某个奇点向壁面附着.角区马蹄涡的传统分离鞍点结构和附着鞍点结构之间存在着一定的演化规律,影响参数包括模型头部钝度、边界层厚度和雷诺数.一定的边界层和雷诺数条件下随着模型头部钝度减小,角区马蹄涡将从传统的分离鞍点结构逐步过渡为附着鞍点结构;对一定的模型若雷诺数越大、边界层越厚则角区流动越趋向于传统的分离鞍点结构,反之则倾向于附着鞍点结构.     

三维抛物化N-S方程全隐格式矢通量分裂求解

采用全隐格式的矢通量分裂算法,在通用曲线坐标系内求解三维无因次抛物化强守恒型N-S(Navier-Stokes)方程.所建立的计算程序,可计算较为复杂边界的三维流场,如求解亚、跨、超音速共存的流场,分区模拟内外流共存流动及其混合剪切尾流.给出了一轴对称亚跨超音速喷管内流场的计算结果,包括计算网格,流矢分布,压力、温度和马赫数等值线分布.也给出了一喷管内外流及其混合尾流的计算结果,包括计算网格,流矢分布,无因次压力、密度和温度等值线分布.