动导数数值预测中的相关问题

动导数是飞行器动态稳定性分析、弹道设计和控制系统设计的重要参数,其预测方法主要有工程近似方法、数值模拟和风洞试验。基于动导数、交叉导数的概念,介绍了强迫振荡法、自由振荡法等预测动导数的数值模拟方法,重点就超声速、高超声速动导数数值预测中存在的问题开展讨论,包括时间步长、子迭代步数的选取;振荡频率对辨识结果的影响;强迫振荡法和自由振荡法辨识动导数可能存在差异的原因分析;交叉导数和交叉耦合导数的辨识问题以及飞行器构型对动导数预测的影响等。并结合算例进行了具体分析,在总结现有研究经验的同时,针对当前研究中存在的困惑和难题,提出了相应的改进建议。     

鸭式布局飞行器的翼体摇滚特性风洞试验

为研究鸭式布局飞行器摇滚特性,设计了一种包括鸭翼、脊型前体、边条翼、主翼和垂尾的模型,进行了自由滚转、扰动滚转、静动态测力和烟线流场显示多种技术手段相结合的风洞试验。通过自由滚转和扰动滚转试验得到了该模型翼体摇滚的时间历程,静态测力和动导数测定验证了非极限环运动形式摇滚的发生。结果表明该鸭式布局模型摇滚不仅同侧存在多个摇滚平衡点,而且在临界俯仰角,摇滚过程中可能出现从一摇滚平衡点跳动至同侧另一摇滚平衡点的突变。通过流场显示技术得到该鸭式布局模型复杂流场的基本形态分布,并对滚转角为0°时的全机涡系干扰和摇滚形成机理进行了简要分析。     

移动粒子半隐式方法GPU加速的双股射流撞击雾化模拟

为实现双股射流撞击雾化过程的高效数值求解并探究射流速度和撞击角度对雾化特性的影响规律,实现了移动粒子半隐式方法（MPS）GPU加速的双股射流撞击雾化模拟。GPU加速程序的最大加速比为16,取得了较好的加速效果。将GPU加速MPS方法应用于典型工况下的双股射流撞击雾化模拟,成功捕捉到了多尺度的液膜形成、液膜破碎成液丝继而破碎成液滴的瞬态过程,模拟得到的液膜破碎长度及雾化角度与试验较为吻合,误差分别为11.7%和0.5%,验证了GPU加速MPS方法在双股射流撞击雾化问题中处理能力。参数化分析了射流速度和撞击角度对液膜破碎长度、雾化角度及一次雾化液滴索尔直径的影响。结果表明撞击角度增加或者射流速度增加均会导致液膜破碎长度减小、雾化角度增加、一次雾化液滴索尔直径减小。     

运输类飞机典型货舱地板下部结构冲击吸能特性

为了研究运输类飞机货舱地板下部结构在冲击载荷作用下的吸能特性,选取三框两段典型货舱地板下部结构试验件开展落重冲击试验,即质量为478.5 kg的落重以3.95 m/s的速度垂直冲击倒置并固定在测力平台上的试验件,分析试验件失效模式及动态响应,同时建立有限元模型进行仿真与试验结果相关性分析及吸能特性研究。结果显示,在此种工况的冲击载荷作用下,中间支撑件发生由32框面向34框面方向的弯曲,并带动机身框发生同向弯曲和扭转,从而导致C型支撑件发生与中间支撑件相反方向的弯曲变形,并最终在机身框与C型支撑件的连接处形成两处塑性铰;紧固件失效以位于中间支撑件附近区域的长桁和剪切角片连接处的22个扁圆头铆钉发生剪切失效为主;试验初始加速度峰值和初始撞击力峰值分别为25.1g和173 kN。仿真与试验获得的结构变形模式吻合较好,仿真获得的最大压缩量与试验结果24.3 mm相差3.7%,仿真获得的压板上初始加速度峰值与试验结果25.1g相差4%。通过仿真分析发现机身框和中间支撑件是主要的吸能部件,吸能贡献分别占总吸能的32.1%和30.4%。     

一种高效的组合动导数分量数值模拟算法

基于传统网格速度法,发展了在流场网格不变的情况下,数值模拟沉浮运动、纯俯仰运动和俯仰振荡运动的算法,提出了一套直接求解旋转阻尼导数和洗流时差导数的方法,该方法不需要实时更新网格,减少了计算时间和所需内存,避免了负体积的出现。首先通过NACA 0006的阵风响应和NACA 0012翼型的强迫俯仰振荡,验证了程序和方法的正确性;然后数值模拟了BFM(basic finner missile)标准模型的沉浮运动、纯俯仰运动和俯仰振荡运动,并计算了其旋转阻尼导数、洗流时差导数和组合动导数,所得结果与动网格方法和风洞试验值吻合,最大误差不超过4%,进一步验证了该方法的正确性。     

基于频响函数的复合材料空间分布模量场识别

针对纤维编织复合材料宏观力学性能的非均匀特性,提出了基于频响函数(FRF)的复合材料梁空间分布弹性模量场的识别方法。采用基于灵敏度分析的方法构造优化问题,以实测和计算加速度频响残差范数最小为目标函数,进而通过迭代求解识别出复合材料梁弹性模量的空间分布。首先,以悬臂梁模型为研究对象进行数值仿真分析,验证识别方法的正确性。进一步开展复合材料梁模态试验研究,将复合材料三点弯曲试验获取的近似均质化弹性模量作为优化问题的初值;利用非接触测量方法获取模态试验中梁上各测点处的动位移响应,并计算得到各测点的加速度频响函数作为优化问题的输入值。结果表明:采用所提出的识别方法获取的模量场计算得到的梁上各处频响函数与试验获取值吻合,且所提方法在实测动响应存在噪声污染工况下是可行的。该方法能够为复合材料等效建模提供更加准确的弹性模量场。     

基于刚性动网格技术的动导数数值模拟

基于刚性动网格技术,建立了俯仰动导数的非定常数值计算方法.首先使用小幅度强迫俯仰振荡方法求解俯仰组合动导数,然后利用小幅度强迫升沉振荡方法求解洗流时差导数,通过两者相减即可得到俯仰阻尼导数.利用国际动导数标准模型Finner导弹进行算例验证,计算得到的俯仰组合动导数与试验值误差为2.76%,洗流时差导数值约为俯仰阻尼导数的11.5%,与文献的结果一致.结论表明:动导数单独模拟方法具有较好的工程实用价值,且可以推广到横向以及航向的动导数数值模拟.      

新型叶片-机匣碰摩模型及其验证

针对航空发动机叶片-机匣碰摩故障,提出了一种新型叶片-机匣碰摩模型,该模型在传统弹性碰摩模型的基础上,考虑了叶片数和转静间隙变化对碰摩力的影响,能够模拟机匣单点、多点、局部及全周,转子的单点、多点、局部和全周的碰摩规律.将所提出的碰摩模型运用于转子-支承-机匣耦合动力学模型中,利用数值积分获取碰摩故障下的机匣加速度响应规律.利用带机匣的航空发动机转子实验器,进行了叶片-机匣的机匣单点-转子全周的碰摩实验,发现了机匣振动加速度信号的碰摩特征具有明显的周期冲击特征,其冲击频率为叶片通过机匣的频率,在数值上等于旋转频率与叶片数的乘积,在频谱高频段出现了叶片通过频率及其倍频,冲击的大小受旋转频率调制,倒频谱具有旋转频率及其倍频的倒频率成分,仿真和实验取得了很好的一致性,验证了所提出的叶片-机匣碰摩新模型的正确有效性.最后,在此基础上,仿真计算了多种碰摩模式下的机匣振动特征和规律.      

考虑多种物理效应的钝锥俯仰稳定性参数影响分析

首先采用二阶ROE格式计算了非定常NS方程下高超声速钝锥体的俯仰阻尼导数,通过与文献中内伏牛顿法和实验结果对比验证了本文所采用非定常计算方法的正确性.同时采用二阶ROE格式及双时间步长、刚性旋转动网格技术等,数值分析了NS方程计算条件下高空高马赫数钝锥体模型的强迫振动非定常流场特征,并运用最小二乘法辨识出其俯仰静稳定性导数与俯仰阻尼导数.通过计算结果对比分析了不同高度(含稀薄气体效应)、Maxwell滑移边界条件及五组元化学非平衡模型对钝锥体模型非定常流动特征及静、动导数的影响.结果分析表明,本文所计算的钝锥模型在所假设强迫振动下的俯仰稳定性参数随高度变化剧烈,且考虑滑移边界条件与化学反应会使得俯仰力矩迟滞曲线形态发生较大程度改变.准确计算高空、高马赫数条件下飞行器俯仰阻尼导数应考虑化学非平衡效应和稀薄气体效应的影响,是否考虑这些物理效应会关系到飞行器静、动稳定的程度.     

加速度冲击对卫星接收机性能影响

基于加速度对晶体振荡器输出频率的影响,研究了加速度冲击对卫星接收机性能的影响。首先,利用ANSYS软件仿真分析了加速度冲击作用对石英晶体形变的影响过程,在此基础上建立加速度冲击环境下卫星接收机参考源晶振的频率偏差数学模型;接着,分析在加速度冲击下由参考晶振频率偏差对接收机射频前端输出中频信号的影响,并结合卫星接收机的工作原理阐述了加速度冲击对其跟踪性能的影响。仿真结果表明,加速度冲击会使接收机跟踪环路的性能恶化,时间较长会引起跟踪环路失锁。