不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

双后掠乘波体设计及性能优势分析

根据密切锥乘波体的设计几何关系,提出双后掠乘波体概念,给出了双后掠乘波体设计的参数与生成乘波体外形之间的关系。使用非均匀有理B样条(NURBS)表达包括圆和直线的激波出口型线辅助设计,研究了钝头区域可控、后掠区域可控的乘波体外形设计方法。使用CFD数值计算方法验证了设计方法的有效性,同时研究了双后掠乘波体外形的性能优势,结果表明在保持高超声速高性能的基础上适当设计外形在低速状态、纵向稳定性和涡效应增升方面具有性能优势,为大空域宽速域高超声速飞行器的研制开拓了新的途径。     

长航时高精度惯导系统姿态阵的四种解算方法

高精度惯性导航系统是保证水下或水面运载体长时间安全航行及搭载的仪器设备正常工作的关键技术之一,姿态阵解算是其核心问题。通过分解姿态阵以及分别采用Euler角和四元数描述姿态参数,建立解算姿态阵的四种方法;分别对这四种方法中的姿态模型进行扰动,并研究相应线性化误差模型的精度及其影响因素;利用数值仿真法对四种姿态解算方法的结果及对应的线性化误差模型的精度进行了比较分析,结果表明:姿态阵解算分解法显著差别于整体法,采用Euler角或四元数描述姿态参数对姿态阵解算精度的影响不显著,分解法对应的线性化误差模型的精度优于整体法。     

自适应渐消Kalman滤波算法在RTK中的应用研究

在动态载波相位差分定位(RTK)中,由于观测环境复杂,会经常发生周跳、卫星信号失锁等情况,严重影响基线解算的连续性和可靠性。针对动态应用环境,提出了一种Kalman滤波算法在RTK技术中的应用方法。该方法可以实时估计模糊度浮点解及其协方差矩阵,在需要重新固定模糊度时可直接用于搜索,起到了周跳修复的作用。此外,采用了自适应渐消Kalman滤波算法提高算法的动态适应性,并引入独立的滑动窗进行新息的收集和处理,解决了由于参考星变化或卫星信号失锁造成观测量中断而无法准确计算新息协方差的难题。仿真结果表明,该算法能够在模糊度发生变化时快速收敛,并且相对于一般Kalman滤波算法在高动态下提高了模糊度浮点解的精度,提高了后续模糊度搜索的效率和固定成功率。     

锥形间隙孔型阻尼密封动力特性分析及抑振机理

应用非定常动网格技术建立了锥形间隙孔型阻尼密封动力特性多频椭圆涡动求解模型,研究了锥形度和涡动频率对孔型阻尼密封泄漏特性与动力特性的影响,分析了锥形间隙孔型阻尼密封的抑振机理。结果表明:锥形间隙增强了孔型阻尼密封的泄漏直通效应,增加了密封泄漏量。在平均间隙相等的情况下,收敛间隙孔型阻尼密封泄漏量小于发散间隙孔型阻尼密封;在同一涡动频率下,刚度系数随着锥形度的增加而增大,阻尼系数随着锥形度的增加而减小;收敛间隙孔型阻尼密封具有较大的正直接刚度,增加了转子系统的固有频率,发散间隙孔型阻尼密封产生较大的负直接刚度;收敛间隙孔型阻尼密封具有较大的有效刚度和穿越频率,较小的有效阻尼,发散间隙孔型阻尼密封具有较大的有效阻尼,较小的有效刚度和穿越频率;在锥形间隙孔型阻尼密封结构设计时,需要同时考虑转子系统固有频率和同心度的影响。      

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

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

基于Kriging模型梯度解析解的改进一次二阶矩方法

改进一次二阶矩（AFOSM）法是一种基于功能函数梯度的结构可靠性分析方法,鉴于其对隐式函数的梯度较难求解,提出了一种基于Kriging模型梯度解析解的AFOSM方法,利用Kriging代理模型的解析表达式推导得到功能函数对输入变量的梯度解析解,为AFOSM中设计点的确定提供高精度的梯度信息。通过Kriging与AFOSM的结合,很好地解决了基于有限元模型的隐式情况下梯度计算量相当大、可靠性分析难的问题。数值与工程算例验证了所提Kriging梯度解析解的较高精确性,同时验证了所提基于Kriging解析解的AFOSM结构可靠性分析方法的正确性与较高精度。     

翼身融合布局低速验证机前缘缝翼设计

翼身融合布局是未来民机最有可能实现的非常规布局形式,其气动布局方案的验证通常采用缩比模型飞行试验的方式进行。以某翼身融合布局低速验证机为研究对象,以数值计算方法为基础,分析了其在飞行试验中存在的纵向和横向不稳定现象,提出了改善的方案——增加前缘缝翼。对此验证机进行前缘缝翼的气动布局设计、典型翼型的二维前缘缝翼设计和机翼三维前缘缝翼的气动设计,利用数值计算方法对设计结果进行纵向和横向分析。结果显示,所设计的前缘缝翼可以明显地增大验证机的失速迎角,改善其纵向力矩特性和横向特性。     

模型自由飞试验中的螺旋桨动力模拟方法研究

开展螺旋桨飞机模型自由飞失速试验需要准确地模拟螺旋桨动力状态。通过对模型自由飞失速试验 动力相似准则和螺旋桨动力模拟相似准则的分析,推导出螺旋桨飞机模型自由飞失速试验中所需要遵循的相 似准则关系;以此开展某四发螺旋桨飞机模型自由飞试验,通过控制螺旋桨转速、设计桨叶角,满足拉力系数和 前进比相等,模拟全尺寸飞机螺旋桨动力特性,获得的失速特性与飞机风洞实验结果一致。结果表明:该方法 能够准确模拟飞机失速过程的螺旋桨动力状态,可为飞机失速试飞提供重要的数据支撑,并为国内后续螺旋桨 类飞机模型自由飞试验失速试飞研究提供理论和试验方法。     

Boundary discontinuous Fourier analysis of thick beams with clamped and simply supported edges via CUF

This paper presents an analytical solution for static analysis of thick rectangular beams with different boundary conditions.Carrera's Unified Formulation (CUF) is used in order to consider shear deformation theories of arbitrary order.The novelty of the present work is that a boundary discontinuous Fourier approach is used to consider clamped boundary conditions in the analytical solution,unlike Navier-type solutions which are restricted to simply supported beams.Governing equations are obtained by employing the principle of virtual work.The numerical accuracy of results is ascertained by studying the convergence of the solution and comparing the results to those of a 3D finite element solution.Beams subjected to bending due to a uniform pressure load and subjected to torsion due to opposite linear forces are considered.Overall,accurate results close to those of 3D finite element solutions are obtained,which can be used to validate finite element results or other approximate methods.