多级注水泵首级叶轮应力分析_毕业论文(编辑修改稿)

多级注水泵首级叶轮应力分析_毕业论文(编辑修改稿)内容摘要：

，周玉丰 [5]以 LB50— 160型离心泵的叶轮为分析对象，利用 ANSYS软件对叶轮在流道流体压力和离心惯性力共同作用下进行应力结果分析，得出应力分布结果及规律，为叶轮盖板的强度计算提供了可靠依据，验证了有限元建模方法和计算方法的正确性，为进行不同类型的有限元分析和优化设计提供了参考。 2020年，王洋 [6]等在 ANSYS Workbench软件中 ，采用单向流固耦合方法对 冲压焊接叶轮在流场中的应力及变形情况 进行仿真计算。 结果表明，各工况流量下， 叶 轮应力分布明显不均，并在局部出现应力集中。 叶轮变形的总位移随半径的增大不断变大，并在叶轮边缘达到最大值。 根据仿真模拟结果对提高叶轮可靠性问题提出解决方案，计算结果为冲压焊接叶轮的结构设计及分析提供有效依据。 有限元分析方法 有限元分析是使用有限元方法来分析静态或动态的物理物体或物理系统。 在这种方法中一个物体或系统被分解为由多个相互联结的、简单、独立的点组成的几何模型。 在这种方法中这些独立的点的数量是有限的，因此被称为有限元。 由实际的物理模型中推导出来得平衡方程式被使用到每个点上，由此产生了一个方程组。 这个方程组可以用线性代数的方法来求解。 有限元分析的精确度无法无限提高。 元的数目到达一定高度后解的精确度不再提高，只有计算时间不断提高 [7]。 在实践中，有限元分析法通常由三个主要步骤组成 [8]： 前处理 ：用户需建立物体待分析部分的 模型 ，在此模型中，该部分的几何形状被分割成若干个离散的子区域 —— 或称为 “ 单元 ”。 各单元在一些称为 “ 结点 ” 的离散点上相互连接。 这些结点中有的有固定的位移，而其余的有给定的载荷。 准备这样的模第 1章 绪论 3 型可能极其耗费时间，所以商用程序之间的相互竞争就在于：如何用最友好的图形化界面的 “ 前处理 模 块 ” ，来帮助用户完成这项繁琐乏味的工作。 有些预处理模块作为计算机化的画图和设计过程的组成部分，可在先前存在的 CAD 文件中覆盖网格，因而可以方便地完成有限元分析。 分析 ：把预处理模块准备好的数据输入到有限元程序中，从而构成并求解用线性或非线性代数方程表示的系统 ijij fuK  式中， u 和 f 分别为各结点的位移和作用的外力。 矩阵 K 的形式取决于求解问题的类 型，本模块将概述桁架与线弹性体应力分析的方法。 商用程序可能带有非常大的单元库，不同类型的单元适用于范围广泛的各类 问题。 有限元法的主要优点之一就是：许多不同类型的问题都可用相同的程序来处理，区别仅在于从单元库中指定适合于不同问题的单元类型。 后处理：分析的早期，用户需仔细地研读程序运算后产生的大量数字，即 列出的模型内各离散位置处的位移和应力。 这种方法容易漏掉重要的趋向与热点，而最新的程序则利用图形显示来帮助用户直接观察运算结果。 典型的后处理模块能显示遍布于模型上的彩色等应力线图，以表示不同的应力水平，显示的整个应力场的图像类似于光弹性法或云纹法的实验结果。 ANSYS 软件简介 ANSYS 软件 概 述 ANSYS 软件 是融结构、流体、电场、磁场、声场分析于一体的大型通用 有限元分析软件。 由世界上最大的有限元分析软件公司之一的美国 ANSYS开发，它能与多数 CAD软件接口 ，实现数据的共享和交换，如 Pro/Engineer, NASTRAN, Alogor, I－ DEAS, AutoCAD 等， 是现代产品设计中的高级 CAE 工具之一。 ANSYS 软件 构成 ANSYS 软件 主要包括三个部分 [9]：前处理模块，分析计算模块和后处理模块。 其中 前处理模块提供了一个强大的实体建模及 网格 划分工具，用户可以方便地构 造有限元模型 ； 分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透（ 11） 第 1章 绪论 4 明显示（可看到结构内部）等图形方式显示出来，也可将计算结果以图表、 曲线 形式 显示或输出。 软件提供了 100 种以上的单元类型，用来模拟工程中的各种结构和材料。 ANSYS 软件 提供的分析类型 结构静力分析 ： 用来求解外载荷引起的位移、应力和力。 静力分析很适合求解惯性和阻尼对结构的影响并不显著的问题。 ANSYS 程序中的静力分析不仅可以进行线性分析，而且也可以进行非线性分析，如塑性、蠕变、膨胀、大变形、大应变及接触分析。 结构动力学分析 ： 结构动力学分析用来求解随时间变化的载荷对结构或部件的影响。 与静力分析不同，动力分析要考虑随时间变化的力载荷以及它对阻尼和惯性的影响。 ANSYS 可进行的结构动力学分析类型包括：瞬态动力学分析、模态分析、谐波响应分析及随机振动响应分析。 结构非线性分析 ： 结构非线性导致结构或部件的响应随外载荷不成比例变化。 ANSYS 程序可求解静态和瞬态非线性问题，包括材料非线性、几何非线性和单元非线性三种。 动力学分析 ： ANSYS 程序可以分析大型三维柔体运动。 当运动的积累影响起主要作用时，可使用这些功能分析复杂结构在空间中的 运动特性 ，并确定结构中由此产生的应力、应变和变形。 热分析 ： 程序可处理热传递的三种基本类型：传导、对流和辐射。 热传递的三种类型均可进行稳态和瞬态、线性和非线性分析。 热分析还具有可以模拟材料固化和熔解过程的相变分析能力以及模拟热与结构应力之间的热－结构耦合分析能力。 电磁场分析 ： 主要用于电磁场问题的分析，如电感、电容、磁通量密度、涡流、电场分布、磁力线分布、力、运动效应、电路和能量损失等。 还可用于螺线管、调节器、发电机、变换器、磁体、加速器、电解槽及无损检测装置等的设计和分析 领域。 流体动力学分析 ： ANSYS 流体单元能进行流体动力学分析，分析类型可以为瞬态或稳态。 分析结果可以是每个节点的压力和通过每个单元的流率。 并且可以利用后处理功能产生压力、流率和温度分布的图形显示。 另外，还可以使用三维表面效应单元和热－流管单元模拟结构的流体绕流并包括对流换热效应。 声场分析 程序的声学功能用来研究在含有流体的介质中声波的传播，或分析浸在流体中的固体结构的动态特性。 这些功能可用来确定音响话筒的频率响应，研究音乐大厅的声场强度分布，或预测水对振动船体的阻尼效应。 第 1章 绪论 5 压电分析 ： 用于分析二维或三维结构对 AC（交流）、 DC（直流）或任意随时间变化的电流或机械载荷的响应。 这种分析类型可用于换热器、振荡器、谐振器、 麦克风等部件及其它电子设备的结构动态性能分析。 可进行四种类型的分析：静态分析、模态分析、谐波响应分析、瞬态响应分析 本文 研究内容 学习 ANSYS 有限元分析软件，掌握模型静力分析方法，熟练应用各种载荷加载方法，能用表面效应单元等方 法解决非法向载荷等问题 ； 在 ANSYS 软件中，分析多级注水泵首级叶轮在 离心惯性力、流道流体压力、前后盖板外侧液体压力、圆盘摩擦力等载荷单独作用及耦合载荷作用下的受力及变形情况 ； 比较各种载荷形式下多级注水泵首级叶轮受力及变形结果并进行总结； 针对应力集中等问题对多级注水泵首级叶轮 安全性 进行校核。 第 2章 多级注水泵首级叶轮三维模型结构 6 第 2章 多级注水泵首级叶轮 三维模型结构 前言 叶轮是水泵的关键过流部件，其制造质量直接影响到水泵性能的好坏以及机组的稳定性。 为了获得性能优良的叶型，传统 的方法是先根据模型换算法或速度系数法计算出流道、叶片截线及木模截线，然后根据图纸制成模型，并在试验台上反复试验，最终得出符合要求的叶片形状。 整个过程费时、费力且成本高。 为了解决这一问题，最有效的途径是对水泵叶轮叶片进行三维造型，并通过流场分析来初步得到叶轮内部的压力、速度等分布状况，做出初步的性能分析，从而为模型制造与生产节省成本和时间。 Pro/E软件 是一套由设计至生产的机械自动化软件，是新一代的产品造型系统，是一个参数化、基于特征的实体造型系统，并且具有单一数据库功能。 Pro/E软件 还提供了实体模型和薄壁 模型的有限元网格自动生成能力，也就是它自动地将实体模型划分成有限元素，以便有限元分析用，所有参数化应力和范围条件可直接在实体模型上指定，即允许设计者定义参数化载荷和边界条件，并自动生成四边形或三角形实体网格。 载荷 /边界条件与网格都直接与基础设计模型相关联，并能像设计时一样进行交互式修改 [10]。 基于 Pro/E软件 的 上述 特点，本文 所分析的多级注水泵首级叶轮三维模型的创建 在Pro/E软件 中实现。 叶轮主要参数机模型的创建 根据对叶轮实物的测绘数据，以 AutoCAD 软件为平台 ，创建叶轮二维模型 [11]， 得到叶片木模图 （如图 21） 及 叶轮结构图 （ 如图 22）。 叶轮的主要设计 参数 为： 叶片进口直径 101mm； 出口直径 308mm； 叶片进口角 10176。 ； 叶片出口角 26176。 ； 叶片厚度 6mm；叶片包角 130176。 ；叶片弦长 ；轮毂厚度 8mm；叶片数 7；扬程 122m；设计流量158m3/h；转速 2950r/min。 图 21 叶轮叶片木模图 图 22 叶轮结构图 第 2章 多级注水泵首级叶轮三维模型结构 7 利用已有的叶片水力木模图数据，结合 Pro/E软件实现叶片、叶轮及流道的三维实体造型如图 23所示。 图 23 叶轮三维模型 第 3章 各种载荷单独作用下叶轮受力分析 8 第 3章 各种载荷单独作用下叶轮受力分析 前言 多级注水泵在正常工作条件下，主要承受离心惯性力、 流道流体压力、前后盖板外侧液体压力、圆盘摩擦力等载荷作用。 在本章中，假设叶轮分别承受上述各种载荷单独作用，在 ANSYS 软件中，分析叶轮在各种载荷作用下的受力及变形情况，校验叶轮受力及设计合理性。 离心惯性力单独作用下受力及变形分析 ANSYS 软件分析过程 对叶轮 施加离心惯性力进行受力及变形分析 ，应用结构分析中的体单元 ，同时叶轮形状较复杂，难以用结构网格进行离散化分析，因此选择与非结构网格相对应的SOLID187 号单元 [12]，该单元为 10 节点四面体单元，具有 2 阶精度，适用于空间结构的弹性、塑性、超弹性、蠕变等结构分析。 在 ANSYS 中没有规定单位 [13]，需要用户自己去定义自己的单位制， 本文中单位统一标准为 ：长度 mm；时间 s；密度 t/mm3；力 N； 弹性模量 MPa；应力 Mpa； 角速度rad/s。 叶轮材质为 ZG1Cr13NiMo， 材料抗拉极限 σb 为 495MPa， 密度取为 109t/mm3，弹性模量 取为 2 105 MPa，泊松比取为。 对叶轮实体模型进行网格划分，设置单元尺寸为 5，采用自由网格划分方式，共将模型离散为单元 84831 个，节点 151260 个，网格划分结果如图 31 所示。 图 31 叶轮网格划分结果 于叶轮 轮毂内侧表面施加约束条件， 本文中 多级注水泵 首级 叶轮额定转速为ω=2900r/min=，按照 ANSYS 转向右手螺旋准则，绕轴向施加上述惯性载荷。 利用软件的通用后处理 (POST1)查看模型的有限元计算结果，求解器类型选为 PCG。 转换坐标系为极坐标系，设 置模型变形放大倍数为 3247倍，查看叶轮在离 心惯性力作用第 3章 各种载荷单独作用下叶轮受力分析 9 下的受力及变形情况。 应力结果分析 径向应力 图 32 为在离心惯性力 单独 作用下的叶轮径向应力分布云图。 从图中可以看出，在离心惯性力单独作用下，叶轮的径向应力主要集中在靠近叶轮入口部分叶片吸力面与后盖板内侧交界处，并由此向叶轮外缘逐渐变小。 最大径向应力为 MPa，远小于叶轮材料的强度极限。 图 32 径向应力分布云图 由上图可以看出， 径向应力在前、后盖板的分布不同， 其中后盖板内侧的径向应力分布明显较前盖板处相应部分大。 周向应力 图 33 为在离心惯性力作用下的叶轮 周向 应力分布云图。 周向应力主要集中在叶片压力面与叶轮前盖板内侧交界处。 最大周向应力为 MPa， 远小于叶轮材料的强度极限。 图 33 周 向应力分布云图 最大 径向应力 最大周向应力 第 3章 各种载荷单独作用下叶轮受力分析 10 由图中可以看出，最大周向应力点不在周向应力 主要分布 区域，而在叶轮外缘处叶片吸力面与后盖板内侧交界部分 ，说明此处为叶轮。