航空发动机强度与振动-课程设计

航空发动机强度与振动课程设计报告姓 名： 肖庭文 专 业： 飞行器动力工程班 级： 080141H 指导教师： 李书明（教授）但敏二0一一年十一月

题目及要求题目基于ANSYS的叶片强度与振动分析叶片模型把叶片简化为根部固装的等截面悬臂梁。叶片模型如右图所示，相关参数如下叶片长度：0.04m叶片宽度：0.008m叶片厚度：0.002m叶根截面距旋转轴的距离为0.16m材料密度：7900kg/m 叶片振动的有限元分析（1） 叶片静频计算与分析要求：给出1到 叶片振动的有限元分析（1） 叶片静频计算与分析要求：给出1到10阶的叶片振型图，并说明其对应振动类型。（2） 叶片动频计算与分析要求：计算出叶片在转速为1000rad/s，2000rad/s,4000rad/s,8000rad/s,10000rad/s下的动频值，用表格形式表示。（3） 共振分析要求：根据前面的计算结果，做出叶片共振图（或称Campbell图），找出叶片的共振点及共振转速。 按要求撰写课程设计报告说明：网格划分必须保证结果具有一定精度。各输出结果图形必须用ANSYS的图片输出功能，不允许截图。弹性模量：2.11011Pa泊松比：0.3叶片的静力分析（1）叶片在转速为5000rad/s下的静力分析。要求：得到vonMises等效应力分布图，并对叶片应力分布进行分析说明。课程设计报告基于ANSYS的叶片强度与振动分析1.ANSYS有限元分析的一般步骤1）前处理前处理的目的是建立一个符合实际情况的结构有限元模型。在Preprocessor处理器中进行。包括：分析环境设置（指定分析工作名称、分析标题、）定义单元类型、定义实常数、定义材料属性（如线弹性材料的弹性模量、泊松比、密度）、建立几何模型（一般用自底向上建模：先定义关键点，由这些点连成线，由线组成面，再由线形成体）、对几何模型进行网格划分（分为三个步骤：赋予单元属性、指定网格划分密度、网格划分）在本课程设计中，先在Preferences中定义了所要研究的对象是structuraK结构），然后在Preprocessor中定义材料的类型为structuralsolid->Brick8node45再设定材料密度为 DENS=7900kg/m3，弹性模量为EX=2.11011P，泊松比为aPRXY=0.3。最后根据叶片在空间的摆放位置创建关键点（Keypoints），然后依次建立面（Areas）->体（V）lumes）。建体的关键是要使简化的长方体叶片模型底面中线正对X轴，并能绕其转动。K2）施加载荷、设置求解选项并求解这些工作通过Solution处理器来实现。指定分析类型（静力分析、模态分析、谐响应分析、瞬态动力分析、谱分析等）设置分析选项（不同分析类型设置不同选项，有非线性选项设置、线性设置和求解器设置）、设置载荷步选项（包括时间、子步数、载荷步、平衡迭代次数和输出控制、加载（ANSYS结构分析的载荷包括位移约束、集中力、面载荷、体载荷、惯性力、耦合场载荷，将其施加于几何模型的关键点、线、面、体上)然后求解。在本课程设计中，静力分析时要固定底面边界，施加5000rad/s绕X轴的转速；模态分析中的静频分析时要固定底面边界，设定10阶最大阶数，然后求解(solve),最后查看结果；模态分析中的动频分析时要固定底面边界，先在static分析类型中第一次求解(solve)出对应转速下的离心拉仲应力，然后再到modal分析中第二次求解(solve)出动频值，求解时要考虑离心拉仲应力的影响。(3)后处理当完成计算以后，通过后处理模块GeneralPostproc查看结果。ANSYS软件的后处理模块包括通用后处理模块( POST1)和时间历程后处理模块(POST26)。可以轻松获得求解计算结果，包括位移、温度、应变、热流等，还可以对结果进行数学运算，然后以图形或者数据列表的形式输出。结构的变形图、内力图(轴力图、弯矩图、剪力图)，各节点的位移、应力、应变，还有位移应力应变云图都可以得出，为我们分析问题提供重要依据。在本课程设计中，主要是通过后处理模块查看叶片变形的位移振动图(DOFsolution)和vonMises等效应力分布图(stress)。算出的动频值结果可以在Resultssummary中查看，另外还可以通过菜单栏中的PlotCtrls->HardCopy->ToFile...中输出白底色图片和PlotCtrls->Animate->modeshape中输出动画。2.叶片的静力分析WC«£jftL5C<LUT.'IC«H^T.E.P=2^UBWC«£jftL5C<LUT.'IC«H^T.E.P=2^UB=1饮]LMIE.=3SEjQV1DIMX=・USE.-03aww=・3aaE.#as酗選=.24SE；ttaaAN212011.3QlE.ttaS .35^E.«0!5 .€7'5E.tt0!5 .SS7E.*aS .132E.*10.151E.*DS .'514E.ttaS .83^E.*QS5000ra.d/s_yiKiglifemxi_^i3o'Ti)a-gjW]n_0S01^.1827分析：图分析：图1转速为5000rad/s时叶片等效应力分布图理论上叶片自上到下应力应该逐渐增大，最小应力MN发生在叶尖部，最大应力MX发生在叶根部。因为在这里叶片可以简化的看成根部完全固装的等截面悬臂杆。把叶片网格划分成有限个微元单元体后，在5000rad/s离心力的作用下，靠近外层的微元单元体所受到的外侧材料的总的离心应力较小，越靠近根部时，截面外侧所有材料的离心力都将加载到该截面上，所以越靠近根部，截面所受到的总的离心应力就越大。用ANSYS软件建模求解后，所得到的叶片应力分布图大致符合理论分析。上面所显示的叶片最底面并不是红色的最大应力区域，这是因为我们在划分网格时，ANSYS软件模拟出的地面固定边界载荷是作用在底面微元单元体上的节点（nodals）上的，而并不是均布在整个底面平面上，所以输出结果不会出现底面是均布的最大载荷的区域。但是从叶片应力分布图来看，总的应力趋势还是从叶片自上到下，应力逐渐增大的，和实际相符。

3.叶片振动的有限元分析叶片静频计算与分析基于ANSYS9.0软件的计算过程：先是建立叶片模型(六面体)，建模过程中要通过关键点确定它在空间中的精确的相对位置；然后是选择模态分析(modal),接着设定10阶模态分析，固定底面边界，然后是求解(solve)；最后通过ReadResults和PlotResults查看1-10阶各阶振动位移图(DOFsolutions),结果如下：(j) 10阶振型图，属于弯曲扭转复合振动(j) 10阶振型图，属于弯曲扭转复合振动图2前1-10阶振型图转(rad/s)阶数1000200040008000100001050.51137.213642026.835934405.74077.24100.24168.54430.65348.55941.26560.76639.36869.67719.28400105699594.49604.89635.89758.51200012023183611844218682196102288724362222152223722305225742361425017292762929829364296273064331371324003240332412324503260032712360093609536350373524106343595503955042950533509435253453678转(rad/s)阶数1000200040008000100001050.51137.213642026.835934405.74077.24100.24168.54430.65348.55941.26560.76639.36869.67719.28400105699594.49604.89635.89758.51200012023183611844218682196102288724362222152223722305225742361425017292762929829364296273064331371324003240332412324503260032712360093609536350373524106343595503955042950533509435253453678总结：除弯曲和扭转振动外，在叶片上还会出现许多其他振型。其中有的是弯曲和扭转的复合振型，有些振型还难以给以名称。其中一阶弯曲振动、二阶弯曲振动、一阶扭转振动„„较为常见，危险性也最大。对于压气机叶片而言，最重要的是一、二弯和一扭振型；对于涡轮叶片，大多是一弯和一扭振型。（2）叶片动频计算与分析基于ANSYS9.0软件的计算过程：首先也是建立叶片模型（六面体）建模过程中也要通过关键点确定它在空间中的精确的相对位置（可以精确地绕X轴旋转）；然后是选择静态分析（tatic），计算前要勾选考虑预应力的影响，把第一次solve求解出的对应转速的离心应力关联到下一步的模态分析中;再然后是选择模态分析（modal）,勾选考虑预应力的影响，设定10阶模态分析，固定底面边界，然后是第二次求解（solve）；最后通过ResultsSummary查看1-10阶的对应转速下的动频值，将结果列表如下：表一各转速下各阶振型的振动频率值10阶

（3）共振分析根据倍频力与转速的关系，将数据列表如下：表二倍频力的频率值与转速的关系6000050000400003000020000100000 1000 2000 3000 4000 5000 6000 7000 8000 90001000011000转速n(rad/s)亠6000050000400003000020000100000 1000 2000 3000 4000 5000 6000 7000 8000 90001000011000转速n(rad/s)亠1阶亠2阶3阶—4阶—5阶亠6阶—7阶一8阶—9阶10阶—K=1__K=2K=3—K=4—K=5—K=6—K=7——K=8-K=16K=30^速^K值(rad/s)\^K=1K=2K=3K=4K=5K=6K=7K=8K=16K=301000159.1318.3477.4636.6795.7954.91114.01273.22546.44774.54000636.61273190925463183381944565092.810185.6190981000015913183477463667957954911140127322546447745图3叶片共振图观察上面的叶片共振图，现对共振现象进行分析：作用在叶片上的一个局部的冲击力可以看作是许多谐力之和，谐力的频率为转速的1、2、3„„倍。或称这些力为1、2、3„„倍频力。倍频力也就是所谓的激振力，它可以分为两类：机械激振力和气动激振力。前者是由于轮盘有振动，因而摇动叶片根部，使叶片发生振动。通常称为“位移激振”或“位移激扰”；后者是由于气流对叶片表面的压强做周期性的变化，激起叶片振动。局部障碍将引起各种频率的