基于CFD动网格技术的飞艇动导数计算方法

由于飞艇和飞机这两种飞行器在飞行原理、气动外形等方面都存在较大的差别,因此现有的飞机动导数工程估算公式和经验数据并不适用于飞艇动导数求解。为了解决上述问题,提出了一种新的基于CFD动网格技术的飞艇动导数计算方法。通过非结构化动网格技术,建立飞艇在流场中的强迫振动模型,采用有限体积法对N—S方程求解,研究飞艇在强迫振动中的气动特性变化规律,并基于数值计算结果求解飞艇的动导数。研究发现,得到的计算结果与风洞试验结果基本一致,表明基于CFD动网格技术进行飞艇的动导数求解是有效可行的。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

基于CFD方法的俯仰静、动导数数值计算

应用时空二阶精度的隐式迭代NND算法,数值模拟了尖锥、钝锥、弹道外形和飞船返回舱等典型再入飞行器作强迫俯仰振荡的粘性动态流场,给出动态气动力和力矩系数的时间历程,在此基础上,运用积分法和奇异分解线性最小二乘法辨识静、动态气动稳定性参数,并与试验结果及半经验理论预测方法进行了比较,表明本文方法具有准确度较高、效率较高及适应性广等优点.     

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

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

基于准定常的飞行器动导数的高效计算方法

结合工程设计经验和工程计算方法,针对飞行器动导数计算方法展开研究。基于求解准定常欧拉方程发展了一种适合于亚跨超声速状态下,任意气动外形气动导数的高效计算方法,并将计算结果与一般非定常方法的求解结果以及国际标准实验数据进行了对比,研究结果表明该方法结果可靠、效率较高。     

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

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

基于CFD的动导数计算与减缩频率影响分析

以CFD流场与气动力计算为基础,分别利用小幅度强迫振动模式和差分两种方法推导了飞行器的纵向组合动导数计算表达式;建立了基于CFD技术的飞行器动导数计算方法;以国际标模Finner导弹为例,利用所推导的两种方法进行该导弹的动导数计算,重点分析了减缩频率对动导数计算的影响,提出了利用CFD方法开展动导数计算时减缩频率的选择...     

飞行器大迎角下俯仰静、动导数的数值计算

 采用一种适用于非定常流动计算的逐次超松弛（SOR）时间推进方法耦合求解非定常Navier-Stokes方程和强迫运动方程,数值模拟了NACA0012翼型、弹道外形(HBS)和有翼导弹等飞行器在不同初始迎角下做强迫俯仰振荡的黏性动态流场,给出动态气动力及力矩系数的时间历程;采用积分法获得飞行器的静导数和动导数,将计算结果与实验结果及半经验理论预测方法进行比较并分析飞行器在大迎角下动态特性的非线性特征;计算结果与实验值相吻合,表明本文采用的非定常数值模拟方法在大迎角下能够较准确地模拟飞行器的动态特性,具有工程应用价值。     

有地面干扰的动导数计算方法

提出了一种有地面干扰情况的飞行器动导数计算方法,其基本思路是从非定常升力面气动力理论出发,利用飞行器的低频运动的气动力,通过振动频率的级数展开,选用适当的气动力部分进行积分,得到所需的气动力导数,为阿了面的影响,将传统的镜像法结合到非定常用功力理论理论方法中,计算结果表明,地面干扰随距地面的高度变化而迅速变化,而且某些动导数的符号都发生了变化,说明地面干扰城飞行动力学的分析是应该加以考虑的。     

超音速、高超音速飞行器动导数的高效计算方法

将超音速、高超音速气动力工程计算方法综合起来,推广到非定常气动力的计算,并用于求解任意外形飞行器作强迫运动时的非定常气动力.在此基础上进一步发展了一种适合于超音速、高超音速大迎角状态的飞行器动导数计算方法和程序.计算结果与国际标准实验数据作了比较,证明了该程序的正确性和精度.研究了减缩频率和基本迎角对动导数的影响,结果表明基本迎角对动稳定性导数影响较大,而减缩频率在小迎角范围内对动导数影响很小.     

基于CFD动网格的内装式空射分离研究

研究了使用飞机发射运载火箭时的机箭分离过程.提出了一种可以解决内装式空射分离问题的整体解决方案,主要基于CAD三维实体建模方法、多体动力学方法和结合六自由度运动的计算流体力学(CFD)动网格方法.以CFD动网格计算为核心,进行了空射分离过程关键问题的研究,得到了火箭在载机气动干扰环境下的运动特性,为轨道卫星和弹道导弹的快速、廉价发射打下了研究基础.     

典型超声速/高超声速动导数计算方法研究

以当地流活塞理论为基础完成超声速/高超声速动导数计算公式推导,发展了一种高效的动导数计算方法.选取1个超声速、2个高超声速动导数算例标模,计算了超声速/高超声速下飞行器动导数变化规律.研究结果表明:BFM超声速流动算例的俯仰组合动导数绝对值随马赫数的增大而减小,计算重心位置后移亦会导致俯仰组合动导数绝对值减小,纵向动稳定性降低;在0° ～20°迎角范围内,高超声速流动中的尖锥和钝头旋成体俯仰组合动导数绝对值随迎角的增大而增大,纵向动稳定性增强;所提方法计算时间约为双时间动导数计算方法的1/27.     

一种高效的CFD/CSD耦合飞行器多学科优化设计方法

建立了一种基于CFD/CSD (computational fluid dynamics/computational structure dynamics)耦合的多学科优化设计方法.该方法采用结构有限元方程进行结构特性分析,同时采用数值求解N-S方程方法获得气动载荷,保证了设计结果的可信性,并采用基于代理模型的优化设计方法解决了采用高精度模型后计算量太大的问题.采用RBF(径向基函数)插值方法进行流固耦合界面数据插值,保证了传递过程中的能量守恒.利用建立的优化设计方法对一典型的翼身组合体进行了气动和结构优化设计,算例表明在满足约束条件的前提下优化后阻力和质量分别减小了5.0%和4.2%,证明了该方法的有效性.      

基于时间谱方法的飞行器动导数高效计算技术

针对动导数计算工程应用亟需的高效、高精度发展目标,提出求解非定常雷诺平均Navier-Stokes(RANS)方程的时间谱方法(TSM)的全隐格式,以改善采样点数较大时的数值稳定性,并将TSM离散推广到Menter剪切应力输运(SST)湍流模型,以提高TSM的工程实用性。将TSM应用于数值模拟NACA0015强迫振荡,所得计算结果与试验数据和双时间步(DTS)方法的计算结果均能较好地吻合,验证了TSM对周期运动的模拟能力。采用发展的TSM对高超声速HBS标模和超声速Finner标模进行计算,并分析研究攻角和马赫数对动导数的影响规律。结果表明:对于周期运动,具有与DTS方法相当的计算精度,但TSM的计算效率会随来流马赫数的增大而提高,其效率优势在高超声速范围时可达一个量级以上。     

机动飞行数值模拟关键技术及初步验证

建立飞机机动飞行仿真及飞行性能评估能力是未来航空CFD发展的重要目标之一,为此,对机动飞行数值模拟的方法、原理、关键技术及潜在的应用方向进行了探讨。采用惯性系下任意拉格朗日-欧拉(ALE)形式描述的Navier-Stokes方程与机体运动学方程相耦合的方法,结合动网格技术建立了通用的机动飞行数值模拟应用框架。从微分方程形式的运动学方程出发,建立了与流动控制方程进行紧耦合和松耦合计算的两套计算方案,并提出了利用预报-校正技术提高耦合计算精度的方法。计算结果表明,所提方案合理、可行,为单自由度/多自由度机动飞行模拟、动导数分析等计算提供了可靠的新工具。