轴流泵叶片有限元剖析.docVIP

轴流泵叶片有限元剖析 　　摘 要：利用UG软件建立轴流泵叶片的三维模型，采用ANSYS CFX软件计算得到叶片上下表面的水压力分布，通过ANSYS Workbench软件强大的耦合分析功能将水压力加载到叶片有限元模型上，最后利用ANSYS软件进行求解计算，得到应力分布结果，根据第三强度理论做应力评定。 　　关键词：轴流泵 UG 流固耦合 应力分析 　　中图分类号：TH312 文献标识码：A 文章编号：1007-3973（2013）003-021-03 　　1 引言 　　轴流泵属于低扬程、大流量型泵，结构简单，使用方便，主要用于农田排灌， 也可用于城市给水排水、热电站中输送循环水、船坞升降水位等。 　　UG的一个最大特点就是混合建模，在一个模型中允许存在无相关性的特征。如在建模过程中，可以通过移动、旋转坐标系创建特征构造的基点，在曲面造型中将会充分体现出来。 　　ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，是现代产品设计中的高级CAE工具之一。ANSYS很早便开始进行流固耦合的研究和应用，目前ANSYS中的流固耦合分析算法和功能已相当成熟，可以通过或者不通过第三方软件（如MPCCI）实现ANSYS Mechanical APDL+CFX、ANSYS Mechanical APDL+FLUENT、ANSYS Mechanical+CFX的流固耦合分析。 　　2 叶轮模型的建立 　　液体在轴流泵叶片内的运动，是一种复杂的空间运动。任何一种空间运动都可以认为是三个互相垂直的运动合成。研究水流在轴流式叶片中的运动时，为了方便起见，我们采用柱坐标系。在分析和设计轴流泵叶片时，提出了圆柱无关性假设。液体质点在以泵轴线为中心的圆柱面上流动，且相邻圆柱面上的液体质点的运动互不相关。即在叶轮的流域中，不存在径向分速度，圆柱面即流面。 　　根据圆柱层无关性假设，可以把叶轮内复杂的运动简化为圆柱流面上的流动。在叶轮内可以作出很多这种圆柱面，圆柱面沿母线割开后，可以展开在平面上。圆柱面和各叶片相交，其截面（翼型剖面或翼型）在平面上构成一组叶栅。 　　由上述分析可知，轴流泵叶轮叶片主要由圆柱面上的型线组成，在水力图上，每条型线都有若干各点的圆柱坐标组成。利用这一点我们可以将叶片各断面木模图导入UG中，然后通过一系列的平移，旋转，投影命令得到轴流泵叶轮圆柱面上的型线，将各型线用样条曲线连接起来，生成沿各型线曲面，并将所有曲面缝合，最后实体化模型得到轴流泵叶片三维模型。图1为叶片翼型空间型线，图2为叶轮三维实体模型。 　　轴流泵具体几何参数如下： 　　轮毂直径d=100mm，叶轮外径D=300mm，叶片数Z=3，转数n=1450r/min，弹性模量E=2.06??2GPa，泊松比 =0.3，密度 =7850Kg/m3，叶片材料为ZG0Cr13NiMo，屈服极限 s=550MPa。 　　图1 叶片翼型空间型线 　　图2 叶轮三维实体模型 　　3 叶轮有限元模型的建立 　　叶片与轮毂连为一体，泵运行过程中叶轮与叶片无相对运动，在约束施加的过程中将轮毂完全固定即可。轴流泵在正常运行时，叶片在泵壳内高速旋转，泵内水流场极其不均匀，叶片表面所承受的水压力在空间上呈无规则变化状态，并且水压力是在ANSYS CFX中计算得到的，可以通过ANSYS Workbench强大的耦合功能加载到叶片有限元模型上。 　　在计算复杂空间扭曲结构的三维模型时，实体单元SOLID92具有良好的适应性，故叶轮单元采用SOLID92四面体单元进行网格划分，在ANSYS Workbench的网格处理器中采用半自动方法生成叶片网格，对叶片根部靠近应力释放槽的过渡区域的网格由手动进行更细网格密度的划分，生成的有限元网格如图3所示。本文共划分四面体单元79685，节点总数127083。采用国际标准单位制。 　　图3 叶轮有限元网格 　　图4 叶轮有限元网格过渡区加密 　　图5 叶片上表面水压力加载结果 　　图6 叶片下表面水压力加载结果 　　在算法上，ANSYS主要采用分离式解法也就是载荷传递法求解流固耦合问题。从数据传递角度出发，流固耦合还可以分为两种：单向流固耦合分析和双向流固耦合分析。轴流泵的结构强度分析属于单向耦合分析，耦合交接面处的数据传递是单向的，CFD分析计算的结果传递给固体结构分析，固体结构分析的结果不用传递给流体分析。也就是说，流体分析对结构分析的影响很大，而结构分析的结果对流体分析的影响很小，甚至可以忽略不计。 　　不同工况点的载荷是由流体和机械载荷综合而成的，对于叶片水压力施加采用的是顺序耦合法。在ANSYS Workbench中，叶片的压强可以通过插值法快速准确的加载到叶片单元表面的节点上。