汽轮机叶片静动特性的随机有限元方法研究

摘要

在叶片的设计、加工、安装、运行各个环节，都存在大量的不可测或不可控的因素，从而导致叶片的静动特性出现随机性。目前叶片静动特性分析主要采用忽略参数随机性的确定性物理模型，确定性模型很难从理论上解释叶片静动特性的分散性，更难解释按确定性方法设计合格的叶片在使用中却发生损坏的现象。本文提出采用随机有限元模型描述叶片静动特性，以反映叶片参数的随机性，并得到叶片静动特性的概率特征，这对于揭示汽轮机叶片静动特性的随机性具有重要的理论意义，对于汽轮机叶片的可靠性设计、制造、运行具有重要的应用价值。 本文分析了影响叶片静动特性的各参数的特点，提出采用随机变量描述转子转速的随机性，采用随机场描述叶片弹性模量、密度、结构尺寸、蒸汽作用力、阻尼、叶根连接条件的随机性，并用收敛性好和精度较高的随机场局部平均法离散随机场；基于摄动技术推导出叶片动力分析的随机Hamilton变分原理，采用20节点等参元建立了叶片动力学随机有限元方程，给出了单元线弹性刚度矩阵、离心刚度矩阵、质量矩阵、离心力向量、蒸汽作用力向量的均值和一阶变异量的表达式；基于摄动技术推导出叶片静态问题随机最小势能变分原理，采用20节点等参元建立了静态特性随机有限元方程；推导了用于叶片动频率计算的初应力刚度矩阵的均值矩阵和一阶变异矩阵的表达式；提出了用随机弹性约束单元来描述叶片根部约束的随机性，并推导了随机弹性约束单元刚度矩阵的均值矩阵和一阶变异矩阵的表达式；给出了叶片静态响应、固有频率和振型、稳态动力响应的均值、灵敏度、方差、变异系数的求解方法。 在理论研究的基础上，采用面向对象随机有限元程序设计思想，编写了叶片随机有限元分析程序。对随机参数等直叶片的静态响应、固有振动特性、稳态动力响应进行了随机有限元分析，通过与解析结果及MonteCarlo模拟结果的对比，验证了本文方法的精度。 本文对随机参数432叶片的静态响应、固有振动特性、稳态动力响应进行了随机有限元分析，详细研究了叶片静动特性对随机参数的灵敏度，给出了参数随机时，叶片静动特性的均值和变异系数，通过分析发现叶片参数的变异将导致叶片静动特性产生显著的变异。利用本文方法和程序可以在叶片设计阶段预测其静动特性的均值和变异，同时还可以获得叶片静动特性对参数的灵敏度，这对于叶片可靠性设计、控制频率分散度和现场调频具有重要意义。

目录

文摘

英文文摘

论文说明：主要符号表

声明及关于学位论文使用授权的说明

第一章前言

1.1课题背景及意义

1.2汽轮机叶片静动特性研究现状

1.2.1确定性力学模型

1.2.2考虑参数不确定性的力学模型

1.2.3叶片振动特性敏感性

1.3随机有限元方法的研究现状

1.3.1随机变分原理

1.3.2随机场描述方法

1.3.3静动力响应随机有限元方程求解方法

1.3.4随机特征值问题求解方法

1.4本文研究内容

第二章叶片参数不确定性的描述方法

2.1随机场局部平均离散方法

2.1.1基本理论

2.1.2规则六面体单元局部平均随机场二阶特性

2.1.3任意六面体单元局部平均随机场二阶特性

2.1.4空间曲面上的局部平均理论

2.1.5算例分析

2.2多随机场离散的局部平均法

2.3 Monte Carlo模拟技术

2.3.1均匀随机数的产生

2.3.2非均匀随机数的产生

2.3.3概率密度函数的级数拟和

2.3.4算例

2.4小结

第三章叶片静动特性三维随机有限元方程

3.1随机Hamilton变分原理

3.2叶片动力学随机有限元方程

3.3叶片动频率计算的初应力刚度矩阵

3.4随机最小势能原理

3.5叶片静态特性随机有限元方程

3.6小结

第四章叶片静态响应的随机有限元分析

4.1摄动随机有限元法

4.2 Neumann展开随机有限元法

4.3叶片稳态汽流作用力

4.4叶根约束随机性处理方法

4.5汽轮机转速的随机性

4.6静态分析摄动随机有限元程序流程

4.7算例

4.7.1随机变量有限元

4.7.2随机场有限元

4.8小结

第五章叶片固有振动特性的随机有限元分析

5.1叶片静动频率的均值及其灵敏度

5.2叶片静动频率的统计特性

5.3叶片静动频率随机有限元分析的程序流程

5.4算例

5.4.1随机变量有限元

5.4.2随机场随机有限元

5.5小结

第六章叶片稳态动力响应的随机有限元分析

6.1叶片动力响应统计特性的求解方法

6.2叶片激振力的处理

6.3叶片阻尼的处理

6.4叶片稳态动力响应随机有限元程序流程

6.5算例

6.5.1随机变量有限元

6.5.2随机场随机有限元

6.6小结

第七章随机参数汽轮机叶片静动特性分析实例

7.1面向对象的随机有限元编程思想

7.2基于Ansys二次开发工具APDL语言的叶片前后处理方法

7.3 432叶片静态特性的随机有限元分析

7.4 432叶片固有振动特性的随机有限元分析

7.5 432叶片动力响应的随机有限元分析

7.6小结

第八章结论

参考文献

致谢

个人简历、在学期间发表的学术论文及研究成果