CN114136573B - 水轮发电机组关键部件的预警振幅的计算方法及相关设备 - Google Patents
水轮发电机组关键部件的预警振幅的计算方法及相关设备 Download PDFInfo
- Publication number
- CN114136573B CN114136573B CN202111437224.5A CN202111437224A CN114136573B CN 114136573 B CN114136573 B CN 114136573B CN 202111437224 A CN202111437224 A CN 202111437224A CN 114136573 B CN114136573 B CN 114136573B
- Authority
- CN
- China
- Prior art keywords
- working condition
- vibration
- amplitude
- fatigue life
- finite element
- 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
- 238000000034 method Methods 0.000 title claims description 45
- 238000004088 simulation Methods 0.000 claims abstract description 60
- 238000001228 spectrum Methods 0.000 claims abstract description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 37
- 230000005284 excitation Effects 0.000 claims abstract description 34
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 230000004044 response Effects 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims description 38
- 230000003068 static effect Effects 0.000 claims description 10
- 230000005484 gravity Effects 0.000 claims description 9
- 238000005215 recombination Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 13
- 238000004891 communication Methods 0.000 description 11
- 238000006073 displacement reaction Methods 0.000 description 11
- 230000003595 spectral effect Effects 0.000 description 11
- 239000000463 material Substances 0.000 description 10
- 230000006870 function Effects 0.000 description 7
- 230000003938 response to stress Effects 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000013016 damping Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/06—Multidirectional test stands
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/025—Measuring arrangements
-
- 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
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本申请提供一种水轮发电机组关键部件的预警振幅的计算方法，包括：建立水轮发电机组关键部件的有限元模型；基于有限元模型，分析关键部件在不同振动频率激励下的振动应力，得到关键部件的频率响应特性；编制典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱；基于频率响应特性，将振动幅值功率谱作为输入激励加载到有限元模型中，疲劳仿真分析得到关键部件的疲劳寿命；设定疲劳寿命阈值，疲劳寿命阈值小于疲劳寿命；将以预设倍数调整振动幅值后的振动幅值功率，作为输入激励加载到有限元模型中，仿真求解得到关键部件的输出疲劳寿命；当输出疲劳寿命接近疲劳寿命阈值时，设定其对应的振动幅值为预警幅值。
Description
技术领域
本申请涉及水轮发电机组技术领域，尤其涉及一种水轮发电机组关键部件的预警振幅的计算方法及相关设备。
背景技术
针对水轮发电机组的关键部件的安全预警，主要通过对在线监测各测点设置预警振幅的方法。因此，预警振幅的设置对于水轮发电机组的安全运行具有重要的意义。
传统的预警振幅的设定方法，在选择安全系数(预警振幅)时多选择较大的安全系数(预警振幅)，或者给设备限定过于苛刻的振动条件。在多数场合这种评估方法偏于保守，且针对不同运行特性的机组，采用单一的振动幅值标准，没有充分考虑机组的运行特性对振动幅值造成的影响，不能适用于运行特性不同的多种机组，因此，不能所有水轮发电机组的运行需求，会降低电厂的发电效益。
发明内容
有鉴于此，本申请的目的在于提出一种水轮发电机组关键部件的预警振幅的计算方法及相关设备。
基于上述目的，本申请提供了一种水轮发电机组关键部件的预警振幅的计算方法，包括：
建立水轮发电机组关键部件的有限元模型；
基于所述有限元模型，分析所述关键部件在不同振动频率激励下的振动应力，得到所述关键部件的频率响应特性；
编制典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱；
基于所述频率响应特性，将所述振动幅值功率谱作为输入激励加载到所述有限元模型中，疲劳仿真分析得到所述关键部件的疲劳寿命；
设定疲劳寿命阈值，所述疲劳寿命阈值小于所述疲劳寿命；
将以预设倍数调整振动幅值后的振动幅值功率，作为输入激励加载到所述有限元模型中，仿真求解得到所述关键部件的输出疲劳寿命；当所述输出疲劳寿命接近所述疲劳寿命阈值时，设定其对应的振动幅值为预警幅值。
在一些实施例中，所述疲劳寿命阈值根据所述疲劳寿命和服役全寿命设定；当所述疲劳寿命呈现无限寿命时，所述输出疲劳寿命与所述疲劳寿命阈值的比值x的数值满足，0.1<x<1。
在一些实施例中，所述编制典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱具体包括：
统计各典型工况工作时间占机组服役寿命的比重；
根据各典型工况的实际比重对测试得到各典型工况下的载荷谱进行重新组合；
将重新组合的结果进行时频变换，编制得到所述典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱。
在一些实施例中，所述根据各典型工况的实际比重对测试得到各典型工况下的载荷谱进行重新组合具体包括：
根据式进行重新组合；其中，n1为静止工况的运行时间；n2为额定工况的运行时间；n8为叠加工况的运行时间；n为各种工况的运行时间之和。
在一些实施例中，所述建立水轮发电机组关键部件的有限元模型具体包括：
通过建立实体模型、网格划分和单元质量控制建立水轮发电机组关键部件的初始有限元模型；
仿真计算：计算所述关键部件的初始有限元模型的仿真频率和仿真振形，得到仿真结果；
现场试验：获取所述关键部件的现场频率和现场振形，得到现场试验结果；
比对结果：比对仿真结果和现场试验结果，当仿真结果与现场结果不一致时，修正所述初始有限元模型，并重复所述仿真计算和现场试验的步骤，直至仿真结果与现场试验结果接近。
在一些实施例中，所述比对结果的步骤之后还包括：
计算所述比对结果所得的有限元模型在典型工况下的振动幅值和振动应力分布，得到典型工况仿真结果；
统计典型工况仿真结果应力值较大的部位，识别易发生失效的潜在危险点；
获取所述关键部件的中心部位在典型工况下的振动幅值，得到典型工况现场试验结果；
统计分析所述典型工况现场试验结果，得到所述关键部件的极端工况；
计算所述比对结果所得的有限元模型的潜在危险点在极端工况下的振动幅值和振动应力分布，得到极端工况仿真结果；
获取所述关键部件的潜在危险点在极端工况下的振动幅值和振动应力分布，得到极端工况现场试验结果；
比对所述极端工况仿真结果和所述极端工况现场试验结果，当不一致时，修正所述比对结果所得的有限元模型，并重复所述仿真计算和现场试验的步骤，直至仿真结果与现场结果接近。
在一些实施例中，所述典型工况包括静止工况、额定工况和叠加工况；所述关键部件包括混流式水轮机组机架、混流式水轮机组顶盖、贯流式水轮机组水导轴承或贯流式水轮机组转轮室。
在一些实施例中，所述关键部件为混流式水轮机组顶盖时，所述典型工况为静止工况、额定工况和叠加工况。
本申请实施例还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如前任意一项所述的方法。
本申请实施例还提供一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令用于使计算机执行如前任一所述方法。
从上面所述可以看出，本申请提供的水轮发电机组关键部件的预警振幅的计算方法，通过建立水轮发电机组关键部件的有限元模型；基于所述有限元模型，分析所述关键部件在不同振动频率激励下的振动应力，得到所述关键部件的频率响应特性；编制典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱；基于所述频率响应特性，将所述振动幅值功率谱作为输入激励加载到所述有限元模型中，疲劳仿真分析得到所述关键部件的疲劳寿命；设定疲劳寿命阈值，所述疲劳寿命阈值小于所述疲劳寿命；将以预设倍数调整振动幅值后的振动幅值功率，作为输入激励加载到所述有限元模型中，仿真求解得到所述关键部件的输出疲劳寿命；当所述输出疲劳寿命接近所述疲劳寿命阈值时，设定其对应的振动幅值为预警幅值。能够基于频响FRF特性，以典型工况下，服役寿命周期内水轮发电机组关键部件的振动幅值功率谱作为输入，结合材料S-N曲线进行仿真，求解得到结构的疲劳寿命。以输出疲劳寿命为目标，不断调整振动幅值利用仿真出来的疲劳寿命是否接近疲劳寿命阈值，从而利用调整出来的振动幅值作为关键部件的振动幅值作为报警值来设置。
附图说明
为了更清楚地说明本申请或相关技术中的技术方案，下面将对实施例或相关技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例的水轮发电机组关键部件的预警振幅的计算方法的流程示意图；
图2为本申请实施例的混流式水轮机组顶盖的监测点的应力随X方向单位振幅激励的变化示意图；
图3为本申请实施例的混流式水轮机组顶盖的监测点的应力随Y方向单位振幅激励的变化示意图；
图4为本申请实施例的混流式水轮机组顶盖的监测点的应力随Z方向单位振幅激励的变化示意图；
图5为本申请实施例的编制典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱的流程示意图；
图6为本申请实施例的各个典型工况的运行时间占总运行时间的比重的示意图；
图7为本申请实施例的重新组合得到的服役周期内关键部件顶盖的振动幅值时域图；
图8a为本申请实施例的机组服役寿命周期内关键部件顶盖在X方向的振动幅值功率谱的密度示意图；
图8b为本申请实施例的机组服役寿命周期内关键部件顶盖在Y方向的振动幅值功率谱的密度示意图；
图8c为本申请实施例的机组服役寿命周期内关键部件顶盖在Z方向的振动幅值功率谱的密度示意图；
图9为本申请实施例的振动疲劳分析的流程示意图；
图10为本申请实施例的振动疲劳分析的具体流程示意图；
图11为本申请实施例的关键部件顶盖的主要材料为Q235时的S-N曲线示意图；
图12为本申请实施例的电子设备的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本申请进一步详细说明。
需要说明的是，除非另外定义，本申请实施例使用的技术术语或者科学术语应当为本申请所属领域内具有一般技能的人士所理解的通常意义。本申请实施例中使用的“包括”或者“包含”等类似的词语意指出现该词前面的元件或者物件涵盖出现在该词后面列举的元件或者物件及其等同，而不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电性的连接，不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也可能相应地改变。
一些水轮机结构关键部件的预警振幅的评估，针对关键部件模型的设计方法主要依赖理论公式以及设计经验，采用静强度对结构安全性进校核，并用动载系统考虑振动对结构安全性的影响。由于传统经验和理论公式难以定量计算额定工作下的应力水平，无法计算结构的动力学特性和振动响应导致的附加应力，更无法计算振动对结构的损伤；且实际运行过程中，由于振动激励源多且特性复杂，系统结构动力学特性难以掌握及控制，因此现有的预警振幅在选择安全系数时多选择较大的安全系数，或者给设备限定过于苛刻的振动条件。
而在多数场合现有的预警振幅的评估方法偏于保守，且根据设备的实际运行情况，以及安全评估后，一些水轮发电机组设备在超报警阈值情况下是可以继续稳定运行的。但是由于目前的预警振幅值，水电厂集控调度会对出现超报警阈值的机组进行负荷限制，对电厂的发电效益造成影响。且现有的预警振幅的评估方法由于未能评定各部件的振动响应及影响，容易对真正的风险部位漏评估。
基于此，本申请实施例提供了基于有限元模型计算水轮机关键部件的额定静工况应力、转动激励下的动力学响应及进行多工况疲劳寿命预测，结合设备实际运行的安全分析，得到预警振幅，能够在一定程度上解决预警振幅的设定偏保守的问题。
请参阅图1，本申请实施例提供一种水轮发电机组关键部件的预警振幅的计算方法，包括：
S100,建立水轮发电机组关键部件的有限元模型；
S200,基于所述有限元模型，分析所述关键部件在不同振动频率激励下的振动应力，得到所述关键部件的频率响应特性；
S300,编制典型工况下，服役寿命周期内水轮发电机组关键部件的振动幅值功率谱；
S400,基于所述频率响应特性，将所述振动幅值功率谱作为输入激励加载到所述有限元模型中，疲劳仿真分析得到所述关键部件的疲劳寿命；
S500,设定疲劳寿命阈值，所述疲劳寿命阈值小于所述疲劳寿命；
S600,将以预设倍数调整振动幅值后的振动幅值功率，作为输入激励加载到所述有限元模型中，仿真求解得到所述关键部件的输出疲劳寿命；当所述输出疲劳寿命与所述疲劳寿命阈值的比值在预设范围x时，设定其对应的振动幅值为预警幅值。
在一些实施例中，步骤S100中，水轮发电机组关键部件可以为，水轮机组长时间运行后会产生疲劳破坏或者容易产生裂纹或者断裂的任意关键部件。例如水轮发电机组的主要承载部件：混流式水轮机组机架、混流式水轮机组顶盖、贯流式水轮机组水导轴承或贯流式水轮机组转轮室等。其中，在混流式水轮机组运行过程中，顶盖主要承受下水流压力、轴系不平衡力以及承受机组附件重量；而机架则承受整个转动部分的重量及水的推力以及自身重量。
在一些实施例中，所述建立水轮发电机组关键部件的有限元模型具体包括：
通过网格划分和单元质量控制建立水轮发电机组关键部件的有限元模型；
仿真计算：计算所述关键部件的初始有限元模型的仿真频率和仿真振形，得到仿真结果；
现场试验：获取所述关键部件的现场频率和现场振形，得到现场试验结果；
比对结果：比对仿真结果和现场试验结果，当仿真结果与现场结果不一致时，修正所述初始有限元模型，并重复所述仿真计算和现场试验的步骤，直至仿真结果与现场试验结果接近。通过将现场试验获取的结果与仿真计算的结果进行比对，能够尽可能提高有限元模型的准确度。
在一些实施例中，可以以关键部件的图纸，结合实际测绘尺寸，尽量模拟实物的机械结构，根据HYPERMESH提供的单元特性预先建立实体模型。单元质量控制可以理解为，对网格划分的单元进行质量检查。网格单元的质量检查可以包括对雅克比系数、翘曲度、长宽比等参数严格按照要求取值。
在一些实施例中，仿真结果可以包括关键部件的振动幅值和应力分布云图。
应当说明的是，仿真计算和现场试验时施加的参数是相同的。
在一些实施例中，在比对结果时，可以根据实际经验和具体的设计需求确定仿真结果与现场试验结果的接近程度。
在一些实施例中，所述比对结果的步骤之后还包括：
计算所述比对结果所得的有限元模型在典型工况下的振动幅值和振动应力分布，得到典型工况仿真结果；
统计典型工况仿真结果振动应力值较大的部位，识别易发生失效的潜在危险点；
获取所述关键部件的中心部位在典型工况下的振动幅值和振动应力分布，得到典型工况现场试验结果；
统计分析所述典型工况现场试验结果，得到所述关键部件的极端工况；
计算所述比对结果所得的有限元模型的潜在危险点在极端工况下的振动幅值和振动应力分布，得到极端工况仿真结果；
获取所述关键部件的潜在危险点在极端工况下的振动幅值和振动应力分布，得到极端工况现场试验结果；
比对所述极端工况仿真结果和所述极端工况现场试验结果，当不一致时，修正所述比对结果所得的有限元模型，并重复所述仿真计算和现场试验的步骤，直至仿真结果与现场结果接近。通过将现场试验识别潜在危险点，并通过对潜在危险点进行获取和仿真，并将结果进行比对，能够尽可能提高有限元模型的精确度。
在一些实施例中，可以通过在关键部件的中心部位布置传感器，对各种典型工况下关键部件的振动幅值进行实时测量。可以通过在关键部件的潜在危险点布置应变片得到通过动态应变采集仪收集各种典型工况下危险点的振动幅值和振动应力波动。
在一些实施例中,步骤S200中，频率响应(FRF)可以理解为通过已知输入和已知响应得到有限元模型的特性，具体测量对象可以包括输入力的频谱和同步的振动响应。具体可以包括，在步骤S100得到的有限元模型的基础上，将关键部件的质量点作为载荷施加激励点，分别施加X方向、Y方向和Z方向的单位位移振动激励。模态频率范围设定为(0～100)Hz，阻尼系数按照测得的动态阻尼曲线输入，以此对有限元模型进行模态频率响应分析，获取有限元模型的激励频率(0-50HZ)与振动应力之间的传递函数。
在应用场景中，对于混流式水轮机组顶盖分析时，选择顶盖上筯板直角处应力较大部位的单元作为监测点，得到其应力随X、Y和Z方向单位振幅激励的变化，请参阅图2、图3和图4。可见，单元的应力响应随着X向振动频率的增加逐渐缓慢升高，而随着Y和Z向振动频率的增加逐渐缓慢降低。但是应力增加或降低都并不明显，说明应力在0-50HZ的位移激励范围内对激励频率并不是很敏感。这种不敏感，是由于结构在50HZ之内只有40HZ一个固有频率点，且40HZ的固有模态振形是一高一低两瓣形态，而激励是整个中心呈现整体上下同步振动，即激励振形与固有模态不一致引起的。
在一些实施例中，如图5所示，步骤S300中，所述编制典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱具体包括：
S310,统计各典型工况工作时间占机组服役寿命的比重；
S320,根据各典型工况的实际比重对测试得到各典型工况下的载荷谱进行重新组合；
S330,将重新组合的结果进行时频变换，编制得到所述典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱。
通过对各典型工况的振幅时域曲线进行重新组成，能够形成整个服役寿命周期内的振幅时域曲线，综合考虑随着机组的全生命周期，叠加水轮发电机组水头等环境的影响，关键部件出现振动、摆度值逐步上升的趋势等，得到全面可靠的振动相关数据，提高后续的疲劳寿命预测的准确性，进而提高最终得到的预警幅值的准确性。
在一些实施例中，典型工况可以包括静止工况、额定工况和叠加工况等。针对不同的关键部件，可以选择不同的典型工况，以提高有限元模型的准确性。
具体地，对于混流式水轮机组顶盖选择静止工况、额定工况和叠加工况等，以得到顶盖只受自身结构及附件重力影响下的应力和应变的分布、顶盖受自身结构及附件重力和轴套反支力及额定水压力作用的应力和应变的分布，以及顶盖在额定工况分析基础上振动源叠加转轮X方向力对顶盖的综合影响下的应力和应变的分布。
在应用场景中，各个典型工况的运行时间占总运行时间的比重可以例如图6所示。
在一些实施例中，在步骤S320中，所述根据各典型工况的实际比重对测试得到各典型工况下的载荷谱进行重新组合具体包括：
根据式进行重新组合。其中，n1为静止工况的运行时间；n2为额定工况的运行时间；n8为叠加工况的运行时间。
在应用场景中，重新组合得到的服役周期内关键部件顶盖的振动幅值时域图如图7所示。
在步骤S330中，时频变换为现有的时频变换(ALT，Alternating Frequency/TimeDomain Method)，本申请不涉及对现有的时频变换的改进。所得振动幅值功率谱为一种功率谱密度(PSD)，能够量度随机变量的均方值，统计工况，为水轮发电机组振动幅值自相关函数的傅里叶变换。功率谱密度(PSD)可以理解为单位频带内的“功率”，是关键部件结构在随机动态载荷激励下响应的统计结果，表现为一条功率谱密度值—频率值的关系曲线。振动幅值功率谱的密度能够描述随机振动，反映了随机振动1秒钟的情况下，每个频率成份的加速度大小及振动次数，通过该位移功率谱密度能够得到主要功率对应的频率。可以理解为，振动幅值功率谱能够反映每个频率成份对应的功率，即相当于知道每秒钟此频率成份的功，由于每次振动的功已知(由加速度和频率可推算出来)，可以推算出每秒钟该频率振动的次数。
在应用场景中，将机组服役寿命周期内关键部件顶盖的振动幅值时域信号进行信号处理，得到其振动幅值功率谱的密度如图8a，图8b和图8c所示。可见，顶盖的X和Y向振动幅值功率集中在1HZ以内的尖点，Z方向的振动幅值功率集中在2.3及16.8HZ的尖点，且其频率分析较宽。
在一些实施例中，在步骤S400中，请参阅图9和图10，疲劳仿真分析的步骤具体可以为：将所述有限元模型进行谐响应分析得到的应力传递函数乘以标准载荷功率谱密度，得到该有限元模型的应力功率谱密度，然后利用Dirlik法模拟得到应力概率密度函数(PSD)，最后根据Miner准则进行疲劳寿命分析。
具体地，假设系统符合线弹性假设，是系统应力响应与激励载荷之间的传递函数，S(f)为激励载荷的标准功率谱密度，由此计算出系统的应力响应功率谱密度G(f)如式(2)，结合i阶频率fi和应力响应功率谱得到i阶功应力率谱密度矩阵mi如式(3)，进而依据应力功率谱密度谱矩m2和m4求出每秒钟的应力循环次数近似值E[P]如式(4)，由此推算出加载时间T内的全部应力循环次数N如式(5)。
G(f)＝|H(f)|2×S(f)(2)，其中，G(f)为系统的应力响应功率谱密度；|H(f)|为系统的应力响应与激励载荷之间的传递函；S(f)为激励载荷的标准功率谱密度。
mi＝∫fi×G(f)df(3)，其中，mi为i阶功率谱密度矩阵；fi为i阶频率；df为对f的微分，f为频率。
其中，E[P]为每秒钟的应力循环次数近似值；m2为2阶应力功率谱密度谱矩；m4为4阶应力功率谱密度谱矩。
N＝E[P]×T(5)，其中，N为时间T内的全部应力循环次数。
某个应力范围所对应的循环次数，可利用式(6)进行计算。ni＝P(Si)×dS×N(6)，其中，Si是第i个应力范围；P(Si)为对应的概率密度函数值；×dS为对S的微分，S为应力。
应力范围的断裂循环次数可根据材料的曲线(也即水轮发电机组关键部件的材料)用式(7)进行计算。
其中，Ni＝为应力范围的断裂循环次数，K和m为材料的常数参数；利用Miner线性损伤法则，将各应力范围水平造成的损伤线性相加。一旦/>则发生疲劳失效，得到疲劳失效时的循环次数，并计算得到对应的应力范围，进而根据功率谱密度得到振动幅值。
在一些实施例中，在步骤S500中，所述疲劳寿命阈值根据所述疲劳寿命和服役全寿命设定。例如，对于混流式水轮机组，由于水轮机每年运行时间为6000小时，以50年作为混流式水轮机组的服役全寿命，混流式水轮机组顶盖的疲劳寿命阈值设置为1.5e9秒；1e10秒作为混流式水轮机组机架的疲劳寿命阈值。
在一些实施例中，在步骤S600中，所述输出疲劳寿命与所述疲劳寿命阈值的接近程度，可以根据输出疲劳寿命与所述疲劳寿命阈值的比值来确定。具体地，当所述疲劳寿命呈现无限寿命时，可以不断以预设倍数调整振动幅值，直至输出疲劳寿命与所述疲劳寿命阈值的比值x的数值满足，0.1<x<1.1。是否呈现无线寿命，可以根据材料的S-N曲线判断。
在应用场景中，对于关键部件混流式水轮机组顶盖的主要材料为Q235，弹性模量为2.10e5MPa，泊松比为0.3，抗拉强度大于为375MPa时，通过ANSYS NCODE疲劳分析软件中的材料参数设置面板，在材料数据管理器中输入Q235的弹性模量E和极限抗拉强度UTS，创建材料的S-N曲线，如图11所示。可见，该S-N曲线为两段线组合形式，其以钢材10E7次为无限寿命次数。因此，呈现出无限寿命的特点。
在应用场景中，经过计算，对于混流式水轮机组的顶盖，其X方向振动允许振幅为现有的6.8倍时，轴向垂直振动允许振幅为现有的7倍。
实施例
对型号为HLA855-LJ，额定功率为250MW的水电站水轮机(以下简称水轮机)的顶盖的预警振幅，进行计算。其主要运行工况参数如下表1所示。
表1水轮机的主要运行工况参数
| 水轮机型号 | HLA855-LJ | 发电机型号 | SF250 |
| 额定水头 | 128m | 额定容量 | 277.8MVA |
| 水头范围 | 97m～156.5m | 额定功率 | 250MW |
| 额定转速 | 166.67r/min | 额定电压 | 15.75KV |
| 上导瓦数 | 12 | 额定电流 | 10183A |
| 水导瓦数 | 10 | 发电机转子 | 560t |
| 固定导叶 | 23 | 推力瓦数 | 24 |
| 转轮叶片 | 15 | 活动导叶数 | 24 |
| 承重机架 | 下机架 | 结构型式 | 全伞式 |
| 转轮名义 | 5.05m | 流量 | 217.81m3/s |
1)水轮机刚强度分析结果
顶盖的位移结果分别如表2所示,应力统计分别如表3。其中安全系数按第四强度理论校核，取Q235屈服强度为235MPa，安全系数按进行计算。其中，σs为极限应力；σmax为许用应力。
表2顶盖位移分布汇总
| 静止工况 | 额定工况 | 叠加工况 | |
| X方向位移/mm | 0.01 | 0.33 | 0.47 |
| Y方向位移/mm | -0.036 | 1.26 | 1.45 |
| Z方向位移/mm | 0.01 | -0.45 | -0.51 |
| 总位移/mm | 0.04 | 1.26 | 1.45 |
表3顶盖刚强度分析较大应力部位汇总
2)水轮机模态分析结果
设置结构阻尼系数为0.025，采用Block Lanczos法提取模态，水轮机顶盖模型提取前6阶模态，前6阶固有频率及对应的振形如下表4所示。
表4顶盖前6阶频率及振形
3)水轮机顶盖振动应力分析结果
现有振动预警系统规定下：顶盖水平振动报警限值：90μm，顶盖垂直振动报警限值：110μm。然而，顶盖振幅最大的工况为变转速100％工况，X方向振幅达到175μm，主要振动频率为0.7HZ。已经超过报警规定值，具体如下表5所示。
表5顶盖振幅最大的工况
| X向振幅 | Y向振幅 | Z向振幅 | 振动频率 | |
| 变转速100％ | -161～175μm | -52～41μm | -12～10μm | 0.7HZ |
鉴于此，很有必要分析上述振动极限工况下顶盖的应力，评估极限振幅对顶盖造成的附加振动应力。分析时，仍然采用建立的顶盖动力学模型，分别将变转速100％工况的振幅作为激励位移源施加到动力学模型中，进行振动响应分析，单独分析振动引起的附加应力。在此基础上，考虑极限振动工况与额定工况同时作用下顶盖的综合应力水平。
分析结果如表6所示。可见：变转速100％工况下的振动应力为40MPa，最大应力部位在矩形板边缘。
表6振幅最大工况下的振动应力
| 总位移 | 最大应力 | 最大应力部位 | |
| 变转速100％ | 0.18mm | 40MPa | 矩形板 |
4)水轮机顶盖疲劳寿命分析结果
该水轮机顶盖疲劳寿命分析是基于有限元模型进行的，而基于前述的水轮机刚强度分析结果、水轮机模态分析结果和水轮机振动应力分析结果，能够从不同维度对有限元模型进行调整，得到准确度和精确度均较高的有限元模型。
经计算，顶盖的疲劳寿命至少为7e22秒，呈现无限寿命特点。
在上基础上，不断增加振幅进行多次计算，分析振幅对疲劳寿命的影响，X方向振幅对疲劳寿命的影响如表7所示。Z方向垂直振动的振幅对疲劳寿命的影响如表8所示。可见，X方向和Z方向振动的振幅对疲劳寿命有显著影响。由于水轮机每年运行时间为6000小时，以50年(约1.1e9秒)作为无限寿命假设，则其X方向振动允许振幅为现有的6.3倍，Z方向垂直振动允许振幅为现有的75倍。
表7 X方向振动对疲劳寿命的影响
| 振幅 | 1倍振幅 | 5倍振幅 | 9倍振幅 | 6倍振幅 | 7倍振幅 | 6.3倍振幅 |
| 寿命/秒 | 7e22 | 3.6e11 | 5e5 | 4.5e9 | 2.7e8 | 1.78e9 |
表8垂直振动对疲劳寿命的影响
| 振幅 | 1倍振幅 | 50倍振幅 | 60倍振幅 | 80倍振幅 | 75倍振幅 |
| 寿命/秒 | 7e22 | 6.8e14 | 6e11 | 4.1e8 | 1.6e9 |
5)评估结论
(1)顶盖满足刚强度要求，顶盖额定工况下最大应力为124MPa；额定工况最大应力为132MPa，现有的预警振幅的安全系数分别为1.9和1.8。
(2)顶盖具有良好的动力学设计。激励频率与顶盖的固有频率相隔较远，不会出现共振现象。测试结果显示，顶盖和的实测振动以4HZ以内的低频为主，其主要振动频率和转子转频显著相关，表现出明显的受迫振动特点，且振动频率与结构自身的固有频率相隔较远。
(3)顶盖振动对应力有较明显影响。顶盖振动最大的变转速100％工况下，其引起的振动应力为40MPa。
(4)由于全寿命周期以额定微小振动工况为主，超预警振动的时长占比很小。现有工况的振幅不会对结构寿命产生实质影响，顶盖满足无限寿命要求。
(5)以50年(约1.1e9秒)作为结构的设计寿命，则顶盖的X方向振动允许振幅为现有的6.3倍。Z方向垂直振动允许振幅为现有的75倍。
需要说明的是，本申请实施例的方法可以由单个设备执行，例如一台计算机或服务器等。本实施例的方法也可以应用于分布式场景下，由多台设备相互配合来完成。在这种分布式场景的情况下，这多台设备中的一台设备可以只执行本申请实施例的方法中的某一个或多个步骤，这多台设备相互之间会进行交互以完成所述的方法。
需要说明的是，上述对本申请的一些实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于上述实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
基于同一发明构思，与上述任意实施例方法相对应的，本申请还提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上任意一实施例所述的水轮发电机组关键部件的预警振幅的计算方法。
图12示出了本实施例所提供的一种更为具体的电子设备硬件结构示意图，该设备可以包括：处理器1010、存储器1020、输入/输出接口1030、通信接口1040和总线1050。其中处理器1010、存储器1020、输入/输出接口1030和通信接口1040通过总线1050实现彼此之间在设备内部的通信连接。
处理器1010可以采用通用的CPU(Central Processing Unit，中央处理器)、微处理器、应用专用集成电路(Application Specific Integrated Circuit，ASIC)、或者一个或多个集成电路等方式实现，用于执行相关程序，以实现本说明书实施例所提供的技术方案。
存储器1020可以采用ROM(Read Only Memory，只读存储器)、RAM(Random AccessMemory，随机存取存储器)、静态存储设备，动态存储设备等形式实现。存储器1020可以存储操作系统和其他应用程序，在通过软件或者固件来实现本说明书实施例所提供的技术方案时，相关的程序代码保存在存储器1020中，并由处理器1010来调用执行。
输入/输出接口1030用于连接输入/输出模块，以实现信息输入及输出。输入输出/模块可以作为组件配置在设备中(图中未示出)，也可以外接于设备以提供相应功能。其中输入设备可以包括键盘、鼠标、触摸屏、麦克风、各类传感器等，输出设备可以包括显示器、扬声器、振动器、指示灯等。
通信接口1040用于连接通信模块(图中未示出)，以实现本设备与其他设备的通信交互。其中通信模块可以通过有线方式(例如USB、网线等)实现通信，也可以通过无线方式(例如移动网络、WIFI、蓝牙等)实现通信。
总线1050包括一通路，在设备的各个组件(例如处理器1010、存储器1020、输入/输出接口1030和通信接口1040)之间传输信息。
需要说明的是，尽管上述设备仅示出了处理器1010、存储器1020、输入/输出接口1030、通信接口1040以及总线1050，但是在具体实施过程中，该设备还可以包括实现正常运行所必需的其他组件。此外，本领域的技术人员可以理解的是，上述设备中也可以仅包含实现本说明书实施例方案所必需的组件，而不必包含图中所示的全部组件。
上述实施例的电子设备用于实现前述任一实施例中相应的水轮发电机组关键部件的预警振幅的计算方法，并且具有相应的方法实施例的有益效果，在此不再赘述。
基于同一发明构思，与上述任意实施例方法相对应的，本申请还提供了一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令用于使所述计算机执行如上任一实施例所述的水轮发电机组关键部件的预警振幅的计算方法。
本实施例的计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。
上述实施例的存储介质存储的计算机指令用于使所述计算机执行如上任一实施例所述的水轮发电机组关键部件的预警振幅的计算方法，并且具有相应的方法实施例的有益效果，在此不再赘述。
所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本申请的范围(包括权利要求)被限于这些例子；在本申请的思路下，以上实施例或者不同实施例中的技术特征之间也可以进行组合，步骤可以以任意顺序实现，并存在如上所述的本申请实施例的不同方面的许多其它变化，为了简明它们没有在细节中提供。
另外，为简化说明和讨论，并且为了不会使本申请实施例难以理解，在所提供的附图中可以示出或可以不示出与集成电路(IC)芯片和其它部件的公知的电源/接地连接。此外，可以以框图的形式示出装置，以便避免使本申请实施例难以理解，并且这也考虑了以下事实，即关于这些框图装置的实施方式的细节是高度取决于将要实施本申请实施例的平台的(即，这些细节应当完全处于本领域技术人员的理解范围内)。在阐述了具体细节(例如，电路)以描述本申请的示例性实施例的情况下，对本领域技术人员来说显而易见的是，可以在没有这些具体细节的情况下或者这些具体细节有变化的情况下实施本申请实施例。因此，这些描述应被认为是说明性的而不是限制性的。
尽管已经结合了本申请的具体实施例对本申请进行了描述，但是根据前面的描述，这些实施例的很多替换、修改和变型对本领域普通技术人员来说将是显而易见的。例如，其它存储器架构(例如，动态RAM(DRAM))可以使用所讨论的实施例。
Claims (9)
1.一种水轮发电机组关键部件的预警振幅的计算方法，其特征在于，包括：
建立水轮发电机组关键部件的有限元模型；
基于所述有限元模型，分析所述关键部件在不同振动频率激励下的振动应力，得到所述关键部件的频率响应特性；
编制典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱，包括：统计各典型工况工作时间占机组服役寿命的比重；根据各典型工况的实际比重对测试得到各典型工况下的载荷谱进行重新组合；将重新组合的结果进行时频变换，编制得到所述典型工况下服役寿命周期内水轮发电机组关键部件的振动幅值功率谱；
基于所述频率响应特性，将所述振动幅值功率谱作为输入激励加载到所述有限元模型中，疲劳仿真分析得到所述关键部件的疲劳寿命；
设定疲劳寿命阈值，所述疲劳寿命阈值小于所述疲劳寿命；
将以预设倍数调整振动幅值后的振动幅值功率，作为输入激励加载到所述有限元模型中，仿真求解得到所述关键部件的输出疲劳寿命；当所述输出疲劳寿命接近所述疲劳寿命阈值时，设定其对应的振动幅值为预警幅值。
2.根据权利要求1所述的水轮发电机组关键部件的预警振幅的计算方法，其特征在于，所述疲劳寿命阈值根据所述疲劳寿命和服役全寿命设定；当所述疲劳寿命呈现无限寿命时，所述输出疲劳寿命与所述疲劳寿命阈值的比值x的数值满足，0.1<x<1。
3.根据权利要求1所述的水轮发电机组关键部件的预警振幅的计算方法，其特征在于，所述根据各典型工况的实际比重对测试得到各典型工况下的载荷谱进行重新组合具体包括：
根据式进行重新组合；其中，n1为静止工况的运行时间；n2为额定工况的运行时间；n8为叠加工况的运行时间；n为各种工况的运行时间之和。
4.根据权利要求1所述的水轮发电机组关键部件的预警振幅的计算方法，其特征在于，所述建立水轮发电机组关键部件的有限元模型具体包括：
通过建立实体模型、网格划分和单元质量控制建立水轮发电机组关键部件的初始有限元模型；
仿真计算：计算所述关键部件的初始有限元模型的仿真频率和仿真振形，得到仿真结果；
现场试验：获取所述关键部件的现场频率和现场振形，得到现场试验结果；
比对结果：比对仿真结果和现场试验结果，当仿真结果与现场结果不一致时，修正所述初始有限元模型，并重复所述仿真计算和现场试验的步骤，直至仿真结果与现场试验结果接近。
5.根据权利要求4所述的水轮发电机组关键部件的预警振幅的计算方法，其特征在于，所述比对结果的步骤之后还包括：
计算所述比对结果所得的有限元模型在典型工况下的振动幅值和振动应力分布，得到典型工况仿真结果；
统计典型工况仿真结果应力值较大的部位，识别易发生失效的潜在危险点；
获取所述关键部件的中心部位在典型工况下的振动幅值，得到典型工况现场试验结果；
统计分析所述典型工况现场试验结果，得到所述关键部件的极端工况；
计算所述比对结果所得的有限元模型的潜在危险点在极端工况下的振动幅值和振动应力分布，得到极端工况仿真结果；
获取所述关键部件的潜在危险点在极端工况下的振动幅值和振动应力分布，得到极端工况现场试验结果；
比对所述极端工况仿真结果和所述极端工况现场试验结果，当不一致时，修正所述比对结果所得的有限元模型，并重复所述仿真计算和现场试验的步骤，直至仿真结果与现场结果接近。
6.根据权利要求1所述的水轮发电机组关键部件的预警振幅的计算方法，其特征在于，所述典型工况包括静止工况、额定工况和叠加工况；所述关键部件包括混流式水轮机组机架、混流式水轮机组顶盖、贯流式水轮机组水导轴承或贯流式水轮机组转轮室。
7.根据权利要求6所述的水轮发电机组关键部件的预警振幅的计算方法，其特征在于，所述关键部件为混流式水轮机组顶盖，所述典型工况为静止工况、额定工况和叠加工况。
8.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如权利要求1至7任意一项所述的方法。
9.一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令用于使计算机执行权利要求1至7任一所述方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111437224.5A CN114136573B (zh) | 2021-11-29 | 2021-11-29 | 水轮发电机组关键部件的预警振幅的计算方法及相关设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111437224.5A CN114136573B (zh) | 2021-11-29 | 2021-11-29 | 水轮发电机组关键部件的预警振幅的计算方法及相关设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114136573A CN114136573A (zh) | 2022-03-04 |
| CN114136573B true CN114136573B (zh) | 2024-01-23 |
Family
ID=80389354
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111437224.5A Active CN114136573B (zh) | 2021-11-29 | 2021-11-29 | 水轮发电机组关键部件的预警振幅的计算方法及相关设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114136573B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115096534B (zh) * | 2022-06-24 | 2023-03-14 | 大连理工大学 | 一种基于无参考点分区测试的柔度面识别方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5388056A (en) * | 1991-07-02 | 1995-02-07 | Hitachi, Ltd. | Method and system for vibration test |
| JPH07248254A (ja) * | 1994-03-11 | 1995-09-26 | Toshiba Corp | 振動監視方法およびその装置 |
| CN103983467A (zh) * | 2014-05-12 | 2014-08-13 | 中国人民解放军空军工程大学 | 一种基于服役状态的单机服役使用寿命监控方法 |
| CN104268335A (zh) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | 微组装组件振动疲劳寿命预测方法和系统 |
| CN107064796A (zh) * | 2017-04-20 | 2017-08-18 | 哈尔滨理工大学 | 智能电机振动在线监测及预警系统 |
| CN110595709A (zh) * | 2019-08-14 | 2019-12-20 | 南京航空航天大学 | 一种透平发动机叶片许用振幅确定方法 |
| CN111365158A (zh) * | 2020-03-02 | 2020-07-03 | 东方电气集团东方电机有限公司 | 水轮机转轮实时状态评估和生命周期管理预测系统 |
| CN112257320A (zh) * | 2020-10-26 | 2021-01-22 | 西安石油大学 | 一种机载机箱振动疲劳特性分析的方法 |
| CN112855408A (zh) * | 2021-01-15 | 2021-05-28 | 湖南五凌电力科技有限公司 | 水轮发电机组预警方法以及预警装置 |
| CN113686378A (zh) * | 2021-07-06 | 2021-11-23 | 国家能源集团科学技术研究院有限公司 | 一种水轮机顶盖运行状态监测方法 |
-
2021
- 2021-11-29 CN CN202111437224.5A patent/CN114136573B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5388056A (en) * | 1991-07-02 | 1995-02-07 | Hitachi, Ltd. | Method and system for vibration test |
| JPH07248254A (ja) * | 1994-03-11 | 1995-09-26 | Toshiba Corp | 振動監視方法およびその装置 |
| CN103983467A (zh) * | 2014-05-12 | 2014-08-13 | 中国人民解放军空军工程大学 | 一种基于服役状态的单机服役使用寿命监控方法 |
| CN104268335A (zh) * | 2014-09-23 | 2015-01-07 | 工业和信息化部电子第五研究所 | 微组装组件振动疲劳寿命预测方法和系统 |
| CN107064796A (zh) * | 2017-04-20 | 2017-08-18 | 哈尔滨理工大学 | 智能电机振动在线监测及预警系统 |
| CN110595709A (zh) * | 2019-08-14 | 2019-12-20 | 南京航空航天大学 | 一种透平发动机叶片许用振幅确定方法 |
| CN111365158A (zh) * | 2020-03-02 | 2020-07-03 | 东方电气集团东方电机有限公司 | 水轮机转轮实时状态评估和生命周期管理预测系统 |
| CN112257320A (zh) * | 2020-10-26 | 2021-01-22 | 西安石油大学 | 一种机载机箱振动疲劳特性分析的方法 |
| CN112855408A (zh) * | 2021-01-15 | 2021-05-28 | 湖南五凌电力科技有限公司 | 水轮发电机组预警方法以及预警装置 |
| CN113686378A (zh) * | 2021-07-06 | 2021-11-23 | 国家能源集团科学技术研究院有限公司 | 一种水轮机顶盖运行状态监测方法 |
Non-Patent Citations (2)
| Title |
|---|
| Turbocharger Radial Turbine Response to Pulse Amplitude;Roberto Mosca;PROCEEDINGS PAPER;全文 * |
| 大型水轮机顶盖动态特性数值模拟与试验研究;谭文胜;万元;夏代雄;李崇仕;;水电能源科学(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114136573A (zh) | 2022-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101835974A (zh) | 用于确定风力涡轮机的动力系中的疲劳损伤的方法 | |
| CN104268335A (zh) | 微组装组件振动疲劳寿命预测方法和系统 | |
| CN107038311A (zh) | 一种风力机叶片共振式疲劳试验载荷配置方法 | |
| CN105205264B (zh) | 预测抽水蓄能发电电动机转子构件使用寿命的试验方法 | |
| CN116577051A (zh) | 一种虑及损伤等效的随机振动疲劳寿命分析方法 | |
| CN103983412A (zh) | 用于振动有限元模型修正航空电子设备工作模态测量方法 | |
| CN104134013B (zh) | 一种风力机叶片模态分析方法 | |
| CN103217282B (zh) | 一种基于等寿命原则的水平轴风力机叶片缩比模型疲劳测试方法 | |
| CN107545110A (zh) | 一种动态应力加速寿命试验剖面编制方法 | |
| CN114136573B (zh) | 水轮发电机组关键部件的预警振幅的计算方法及相关设备 | |
| CN102564750A (zh) | 一种水平轴风力机叶片模态测试结果修正方法 | |
| CN109614640A (zh) | 一种大型风电机组轮毂疲劳寿命预测方法及系统 | |
| CN109115408B (zh) | 一种基于离心力方程的大型水电机组动平衡试验方法 | |
| Bouajila et al. | Hydraulic phenomena frequency signature of Francis turbines operating in part load conditions | |
| KR100946278B1 (ko) | 풍력 발전 단지의 모델링 및 해석 파라미터 시험 방법 | |
| CN103196689A (zh) | 一种水平轴风力机叶片静力测试结果修正方法 | |
| CN115422649A (zh) | 吊挂结构的腐蚀疲劳寿命预测方法及相关设备 | |
| CN108535558B (zh) | 一种电力稳定器的测试方法及系统 | |
| CN109085375A (zh) | 一种旋转机械的键相脉冲信号监测方法及装置 | |
| Kang et al. | Tacholess estimation of time-varying dynamic coefficients of journal bearing based on the square-root cubature Kalman filter | |
| KR102509719B1 (ko) | 터빈-발전기 축계 비틀림 응력 평가 시스템 및 방법 | |
| CN106777654B (zh) | 干摩擦阻尼隔振器等效阻尼确定方法 | |
| CN103399234A (zh) | 一种电网功角暂态稳定在线预测方法 | |
| Xie et al. | Vibration source analysis and structural optimization design of rotary table based on OTPA | |
| CN114201810A (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 |