叶轮机械失速颤振的全三维数值分析方法

给出了建立在全三维N-S方程求解基础上的叶轮机械失速颤振分析的数值方法。包括用定常N-S方程求解三维粘性失速流和用三维N-S振幅方程求解三维粘性失速振荡流。为了验证方法的准确性, 用压气机环形叶栅的失速流实验和振荡叶片振荡流场的实验作了对比, 同时应用本方法分析了某30万千瓦汽轮机末级叶片的失速颤振。      

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

跨音弯掠风扇转子叶片的气动弹性剪裁

跨音风扇进口级增压比的进一步提高,主要受两方面的制约,一是效率,二是气动弹性失稳。这两个因素与转子流场中的激波和激波诱导的大尺度分离紧密相关。激波结构又和转子叶片前缘空间曲线形状直接联系。因此,如何精心设计转子叶片前缘空间曲线形状来控制激波结构,就成了当代风扇气动力学的前沿,并导致风扇弯掠空气动力学概念的出现。在具体设计中,此问题可称之为气动与气动弹性综合剪裁。本文简述综合剪裁中的一个组成部分,即气动弹性剪裁问题,并对一个单级风扇转子叶片给出实例。     

考虑顶隙的压气机单转子内三维跨音速粘性流场结构的数值模拟

借助于考虑叶片顶部间隙的压气机转子内跨音速湍流流场的三维数值解,得到了该转子流场中的复杂涡系与波系结构。给出了叶片表面及轮毂面上的壁面流型,分析并揭示了转子叶道内的马蹄涡、通道涡、角涡、尾涡等复杂涡系及其与激波的相互干涉;特别是,分析并揭示了转子叶道内顶隙射流的特点,得到了顶隙附近激波与泄漏流动的干涉形态及泄漏流动的发展,并着重给出了其与激波的干涉形态。     

压电复合材料壁板颤振的控制

针对考虑热效应的壁板颤振问题,利用von Karman大变形应变-位移关系、气动力活塞理论和准定常热应力理论建立了埋入压电材料的复合材料板颤振的气动弹性力学模型,使用Bogner-Fox-Schmit单元推导出壁板颤振的非线性有限元方程.将该动力学方程转换到状态空间模型,给出了基于这个非线性模型的最优控制.应用Runge-Kutta方法在时域内对一个四边简支的含压电材料的复合材料壁板颤振的控制进行了仿真,数值结果证明了该方法的有效性.     

航空航天薄壁件铣削过程加工状态监测研究进展

对切削加工状态进行精准监测,是实现航空航天薄壁件加工变形控制的重要保障。围绕航空航天薄壁件铣削加工状态监测的最新研究进展进行了评述,详细介绍了建立加工状态监测模型的关键技术与方法,包括加工信息采集处理、特征提取和特征融合。归纳了学者们在薄壁件加工过程中对刀具磨破损、铣削颤振、铣削变形等具体状态监测的研究进展。基于数字孪生技术,构建了面向薄壁件铣削加工状态监测的优化系统。最后,根据现阶段本领域发展状况对薄壁件铣削加工状态监测进行了展望。     

航发叶片机器人磨削颤振检测方法

针对机器人磨削航发叶片过程中末端执行器颤振检测不及时的问题,提出一种基于排列熵的磨削颤振检测方法。利用滑动窗口算法计算颤振原始信号的连续排列熵值,通过排列熵经验阈值0.95判断机器人末端是否发生颤振。其中,采用基于序数模式的排列熵算法,很大程度上提高了排列熵特征的提取效率。当信号长度为10000时,所提出算法的排列熵的计算时间减小到约0.25 s,比传统排列熵算法的计算时间减小了一个数量级。数值模拟及试验结果表明,所提出的方法可以在颤振爆发前约0.48s检测出颤振,为采取颤振抑制措施争取了更多的反应时间。     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

考虑光压力和热载荷的太阳帆薄膜力学行为分析

基于冯卡门板理论描述薄膜结构大变形,耦合光压力与热载荷,构建10μm量级薄膜非线性动力学主控方程。采用伽辽金离散方法和龙格-库塔推进方法数值求解方程时域解,并引入最大李雅普诺夫指数判别薄膜混沌运动。分析表明,法向光压力远大于切向光压力并抑制薄膜动力学颤振,仅计入光压力时,太阳帆膜呈静变形特征。当计入的光辐射热载荷达到临界值后,考虑切向光压力与法向光压力的薄膜都呈现动力学混沌运动。综上,光压力抑制太阳帆薄膜颤振,热载荷导致薄膜进入混沌运动。