理论解用于复杂构型飞行器气动阻力方法研究

针对现有的理论解方法尚不适用于复杂构型飞行器气动力求解的问题,从稀薄气体分子的动力学模型出发,制作等效平面以分析各部分来流的宏观特性,并将入射流和反射流作用分开进行考虑,从而计算出几种典型情况下的气动阻力值。后续验证表明,在分子-壁面间多次碰撞不显著时,此方法能够获得较准确的气动力结果。所提出的气动力计算方法对低轨道飞行器的气动力计算及气动构型优化等工作均具有一定的参考价值。     

直升机飞行力学线性变参数建模方法研究

针对直升机飞行力学建模问题,分别基于模型缝合及模糊逼近提出了两种线性变参数(LPV)建模方法。在模型缝合方法中,基于经过滤波加权的飞行速度进行气动导数和操纵导数的插值,完成气动力的实时计算,结合非线性欧拉运动方程完成LPV模型的建立。在模型逼近方法中,建立了基于T-S模糊模型的非线性状态空间逼近模型,采用anfis结构完成了模糊模型训练。利用UH-60直升机非线性模型生成了不同飞行速度的线性状态空间模型,并分别通过两种方法建立了两套LPV模型。最后,基于频率响应和泰尔不等系数对典型LPV建模参数对模型精度的影响特性进行了分析。     

多故障转子-滚动轴承系统的动力学特性

为深入研究多故障转子-滚动轴承系统的机理,充分考虑了旋转部件的扭矩、剪切力、转动惯量及陀螺力矩等效应,采用Timoshenko梁单元和转盘单元建立了转子的有限元模型,集成了滚动轴承动力学模型,及不平衡、不对中、碰摩、裂纹及松动等故障模型,得到了多故障转子-滚动轴承系统的动力学方程。采用变步长龙格库塔法,获得了转子-滚动轴承系统的位移响应,研究了系统位移响应随转速变化的分岔图、庞加莱截面,以及故障特征处的时域波形图、轴心轨迹图、相平面图和频谱图,探究了多故障转子-滚动轴承系统的故障特征及其动力学特性和规律。仿真分析结果表明：多故障转子-滚动轴承系统在运行过程中存在周期运动,拟周期运动和复杂的类混沌运动等丰富的非线性行为,且在各非线性行为转速区间内,会呈现不平衡、不对中、碰摩、裂纹和松动等单一故障或耦合故障的特征和现象。     

基于本征正交分解的叶片碰摩系统降阶方法

针对叶片碰摩响应求解问题,提出了采用本征正交分解进行动力学降阶的方法。通过对快照矩阵进行本征正交分解生成投影子空间,将系统动力学方程投影到子空间进行模型降阶,并结合数值积分方法进行了碰摩响应求解。基于降阶模型分析了不同转速及侵入量参数下的叶片碰摩响应,并与全阶模型进行了对比。结果表明:降阶模型的时域响应幅值偏差小于5%,计算效率了提升98.4%;通过改变叶片转速、侵入量参数验证了降阶模型的鲁棒性,并且发现随着转速、侵入量的增加,本征正交模态能量在低阶与高阶之间发生转移,并呈现出不同的传递规律,由转速引起的模态能量转移与结构的固有频率存在一定关系。该方法及结论可为叶片碰摩分析及故障诊断提供依据。      

基于叶端定时的转子叶片动应变重构不确定性量化

基于叶片非接触式动应变重构理论,开展动应变重构不确定性量化方法研究。基于方差合成定理建立重构叶片动应变不确定性量化分析模型;以模拟转子叶片为研究对象,开展旋转叶片叶端定时试验,利用周向傅里叶算法获取不同叶端定时传感器布局下的测点振幅,通过最大熵方法拟合振幅分布概率密度函数,确定叶端定时测振的不确定性参数;结合Kriging代理模型和试验共振频率数据对叶片有限元模型进行修正,获取关键测点位移-应变转换因子,并获取考虑共振转速以及测点位置不确定性的转换因子不确定性参数;获取重构动应变的均值、标准不确定度和包含区间,与应变片测量数据作对比。结果表明,除5号叶片A测点外,测量动应变均位于重构动应变的95%置信度下的包含区间内,且所有叶片的应变片测点动应变重构误差不超过15%。     

基于Tanaka-Mura模型的电子束焊接头疲劳裂纹萌生模拟

为揭示高温合金电子束焊接头的疲劳特性,对其开展了疲劳裂纹萌生数值模拟研究。考虑焊缝区微观组织特性,对Voronoi图法进行改进,建立了焊缝区包含柱状晶、细等轴晶及粗等轴晶的混合晶区微观组织模型;对ABAQUS进行二次开发,考虑晶粒随机取向,生成晶粒多滑移带模型。基于Tanaka-Mura位错滑移模型,编写了疲劳裂纹萌生算法,考虑晶界处裂纹的连接与合并,对算法进行了改进,并结合有限元计算建立了电子束焊接头疲劳裂纹萌生数值模拟方法。基于上述方法对GH4169电子束焊接头不同载荷大小的疲劳裂纹萌生进行数值模拟,分析了裂纹萌生过程及萌生寿命,并与试验结果进行对比验证;还探讨了不同热影响区晶粒尺寸对焊接接头疲劳裂纹萌生的影响规律。结果表明,电子束焊接头疲劳裂纹均萌生于热影响区,但随着载荷水平的提高,萌生位置向熔合区一侧靠近;当热影响区晶粒尺寸与母材区晶粒尺寸越接近时,接头疲劳寿命越长。     

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

2.5维编织树脂基复合材料平板多尺度动力学模型修正

在2.5维编织树脂基复合材料多尺度固有振动分析的基础上,开展基于实测结构模态数据的复合材料平板动力学模型修正研究,以获取精准反映其动力学特性的模型。首先,建立复合材料平板固有频率对弹性参数的灵敏度分析方法,分析了平板固有振动特性的主要影响因素;其次,建立基于灵敏度分析的复合材料多尺度动力学模型修正方法,形成了基于MSC.NASTRAN平台的模型修正程序。最后,利用实测结构模态测试数据,实现对复合材料平板动力学模型的修正。结果表明,对关键弹性参数进行修正后,不同边界条件和纱线走向的复合材料平板试验件固有频率实测数据与仿真数据的最大误差由10.44%下降至2.39%,相关性得到了大幅改善。此外,所提出的自由模态测试方案结合动力学模型修正方法,为复合材料等效弹性参数的试验辨识提供了新的技术途径。     

基于显式预测控制的涡扇发动机控制器设计

模型预测控制因其能够较好地处理涡扇发动机约束问题且实现多变量控制而具有较大的应用潜力。为了解决传统模型预测控制算法在实际应用中存在运算量大、实时性较差的问题,基于显式模型预测设计了航空发动机多变量控制器。在控制结构上采用多速率双闭环系统,保证了控制精度;在计算上通过多参数规划将在线优化问题转化为线性函数计算问题,大幅度减少了计算量;在发动机从开环切闭环的过程提出一种增量式切换方法,实现无扰切换。数值仿真和硬件在环仿真结果表明,转速和压比稳态误差分别不超过±0.25%和±1%,控制器在25ms控制周期内能完成计算,满足嵌入式系统实时性要求。     

亚格子模型对泵喷推进器宽带非定常力预测的影响研究

为了研究亚格子模型对泵喷推进器非定常流动与宽带非定常力预报结果的影响,采用分块结构化网格建立了模型尺度下艇后泵喷推进器的计算模型,并进行了大涡模拟数值计算。从艇尾非定常流场特征量和推进器转子脉动载荷两个方面对比了三种不同亚格子模型计算结果的差异,并分析了泵喷内部流动与转子非定常力间的内在联系。研究结果表明：三种亚格子模型均能得到含有叶频宽带峰的轴向推力谱,且整体趋势相近。泵喷转子上游的湍流强度和尺度的分布对亚格子模型较敏感,其中Smargrinsky-Lilly模型得到的湍流强度较强,尺度较大,该模型下的湍流谱在非平衡区量级较大,但由低频向高频的衰减较快,并且预测到的分离区范围大,导叶尾缘脱落涡分散,导致转子上游来流空间分布不均程度较强。对于转子叶片上的载荷脉动,Smargrinsky-Lilly模型预测的推力谱中线谱成分明显,并且叶频处宽带谱峰“陡峭”,而WALE模型和KET模型的结果宽带特性较强,对比非定常推力测试结果可知,WALE模型和KET模型更适于该问题宽带非定常力预测。