孤立气泡生长过程的短时微重力落塔实验研究

利用中国科学院国家微重力实验室北京落塔提供的3.6s微重力时间开展了短时微重力条件下的池沸腾实验研究, 分析了微重力条件下孤立的单个气泡生长过程特征. 实验中采用掺杂磷的N型光滑硅片作为加热面(加热片尺寸10mm×10mm×0.5mm), 以含气率0.0046 (气液摩尔分数比)的FC-72作为工质, 利用恒流源对加热片通电加热. 通过对实验观测到的单个气泡生长图像及相应传热数据分析可知, 经典传热机制控制的气泡生长模型可以描述其早期特征. 相关模型中经验参数的拟合结果在文献报道的数值范围内, 表明重力对气泡生长早期影响较小, 但较大的气泡尺寸可以提供更准确的数值结果.      

升力浮力复合型飞艇动导数分析

计算了球体和椭球体的平移、转动加速度导数,数值模拟了NACA0015翼型的俯仰振荡,验证了计算流体动力学(CFD)方法计算动导数的有效性.模拟升力浮力复合型飞艇的强迫振荡运动,得到了全艇以及机翼、艇身、平尾和垂尾的动导数,并与常规飞艇进行了比较分析.机翼对速度导数的贡献导致复合型飞艇升沉速度导数、滚转角速度导数较常规飞艇显著增大.常规飞艇加速度导数以惯性力成分为主,但黏性力成分以及机翼对艇身和尾翼的干扰对复合型飞艇升沉、滚转和俯仰3个方向的加速度导数有重要影响,不能将加速度导数简化为附加质量系数.      

椭圆轨道相对动力学状态转移矩阵

针对已有研究的不足,给出了更加便于理论分析和工程应用,并具有普遍适用性的状态转移矩阵的一种新形式。基于Lawden方程推导了以三维空间中位置-速度表示的状态转移矩阵元素的具体表达式,并将其由角域变换到时域,得到了表示真实位置速度的状态转移矩阵,通过初等变换可以从该矩阵得到任意需要的形式。仿真结果表明了状态转移矩阵的有效性,可将其用于任意偏心率椭圆轨道的三维轨道交会、轨道修正、编队队形设计及控制问题。     

钢球碰撞柔性杆动力分析的有效数值方法(英文)

柔性物体的碰撞动力学是工程实际中的一个重要问题。获得问题的动力响应是一个有挑战性的任务,通常通过数值计算来完成。首次采用离散奇异卷积法分析碰撞动力学问题,并通过钢球碰撞柔性杆的算例来研究几种数值方法的有效性。数值方法包括常规的有限元法、离散奇异卷积法、谱有限元法和提出的改进谱有限元法,将各种方法获得的算例的数值结果与理论解比较来研究方法的有效性。研究结果表明,采用恰当的常规有限元长度后,提出的改进谱有限元法给出的接触应力的结果最精确。另外还发现,离散奇异卷积法可以用于求解碰撞问题。     

新型重叠网格洞面优化方法及其应用

高度自动化的洞面优化方法是重叠网格关键技术之一。通过对割补法的分析,针对其洞面切割后洞边界位置不确定性及不可预知性问题,提出将其与物面距离优化准则相结合的新型混合洞面优化方法。该方法针对各物体间体网格,采用物面距离优化准则,将洞边界位置控制在物体中间的位置,使重叠区域更合理、可预测;针对背景网格,将其物面距离设置为较大值,采用物面距离优化准则将其落入物体网格内部的单元一并挖去;针对物面处重叠的网格,采用割补法对其进行切割和填补,保证物面及近壁面网格得以有效的重叠。该方法在保证高适用性的前提下提高了网格重叠质量,自动化程度高,无需人工干预。3个典型复杂流动算例计算结果与实验结果吻合良好,网格重叠区域流场变量传递正确,等值线过渡光滑,流场刻画准确,证明该方法准确可靠。     

一种挠性航天器的对偶四元数姿轨耦合控制方法

为解决复杂的挠性航天器的姿轨控制问题,对于挠性航天器的姿轨耦合动力学建模与控制展开研究。基于对偶四元数原理,推导给出一套挠性航天器的姿轨一体化动力学模型。此种模型能够紧凑描述航天器的轨道和姿态,且能够自动引入航天器平动、转动与挠性附件振动三者之间的关联耦合作用。基于此模型设计了一种自适应位置姿态跟踪控制器,该控制器能够在航天器质量特性参数未知的情况下,对其位置和姿态进行轨迹跟踪控制,并使位置和姿态误差收敛。该自适应控制器还可对航天器上挠性附件对系统的耦合作用进行估计,进而在控制输出中对其进行补偿,提高卫星控制系统的稳定性。通过仿真对控制律进行校验,结果表明该控制律对挠性航天器控制效果良好,具有一定的工程应用参考价值。     

气动热数值模拟中壁面热流后处理方法

针对气动热数值模拟中壁面热流后处理方法进行了系统的数值研究,提出了一种有别于常规微分法和积分法的壁面热流后处理方法.通过3种算例进行了对比分析,结果表明:即使在流场参数相同情况下,采用不同的壁面热流后处理方法计算的热流结果也会有较大差异,且表现出不同的网格依赖性.微分法对网格的依赖性较大,单纯提高差分精度不一定能够减弱网格依赖性和提高计算精度,还与近壁处网格分布密切相关;积分法从能量平衡积分方程出发,极大地减弱了网格依赖性,具有较高的计算精度.所提出方法具有与积分法相当的网格依赖性和计算精度,尤其在驻点区域,其网格依赖性更弱于积分法,且计算简单,提高了计算效率.      

V形兜孔圆柱滚子轴承的高速动态性能

保持架打滑和动态不稳定是航空发动机主轴圆柱滚子轴承面临的难题,为此提出一种V形兜孔圆柱滚子轴承,并对其保持架打滑及稳定性展开研究。运用Hertz接触、弹流和流体润滑及牛顿-欧拉动力学理论,建立径向刚性加载下轴承的动力学模型,在此基础上,利用变步距龙格库塔数值积分法进行动力学数值仿真,探讨了V形兜孔的几何参数对保持架打滑及稳定性的影响,分析了兜孔优化后在不同转速下保持架的打滑特性、稳定性及兜孔/滚子碰撞特性,结果表明,V形兜孔的几何参数对保持架打滑及稳定性的影响显著,在30 000~60 000 r/min的转速范围内,优化轴承保持架的打滑率明显低于普通轴承,涡动半径明显小于普通轴承,滚子对兜孔局部碰撞的力幅值和频率也明显小于普通轴承。     

基于广义Richardson外插方法的阻力预测精度分析

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

固定翼双旋弹动力学分岔特性分析

针对一种滚转稳定的固定翼双旋弹,对其非线性动力学进行了分岔特性分析,并在此基础上研究了各系统参数对其动力学分岔特性的影响。根据固定翼双旋弹非对称的特点,通过数值计算方法研究其飞行过程中平衡点随同向鸭翼安装角的变化规律,通过系统的分岔图得知系统具有三组稳定平衡点,其中只有一组平衡点为理想可行的稳定平衡点,因此需限定同向鸭翼安装角的范围以使固定翼双旋弹保持稳定飞行。在此基础上针对固定翼双旋弹弹道修正组件周期旋转和转角固定两种工作模式,通过各系统参数下的系统分岔图总结了固定翼双旋弹结构及气动力参数对其动力学系统分岔特性的影响。仿真结果表明,固定翼双旋弹的各气动力参数及飞行速度均对系统的分岔特性具有较大影响,应合理选定这些系统参数以使其具有良好的气动特性。