水轮机叶片

一、工程背景及水轮机叶片简介

图 1、为某型水轮机叶片的CAD模型。在发电工作工程中水流由进水口流向出水口，叶片承受水流的冲刷从而开始运动，这种运动通过传动轴传递到发电机，从而带动发电机工作发电。但是水轮机在工作仅仅一年多时间以后，就有数片叶片发生了疲劳断裂事故，使得水轮机不能正常工作发电，造成了一定的经济损失，同时也说明水轮机叶片在结构的设计方面确实存在不完善之处。然而，由于水轮机在水下进行工作，很难通过测量得方法获得叶片上应力和位移的分布情况，也就无法知道叶片为何会断裂，无法有效的改善叶片的几何结构。在这种情况下，长江水利委员会陆水枢纽局的委托我们对LS591水轮机叶片的进行Ansys有限元模拟计算，获得叶片的应力场和位移场的分布，从而为叶片断裂事故分析提供技术支持，并对叶片结构的改进提供具体方案。

传动轴

进水口出水口

图1、CAD模型

二、ANSYS简介及解题步骤

1、ANSYS简介

对于大多数工程技术问题，由于物体的几何结构比较复杂或则问题的某些特征是非线性的，我们很难求得其解析解。这类问题的解决通常具有两种途径：一是引入简化假设，但这种方法只是在有限的情况下是可行的。也正是因为这样，有限元数值模拟的技术产生了。有限元方法通过计算机程序在工程中得到了广泛的应用。到80年代初期，国际上较大型的面向工程的有限元通用软件达到了几百种，其中著名的有：ANSYS,NASTRAN,ASKA, ADINA,SAP等。其中，以ANSYS为代表的工程数值模拟软件，即有限元分析软件，不断的吸取计算方法和计算机技术的最新进展，将有限元分析、计算机图形学和优化技术相结合，已成为解决现代工程问题必不可少的有力工具。尤其是在某些环境中，样机试验是不方便的或者不可能的，而利用ANSYS软件，对这个问题有了很好的解决。本文中水轮机叶片是在水下的环境进行工作，测量很难进行，利用有限元软件ANSYS这个问题得到了很好的解决。

2、ANSYS分析步骤

ANSYS分析可以分为三个步骤：

a、创建有限元模型

（1）创建或读入几何模型

根据实体模型按照给定的尺寸建立模型或者直接导

入已经生成的几何模型，并对其进行一定程度的修

复、简化等。

（2）定义单元类型，设定实常数、定义材料的属性定义单元类型：对于任何分析，必须在单元类型

库中选择一个或几个合适的单元类型，单元的类型决定了附加的自由度（位移、转角、温度等）。许多单

元还要设置一些单元的选项，诸如单元特性和假设

等。

设定实常数：有些单元的几何特性，不能仅用

其节点的位置充分表示出来时，就需要提供一些实

常数来补充几何信息。

定义材料属性：材料属性是与几何模型无关的本

构属性，例如杨氏模量、密度等。虽然材料属性并

不与单元类型联系在一起，但由于计算单元刚度矩

阵时需要材料属性，所以在此我们要对材料的属性

进行相关的定义。

（3）划分网格（节点及单元）

在做好上述的所有工作后，接下来就是对实体模型进

行网格划分，此步尤为关键，因为网格划分的好坏将

直接影响到计算结果的精确度与收敛性。根据模型的

拓扑结构决定采用映射网格还是自由网格，之后对网

格的尺寸进行设定，对关心的部位或者危险部位进行

必要的网格细化。总之,网格要足够细，才能保证结

果的精确性。

b、施加载荷并求解

（1）施加载荷及载荷选项、设定约束条件

施加约束根据具体情况对有限元模型进行约束设定

施加载荷包括集中载荷、面载荷、体载荷、惯性载荷

等。

（2）求解

选择求解器 ANSYS提供了两个直接求解器：波前求解器、稀疏矩阵求解器，同时还提供了三个迭代求解器：PCG、JCG、ICCG。因此，在前根据具体情况选择合适的求解器，这样直接影响求解的速度和结果的精确度。

进行求解

c、后处理

（1）查看结果

静态分析的结果写入结果文件，结果由以下数据构成：

基本数据——节点位移（UX,UY,UZ,ROTX,ROTY,ROTZ）

导出数据——节点单元应力、单元应变、单元集中力、节点反力等。

（2）分析结果

可直接LIST结果数据，也可通过等值线、矢量图等形式

对结果进行观察分析。

三、叶片几何模型的建立

叶片的几何模型根据相关的图纸建立，CAD模型如图1、图2所示。

图2、几何模型

我们把整个叶片离散成474个小六面体分别建立，其目的有三：

●建立模型的需要

从图2、知道叶片的几何结构不是很规则，其上下两个表面都是形式比较复杂的超曲面。基于这种情况，很难建立一个和实际叶片一模一样、丝毫不差的模型，只要把叶片离散成474个足够小的小六面体逼近实际模型，这个问题就得以解决。为了保证叶

片几何模型上下两个表面光滑，我们在建模时采用了以下的方案：

A

B

C

D

E

F

G

H

图3、建模方案

其中A,B,C,D,E,F,G,H分别为相邻小六面体的顶点，首先通过ABCD,EFGH建立两条三次样条插值曲线，然后建立曲线ED和曲线ABCD,EFGH相切，这样就使得ABCDEFGH成为一条光滑的曲线。使用这种方法，也就保证了叶片的上下两个表面都是光滑的曲面。

由于叶片的上下两个面是通过叶片上的一些离散点三次样条插值得到，所以几何模型和实际模型的逼近效果较好。

●划分高精度六面体映射网格的需要

映射网格比自由网格具有更高的计算精度。然而划分六面体映射网格对模型的拓扑结构有严格的限制：只有形状较规则的六面体和三棱柱才能划分映射网格。我们把叶片分成474个小六面体，这些小体形状不很奇异、比较规则，正好可以满足划分映射网格的要求。

●在叶片上加载的需要