基于本征方程的柔性机翼非线性颤振分析

柔性机翼在气动载荷作用下常常会产生较大的变形,颤振特性会随之发生变化,针对此问题线性理论常常难以进行合理的预测。以几何精确本征梁模型建立了机翼的运动方程,耦合ONERA-EDlin非线性气动模型,建立了柔性机翼的非线性气动弹性分析模型。利用Newton-Raphson和Backward-Differentiation-Formula(BDF)分别求解机翼的静态变形和动态响应,基于机翼平衡位置附近的线性化方程来判断系统的稳定性,进而确定颤振临界速度。通过算例验证了模型的准确性,并分析了不同刚度、后掠角、机翼安装角等参数对颤振速度的影响。     

基于灰色补偿反步法的机动指令跟踪控制

基于三自由度飞行动力学模型和运动学模型,提出了一种带有灰色补偿的反步(Back-stepping)控制方法,以实现某些预定的战术机动动作。首先,将飞行动力学模型和运动学模型分解成3个轴方向的相对简单的子系统;然后,针对3个子系统采用反步法分别设计出模型中没有干扰时的控制律,并采用GM(0,N)预测模型对不确定部分的模型参数进行辨识,进一步根据估计出的参数设计灰色补偿控制律;最后,根据3个子系统设计出的控制律解算出同时满足3个子系统渐进稳定的控制输入,这样使得整个系统渐进稳定。仿真结果表明,GM(0,N)预测模型能精确预测不确定模型参数,该算法能够高效实现预定的战术机动动作。     

三支点柔性转子系统支承不同心激励特征及振动响应分析

针对航空发动机三支点柔性转子系统的支承不同心问题,充分考虑转子结构特征和载荷特征,首次将当量刚度引入多支点柔性转子不同心问题的动力学分析,定量描述转子系统各支承间不同心度带来的转子轴段刚度非线性,并提出了多跨度柔性转子系统支承不同心激励的数学描述,建立了不同心激励下多跨度柔性转子系统的力学模型。基于Lagrange能量法,给出了转子系统动力学方程的求解方法,研究得到了支承不同心转子系统的动力响应特征。结果表明:支承不同心不仅引起转子过渡轴的刚度非线性,产生2倍频激励,还会给转子系统带来附加不平衡激励;对于三支点柔性转子系统而言,2倍频分量同样是支承不同心下转子系统振动响应的典型特征之一。转子系统2倍频分量随不同心量的增加而迅速增加,而1倍频分量基本保持不变。同时转子振动响应呈现"缓增速降"趋势,且随非线性刚度、不平衡量的增大愈加明显。     

基于Gerber模型的DFR法与结构细节效应

在实际工程应用中,基于Goodman模型估算结构细节疲劳额定强度(DFR)的方法偏于保守,不够经济.本工作推导基于Gerber模型的DFR值D(应力比R=0.06)的计算公式,引入细节效应系数以考虑多细节对置信度系数和D估算结果的影响;与已有结果的比较发现:基于Gerber模型的D的估算方法能够较大程度地发挥材料的潜能,降低研发成本.利用构件细节额定值系数法,通过单个细节结构的D估算了含多个细节的结构的D.估算结果与试验结果的比较验证了基于Gerber模型的DFR法的准确性和适用性.     

基于不可逆热力学的宏细观粘弹性损伤本构模型

为建立复合固体推进剂的损伤本构模型,基于不可逆热力学和叠加原理,通过引入宏观损伤效应张量,推导出一个通过有效应力表征损伤的蠕变型损伤本构方程。假设材料为初始各向同性,进一步引入细观标量损伤效应函数表征材料对称性的改变,进而得到由细观损伤效应函数表征宏观损伤效应张量的一般表达式。通过选取合适的细观损伤效应函数,文中建立的本构方程可以用于表征材料的正交各向异性损伤、横观各向同性损伤和各向同性损伤。最后,基于Schapery粘弹性微裂纹扩展模型,选取相对微裂纹密度为损伤内变量,建立了一个由损伤热力学对偶力表征的幂律型损伤演化方程。数值结果表明,建立的模型能够较好地反映材料损伤的率相关性和温度依赖性,具有良好的预测精度。     

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

基于连续损伤力学的轮盘超转破裂行为预测与验证

为有效预测轮盘的超转破裂行为,建立了1种基于连续损伤力学的预测方法。给出了该方法下基于光滑拉伸试验数据的直接时效GH4169合金弹塑性本构模型参数与损伤模型参数确定方法。对缺口半径分别为0.5、2和5 mm的3种缺口圆棒试样的拉伸响应曲线进行了有限元计算,并与试验结果进行了对比。结果表明:该模型可以较准确地预测缺口试样的极限拉伸载荷和拉伸断裂变形。采用弧长法对模拟盘的超转破裂过程进行了非线性有限元计算,使用模拟盘超转破裂试验进行了验证。该模型不仅可以较准确地预测模拟盘的破裂转速,还能有效预测模拟盘的超转破裂模式。对于实际轮盘破裂转速预测具有一定的指导意义。     

基于在线滚动序列核极限学习机的涡轴发动机非线性模型预测控制

针对涡轴发动机控制系统设计,提出了1种基于在线滚动序列核极限学习机的非线性模型预测控制方法。综合考虑直升机旋翼扭矩、燃气涡轮转速、动力涡轮转速、涡轮级间温度和压气机喘振裕度等信息,设计具有较好实时性、精度和泛化能力的多输出在线滚动序列核极限学习机作为预测模型,引入预测模型输出与发动机输出的误差进行反馈校正,利用序列二次规化算法在线求解包含限制约束的预测控制问题。在某型直升机/涡轴发动机综合平台的仿真环境中进行了直升机大幅度机动飞行仿真验证,结果表明:该模型预测控制器相比于传统串级控制具有更好的控制品质,可显著降低动力涡轮转速超调/下垂量。     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。