离心压缩机三元叶轮形状优化设计.doc

离心压缩机三元叶轮形状优化设计 刘万青 王跃方/大连理工大学工业装备结构分析国家重点实验室 张 勇 杨树华/沈阳鼓风机集团有限公司 摘要：在三维参数化建模的基础上，对某离心三元叶轮进口叶片的前缘进行了形状优化设计，在叶轮质量和变形满足要求的条件下，得到了叶轮进口叶片前缘的最优结构尺寸参数，降低了叶轮最大应力并改善了应力分布状态，为离心压缩机三元叶轮设计提供了重要的依据。关键词：离心式压缩机；三元叶轮；优化设计 中图分类号：TH452 文献标识码：B文章编号：1006-8155（2008）04-0031-04Optimization Design for the Shape of 3-D Impeller of Centrifugal CompressorsAbstract：The shape optimization design for the front edge of the inlet blade of some centrifugal 3-D impeller is presented based on the three-dimensional parametric modeling. The optimum structural parameters are obtained for the front edge blade of inlet impeller under the requirements of impeller quality and distortion，the maximum stress of the impeller is brought down and the distribution of stress is improved. The shape optimization design provides important foundation for designing centrifugal compressor 3-D impeller.Key words ：centrifugal compressor；3-D impeller；optimization design0 引言离心压缩机叶轮叶片沿长度和高度方向是变化的，模型较为复杂，在转动中承受着多种载荷作用，故在许多计算中都采用简化为等厚度叶片，采用轴对称有限元模型1进行计算，此时叶片被简化成等效的盘形夹层，蔡显新等2利用叶片轮盘藕合分析方法进行求解 。这些方法都不能真实地反映叶片的强度和刚度。而在工程应用中叶片正是薄弱环节，许多叶轮的破坏形式正是叶片撕裂。只有对叶轮进行真实的三维有限元模拟才能得到叶轮的真实应力分布3-4。随着技术的进步，压缩机的转速、流量朝着更高的方向发展，叶轮的直径和出口宽度也变大，人们对叶轮强度方面的要求也越来越高。因此，对离心叶轮进行有限元分析、强度校核和优化设计越来越受到人们的重视。本文研究的高效大三元叶轮，其出口宽度和叶片的空间扭曲较大，强度问题非常突出。在叶轮的原始设计中，经过有限元结构分析发现叶片进口前缘的应力值很大，属于危险状态，在复杂工况及交变载荷的作用下很容易引起叶片的断裂破坏，在实际工程应用中也发生过类似的破坏情况。通过力学分析和优化设计，降低了其最大应力，提高了叶轮的刚度。1 叶轮形状优化设计以有限元分析软件ANSYS为工作平台，使用其二次开发语言APDL，采用参数化设计的思想，使有限元建模与分析实现参数化。创建了参数化的分析流程后利用ANSYS的优化模块对叶轮的设计参数进行优化设计，从而达到最优化设计目标。收稿日期：2008-05-20 大连市 116024叶轮运转时所受载荷包括: 离心力、轴向力、气动力及温度载荷,并以离心力为主，笔者考虑的载荷为离心力，在叶轮的不同半径和轴向位置会产生不同的切向应力和径向应力，叶片的最大应力由于叶型和厚度分布不同, 其发生部位也不同。以往对叶轮进行形状优化主要针对叶轮轮盘和盖盘的外边界，同时，对于离心叶轮的优化设计，均以减轻叶轮质量或者降低叶轮应力作为优化目标2,4-5。然而在三元叶轮的原始结构设计中，经过有限元结构分析发现叶片进口前缘叶片与盖盘连接处及叶片与轴盘连接处应力值很大，属于危险状态，在复杂工况及交变载荷的作用下很容易引起叶片的断裂破坏。因此，针对叶轮进口叶片前缘形状作了优化设计，选取最大应力作为目标函数。图1 为高效大三元叶轮的剖面图,选取叶片前缘边界B,C,D,E,F,G 6点控制的样条曲线及叶片厚度为设计变量。由于叶轮属于循环对称结构，受离心力作用时，圆周方向上每个叶片相同的位置的应力基本相等，因此设计变量选取一个叶片上的坐标即可，而目标函数选取为所有节点最大应力，所作的优化为降低最大应力。图1 离心叶轮截面图边界约束要考虑到零件之间的装配关系、加工等因素并结合经验给出，与性态约束与目标函数的选取有关。优化问题统一表示为subjected to （1） 其中表示叶轮某关键部位的径向位移和轴向位移；表示容许位移；(X,Y,Z)，t为设计变量，分别为B,C,D,E,F,G点的坐标值及叶片厚度；为目标函数，为B,C,D,E,F,G点的坐标值的函数，选取为,表示叶轮最大von Mises等效应力。2 参数化建模以有限元分析软件ANSYS为工作平台，使用其二次开发语言APDL6-7语言实现叶轮的三维参数建模，具体建立模型方法如下。（1）将轴盘及盖盘截面关键点坐标数据整理为ANSYS规定三格式文件，利用*VREAD命令读取数据文本文件，将轴盘数据存储为数组或者表参数，利用APDL宏文件循环结构，生成轴盘及盖盘截面的关键点，由这些点连成线，再生成轮盘及盖盘截面，旋转成轮盘和盖盘。（2）将叶片的坐标数据整理为ANSYS规定格式文件，利用*VREAD命令读取数据文本文件，将叶片数据存为数组或者表参数，以便随时调用。（3）以变量形式定义B,C,D,E,F,G坐标，利用APDL的循环结构，生成一系列叶片关键点，生成参数化叶片前缘曲线,由叶片数据生成叶片中面，再利用叶片厚度变量生成厚度参数化的叶片。 （4）将叶片旋转复制生成全部叶片，并通过布尔运算得到叶轮实体。图2为某离心叶轮实体模型图。参数化模型建立后，即可施加约束以及离心力载荷，对结构进行有限元分析。图2 离心叶轮实体3 计算及结果分析三元叶轮由轴盘、盖盘和叶片构成，轮径为1100mm，轴径为337mm，轴孔长度145mm。叶轮共有叶片19个。叶轮材料为某型号钢，弹性模量为2.1E5MPa，泊松比为0.3，密度为7860kg/m3，材料许用应力为1100MPa。单元选择10节点四面体三维实体单元,单元类型为SOLID187，应用软件自动划分单元功能，选取网格尺寸为20mm，进行网格划分。对轮盘的轴孔进行周向以及轴向约束。作用载荷离心力载荷，转速为n=491.9r/s。图3是叶轮初始设计的von Mises等效应力云图，图4和图5分别为叶轮径向和轴向位移云图。从图中看出，最大应力出现在叶片进口前缘，最大应力值为1272MPa，大于材料容许应力1100MPa，处于危险状态，极易破坏，在工程应用中也发生过类似的破坏事故。故有必要对该叶轮进行优化设计。图3 初始设计叶轮von Mises应力云图 图4 初始设计叶轮径向位移云图图5 初始设计叶轮轴向位移云图经过迭代计算，得到优化结果，优化后应力位移云图分别见图6图8。优化前后设计变量值如表1所示，表2是优化前后计算结果的比较。可以看出，优化目标函数即所有节点处的最大von Mises等效应力，由初始的1272 MPa降低到1103MPa，下降率为13.29%；优化效果很明显。优化后的最大应力点的应力计算值为1103MPa，局部最大应力区域的应力水平约为900MPa。图6 优化后叶轮von Mises应力云图图7 优化后叶轮径向位移云图图8 优化后叶轮轴向位移云图 表1 初始设计和优化设计中设计变量值 mm 坐标BCDEFGt初始设计X214.38217.56213.89210.22207.17206.198Y304.33259.11226.11204.84185.53168.67Z234.67210.22185.78161.33136.89112.44优化设计X215.43219.00214.94211.14208.35 206.079.3804Y310.48263.69230.04207.46185.32 168.47Z237.68213.36189.05164.34136.58112.72 表2 优化前后计算结果比较初始设计优化设计减少百分比/%最大应力/MPa1272110313.29叶轮最大径向位移/mm3.4862.73721.49叶轮最大轴向位移/mm1.9441.72211.42 4 结论以有限元分析软件ANSYS为工作平台，使用其二次开发语言APDL语言实现了叶轮的三维参数建模。在参数化建模的基础上，对叶轮进口叶片的前缘进行了形状优化设计，以叶轮进口叶片前缘的几何尺寸为设计变量，以降低最大应力为目标，使叶轮的质量以及变形尽可能地小，大大减少了工作量，提高了工作效率。参 考 文 献1 刘士学，方先清.透平压缩机强度与振动M.北京: 机械工业出版社, 1997.2 蔡显新，尹泽勇，高德平.整体离心叶轮的形状优化设计J .航空动力学报,2000,15(2):219-221. 3 邱凯,徐自力,王尚锦,等.离心压缩机叶轮三维有限元强度分析系统J .应用力学学报,19