基于有限元分析的大型水轮机轴联销钉拔销强度及寿命预测方法.pdf

第4 6 卷第1 1 期2023年1 1 月水电站机电技术Mechanical&Electrical Technique of Hydropower StationVol.46 No.11Nov.20231基于有限元分析的大型水轮机轴联销钉拔销强度及寿命预测方法钱祥1，何亚文1，张宝庆1，李友才1，吴秩嵩1，周海林1，黄俊玮，刘桂英（1.湖南五凌电力工程有限公司，湖南长沙4 1 0 0 0 4；2.长沙理工大学电气与信息工程学院，湖南长沙4 1 0 1 1 4）摘要：销钉拔销系统是大型水轮机轴联检修中拆解联轴的必要设备，其强度和寿命影响使用过程的安全。本文针对大型水轮机轴联销钉拔销系统强度和寿命预测问题，提出了一种基于有限元分析的销钉拔销系统强度和寿命预测新方法。首先，基于有限元分析构建起销钉拔销系统有限元模型；然后，根据各有限单元的应力和应变分布特征，结合SWT模型，建立起销钉拔销系统的强度与寿命预测模型。在销钉拔销系统有限元模型中，简化了非线性预测问题，将其转化为线性预测问题，以提高基于有限元模型的强度与寿命预测模型的求解速度。最后，通过计算机仿真，验证了所提出的预测模型可以实现快速、高精度的大型水轮机销钉拔销系统强度与寿命预测。关键词：销钉；拔销；轴联；有限元；强度；预测中图分类号：TK730.2DOl:10.13599/ki.11-5130.2023.11.001文献标识码：A文章编号：1 6 7 2-5 38 7(2 0 2 3)1 1-0 0 0 1-0 5分布的方式1 0.。在工程应用可靠性分析中,该分0引言布在产品寿命以及可靠性分析预测中应用广泛2 。销钉拔销系统是大型水轮机轴联检修中拆解联轴的必要设备，它的状态对大型水轮机的安全检修起着举足轻重的作用,2 1。疲劳失效是机械元件失效中的一种主要形式，销钉拔销系统断裂主要发生在高应力区域1 3.4。为了提高大型水轮机检修拔出销钉的安全性，降低机械元件的运维成本，减少检修事故的发生次数，对于销钉拔销系统强度疲劳程度及寿命的准确预测就具备了至关重要的工程意义。已有大量学者对于销钉损伤程度及寿命预测进行了研究,姜丰等通过计算不同负荷条件下剩余寿命的分布情况，建立了载荷数疲劳寿命的干涉模型。刘俭辉等口采用剩余强度作为控制参数的传统应力强度干涉模型来对机械元件的疲劳可靠性进行分析。余天堂等8 提出了一种含有位置参数的对数正态分布来进行销钉剩余寿命预测的方法，同时分析了在恒幅值载荷和变幅值载荷作用下销钉的剩余寿命分布形式。杨嘉顺9提出了一种在复杂应力状态下仍能适用的剩余寿命预测模型，该模型是基于临界平面法搭建的。现有的研究中,对于疲劳寿命的描述主要采用正态分布以及对数正态但是，专门针对销钉拔销系统的研究文献极少。有鉴于此，本文针对销钉拔销系统的强度和寿命预测问题进行研究，基本思想是：根据有限元理论构建起销钉拔销系统的有限元模型，在有限元模型基础上，根据各有限单元的应力和应变分布特征，运用Smith-Watson-Topper(SWT)模型理论，建立起包含销钉拔销系统的材料属性及其疲劳参数、损伤参数的销钉疲劳强度与剩余寿命的SWT预测模型在销钉拔销系统有限元模型中，简化了非线性预测问题,将其转化为线性预测问题,以提高基于有限元模型的强度与寿命预测模型的求解速度。1销钉拔销的有限元模型有限元分析(Finite Element Analysis，F EA)是一种利用数学近似对真实物理系统进行模拟的方法。基本思想是用较简单的问题代替复杂问题后再求解，利用简单而又相互作用的有限微小单元及其有收稿日期：2 0 2 3-0 9-0 8作者简介：钱祥（1 990-），男，工程师，从事大型水轮发电运行与故障诊断工作。2限数量的未知量去逼近无限未知量的真实系统，将整体系统的求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的近似解，然后推导求解这个域总的满足条件（如结构的平衡条件),从而得到整体系统的解1 3。销钉拔销系统的有限元建模原理如下：(1）建立单元位移模式模型将销钉拔销系统中的销钉划分为有限数量的N个微小单元u，微小单元u为刚性六面体，每个单元有8 个节点，每个节点有6 个自由度的位移,其中3个距离位移,3个相位位移。因此,第n(n=1,2,N。下同）单元u，的第i节点的未知量包括6 个应力变量、6 个应变变量、3个位移变量，分别构成应力向量、应变向量、位移向量。位移未知量也称基本未知量。第n单元u.的第i节点的应力向量可以表示为(1)第n单元u,的第i节点的应变向量可以表示为(2)第n单元u,的第i节点的位移向量可以表示为(3)第n单元的全部节点的位移向量构成的节点位移扩展向量可以表示为(4)将微小单元u看成一个质点，则第n单元u有6个应力变量、6 个应变变量、3个位移变量，因此，第n单元u.的应力向量可以表示为(5)第n单元u.的应变向量可以表示为Cu=leux Euy Cua Cuxy Cuy Cual第n单元u.的位移模式Qu表示为设第n单元的位移模式使用线性差值表示，则.可以表示为m,=0,(x,y,z)=a.jo+aj/,+amj2y.+amj3m,j=x,z(8)写成矩阵形式,有anxanx2anx3anylamy2amy3二unLanzil对于任意单元n的任意节点i的位移变量2应满足：Anij=nj(x,y,z)=anjo+amjimi+anj2/mi+anj3mij=x,J,z写成矩阵形式，有水电站机电技术anlanx2anx3Xanxo.i=anylaay2amy3Lamzlanmz2anz3/z;JanzoJ其中，Xai、Ja i v、Za 为位移向量a在xyz坐标系中的x、y、z 轴上的坐标。根据式(7)(1 1),可以得到,=Aa.其中，A00A200A：0 0A=0A100ALO0A100称为形态函数(形函数)矩阵。(2)建立单元应力向量与应变向量关系模型第n单元u,的单元应力向量与应变向量关系.可以写成.-Deu.式(1 3)中的矩阵D为1-u1-u10001-u1-u1000E(1-)1-u1-uD=(1+)(1-2/1-2u0002(1-M)1-2u00 002(1-)1-2000002(1-M)其中，E为杨氏弹性模量，为柏松比。(6)(3)建立单元应变向量与位移向量关系模型第n单元u,的单元应变向量与节点位移扩展(7)向量关系可以写成u.=BA,式(1 5)中,B=Bn1A00axa.nxoayoanz2anz3Jz第4 6 卷aayo,i=1,2,8(11)(12)0 A:0.00AJ0Bm2.Bmsl,其中,0aA.0ay(9)OA.00B.=QAA0yaxaAA0(10)A0ALozay(13)(14)00(15)i=1,2,8。ay第1 1 期钱祥，等：基于有限元分析的大型水轮机轴联销钉拔销强度及寿命预测方法3(4)建立单元刚度矩阵模型与载荷向量模型根据式(1 3)、(1 5),可以得到单元应力向量u.与节点位移扩展向量的关系为.u.=Deu.=DBa,=San其中,S为应力矩阵,S=DB。基于虚功原理，可得到单元的节点力扩展向量（即为单元载荷向量）与节点位移扩展向量的关系为F,=K,an(17)其中，F，也称为节点力扩展向量，F,=F.1 Fa2式（1 7)也称单元特性方程；K,称为单元刚度矩阵,K,=Jv,BDBdV=BDBV。V,表示单元的体积。(5)建立整体的总刚度矩阵模型与总载荷向量模型整体结构的总刚度矩阵的总载荷向量分别是单元刚度矩阵和单元节点力扩展向量的叠加，因此有NNK=Z Kn,F=Z F.n=1F=Kan,n=1,2,N如果单元载荷比较复杂，按如下表达式求取等效载荷作为单元载荷向量。F,=AF.+Js,AF,ds=Jv.AFvdV其中，F。F、F v 分别表示作用于单元上的集中力、面力和体力。(6)建立节点载荷平衡方程在整体结构中，一个节点为几个单元所共有,因此,对每个节点i,均有E.Tai=R,或,2.Fri=0其中，u表示共有节点i(与节点i相连）的单元的集合，Ti、R,分别表示第n单元的第i节点（假设共有的节点在单元中的编号相同）的节点力（也称内载荷）和外力（也称外载荷）,因此,前述的F是T(内力)和R(外力)之和。(7)建立边界条件1)力条件边界上的单位面积的外载荷的面力R。等于内载荷的面力T。2)位移条件位移条件也称几何条件，在边界上的单元的节点位移已知为元，则入.=元m。(8)计算节点应力根据节点载荷平衡方程和边界条件，首先，求解式(1 7)得到位移a，然后，根据式(1 8)计算各单元的节点力扩展向量F,，根据式(1 6),计算各单元的单元应力向量u，从而预测出销钉拔销的受力分布。(16)2基于有限分析模型的寿命预测SWT模型销钉拔销系统的疲劳主要是以销钉的微动疲劳为主。微动疲劳是一种常见的接触疲劳现象，它发生在两个表面紧密配合并有微小的相对运动时。这种运动会在接触区域产生高应力和高应变，从而引起表面的损伤和裂纹。销钉拔销系统作为销钉拔出过程使用的一种重要器件,反复使用受到微动疲劳的影响,其寿命会因此降低。本文采用临界平面法中的SWT模型来建立销钉拔销疲劳受损的剩余寿命预测模型0 。临界平面法是一种分析裂纹萌生的方法4,基本思想是：将裂纹萌生问题转化为寻找临界平面和综合参数的问题,在复杂应力状态下,存在一个最有利于裂纹萌生的平面,称为临界平面，在有限分析模n=1型基础上，通过在有限单元积分点计算一个综合参(18)数，可以比较不同截面的裂纹萌生倾向，从而预测工程构件的疲劳寿命。2.1基于有限分析模型的微动损伤预测模型(19)裂纹的萌生通常发生在一个特定的危险平面上，这个平面可以用一个综合参数来描述，称为临界平面参数。临界平面参数可以反映材料在复杂应力状态下的疲劳性能,有多种形式,如SWT模型、FP模型、SSR模型、MSSR模型等1 5 。本文采用SWT模型建立销钉拔销系统微动损伤预测模型，如式(20)(21)。根据销钉拔销系统的销钉微动损伤度，从而预测出销钉拔销系统的耐力强度。D=omaxca式中：D为微动损伤参量，max为最大应力，ea为应变幅值。2.2基于微动损伤模型的疲劳寿命预测模型微动疲劳寿命预测模型的核心在于解决接触问题1 3。接触表面的应力和位移是微动疲劳的两个主要控制因素。为了提高应力和位移计算的精度，本文采用有限元分析模型进行计算，然后根据临界平面法的SWT模型，建立了微动疲劳寿命预测模型。微动损伤疲劳是一种由于微小的位移或变形而导致材料失效的现象。微动损伤疲劳的影响因素众多，包括材料的性质、微动损伤幅度、作用力幅值及(21)4其频率、温度、环境等。基于微动损伤的疲劳寿命预测模型可以表示为：T=AVDf式中:A和m为材料特性相关的系数,8 为应力，厂为力的频率。2.3强度与寿命预测流程销钉拔销系统的强度与寿命预测流程如图1 所示，通过建立销钉拔销系统的有限元分析模型，计算出销钉拔销系统的节点力和单元应力分布，在此基础上，计算销钉拔销系统微动损伤度和疲劳寿命。开始建立单元节点力和整体节点力计算模型选择单元类型和单元划分数量建立节点应力、应变、位移变量建立单元里建立位移模式构建单元应变矩阵和应力矩阵建立单元和整体刚矩阵计算模型图1 强度与寿命预测流程图3仿真与试验验证3.1拔销过程的销钉应力分布仿真现在常用的有限元软件有Ansys、A b a q u s、Cosmos和Marc等。Abaqus是一款专业的有限元分析软件，在国际上的知名度都是首屈一指的，完美地融合了同类软件中的优点，拥有各种类型的材料模型库，可以模拟典型工程材料的性能,适用于相对简单的线性分析和复杂的非线性分析。本文使用Abaqus2022对销钉拔销系统的销钉拔销过程进行有限元分析仿真计算。仿真的销钉拔销系统的结构示意图如图2 所示，销钉是构成销钉拔销系统的核心元件，销钉为圆柱体金属构件，仿真的大型水轮机轴联的销钉拔销水电站机电技术系统的销钉为6 个，呈等间隔圆环设置，使用6 个销钉联接的大型水轮机轴联系统的Abaqus2022仿真模型结构图如图3所示，仿真的销钉拔销系统的销(22)钉长度4 0 0 mm，直径1 2 0 mm，弹性模量为2 0 6 MPa，泊松比为0.3,销钉拔销系统在拔销过程中的外载荷(外力)90 0 MPa,属于面力。销钉拔销系统拔销过程的销钉云图如图4 所示。由图4 可以看出，在销钉的4 处出现较高应力，此应力集中现象与环带的微动疲劳问题有关，符合实际构件破坏情况，建立节点力平衡方程式和边界条件计算各单元的节点位移、节点力和应力计算微动损伤度计算疲劳寿命输出第4 6 卷251、销钉；2、注油孔；3、高压接头；4、高压油泵管路；5、销孔；6、油泵图2 销钉拔销系统的结构示意图图3销钉联接的大型水轮机轴联系统的Abaqus模型图4 销钉拔销过程的销钉应力分布仿真模型3.2销钉拔销寿命的仿真和实验销钉拔销系统仿真与实验的销钉采用锻钢35CrMo材质，销钉的化学成分及机械性能参数见表1。销钉拔销系统仿真，采用Abaqus2022+Maltlab仿真软件。销钉拔销系统的实验过程，采用电子式万能试验机模拟销钉拔销，采用直读光谱仪进行监测，试验在室温空气介质中进行，仿真和实验结果如图5和图6 所示。图5 为销钉拔销系统的销钉SWT参数与疲劳寿命的关系曲线,图6 为销钉疲劳寿命的仿真和实验预测结果。由图5 可以看出，SWT参数的销钉疲劳数据拟合的效果较好。由图6 可以看出，实验预测结果基本上均布在试验结果的两侧,效果较为理想，销钉疲第1 1 期劳寿命实验与仿真预测的误差因子小于1.5，仿真预测结果小于实验结果，仿真预测结果偏于保守，但是与实验结果的偏差较小，由此可见，本文所研究的寿命预测模型具有较高的疲劳寿命预测精度。表1 销钉的化学成分与机械性能参数化学成分1%力学性能参数C0.32Si0.24Mn0.44P0.016S0.002Cr0.89Mo0.17Cu0.006100sMS10102图5 销钉SWT参数与疲劳寿命的关系106104103102102图6 销钉疲劳寿命的仿真和实验预测结果4结论本文旨在利用有限元法和SWT对大型水轮机轴联销钉拔销系统的强度和寿命进行预测。通过理论分析和仿真与实验，可以得到以下结论：(1)基于有限元分析，可以确定销钉拔销系统的销钉是否在承受负载时会发生塑性变形或破坏。如果超过了材料的屈服强度，销钉拔销系统可能会发生塑性变形。如果应力达到了销钉拔销系统的破坏强度，销钉可能会发生破坏。(2)建立的销钉拔销系统疲劳寿命预测SWT模型，在估算疲劳寿命时，以临界平面上的损伤为主，钱祥，等：基于有限元分析的大型水轮机轴联销钉拔销强度及寿命预测方法的预测精度,误差因子小于1.5。(3)基于有限元法和SWT的预测方法，可以较Rro2/MPa728高精度地评估销钉拔销系统的强度和寿命。该方法R./MPa899A/%12.0ZI%59.0HBW262KV./J69103104预测寿命（循环次数）10310实验寿命（循环次数）5同时了考虑材料的总体损伤程度，适用于变幅或随机多轴载荷情况下的多损伤平面情况的寿命预测。以临界平面上最大剪应变幅、正应变幅为参数建立微动损伤和疲劳寿命预测SWT模型，能够得到较高为大型水轮机轴联销钉拔销设计和工程实践提供了重要的技术支持,能够提高水轮机的可靠性,降低运维成本，并减少检修事故的发生。参考文献：1熊联英.浅谈水轮机振动故障诊断分析方法.水电站机电技术,2 0 2 1,4 4(4):1 0-1 2.2陈炼红.水轮机剪断销剪断处置流程的思考.水电站机电技术,2 0 2 0,4 3(1 1):1-2.3高婧，罗城.基于SWT模型和威布尔分布的CFRP环带微动疲劳寿命预测.材料导报，2 0 2 1,35(1 6)：1 6 0 1 5-1 6 0 2 0.4 JI S Y,MOU Q Y,YUAN G M,et al.Dowel-bearingbehavior of bamboo scrimber for bolted-type joint.In-10510O模型预测结果实验预测结果105dustrial Crops&Products,2023:193.5刘宇航,石宇强,王俊佳.基于FCM-LSTM 的滚动轴承多阶段寿命预测.机械设计，2 0 2 3,4 0(5)：4 3-5 0.6姜丰，吴华丽.某型地铁车辆车下帽型梁裂纹扩展寿命预测.现代制造技术与装备，2 0 1 7(7)：35-38.7刘俭辉，薛文卓，魏泰.含冲击损伤复合材料层合板疲劳寿命预测D).西安理工大学学报，2 0 2 2,38(2)：1 6 4-1 7 0.8余天堂，万林林，任青文.不连续岩体的扩展有限单元法分析.水力发电学报,2 0 1 0,2 9(5):92-97.9杨嘉顺.基于概率晶体塑性有限元方法的疲劳弹性寿命预测研究D.哈尔滨：哈尔滨工业大学,2 0 2 110KARGE C,STROETMANN R,WISNEWSKI J.106Experimental and numerical investigation of dowelstrips forlongitudinal and transverse loadingJ.ce/papers,2023,6(1).11GEORGIOS C,JOSEF H,MARTIN C.Multi-spancomposite girders with composite dowels-Experimentalinvestigations and design approach to account for thebeneficial effect of shear force redistribution along thecomposite interfaceJ.Engineering Structures,2023:285.12赵文峻.基于有限元方法的某变速器箱体的疲劳寿命预测与研究D.长春：吉林大学,2 0 1 2.13武洋洋，王莹.螺栓连接件微动疲劳特性分析.力学季刊,2 0 2 1,4 2(4):6 7 4-6 8 4.14周维,刘义伦,李松柏，等，一种基于临界平面法的多轴疲劳寿命预测模型.材料导报,2 0 1 5,2 9(1 6)：1 4 7-1 5 0.15张世越,吴昊.基于顺序雨流法和临界平面法的多轴变幅低周疲劳寿命评价D).力学季刊,2 0 2 0,4 1(3)：4 6 5-4 7 6.(Vol.46 No.11)QIAN Xiang,HE Yawen,ZHANG Baoqing,LI Youcai,WU Zhisong,ZHOU Hai-lin,HUANG Junwei,LIU Gui-ying2.School of Electrical and Information Engineering of Changsha University of Science&Technology,Changsha 410114,China)Abstract:The pin pulling system is essential for disassembling the coupling in maintenance of the shaft coupling fora large hydraulic turbine,and its strength and service life affect the safety during its service.This paper proposes a newmethod of predicting the system strength and life based on the finite element analysis.Firstly,a finite element model ofthe system is established based on the finite element analysis.Next,a strength and life prediction model of the system isestablished based on the stress and strain distribution characteristics of each finite element and in combination with theSWT model.In the finite element model of the system,the nonlinear prediction problem is simplified and transformed into alinear prediction problem to improve the solving speed of the strength and life prediction model based on the finite elementmodel.Finally,it is verified through computer simulation that the prediction model proposed in this paper can realize fastand high-precision strength and life prediction of the pin pulling system for the large hydraulic turbine.Keywords:pin;pin pull;shaft coupling;finite element;strength;predictionWANG Xiao-lan,HU Si-yu,GUAN Yu-yao,LIU Yang-cheng,WEI Zong-kai(Datang Hydropower Science&Technology Research Institute Co.,Ltd,Nanning 530007,China)Abstract:The safe and stable operation of hydraulic turbine-generator units is very important for hydropower stations.For multiple factors,such as mechanical,hydraulic and electromagnetic factors and aging of mechanical components,mostfaults of these units are manifested in the form of runout which can thus directly characterize the operating state of units.Inthis paper,a runout analysis method based on multi-model fusion is studied.The model integrates three machine learningalgorithms:SVR,LightGBM and XGBoost through the least square method,and the R2 coefficient of determination is usedto evaluate the model.Experiments show that the fusion model can stably predict the runout value and get a R2 score above0.98,more than 0.3 higher than that of a single model.This validates the effectiveness of the model.Keywords:hydraulic turbine;runout analysis;fusion model;machine learning(1.Guangzhou Branch of Guangdong Yunshun Integrated Energy Technology Co.,Ltd.,Guangzhou 510540,China;Abstract:The stability test is carried out for Unit 1 of a power station under four typical heads to test the vibration,swingand water pressure pulsation of the unit under no load and different loads.The causes of unit vibration are analyzed fromhydraulic,mechanical and electrical factors respectively.Based on the operation stability of the unit,the operating zone isdivided into vibration zone under low load,working condition zone under vortex strip and stable operating zone under highload.This provides accurate basis for safe operation,technical transformation and scientific dispatching of the unit.Keywords:Francis turbine generator unit;stability;vibration;operating zoneXIN Jian-zhong,MA Dai-bin(Dongfeng Electric Machinery Co.,Ltd.,Deyang 618000,China)Abstract:This paper elaborates the key points and difficulties about adjustment,fixing,welding and quality control of thenew stator elastic element during its installation by taking the installation and quality control of a new type of stator elasticelement for a 350MW pumped storage unit in China as an example and analyzing the basic principles of the element.Incombination with the magnetization test of stator core and the measured data of unit operation after the installation,thefeasibility and effectiveness of the installation process of the element are further verified.This provides a reference for theinstallation of elastic elements in similar units.Keywords:pumped storage power station;generator motor;stator;new elastic elementMechanical&Electrical Technique of Hydropower StationFEA-based pin pull strength and life prediction of shaft coupling for large hydraulic turbine(1.Hunan Wuling Power Corporation Co.,Ltd.;Changsha 410114,China;Multi-model fusion-basedrunout rationalityanalysis and prediction of hydraulic turbineStability analysis andoperating zone division of Francis turbine generator unit2.Southern Power Grid Technology Co.,Ltd.,Guangzhou 510080,China)Installation of new stator elastic element for generator motor in pumped storage power station131ABSTRACTSCAO Bin,YAO Ze,LI Xi?