一种螺旋桨发动机颤振工程计算方法与流程

本发明属于航空气弹分析技术，涉及螺旋桨发动机系统颤振分析方法。
背景技术：
：装有螺旋桨发动机的飞机，存在螺旋桨与发动机安装结构组合的螺旋颤振稳定性问题。对该问题的研究方法常采用螺旋桨－发动机系统模型，考虑俯仰和偏航两个自由度，计算螺旋桨气动力和陀螺力矩，最后推导出螺旋颤振的动力学方程，采用数值方法求解复矩阵广义特征值问题，获得颤振稳定性计算结果。参见论文：金一忱，《螺旋颤振分析》，航空学报，1993年4月第14卷第4期；姚一龙、施辉，《涡桨发动机飞机的螺旋颤振分析》，民用飞机设计与研究，1996年第4期。技术实现要素：本发明的目的是：现有方法计算螺旋颤振分析主要采用数值计算方法，关键的螺旋桨气动力计算和颤振方程求解采用设计人员个人自编程序进行计算。计算的准确度因人而异，颤振计算程序不对外发布，不利于计算方法的推广应用。另外螺旋颤振计算需要考虑螺旋桨几何参数、转速、飞行速度、飞行高度、空气密度和发动机系统质量、刚度这些参数，并进行参数的组合变化，计算状态较多，计算量大，自编程序不适合工程应用。本发明方法采用国际著名软件msc.nastran进行计算，计算精度可靠，计算方便，计算模型可视化，利于计算方法推广应用。本发明的技术方案是：一种螺旋桨发动机颤振工程计算方法，包括以下步骤：1.采用有限元方法计算螺旋桨－发动机系统质量矩阵和刚度矩阵(1)螺旋桨－发动机系统有限元模型简化方法：飞机螺旋桨－发动机系统可简化概括为三部分结构，一是螺旋桨，二是发动机主结构，三是发动机舱。结构上螺旋桨和发动机通过旋转轴连接在一起，发动机再安装在发动机舱上，发动机舱最后与机翼连接。螺旋桨和发动机简化为集中质量单元，发动机舱结构简化为杆、板、梁单元。螺旋桨和发动机的质量和转动惯量采用计算或者测量的方法可以得到。(2)创建螺旋桨－发动机系统有限元模型：将发动机舱三维几何模型导入msc.nastran有限元软件。发动机舱由杆、板、梁单元组成。螺旋桨和发动机简化为集中质量单元，通过rbe2单元与发动机舱模型连接，发动机安装点和发动机安装接头之间的减震器用弹簧元来模拟，模拟发动机的安装刚度。模型在发动机舱悬挂接头处和与机翼连接处约束。螺旋桨发动机系统有限元模型创建完毕。(3)计算螺旋桨－发动机系统质量矩阵和刚度矩阵采用1－2中建立的螺旋桨－发动机系统有限元模型，用msc.nastran有限元软件进行模态分析。通过模态分析可以获得系统偏航模态和俯仰模态及频率值，通过模态分析同时可以获得螺旋桨－发动机系统的质量矩阵和刚度矩阵。2.采用msc公司propf.for程序计算螺旋桨的气动力系数(1)输入数据准备：输入数据为螺旋桨几何参数、转动惯量、质量，发动机转动惯量、质量、螺旋桨转速、飞机飞行速度、空气密度等，将这些数据放入一个＊.inp文本文件，文件名称任意。(2)运行propf.for：在fortran编译环境下运行propf.for，生成propf.exe可执行程序，运行propf.exe可执行程序后，程序以dos界面形式出现，并在界面中提示输入＊.inp文本文件的文件名。输入(1)中保存的＊.inp文件名后，程序自动开始计算，计算完毕将程序自动关闭，计算结果以＊.out和aerodmig.dat文本文件保存在当前程序文件目录下。(3)气动力系数数据获取打开＊.out文本文件，从中可以读取到以下气动力系数结果：cyq——侧向力系数对俯仰速率的导数cmq——俯仰力矩系数对俯仰速率的导数cmpsi——俯仰力矩系数对偏航角的导数czpsi——垂向力系数对偏航角的导数cypsi——侧向力系数对偏航角的导数打开aerodmig.dat文件，可以获得由螺旋桨气动力pz，py，my,p，mz,p分解成的msc.nastran软件可以识别的阻尼矩阵[b2pp]和刚度矩阵[k2pp]dmig卡片。根据不同飞机飞行速度和高度，可以计算一系列阻尼矩阵[b2pp]和刚度矩阵[k2pp]dmig卡片数据。3.采用msc.nastran软件的结构化sol145流程计算颤振最后将由1计算得到的螺旋桨－发动机系统的质量矩阵和刚度矩阵数据和由2计算得到的一系列阻尼矩阵[b2pp]和刚度矩阵[k2pp]dmig卡片组合到msc.nastran软件的结构化sol145流程计算卡片中，就可以直接计算螺旋颤振了。计算结果可以绘制螺旋桨－发动机系统飞行速度随阻尼和频率的变化曲线，从中可以识别螺旋桨－发动机系统颤振临界速度，用于判别系统颤振稳定性。本发明优点是：1.螺旋桨－发动机系统计算模型采用有限元创建可视化好，动力学特性采用msc.nastran计算，计算结果精确，与试验结果符合性好；2.螺旋桨气动力和颤振计算采用国际知名软件计算，计算的螺旋桨力矩系数随速度的变化曲线与发动机厂家计算的曲线符合性好，螺旋桨－发动机系统颤振稳定性通过了飞行试验验证，计算结果可靠；3.可以一次性计算不同参数变化情况下的螺旋颤振，减少了计算时间，减轻了设计人员工作量，提高了工作效率。4.采用本发明方法计算螺旋颤振，适合工程应用。附图说明图1发动机－螺旋桨系统模型图2发动机舱俯仰模态图3发动机舱偏航模态图4力矩系数对比图5阻尼矩阵[b2pp]和刚度矩阵[k2pp]dmig卡片图6海平面v－g图图7海平面v－f图图8不同螺旋桨转速v－g图图9不同螺旋桨转速v－f图具体实施方式实施例一：飞机螺旋颤振分析1.螺旋桨－发动机系统结构动力学特性计算螺旋颤振研究建立了发动机舱、发动机、螺旋桨系统的有限元模型如图1所示，模型在发动机舱悬挂接头处和与机翼连接处简支约束。通过模态分析获得系统偏航模态和俯仰模态及频率值。螺旋桨－发动机系统俯仰和偏航模态振型见图2、图4。修正分析模型使发动机舱系统俯仰和偏航模态频率与全机地面振动试验获得的结果误差在5％以内，保证了分析模型的刚度和实际结构一致，模型可以用于螺旋颤振分析。分析模型修正后的振动模态与地面振动试验结果的对比见表1。表1发动机舱系统模态频率模态分析频率(hz)试验频率(hz)误差俯仰6.697.01－4.5％偏航9.118.852.9％2.螺旋桨空气动力系数计算msc公司编制的fortran计算程序propf.for，可以用来计算气动力力矩系数。此程序需要的输入数据为螺旋桨几何参数、转动惯量、质量，发动机转动惯量、质量等。螺旋桨叶片类型为m10876an(s)(k)，将螺旋桨叶片分割为14个剖面，螺旋桨参数见表2，每个剖面半径和宽度示例见下表3。运行propf.for程序，输入这些参数后，程序计算完毕输出可以螺旋桨叶片气动力系数，也可以输出阻尼矩阵[b2pp]和刚度矩阵[k2pp]。通过将将msc公司propf程序和螺旋桨厂家叶素理论计算的气动力系数做了对比，见图4至图5，图中programe为程序计算结果，factory为厂家计算结果。从力矩系数随速度的变化曲线可见，在飞行速度范围内两者的符合性较好。计算得到的阻尼矩阵[b2pp]和刚度矩阵[k2pp]dmig卡片见图5所示。表2螺旋桨参数叶片数转速rpm半径inch51700111.2表3叶片参数半径inch宽度inch7.853.562略略55.602.83.螺旋颤振计算螺旋颤振方程的求解采用msc/nastran软件的结构化sol145流程求解颤振方程。螺旋桨气动力pz，py，my,p，mz,p分解成阻尼矩阵[b2pp]和刚度矩阵[k2pp]，[b2pp]阵和[k2pp]阵是关于螺旋桨几何参数、转速、飞行速度、飞行高度、空气密度等参数的函数，可以用msc公司计算程序propf.for按一系列不同的飞行速度(等高度)直接生成一系列代表阻尼阵和刚度阵的dmig卡片。将它们加入到带有螺旋桨－发动机动特性数据的sol145流程控制卡中。然后合并到数据卡片中，最后用pk法求解颤振稳定性。螺旋颤振计算得到的速度－阻尼(v-g)图和速度－频率(v-ω)图的见下图所示。y12f飞机已经取得faa和caac适航证，螺旋颤振分析方法已经获得适航认可，并且y12f飞机在颤振试飞过程中未发生螺旋颤振问题，表明计算结果可信。当前第1页12