[发明专利]一种基于Grover‑Manson准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法

摘要

一种基于Grover‑Manson准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，根据抽水蓄能电厂提供实际运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算，分别得到发电电动机转子鸽尾部在发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下离心力所引起的应力随时间变化的分布情况；绘制在不同工况时最容易出现疲劳破坏点处的应力随时间变化曲线，计算此处可使用的疲劳寿命次数；将已运行m年的发电电动机实际各种工况发生次数，及各种工况下疲劳寿命使用次数代入Grover‑Manson准则计算公式，可分别考虑疲劳裂纹形成和疲劳裂纹扩展，并根据裂纹形成及扩展过程计算得到该发电电动机转子鸽尾部的疲劳寿命使用年限。

主权项

一种基于Grover‑Manson准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，其特征在于包括以下步骤：步骤1)：建立电机1/2周期模型，采用有限元法进行电磁场‑温度场‑结构场耦合数值计算，通过对电磁场控制方程(1)‑(3)和温度场控制方程(4)(5)进行有限元数值计算得到由电磁损耗引起的温度稳态分布，再由温度相比初始温度的变化对方程(6)进行求解得到热应力分布情况；▿×(1μ▿×A·)-▿(1μ▿·A·)=J·sinV2---(2)]]>Q=∫ΩJ2σdΩ---(3)]]>式中，V1是转子绕组区域，V2是定子绕组区域，σ为电导率，μ为相对磁导率，为源电流密度，Q为电磁损耗，包括源电流及涡流引起的损耗；是哈密顿算子，为矢量磁位，ω为角频率，j为虚数单位，为电位；J为总电流密度，即源电流密度和涡流密度之和；Ω代表存在电磁能量损耗的区域；∂∂x(kx∂T∂x)+∂∂y(ky∂T∂y)+∂∂z(kz∂T∂z)=-Q---(4)]]>k∂T∂n+h(T-T0)=0---(5)]]>式中，Q为能量损耗；kx,ky,kz分别表示热导率的各向异性参数；h为传热系数；T为求解温度；T0为环境温度；σij,i+Fi=0ϵij=1+vEσij-vEσkkδij+βΔT12(ui,j+uj,i)=ϵij---(6)]]>式中，i,j,k＝1,2,3；εij为应变张量；σij为应力张量；σij,j为应力张量对坐标的偏导数；E为弹性模量；ν为泊松比；β为热膨胀系数；ΔT为温度相比初始温度的变化量；Fi为外力的分量；ui,j为位移对坐标的偏导数；δij为应力因子，i＝j时为1，i≠j时为0；步骤2)：根据抽水蓄能电厂提供实际发电启动运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到发电启动工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；σij,i+Fi=ρ∂2ui(t)∂t2ϵij=1+vEσij-vEσkkδij12(ui,j+uj,i)=ϵij---(7)]]>式中，ρ为密度，ui为位移，其余参数与公式(6)相同；步骤3)：根据抽水蓄能电厂提供实际发电停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到发电停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；步骤4)：根据抽水蓄能电厂提供实际电动启动运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动启动工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；步骤5)：根据抽水蓄能电厂提供实际电动停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；步骤6)：根据抽水蓄能电厂提供实际甩负荷运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到甩负荷工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；步骤7)：根据抽水蓄能电厂提供实际飞逸运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到飞逸工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；步骤8)：根据步骤2)计算得到的发电启动工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出B点在此种工况运行情况下可使用的疲劳寿命次数N1；N=Cσ-aσ=Kσϵσβσσmax-σmin2+ψaσmax+σmin2---(8)]]>式中：C和a为材料疲劳系数；σmax为总应力变化曲线中的最大值；σmin为总应力变化曲线中的最小值；Kσ，εσ，βσ和ψa分别为有效应力集中系数、零件尺寸系数、表面系数和平均应力系数；步骤9)：根据步骤3)计算得到的发电停机工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出转子鸽尾部在此种工况运行情况下可使用的疲劳寿命次数N2；步骤10)：根据步骤4)计算得到的电动启动工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出转子鸽尾部在此种工况运行情况下可使用的疲劳寿命次数N3；步骤11)：根据步骤5)计算得到的电动停机工况时转子转子3号鸽尾B点处应力分布曲线，由公式(8)计算出转子鸽尾部在此种工况运行情况下可使用的疲劳寿命次数N4；步骤12)：根据步骤6)计算得到的甩负荷工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出转子鸽尾部在此种工况运行情况下可使用的疲劳寿命次数N5；步骤13)：根据步骤7)计算得到的飞逸工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出转子鸽尾部在此种工况运行情况下可使用的疲劳寿命次数N6；步骤14)：将步骤8)‑13)计算出的结果Ni(i＝1,2,3,…6)分别代入Grover‑Manson准则公式(9)和(10)，分别计算出发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况共六种工况下各自对应的裂纹扩展寿命ΔNi和裂纹形成寿命N0i(i＝1,2,3,…6)；△Ni＝PNi0.6          (9)式中，Ni为六种工况单独作用时的转子鸽尾部疲劳寿命次数；P为材料系数，根据《疲劳强度》手册，“1988年第一版，作者：徐灏”，P值取14；N0i＝Ni‑△Ni          (10)步骤15)：结合运行m年的各种工况实际出现次数，将步骤14)计算结果代入公式(11)，得到裂纹形成及扩展的评估系数DG‑M；DG-M=ΣiniN0i---(11)]]>式中，ni(i＝1,2,3,…6)分别为运行m年出现的发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况次数；步骤16)：假设在运行xm年后裂纹形成，则有公式(12)可解得x的值；xDG‑M＝1‑DG‑M                                        (12)步骤17)：根据公式(13)求解裂纹的扩展寿命系数y；yΣiniΔN0i=1---(13)]]>步骤18)：根据公式(14)计算求得累积疲劳寿命年限NGrover‑Manson，NGrover‑Manson＝m(1+x)+my                                  (14)。