对称交叉激波和湍流边界层相互作用的数值研究

针对两种双尖鳍外形的对称交叉激波与湍流边界层相互作用,采用N-S方程和两种湍流模型进行了计算。研究了网格收敛性、鳍的角度和湍流模型对壁面压强、Stanton数和壁面摩擦力线的影响。弱相互作用的计算结果较好,强相互作用的壁面压强和摩擦力线的计算结果与试验吻合较好,而Stanton数的计算结果较差,峰值高达试验的2.5倍左右。随着鳍的角度的增加,壁面压强和Stanton数的分布从单调分布发展为M型分布,两者的峰值不在相同的位置。湍流模型对壁面压强和壁面摩擦力线影响很小,对Stanton数计算的影响很大,SST模型比BSL模型表现好一些。     

作大范围运动刚柔混合体的动力学建模与仿真

考虑耦合效应的柔性空间结构耦合动力学研究在航空航天领域具有重要的应用价值.用有限元法对柔性结构进行离散化处理,得到了柔性结构的模态阵.以模态展开方法以及Kane方程为理论基础,建立了作大范围运动刚柔混合体的动力学模型,并通过Matlab软件编写程序进行求解.为了检验该动力学模型及其程序的正确性,将Matlab仿真结果与商用软件仿真结果进行对照,进而验证了该仿真的正确性.本研究解决了中心刚体加柔性附件这一类航天器的动力学仿真问题.      

基于虚网格的格心ADT搜索法

针对重叠网格寻点问题,提出了一种直接在格心网格下操作,基于虚网格的ADT(Alternating Digital Tree)搜索方法.通过对计算网格进行扩展,建立虚边界网格,解决了格心网格因覆盖区域不完整而无法直接建立ADT数据结构的问题.搜索结果即为包含解的格心单元集合,可直接对结果列表遍历以得到合理贡献单元,故完全摒弃了可靠性差的Stencil Walk方法.由虚网格的定义,使寻点在边界附近的处理更为灵活,可以准确给出边界附近贡献单元的有效信息,同时简化了虚网格体系的构建.扩展了重叠边界类型,构建的搜索空间能完整覆盖网格范围,解决了对称面重叠问题.算例研究表明:该方法可靠性好,边界处理能力强,有效提高了重叠网格方法对大型复杂网格的解算能力.     

基于多体动力学的旋翼模型与气弹稳定性

针对无轴承等先进旋翼结构特点,建立了基于柔性多体系统动力学方法的旋翼气弹分析模型.旋翼各构件建模相互独立,便于组成不同构型旋翼,适合于新构型旋翼气弹分析.分析模型集成了适合于旋翼气弹分析的隐式大变形桨叶模型,构件动力学方程增加了铰链动力学方程,统一约束方程形式,改进了角运动约束方程,避免连续旋转引起的奇点问题.应用分析模型计算无铰式旋翼和无轴承旋翼的气弹稳定性,分析结构参数对旋翼气弹稳定性的影响.分析结果表明模型能准确计算结构弹性变形的耦合及非线性,提高旋翼气弹稳定性分析精度.      

高转速袋型阻尼密封泄漏的特性

采用有限元分析和计算流体动力学耦合数值分析方法(FEA/CFD),研究了考虑转子系统由于高速旋转造成转子径向伸长效应时袋型阻尼密封的泄漏特性和流场形态.计算了密封间隙为0.13mm时,袋型阻尼密封试验件在3种压比、3种转速下的泄漏量,并与试验值和不考虑转子伸长效应的CFD数值模拟结果进行了比较,验证了所采用的FEA/CFD数值方法的可靠性(误差小于1.3%)和高转速下考虑转子伸长的必要性.研究了6种压比、6种转速下密封间隙为0.25mm的袋型阻尼密封的流场和泄漏量,分析了压比和转速对袋型阻尼密封泄漏特性的影响规律.结果表明:当压比大于0.26时,随压比的减小,袋型阻尼密封泄漏量逐渐增大;当压比减小到一定值时(π<0.26),泄漏流体的马赫数在最后一个密封齿间隙处达到1.0,即发生了堵塞,此时泄漏量达到与密封进口总压相对应的最大值;在高转速下,考虑和不考虑转子半径的伸长,密封泄漏量均随转速的增大而减小;在转子面周向马赫数大于0.35时,需要考虑转子半径伸长对泄漏量的影响.      

考虑执行器动态补偿的动中通系统非线性鲁棒控制

研究用于移动载体卫星通信的动中通系统的稳定跟踪问题。选取动中通系统的三轴非线性模型为被控对象,较单轴模型更能准确地描述子系统之间的运动学与动力学特性。在执行器动态特性的基础上,设计了一种具有执行器动态补偿的非线性鲁棒控制器。该控制器不仅能够使系统在载体移动过程中隔离载体扰动的影响,精确跟踪目标卫星,而且还进一步解决了执行器动态特性对系统的影响,更符合实际动中通控制系统设计的需要。最后通过三组对比仿真结果验证了所设计控制器在不同扰动条件下的稳定跟踪性能。     

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

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

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

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

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

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

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

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