CN108228979B - 一种用于发电机结构的低周疲劳寿命分析方法 - Google Patents
一种用于发电机结构的低周疲劳寿命分析方法 Download PDFInfo
- Publication number
- CN108228979B CN108228979B CN201711360071.2A CN201711360071A CN108228979B CN 108228979 B CN108228979 B CN 108228979B CN 201711360071 A CN201711360071 A CN 201711360071A CN 108228979 B CN108228979 B CN 108228979B
- Authority
- CN
- China
- Prior art keywords
- strain
- stress
- generator structure
- node
- epsilon
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种用于发电机结构的低周疲劳寿命分析方法,其特征在于,包括:步骤1:将材料的循环应力‑应变曲线输入发电机结构的有限元数值仿真模型,导出发电机结构所有节点在额定工况下的应变分量、在静止工况下的应变分量、额定工况下的等效应力、静止工况下的等效应力、第一主应力和第三主应力;步骤2:计算应变幅和平均应力;步骤3:计算发电机结构在任一节点处的低周疲劳寿命N。本发明提供了一种适用于多种发电机结构、计算精准、操作简单、可实现计算方法与流程标准化的发电机结构低周疲劳寿命分析方法。
Description
技术领域
本发明涉及一种用于发电机结构的低周疲劳寿命分析方法,适用于发电机与大型电机,属于发电机的结构分析技术领域。
背景技术
电网对电厂发电能力的动态调整,发电机设备定期或非定期维护,都需要对设备进行起停,这使得护环、转子端部齿、磁极引线和引水拐脚等很多部件将承受由频繁起停机引起的低周疲劳载荷。如果这些部件在设计安全寿命内发生疲劳断裂破坏将会给安全生产带来不可估量的损失。虽然疲劳损伤是一个连续的过程,但对损伤的在线监测或维修检测却都非常困难或不可能。此外,疲劳的发生经常具有突然性。因此,在发电机结构设计时,必须对发电机转子各部件进行低周疲劳寿命分析,以保证发电机的正常稳定运行。
以往的发电机结构低周疲劳寿命计算采用的都是疲劳理论的经典手工公式算法。即,采用解析公式或有限元数值计算得到结构在各工况下的应力、应变,再采用局部应力-应变法对结构危险部位的局部应力-应变响应进行弹塑性修正,最后利用应变-寿命曲线对结构的低周疲劳寿命进行评估。此种计算方法存在以下不足:(1)局部应力-应变法是以单轴应力为假设条件的,无法有效处理多轴应力问题;(2)结构的疲劳破坏危险点需要人为预判选择,容易导致疲劳分析的偏差;(3)公式计算结构寿命未实现流程化分析,效率较低,且结果不够显性化。
发明内容
本发明要解决的技术问题是提供一种适用于多种发电机结构、操作简单、可实现计算方法与流程标准化的发电机结构低周疲劳寿命分析方法。
为了解决上述技术问题,本发明提供了一种用于发电机结构的低周疲劳寿命分析方法,其特征在于,包括:
步骤1:根据发电机的实际结构建立有限元数值仿真模型,将材料的循环应力-应变曲线输入发电机结构的有限元数值仿真模型,导出发电机结构所有节点在额定工况下的6个应变分量(包括所有节点在额定工况下x向主应变分量εx、y向主应变分量εy、z向主应变分量εz、xy向剪应变分量εxy、yz向剪应变分量εyz和zx向剪应变分量εzx)、所有节点在静止工况下的6个应变分量(包括所有节点在静止工况下x向主应变分量εx'、y向主应变分量εy'、z向主应变分量εz'、xy向剪应变分量εxy'、yz向剪应变分量εyz'和zx向剪应变分量εzx')、额定工况下的等效应力(σeqv)额定、静止工况下的等效应力(σeqv)静止、第一主应力σ1和第三主应力σ2;
步骤2:计算应变幅和平均应力:
发电机结构在任一节点处的等效应变幅Δε的计算公式为:
式中,Δε为等效应变幅;dεx为x向的主应变幅,dεx=εx-εx';dεy为y向的主应变幅,dεy=εy-εy';dεz为z向的主应变幅,dεz=εz-εz';dεxy、dεyz和dεzx分别为xy、yx、zx三个方向的剪应变幅,dεxy=εxy-εxy',dεyz=εyz-εyz',dεzx=εzx-εzx';
发电机结构在循环交变工况下的平均应力σm的计算公式为:
式中,μ为修正系数,当第一主应力σ1的绝对值大于第三主应力σ2的绝对值,则节点处于受拉状态,μ取值为1;反之,则节点处于受压状态,μ取值为-1;
步骤3:计算发电机结构在任一节点处的低周疲劳寿命N,计算公式(应变-寿命公式)为:
式中,E为材料弹性模量;σf′为疲劳强度系数;b疲劳强度指数;εf′疲劳塑性系数;c疲劳塑性指数;N为发电机结构在任一节点处的低周疲劳寿命。
优选地,所述的用于发电机结构的低周疲劳寿命分析方法还包括:
步骤4:将发电机结构在每一节点处的低周疲劳寿命N的对数导入到步骤1中建立的有限元数值仿真模型中,并作为温度载荷施加到对应的节点上,显示结构的温度载荷分布即可显性化表示出结构的疲劳寿命分布。
所述的发电机结构为护环、转子端部齿、磁极引线、引水拐脚、平板、梁等。
为保证计算的精度,所述的有限元数值仿真模型采用2维的8节点等参单元或3维的20节点等参单元。
所述的步骤1中,由于低周疲劳分析中,结构承受的应力水平较高并已经入弹塑性范围。建模时需要输入材料的循环应力-应变曲线。采用弹塑性有限元分析计算出结构在额定和静止工况下的应力和应变。
所述的步骤2中,结构的低周疲劳寿命主要由应变决定,需要导出结构应变进行后续的疲劳分析。由于发电机结构通常处于多轴应力、应变状态,因此导出结构在静止和额定工况下的所有应变分量以组成应变张量。在三维笛卡尔坐标系下,应变张量可表示为
应变张量是一个对称张量,需要有限元分析导出的是其中6个独立的分量,包括:x向主应变εx,y向主应变εy,z向主应变εz以及三个剪应变分量εxy、εyz和εzx。
此外,结构的寿命还受到平均应力的影响,因此需要导出结构的等效应力的结果。等效应力为正值,无法考虑节点是处于受压还是受拉状态。因此,还要导出结构的第一和第三主应力来判断节点处的等效应力状态。
所述的步骤3中,结构在任一节点处由额定到静止工况下的应变变化范围可由下式计算得到
[dε]=[ε]额定-[ε]静止 (1)
式中,[ε]额定和[ε]静止分别为步骤2中从有限元分析导出的额定和静止工况下应变张量(3×3矩阵)。
由于通常情况下应变-寿命曲线是在对称循环条件下试验得出的。对于复杂载荷历程作用下的疲劳问题,平均应力的存在是不可避免的。因此,疲劳寿命分析时应计入平均应力的影响。结构在循环交变工况下的平均应力。可由有限元计算出的额定和静止工况下的等效应力取平均值得到。由于等效应力均为正值,无法考虑节点是处于受压还是受拉状态。因此,还需根据第一和第三主应力来判断节点处的等效应力状态。当第一主应力的绝对值大于第三主应力的绝对值,则节点处于受拉状态,平均应力为正值;反之节点处于受压状态,平均应力为负值。
与现有技术相比,本发明的有益效果是:
本发明提供了一种适用于多种发电机结构、计算精准、操作简单、可实现计算方法与流程标准化的发电机结构低周疲劳寿命分析方法。
本发明提供的方法克服了现有技术的不足,可以实现发电机结构低周疲劳寿命的精确、快速分析。同时将疲劳寿命结果做为温度载荷加载至有限元数值仿真模型上可以实现结构寿命的显性化显示。
本发明的优点在于可以完成多种发电机结构的低周疲劳寿命计算,操作简单,实现了计算方法与流程的标准化,大大提升了工作效率,使得原本1天的计算工作量大幅缩短为5分钟,一定程度上降低了设计开发成本。此外,通过结构疲劳寿命的可视化可以帮助研发人员快速准确判断结构的危险位置,避免了人工公式计算导致的误判断。
本发明通过有限元分析得到结构的应力、应变,对应力应变结果进行低周疲劳寿命分析,最后将结构的寿命结果导入到有限元数值仿真模型中,可以实现结构寿命的可视化。本发明的优点是将结构疲劳寿命分析流程化、可视化。相比之前的传统理论公式分析方法,大大提高了数据处理的效率,同时能将结构各位置的寿命直观地显示在有限元数值仿真模型上。
附图说明
图1为一悬臂梁结构的有限元数值仿真模型;
图2为悬臂梁疲劳寿命结果分布;1—寿命分布云图;2—寿命对数;3—最小寿命说明;
图3为低周疲劳分析流程。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
实施例
如图1所示,为一悬臂梁结构的有限元数值仿真模型,具体尺寸参数为长120mm,截面宽20mm、高20mm,圆角半径10mm,采用铜材料,施加载荷为自由端的±5kN的交变力。将其作为本实施例的分析对象。
如图3所示,首先由有限元弹塑性分析得到此悬臂梁分别5kN和-5kN载荷下应力和应变(步骤1)。再将应力和应变结果导入到自编的MATLAB程序(步骤2-3)进行低周疲劳分析,计算出悬臂梁有限元模型在各节点处的疲劳寿命。最后再将各节点处的疲劳寿命作为温度载荷加载于每个节点上并显示出来(步骤4)。
所述的用于发电机结构的低周疲劳寿命分析方法为:
步骤1:根据发电机的实际结构建立有限元数值仿真模型,如图1所示,建立图1所示的有限元数值仿真模型,为保证计算的精度建模采用3维的20节点等参单元。将材料的循环应力-应变曲线输入有限元数值仿真模型,导出发电机结构所有节点在额定工况下的6个应变分量(包括所有节点在额定工况下x向主应变分量εx、y向主应变分量εy、z向主应变分量εz、xy向剪应变分量εxy、yz向剪应变分量εyz和zx向剪应变分量εzx)、所有节点在静止工况下的6个应变分量(包括所有节点在静止工况下x向主应变分量εx'、y向主应变分量εy'、z向主应变分量εz'、xy向剪应变分量εxy'、yz向剪应变分量εyz'和zx向剪应变分量εzx')、额定工况下的等效应力(σeqv)额定、静止工况下的等效应力(σeqv)静止、第一主应力σ1,和第三主应力σ2。
步骤2:计算应变幅和平均应力:
发电机结构在任一节点处的等效应变幅Δε的计算公式为:
式中,Δε为等效应变幅;dεx为x向的主应变幅,dεx=εx-εx';dεy为y向的主应变幅,dεy=εy-εy';dεz为z向的主应变幅,dεz=εz-εz';dεxy、dεyz和dεzx分别为xy、yx、zx三个方向的剪应变幅,dεxy=εxy-εxy',dεyz=εyz-εyz',dεzx=εzx-εzx';
发电机结构在循环交变工况下的平均应力σm的计算公式为:
式中,μ为修正系数,当第一主应力σ1的绝对值大于第三主应力σ2的绝对值,则节点处于受拉状态,μ取值为1;反之,则节点处于受压状态,μ取值为-1。
步骤3:计算发电机结构在任一节点处的低周疲劳寿命N,计算公式为:
式中,E为材料弹性模量,取值为1.65×105MPa;σf′为疲劳强度系数,取值为506.605MPa;b疲劳强度指数,取值为-0.0508;εf′疲劳塑性系数,取值为0.0939;c疲劳塑性指数,取值为-0.3795;N为发电机结构在任一节点处的低周疲劳寿命。
步骤4:将发电机结构在每一节点处的低周疲劳寿命N的对数(底为10)导入到步骤1中建立的有限元数值仿真模型中,并作为温度载荷施加到对应的节点上,显示结构的温度载荷分布即可显性化表示出结构的疲劳寿命分布,如图2所示。
计算结果直观地显示出悬臂梁的较小寿命主要集中于固定端得圆角区域。最小寿命为1349次载荷交变,且位于3546号节点位置。
Claims (3)
1.一种用于发电机结构的低周疲劳寿命分析方法,其特征在于,包括:
步骤1:根据发电机的实际结构建立有限元数值仿真模型,将材料的循环应力-应变曲线输入发电机结构的有限元数值仿真模型,导出发电机结构所有节点在额定工况下的6个应变分量,包括所有节点在额定工况下x向主应变分量εx、y向主应变分量εy、z向主应变分量εz、xy向剪应变分量εxy、yz向剪应变分量εyz和zx向剪应变分量εzx,所有节点在静止工况下的6个应变分量,包括所有节点在静止工况下x向主应变分量εx'、y向主应变分量εy'、z向主应变分量εz'、xy向剪应变分量εxy'、yz向剪应变分量εyz'和zx向剪应变分量εzx',额定工况下的等效应力(σeqv)额定,静止工况下的等效应力(σeqv)静止,第一主应力σ1和第三主应力σ2;
步骤2:计算应变幅和平均应力:
发电机结构在任一节点处的等效应变幅Δε的计算公式为:
式中,Δε为等效应变幅;dεx为x向的主应变幅,dεx=εx-εx';dεy为y向的主应变幅,dεy=εy-εy';dεz为z向的主应变幅,dεz=εz-εz';dεxy、dεyz和dεzx分别为xy、yx、zx三个方向的剪应变幅,dεxy=εxy-εxy',dεyz=εyz-εyz',dεzx=εzx-εzx';
发电机结构在循环交变工况下的平均应力σm的计算公式为:
式中,μ为修正系数,当第一主应力σ1的绝对值大于第三主应力σ2的绝对值,则节点处于受拉状态,μ取值为1;反之,则节点处于受压状态,μ取值为-1;
步骤3:计算发电机结构在任一节点处的低周疲劳寿命N,计算公式为:
式中,E为材料弹性模量;σf′为疲劳强度系数;b疲劳强度指数;εf′疲劳塑性系数;c疲劳塑性指数;N为发电机结构在任一节点处的低周疲劳寿命。
2.如权利要求1所述的用于发电机结构的低周疲劳寿命分析方法,其特征在于,还包括:
步骤4:将发电机结构在每一节点处的低周疲劳寿命N的对数导入到步骤1中建立的有限元数值仿真模型中,并作为温度载荷施加到对应的节点上,显示结构的温度载荷分布即可显性化表示出结构的疲劳寿命分布。
3.如权利要求1所述的用于发电机结构的低周疲劳寿命分析方法,其特征在于,还包括:所述的有限元数值仿真模型采用2维的8节点等参单元或3维的20节点等参单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711360071.2A CN108228979B (zh) | 2017-12-15 | 2017-12-15 | 一种用于发电机结构的低周疲劳寿命分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711360071.2A CN108228979B (zh) | 2017-12-15 | 2017-12-15 | 一种用于发电机结构的低周疲劳寿命分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108228979A CN108228979A (zh) | 2018-06-29 |
CN108228979B true CN108228979B (zh) | 2021-04-13 |
Family
ID=62652345
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711360071.2A Active CN108228979B (zh) | 2017-12-15 | 2017-12-15 | 一种用于发电机结构的低周疲劳寿命分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108228979B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111695214B (zh) * | 2020-05-26 | 2023-09-08 | 上海拜安传感技术有限公司 | 基于统计模型的岸桥起重机疲劳损伤的确定方法 |
CN112115540B (zh) * | 2020-09-27 | 2022-12-20 | 武汉第二船舶设计研究所(中国船舶重工集团公司第七一九研究所) | 一种融合实测数据的海工支撑结构疲劳分析方法及系统 |
CN112378652A (zh) * | 2020-11-10 | 2021-02-19 | 北京航空航天大学 | 一种用于航空膜盘联轴器寿命分析方法 |
CN113392558B (zh) * | 2021-06-24 | 2024-02-23 | 江南大学 | 一种基于最大剪应变与主拉应变的岩石损伤确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101750216A (zh) * | 2010-01-28 | 2010-06-23 | 清华大学 | 次同步振荡引起的汽轮发电机轴系疲劳损伤在线分析方法 |
CN106484954A (zh) * | 2016-09-14 | 2017-03-08 | 华北电力大学 | 一种燃气蒸汽联合循环机组启动曲线在线优化调整方法 |
CN107389471A (zh) * | 2017-07-13 | 2017-11-24 | 电子科技大学 | 裂纹失效模式判定方法及基于该方法的疲劳寿命预测方法 |
CN107423540A (zh) * | 2017-04-09 | 2017-12-01 | 北京工业大学 | 一种基于权平均最大剪切应变幅平面的多轴疲劳寿命预测方法 |
-
2017
- 2017-12-15 CN CN201711360071.2A patent/CN108228979B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101750216A (zh) * | 2010-01-28 | 2010-06-23 | 清华大学 | 次同步振荡引起的汽轮发电机轴系疲劳损伤在线分析方法 |
CN106484954A (zh) * | 2016-09-14 | 2017-03-08 | 华北电力大学 | 一种燃气蒸汽联合循环机组启动曲线在线优化调整方法 |
CN107423540A (zh) * | 2017-04-09 | 2017-12-01 | 北京工业大学 | 一种基于权平均最大剪切应变幅平面的多轴疲劳寿命预测方法 |
CN107389471A (zh) * | 2017-07-13 | 2017-11-24 | 电子科技大学 | 裂纹失效模式判定方法及基于该方法的疲劳寿命预测方法 |
Non-Patent Citations (2)
Title |
---|
QFS型180MW级双水内冷发电机试验分析;吴新亚 等;《电机技术》;20171031;第36-39页 * |
大型汽轮发电机转子低周疲劳寿命分析;杨昔科 等;《电机技术》;20150630;第5-7页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108228979A (zh) | 2018-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108228979B (zh) | 一种用于发电机结构的低周疲劳寿命分析方法 | |
CN112710539B (zh) | 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法 | |
CN108204925B (zh) | 复合材料的疲劳寿命预测方法及预测系统 | |
Ye et al. | A multi-scale model for studying failure mechanisms of composite wind turbine blades | |
CN109614715B (zh) | 一种多轴载荷作用下考虑缺口效应的能量场强法及其应用 | |
CN106202630A (zh) | 汽轮机转子超高周疲劳疲劳强度与疲劳寿命评估方法 | |
CN102589891A (zh) | 一种车用电机驱动系统耐久寿命的估计方法 | |
CN105008887A (zh) | 使用非接触测量和动态响应重构技术的涡轮机叶片疲劳寿命分析 | |
CN109828033B (zh) | 基于振动响应相似度分析的损伤识别方法和系统 | |
CN103217282B (zh) | 一种基于等寿命原则的水平轴风力机叶片缩比模型疲劳测试方法 | |
WO2009071817A3 (fr) | Procede de determination de la valeur d'un parametre representatif de la manoeuvrabilite d'un reacteur nucleaire | |
CN109855991B (zh) | 一种预测复合材料疲劳寿命的方法 | |
CN109858152B (zh) | 一种预测含初始分层损伤复合材料的剩余疲劳寿命的方法 | |
CN109614640A (zh) | 一种大型风电机组轮毂疲劳寿命预测方法及系统 | |
Xie et al. | Comprehensive fatigue estimation and fault diagnosis based on Refined Generalized Multi-Scale Entropy method of centrifugal fan blades | |
CN106596255A (zh) | 一种基于光纤光谱图像分析的孔边裂纹诊断方法 | |
Liu et al. | Crack growth analysis and fatigue life estimation in the piston rod of a Kaplan hydro turbine | |
CN117807659A (zh) | 一种建筑结构缺陷检测和纠正方法及系统 | |
CN103991556B (zh) | 一种载荷测量方法 | |
CN109948216B (zh) | 总应变能密度修正的缺口件低周疲劳预测方法 | |
CN114036792B (zh) | 一种类裂纹不连续区结构的疲劳损伤系数获取方法及装置 | |
Singh et al. | A multiscale approach to modeling of composite damage | |
CN112906237A (zh) | 一种发动机部件故障分析方法及系统 | |
Vettori et al. | A virtual sensing approach to operational modal analysis of wind turbine blades | |
CN102879185A (zh) | 阀门信号测试方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |