《2023年基于ANSYS的水轮机推力轴承动力学分析水轮机推力轴承.docx》由会员分享，可在线阅读，更多相关《2023年基于ANSYS的水轮机推力轴承动力学分析水轮机推力轴承.docx（16页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、2023年基于ANSYS的水轮机推力轴承动力学分析水轮机推力轴承 机械制造匡涛，等基于ANSYS 的水轮机推力轴承动力学分析 基于ANSYS 的水轮机推力轴承动力学分析 匡涛，秦战生 （河海高校能源与电气学院，江苏南京211106） 要：针对某一水轮机结构，推力轴承作为其重要部件，其运转状态的好坏对水轮机运行 的平安性和平稳性产生重大影响。因此对水轮机推力轴承进行动力学仿真及结构分析，研摘 究推力轴承运动时的特征，为该推力轴承实时监控供应理论支撑，具有非常重要的意义。介绍了该水轮机轴承内部的运动学关系，建立了该轴承的有限元模型，利用ANSYS /LSDYNA 显式动力学模块对该轴承的运动过程进
2、行了仿真，并用LS PEPOST进行后处理，分析了滚动轴承在正常状况下，轴承各部件的应力、位移分布规律，为水轮机推力轴承的正常运转供应了理论依据。 关键词：推力轴承；有限元分析；显式动力学；动态仿真 5276（2023）01-0038-04中图分类号：TH13333+1；TP3919文献标记码：A 文章编号：1671- Dynamic Analysis of Turbine Thrust Bearing Based on ANSYS KUANG Tao ，QIN Zhan-sheng （Hohai University ，College of Energy and Electrical Eng
3、ineering ，Nanjing 211106，China ） Abstract ：Thrust bearing is an important component of the turbine structure ，its quality affects the safety and stability of the turbine operationTherefore ，its dynamics simulation and structural analysis and investigation on its movement characteristics are of great
4、 significance to theoretical support for real time monitoringThis paper describes the kinematic relations inside the turbine bearing ，builds its finite element model and simulases motion of rolling bearings by using ANSYS /LSDYNA explicit dynamics moduleIt also shows the stress and displacement dist
5、ribution of the bearing operating under normal conditions by processing the result using LS PEPOSTThe theoretical basis is provided for smoothly operating the turbine Keywords ：thrust bearings ；finite element analysis ；explicit dynamics ；dynamic simulation 1 敛性问题。 ANSYS /LS DYNA 是以非线性动力分析为主，应用 0引言 在
6、水轮机运行中，轴承是平安运行的关键部件之一，其运行过程中状态的好坏对水轮机平安构成重大影响，轴承运行是否稳定关系到水力发电的平安生产。 本文在滚动轴承运动仿真的基础上，依据某实际水轮机推力轴承（推力球轴承）运行状况，采纳LS DYNA 显示动力学对水轮机推力轴承进行分析，获得该类水轮机推力轴承运行的基本特征。 显式时间积分的大型有限元分析程序在动力学有限元分析中，系统的求解方程为： M a t +Ca t +Ka t =Q t （1）式中：a t 为系统节点的加速度向量；a t 为系统节点的速 度向量；a t 为系统节点的位移向量；M 为系统的质量矩阵；C 为系统的阻尼矩阵；K 为系统的刚度矩
7、阵；Q t 为系统节点载荷向量。 在显式动力学中心差分法中，假定时间t =0的位移 a 0，速度a 0和加速度a 0已知。假设时间求解域被等分为n 等时间间隔t ，并且0t 时刻内的全部解已经求得，计算目的在于求解t +t 时刻的解。 利用泰勒绽开式，将a t +t 在时刻点t 展成泰勒多项式，并取有限项作为a t +t 的近似值： 22 t +t （p ） =+a t a a t +t a t ta t 2p t 1显式动力学基本算法 有限元软件包含显式算法与隐式算法两个核心算 法。隐式算法是基于虚功原理，每一增量步都须要对静态平衡方程进行迭代求解，并且在迭代的过程中须要求解大型线性方程组，
8、这样就占用相当多的计算资源（包括磁盘空间和内存等）；显式算法采纳动力学方程的一 线性加速度法些差分格式（如广泛运用的中心差分法、 和Newmark 法等），同时也不不须要干脆求解切线刚度， 须要进行平衡迭代，计算速度快，可节约大量计算时间，只要时间步取的足够小，对于非线性问题一般不存在收 （2） 式中：a t 为a t 的p 次微分。 依据系统求解式（1）为二次微分，取： a t +t =a t +ta t + 2 t a t 2 （p ） （3） 对式（3）进行求导，省略高阶得： 硕士探讨生，主要探讨方向为流体机械结构。男，四川绵阳人，作者简介：匡涛（1988）， 38http ：ZZHDc
9、hinajournalnetcn E-mail ：ZZHDchainajournalnetcn 机械制造与自动化 机械制造 匡涛，等基于ANSYS 的水轮机推力轴承动力学分析 a t +t =a t +t a t （4） t ）与（t ，t +t ）内，又在时间区间（t t ，速度可用位移 近似表达，有： 1 （a t a t t ）t 1 a t +t =（a t +a t +t ） t a t =（5）（6） 式（6）代入到式（3）里有：将式（5）、1 a t =（a t t 2a t +a t +t ） t t +t ）有：同样在时间区间（t t ， （7） 图1 为全约束。 推力轴承所
10、在位置如图2所示：轴承在工作中受压力 载荷和转速载荷，推力轴承承受质量约为2t ，转子转速为136r /min，压力和转速施加在相应位置的刚性面上。为模拟真实水轮机启动过程的转数和应力加载，将两参数加载历程设置如图3和图4所示 ： 轴承有限元模型图 a t = 1 （a +a ）2t t t t t （8） 式（8）代入式（1），系统求解方程可得求解将式（7）、 各个离散时间点位移值的递推公式： ( Q t K 11 +M C a t +t =2 2t t ) 211M a t M C a t t 22 2t t t 211 。即为经典的中心差分格式 ()() 2水轮机推力轴承仿真模型的建立
11、a ）有限元模型的建立 以某低水头水轮机（该水轮机不用于发电，用于回收 冷却塔能量回收）推力轴承51322为仿真对象，其结构尺寸如表1所示。采纳自下而上的方式，建立推力轴承的几何模型，并对模型做了如下简化： 1）在仿真过程中没有考虑滚动轴承径向游隙和轴向游隙以及油膜的影响。 2）由于滚动轴承在工作过程中塑性变形很小，所以仿真过程中轴承材料均为线弹性材料 表1 1，3 图2 低水头水轮机局部结构图 。 推力球轴承51322结构参数 数值13 符号物理含义数值符号物理含义数值符号物理含义d 轴承内径110mm d 1座圈内径113mm z 滚子数量 D 轴承外径190mm D 1轴圈外径187mm
12、 d z T 轴承高度63mm r 圆角半径2mm d r 滚子直径1485mm 滚道直径 15mm 图3 转数加载过程图 有限元单元选用3D Solid 164体单在分网过程中，元，采纳扫略分网和映射分网相结合的方式划分网格，采纳六面体单元。为模拟边界条件，将与其他刚结构相作用的面假设为刚性面。由于3D Solid 164单元没有旋转自由度，无法施加转速，因此将轴承轴圈刚性面设为Thin Shell 163单元以便施加径向力和转速。滚动轴承有限元模型如图1所示。 b ）边界条件、材料参数和载荷参数设置 轴圈、座圈和滚动体材料均为GCr15钢，其密度 211 为7830kg /m，弹性模量为2
13、06 10Pa ，泊松比为0 23；保持架材料为冷轧钢板，密度为7830kg /m，弹性模 11 量为196 10Pa ，泊松比为024。为了模拟轴承固 图4应力加载过程图 定在轴承座中的工况，将轴承座圈下表面的刚性面 设 Mac hine Building A utomation ，Feb 2023，44（1）：38 41 39 转数由0加载至136r /min，在05s 内完成（真实水 轮机启动过程时间较长，计算量太大，仿真效果基本一），致，在此设置为05s 。应力在整个仿真过程中不改变。 c ）接触模型的设定 在显式分析中没有接触单元，与隐式分析不同的是， 须要定义接触面以及接触类型和接
14、触有关的参数来仿真 实际接触的状况。 滚动轴承在工作中存在三种接触，即滚动体与轴圈、座圈滚道之间的接触和滚动体与保持架兜孔之间的接 1，4 。因此分析选用自动面面接触类触，均为面面接触型，依据接触面与目标面的指定原则，定义轴圈滚道表 面、座圈滚道表面和保持架兜孔表面为目标面，滚动体的外表面为接触面的39组接触对，依据材料性质等参数将滚动体与座圈、轴圈和保持架之间的静摩擦系数分别035和02，016设为035、动摩擦系数分别设为016、和01。 3轴承的显式动力学分析 图5 滚动轴承等效应力云图 依据上述有限元模型和边界、载荷条件，将计算时间取为100ms ，输出步数设为1000进行求解，采纳L
15、S PE-POST 进行后处理。计算结果如下。 1）滚动轴承元件的等效应力分布状况分析 在计算结果中任取三个时刻进行分析，为体现整个过程如图5（a ）图5（d ）所示取51322轴承各部分在较有代70ms 、85ms 和100ms 时刻的等效应力云图，表性40ms 、其 他时刻的等效应力云图类似。在40ms 时刻正常滚动轴承70ms 时刻轴承最大等效应力最大等效应力为1934MPa ， 85ms 时刻轴承最大等效应力为3499MPa ，为3368MPa ， 100ms 时刻最大等效应力为3476MPa 。比较三幅图可以看出，轴承较大的应力出现在滚动体与座圈、轴圈的接触 区域处，即图中滚动体与轴
16、圈、座圈接触区域颜色较深处。在轴承运转过程中，随着滚动体的运动最大应力值和最大应力位置都发生改变 。 以下为了能对轴承系统的应力状况进行全面的分析， 给出70ms 时滚动轴承各元件等效应力云图，其他时刻元件应力云图类似，如图6（a ）图6（d ）所示 。 40http ：ZZHDchinajournalnetcn E-mail ：ZZHDchainajournalnetcn 机械制造与自动化 些必要的简化，选取适当的单元类型和材料模型，对轴承 建立几何模型并划分网格生成有限元模型，然后依据该水轮机推力轴承的实际工作状况设定有效的约束和载荷，并设立接触对，建立好轴承各部件之间的联系，最终在AN-
17、SYS /LSDYNA 中完成了该轴承的动力学仿真，在LS PEPOST进行结果的后处理，并对仿真结果进行分析，结果表明： 1）轴承比较大的等效应力出现在滚动体与轴圈、座 图6轴承各元件等效应力云图 圈的接触区域处，应力最大值出现在接触表面以下肯定深度区域处，并渐渐向外衰减，并且球轴承接触应力在接触表面的形态为椭圆形。 2）轴承在工作过程中，应力峰值均发生在节点与轴圈、座圈接触位置，位于滚子表面。3）在轴承运转过程中，随着滚动体的运动最大应力值和最大应力位置都发生改变，接触应力，应力水平从大到小依次是滚动体、轴圈、座圈和保持架。 4）在仿真基础上，分析了轴承等效应力的改变特点。运用该方法可以较
18、精确地为正常运行的轴承供应一些数据参考，同时在此基础上可进一步作轴承的故障分析，为轴承平安运行供应更加牢靠的数据与理论。参考文献： 1高春良滚动轴承故障的显式动力学仿真与振动特征分析 D 成都：电子科技高校，2023：22-29 2张乐乐，谭南林，樊莉滚动轴承故障的显式动力学仿真与 J 上海交通高校学报，2023：1506-1509分析 3马士垚，J 机械张进国滚动轴承接触问题的有限元分析 2023设计与制造， 4樊莉机车滚动轴承运动和故障状态仿真与试验探讨D 2023：5-30北京：北京交通高校， 5万长森滚动轴承的分析方法M 北京：机械工业出版社， 1987：72-74 6沈成武滚动轴承接
19、触问题的有限元分析J 固体力学学 1982（1）：106-111报， 7李国超，彭炜，李勇才，等滚动轴承外圈故障的显式有限 2023，23（23）：元动态仿真分析J 中国机械工程，2825-2829 机8高春良，王成栋，苗强滚动轴承动力学仿真与分析J 由图5图6可以看出，滚动体、轴圈、座圈和保持架 在同一时刻最大等效应力并不相同，其中，滚动体应力最大，为3368MPa ，其次为轴圈2365MPa ，座圈2229MPa ，应力最小的为保持架618MPa 。滚动体比较大的应力出现在与轴圈、座圈接触区域处，且滚动体与轴圈接触处应力大于滚动体与座圈接触处的应力。轴圈和座圈是轴承中与滚动体干脆接触并承载
20、的主要元件，其应力分布和滚动体的应力分布有紧密的联系，因此，座圈、轴圈比较大的应力也是分布在与滚动体接触区域，应力分布区域在接触表面呈现椭圆形。保持架在各个兜孔之间的应力分布比较匀称，较大的应力出现在兜孔与滚动体接触区域。同时，由图6（a ）的应力分布可知，应力峰值均出现在滚动体与轴圈、座圈相接触位置，位于滚动体之上。 2）滚动轴承滚动体与座圈、轴圈的接触状况分析 为了清晰地体现出滚子与座圈轴圈接触状况，在LS PEPOST后处理中作出某些节点的位移曲线图，由曲线可以清晰看出在轴承转动过程中滚子与座圈轴圈的接触状况。如图7为座圈节点27793与轴圈节点7149的z 方向位移曲线 。 图7轴圈上
21、节点的z 方向位移曲线座圈、 2023：193-195械设计与制造， 9黄乾贵，邓四二，腾弘飞滚动轴承系统仿真技术的现状及发 J 轴承，2023（4）：34-37展 10高红斌基于有限元法的滚动轴承结构和模态分析与探讨 J 机械工程与自动化，2023（4）：90-92 11熊小晋，张晓，熊晓燕滚动轴承接触的非线性有限元分析 J 测试技术学报，2023，23（1）：23-2709-29收稿日期：2023- 可以清晰看出该位移呈周期性，由于模型建立时采纳 的是z 方向施加应力，对于轴圈节点有波峰为与滚子接触瞬间的位移，波谷为与滚子距离最远时的位移状况，由于运动过程中以间隙改变及应变影响，曲线整体有向下偏移的趋势。对于座圈节点的状况同样是在波峰为接触瞬间，波谷为距滚子最远时刻的状况。 4结语 首先依据仿真的实际条件对该水轮机推力轴承做一 Mac hine Building A utomation ，Feb 2023，44（1）：38 41 41