【发明公布】一种卧式半调节轴流泵叶轮叶片固定螺栓疲劳寿命预测方法_扬州大学_201911105891.6 

申请/专利权人:扬州大学

申请日:2019-11-13

发明/设计人:仇宝云;苏保铭;周钰林;张金龙;盛林华;黄先北;董媛媛;刘木秀

公开(公告)日:2020-07-24

代理机构:扬州苏中专利事务所(普通合伙)

公开(公告)号:CN111444639A

代理人:许必元

主分类号:G06F30/23(20200101)

地址:225009 江苏省扬州市大学南路88号

分类号:G06F30/23(20200101);G01M13/00(20190101);G06F111/10(20200101);G06F119/14(20200101)

优先权:

专利状态码:在审-公开

法律状态:2020.07.24#公开

摘要:一种卧式半调节轴流泵叶轮叶片固定螺栓疲劳寿命预测方法,属于动力机械设备可靠耐久性技术领域,当卧式半调节轴流泵运行时,水泵叶轮叶片在受到水的周向阻力、径向离心力、轴向水压力以及叶片重力的共同作用下,以预测叶片与轮毂的固定螺栓疲劳寿命。包括泵装置内流场数值模拟、叶轮叶片受力计算;根据叶片固定螺栓、叶片根部法兰、叶轮轮毂及压板变形协调条件,计算叶片固定螺栓危险截面应力,确定疲劳载荷谱;最终计算确定各固定螺栓疲劳寿命。本发明能够准确预测卧式半调节轴流泵叶轮叶片固定螺栓的疲劳寿命,对改进叶片根部结构与连接固定设计,合理确定固定螺栓更换周期,提高水泵运行可靠性,具有重要的理论意义和实际应用价值。

主权项:1.一种卧式半调节轴流泵叶轮叶片固定螺栓疲劳寿命预测方法,其特征在于,包括如下操作步骤:A.卧式半调节轴流泵装置流场计算区域三维建模与网格划分;B.卧式半调节轴流泵装置流场数值模拟与叶轮叶片受力计算确定;C.叶片固定螺栓危险截面应力计算;D.叶片固定螺栓疲劳载荷谱确定;E.叶片固定螺栓应力循环S-N曲线计算确定;F.叶片固定螺栓应力循环S-N曲线修正;G.叶片固定螺栓疲劳寿命计算;步骤A中所述卧式半调节轴流泵装置流场计算区域三维建模与网格划分的方法如下:通过UG对泵装置流场计算区域进行三维建模,包括进水流道、叶轮、导叶、出水弯管、出水流道,叶轮部分构建不同旋转角度的叶轮;应用ANSYSICEMCFD对各部分进行网格划分,并在CFX–Pre进行拼接;步骤B中所述卧式半调节轴流泵装置流场数值模拟与叶轮叶片受力计算确定的方法如下:流动控制方程采用连续方程和动量方程: 式中:ρ—水体密度;t—时间;U—速度矢量;B—体积力总和;μeff—有效粘度;p'—修正压力;—散度;·—矢量积;—叉乘,T—转置;其中:μeff=μ+μt3 式中:μt—湍流粘度;k—湍动能;p—压力;μ—水的动力粘度;k-ε湍流模型假设湍流粘度与湍动能和湍动能耗散有关,即 式中:ε—湍动能耗散率;Cμ—k-ε湍流模型常数;采用k-ε紊流模型: 式中:σk、σε—常数;Pk—湍流剪切产出项;边界条件:水泵装置流场计算区域采用压力进口边界条件,质量流量出口边界条件;叶片受力理论分析计算确定:水泵运行时,叶轮叶片受到水流作用力、叶片重力及离心力,其中应用CFX流体计算软件计算水流对叶片的作用力;叶片重力G作用于重心,方向垂直向下;重力分量为: 式中,Gu—重力的周向分量,Gr—重力的径向分量,θ—顺水流方向观察叶轮逆时针转动相对于x轴正向转过的角度;叶片质量惯性离心力: 式中,G—叶片重量;RG—叶片重心所在半径;叶轮叶片处于不同转角和出口两导叶叶片间不同位置角下的横断面水平方向x轴、水泵轴向y轴、横断面竖直方向z轴各方向的水体作用力Fx,Fy,Fz,应用CFX流体计算软件计算得到;通过式8、9计算出叶片的重力分量、离心力,然后应用式10计算出叶片受到的轴向力Fa、径向力Fu和周向力Fr; 步骤C中所述叶片固定螺栓危险截面应力计算的方法如下:在叶片根部法兰处以其中心、即叶片轴中心Oˊ为原点建立相对坐标系rua,a轴与水泵轴线平行、正向垂直于纸面向内;通过计算相对坐标系中法兰受到的对于各轴的力矩Mr,Mu,Ma,判断出受拉力最大的螺栓; 根据水泵叶轮的转动方向,可以分析得出在相对坐标系rua中径向力Fr、周向力Fu大于零,轴向力Fa小于零,将叶片质心力系作用点C坐标l,0,0代入式11,得到Mr=0,即无扭矩影响;Mu,Ma为正值;叶片在所受力矩Ma、Mu合成的合力矩M的作用下,叶片法兰受到的固定螺栓1~6的拉力与B点轮毂对叶片法兰的向外的反力组成反力矩,该反力矩与主动力作用的合力矩M平衡;过点B作叶片法兰切线b,分析计算叶片受主动力的周向力、径向力、轴向力分量对切线b的合力矩Mb;在坐标系uO’a中,第i只固定螺栓的坐标为ui,ai,叶片根部法兰受轮毂反力的作用点B点坐标为uB,aB,其中第i只固定螺栓所在圆心角αi以u轴正向为起点逆时针转过为正,坐标由下列各式计算得到: 在坐标系中,法兰轮廓为圆,其方程为u2+a2=R12,根据点B坐标得到通过B点的切线b的直线方程:uBu+aBa-R12=014根据点到直线距离公式,求出第i只固定螺栓到切线b的距离: 式中,lib—编号为i的固定螺栓到轴b的距离;当叶片不受主动力作用时,法兰上无工作载荷作用,法兰受到固定螺栓的预紧力F0i和轮毂的反作用力F0hubi,二力平衡;当水泵运行时,叶片还受到水流的轴向水压力、周向水阻力、旋转离心力和叶片自身重力共4个主动力的作用,法兰受到工作载荷Fri;无工作载荷时:F0i=F0hubi=Fpre16式中,F0i—编号为i的固定螺栓对法兰的作用力;F0hubi—编号为i的固定螺栓处的轮毂对法兰的作用力;Fpre—固定螺栓预紧力;当有工作载荷时,编号为i的固定螺栓处的叶片根部法兰受叶片向右的拉力Fri、固定螺栓的拉力Fi和轮毂的反作用力Fhubi,考虑法兰在三力作用下的平衡:Fhubi+Fri-Fi=017由力矩平衡原理可知,6只固定螺栓处法兰的受反力Fri对b轴的总力矩与叶片所受周向力、径向力、轴向力三个方向主动力对轴b的力矩Mb相互平衡: 考虑到螺栓、法兰、轮毂及压板的形变的影响,结合式15、16、17及18计算得到固定螺栓受拉力Fi;固定螺栓的受拉应力满足下式: 式中,S—螺栓受拉应力;A—螺栓截面积;步骤D中所述叶片固定螺栓疲劳载荷谱确定的方法如下:由于在叶轮旋转过程中后导叶对叶轮的动静干涉作用及叶片重力相对于叶片轴的作用方向的改变,水泵叶轮叶片受到的轴向水压力和周向水阻力及重力为交变荷载,叶片固定螺栓所受拉应力为交变应力,对螺栓的疲劳寿命进行预测,首先应确定其载荷谱;本发明采用雨流计数法确定载荷谱;利用ANSYSnCodeDesignLife的rainflow模块计算;雨流法的计数规则如下:1重新安排载荷历程,以最高峰值或最低谷值为起点;2雨流依次从每个峰谷的内侧向下流,在下一个谷峰处落下,直到对面有一个比其起点更高的峰值或更低的谷值停止;3当雨流遇到自上面屋顶流下的雨流时即行停止;4取出所有的全循环,并记录各自的范围或幅值和均值;幅值: 均值: 其中,Smax—应力最大值;Smin—应力最小值;Sa—应力幅值;Sm—应力均值;步骤E中所述叶片固定螺栓应力循环S-N曲线计算确定的方法如下:根据已知的材料性能估算应力循环S-N曲线;一般应力循环S-N曲线采用幂函数的双对数形式:SmN=C22对式22两边取对数: 其中,S—循环应力值;m、C—材料系数;N—循环次数;一般低周循环N1﹤104,当N1=104时,疲劳极限较高,接近于屈服强度,一般取S1=0.9σ0.2;N=Nf的疲劳极限Sf根据材料加工工艺查抗疲劳设计手册得到,水泵叶片固定螺栓材料使用的是淬火并回火的合金结构钢;S1=0.9σ0.224σ-1=0.383σb+94MPa25Sf=σ-126 式中,σ0.2—材料屈服强度残余塑性形变等于原长度的0.2%时的应力;σb—材料抗拉强度;σ-1—对称循环的循环次数为Nf的疲劳极限;B—高周循环S-N对数曲线斜率;步骤F中所述叶片固定螺栓应力循环S-N曲线修正的方法如下:材料的S-N曲线只能代表光滑试样的疲劳性能,以合金钢为材料的固定螺栓的尺寸、形状以及表面情况会对S-N曲线造成影响,由于固定螺栓的疲劳极限较高,所以应对高周循环N104的S-N曲线进行修正;查阅抗疲劳设计手册得到表面加工系数β1,螺纹联接部分的疲劳缺口系数kσ,螺纹联接的尺寸系数εσ,最后得到S-N曲线修正系数kσD: 得到修正系数kσD后,对循环次数为Nf时的疲劳极限进行修正: 同时考虑到平均应力Sm的影响,采用Goodman修正方法对应力幅值Sa进行修正: 式中,Sa-1—对称应力循环下的应力幅值;Su—抗拉强度;步骤G中所述叶片固定螺栓疲劳寿命计算的方法如下:计算螺栓疲劳寿命时,水泵运行时,叶轮转速较高,叶轮叶片流动受导叶干扰的频率高,叶片每转1圈为1个应力循环,叶片固定螺栓应力循环次数高,远远大于Nf,所以应该将修正后的应力循环S-N曲线的双对数形式转换为幂函数形式,即由式23变为式22;将修正的应力幅值Sa-1代入式22得到各应力水平下应力循环次数Nj,各应力水平下的疲劳损伤Dj=1Nj;根据Miner法则的线性累积损伤理论,得到螺栓的总应力循环次数Nsum: 最终螺栓的疲劳寿命由式32计算: 式中,T—疲劳寿命,h;n—水泵转速,rmin。

全文数据:

权利要求:

百度查询: 扬州大学 一种卧式半调节轴流泵叶轮叶片固定螺栓疲劳寿命预测方法