基于部件特性的螺旋桨数学模型通用建模算法

为了有效地建立螺旋桨数学模型和涡桨发动机数学模型,以及为涡桨发动机控制规律设计奠定基础,借鉴压气机部件特性缩比方法,提出了一种适用于螺旋桨部件特性的修正缩比方程;基于螺旋桨静态特性,提出了静拉力状态下的螺旋桨建模优化算法;针对螺旋桨低速前进状态下和高速前进状态下的不同工作特点,提出了两种不同的螺旋桨建模算法,以实现全飞行包线内的螺旋桨数学建模。通过与GSP (Gas turbine Simulation Program)软件仿真数据对比验证,其结果表明,基于所提出的算法建立的螺旋桨数学模型输出拉力、功率和效率的最大相对误差分别不超过3×10-6,3×10-6和6×10-5,同时,验证了算法有效性和通用性。     

CRM-WBPN风洞试验模型数值模拟

基于雷诺平均Navier-Stokes（RANS）方程和多块对接结构网格技术,采用亚跨超CFD软件平台（TRIP 3.0）,研究了共同研究模型（CRM）风洞模型支撑装置对CRM翼/身/架/舱（CRM-WBPN）模型压力分布及气动特性的影响。构造了一个包含小、粗、中、细网格的4套网格,并进行了网格收敛性研究,在升力系数为0.50条件下,开展了固定升力系数下模型支撑对CRM-WBPN模型和CRM-WB模型压力分布及挂架短舱阻力增量的影响,以及不同来流迎角下模型支撑对CRM-WBPN模型气动特性的影响。通过与NASA国家跨声速风洞设备（NTF）试验结果和CRM-WB模型数值模拟结果的对比分析,表明模型支撑对固定升力系数下挂架短舱阻力增量影响较小;在计算构型中考虑模型支撑的影响,引起升力系数和阻力系数下降、俯仰力矩系数增加,数值模拟结果更加接近相应的试验结果。     

无轴涵道旋翼气动特性数值研究

采用基于非结构网格的滑移网格技术对无轴涵道旋翼与传统涵道旋翼的气动特性进行了非定常Euler方程数值模拟。分别考察了无轴涵道旋翼与传统涵道旋翼气动特性差异,中心孔径、涵道扩散角、涵道翼型对无轴涵道旋翼气动特性的影响。涵道扩散角变化范围为-6°~10°,涵道翼型选择NACA66、NACA0018及NACA4415三种翼型进行研究。研究发现:涵道翼型类型对无轴涵道旋翼的拉力分配影响较大,对称翼型能减弱旋翼上方低压涡从而涵道能产生更大的拉力;无轴涵道旋翼比传统涵道旋翼具有更优的拉力性能,转速为18 000 r/min时其总拉力是后者的1.185倍;减小无轴涵道旋翼的中心孔径能提高总拉力值,但整体耗能也将随之提高;涵道拉力占比越高的翼型,其最佳涵道扩散角越大,该状态下涵道拉力占比最高。      

共轴双旋翼非定常流场干扰特性

直升机共轴双旋翼相遇过程的非定常气动干扰在其周期性干扰中最为强烈。为了研究该过程中的干扰特征,建立了一个基于非定常雷诺平均Navier-Stokes方程的气动干扰数值方法,双旋翼的反转运动采用运动嵌套网格进行模拟。为了探究双旋翼相遇时的气动干扰机理,以两个相对运动的双翼型系统来模拟双旋翼特定展向截面的非定常相遇过程。分析了双旋翼和双翼型相遇时的气动特性和流场特征,并对双翼型系统进行了参数影响研究,结果表明:双旋翼上下桨叶相遇时,上下旋翼拉力均会出现先增后减的变化趋势且波动幅值分别为其对应总拉力的30%和22%,双翼型系统的升力波动趋势与之相似;间距增大使上下翼型间的气动干扰减小,且上下翼型升力波动对间距变化的敏感时间不同;翼型相对厚度增大使双翼型升力波动幅值增大,并出现二次波动;有附加来流时,相对来流速度较小的翼型升力波动的幅值与范围更大。      

超声速绕二维波纹壁的气动与传热问题数值研究

本文采用Ｎ－Ｓ方程、代数湍流模式，对超声速绕二维波纹壁的气动与传热问题进行了数值研究，得到平缓压力分布与剧烈热流变的结果，与国外相同条件下的实验结果一致。这对于进一步深入研究小尺度旋涡的特性具有重要的意义。     

带有分离的二维可压粘性喷管流动计算

采用二维、时变、可压纳维-斯托克斯(N-S)方程计算喷管流动。流动可为无粘、层流或湍流。湍流模型可采用代数混合长模型,一方程和二方程模型。喷管形状为二维(平面或轴对称)。计算网格内点采用MacCormack显式格式。采用局部细网格,快速解法以及时间间隔光滑方法,以加速高雷诺数下收敛。边界点采用参考平面特征线格式计算,粘性项作为源项处理。给出了初步算例结果,与其它计算和实验数据比较,符合良好。     

飞机在变化风场中三维运动的数值仿真

本文在引入和定义了必要的坐标系和Euler角以后,推导了计及任意变化风场的飞机一般运动方程.所得的两组运动方程均适宜于在大气紊流和风切变中飞机三维飞行运动的数值仿真。还讨论了风梯度的气动力效应。     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

双框架控制力矩陀螺群的建模与分析

 研究采用双框架控制力矩陀螺群在卫星高精度姿态控制中的建模问题。对于中心刚体带有多个双框架控制力矩陀螺,构造了该系统的拉格朗日函数。采用拟拉格朗日方程建立了卫星姿态动力学模型。在建模中,陀螺的内、外框架未进行无惯量的假设,使得双框架运动带来的动力学特性得到了准确的反映;设计了适宜的标记,使得方程形式简洁并且统一。该模型为高精度的数学仿真和分析提供了基础。基于系统的拉格朗日函数,应用拉格朗日方程,分别建立了沿内、外框架轴和转子轴的轴向动力学方程,为执行机构伺服系统的动力学特性分析和控制器设计提供模型的依据。     

直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.