基于VTPR方程计算N2/CF3I充填特性

采用比容平移法修正原始的Peng-Robinson状态（VTPR）方程,提高三氟碘甲烷（CF₃I）液相饱和区密度的计算精度。基于VTPR方程,结合经典的范德瓦尔混合规则,计算了N₂作为增压介质时,CF₃I充填1/2灭火瓶和2/3灭火瓶所需的N₂质量,并与试验值和文献值进行对比。结果表明,不同充填工况的计算值与试验值基本吻合,并且优于PROFISSY软件的结果。获得了充填压力分别为2.5 MPa和4.2 MPa、CF₃I不同充填密度下的总压力与温度的关系式,并计算了灭火剂热膨胀充满灭火瓶的临界温度和临界压力,该压力-温度关系式可用于以CF₃I为灭火剂的机载灭火系统工程设计。     

基于LES与FW-H方程的圆柱-翼型干涉噪声数值研究

采用大涡模拟（LES）结合Ffowcs Williams-Hawkings（FW-H）方程研究了圆柱-翼型干涉噪声的频谱特征和四极子噪声的影响.流场结果与实验吻合较好,声场获得了与实验相一致的频谱形状,翼型前缘30%弦长为主要干涉噪声源.针对计算与实验展向高度不同的问题,首先基于相关函数为三角函数的假定推导了展向修正公式,并与相关文献采用矩形函数、指数函数及高斯函数所得结果进行了对比.结果表明:不同的修正方法修正量差别高达7.6dB,三角函数与矩形函数修正效果类似,优于高斯函数,指数函数修正效果最差.分别采用固体壁面和可穿透数据面作为积分面,研究了体声源的影响,结果表明低马赫数下四极子噪声的影响可以忽略不计.      

自适应笛卡尔网格Ghost Cell方法研究

在笛卡尔网格中,利用Ghost Cell方法处理浸入边界,模拟二维无粘可压缩流。针对静止物体,比较了各种不同的Ghost Cell边界条件下的熵误差,总压误差,以及阻力系数。此外为提高激波分辨率,将Ghost Cell方法与基于叉树结构的自适应笛卡尔网格算法相结合,数值结果显示本文的方法是切实可行的。     

盘形铣刀加工等基圆锥齿轮的原理及齿面求解

为了实现等基圆锥齿轮的高速、高效加工,基于等基圆锥齿轮理论,分析了盘形铣刀加工等基圆锥齿轮的可行性,建立了盘形铣刀数学模型、盘形铣刀加工坐标系,分析并确立了各坐标系间的转换关系.通过刀具曲面与齿面在共轭接触点处啮合方程的建立及求解,得到了盘形铣刀加工的等基圆齿面表达式.对盘形铣刀加工理论下的等基圆锥齿轮齿面与理论等基圆齿面进行了分析比较,进行了盘形铣刀的仿真加工.结果表明:通过盘形铣刀截形、盘形铣刀加工运动中的回转角度、倾斜角度及倾斜中心点的合理设计,盘形铣刀加工而成的实际齿面与理论齿面的接近程度达到工程要求,使用盘形铣刀可以加工等基圆锥齿轮.      

旋翼翼型非定常动态失速特性的CFD模拟及参数分析

构建了一套基于运动嵌套网格技术和可压缩RANS方程的旋翼翼型非定常流动特性模拟的高效、高精度的CFD方法。首先,发展了基于Poisson方程求解的围绕翼型的粘性贴体正交网格生成方法,并提出了基于最小距离法(MDM)改进策略的运动嵌套网格生成方法,克服了弹簧法可能导致网格畸变的不足;其次,为准确模拟由湍流分离和气流再附引起的气动力的迟滞效应,基于RANS方程、双时间方法和高阶插值格式,建立了旋翼翼型非定常气动特性分析的高精度数值方法,并采用能够较好捕捉气流分离现象的S-A湍流模型;再次,针对旋翼后行桨叶动态失速时桨叶剖面来流速度较低、迎角较大的特点,为解决低来流速度时L-B半经验模型在旋翼翼型非定常动态失速计算中的局限性,并克服可压缩方程对低速流场计算收敛困难和精度低的问题,建立了基于Pletcher-Chen低速预处理方法、FAS多重网格法和隐式LU-SGS方法相结合的高效数值方法。应用发展的方法,分别针对NACA0012、SC1095旋翼翼型静态和轻度、深度动态失速进行计算,精确捕捉了气动力迟滞效应以及翼型前缘脱体涡的产生、对流和脱落过程,验证了本文方法的有效性;最后,着重针对NACA0012动态失速状态,开展了振荡参数对旋翼翼型非定常动态失速特性影响的分析,研究结果表明翼型迎角平均值、振幅及减缩频率的变化均能引起迟滞效应的改变并使得气动力峰值发生有规律的前、后移现象等。     

复杂外形网格生成技术

复杂外形网格生成是进行流场数值模拟的一个关键技术及前提。本文论述了采用保角变换方法,求解偏微分方程方法以及它们在代数方法结合生成二维及三维具有多块网格结构的贴体且与边界正交网格的技术。     

基于应变速率循环法的TA15钛合金超塑性本构方程

采用应变速率循环法对TA15钛合金进行三组高温超塑性拉伸试验,变形温度区间为850~950℃,应变速率循环区间为5×10-6~5×10-4s-1。分析拉伸试验数据后,计算出TA15钛合金动态再结晶激活能Q,结合金相组织分析得出其热变形过程中发生了动态再结晶的结论;并利用Arrhenius模型构建超塑性本构方程,应用origin数据处理软件进行数据分析,求得TA15钛合金高温条件下的超塑性本构方程。运用1stopt软件修正了该本构方程,使其精度达到99.3%。结果表明,TA15钛合金的流动应力对变形温度较为敏感,随着温度的升高,流变应力逐渐减小,软化机制愈发明显,且在900℃附近的超塑性较好,伸长率达到了846%。     

带控制律导弹投放数值模拟

基于非结构重叠网格技术,耦合舵面控制律模块,采用非定常雷诺平均NS方程,建立了模拟控制律作用下导弹投放分离动态特性的CFD数值方法.非结构重叠网格生成技术以“物面相交”方法确定洞边界,在多套网格重叠部分采用算术平均进行物理量的计算传递.采用机翼/挂架/带舵外挂物模型,以及不同的舵面控制规律展开了对比计算,分析了结果.并比较了某导弹有无控制律作用时从某战斗机分离的动态特性,研究了增加控制律的必要性.     

空间飞行器抓捕非合作目标后的复合控制

随着空间技术的不断发展，空间飞行器抓取非合作目标的稳定控制技术变得日益重要。基于拉格朗日动力学原则对空间飞行器进行动力学建模，考虑了非合作目标的未知参数、空间飞行器模型的不确定性及存在外部干扰的情况，理论分析了自适应控制及基于非线性干扰观测控制的可行性，并对比了二者的优缺点。在此基础上提出了一种可以在线精确估计非合作目标未知参数的复合控制方法，并进行了仿真验证。仿真结果表明，所提方法能够将所有参数的估计误差限制在很小的范围内，相对于自适应控制方法有了显著改善。     

纵列式直升机双旋翼气动特性分析

由于纵列式直升机前后旋翼之间会产生复杂的干扰现象,本文针对两旋翼之间的气动干扰问题进行了流场分析。首先,建立了基于动量源模型和N-S（Navier-Stokes）方程的数值研究方法,选用k?ω SST两方程湍流模型,采用基于压力隐式求解器,并对所建立求解方法进行了算例验证。随后,对悬停状态,前向来流状态,后向来流状态的纵列式双旋翼直升机前后旋翼干扰流场分别进行了数值模拟,并结合计算结果与流场特性,对干扰流场所产生的气动现象进行了分析。结果表明,针对本文算例,悬停状态前后旋翼性能均降低,后旋翼对前旋翼的性能影响较大。前向来流状态后旋翼升力损失大于悬停状态,后向来流状态前旋翼升力损失更为明显,最低处旋翼效率损失近半。