Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

民用飞机全机静力试验前机身工况配平载荷设计

民用飞机全机静力试验中由于加载方案复杂,为了简化试验加载数量,传统的配平方案主要保证考核区的加载准确性,对于非考核区的载荷简化到若干加载点上。对于飞机机身这种薄壁结构,非考核区的配平对考核区影响较大。以某飞机前机身严重工况为例,基于飞机全机有限元模型,研究静力试验的载荷配平方案;通过对比分析态载荷与传统试验态配平载荷方案下的内力解,对传统全机静力试验载荷配平方案进行优化。结果表明：优化后的试验态内力解与分析态内力解基本吻合;对于飞机机身薄壁结构,圣维南原理并不适用,全机静力试验的载荷配平要与真实工况受力情况吻合,才能达到试验验证目的。     

一种基于随机优化算法的民用飞机突风动载荷重建方法

民用飞机突风载荷的准确重建对于飞行安全具有十分重要的意义。基于随机优化算法,提出一种突风动载荷重建方法。首先,基于频域气弹耦合方程,建立一个突风动载荷时域分析模型,其次,引入Hicks-Henne 型函数进行突风速度场的参数化建模,在此基础上,基于随机优化算法对各型函数参数进行迭代寻优,使得优化目标值与测量值一致,最终得到重建后的突风速度及对应的动态响应及载荷。最后,以GTA 模型经历1-cos 形状突风为例进行数值仿真。结果表明：本文提出的突风动载荷重建方法可有效实现突风速度场的识别,进而完成突风动载荷的重建;针对本文算例,过载最大误差为3.5%,载荷最大误差不超过1%。     

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

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

航空发动机低压转子系统的“可容模态”设计及实验验证

建立了带弹性支承和阻尼器的航空发动机低压柔性转子动力学模型,进行了模态计算以及临界转速处的响应计算。在考虑临界转速约束与“临界跟随”现象约束的条件下,综合模态不平衡影响因子、弹支应变能占比以及套齿连接结构稳定性构造了可容度评价函数,建立了低压转子系统的“可容模态”优化设计方法。设计并搭建了低压转子实验系统,从模态测试实验、阻尼器减振实验以及长时间“共振”实验,验证了设计方法的可靠性。研究结果为计算的临界转速与实际测量的临界转速最大误差为3.86%,阻尼器在1阶与2阶临界转速处的减振比最大可达45.6%,实验转子系统在1阶与2阶临界转速处各完成了长达412.5 s和429.8 s的“共振”实验,共振过程中转子系统各通道振动单峰值稳定在100μm以内,且无次谐波产生。表明了所建立的航空发动机低压转子系统“可容模态”优化设计方法是可行的。     

C/C复合材料喷管烧蚀壁面退移流固耦合仿真研究

为了保障固体火箭发动机C/C喷管的可靠性,建立了一套正确反映发动机喷管烧蚀过程的流固耦合计算模型,以实现对喷管烧蚀率的高精度预估。依据热化学烧蚀理论以及喷管内燃气与喷管结构体界面的质量平衡和能量平衡关系,建立并验证了考虑壁面退移的C/C喷管流固耦合方法,实现了燃气流动、异相化学反应、结构体传热三者间的耦合。通过实验发动机喷管的烧蚀计算,论证了模型的正确性,并分析了不同金属铝含量对烧蚀率的影响,计算所得的烧蚀率与实验值最大相对误差为4.3%,与不考虑壁面退移的耦合算法计算结果对比,计算精度最高可提升46%。计算结果表明：C/C喷管在喉部附近烧蚀最为严重;推进剂中Al含量的增加导致燃气中氧化组分浓度降低,进而减少了烧蚀速率,这些结论与C/C喷管烧蚀相关研究结果一致。     

温度参数下轴编C/C复合材料的超高温热结构性能分析

为了研究轴编碳/碳（C/C）复合材料的超高温热结构性能,开展了2800℃超高温载荷下复合材料的拉伸实验,表征了复合材料组分材料的高温失效形貌;提出了立方体单胞和正六棱柱单胞的温度周期性边界条件的算法和预测复合材料等效热结构参数的方法,在材料属性各向同性正六棱柱单胞模型上验证了方法的正确性;最后建立了轴编C/C复合材料的代表性体积单元,计算了复合材料随温度载荷变化下的等效热结构参数,与实验数据进行了对比分析。分析表明：开展的轴编C/C复合材料的超高温拉伸实验,可以得到复合材料在2800℃载荷下复合材料的等效刚度值,同时分析了轴向试件和径向试件的高温失效机理,径向试件的纤维束拉伸强度对复合材料的径向拉伸强度起到决定性的作用。在应用位移以及温度周期性边界条件的基础上,提出的方法可以得到随温度载荷变化下复合材料的等效热结构参数,得出复合材料在横向满足近似各向同性的特性,预测结果表明只改变温度参数下复合材料的等效热导率不会发生改变。提出的方法同样也适用于其他编织复合材料的热结构性能研究。     

三维五向编织复合材料的冲击压缩特性及破坏机制

利用分离式霍普金森压杆(split Hopkinson pressure bar,SHPB)实验系统和高速摄像技术对三维五向碳/环氧编织复合材料的动态压缩特性进行研究。通过对编织角为22.3°的试样分别进行沿纵向和横向方向的冲击压缩实验,得到材料在200~1200 s^-1应变率范围内的应力应变曲线,并结合高速摄像记录的动态压缩过程,对不同应变率下材料在高速变形下的渐进破坏规律进行分析。同时,综合试样的宏观破坏特征和微观断口形貌特征,进一步分析材料的破坏模式及破坏机理。结果表明:随着应变率的增加,材料在纵向和横向均具有一定的应变率强化效应,在横向方向的应变率强化效应更为显著;不同加载方向下材料的渐进破坏过程、应力应变曲线特征以及破坏方式均具有明显差异,且随着应变率的增加而发生改变。     

螺旋桨向量平衡

螺旋桨旋转时各桨叶的离心力可分解为在桨叶轴线旋转平面内的静不平衡力和垂直于此旋转平面的动不平衡力偶。由于各桨叶不完全相同,这两种不平衡度在不同矩角位置会有差异。但所加静平衡力和动平衡力偶则不随变矩而变化,它们的平衡作用在所有矩角位置都是一样的,可以采用本文所述的向量分析方法,在特制的向量平衡图上预测对各矩角平衡的结果。这样就可以只需一次平衡即使各矩角都平衡合格并达到最佳平衡状态,而用常规平衡方法需要反复凑试甚至无法合格。     

副翼配重的着陆冲击响应

提供了一个计算副翼配重着陆冲击响应的方法。由于主起落架为纯摇臂式,考虑机体的弹性,机体和主起落架运动联立求解的加速度不是显式,因此必须用高斯方法解加速度的线性方程组,再用吉尔方法解微分方程。