A method, system and computer device for updating the fatigue life of a dangerous part of a runner

Abstract

The present application belongs to the technical field of hydroelectric generating set state monitoring and life evaluation, and discloses a runner dangerous position fatigue life updating method, system and computer equipment, the method comprising: acquiring real-time monitoring response value and predicted response value in the current time window of the runner dangerous position respectively, and constructing a deviation index; judging whether the runner dangerous position triggers finite element directional calibration; if the runner dangerous position triggers finite element directional calibration, combining the calibration fatigue damage increment with the pre-calculated predicted fatigue damage increment, determining the effective fatigue damage increment of the current time window; if the runner dangerous position does not trigger finite element directional calibration, taking the pre-calculated predicted fatigue damage increment as the effective fatigue damage increment of the current time window; and calculating the updated residual life and life risk state of the runner dangerous position. The present application can continuously update the cumulative fatigue damage and residual life of the runner dangerous position under complex working conditions.

Description

Technical Field

[0001] This invention relates to the field of hydropower unit condition monitoring and life assessment technology, and in particular to a method, system and computer equipment for updating the fatigue life of dangerous parts of the turbine runner. Background Technology

[0002] The turbine runner is a critical load-bearing and energy conversion component in a pumped storage unit. Under complex operating conditions such as start-up and shutdown switching, load fluctuations, partial load operation, and hydraulic pulsation, local high-stress areas of the runner are prone to alternating load concentration, leading to fatigue damage accumulation in critical locations. As operating time increases, fatigue damage in these critical locations will continue to evolve, potentially causing crack initiation and propagation, thus affecting the safe and stable operation of the unit. Therefore, identifying fatigue damage conditions and updating remaining service life for critical parts of the runner is of great significance for condition-based maintenance and safety early warning of the unit.

[0003] Currently, the analysis of turbine runner fatigue life typically employs the finite element method (FEM) to establish a turbine runner structural model. Under several typical operating conditions, the stress response of critical components is calculated, and fatigue life is estimated using the SN curve method, local strain method, or Miner's linear cumulative damage theory. While these methods can effectively reflect the stress level and life characteristics of critical components under preset operating conditions, their analysis usually relies on offline modeling and preset boundary conditions. This makes it difficult to reflect the response deviations of critical components caused by sudden load changes, flow field pulsations, start-up / shutdown transitions, or abnormal vibrations during actual operation, thus hindering continuous updates to the fatigue life of critical components. For example, existing technology has conducted fatigue life prediction studies on pump-turbine runners in pumped storage power stations. The stress spectrum of the pump-turbine under various operating conditions is analyzed using the nominal stress method combined with a finite element model with preset boundaries, and the fatigue life is estimated using Miner's cumulative damage theory. The results show that the fatigue life of the runner under variable load conditions is significantly lower than that under steady-state conditions. Existing technologies simulate the dynamic stress distribution of the runner of a Francis-type pumped storage power station based on the finite element method. After offline modeling, the local strain method combined with fluid-structure interaction is used to estimate the initiation life of the runner crack. However, it still mainly relies on simulation boundary conditions and empirical assumptions.

[0004] To improve online assessment capabilities, existing technologies utilize real-time monitoring data, empirical models, or neural network models to predict equipment response or life parameters. These methods offer the advantage of fast online calculation speeds; however, when operating conditions change significantly, monitoring data distribution shifts, or abnormal responses occur in localized hazardous areas, relying solely on predictions can easily lead to discrepancies with the actual stress state. Particularly during fatigue life updates for hazardous areas, the lack of a re-analysis mechanism based on the degree of deviation makes it difficult to determine when finite element analysis needs to be reintroduced to correct the response of hazardous areas, thus affecting the accuracy and stability of life update results.

[0005] In summary, the existing technology has the following shortcomings: First, it mainly relies on offline finite element analysis results, making it difficult to dynamically correct the fatigue life of the dangerous parts of the turbine based on the real-time operating status; second, the existing online prediction methods are prone to deviation under complex and variable working conditions, and lack a reliable correction mechanism for updating the life of dangerous parts; third, there is a lack of a method to trigger finite element reanalysis based on the degree of monitoring deviation and to update the fatigue damage and remaining life of the dangerous parts of the turbine using the calibration results.

[0006] Therefore, how to provide a method, system, and computer equipment for updating the fatigue life of dangerous parts of a turbine is an urgent problem to be solved. Summary of the Invention

[0007] This invention provides a method, system, and computer device for updating the fatigue life of critical parts of a turbine runner, in order to solve the problems mentioned above in the prior art.

[0008] According to a first aspect of the present invention, a method for updating the fatigue life of a critical part of a turbine rotor is provided.

[0009] In one embodiment, the method for updating the fatigue life of the critical parts of the turbine runner includes:

[0010] A finite element reference model of the runner is established based on its characteristic parameters, and the dangerous parts of the runner are identified.

[0011] The real-time monitoring response value and the predicted response value of the dangerous parts of the turbine runner within the current time window are obtained respectively, and a deviation index is constructed based on the difference between the real-time monitoring response value and the predicted response value.

[0012] The deviation index is compared with a preset threshold, and the comparison result is used to determine whether the dangerous part of the runner triggers finite element orientation calibration.

[0013] If the dangerous part of the runner triggers finite element orientation calibration, obtain the calibration response result of the dangerous part of the runner within the current time window and calculate the calibration fatigue damage increment. Combine the calibration fatigue damage increment with the pre-calculated predicted fatigue damage increment to determine the effective fatigue damage increment of the current time window.

[0014] If the dangerous part of the runner does not trigger the finite element orientation calibration, the pre-calculated predicted fatigue damage increment will be used as the effective fatigue damage increment for the current time window.

[0015] The effective fatigue damage increment is added to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the dangerous parts of the turbine runner. Based on the updated cumulative fatigue damage value, the remaining life and life risk status of the dangerous parts of the turbine runner after the update are calculated.

[0016] According to a second aspect of the present invention, a fatigue life renewal system for critical parts of a turbine is provided.

[0017] In one embodiment, the fatigue life update system for critical parts of the turbine runner includes:

[0018] The benchmark modeling and positioning module is used to establish a finite element benchmark model of the runner based on its characteristic parameters and to determine the dangerous parts of the runner.

[0019] The response deviation assessment module is used to obtain the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and to construct a deviation index based on the difference between the real-time monitoring response value and the predicted response value.

[0020] The calibration trigger determination module is used to compare the deviation index with the preset threshold and determine whether the dangerous part of the runner triggers finite element orientation calibration based on the comparison result;

[0021] The fusion damage calculation module is used to obtain the calibration response results of the dangerous part of the runner within the current time window and calculate the calibration fatigue damage increment if the dangerous part of the runner triggers finite element orientation calibration. The effective fatigue damage increment of the current time window is determined by combining the calibration fatigue damage increment with the pre-calculated predicted fatigue damage increment.

[0022] The predicted damage carry-over module is used to use the pre-calculated predicted fatigue damage increment as the effective fatigue damage increment for the current time window if the dangerous part of the runner does not trigger the finite element orientation calibration.

[0023] The life risk status update module is used to add the effective fatigue damage increment to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the dangerous parts of the turbine. Based on the updated cumulative fatigue damage value, the remaining life and life risk status of the dangerous parts of the turbine are calculated.

[0024] According to a third aspect of the present invention, a computer device is provided.

[0025] In some embodiments, the computer device includes a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement the steps of the method described above.

[0026] According to a fourth aspect of the present invention, a computer program product is provided.

[0027] In one embodiment, the computer program product stores a computer program that, when executed by a processor, implements the steps of the method.

[0028] According to a fifth aspect of the present invention, a computer-readable storage medium is provided.

[0029] In one embodiment, a computer program is stored on the computer-readable storage medium, which, when executed by a processor, implements the steps of the above method.

[0030] The technical solutions provided by the embodiments of the present invention may include the following beneficial effects:

[0031] 1. The fatigue life update method for critical parts of a turbine rotor provided by this invention addresses the problems in existing technologies, such as the difficulty of dynamically correcting the fatigue life of critical parts of a turbine rotor based on real-time operating conditions in offline finite element analysis, the tendency for online prediction results to deviate from the actual response of critical parts under complex variable operating conditions, and the lack of a mechanism to trigger finite element re-analysis based on the degree of deviation and use calibration results for life update. By constructing a deviation index between the real-time monitoring value and the predicted value of the critical part, finite element directional calibration is triggered for the critical part based on the degree of deviation, and the fatigue damage increment in the current time window is corrected using the calibration results. This achieves dynamic updating of the cumulative fatigue damage and remaining life of the critical parts of the turbine rotor, improving the accuracy and stability of the life update results.

[0032] 2. This invention focuses on the fatigue life update of the dangerous parts of the turbine, which can more accurately reflect the life evolution process of the areas where fatigue damage is prone to occur and expand, and improve the pertinence of life update. By constructing a deviation index between the real-time monitoring value and the predicted value of the dangerous parts, and judging whether to trigger finite element directional calibration based on the degree of deviation, a calibration mechanism can be introduced in time when the predicted result deviates from the actual operating state, thereby improving the adaptability of life update under complex and variable working conditions.

[0033] 3. This invention adopts a finite element-oriented calibration method for hazardous parts, updating and solving only the local loads, boundary conditions, or local sub-models related to the fatigue damage evolution of hazardous parts. This helps reduce the computational scale and improve the efficiency of life update. Furthermore, based on the degree of deviation, the predicted damage increment and the finite element-calibrated damage increment are segmented or weighted, so that the effective fatigue damage increment of the current time window can be dynamically adjusted with the changes in the operating status of the hazardous parts. This improves the accuracy and stability of the cumulative fatigue damage and remaining life update results. As a result, it can continuously update the cumulative fatigue damage and remaining life of the turbine runner under complex operating conditions such as load fluctuations, start-stop switching, hydraulic pulsation changes, and response drift, providing a basis for unit status early warning and maintenance decisions.

[0034] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and are not intended to limit the invention. Attached Figure Description

[0035] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0036] Figure 1 This is a flowchart illustrating a method according to an exemplary embodiment;

[0037] Figure 2 This is a system schematic diagram illustrated according to an exemplary embodiment;

[0038] Figure 3 This is a schematic diagram of the structure of a computer device according to an exemplary embodiment;

[0039] Figure 4 This is a general flowchart illustrated according to an exemplary embodiment;

[0040] Figure 5 This is a schematic diagram illustrating the wheel structure and the location of dangerous parts according to an exemplary embodiment;

[0041] Figure 6 This is a schematic diagram illustrating a finite element orientation calibration triggering determination based on a deviation index, according to an exemplary embodiment.

[0042] Figure 7 This is a schematic diagram illustrating the incremental correction of fatigue damage and the update of remaining life at dangerous locations according to an exemplary embodiment. Detailed Implementation

[0043] The following description and accompanying drawings fully illustrate specific embodiments described herein to enable those skilled in the art to practice them. Some portions and features of certain embodiments may be included in or replace portions and features of other embodiments. The scope of the embodiments herein includes the entire scope of the claims and all available equivalents thereof. The various embodiments described herein are presented in a progressive manner, with each embodiment focusing on its differences from other embodiments; similar or identical parts between embodiments can be referred to interchangeably.

[0044] The modules in the apparatus or system of this application can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device in hardware form, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each module.

[0045] Where there is no conflict, the embodiments and features in the embodiments of the present invention can be combined with each other.

[0046] Figure 1 An embodiment of a method for updating the fatigue life of a critical part of a turbine according to the present invention is shown.

[0047] In this optional embodiment, the fatigue life update method for the critical parts of the turbine runner includes:

[0048] S101. Establish a finite element reference model of the impeller based on its characteristic parameters, and determine the dangerous parts of the impeller;

[0049] S102. Obtain the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and construct a deviation index based on the difference between the real-time monitoring response value and the predicted response value.

[0050] S103. Compare the deviation index with the preset threshold, and determine whether the dangerous part of the runner triggers finite element orientation calibration based on the comparison result;

[0051] S104. If the dangerous part of the runner triggers finite element orientation calibration, obtain the calibration response result of the dangerous part of the runner within the current time window and calculate the calibration fatigue damage increment. Combine the calibration fatigue damage increment with the pre-calculated predicted fatigue damage increment to determine the effective fatigue damage increment of the current time window.

[0052] S105. If the dangerous part of the runner does not trigger the finite element orientation calibration, the pre-calculated predicted fatigue damage increment will be used as the effective fatigue damage increment of the current time window.

[0053] S106. The effective fatigue damage increment is added to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the dangerous parts of the turbine runner. The remaining life and life risk status of the dangerous parts of the turbine runner are calculated based on the updated cumulative fatigue damage value.

[0054] In this optional embodiment, the step of establishing a finite element reference model of the runner based on its characteristic parameters and determining the dangerous parts of the runner includes:

[0055] Based on the structural characteristic parameters, material parameters and installation constraints of the runner, a finite element reference model of the runner is established;

[0056] Based on the structural stress analysis results, historical failure locations, and fatigue sensitivity assessment results, the transition zone connecting the blade root and the upper crown in the finite element reference model of the runner is identified as the critical part of the runner.

[0057] It should be noted that the dangerous areas of the runner are the transition zone between the blade root and the upper crown, the transition zone between the blade root and the lower ring, the area adjacent to the weld, or the area with a high incidence of historical cracks.

[0058] The baseline response of the critical parts of the runner under various typical operating conditions was calculated based on the finite element baseline model.

[0059] In this optional embodiment, the reference response includes stress response, strain response, displacement response, and vibration response.

[0060] In this optional embodiment, the structural feature parameters include the runner structure dimensions, blade geometry parameters, crown structure, and ring structure;

[0061] The material parameters include elastic modulus, Poisson's ratio, density, and stress-life curve parameters;

[0062] The typical operating conditions include rated operating conditions, partial load operating conditions, start-stop transition operating conditions, load fluctuation operating conditions, and flow field pulsation operating conditions.

[0063] In this optional embodiment, the step of obtaining the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and constructing a deviation index based on the difference between the real-time monitoring response value and the predicted response value, includes:

[0064] Obtain the real-time monitoring response value of the dangerous parts of the turbine within the current time window;

[0065] The current typical operating condition parameters and historical response sequences are used as inputs to the preset response prediction model. The predicted response values ​​of the dangerous parts of the turbine runner within the current time window are output, and the predicted fatigue damage increment is obtained through fatigue damage calculation.

[0066] The deviation index for the current time window is calculated based on the difference between the real-time monitoring response value and the predicted response value of the dangerous parts of the runner within the current time window.

[0067] In this optional embodiment, the formula for calculating the deviation index of the current time window is:

[0068] ;

[0069] In the formula, i is the sampling point number within the current time window. This represents the real-time monitoring response value of the dangerous part of the rotor at the i-th sampling point within the current time window. This is the predicted response value corresponding to the i-th sampling point within the current time window. This represents the number of sampling points within the current time window. To prevent tiny positive numbers with a denominator of zero, This is the deviation indicator for the current time window.

[0070] In this optional embodiment, comparing the deviation index with a preset threshold and determining whether the dangerous part of the runner triggers finite element orientation calibration based on the comparison result includes:

[0071] The deviation index of the current time window is compared with the preset first calibration trigger threshold and the second forced calibration threshold, respectively;

[0072] When the deviation index is less than or equal to the first calibration trigger threshold, the finite element orientation calibration of the dangerous part of the runner will not be triggered.

[0073] When the deviation index is greater than the first calibration trigger threshold and less than or equal to the second mandatory calibration threshold, finite element orientation calibration is performed on the dangerous parts of the runner.

[0074] When the deviation index exceeds the second mandatory calibration threshold, finite element orientation calibration is performed on the dangerous parts of the runner.

[0075] In this optional embodiment, if the dangerous part of the turbine triggers finite element orientation calibration, obtaining the calibration response result of the dangerous part of the turbine within the current time window and calculating the calibration fatigue damage increment includes:

[0076] Based on the finite element reference model of the turbine, which includes the dangerous part of the turbine and the surrounding transition structure, a local sub-model of the dangerous part is constructed. The geometric parameters, material parameters and initial boundary conditions of the local sub-model of the dangerous part are derived from the finite element reference model of the turbine, and the mesh of the local sub-model of the dangerous part is refined in the dangerous part of the turbine and its adjacent area.

[0077] Based on the global response results of the finite element reference model of the runner under the corresponding operating conditions, the boundary transfer conditions of the local sub-model of the dangerous part are determined by at least one of the following methods: displacement constraint at the truncation boundary, equivalent nodal force or interpolated load.

[0078] Based on the real-time operating conditions within the current time window, update the local load spectrum parameters, local boundary conditions, and hydraulic pulsation loads acting on the local sub-model of the dangerous area.

[0079] Based on the updated local load spectrum parameters, local boundary conditions, and hydraulic pulsation loads, the local sub-model of the dangerous part is solved again by finite element method, and the calibration response of the dangerous part of the runner within the current time window is output. The calibration response includes calibration stress response and calibration strain response.

[0080] The calibration response is sequentially subjected to equivalent stress conversion, rainflow counting, stress-life curve conversion, and linear cumulative damage calculation to obtain the calibration fatigue damage increment.

[0081] In this optional embodiment, the step of adding the effective fatigue damage increment to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the critical parts of the turbine runner, and calculating the updated remaining life and life risk status of the critical parts of the turbine runner based on the updated cumulative fatigue damage value includes:

[0082] The effective fatigue damage increment within the current time window is added to the cumulative fatigue damage value at the end of the previous time window to obtain and update the cumulative fatigue damage value at the end of the current time window.

[0083] Based on the stress-life curve method and linear cumulative damage theory, the remaining life of the critical parts of the turbine after the update is calculated based on the updated cumulative fatigue damage value;

[0084] When the updated cumulative fatigue damage value reaches the preset damage threshold, or the updated remaining lifespan is lower than the preset safety threshold, the life risk status of the dangerous parts of the rotor is output.

[0085] In this optional embodiment, the formula for calculating the remaining lifespan of the turbine after the critical parts are updated is as follows:

[0086] RUL=(D lim -D k ) / v D ;

[0087] In the formula, D lim As a preset damage threshold, v D D represents the average damage growth rate within a preset number of continuous time windows. k RUL is the cumulative fatigue damage value updated at the end of the current time window, and RUL is the remaining life of the critical parts of the turbine after the update.

[0088] Figure 2 An embodiment of a fatigue life renewal system for critical parts of a rotor according to the present invention is shown.

[0089] In this optional embodiment, the fatigue life update system for the critical parts of the turbine runner includes:

[0090] The benchmark modeling and positioning module 201 is used to establish a finite element benchmark model of the runner based on its characteristic parameters and to determine the dangerous parts of the runner.

[0091] The response deviation assessment module 202 is used to obtain the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and to construct a deviation index based on the difference between the real-time monitoring response value and the predicted response value.

[0092] The calibration trigger determination module 203 is used to compare the deviation index with the preset threshold and determine whether the dangerous part of the runner triggers finite element orientation calibration based on the comparison result.

[0093] The fusion damage calculation module 204 is used to obtain the calibration response result of the dangerous part of the runner within the current time window and calculate the calibration fatigue damage increment if the dangerous part of the runner triggers finite element orientation calibration. It combines the calibration fatigue damage increment with the pre-calculated predicted fatigue damage increment to determine the effective fatigue damage increment of the current time window.

[0094] The predicted damage module 205 is used to use the pre-calculated predicted fatigue damage increment as the effective fatigue damage increment of the current time window if the dangerous part of the runner does not trigger the finite element orientation calibration.

[0095] The life risk status update module 206 is used to add the effective fatigue damage increment to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the dangerous parts of the runner, and calculate the updated remaining life and life risk status of the dangerous parts of the runner based on the updated cumulative fatigue damage value.

[0096] In this optional embodiment, the step of establishing a finite element reference model of the runner based on its characteristic parameters and determining the dangerous parts of the runner includes:

[0097] Based on the structural characteristic parameters, material parameters and installation constraints of the runner, a finite element reference model of the runner is established;

[0098] Based on the structural stress analysis results, historical failure locations, and fatigue sensitivity assessment results, the transition zone connecting the blade root and the upper crown in the finite element reference model of the runner was identified as the critical part of the runner.

[0099] The baseline response of the critical parts of the runner under various typical operating conditions was calculated based on the finite element baseline model.

[0100] In this optional embodiment, the reference response includes stress response, strain response, displacement response, and vibration response.

[0101] In this optional embodiment, the structural feature parameters include the runner structure dimensions, blade geometry parameters, crown structure, and ring structure;

[0102] The material parameters include elastic modulus, Poisson's ratio, density, and stress-life curve parameters;

[0103] The typical operating conditions include rated operating conditions, partial load operating conditions, start-stop transition operating conditions, load fluctuation operating conditions, and flow field pulsation operating conditions.

[0104] In this optional embodiment, the step of obtaining the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and constructing a deviation index based on the difference between the real-time monitoring response value and the predicted response value, includes:

[0105] Obtain the real-time monitoring response value of the dangerous parts of the turbine within the current time window;

[0106] The current typical operating condition parameters and historical response sequences are used as inputs to the preset response prediction model. The predicted response values ​​of the dangerous parts of the turbine runner within the current time window are output, and the predicted fatigue damage increment is obtained through fatigue damage calculation.

[0107] The deviation index for the current time window is calculated based on the difference between the real-time monitoring response value and the predicted response value of the dangerous parts of the runner within the current time window.

[0108] To facilitate understanding of the above technical solutions of the present invention, the following further explains the above technical solutions of the present invention from the perspective of architecture and principle, as follows:

[0109] This invention proposes a method for updating the fatigue life of critical parts of a turbine runner. Addressing the shortcomings of existing technologies that primarily rely on offline finite element analysis (FEM) for updating the fatigue life of critical parts of turbine runners, which struggles to reflect complex operating conditions in a timely manner and lacks an effective correction mechanism for sudden changes in operating conditions or response drift, this invention constructs a dynamic updating method encompassing critical part identification, monitoring and prediction deviation calculation, finite element-based directional trigger calibration, damage increment correction, and remaining life updating. This method enables continuous updating of the fatigue damage state and remaining life of critical parts of the turbine runner, specifically including:

[0110] Based on the runner's structural parameters, material parameters, connection boundary conditions, and typical operating conditions, a finite element reference model of the runner was established. Specifically, the runner's finite element reference model includes the upper crown, lower ring, blades, hub, and connection transition zone; material parameters include elastic modulus, Poisson's ratio, density, and SN curve (stress-life curve) parameters; connection boundary conditions include runner-spindle connection constraints, upper crown constraints, lower ring constraints, and equivalent constraints for welded connections; typical operating conditions include rated operating conditions, partial load operating conditions, start-stop transition operating conditions, and load fluctuation operating conditions. After local mesh refinement of the candidate critical areas, centrifugal loads, hydrostatic pressure, operating pressure field, and hydraulic pulsation loads were applied to solve for the reference stress response, reference strain response, displacement response, and vibration response of the critical parts. Based on historical failure locations, stress concentration analysis results, and fatigue sensitivity analysis results, the critical parts of the runner were determined. Critical parts were identified as the transition zone between the blade root and the upper crown, the transition zone between the blade root and the lower ring, the area adjacent to the weld, or areas with a high incidence of historical cracks. Based on the finite element baseline model, the baseline response of the critical part (the critical part is the critical part of the turbine, hereinafter referred to as the critical part) under preset typical working conditions is calculated. The baseline response includes stress response, strain response, displacement response and vibration response, which are used for subsequent fatigue damage updates of the critical part.

[0111] Acquire real-time monitoring and predicted values ​​for hazardous locations within the current time window. Within the current time window, collect data on surface strain, unit vibration response, flow field pulsation pressure, speed, load, load change rate, and start-up / shutdown status for the hazardous location. The prediction model (a response prediction model, or simply prediction model) takes current operating parameters, historical response sequences, and baseline responses as inputs and outputs predicted response values ​​for the hazardous location within the current time window. Predict fatigue damage increment ΔD. p The predicted response value is obtained through equivalent stress conversion, rainflow counting, SN curve conversion, and Miner linear cumulative damage calculation. The current time window is a fixed duration window of 5 to 60 seconds, or an integer multiple of one blade frequency cycle of the impeller.

[0112] The prediction model is a time series regression model, which is trained using historical operating data. The training samples include operating parameters, historical response sequences, finite element reference responses, and corresponding measured responses. The training objective is to reduce the error between the predicted response and the measured response.

[0113] Construct a monitoring-prediction deviation index for hazardous areas. Based on the difference between the real-time monitoring values ​​and predicted values ​​of hazardous areas within the current time window, a monitoring-prediction deviation index E is constructed to characterize the degree of deviation between the predicted results and the actual operating status of the hazardous areas. The deviation index E is calculated using the following formula:

[0114] ;

[0115] In the formula, i is the sampling point number within the current time window. This represents the real-time monitoring response value of the dangerous part of the rotor at the i-th sampling point within the current time window. This is the predicted response value corresponding to the i-th sampling point within the current time window. This represents the number of sampling points within the current time window. To prevent tiny positive numbers with a denominator of zero, The deviation index is the current time window; the real-time monitoring response value and the predicted response value are taken as the equivalent variation amplitude value, equivalent stress amplitude value, and dominant frequency vibration response amplitude value of the dangerous part.

[0116] The determination of whether to trigger finite element orientation calibration is based on the deviation index. A first calibration trigger threshold T1 and a second mandatory calibration threshold T2 are set, with T2 > T1. The deviation index E is compared with thresholds T1 and T2, and the result determines whether finite element orientation calibration is triggered for the hazardous area. The specific judgment rules are as follows: When E ≤ T1, the deviation between the current hazardous area prediction result and the actual operating state is considered small, and finite element orientation calibration is not triggered; the prediction result is directly used as the input for fatigue damage calculation in the current time window. When T1 < E ≤ T2, the deviation between the current hazardous area prediction result and the actual operating state is considered moderate, and finite element orientation calibration is triggered; the predicted damage increment is corrected based on the calibration result. When E > T2, the deviation between the current hazardous area prediction result and the actual operating state is considered large, and finite element orientation calibration is triggered; the calibrated finite element result is used as the input for fatigue damage calculation in the current time window. Thresholds T1 and T2 are determined based on the statistical results of the deviation distribution of hazardous areas under historical normal and abnormal operating conditions.

[0117] Finite element orientation calibration is performed on hazardous areas. When the deviation index meets the triggering condition, finite element orientation calibration is performed on the hazardous areas based on the real-time operating condition parameters within the current time window to obtain the calibration response results of the hazardous areas within the current time window, and the calibration fatigue damage increment is calculated based on the calibration response results. Finite element orientation calibration includes the following:

[0118] 1. Based on the finite element reference model of the turbine wheel, a local sub-model of the hazardous area and its surrounding transition structure is constructed. The geometric parameters, material parameters, element types, and initial boundary conditions of the local sub-model of the hazardous area are derived from the finite element reference model of the turbine wheel, and local mesh refinement is performed in the hazardous area and its adjacent region.

[0119] 2. Based on the global displacement response, stress response, or load mapping results of the finite element reference model of the runner under the corresponding operating conditions, determine the displacement constraints, equivalent nodal forces, or interpolated loads at the cut-off boundary of the local sub-model of the dangerous part, so that a boundary and load transfer relationship is formed between the local sub-model of the dangerous part and the finite element reference model of the runner;

[0120] 3. Based on the rotational speed, load, load change rate, vibration response, and flow field pulsating pressure within the current time window, update the local load spectrum parameters, local boundary conditions, and hydraulic pulsating loads acting on the local sub-model of the dangerous area;

[0121] 4. Based on the updated local load spectrum parameters, local boundary conditions, and hydraulic pulsation loads, the local sub-model of the dangerous area is solved again using finite element method.

[0122] 5. Output the calibration stress response, calibration strain response, calibration equivalent fatigue characterization quantity, and calibration fatigue damage increment of the critical part within the current time window.

[0123] The calibration response results include the calibration stress response or calibration strain response within the current time window of the critical area. When the output is a calibration strain response, it is converted into a calibration equivalent stress response; when the output is a calibration stress response, it is directly used for fatigue damage calculation. Rainflow counting is performed on the calibration equivalent stress response to obtain the stress amplitude and cycle number of each stress cycle. Then, the corresponding allowable number of cycles is determined according to the SN curve, and the calibration fatigue damage increment ΔD is calculated according to Miner's linear cumulative damage theory. f .

[0124] Determine the effective fatigue damage increment within the current time window. This is done by combining the fatigue damage increment ΔD obtained from the prediction model and the fatigue damage increment obtained from the finite element orientation calibration within the current time window. Determine the effective fatigue damage increment within the current time window. .

[0125] When E≤T1: When T1 < E ≤ T2: The weighting coefficient α is calculated using the following formula: When E > T2: ;in: To predict the increment of fatigue damage, To calibrate the fatigue damage increment, This represents the effective increase in fatigue damage.

[0126] Using the above method, the effective fatigue damage increment within the current time window is determined based on the degree of deviation between the predicted results of hazardous locations and the actual operating conditions. The effective fatigue damage increment determined within the current time window is then used as the basis for this determination. This value is added to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value for the critical area. : ;in, This represents the cumulative fatigue damage value of the hazardous component at the end of the previous time window. Based on the SN curve method and Miner's linear cumulative damage theory, the remaining life (RUL) of the hazardous component after the update is calculated. When the updated cumulative fatigue damage reaches a preset damage threshold, or the remaining life is lower than a preset safety threshold, the life risk status of the hazardous component is output for unit status early warning and maintenance decision-making. The formula for calculating the remaining life (RUL) of the hazardous component after the update is: RUL = (D lim -D k ) / v D， Among them, D lim As a preset damage threshold, v D The average damage growth rate is defined as the average damage growth rate within a recently preset number of consecutive time windows. This preset number of consecutive time windows includes the current time window and one or more consecutive time windows preceding it. The preset number can be pre-set based on the time window length, monitoring sampling frequency, and lifetime update accuracy, thus clearly defining the range of the average damage growth rate. k RUL is the cumulative fatigue damage value updated at the end of the current time window, and RUL is the remaining life of the critical parts of the turbine after the update.

[0127] The present invention will be further described in detail below with reference to specific embodiments.

[0128] This embodiment takes the runner of a pumped storage unit as the object and adopts the fatigue life update method for dangerous parts of the runner proposed in this invention to dynamically update the fatigue damage state and remaining life of dangerous parts of the runner during actual operation.

[0129] like Figure 4 As shown, the method in this embodiment includes: identifying the dangerous parts of the turbine, establishing a finite element reference model of the dangerous parts, obtaining the real-time monitoring value and predicted value of the dangerous parts within the current time window, calculating the monitoring-prediction deviation index E, determining whether to trigger finite element orientation calibration based on the deviation index, determining the effective fatigue damage increment of the current time window, and updating the cumulative fatigue damage and remaining life of the dangerous parts.

[0130] Example 1: The method of the present invention is described using the transition area between the root of the rotor blade and the upper crown as the dangerous part.

[0131] The first step is to identify the critical areas and establish a finite element baseline model. Based on the runner's structural dimensions, blade geometry, crown structure, annular structure, material parameters, and installation constraints, a three-dimensional finite element baseline model of the runner is established. Material parameters include elastic modulus, Poisson's ratio, density, and the SN curve parameters required for fatigue life analysis. In this embodiment, the elastic modulus of the runner material is taken as 2.06 × 10⁻⁶. 5 MPa, Poisson's ratio is taken as 0.30, and density is taken as 7850 kg / m³ 3 Operating conditions include rated operating conditions, partial load operating conditions, start-stop transition operating conditions, load fluctuation operating conditions, and flow field pulsation operating conditions.

[0132] Based on the structural stress analysis results, historical failure locations, and fatigue sensitivity assessment results, the transition zone connecting the blade root and the upper crown was identified as a critical area. Specifically, as follows... Figure 5 As shown, the equivalent stress, principal stress, displacement response, and vibration response of the hazardous part under the above typical operating conditions are calculated to form finite element reference response data for the hazardous part, which is used for subsequent directional calibration and life update. In this embodiment, the finite element reference model calculates that the equivalent stress amplitude of the hazardous part is 92 MPa under rated operating conditions, 118 MPa under partial load conditions, 135 MPa under start-stop transition conditions, and 126 MPa under load fluctuation conditions.

[0133] The second step involves acquiring real-time monitoring and predicted values ​​for the hazardous area within the current time window. During unit operation, a 10-second time window is used as the current time window. Within each time window, surface strain, unit vibration response, speed, and load signals for the hazardous area are collected. Surface strain at the hazardous area is acquired by strain sensors positioned near the transition zone between the blade root and the upper crown. The unit vibration response is acquired by a vibration monitoring device, and the speed and load are acquired by the unit operation monitoring system. Fatigue damage at the hazardous area is related not only to structural stress but also to the combined effects of flow field pulsation pressure, load fluctuation frequency, start-stop cycles, and vibration response. In this embodiment, the unit load increases from 80% to 92% of the rated load within the current time window, and the equivalent amplitude value for the hazardous area is 168 με.

[0134] Based on the current operating parameters and historical response sequences, a pre-defined hazardous part response prediction model is used to output the predicted response value of the hazardous part within the current time window. The predicted fatigue damage increment ΔD is then calculated using fatigue damage analysis. pIn this embodiment, the prediction model, based on the current rotational speed, load, load change rate, historical strain response sequence, and finite element reference response, obtains a predicted equivalent amplitude value of 142με for the critical part. The predicted response is then subjected to equivalent stress conversion, rainflow counting, SN curve conversion, and Miner linear cumulative damage calculation to obtain the predicted fatigue damage increment ΔD. p =1.6×10 -7 .

[0135] The third step is to calculate the monitoring-prediction deviation index for hazardous areas. Based on the difference between the real-time monitoring response value and the predicted response value of the hazardous area within the current time window, the monitoring-prediction deviation index E is calculated to characterize the degree of deviation between the current prediction result and the actual operating state. In this embodiment, the real-time monitoring response value and the predicted response value are taken as the equivalent amplitude values ​​of the hazardous area. The deviation index E is calculated using the normalized error average. Based on the real-time monitoring equivalent amplitude value of 168 με and the predicted equivalent amplitude value of 142 με, the deviation index E for the current time window is calculated to be 0.154.

[0136] The fourth step is to determine whether to trigger finite element orientation calibration based on the deviation index. A first calibration trigger threshold T1 and a second forced calibration threshold T2 are set, where T2 > T1. In this embodiment, thresholds T1 and T2 are determined based on the statistical results of the deviation distribution during the unit's historical normal operation and abnormal operation phases, where T1 = 0.08 and T2 = 0.20, as detailed below. Figure 6 As shown.

[0137] The deviation index E obtained in the current time window is compared with the thresholds T1 and T2, and processed according to the following rules: when E≤T1, it is considered that the deviation between the prediction result of the current dangerous part and the actual operating state is small, the finite element orientation calibration is not triggered, and the predicted fatigue damage increment ΔD is directly used. p The effective fatigue damage increment is used as the current time window. When T1 < E ≤ T2, it is determined that there is a moderate deviation between the predicted result of the current dangerous part and the actual operating state, triggering finite element orientation calibration, and the predicted damage increment is corrected according to the calibration result. When E > T2, it is determined that there is a large deviation between the predicted result of the current dangerous part and the actual operating state, triggering finite element orientation calibration, and the calibrated fatigue damage increment ΔD is directly used. f This serves as the effective fatigue damage increment for the current time window. In this embodiment, E = 0.154, satisfying T1 < E ≤ T2, thus triggering finite element orientation calibration and correcting the predicted fatigue damage increment.

[0138] The fifth step is to perform finite element orientation calibration of the hazardous parts. When the deviation index meets the triggering conditions, finite element orientation calibration is performed on the hazardous parts based on the real-time operating parameters within the current time window. This calibration does not recalculate the entire finite element model, but only updates and resolves the local sub-model containing the hazardous parts.

[0139] Specifically, based on the rotational speed, load, load change rate, vibration response, and flow field pulsating pressure within the current time window, the local load spectrum parameters, local boundary conditions, and hydraulic pulsating loads related to the critical area are updated. Then, the local sub-model of the critical area is solved again using the finite element method, outputting the calibration stress response and calibration strain response of the critical area within the current time window. The calibration response is then subjected to equivalent stress conversion, rainflow counting, SN curve conversion, and Miner linear cumulative damage calculation to obtain the calibration fatigue damage increment ΔD. f In this embodiment, the calibration fatigue damage increment ΔD obtained from finite element orientation calibration is... f =2.8×10 -7 .

[0140] Step 6: Determine the effective fatigue damage increment within the current time window. This is based on the predicted fatigue damage increment ΔD within the current time window. p and calibration fatigue damage increment ΔD f Determine the effective fatigue damage increment ΔD for the current time window. c Since T1 < E ≤ T2 in this embodiment, the formula is used to calculate α = (T2 - E) / (T2 - T1). Substituting T1 = 0.08, T2 = 0.20, and E = 0.154, we get α = (0.20 - 0.154) / (0.20 - 0.08) = 0.383. Therefore, the effective fatigue damage increment in the current time window is ΔD. c =αΔD p +(1-α)ΔD f =0.383×1.6×10 -7 +0.617×2.8×10 -7 =2.34×10 -7 .

[0141] Therefore, within the current time window, if only the predicted fatigue damage increment is used, the damage increment for the current window is 1.6 × 10⁻⁶. -7 After calibration and correction using the finite element method with directional triggering as described in this invention, the effective fatigue damage increment is 2.34 × 10⁻⁶. -7 This avoids the underestimation of fatigue damage in dangerous areas by the prediction model under complex and variable working conditions.

[0142] Step 7, as Figure 7As shown, update the cumulative fatigue damage and remaining life of the critical parts. This is based on the effective fatigue damage increment ΔD determined within the current time window. c This is added to the cumulative fatigue damage value D at the end of the previous time window. k-1 The cumulative fatigue damage value D at the end of the current time window is obtained. k If the cumulative fatigue damage value D of the dangerous area at the end of the previous time window... k-1 =0.362000000, then the cumulative fatigue damage value at the end of the current time window is: D k =D k-1 +ΔD c =0.362000000+2.34×10 -7 =0.362000234.

[0143] Based on the SN curve method and Miner's linear cumulative damage theory, and combined with a preset damage threshold and the damage growth rate within the most recent m consecutive time windows, the remaining life (RUL) of the hazardous part is updated. When the updated cumulative fatigue damage reaches the preset damage threshold, or the remaining life is lower than the preset safety threshold, the life risk status of the hazardous part is output for condition warning and maintenance scheduling.

[0144] In one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as follows: Figure 3 As shown, the computer device includes a processor, memory, and a network interface connected via a system bus. The processor provides computing and control capabilities. The memory includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The database stores static and dynamic information data. The network interface communicates with external terminals via a network connection. When the computer program is executed by the processor, it implements the steps in the above method embodiments.

[0145] Those skilled in the art will understand that Figure 3 The structure shown is merely a block diagram of a portion of the structure related to the present invention and does not constitute a limitation on the computer device to which the present invention is applied. A specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0146] In addition, the present invention also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps in the above method embodiments.

[0147] In addition, the present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps in the above method embodiments.

[0148] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the methods described above. Any references to memory, storage, databases, or other media used in the embodiments provided by this invention can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, or optical storage, etc. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0149] This invention is not limited to the structures described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of this invention is limited only by the appended claims.

Claims

1. A method for extending the fatigue life of critical parts of a rotor, characterized in that, The method includes: A finite element reference model of the runner is established based on its characteristic parameters, and the dangerous parts of the runner are identified. The real-time monitoring response value and the predicted response value of the dangerous parts of the turbine runner within the current time window are obtained respectively, and a deviation index is constructed based on the difference between the real-time monitoring response value and the predicted response value. The deviation index is compared with a preset threshold, and the comparison result is used to determine whether the dangerous part of the runner triggers finite element orientation calibration. If the dangerous part of the runner triggers finite element orientation calibration, obtain the calibration response result of the dangerous part of the runner within the current time window and calculate the calibration fatigue damage increment. Combine the calibration fatigue damage increment with the pre-calculated predicted fatigue damage increment to determine the effective fatigue damage increment of the current time window. If the dangerous part of the runner does not trigger the finite element orientation calibration, the pre-calculated predicted fatigue damage increment will be used as the effective fatigue damage increment for the current time window. The effective fatigue damage increment is added to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the dangerous parts of the turbine runner. Based on the updated cumulative fatigue damage value, the remaining life and life risk status of the dangerous parts of the turbine runner after the update are calculated.

2. The method for updating the fatigue life of critical parts of a turbine runner according to claim 1, characterized in that, The process of establishing a finite element reference model of the runner based on its characteristic parameters and determining the dangerous parts of the runner includes: Based on the structural characteristic parameters, material parameters and installation constraints of the runner, a finite element reference model of the runner is established; Based on the structural stress analysis results, historical failure locations, and fatigue sensitivity assessment results, the transition zone between the blade root and the upper crown in the finite element reference model of the runner was identified as the critical part of the runner. The baseline response of the critical parts of the runner under various typical operating conditions was calculated based on the finite element baseline model.

3. The method for updating the fatigue life of critical parts of a turbine runner according to claim 2, characterized in that, The reference response includes stress response, strain response, displacement response, and vibration response.

4. The method for updating the fatigue life of critical parts of a turbine runner according to claim 3, characterized in that, The structural feature parameters include the runner structure dimensions, blade geometric parameters, crown structure, and ring structure; The material parameters include elastic modulus, Poisson's ratio, density, and stress-life curve parameters; The typical operating conditions include rated operating conditions, partial load operating conditions, start-stop transition operating conditions, load fluctuation operating conditions, and flow field pulsation operating conditions.

5. The method for updating the fatigue life of critical parts of a turbine runner according to claim 1, characterized in that, The process of acquiring the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and constructing a deviation index based on the difference between the real-time monitoring response value and the predicted response value, includes: Obtain the real-time monitoring response value of the dangerous parts of the turbine within the current time window; The current typical operating condition parameters and historical response sequences are used as inputs to the preset response prediction model. The predicted response values ​​of the dangerous parts of the turbine runner within the current time window are output, and the predicted fatigue damage increment is obtained through fatigue damage calculation. The deviation index for the current time window is calculated based on the difference between the real-time monitoring response value and the predicted response value of the dangerous parts of the runner within the current time window.

6. The method for updating the fatigue life of critical parts of a turbine runner according to claim 5, characterized in that, The formula for calculating the deviation index of the current time window is: ； In the formula, i is the sampling point number within the current time window. This represents the real-time monitoring response value of the dangerous part of the rotor at the i-th sampling point within the current time window. This is the predicted response value corresponding to the i-th sampling point within the current time window. This represents the number of sampling points within the current time window. To prevent tiny positive numbers with a denominator of zero, This is the deviation indicator for the current time window.

7. The method for updating the fatigue life of critical parts of a turbine runner according to claim 1, characterized in that, The step of comparing the deviation index with a preset threshold and determining whether the dangerous part of the runner triggers finite element orientation calibration based on the comparison result includes: The deviation index of the current time window is compared with the preset first calibration trigger threshold and the second forced calibration threshold, respectively; When the deviation index is less than or equal to the first calibration trigger threshold, the finite element orientation calibration of the dangerous part of the runner will not be triggered. When the deviation index is greater than the first calibration trigger threshold and less than or equal to the second mandatory calibration threshold, finite element orientation calibration is performed on the dangerous parts of the runner. When the deviation index exceeds the second mandatory calibration threshold, finite element orientation calibration is performed on the dangerous parts of the runner.

8. The method for updating the fatigue life of critical parts of a turbine runner according to claim 1, characterized in that, If the critical part of the runner triggers finite element orientation calibration, obtaining the calibration response result of the critical part of the runner within the current time window and calculating the calibration fatigue damage increment includes: Based on the finite element reference model of the turbine, which includes the dangerous part of the turbine and the surrounding transition structure, a local sub-model of the dangerous part is constructed. The geometric parameters, material parameters and initial boundary conditions of the local sub-model of the dangerous part are derived from the finite element reference model of the turbine, and the mesh of the local sub-model of the dangerous part is refined in the dangerous part of the turbine and its adjacent area. Based on the global response results of the finite element reference model of the runner under the corresponding operating conditions, the boundary transfer conditions of the local sub-model of the dangerous part are determined by at least one of the following methods: displacement constraint at the truncation boundary, equivalent nodal force or interpolated load. Based on the real-time operating conditions within the current time window, update the local load spectrum parameters, local boundary conditions, and hydraulic pulsation loads acting on the local sub-model of the dangerous area. Based on the updated local load spectrum parameters, local boundary conditions, and hydraulic pulsation loads, the local sub-model of the dangerous part is solved again by finite element method, and the calibration response of the dangerous part of the runner within the current time window is output. The calibration response includes calibration stress response and calibration strain response. The calibration response is sequentially subjected to equivalent stress conversion, rainflow counting, stress-life curve conversion, and linear cumulative damage calculation to obtain the calibration fatigue damage increment.

9. The method for updating the fatigue life of critical parts of a turbine runner according to claim 1, characterized in that, The process of adding the effective fatigue damage increment to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the critical parts of the runner, and calculating the updated remaining life and life risk status of the critical parts of the runner based on the updated cumulative fatigue damage value, includes: The effective fatigue damage increment within the current time window is added to the cumulative fatigue damage value at the end of the previous time window to obtain and update the cumulative fatigue damage value at the end of the current time window. Based on the stress-life curve method and linear cumulative damage theory, the remaining life of the critical parts of the turbine after the update is calculated based on the updated cumulative fatigue damage value; When the updated cumulative fatigue damage value reaches the preset damage threshold, or the updated remaining lifespan is lower than the preset safety threshold, the life risk status of the dangerous parts of the rotor is output.

10. The method for updating the fatigue life of critical parts of a turbine runner according to claim 9, characterized in that, The formula for calculating the remaining lifespan of the critical parts of the turbine after replacement is as follows: RUL=(D lim d k ) / v D ; In the formula, D lim As a preset damage threshold, v D D represents the average damage growth rate within a preset number of continuous time windows. k RUL is the cumulative fatigue damage value updated at the end of the current time window, and RUL is the remaining life of the critical parts of the turbine after the update.

11. A fatigue life renewal system for critical parts of a rotor, characterized in that, The system includes: The benchmark modeling and positioning module is used to establish a finite element benchmark model of the runner based on its characteristic parameters and to determine the dangerous parts of the runner. The response deviation assessment module is used to obtain the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and to construct a deviation index based on the difference between the real-time monitoring response value and the predicted response value. The calibration trigger determination module is used to compare the deviation index with the preset threshold and determine whether the dangerous part of the runner triggers finite element orientation calibration based on the comparison result; The fusion damage calculation module is used to obtain the calibration response results of the dangerous part of the runner within the current time window and calculate the calibration fatigue damage increment if the dangerous part of the runner triggers finite element orientation calibration. The effective fatigue damage increment of the current time window is determined by combining the calibration fatigue damage increment with the pre-calculated predicted fatigue damage increment. The predicted damage carry-over module is used to use the pre-calculated predicted fatigue damage increment as the effective fatigue damage increment for the current time window if the dangerous part of the runner does not trigger the finite element orientation calibration. The life risk status update module is used to add the effective fatigue damage increment to the historical cumulative fatigue damage to obtain the updated cumulative fatigue damage value of the dangerous parts of the turbine. Based on the updated cumulative fatigue damage value, the remaining life and life risk status of the dangerous parts of the turbine are calculated.

12. The fatigue life renewal system for critical parts of a turbine runner according to claim 11, characterized in that, The process of establishing a finite element reference model of the runner based on its characteristic parameters and determining the dangerous parts of the runner includes: Based on the structural characteristic parameters, material parameters and installation constraints of the runner, a finite element reference model of the runner is established; Based on the structural stress analysis results, historical failure locations, and fatigue sensitivity assessment results, the transition zone connecting the blade root and the upper crown in the finite element reference model of the runner was identified as the critical part of the runner. The baseline response of the critical parts of the runner under various typical operating conditions was calculated based on the finite element baseline model.

13. The fatigue life renewal system for critical parts of a turbine runner according to claim 12, characterized in that, The reference response includes stress response, strain response, displacement response, and vibration response.

14. The fatigue life renewal system for critical parts of a turbine runner according to claim 13, characterized in that, The structural feature parameters include the runner structure dimensions, blade geometric parameters, crown structure, and ring structure; The material parameters include elastic modulus, Poisson's ratio, density, and stress-life curve parameters; The typical operating conditions include rated operating conditions, partial load operating conditions, start-stop transition operating conditions, load fluctuation operating conditions, and flow field pulsation operating conditions.

15. The fatigue life renewal system for critical parts of a turbine runner according to claim 11, characterized in that, The process of acquiring the real-time monitoring response value and the predicted response value of the dangerous part of the runner within the current time window, and constructing a deviation index based on the difference between the real-time monitoring response value and the predicted response value, includes: Obtain the real-time monitoring response value of the dangerous parts of the turbine within the current time window; The current typical operating condition parameters and historical response sequences are used as inputs to the preset response prediction model. The predicted response values ​​of the dangerous parts of the turbine runner within the current time window are output, and the predicted fatigue damage increment is obtained through fatigue damage calculation. The deviation index for the current time window is calculated based on the difference between the real-time monitoring response value and the predicted response value of the dangerous parts of the runner within the current time window.

16. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 10.

Images