汽轮机的结构静力学分析.docxVIP

汽轮机的结构静力学分析 主要汽油公司为双吸式离心。该泵是汽轮机油系统的关键设备之一,主油泵运行正常与否将直接影响机组的安全可靠性。叶轮是主油泵内的旋转部件,其工作特点是转速高,达到3 000 r/min以上,工作压力等级高,高达1.8 MPa,因此该叶轮是主油泵结构分析中的重点。本文采用有限元软件ANSYS对主油泵叶轮进行结构静力学分析,通过这一分析可以掌握叶轮部件的应力、变形等力学特性,并为以后的模态分析和瞬态分析打下良好的基础。 1 叶片静力学分析及计算过程 1.1 叶轮网格划分 由于离心泵叶轮具有复杂的空间几何形状,其叶面是由多个复杂的雕塑曲面组成,因此在三维CAD软件UG中建立起了它的实体模型。然后将叶轮实体模型存储为Parasolid格式的文件,ANSYS有Parasolid格式文件的接口,导入该文件并显示叶轮的实体模型(如图1所示)。 在ANSYS中检查叶轮实体模型并确认无误,然后设置叶轮的材料属性(如表1所示)。接下来选择单元类型,由于ANSYS是通用有限元软件,可以进行多种类型的分析,针对不同的分析类型它提供了二百多种单元类型。本文是对叶轮进行结构静力学分析,需要选择与结构分析相关的单元类型。首先叶轮体上承受惯性载荷,这样就需要应用结构分析中的体单元,同时由于叶轮体形状复杂,难以用结构网格对其进行离散,因此我们选择与非结构网格相对应的solid187号单元。该单元为10节点四面体单元,具有2阶精度,适用于空间结构的弹性、塑性、超弹性、蠕变等结构分析。然后叶轮表面还要承载油压力,以及键槽表面要承受方键的压力,这些载荷都可以选用surf154号单元来加载。该单元为空间结构表面效应单元,适用于三维表面变载荷或表面效应的加载分析,并且该单元可以与solid187号单元相适应产生三角形的表面单元。 在选择好单元后,就要用相应的单元对叶轮划分网格,ANSYS提供了自由网格与映射网格的划分方式。由于叶轮形状较复杂,难以使用映射网格划分该几何体,所以我们选择使用自由网格划分方式。另外,ANSYS还提供了智能网格划分功能,智能网格划分功能是一种比较高效的自由网格划分方法,它考虑几何图形的曲率以及线与线的接近程度自动进行网格划分,共分为10个级别,级别越高划分的越粗,反之越细。经过多次划分,在使用6以下级别划分叶轮网格时失败,所以使用级别为7智能网格划分功能对叶轮进行网格划分。划分好网格后的叶轮的单元个数为226 553个,节点个数为305 704个,叶轮网格如图2所示。至此完成了叶轮有限元模型的建立。 1.2 叶轮的压力状态 计算工况。 叶轮的设计工况为3 000 r/min,设计的最大转速工况为3 600 r/min。针对设计要求,我们将计算工况也设置为3 000 r/min和3 600 r/min两种工况。 受力分析。叶轮工作时是由主轴驱动做高速旋转,主轴通过方键将扭矩传递给叶轮,因此叶轮的键槽表面要承受的由主轴传来的扭矩。同时在泵腔内充满了透平油,由叶轮旋转将油泵到油道中,因此叶轮的叶面上还要承受很大的油的压力。此外由于叶轮的转速极高,因此在叶轮旋转时还要产生极大的惯性力。还有叶轮自身的重力,但由于叶轮的转速高,惯性大,而且叶面承受的油压也较大,所以与上述两种力相比叶轮自身的重力可以被忽略不计。 施加载荷。首先我们分析叶轮在空气中旋转时的情况,这时只在叶轮上加载由其自身旋转所产生的惯性力,此时只用187号单元即可。由于叶轮并非完全对称,因此需要在轮毂内的四个键槽底面上均匀选取8个节点,在这些节点上施加周向和轴向位移约束。这样才能平衡由于叶轮的不对称性所产生的轴向位移,并得到完全的约束方程使得计算收敛,同时还不会影响计算结果的真实性。 其次我们分析叶轮在正常工作时的情况,这时除了加载上面所有的载荷和约束外,还要在叶片表面加载油压力,压力数据由FLUENT中提取并读入ANSYS中,加载压力载荷后的叶片表面单元如图3和图4所示。另外在键槽向着旋转方向一侧的表面上要加载由驱动轴传来的压力,这时除了用187号单元外还需要使用154号单元来加载上面的两种压力载荷。通过有无油压两种载荷状态的计算,可以分析比较得出油压力和键槽受力对叶轮的应力及变形的影响。 求解计算。进入ANSYS求解器,设置分析类型为结构分析,求解器类型为PCG,其他均取默认值,分别对4种载荷情况进行求解计算。 1.3 结构分析方法 计算完成后,进入ANSYS的后处理器post1,在这里ANSYS提供了关于结构分析的多种计算结果,而且既能用云图表示也可进行数据列表,读取叶轮的应力及变形数据结果,提取各工况下的应力和变形图(如图5-图12所示)。 1.4 叶轮应力及变形的特性 首先对比两种转速工况下的应力及变形。由图5-12看出叶轮在3 600 r/min时的应力及变形