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

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

钝锥在添质作用下的稳定性导数的分析

本文在国外有关工作的基础上进一步分析了烧蚀引起的添质作用的力学模型及其对球钝锥俯仰力矩稳定性导数的影响机理,一起计入了热—化学延迟效应和气流延迟效应,考虑了边界层为层流以及转捩为湍流的各种情况,探讨了稳定性导数与添质率、边界层转捩位置、热—化学延迟相位角、减缩频率以及攻角等诸因素的关系。并且提供了计算具有添质作用的球钝锥在小攻角下以及在不同边界层流态下的俯仰力矩稳定性导数的一套方法,较前扩展了应用范围。分析表明,当添质率较大,而在重心附近又发生边界层转捩时,钝锥的俯仰阻尼力矩导数可以出现严重的负阻尼情况。本文的计算结果与某些弹道靶试验数据以及飞行试验遥测换算数据相当接近。     

再入体动导数试验技术

外形短粗，飞行姿态复杂，其动导数量级很小的再入体模型，安装在新研制的尺寸小，厚度薄的整体结构之天平元件上，天平与尾支杆相连，由气活动塞推动拨杆，经过拨螺钉撞击模型上的楔块并驱使模型作自由振动，从而测量了大攻角俯仰，偏航及滚转动导数。     

侧向喷射受控导弹的数据分析和设计优化

本文叙述的是侧向喷射受导弹的数据研究，经过数据分析研究了导弹的法向力和俯仰力矩的特征，数据来自不同的喷射流动态、攻角、喷嘴在弹体上的位置以及喷射角度。以侧向喷射受控导弹四周各种喷射气流为对象，用空气动力分析得出超音速气流的结果，把俯仰力矩和法向力作为目标函数和约束函数，把飞行马赫数，攻角和侧向喷射流体角度定为设计变量，通过采用全面优化的初始算法和响应面方法，完成侧向喷射受控导弹优化设计，得出导弹控制的最佳飞行条件。     

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.     

浅析动导数对飞机大攻角动态特性的影响

在分析飞机大攻角动态特性时必须考虑非线性空气动力影响,因为大攻角下的空气动力现象比较复杂,难以定量预测。近年来由于对飞行品质要求的提高及扩大飞行包线使用范围的需要,对大攻角的动态特性提出了更高的要求。本文简要介绍了大攻角动态特性、动导数研究方法,并就动导数对大攻角飞行特性的影响作了一些定性分析。     

某型导弹动导数风洞试验和数值模拟

动导数的精确预测对于弹箭类武器的成功研制至关重要。以某型导弹动态特性的研究为对象,在FD06风洞中采用自由振动法完成动态特性试验,获取该型导弹俯仰和滚转两个方向的典型振幅衰减曲线和动导数结果;针对该型导弹,基于双时间推进法数值求解非定常N-S方程和强迫运动方程,获取该型导弹俯仰动导数和滚转动导数。结果表明:CFD计算结果与风洞试验结果有较好的一致性;在给定马赫数和攻角下CFD数值模拟结果可以为后续工程型号研制提供参考。     

载人飞船返回舱的动稳定性

载人飞船返回舱动稳定性的地面模拟试验方法有强迫振动法。有限自由振动法，自由翻滚法和模型自由飞试验法。试验结果表明，在配平攻角区阿波罗返回舱除马赫数Ｍ〈０．７以外，一般都具有俯仰正阻尼特性。双子星座返回舱在配平攻角区，当Ｍ〈０．５时存在俯仰负阻尼特性。钝头气流分离效应，后体气再附效应，船尾近尾流交应和动态时滞效应等对静，动稳定性都有相反效应。这些效应会使返回舱的静稳定性增加，而使动稳定降低。角振幅增     

气动导数非定常数值模拟研究

非定常数值模拟飞行器的强迫振动过程,获得任意时刻的气动参数,可以解决复杂外形以及多种飞行条件的气动导数问题。通过AVICFD-X软件计算0?~80?初始攻角下的BFM带翼导弹的俯仰动导数,验证该方法的可行性。计算结果与实验值具有很好的一致性,反映了大攻角流场非线性特征。计算中充分考虑了工程实际需求,尽可能地减小计算代价,获得了理想结果,表明该方法适用于复杂外形气动导数计算。     

小迎角高超声速非旋转钝锥非线性运动分析

采用分段拟合技术对小迎角高超声速(马赫数Ma=6)非旋转钝锥双平面拍摄风洞自由飞试验结果进行了分析,获得了非旋转钝锥在小迎角高超声速下气动导数非线性的具体形式,分析了非旋转钝锥在小迎角下的非线性运动特点。研究发现在小迎角范围内钝锥的动导数系数均呈现明显的非线性,而静导数系数非线性较弱,可近似为线性。各组试验的静、动导数系数的非线性形式表明,除模型Ⅱ外的其余模型均为两方向振幅不同的极限圆锥运动,模型Ⅱ为极限平面运动。每组试验两个平面的静导数系数在小迎角范围内保持基本相等。而无论出现极限圆锥运动的四组试验还是出现极限平面运动的一组试验,模型在俯仰和偏航两个方向的动导数系数均存在较大差异。     

基于全局灵敏度分析的侧向气动导数不确定性对侧向飞行载荷的影响

考虑飞行载荷计算中使用的气动导数存在不确定性,利用基于方差的全局灵敏度分析(GSA)方法,结合偏航机动的动力学模型,分析了侧向气动导数不确定性对侧向飞行载荷的影响。以某型飞机为例,运用该方法得到侧向气动导数的全局灵敏度排序。结果表明:侧滑角、阻尼贡献的垂尾侧向载荷及垂尾侧向总载荷受全机侧力系数对侧滑角导数的影响最大,受航向静稳定导数及方向舵操纵效能的影响次之;方向舵偏度贡献的垂尾侧向载荷只受全机侧力系数对方向舵偏度导数的影响;无尾飞机侧向载荷主要受航向静稳定导数、方向舵操纵效能及无尾飞机侧力系数对侧滑角导数的影响;偏航阻尼导数基本不影响各侧向飞行载荷。同时也验证了方法的有效性,对提高飞行载荷的计算精度有一定的指导意义。     

考虑进排气的飞翼静、动态气动特性数值模拟

采用计算流体力学（CFD）数值模拟技术,以二维模型为例,建立了基于小幅度非定常运动的动导数计算方法,构建了考虑喷流引射效应的进排气模型、考虑管道效应的通气模型以及常规的保形模型,分别对其静态以及动态气动特性进行了数值模拟计算,并辨识了其纵向俯仰力矩系数组合动导数。研究表明:非定常动导数辨识方法较为准确可靠,相比保形模型,进排气以及通气模型的静态失速攻角（AOAs）增大,升力系数、阻力系数及俯仰力矩系数也增大,但力矩系数斜率基本不变,说明静态稳定性差异不大。而相同攻角下,进排气影响下的俯仰力矩系数组合动导数绝对值最大,表明了进排气模型具有最大的动态阻尼,而通气模型次之,保形模型的动导数绝对值最小。      

基于绳系并联机器人支撑系统的SDM动导数试验可行性研究

详细给出了在低速风洞中,采用绳系并联机器人(WDPR)支撑模型,用强迫振荡法进行标准动态模型(SDM)动导数试验可行性的研究。试验中将杆式六分量应变天平内置入模型中以测量模型的气动力和气动力矩,建立了适用于绳系并联机器人支撑系统的模型运动控制子系统和数据采集子系统。采用绳拉力作为参考信号,对气动力矩信号与位姿信号进行数据的同步处理,解决了绳系并联机器人支撑系统应用于动导数试验时所测力矩信号与位姿信号之间的相位差确定问题,给出了WDPR支撑下模型动导数的计算方法。整个试验样机置于某开口式低速直流风洞中进行了俯仰、带偏航角的俯仰以及升沉的动导数试验,通过测量和计算得到各动导数。试验结果与参考文献相比较具有合理的一致性。研究结果表明,采用绳系并联机器人支撑模型进行动导数试验是可行的,至少对于SDM是这样的结果;使用一套绳系并联机器人支撑系统,可以完成多套硬式支撑系统才能完成的动导数试验,从而提高试验效率,降低试验成本。     

三角翼受迫俯仰滚转耦合运动的气动特性研究

基于刚性动网格的技术,选用B-L湍流模型,利用有限控制体积法对N-S方程进行数值离散,对76°大后掠三角翼的受迫俯仰滚转耦合运动进行了数值模拟,在此基础上,对俯仰滚转耦合运动的气动力特性和流场结构进行了分析。计算结果表明:俯仰滚转耦合运动时,三角翼上表面的涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩,滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大,但对法向力影响不大。     

大攻角风洞自由飞动导数实验研究

实验是在0.2×0.2m~2超音速风洞中进行的,其来流Mach数为2.0。模型为带翼导弹模型。实验时用气体发射装置将模型向实验段上游发射;并以另一气流横吹模型使其呈大攻角状态飞行。最大攻角可达100°。采用高速鼓轮相机记录模型飞行姿态随时间的变化。用改进的数值积分法进行数据处理,求得俯仰阻尼导数随攻角的变化。     

三角翼受迫俯仰滚转运动气动特性的数值模拟研究

基于动网格的技术,利用有限体积法对N—S方程进行数值离散,选用湍流模型,对65°大后掠三角翼的受迫俯仰滚转运动进行了数值模拟,在此基础上,对俯仰滚转耦合运动的气动力特性和流场结构进行了分析。计算结果表明：俯仰滚转耦合运动时,三角翼上表面涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩,滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大,但对法向力影响不大。     

基于非定常CFD方法跨声速民用飞机标模俯仰气动导数计算研究

采用非定常计算流体力学(CFD)方法和动网格技术,计算了跨声速民用飞机全机构型TCR气动标准模型在0°~30°迎角范围内气动力系数的时均值、俯仰静导数和俯仰动导数;通过与试验结果对比,研究了非定常CFD方法在较大迎角范围内对包含鸭翼、垂尾等部件的复杂外形气动导数的计算精度。与试验结果对比表明,在16°~24°迎角范围内,受鸭翼和机翼上流动分离再附以及内外翼分离涡相互作用等复杂流动结构的影响,俯仰力矩系数的动导数与试验数据有相对较大的差异。     

平头增阻再入体俯仰动态特性计算与流动机理分析

采用自由振荡法数值模拟了一种平头增阻再入体的高超声速俯仰动态特性,沿用稳定性导数模型作为动态特性的指标,通过把稳定性导数视为攻角的函数来考虑非线性影响,并应用奇异分解线性最小二乘法辨识稳定性导数。对于本文研究对象,计算结果表明,平头外形在较高马赫数和小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的。文中对平头外形出现这种动态特性现象的空气动力学原因进行了分析,有必要作进一步的确证与研究。     

由稳态飞行确定纵向气动导数

研究了从稳态对称飞行数据估计气动导数的一种方法，要估计的导数是纵向静稳定性和操纵性导数，尾翼产生的阻尼导数、速度对升力和俯仰力矩系数影响的导数，其方法是基于扩展的纵向静稳定性和操纵性理论以及相应的飞行试验数据表示方法，对于大多数气导数，假设飞行测量数据是以平衡曲线和ＣＬ（α）曲线形式给出的，导数估计表达式是由已知常数，直接或非直接测量量的代数关系式构成，估计精度用结果标准误差或系统误差带表示。所研     

平／钝头外形俯仰动态特性对落点姿态的影响分析

采用自由振荡法数值模拟了平头、钝头外形的超声速俯仰振荡时间历程,并应用奇异分解线性最小二乘法辨识稳定性导数,得到动导数随马赫数和攻角非线性变化的规律.平头外形在较高马赫数和中小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的;小钝头外形则没有动不稳定现象.蒙特卡洛弹道仿真表明,平头外形初始时刻俯仰角的扰动会经历一个先增长然后衰减的过程,最终落地时俯仰角小于1°,满足设计要求,但其落地速度有所降低,可能影响该外形对靶标的侵彻效果.