撰文/河北工程大学 水电学院 范江峰 孙中伟 李丙尧
SolidWorks Simulation在水轮机转轮模态分析中的应用
撰文/河北工程大学 水电学院 范江峰 孙中伟 李丙尧
SolidWorks Simulation作为一款强大的有限元分析软件,其便捷性与易用性已得到公认。本文利用SolidWorks Simulation对某模型水轮机的转轮分别进行了零件状态下的干模态分析,并试算转轮与水体装配体状态下的湿模态分析。分析证明SolidWorks Simulation可以方便、快速的进行零件相关的设置和分析,计算结果可以为转轮的优化设计提供理论依据,但在湿模态分析中还存在不足。
在水轮发电机运行时,由水力因素,机械因素和电磁因素引起的各种振动中,对机组运行影响最大的是振动频率与水轮机固有频率接近的振动。在该振动频率下,往往引起机组的共振,从而导致机组部件金属和焊缝疲劳破坏区的形成和扩大,从而产生裂纹甚至断裂而导致部件的报废。因此,确定水轮机的固有频率,并使机组运行避开这个振动区域对机组的安全运行有重大的意义。
SolidWorks Simulation是由COSMOSWorks软件发展而来的基于有限元分析技术的设计分析模块。COSMOSWorks软件是由SRAC公司与SolidWorks公司合作开发的插入式分析软件,可以与SolidWorks的CAD系统直接对接。软件提供了丰富的机械有限元分析功能,包括静应力分析、热力、频率、屈服和跌落测试等。软件界面十分友好,使分析人员可以方便快捷的对各种机械进行有限元分析。
本文通对利用SolidWorks Simulation对某模型水轮机转轮进行模态分析,得出使用SolidWorks Simulation进行转轮干模态分析的方法以及在进行湿模态分析时上辈子在的的不足。
在SolidWorks中建立转轮模型(图1),转轮直径6.9m,叶片数17片,由于在实验中存在频率分析后出现叶片、转轮上冠以及下环分离的情况,在分析前将这三部分进行组合运算,使其成为整体。设置转轮材料为0Cr13Ni5Mo,材料基本属性为:弹性模量2.1×1011N/m2,中泊松比0.3,中抗剪模量3.189×108N/m2,质量密度7850kg/m3。在SolidWorks Simulation中新建算例,选择新算例,选择频率。
由于水轮机转轮上冠与主轴是通过法兰相连接,所以将此处设置为夹具,选择固定几何体。SolidWorks Simulation中提供了方便的四面体网格划分方式。四面体单元由于纵横比小,相比于六面体单元,往往需要划分更多的网格才能描述同一个模型。但由于SolidWorks Simulation的内置算法为FFE(Fast finite Element)快速有限元算法,因此在网格数比较多的情况下,SolidWorks Simulation依然保持了较快的运算速度。本例中利用网格工具自动生成网格,并对叶片、上冠、下环部分应用网格控制实现局部网格加密(图2)。在网格划分过程中偶尔会出现网格划分失败的情况,通过调整网格尺寸,修改网格控制选项以及雅可比点数等,一般可以解决划分网格失败的问题。
图1 在SolidWorks中建立的转轮模型
图2 加密后生成的网格
应用网格细节可以清楚的看到生成网格的详细信息,共生成118734个节点,70379个网格。在有限元分析中网格的质量与分析结果的精度有着直接的联系,因此划分网格时一定要确定网格控制的策略,保证生成网格的质量。
干模态分析是指在空气环境下,对分析对象进行模态分析。运行次算例后,得出转轮在空气中振动的前五阶模态。
由转轮在空气中的前五阶模态图可以看出,位移最大的部分总出现在转轮下环处,这里在运行状态下还要承受水对转轮冲击所带来的载荷,分析结果与实际一致,频率值在可接受的范围内。
图3 第一阶模态(干)
图4 第二阶模态(干)
图5 第三阶模态(干)
图6 第四阶模态(干)
图7 第五阶模态(干)
图8 列举频率(干)
由于水轮机转轮是运行在水体中,周围水体的存在必然影响转轮的振动频率,因此进行水环境下的湿模态分析是十分必要的。SolidWorks Simulation中并未提供单独的湿模态分析选项。鉴于湿模态分析是水体与转轮之间的振动频率分析,本文利用装配体频率分析进行了如下试算。
在SolidWorks中建立一个圆柱型水体,插入转轮零件进行组合删减运算,得到包含转轮型空腔的水体零件,在SolidWorks环境下可以方便的更改圆柱水体的半径及高度参数。在装配体模式下将水体与转轮进行装配,得到被水体包围的转轮。模型如图9所示,设置水体材料为水。由于SolidWorks Simulation下进行频率分析要求材料弹性模量不为0,参考鲁丽等发表在《西南交通大学学报》中的研究结果《反应堆吊篮在空气和静水中的振动特性分析》,取水体的材料弹性模量为转轮材料的1%,即2.1×109N/m2,泊松比取值相对较小,取为0.3×10-4,同样设置转轮上冠与转轴法兰相连处为固定约束,分别对转轮与水体应用网格控制进行网格划分,转轮叶片处网格应适当加密。
图9 水体与转轮的装配体模型
表 不同尺寸水体装配体对应频率表
设置不同的水体尺寸,得到不同水体尺寸下装配体的固有频率如表所示。结果显示这种方法下进行的湿模态分析频率,由于水体尺寸的不同,装配体的固有频率也不同,参考梁权伟等发表在《水力发电学报》中的研究结果《考虑流固耦合的混流式水轮机转轮模态分析》给出的数据,可知干湿模态下的频率影响系数分辨值介于0.671~0.904间相差较大。原因是所得的固有频率不是在水体阻碍作用下转轮的固有频率,而是水体与转轮所组成装配体的固有频率。且由结果看出不同水体形状会产生不同的固有频率。
(1)通过在SolidWorks Simulation中进行水轮机转轮的干模态分析,可以发现Simulation插件可以快捷方便的进行不考虑外部环境的频率分析,在进行干模态分析时利用Simulation可以提高工作效率。
(2)由于Simulation中没有湿模态分析的相关单元,其中的外部载荷选项中提供的流动效应是针对流水载荷,并不能进行静水下的频率分析。本文提出了一种通过构建水体装配体模型的方法,结果证明同样不能有效模拟静水环境下转轮的振动情况。
(3)Simulation作为一款方便的有限元分析工具,其简便性深受好评,但由于软件功能相对简单,在今后的发展中依然要不断更新升级自身的功能,才能不断在工业制造中发展壮大