CN106021714A - 一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法 - Google Patents

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

Info

Publication number
CN106021714A
CN106021714A CN201610331418.XA CN201610331418A CN106021714A CN 106021714 A CN106021714 A CN 106021714A CN 201610331418 A CN201610331418 A CN 201610331418A CN 106021714 A CN106021714 A CN 106021714A
Authority
CN
China
Prior art keywords
operating mode
stress
rotor
under
fatigue life
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610331418.XA
Other languages
English (en)
Other versions
CN106021714B (zh
Inventor
张宇娇
刘东圆
黄雄峰
徐彬昭
吴刚梁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN201610331418.XA priority Critical patent/CN106021714B/zh
Publication of CN106021714A publication Critical patent/CN106021714A/zh
Application granted granted Critical
Publication of CN106021714B publication Critical patent/CN106021714B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04—Ageing analysis or optimisation against ageing

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Geometry (AREA)
  • Game Theory and Decision Science (AREA)
  • Evolutionary Computation (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
  • Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)

Abstract

一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，根据抽水蓄能电厂提供实际不同运行工况下的转子速度曲线得到加速度分布作为载荷，得到发电电动机转子鸽尾部在发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下离心力所引起的应力随时间变化的分布情况；将由电磁损耗引起的电机温度变化，而由此产生的热应力与离心力引起的应力分布进行矢量求和，得到总应力分布；将已运行m年的发电电动机实际各种工况发生次数，及各种工况下疲劳寿命使用次数代入Miner准则计算公式，根据寿命预测系数及实际运行时间年数m计算得到该发电电动机转子鸽尾部的疲劳寿命使用年限。

Description

一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法
技术领域
本发明一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，涉及抽水蓄能发电电动机转子鸽尾部疲劳寿命预测领域。
背景技术
抽水蓄能发电电动机有发电和电动两种主要工况，其转速高，单机容量大，且机组启停频繁、转子存在正反转情况，运行工况复杂，运行环境较一般水轮发电机更恶劣，同时随着抽水蓄能电站运行年限的增长，使得发电电动机的故障和事故频发，安全稳定问题日益突出。
目前对于抽水蓄能发电电动机转子鸽尾部疲劳寿命预测的研究较少，由于发电电动机存在发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况，各种工况对转子鸽尾部的综合作用形成应力载荷谱作为累积疲劳分析的载荷，再结合一定的计算准则进行疲劳寿命预测。已有研究方法是通过静力学计算得到每种工况的最大应力分布，再将其作为载荷进行疲劳寿命预测，这种方法无法反映出应力分布随时间变化的情况。此外，也没有通过计算准则对抽水蓄能发电电动机转子鸽尾部进行疲劳寿命预测的研究。
发明内容
为解决上述技术问题，本发明提供一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，根据抽水蓄能电厂提供实际不同运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算得到不同工况下离心力所引起的应力随时间变化的分布情况，并结合热应力的作用，分别对各种工况单独进行疲劳寿命使用次数计算，根据已运行m年的发电电动机各种工况实际发生次数，应用Miner准则，实现了考虑各种工况共同影响的累积疲劳寿命预测。
本发明所采用的技术方案是：
一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，包括以 下步骤：
1)、根据抽水蓄能电厂提供实际运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算，分别得到发电电动机转子鸽尾部在发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下离心力所引起的应力随时间变化的分布情况；
2)、由电磁损耗引起电机的温度变化，而由此产生的热应力与离心力引起的应力进行矢量求和；
3)、根据发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下的总应力分布情况，绘制在不同工况时最容易出现疲劳破坏点处的应力随时间变化曲线，并根据曲线计算此处可使用的疲劳寿命次数；
4)、将已运行m年的发电电动机实际各种工况发生次数，及各种工况下疲劳寿命使用次数代入Miner准则计算公式，可得到寿命预测系数，并根据寿命预测系数及已知运行年数m计算得到该发电电动机转子鸽尾部的累积疲劳寿命使用年限。
一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，包括以下步骤：
步骤1)：建立电机1/2周期模型，采用有限元法进行电磁场-温度场-结构场耦合数值计算，通过对电磁场控制方程(1)-(3)和温度场控制方程(4)(5)进行有限元数值计算得到由电磁损耗引起的温度稳态分布，再由温度相比初始温度的变化对方程(6)进行求解得到热应力分布情况；
式中，V1是涡流区(转子绕组)，V2为源电流区(定子绕组)，σ为电导率，μ为相对磁导率，为源电流密度，Q为电磁损耗(包括源电流及涡流引起的损耗)。
式中，Q为能量损耗；kx,ky,kz分别表示热导率的各向异性参数；h为传热系数；T为求解温度；T0为环境温度。
式中，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)计算得到的热应力进行矢量求和，得到该工况下总应力；
式中，ρ为密度，ui为位移，其余参数与公式(6)相同。
步骤3)：根据抽水蓄能电厂提供实际发电停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到发电停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤4)：根据抽水蓄能电厂提供实际电动启动运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动启动工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤5)：根据抽水蓄能电厂提供实际电动停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤6)：根据抽水蓄能电厂提供实际甩负荷运行工况下的转子速度曲线得到加速度分 布作为载荷，通过动力学计算公式(7)，求解得到甩负荷工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤7)：根据抽水蓄能电厂提供实际飞逸运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到飞逸工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤8)：根据步骤2)计算得到的发电启动工况时转子3号鸽尾B点(如图2所示)处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N1；
式中：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)：根据某台发电电动机实际运行m年中各种工况发生次数，将步骤8)-13)计算出的结果代入Miner准则公式(9)，计算出寿命预测系数Dminer；
式中，ni(i＝1,2,3,…6)分别为运行m年出现的发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况次数；Ni(i＝1,2,3,…6)为步骤8)-13)中计算得到的每种工况下转子鸽尾部可使用疲劳寿命次数；
步骤15)：根据电机已知运行时间m年数及步骤14)计算得到的Dminer，代入公式(10)可求出转子鸽尾部疲劳寿命预测年限Nminer。
本发明一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，优点在于：
1)、根据抽水蓄能电厂提供实际运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算得到不同工况下离心力所引起的应力随时间变化的分布情况，将不同工况的应力变化以应力载荷谱的形式作为疲劳分析的激励，可以考虑长期运行在不同工况下的累积损伤对疲劳的综合作用。比常规方法，通过静力学计算出一个最大应力作为疲劳寿命计算载荷更符合实际情况，更准确。
2)、与电磁损耗引起的温度变化及由此产生热应力进行矢量求和，考虑总应力对疲劳的影响，也更符合实际工况。
3)、由于在电机设计阶段，无法预知每年实际发生的六种工况(即发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况)的具体次数，只能对其关键部件的疲劳寿命做简单预测，本方法将电机实际运行时间中的已知各种工况发生次数代入计算中，再对鸽尾部的疲劳寿命进行预测，得到的结果更合理更准确。
4)、Miner准则相较于其他疲劳寿命预测准则，只需计及实际运行时间中发生的各种工况次数与每种工况单独作用时的疲劳寿命次数，无需加入更多参数进行计算，既考虑了多种工况的累积损伤影响又结合了工程实际情况，同时计算方法更为简便，对于生产实际一线技术人员来说更易于掌握。
附图说明
图1是发电电动机1/2周期模型结构示意图。
图2是转子绕组鸽尾部局部放大示意图。
图3是电机转子稳态温度分布图。
图4是电机转子热应力分布图。
图5是发电启动工况时转子3号鸽尾B点处应力随时间变化曲线。
图6是发电停机工况时转子3号鸽尾B点处应力随时间变化曲线。
图7是电动启动工况时转子3号鸽尾B点处应力随时间变化曲线。
图8是电动停机工况时转子3号鸽尾B点处应力随时间变化曲线。
图9是甩负荷工况时转子3号鸽尾B点处应力随时间变化曲线。
图10是飞逸工况时转子3号鸽尾B点处应力随时间变化曲线。
具体实施方式
一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，包括以下步骤：
1)、根据抽水蓄能电厂提供实际运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算，分别得到发电电动机转子鸽尾部在发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下离心力所引起的应力随时间变化的分布情况；
2)、由电磁损耗引起电机的温度变化，而由此产生的热应力与离心力引起的应力进行矢量求和；
3)、根据发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下的总应力分布情况，绘制在不同工况时最容易出现疲劳破坏点处的应力随时间变化曲线，并根据曲线计算此处可使用的疲劳寿命次数；
4)、将已运行m年的发电电动机实际各种工况发生次数，及各种工况下疲劳寿命使用次数代入Miner准则计算公式，可得到寿命预测系数，并根据寿命预测系数及已知运行年数m计算得到该发电电动机转子鸽尾部的累积疲劳寿命使用年限。
一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，包括以下步骤：
步骤1)：建立电机1/2周期模型，采用有限元法进行电磁场-温度场-结构场耦合数值计算，通过对电磁场控制方程(1)-(3)和温度场控制方程(4)(5)进行有限元数值计算得到由电磁损耗引起的温度稳态分布，再由温度相比初始温度的变化对方程(6)进行求解得到热应力分布情况；
式中，V1是涡流区(转子绕组)，V2为源电流区(定子绕组)，σ为电导率，μ为相对磁导 率，为源电流密度，Q为电磁损耗(包括源电流及涡流引起的损耗)。
式中，Q为能量损耗；kx,ky,kz分别表示热导率的各向异性参数；h为传热系数；T为求解温度；T0为环境温度。
式中，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)计算得到的热应力进行矢量求和，得到该工况下总应力；
式中，ρ为密度，ui为位移，其余参数与公式(6)相同。
步骤3)：根据抽水蓄能电厂提供实际发电停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到发电停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤4)：根据抽水蓄能电厂提供实际电动启动运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动启动工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤5)：根据抽水蓄能电厂提供实际电动停机运行工况下的转子速度曲线得到加速度 分布作为载荷，通过动力学计算公式(7)，求解得到电动停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤6)：根据抽水蓄能电厂提供实际甩负荷运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到甩负荷工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤7)：根据抽水蓄能电厂提供实际飞逸运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到飞逸工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤8)：根据步骤2)计算得到的发电启动工况时转子3号鸽尾B点(如图2所示)处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N1；
式中：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)：根据某台发电电动机实际运行m年中各种工况发生次数，将步骤8)-13)计算出的结果代入Miner准则公式(9)，计算出寿命预测系数Dminer；
式中，ni(i＝1,2,3,…6)分别为运行m年出现的发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况次数；Ni(i＝1,2,3,…6)为步骤8)-13)中计算得到的每种工况下转子鸽尾部可使用疲劳寿命次数；
步骤15)：根据电机已知运行时间m年数及步骤14)计算得到的Dminer，代入公式(10)可求出转子鸽尾部疲劳寿命预测年限Nminer。
图3是电机转子稳态温度分布图。在步骤1)中进行了电磁场-温度场的耦合计算后得到的温度分布云图；
图4是电机转子热应力分布图。在步骤1)中得到温度分布(如图3)后，作为载荷加载到结构场进行计算，得到温度引起的形变情况，即热应力分布情况；
图5是发电启动工况时转子3号鸽尾B点处应力随时间变化曲线。步骤2)中计算得到电机在发电启动工况下的转子离心力引起的应力与热应力进行矢量求和，得到总应力分布，绘制出图2中3号鸽尾B点处应力随时间变化曲线；
图6是发电停机工况时转子3号鸽尾B点处应力随时间变化曲线。步骤3)中计算得到电机在发电停机工况下的转子离心力引起的应力与热应力进行矢量求和，得到总应力分布，绘制出图2中3号鸽尾B点处应力随时间变化曲线；
图7是电动启动工况时转子3号鸽尾B点处应力随时间变化曲线。步骤4)中计算得到电机在电动启动工况下的转子离心力引起的应力与热应力进行矢量求和，得到总应力分布，绘制出图2中3号鸽尾B点处应力随时间变化曲线；
图8是电动停机工况时转子3号鸽尾B点处应力随时间变化曲线。步骤5)中计算得到电机在电动停机工况下的转子离心力引起的应力与热应力进行矢量求和，得到总应力分布，绘制出图2中3号鸽尾B点处应力随时间变化曲线；
图9是甩负荷工况时转子3号鸽尾B点处应力随时间变化曲线。步骤6)中计算得到电机在甩负荷工况下的转子离心力引起的应力与热应力进行矢量求和，得到总应力分布，绘制出图2中3号鸽尾B点处应力随时间变化曲线；
图10是飞逸工况时转子3号鸽尾B点处应力随时间变化曲线。步骤7)中计算得到电机在飞逸工况下的转子离心力引起的应力与热应力进行矢量求和，得到总应力分布，绘制出图2中3号鸽尾B点处应力随时间变化曲线。
具体算例：以广蓄B厂8#机组为例
首先建立电机的1/2周期模型如图1、2所示，按照步骤1)对电机电磁场、温度场、结构场进行有限元数值计算，得到其转子的稳态温度分布如图3所示，再得到温度引起的热应力分布如图4所示。
根据步骤2)-步骤7)分别将转子鸽尾部在发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下离心力引起的应力与热应力进行矢量求和，得到各种工况的总应力分布，并绘制出鸽尾部总应力最大点处的应力分布曲线，如图5-10所示。再根据以上各工况应力分布曲线分别计算出B点在发电启动工况疲劳寿命使用次数N1为45320次，发电停机工况疲劳寿命使用次数为N2为149300次，电动启动工况疲劳寿命使用次数N3为217500次，电动停机工况疲劳寿命使用次数N4为362000次，甩负荷工况疲劳寿命使用次数N5为8991次，飞逸工况疲劳寿命使用次数N6为781次。
广蓄#8机组历年运行统计数据，自2000年3月14日#8机组投运起，至2012年2月22日机组开始大修为止，#8机组运行13年间，发电启停各9543次，电动启停各8286次，甩负荷26次，飞逸0次。根据步骤14)中的公式(9)计算出寿命预测系数Dminer：
根据步骤15)中的公式(10)计算得到转子鸽尾部疲劳寿命预测年限Nminer：

Claims (2)

1.一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，其特征在于包括以下步骤：
1)、根据抽水蓄能电厂提供实际运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算，分别得到发电电动机转子鸽尾部在发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下离心力所引起的应力随时间变化的分布情况；
2)、由电磁损耗引起电机的温度变化，而由此产生的热应力与离心力引起的应力进行矢量求和；
3)、根据发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况下的总应力分布情况，绘制在不同工况时最容易出现疲劳破坏点处的应力随时间变化曲线，并根据曲线计算此处可使用的疲劳寿命次数；
4)、将已运行m年的发电电动机实际各种工况发生次数，及各种工况下疲劳寿命使用次数代入Miner准则计算公式，可得到寿命预测系数，并根据寿命预测系数及已知运行年数m计算得到该发电电动机转子鸽尾部的累积疲劳寿命使用年限。
2.根据权利要求1所述一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法，其特征在于包括以下步骤：
步骤1)：建立电机1/2周期模型，采用有限元法进行电磁场-温度场-结构场耦合数值计算，通过对电磁场控制方程(1)-(3)和温度场控制方程(4)(5)进行有限元数值计算得到由电磁损耗引起的温度稳态分布，再由温度相比初始温度的变化对方程(6)进行求解得到热应力分布情况；
式中，V1是涡流区(转子绕组)，V2为源电流区(定子绕组)，σ为电导率，μ为相对磁导率，为源电流密度，Q为电磁损耗，包括源电流及涡流引起的损耗；
式中，Q为能量损耗；kx,ky,kz分别表示热导率的各向异性参数；h为传热系数；T为求解温度；T0为环境温度，
式中，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)计算得到的热应力进行矢量求和，得到该工况下总应力；
式中，ρ为密度，ui为位移，其余参数与公式(6)相同；
步骤3)：根据抽水蓄能电厂提供实际发电停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到发电停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤4)：根据抽水蓄能电厂提供实际电动启动运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动启动工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤5)：根据抽水蓄能电厂提供实际电动停机运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到电动停机工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤6)：根据抽水蓄能电厂提供实际甩负荷运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到甩负荷工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤7)：根据抽水蓄能电厂提供实际飞逸运行工况下的转子速度曲线得到加速度分布作为载荷，通过动力学计算公式(7)，求解得到飞逸工况下由转子离心力所引起的应力随时间分布，并与步骤1)计算得到的热应力进行矢量求和，得到该工况下总应力；
步骤8)：根据步骤2)计算得到的发电启动工况时转子3号鸽尾B点，处应力分布曲线，由公式(8)计算出在此种工况运行情况下可使用的疲劳寿命次数N1；
式中：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)：根据某台发电电动机实际运行m年中各种工况发生次数，将步骤8)-13)计算出的结果代入Miner准则公式(9)，计算出寿命预测系数Dminer；
式中，ni(i＝1,2,3,…6)分别为运行m年出现的发电启动工况、发电停机工况、电动启动工况、电动停机工况、甩负荷工况、飞逸工况次数；Ni(i＝1,2,3,…6)为步骤8)-13)中计算得到的每种工况下转子鸽尾部可使用疲劳寿命次数；
步骤15)：根据电机已知运行时间m年数及步骤14)计算得到的Dminer，代入公式(10)可求出转子鸽尾部疲劳寿命预测年限Nminer，
CN201610331418.XA 2016-05-18 2016-05-18 一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法 Active CN106021714B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610331418.XA CN106021714B (zh) 2016-05-18 2016-05-18 一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610331418.XA CN106021714B (zh) 2016-05-18 2016-05-18 一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法

Publications (2)

Publication Number Publication Date
CN106021714A true CN106021714A (zh) 2016-10-12
CN106021714B CN106021714B (zh) 2017-09-15

Family

ID=57097784

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610331418.XA Active CN106021714B (zh) 2016-05-18 2016-05-18 一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法

Country Status (1)

Country Link
CN (1) CN106021714B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106990002A (zh) * 2017-02-22 2017-07-28 国家电网公司 混流式水轮机承重支架疲劳寿命测试分析方法
CN109033488A (zh) * 2018-05-25 2018-12-18 华中科技大学 一种基于水机电耦合模型的抽水蓄能机组寿命分析方法
CN112594015A (zh) * 2020-11-23 2021-04-02 华能国际电力股份有限公司大连电厂 一种汽轮机寿命预测方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200057A (zh) * 2014-07-02 2014-12-10 三峡大学 抽水蓄能发电电动机转子t尾部剩余使用寿命的预测方法
WO2014204232A1 (ko) * 2013-06-20 2014-12-24 Lee Jung Chul 세라믹 타격절단장치
CN105205264A (zh) * 2015-09-23 2015-12-30 中国南方电网有限责任公司调峰调频发电公司 预测抽水蓄能发电电动机转子构件使用寿命的试验方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014204232A1 (ko) * 2013-06-20 2014-12-24 Lee Jung Chul 세라믹 타격절단장치
CN104200057A (zh) * 2014-07-02 2014-12-10 三峡大学 抽水蓄能发电电动机转子t尾部剩余使用寿命的预测方法
CN105205264A (zh) * 2015-09-23 2015-12-30 中国南方电网有限责任公司调峰调频发电公司 预测抽水蓄能发电电动机转子构件使用寿命的试验方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106990002A (zh) * 2017-02-22 2017-07-28 国家电网公司 混流式水轮机承重支架疲劳寿命测试分析方法
CN106990002B (zh) * 2017-02-22 2019-06-14 国家电网公司 混流式水轮机承重支架疲劳寿命测试分析方法
CN109033488A (zh) * 2018-05-25 2018-12-18 华中科技大学 一种基于水机电耦合模型的抽水蓄能机组寿命分析方法
CN109033488B (zh) * 2018-05-25 2020-05-19 华中科技大学 一种基于水机电耦合模型的抽水蓄能机组寿命分析方法
CN112594015A (zh) * 2020-11-23 2021-04-02 华能国际电力股份有限公司大连电厂 一种汽轮机寿命预测方法及系统

Also Published As

Publication number Publication date
CN106021714B (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
Presas et al. Fatigue life estimation of Francis turbines based on experimental strain measurements: Review of the actual data and future trends
CN106021713A (zh) 一种基于Corten-Dolan准则的抽水蓄能发电电动机转子鸽尾部疲劳寿命预测方法
CN106021714A (zh) 一种基于Miner准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法
CN114412589A (zh) 核电汽轮机的多目标多维度多工况的监控方法
US20170292401A1 (en) Method for determination of fatigue lifetime consumption of an engine component
CN104034480B (zh) 渐进式弯曲转子的现场动平衡方法
CN105893713B (zh) 一种基于Grover‑Manson准则的抽水蓄能发电电动机转子鸽尾部累积疲劳寿命预测方法
Connor et al. CFD modelling of an entire synchronous generator for improved thermal management
Liu et al. Coupled torsional vibration and fatigue damage of turbine generator due to grid disturbance
Dhopade et al. Aeromechanical modeling of rotating fan blades to investigate high-cycle and low-cycle fatigue interaction
CN204532445U (zh) 超临界机组汽轮机转子热应力在线监测系统
Zarifakis et al. Models for the transient stability of conventional power generating stations connected to low inertia systems
Gazizova et al. Analysis of short circuit transients with separate operation of iron and steel industry power plant with account of dynamic characteristics of industrial load
Blocher et al. Time-Linearized Forced Response Analysis of a Counter Rotating Fan: Part I—Theoretical Concept of a Fully Time-Linear Forced Response Analysis
Liao et al. Finite element analysis and Lightweight design of hydro generator lower bracket
CN106021716B (zh) 一种抽水蓄能发电电动机转子鸽尾部力学性能检测方法
Peng Tip running clearances effects on tip vortices induced axial compressor rotor flutter
Daroogheh et al. Engine life evaluation based on a probabilistic approach
Ma et al. Interval analysis method for rotordynamic with uncertain parameters
Tsypkaykin et al. Full rotortrain modelling approach for assessment of dynamic stresses on rotor and blades due to torsional disturbances
太兴宇 et al. A new representation model of rotating blade-casing rub-impact
Alendar et al. The problems of creating a propulsion system of a new generation supersonic passenger aircraft
Raza et al. Analysis of material coating for damping in beam structures
Curtin et al. Evaluating splitter blade designs to avoid impeller failure in a high-speed unshrouded centrifugal compressor
Srinivasa Rao et al. Creep Life Consumption Monitoring of a Turbine Rotor Blade

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant