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

权利要求书

1.一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，其特征在于包括以下步骤：

步骤1)：建立电机1/2周期模型，采用有限元法进行电磁场-温度场-结构场耦合数值计算，通过对电磁场控制方程(1)-(3)和温度场控制方程(4)(5)进行有限元数值计算得到由电磁损耗引起的温度稳态分布，再由温度相比初始温度的变化对方程(6)进行求解得到热应力分布情况；

▿×(1μ▿×A·)-▿(1μ▿·A·)=J·sinV2---(2)]]>

Q=∫ΩJ2σdΩ---(3)]]>

式中，为源电流密度，电磁场数值计算公式中为矢量磁位的通用表达参数，ω为频率的通用表达参数；V₁是转子绕组，V₂为定子绕组，σ为电导率，μ为相对磁导率，Q为电磁损耗，包括源电流及涡流引起的损耗，J为电流密度的有效值；

∂∂x(kx∂T∂x)+∂∂y(ky∂T∂y)+∂∂z(kz∂T∂z)=-Q0---(4)]]>

k∂T∂n+h(T-T0)=0---(5)]]>

式中，Q₀为能量损耗；k_x,k_y,k_z分别表示热导率的各向异性参数；h为传热系数；T为求解温度；T₀为环境温度，

σ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为温度相比初始温度的变化量；F_i为外力的分量；u_i,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)]]>

式中，ρ为密度，u_i为位移，其余参数与公式(6)相同；

步骤3)：根据抽水蓄能电厂提供实际发电停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到发电停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；

步骤4)：根据抽水蓄能电厂提供实际电动启动运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动启动工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；

步骤5)：根据抽水蓄能电厂提供实际电动停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；

步骤6)：根据抽水蓄能电厂提供实际甩负荷运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到甩负荷工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；

步骤7)：根据抽水蓄能电厂提供实际飞逸运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到飞逸工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；

步骤8)：根据步骤2)计算得到的发电启动工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N₁；

N=Cσ-aσ=Kσϵσβσσmax-σmin2+ψaσmax+σmin2---(8)]]>

式中：C和a为材料疲劳系数；σ_max为总应力变化曲线中的最大值；σ_min为总应力变化曲线中的最小值；K_σ，ε_σ，β_σ和ψ_a分别为有效应力集中系数、零件尺寸系数、表面系数和平均应力系数；

步骤9)：根据步骤3)计算得到的发电停机工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N₂；

步骤10)：根据步骤4)计算得到的电动启动工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N₃；

步骤11)：根据步骤5)计算得到的电动停机工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N₄；

步骤12)：根据步骤6)计算得到的甩负荷工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N₅；

步骤13)：根据步骤7)计算得到的飞逸工况时转子3号鸽尾B点处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N₆；

步骤14)：根据某台发电电动机实际运行m年中各种工况发生次数，将步骤8)-13)计算出的结果代入Miner准则公式(9)，计算出寿命预测系数D_miner；

Dminer=ΣiniNi---(9)]]>

式中，n_i，i＝1,2,3,…6，分别为运行m年出现的发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况次数；N_i，i＝1,2,3,…6为步骤8)-13)中计算得到的每种工况下转子鸽尾部可使用疲劳寿命次数；

步骤15)：根据电机已知运行时间m年数及步骤14)计算得到的D_miner，代入公式(10)可求出转子鸽尾部疲劳寿命预测年限N_miner，

Nminer=mDminer---(10).]]> 3