发动机静子尾流激励对叶盘结构高阶共振可激性的研究

针对航空发动机中转静干涉流场诱发的叶盘结构强迫响应问题，分析由导流叶片和动叶构成的转静干涉流场的尾流激励所激起的叶盘结构高阶共振的可能性和危险性.基于非定常转静干涉流场数值模拟结果，通过时域与空间傅里叶分析得到叶盘所受气动压力的频谱和空间波谱，同时定义投影于叶盘结构模态空间的激励（即模态力）为尾流激励的可激性，用于对叶盘结构共振响应强度的评估.研究结果表明:尾流激励对非主要共振点的可激性远远小于其对主要共振点的可激性；在发动机工作转速范围内允许那些模态位移场与激振力场的夹角较大（余角较小）的共振存在，即尾流激励可激性很小的高阶共振的存在是合理的.      

临近空间太阳能飞行器横航向稳定性

临近空间太阳能飞行器的横航向气动导数和质量特性与常规飞行器存在显著差异,其横航向模态具有不同于常规飞行器的特点。理论分析表明,临近空间太阳能飞行器滚转交感动导数较大,偏航阻尼动导数较小,导致螺旋模态发散;航向气动阻尼力矩和惯性力矩之比较大,因而荷兰滚模态阻尼比较大。为解决常规线性化动力学模型无法体现风场影响的问题,以地速在机体坐标系中的投影作为状态变量,建立了考虑风场影响、适用于地速为0 m/s状态的线性化横航向动力学模型。利用此模型分析了临近空间太阳能飞行器在水平风和垂直风中的横航向稳定性。分析结果表明,稳定风场对横航向模态特征根无影响,但会导致横航向模态特征矢量发生改变。建立的动力学模型可用于此类飞行器的动力学分析和仿真。     

基于CFD动网格技术的飞艇动导数计算方法

由于飞艇和飞机这两种飞行器在飞行原理、气动外形等方面都存在较大的差别,因此现有的飞机动导数工程估算公式和经验数据并不适用于飞艇动导数求解。为了解决上述问题,提出了一种新的基于CFD动网格技术的飞艇动导数计算方法。通过非结构化动网格技术,建立飞艇在流场中的强迫振动模型,采用有限体积法对N—S方程求解,研究飞艇在强迫振动中的气动特性变化规律,并基于数值计算结果求解飞艇的动导数。研究发现,得到的计算结果与风洞试验结果基本一致,表明基于CFD动网格技术进行飞艇的动导数求解是有效可行的。     

扑翼飞行器机翼平面形状设计方法研究

通过均匀试验法在设计变量状态空间内寻找分布更均匀、更能反映系统性能的状态点,使用改进的片条理论计算机翼的气动特性,建立了扑翼飞行器机翼面积分布(尖削比)、扑动迎角与扑动频率的响应面模型,得到了尖削比对机翼气动特性的影响规律和机翼升力系数关于尖削比、频率与迎角的估算公式。     

炮舱结构段的动特性与动响应分析

对B型机原炮舱结构进行了动特性与动响应分析,并以此为基础,对其结构进行了加装阻尼材料的最佳设计。本计算采用了空间迭代法求结构动特性,采用模态迭加法求结构动响应。经试验验证,计算结果与试验结果基本吻合。     

中介传动齿轮的耦联共振

 本文研究了某减速器的振动故障,发现并研究了其中介传动齿轮的耦联共振。用试验模态分析结果、模态综合技术、瑞利法相结合的办法,提高了试验模态分析所得振型的精度。对多种结构修改方案计算求得了它们的各阶动频,从中选定了最佳的排振修改方案。     

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

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

风切变中飞行能量高度分析及反馈

本文采用模态集结法,导出了能量高度响应的近似公式并分析了高度模态与能量高度的对应关系,由此给出了采用能量高度反馈的设计方法。计算表明近似公式反映了飞机在穿越低空风切变时的响应特性,模拟结果表明本文中的反馈设计明显提高了飞机穿越低空风切变的能力。     

基于响应传递比的转子叶片动应变重构

提出基于响应传递比的转子叶片动应变重构方法,从系统应变频响特性出发,在频域获得应变-应变响应传递比关于模态振型的解析表达式,建立叶片已知测点应变与不可测位置应变的映射关系;开展高速旋转叶片动应变测量实验,建立叶片有限元模型,并采用转子叶片前4阶固有频率的实验结果修正有限元模型;开展考虑转速影响的转子叶片模态分析,提取对应的应变模态振型,计算得到单模态振动下应变-应变响应传递比;综合应变响应传递比和实测动应变重构叶片不可测位置的动应变。实验结果表明:与应变片实测结果相比,基于响应传递比的转子叶片动应变重构结果相对误差小于10%。      

非定常超音速侧滑机翼的升力、力矩系数及动导数的数值计算

 本文利用分布在机翼平面上的非定常超音速点源的特征线网格法给出非定常超音速侧滑机翼的数值解。本文给出侧滑机翼作沉浮、俯仰、滚转谐振时的振型以及计算升力、力矩和动导数（包括交叉导数）的计算方法。按本文用非定常理论数值解对典型的无侧滑机翼计算出的动导数与常用的按准定常理论解析公式计算结果进行比较和鉴定,结果表明本文方法是令人满意的。     

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

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

非定常耦合流动对叶片动力响应的影响

针对"在轴流压气机中实现非定常耦合流动状态时,是否会引发下游转子产生结构振动问题"进行了研究.在已有的三维非定常CFD(计算流体动力学)计算频谱分析结果的基础上,使用有限元方法进行模态分析,通过调整上游静叶数的手段,研究实现耦合流动后产生的动力响应问题,并对调整前后压气机的气动特性进行比较.结果证明,通过适当的调整上游静子叶片数,可以在最大化发掘非定常耦合流动的潜力的同时避开转子的共振频率.      

耦合六自由度运动弹性飞机阵风响应数值模拟

耦合飞机刚体六自由度运动,考虑飞机结构气动弹性变形,基于全湍流Navier—Stokes方程流场数值模拟方法．建立了离散阵风作用下弹性飞行器气动载荷与飞行姿态响应的数值模拟技术。其中模态空间中结构动力学方程以及六自由度运动方程分别采用四阶R—K进行时间推进求解,采用RadialBasicFunction（RBF）技术进行气动／结构数据耦合。对应飞机姿态以及弹性变形采用RBF与TFI组合模式的动网格方法进行网格更新;以飞翼布局弹性飞机为数值研究对象,研究了离散阵风作用下飞行器的气动载荷响应以及飞行姿态的变化,对比分析了刚性飞行器与弹性飞行器对离散阵风响应的区别,为进一步系统地分析真实飞机飞行品质提供数值平台。     

Helicopter aeroelastic analysis with spatially uncertain rotor blade properties

This paper investigates the effects of spatially uncertain material properties on the aeroelastic response predictions (e.g., rotating frequencies, vibratory loads, etc.) of composite helicopter rotor. Initially, the spatial uncertainty is modeled as discrete random variables along the blade span and uncertainty analysis is performed with direct Monte Carlo simulations (MCS). Uncertainty effects on the rotating frequencies vary with the higher-order modes in a non-linear way. Each modal frequency is found to be more sensitive to the uncertainty at certain sections of the rotor blade than uncertainty at other sections. Uncertainty effects on the vibratory hub load predictions are studied in the next stage. To reduce the computational expense of stochastic aeroelastic analysis, a high-dimensional model representation (HDMR) method is developed to approximate the aeroelastic response as functions of blade stiffness properties which are modeled as random fields. Karhunen–Loève expansion and a lower-order expansion are used to represent the input and outputs, respectively, in the HDMR formulation which is similar to the spectral stochastic finite element method. The proposed method involves the approximation of the system response with lower-dimensional HDMR, the response surface generation of HDMR component functions, and Monte Carlo simulation. The proposed approach decouples the computationally expensive aeroelastic simulations and uncertainty analysis. MCS, performed with computationally less expensive HDMR models, shows that spatial uncertainty has considerable influence on the vibratory hub load predictions.     

阻尼不相似动力学模型的时域响应修正方法

 动力学相似的缩比模型试验中,阻尼往往难以保证相似,使得缩比模型的动响应测试数据难以直接换算到结构原型上,因此提出一种阻尼不相似动力学模型的时域响应修正方法来解决这一问题。该修正方法假定缩比模型的阻尼可用比例阻尼模型近似表征,并直接从实际缩比模型与理想缩比模型的一般强迫振动响应计算方程出发,利用线性系统的叠加原理和模态叠加法,将动响应修正量的求解转换为理想模型在一个等效附加激励力作用下的时域响应求解。同时,针对实际工程中响应测点数目有限的问题,利用模态缩聚法进行了未测量点的响应反演。该修正方法仅需已知结构原型和实际缩比模型的模态阻尼比,以及准确建立的实际缩比模型的质量矩阵和刚度矩阵,即可实现实际缩比模型在任意激励工况下的测试动响应修正。以某型飞机的翼梁缩比模型为研究对象,对所提出的响应修正方法进行了验证。试验和计算结果对比分析表明,修正后的响应功率谱密度(PSD)和响应的均方根(RMS)值与理想模型基本一致,表明了本文方法的可行性和有效性。     

A novel stability prediction approach for thin-walled component milling considering material removing process

The milling stability of thin-walled components is an important issue in the aviation manufacturing industry,which greatly limits the removal rate of a workpiece.However,for a thin-walled workpiece,the dynamic characteristics vary at different positions.In addition,the removed part also has influence on determining the modal parameters of the workpiece.Thus,the milling stability is also time-variant.In this work,in order to investigate the time variation of a workpiece's dynamic characteristics,a new computational model is firstly derived by dividing the workpiece into a removed part and a remaining part with the Ritz method.Then,an updated frequency response function is obtained by Lagrange's equation and the corresponding modal parameters are extracted.Finally,multi-mode stability lobes are plotted by the different quadrature method and its accuracy is verified by experiments.The proposed method improves the computational efficiency to predict the time-varying characteristics of a thin-walled workpiece.     

基于系统响应特征的指尖密封泄漏特性分析

结合指尖密封动态工作特点,通过研究转速和转子不平衡力与转子跳动之间的关联规律,得到转子的位移激励,构建了指尖密封系统的动态计算模型,依据指尖密封的系统响应特征研究其泄漏间隙的获取办法,进而基于泄漏间隙特征建立了指尖密封动态泄漏量计算方法.针对某型发动机转子进行了动态泄漏量分析,结果表明:在转子每个转动周期内,指尖密封的动态响应跟随转子激励周期性变化,由此产生的动态迟滞泄漏间隙也随时间而周期性变化,迟滞泄漏的大小随着密封上下游压差的增加而增大,随着转子转速的升高先增加后减小;设置一定的装配过盈量能够一定程度地减小迟滞泄漏;指尖密封磨损后迟滞泄漏规律表现出先减小后增加再降低的变化趋势;与国内外试验结果的对比分析初步验证了提出的计算方法,为指尖密封动态性能设计方法研究提供了参考.      

基于ERSM 的涡轮机匣径向变形动态概率分析

针对复杂结构动态概率分析,提出了高效率动态概率分析的极值响应面方法(ERSM)。首先介绍了动态概率分析的极值响应面方法基本原理,并推导出了其数学模型;其次,考虑机匣材料属性与边界条件的非线性、热载荷动态性以及多个随机变量,将该方法应用到高压涡轮机匣径向变形概率分析中,得出了各输入输出随机变量的分布特征和影响机匣径向动态变形的主要因素。与Monte Carlo法和响应面法进行比较,结果表明:ERSM在保证计算精度的前提下,既能大大缩短计算时间,又能提高计算效率,同时也验证了该方法的可行性和有效性。     

多目标自然层流翼型反设计方法

 进行了基于扰动放大N-因子的目标压力分布设计方法的多目标自然层流（NLF）翼型反设计方法研究。流场分析和转捩位置计算用XFOIL程序,大大减少了计算花费。用N-因子设计方法进行有NLF范围要求和满足气动约束的目标压力计算,压力恢复段的压力分布用Stratford分离准则来进行设计。用基于响应面方法的优化方法来进行反设计计算,使用不含二阶交叉项的二次多项式模型的响应面模型,大大减少了构造模型所需的试验次数;设计空间内试验点的选取满足D-优化准则。根据设计目标的设计状态,进行了多目标翼型反设计。计算结果表明,设计结果的层流范围和设计目标基本吻合,该方法可以用于NLF翼型的多目标反设计中。     

热结构动力优化设计

 考虑温度沿飞行器表面的分布及材料物性系数随温度的变化,按Ｋｕｈｎ－Ｔｕｃｋｅｒ条件建立迭代式,研究具有频率约束、振型节点约束以及变量上、下限约束的热结构动力优化设计。算例讨论了前３阶固有特性及最小目标随温度场的变化规律,证实本方法有效、可行。