Power station pipeline vibration transmission path analysis method and system based on finite element simulation

Abstract

The invention discloses a power station pipeline vibration transmission path analysis method and system based on finite element simulation. The method comprises the following steps: constructing a pipeline finite element model; based on a pipeline finite element model, determining excitation source positions and target point positions of pipeline vibration, and calculating to obtain frequency response functions from each excitation source to a target point; based on the frequency response function from each excitation source to the target point, carrying out load identification on the load borne by the target point; and on the basis of the identified load, pipeline key transmission path analysis is carried out, the total vibration response of the target point # imgabs0 # and the vibration contribution value and proportion of each excitation source point to the target point # imgabs1 # are obtained, namely, the reason for the vibration problem at the target point position is found, and then pipeline vibration is processed. The system comprises a model construction module, a calculation module, an identification module and an analysis processing module. According to the invention, system evaluation of pipeline vibration is realized.

Description

Technical Field

[0001] The present invention belongs to the technical field of pipeline vibration, and in particular relates to a method and system for analyzing the vibration transmission path of a power station pipeline based on finite element simulation. Background Art

[0002] The vibration problem of power plant pipelines is a very complex problem involving many factors. The force that causes vibration is called the exciting force. According to the source of the exciting force, it can be summarized into several types such as mechanical vibration, fluid vibration, valve self-excited vibration, earthquake, etc. Among them, mechanical vibration includes direct and indirect. Direct mechanical vibration is manifested as the vibration of the equipment directly transmitted to the pipeline through the nozzle; indirect mechanical vibration is manifested as the vibration of the structure or foundation caused by the vibration of the equipment, and then transmitted to the pipeline through the support. For pipelines with smooth medium flow, the main type of pipeline vibration is mechanical vibration. All pipeline vibration problems start from the excitation source and then transmit to a perceptible response position through different transmission paths such as structures. For some important pipelines in power plants, the vibration phenomenon at a specific location is often caused by one or more distant excitation sources. Each excitation source transmits the energy to the response position where the vibration problem occurs after amplifying or attenuating it through different paths. In order to better diagnose and optimize vibration problems, it is necessary to comprehensively consider the conditions of each vibration source and transmission path. The transfer path analysis method (TPA) is an effective solution. The transfer path analysis method (TPA) is a system-level solution based on excitation source-path-response. The TPA method solves specific vibration problems by analyzing the impact of the excitation source and the transfer path on the response and adjusting one or several of these factors. The goal of the transfer path analysis is to calculate the contribution of each path from the source to the response, identify the vibration transfer characteristics of each component in the transfer path, and solve specific problems by adjusting them. At present, there is no systematic vibration analysis of this type of pipeline vibration from each excitation source, and it is impossible to comprehensively evaluate the impact of excitation sources on pipeline vibration in various places. Summary of the invention

[0003] The purpose of the present invention is to provide a method and system for analyzing the vibration transfer path of a power station pipeline based on finite element simulation, so as to realize a systematic evaluation of the pipeline vibration.

[0004] In order to achieve the above object, the present invention adopts the following technical solution: The power station pipeline vibration transfer path analysis method based on finite element simulation includes: Step 1: Construct a pipeline finite element model; Step 2: Based on the pipeline finite element model, determine the excitation source position and target point position of the pipeline vibration, and calculate the frequency response function from each excitation source to the target point; Step 3: Based on the frequency response function from each excitation source to the target point, load identification is performed on the load received by the target point; Step 4: Based on the identified load, perform pipeline critical transfer path analysis to obtain the target point The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed.

[0005] A further improvement of the present invention is that step 1, constructing a pipeline finite element model, comprises: First, three-dimensional modeling is performed according to the design parameters of the pipeline, and the three-dimensional model is pre-processed. Then, the pipeline finite element model is obtained through the processed three-dimensional model. Finally, the pipeline is modally corrected and parameter calculated through finite element modal analysis and experimental modal analysis. The modal correction is to determine the real boundary conditions of the pipeline, and the parameter calculation is to obtain the mass damping coefficient of the pipeline system. and stiffness damping coefficient .

[0006] The present invention is further improved in that, in step 2, based on the pipeline finite element model, the excitation source position and the target point position of the pipeline vibration are determined, and the frequency response function from each excitation source to the target point is calculated, including: Through the pipeline finite element model obtained in step 1, the excitation source of pipeline vibration is analyzed to determine the excitation source point causing pipeline vibration. , and determine the target point of pipeline vibration , by performing harmonic response analysis on the model to obtain the excitation source point To the destination The frequency response function , then the excitation source point i To target point T Contribution of vibration Obtained by formula (1); (1) Where: ——Source of motivation i To target point T The contribution of vibration; ——Source of motivation i Frequency response function to the target point; ——path load on the excitation source point i; When the excitation source has Hours, target point The total amount of vibration is obtained by the following formula: (2).

[0007] A further improvement of the present invention is that step 3, based on the frequency response function from each excitation source to the target point, load identification is performed on the load received by the target point, including: The measured response signal is used in combination with the frequency response function matrix of the system to reversely solve the load, that is, to identify the load; Assume the measured response is X , the frequency response function is H , limited by the measurement method, X and H There will be errors. To control and measure the errors, the overall absolute error is e , which is obtained by the following formula (3) Where: F ——The external load on the system is a matrix composed of ; To make the error e As small as possible, introduce penalty function: (4) In formula (4), when J right F When the first derivative of is zero, the error e If there is a minimum value, then we get F ; (5) Where: I ——the identity matrix; ——Regularization parameter.

[0008] The regularization method is used to solve the inverse of the load. The specific formula is as follows: (6) Where: ——Euclidean norm; X is Measuring point response, is the rate response function, - number of system response points; ; ——Diagonal matrix, its diagonal terms are obtained by the following formula: (7) Where: --matrix The diagonal terms of ; --matrix The diagonal terms of ; When formula (7) obtains the minimum value , which is the regularization parameter determined by the OCV method.

[0009] A further improvement of the present invention is that, in step 4, based on the identified load, a critical transfer path analysis of the pipeline is performed to obtain the target point The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed, including: The frequency response function of each excitation point to the target point is calculated by step 2 , obtained by load identification in step 3 Finally, the target point is obtained through equations (1) and (2) The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point are obtained, that is, the cause of the vibration problem at the target point is found; thus, the pipeline vibration is processed.

[0010] The power station pipeline vibration transfer path analysis system based on finite element simulation includes: Model building module, building pipeline finite element model; The calculation module determines the excitation source position and target point position of the pipeline vibration based on the pipeline finite element model, and calculates the frequency response function from each excitation source to the target point; An identification module performs load identification on the target point based on the frequency response function from each excitation source to the target point; The analysis and processing module performs pipeline critical transfer path analysis based on the identified loads to obtain the target points The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed.

[0011] A further improvement of the present invention is that, in the model building module, a pipeline finite element model is built, including: First, three-dimensional modeling is performed according to the design parameters of the pipeline, and the three-dimensional model is pre-processed. Then, the pipeline finite element model is obtained through the processed three-dimensional model. Finally, the pipeline is modally corrected and parameter calculated through finite element modal analysis and experimental modal analysis. The modal correction is to determine the real boundary conditions of the pipeline, and the parameter calculation is to obtain the mass damping coefficient of the pipeline system. and stiffness damping coefficient .

[0012] A further improvement of the present invention is that, in the calculation module, based on the pipeline finite element model, the excitation source position and the target point position of the pipeline vibration are determined, and the frequency response function from each excitation source to the target point is calculated, including: The pipeline finite element model obtained by the model building module is used to analyze the excitation source of pipeline vibration and determine the excitation source point causing pipeline vibration. , and determine the target point of pipeline vibration , by performing harmonic response analysis on the model to obtain the excitation source point To the destination The frequency response function , then the excitation source point i To target point T Contribution of vibration Obtained by formula (1); (1) Where: ——Source of motivation i To target point T The contribution of vibration; ——Source of motivation i Frequency response function to the target point; ——path load on the excitation source point i; When the excitation source has Hours, target point The total amount of vibration is obtained by the following formula: (2).

[0013] A further improvement of the present invention is that, in the identification module, based on the frequency response function from each excitation source to the target point, load identification is performed on the load received by the target point, including: The measured response signal is used in combination with the frequency response function matrix of the system to reversely solve the load, that is, to identify the load; Assume the measured response is X , the frequency response function is H , limited by the measurement method, X and H There will be errors. To control and measure the errors, the overall absolute error is e , which is obtained by the following formula (3) Where: F ——The external load on the system is a matrix composed of ; To make the error e As small as possible, introduce penalty function: (4) In formula (4), when J right F When the first derivative of is zero, the error e If there is a minimum value, then we getF ; (5) Where: I ——the identity matrix; ——Regularization parameter.

[0014] The regularization method is used to solve the inverse of the load. The specific formula is as follows: (6) Where: ——Euclidean norm; X is Measuring point response, is the rate response function, - number of system response points; ; ——Diagonal matrix, its diagonal terms are obtained by the following formula: (7) Where: --matrix The diagonal terms of ; --matrix The diagonal terms of ; When formula (7) obtains the minimum value , which is the regularization parameter determined by the OCV method.

[0015] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of a method for analyzing vibration transfer paths of power plant pipelines based on finite element simulation.

[0016] Compared with the prior art, the present invention has at least the following beneficial technical effects: The power station pipeline vibration transfer path analysis method and system based on finite element simulation provided by the present invention can systematically perform vibration analysis on pipeline vibration problems from various excitation sources, and comprehensively conduct vibration assessment on the influence of excitation sources in various places on pipeline vibration. Compared with the frequency response function method obtained by traditional test measurement methods, obtaining the frequency response function through finite element simulation can greatly save working time and cost, especially for the acquisition of frequency response functions of complex structures, the traditional measurement method is very time-consuming and labor-intensive. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 It is a flow chart of the present invention.

[0018] Figure 2 Schematic diagram of vibration transmission.

[0019] Figure 3 This is the finite element model diagram.

[0020] Figure 4 It is the frequency amplitude diagram of the signal transmitted from the excitation source point 1 to the target point through path 1 and the measured signal at the target point.

[0021] Figure 5 It is the frequency amplitude diagram of the signal transmitted from the excitation source point 2 to the target point through path 2 and the measured signal at the target point.

[0022] Figure 6 It is a structural block diagram of the power station pipeline vibration transfer path analysis system based on finite element simulation of the present invention. DETAILED DESCRIPTION

[0023] In the following, only some exemplary embodiments are briefly described. As those skilled in the art will appreciate, the described embodiments may be modified in various ways without departing from the spirit or scope of the present invention. Therefore, the drawings and descriptions are considered to be exemplary and non-restrictive in nature.

[0024] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0025] It should also be understood that the terms used in the present specification are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the present specification and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include plural forms.

[0026] It should be further understood that the term "and / or" used in the present description and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0027] Various structural schematic diagrams of the embodiments disclosed in the present invention are shown in the accompanying drawings. These figures are not drawn to scale, and some details are magnified and some details may be omitted for the purpose of clear expression. The shapes of various regions and layers shown in the figures and the relative sizes and positional relationships therebetween are only exemplary, and may deviate in practice due to manufacturing tolerances or technical limitations, and those skilled in the art may additionally design regions / layers with different shapes, sizes, and relative positions according to actual needs.

[0028] The embodiments of the present invention are described in detail below with reference to the accompanying drawings.

[0029] Example 1 The present invention provides a method for analyzing the vibration transfer path of a power station pipeline based on finite element simulation, comprising: Step 1: Construct a pipeline finite element model; Step 2: Based on the pipeline finite element model, determine the excitation source position and target point position of the pipeline vibration, and calculate the frequency response function from each excitation source to the target point; Step 3: Based on the frequency response function from each excitation source to the target point, load identification is performed on the load received by the target point; Step 4: Based on the identified load, perform pipeline critical transfer path analysis to obtain the target point The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed.

[0030] Example 2 refer to Figure 1 , Figure 2 The power station pipeline vibration transfer path analysis method based on finite element simulation provided by the present invention comprises the following steps: Step 1: Construct a finite element model of the pipeline. First, perform high-quality 3D modeling based on the design parameters of the pipeline, and perform certain pre-processing on the 3D model. Then, obtain a relatively accurate finite element model of the pipeline through the processed 3D model. Finally, perform modal correction and parameter calculation on the pipeline through finite element modal analysis and experimental modal analysis. Modal correction is to determine the true boundary conditions of the pipeline, and parameter calculation is to obtain the mass damping coefficient of the pipeline system. and stiffness damping coefficient .

[0031] Step 2: Determine the excitation source position and target point position of the pipeline vibration, and calculate the frequency response function from each excitation source to the target point. Through the pipeline finite element model obtained in step 1, analyze the excitation source of pipeline vibration and determine the excitation source point causing pipeline vibration. , and determine the target point of pipeline vibration , by performing harmonic response analysis on the model to obtain the excitation source point To the destination The frequency response function , then the excitation source point i To target point T Contribution of vibration Obtained by formula (1).

[0032] (1) Where: ——Source of motivation i To target pointT The contribution of vibration, usually an acceleration signal; ——Source of motivation i Frequency response function to the target point; ——The path load on the excitation source point i, which is usually the load force in vibration analysis.

[0033] When the excitation source has Hours, target point The total amount of vibration can be obtained by the following formula: (2) Step 3: Load identification of the target point. According to the TPA theoretical calculation formula (1), the external load It is an important factor affecting the accuracy of TPA analysis. It cannot be obtained directly through measurement, but the measurement of the structural response caused by the load is relatively easy and accurate. Therefore, the measured response signal can be used in combination with the frequency response function matrix of the system to reversely solve the load, that is, to identify the load.

[0034] Assume the measured response is X , the frequency response function is H , limited by the measurement method, X and H There will be errors. To control and measure the errors, the overall absolute error is e , which can be obtained by the following formula (3) Where: F ——The external load on the system is a matrix composed of .

[0035] To make the error e As small as possible, introduce penalty function: (4) In formula (4), when J right F When the first derivative of is zero, the error e If there is a minimum value, we can get F (5) Where: I ——the identity matrix; ——Regularization parameter.

[0036] Load identification is an inverse problem. Load identification in actual engineering is often ill-posed. The ill-conditioned system frequency response function matrix will seriously amplify the smaller measurement error in the response signal during the direct inversion process, causing the inverse load to deviate seriously from the actual load, making load identification meaningless. In order to solve this problem caused by the ill-conditioned system, a regularization method (ordinary cross validation method (OCV)) is used to solve the inverse load. The specific formula is as follows: (6) Where: ——Euclidean norm; X is Measuring point response, is the rate response function, - number of system response points; ; ——Diagonal matrix, its diagonal terms are obtained by the following formula (7) Where: --matrix The diagonal terms of ; --matrix The diagonal terms of .

[0037] When formula (7) obtains the minimum value , which is the regularization parameter determined by the OCV method.

[0038] Step 4: Perform pipeline critical transfer path analysis and obtain the frequency response function of each excitation point to the target point through calculation in step 2. , obtained by load identification in step 3 Finally, the target point can be obtained through equations (1) and (2) The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the main cause of the vibration problem at the target point can be found. Therefore, the pipeline vibration can be treated by corresponding means.

[0039] Example 3 The research object is a cantilever shell structure with a length of 500 mm, a thickness of 5 mm, an arc length of 290 mm, a diameter of 40 mm, and a material of 45 steel with an elastic modulus of N / m 2 , Poisson's ratio , density 7890kg / m 3 , the structural damping is unknown. Figure 3It is a finite element model of a cantilever shell structure, and the mesh is divided by shell elements, with a shell element thickness of 5 mm. Two exciters are arranged on the cantilever shell. The two exciters serve as excitation sources and two vibration transmission paths. The contact position between the exciter and the shell is the two path points. The two excitation sources are excitation source 1 and excitation source 2 respectively. The transmission path from excitation source 1 to the target point is path 1, and the transmission path from excitation source 2 to the target point is path 2. Among the five acceleration sensors arranged on the shell, four are used as calculation reference points for load identification; one is used as the target point for transmission path analysis to collect the response signal of the target point. In this way, a vibration system with two independent excitation sources is formed, which transmits vibration energy to the target position through two paths. Figure 4 It is the frequency amplitude diagram of the signal transmitted from the excitation source point 1 to the target point through path 1 and the measured signal at the target point. Figure 5 It is the frequency amplitude diagram of the signal transmitted from the excitation source point 2 to the target point through path 2 and the measured signal at the target point. It obtains the situation that different paths of the structure are transmitted to the target point in different frequency ranges, that is, the contribution of each path at different frequencies, and finds the cause and source of the specific vibration problem. Therefore, the structure can be optimized and adjusted according to the TPA analysis results to avoid vibration problems at the target point within the specific frequency band of concern. Compared with the frequency response function obtained by traditional experimental measurement methods, obtaining the frequency response function through finite element simulation can greatly save working time and cost, especially for the acquisition of frequency response functions of complex structures. Traditional measurement methods are very time-consuming and labor-intensive.

[0040] Example 4 refer to Figure 6 The power station pipeline vibration transfer path analysis system based on finite element simulation provided by the present invention includes: Model building module, building pipeline finite element model; The calculation module determines the excitation source position and target point position of the pipeline vibration based on the pipeline finite element model, and calculates the frequency response function from each excitation source to the target point; An identification module performs load identification on the target point based on the frequency response function from each excitation source to the target point; The analysis and processing module performs pipeline critical transfer path analysis based on the identified loads to obtain the target points The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed.

[0041] Example 5 The present invention provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the steps of the power plant pipeline vibration transfer path analysis method based on finite element simulation are implemented.

[0042] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0043] The present application is described with reference to the flowcharts and / or block diagrams of the methods, systems and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A system that specifies the functions of a box or boxes.

[0044] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0045] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0046] The above shows and describes the basic principles and main features of the present invention and the advantages of the present invention. It is obvious to those skilled in the art that the present invention is not limited to the details of the above exemplary embodiments, and the present invention can be implemented in other specific forms without departing from the spirit or basic features of the present invention. Therefore, no matter from which point of view, the embodiments should be regarded as exemplary and non-restrictive. The scope of the present invention is defined by the attached claims rather than the above description, and it is intended that all changes falling within the meaning and scope of the equivalent elements of the claims are included in the present invention. Any figure mark in the claims should not be regarded as limiting the claims involved.

[0047] In addition, it should be understood that although this specification is described in accordance with the implementation modes, not every implementation mode contains only one independent technical solution. This description of the specification is only for the sake of clarity. Those skilled in the art should regard the specification as a whole. The technical solutions in each embodiment can also be appropriately combined to form other implementation modes that can be understood by those skilled in the art. The above content is only to illustrate the technical idea of ​​the present invention, and cannot be used to limit the protection scope of the present invention. Any changes made on the basis of the technical solution according to the technical idea proposed by the present invention shall fall within the protection scope of the claims of the present invention.

Claims

1. A power station pipeline vibration transfer path analysis method based on finite element simulation, characterized in that: include: Step 1: Construct a pipeline finite element model; Step 2: Based on the pipeline finite element model, determine the excitation source position and target point position of the pipeline vibration, and calculate the frequency response function from each excitation source to the target point; Step 3: Based on the frequency response function from each excitation source to the target point, load identification is performed on the load received by the target point; Step 4: Based on the identified load, perform pipeline critical transfer path analysis to obtain the target point The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed.

2. The power station pipeline vibration transfer path analysis method based on finite element simulation according to claim 1 is characterized in that: Step 1: Construct a pipeline finite element model, including: First, three-dimensional modeling is performed according to the design parameters of the pipeline, and the three-dimensional model is pre-processed. Then, the pipeline finite element model is obtained through the processed three-dimensional model. Finally, the pipeline is modally corrected and parameter calculated through finite element modal analysis and experimental modal analysis. The modal correction is to determine the real boundary conditions of the pipeline, and the parameter calculation is to obtain the mass damping coefficient of the pipeline system. and stiffness damping coefficient .

3. The power station pipeline vibration transfer path analysis method based on finite element simulation according to claim 2 is characterized in that: Step 2: Based on the pipeline finite element model, determine the excitation source position and target point position of the pipeline vibration, and calculate the frequency response function from each excitation source to the target point, including: Through the pipeline finite element model obtained in step 1, the excitation source of pipeline vibration is analyzed to determine the excitation source point causing pipeline vibration. , and determine the target point of pipeline vibration , by performing harmonic response analysis on the model to obtain the excitation source point To the destination The frequency response function , then the excitation source point i To target point T Contribution of vibration Obtained by formula (1); （1） Where: ——Source of motivation i To target point T The contribution of vibration; ——Source of motivation i Frequency response function to the target point; ——path load on the excitation source point i; When the excitation source has Hours, target point The total amount of vibration is obtained by the following formula: （2）。 4. The power station pipeline vibration transfer path analysis method based on finite element simulation according to claim 3 is characterized in that: Step 3: Based on the frequency response function from each excitation source to the target point, load identification is performed on the load on the target point, including: The measured response signal is used in combination with the frequency response function matrix of the system to reversely solve the load, that is, to identify the load; Assume the measured response is X , the frequency response function is H , limited by the measurement method, X and H There will be errors. To control and measure the errors, the overall absolute error is e , which is obtained by the following formula （3） Where: F ——The external load on the system is a matrix composed of ; To make the error e As small as possible, introduce penalty function: （4） In formula (4), when J right F When the first derivative of is zero, the error e If there is a minimum value, then we get F ; （5） Where: I ——the identity matrix; ——Regularization parameter; The regularization method is used to solve the inverse of the load. The specific formula is as follows: （6） Where: ——Euclidean norm; X is Measuring point response, is the rate response function, - number of system response points; ; ——Diagonal matrix, its diagonal terms are obtained by the following formula: （7） Where: --matrix The diagonal terms of ; --matrix The diagonal terms of ; When formula (7) obtains the minimum value , which is the regularization parameter determined by the OCV method.

5. The power station pipeline vibration transfer path analysis method based on finite element simulation according to claim 4 is characterized in that: Step 4: Based on the identified load, perform pipeline critical transfer path analysis to obtain the target point The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed, including: The frequency response function of each excitation point to the target point is calculated by step 2 , obtained by load identification in step 3 Finally, the target point is obtained through equations (1) and (2) The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point are obtained, that is, the cause of the vibration problem at the target point is found; thus, the pipeline vibration is processed.

6. The power station pipeline vibration transfer path analysis system based on finite element simulation is characterized by: include: Model building module, building pipeline finite element model; The calculation module determines the excitation source position and target point position of the pipeline vibration based on the pipeline finite element model, and calculates the frequency response function from each excitation source to the target point; An identification module performs load identification on the target point based on the frequency response function from each excitation source to the target point; The analysis and processing module performs pipeline critical transfer path analysis based on the identified loads to obtain the target points The total vibration response and each excitation source point to the target point The vibration contribution value and proportion of the target point can be found, that is, the cause of the vibration problem at the target point can be found, and then the pipeline vibration can be processed.

7. The power station pipeline vibration transfer path analysis system based on finite element simulation according to claim 6 is characterized in that: In the model building module, the pipeline finite element model is constructed, including: First, three-dimensional modeling is performed according to the design parameters of the pipeline, and the three-dimensional model is pre-processed. Then, the pipeline finite element model is obtained through the processed three-dimensional model. Finally, the pipeline is modally corrected and parameter calculated through finite element modal analysis and experimental modal analysis. The modal correction is to determine the real boundary conditions of the pipeline, and the parameter calculation is to obtain the mass damping coefficient of the pipeline system. and stiffness damping coefficient .

8. The power station pipeline vibration transfer path analysis system based on finite element simulation according to claim 7 is characterized in that: In the calculation module, based on the pipeline finite element model, the excitation source position and target point position of the pipeline vibration are determined, and the frequency response function from each excitation source to the target point is calculated, including: The pipeline finite element model obtained by the model building module is used to analyze the excitation source of pipeline vibration and determine the excitation source point causing pipeline vibration. , and determine the target point of pipeline vibration , by performing harmonic response analysis on the model to obtain the excitation source point To the destination The frequency response function , then the excitation source point i To target point T Contribution of vibration Obtained by formula (1); （1） Where: ——Source of motivation i To target point T The contribution of vibration; ——Source of motivation i Frequency response function to the target point; ——path load on the excitation source point i; When the excitation source has Hours, target point The total amount of vibration is obtained by the following formula: （2）。 9. The power station pipeline vibration transfer path analysis system based on finite element simulation according to claim 8, characterized in that: In the identification module, based on the frequency response function from each excitation source to the target point, the load on the target point is identified, including: Using the measured response signal and combining it with the system's frequency response function matrix, the load is inversely solved, that is, the load is identified; Assume the measured response is X , the frequency response function is H , limited by the measurement method, X and H There will be errors. To control and measure the errors, the overall absolute error is e , which is obtained by the following formula （3） Where: F ——The external load on the system is a matrix composed of ; To make the error e As small as possible, introduce penalty function: （4） In formula (4), when J right F When the first derivative of is zero, the error e If there is a minimum value, then we get F ; （5） Where: I ——the identity matrix; ——Regularization parameter; The regularization method is used to solve the inverse of the load. The specific formula is as follows: （6） Where: ——Euclidean norm; X is Measuring point response, is the rate response function, - number of system response points; ; ——Diagonal matrix, its diagonal terms are obtained by the following formula: （7） Where: --matrix The diagonal terms of ; --matrix The diagonal terms of ; When formula (7) obtains the minimum value , which is the regularization parameter determined by the OCV method.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the power plant pipeline vibration transfer path analysis method based on finite element simulation described in any one of claims 1 to 5 are implemented.

Images