A Method and System for Identifying Dynamic Loads of Current-Carrying Pipelines by Fusing Data and Models

Abstract

The present invention belongs to the technical field of load identification, and discloses a method and system for identifying dynamic loads of a current-carrying pipeline by fusing data and models. The present invention constructs an initial dynamic model of the current-carrying pipeline; obtains the frequency response function of a local position of the current-carrying pipeline in the pump source off state through testing; uses the frequency response function to correct the initial dynamic model to obtain a corrected dynamic model; obtains the translational vibration response of a local position of the current-carrying pipeline in the pump source on state through testing; establishes a dynamic load identification model of the current-carrying pipeline based on the corrected dynamic model and the translational vibration response; solves the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline. The present invention can improve the accuracy of dynamic load identification of the current-carrying pipeline and reduce costs.

Description

Technical Field

[0001] The present invention belongs to the technical field of load identification, and more specifically, relates to a method and system for identifying dynamic loads of current-carrying pipelines by fusing data and models. Background Art

[0002] In a pipeline system, vibrations caused by internal and external excitations will affect its normal operation and safety. Due to the complex structure and harsh environmental conditions of the pipeline system, studying its load distribution and identification methods can not only help improve the design of the pipeline system, ensure its reliability and safety during use, but also be used for fault diagnosis, vibration reduction and noise reduction, and optimization design.

[0003] For the identification of distributed dynamic loads of current-carrying pipelines, the following main problems currently exist: (1) To identify the excitation load of a current-carrying pipeline, it is first necessary to establish a dynamic model of the pipeline. Currently, in engineering, the dynamic model is often obtained by testing methods. When establishing the dynamic model of a current-carrying pipeline, it is necessary to apply excitation to the pipeline at all measuring points using a force hammer or a shaker, and test the frequency response function of the pipeline. Due to the inevitable influence of measurement errors during cyclic excitation, the established dynamic model of the pipeline has problems such as poor orthogonality and local singularity, resulting in low accuracy of the load identified using this dynamic model. (2) Generally speaking, for the identification of distributed pipeline loads, the dynamic characteristics at each position of the pipeline need to be known. If the traditional direct test method is used to identify the distributed load of the pipeline, it is necessary to perform intensive grid division and large-scale sensor arrangement on the test pipeline. This method has relatively high time and economic costs and is difficult to implement in actual engineering. In summary, the existing solutions have problems of low identification accuracy and high costs. How to solve this problem to meet the actual requirements of dynamic load identification of current-carrying pipelines is a topic of concern and research in this field. Summary of the Invention

[0004] The present invention provides a method and system for identifying dynamic loads of current-carrying pipelines by fusing data and models, which solves the problems of low accuracy and high cost in the identification of dynamic loads of current-carrying pipelines in the prior art.

[0005] The present invention provides a method for identifying dynamic loads of current-carrying pipelines by fusing data and models, including the following steps:

[0006] Construct an initial dynamic model of the current-carrying pipeline; obtain the frequency response function of a local position of the current-carrying pipeline in the pump source off state through testing; use the frequency response function to correct the initial dynamic model to obtain a corrected dynamic model;

[0007] Obtain the translational vibration response of a local position of the current-carrying pipeline in the pump source on state through testing;

[0008] Based on the corrected kinetic model and the translational and vibrational responses, establish a dynamic load identification model for the current-carrying pipeline;

[0009] Solve the dynamic load identification model for the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline.

[0010] Preferably, use the finite element method to construct the initial kinetic model according to the geometric parameters and material parameters of the current-carrying pipeline.

[0011] Preferably, the initial kinetic model is expressed as follows:

[0012]

[0013] In the formula, is the dynamic stiffness matrix of the i th order of the current-carrying pipeline, is the modal vibration mode vector of the i th order of the current-carrying pipeline, is the natural frequency of the i th order of the current-carrying pipeline, is the mass matrix of the current-carrying pipeline, is the damping matrix of the current-carrying pipeline, is the stiffness matrix of the current-carrying pipeline, j is the imaginary unit.

[0014] Preferably, the establishment of the dynamic load identification model for the current-carrying pipeline includes: establishing a vibration transfer matrix, and establishing an objective function according to the vibration transfer matrix and the translational and vibrational responses;

[0015] Among them, the establishment of the vibration transfer matrix includes the following sub-steps: based on the corrected kinetic model, obtain the dynamic flexibility matrix in the translational direction of the current-carrying pipeline; based on the Chebyshev polynomial basis function, establish the dynamic load basis function vector of the current-carrying pipeline; based on the dynamic flexibility matrix in the translational direction of the current-carrying pipeline and the dynamic load basis function vector, construct the vibration transfer matrix.

[0016] Preferably, based on the corrected kinetic model, obtain the dynamic stiffness matrix of all degrees of freedom of the current-carrying pipeline; use matrix elementary transformation to block the dynamic stiffness matrix in two dimensions of translation and rotation; based on the dynamically blocked stiffness matrix obtained by blocking, calculate the dynamic flexibility matrix in the translational direction of the current-carrying pipeline.

[0017] Preferably, the block representation of the dynamic stiffness matrix in two dimensions of translation and rotation is as follows:

[0018]

[0019] In the formula, is the dynamic stiffness matrix of all degrees of freedom of the current-carrying pipeline, , , , are all dynamic stiffness sub - matrices, Ω is the translational direction identifier, and Φ is the rotational direction identifier;

[0020] The dynamic flexibility matrix in the translational direction of the current - carrying pipeline is:

[0021]

[0022] In the formula, is the dynamic flexibility matrix in the translational direction of the current - carrying pipeline, represents the inverse matrix of

[0023] Preferably, the dynamic load basis function vector is expressed as follows:

[0024]

[0025] In the formula, is the dynamic load basis function vector, corresponds to the s th - order Chebyshev polynomial basis function, representing the s th element in the dynamic load basis function vector, s = 1, 2, ……, n , n is the maximum order of the Chebyshev polynomial basis function participating in the load identification calculation;

[0026] The vibration transfer matrix is expressed as follows:

[0027]

[0028] In the formula, Q is the vibration transfer matrix, is the length of the current - carrying pipeline; is the real part of the vibration response of the x th measurement point under a unit excitation applied at the q point, is the imaginary part of the vibration response of the x th measurement point under a unit excitation applied at the q point, q = 1, 2, ……, m , m is the number of measurement points arranged on the current - carrying pipeline; represents the real part of the vibration response of the q th measurement point when is loaded onto the current - carrying pipeline, q ​​

[0029] Preferably, the objective function is expressed as:

[0030]

[0031] Among them, the translational vibration response is expressed as a response vector to obtain a translational vibration response vector, which is expressed as follows:

[0032]

[0033] In the formula, represents the translational vibration response vector, and are respectively the real part and the imaginary part of the translational vibration response of the k th measurement point, where k = 1, 2,..., m , m is the number of measurement points arranged on the current-carrying pipeline; Q is the vibration transfer matrix; represents the coefficient vector to be solved, is the penalty factor.

[0034] Preferably, solving the dynamic load identification model of the current-carrying pipeline includes:

[0035] Using a genetic algorithm to solve the objective function to obtain the coefficient vector A;

[0036] Using the following formula to obtain the dynamic load applied to the current-carrying pipeline:

[0037]

[0038] In the formula, is the dynamic load applied to the current-carrying pipeline, and are respectively the real part coefficient vector and the imaginary part coefficient vector.

[0039] On the other hand, the present invention provides a current-carrying pipeline dynamic load identification system integrating data and model, including:

[0040] A testing unit for testing and obtaining the frequency response function of a local position of the current-carrying pipeline in the pump source off state, and for testing and obtaining the translational vibration response of a local position of the current-carrying pipeline in the pump source on state;

[0041] An initial dynamic model construction unit for constructing an initial dynamic model of the current-carrying pipeline;

[0042] A model correction unit for correcting the initial dynamic model by using the frequency response function to obtain a corrected dynamic model;

[0043] An identification unit, configured to establish a dynamic load identification model of a current-carrying pipeline based on the corrected dynamic model and the translational and vibrational responses, and to solve the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline;

[0044] The current-carrying pipeline dynamic load identification system with data and model fusion is used to execute the steps in the above-mentioned current-carrying pipeline dynamic load identification method with data and model fusion.

[0045] One or more technical solutions provided in the present invention have at least the following technical effects or advantages:

[0046] For the research object of the dynamic load of the current-carrying pipeline, the present invention proposes an identification scheme with data and model fusion. In terms of data, the present invention obtains the frequency response function of a local position of the current-carrying pipeline in the pump source off state through testing, and obtains the translational and vibrational responses of a local position of the current-carrying pipeline in the pump source on state through testing. In terms of model, it includes constructing an initial dynamic model of the current-carrying pipeline. The data and model fusion includes: using the frequency response function of a local position of the current-carrying pipeline in the pump source off state obtained through testing to correct the initial dynamic model to obtain a corrected dynamic model; based on the corrected dynamic model and the translational and vibrational responses of a local position of the current-carrying pipeline in the pump source on state obtained through testing, establishing a dynamic load identification model of the current-carrying pipeline. Finally, the present invention solves the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline.

[0047] The identification scheme with data and model fusion proposed by the present invention uses measured data to update and correct the initial dynamic model of the current-carrying pipeline, obtaining a high-precision dynamic model of the current-carrying pipeline with all degrees of freedom, thereby improving the identification accuracy. At the same time, relying on the established high-precision dynamic model and measured vibration responses, the present invention establishes a dynamic load identification model of the current-carrying pipeline to describe the dynamic load of the current-carrying pipeline, transforming the dynamic load identification problem into an identification model solving problem. It can obtain the global load of the pipeline through local vibration responses with a small number of sensors arranged, accurately identify the dynamic load applied to the current-carrying pipeline, that is, the present invention does not require a large number of sensors to be arranged, effectively reducing the time and economic costs. In summary, the present invention can improve the accuracy of dynamic load identification of the current-carrying pipeline and reduce costs. Description of the Drawings

[0048] Figure 1 It is a flowchart of a method for identifying the dynamic load of a current-carrying pipeline with data and model fusion provided in Embodiment 1 of the present invention. Detailed Embodiment

[0049] In order to better understand the above technical solutions, the above technical solutions will be described in detail below in conjunction with the accompanying drawings of the specification and specific embodiments.

[0050] Example 1:

[0051] Example 1 provides a method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model. Refer to Figure 1 , which includes the following steps:

[0052] Construct an initial dynamic model of the current-carrying pipeline; obtain the frequency response function of a local position of the current-carrying pipeline in the pump source closed state through testing; use the frequency response function to correct the initial dynamic model to obtain a corrected dynamic model;

[0053] Obtain the translational vibration response of a local position of the current-carrying pipeline in the pump source open state through testing;

[0054] Based on the corrected dynamic model and the translational vibration response, establish a dynamic load identification model for the current-carrying pipeline;

[0055] Solve the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline.

[0056] The following will separately elaborate on each part included in Example 1 in detail.

[0057] (1) Construct an initial dynamic model of the current-carrying pipeline.

[0058] Using the finite element method, construct the initial dynamic model according to the geometric parameters and material parameters of the current-carrying pipeline. That is, establish a dynamic finite element simulation beam model of the actual pipeline to be load-identified.

[0059] Specifically, the initial dynamic model is expressed as follows:

[0060]

[0061] In the formula, is the dynamic stiffness matrix of the i th order of the current-carrying pipeline, is the modal vibration mode vector of the i th order of the current-carrying pipeline, is the natural frequency of the i th order of the current-carrying pipeline, is the mass matrix of the current-carrying pipeline, is the damping matrix of the current-carrying pipeline, is the stiffness matrix of the current-carrying pipeline, j is the imaginary unit.

[0062] Assume that the established dynamic finite element beam model of the pipeline has a total of N nodes, then is the mass matrix of the pipeline 6( N +1)×6( N +1), is a damping matrix of dimension 6( N +1)×6( N +1); is a stiffness matrix of dimension 6( N +1)×6( N +1); represents the dynamic stiffness matrix of dimension 6( N +1)×6( N +1); represents the i -th 6( N +1)-dimensional modal vibration mode vector.

[0063] (2) Obtain the frequency response function at a local position of the flow-carrying pipeline in the pump source closed state through testing; use the frequency response function to correct the initial dynamic model to obtain a corrected dynamic model.

[0064] That is, use the frequency response function measured when the pump source is closed to update the system parameters of the established numerical simulation model, construct a simulation model consistent with the dynamic characteristics of the test pipeline, and obtain the high-precision dynamic stiffness matrix of the pump source pipeline, expressed as follows:

[0065]

[0066] (3) Obtain the translational vibration response at a local position of the flow-carrying pipeline in the pump source open state through testing.

[0067] Assume that the number of measurement points arranged on the flow-carrying pipeline is m , m <3( N +1), and express the measured translational vibration response of the flow-carrying pipeline as a response vector.

[0068]

[0069] In the formula, represents the translational vibration response vector, represents the k -th measurement point's translational vibration response, k = 1, 2,..., m , m is the number of measurement points arranged on the flow-carrying pipeline.

[0070] Considering that is a complex number and not convenient for later solution, so expand Equation (3) to:

[0071]

[0072] In the formula, and are respectively The real and imaginary parts.

[0073] (4) Based on the corrected dynamic model and the translational and vibrational responses, establish a dynamic load identification model for the current-carrying pipeline.

[0074] Among them, establishing the dynamic load identification model for the current-carrying pipeline includes:

[0075] (4.1) Establish a vibration transfer matrix, and (4.2) establish an objective function according to the vibration transfer matrix and the translational and vibrational responses.

[0076] Specifically, (4.1) Establishing the vibration transfer matrix includes the following sub-steps:

[0077] (4.1.1) Based on the corrected dynamic model, obtain the dynamic flexibility matrix in the translational direction of the current-carrying pipeline.

[0078] (a) Based on the corrected dynamic model, obtain the dynamic stiffness matrix of the current-carrying pipeline with all degrees of freedom.

[0079] (b) Using elementary matrix transformation, block the dynamic stiffness matrix in two dimensions of translation and rotation.

[0080] The block representation of the dynamic stiffness matrix in two dimensions of translation and rotation is as follows:

[0081]

[0082] In the formula, is the dynamic stiffness matrix of the current-carrying pipeline with all degrees of freedom, , , , are all dynamic stiffness block matrices, Ω is the translational direction identifier, and Φ is the rotational direction identifier.

[0083] (c) Based on the obtained dynamic stiffness block matrices by blocking, calculate the dynamic flexibility matrix in the translational direction of the current-carrying pipeline.

[0084] The dynamic flexibility matrix in the translational direction of the current-carrying pipeline is:

[0085]

[0086] In the formula, is the dynamic flexibility matrix in the translational direction of the current-carrying pipeline, represents the inverse matrix of.

[0087] (4.1.2) Based on the Chebyshev polynomial basis function, establish the dynamic load basis function vector of the current-carrying pipeline.

[0088] First, construct the Chebyshev polynomial basis functions

[0089]

[0090] wherein is the s -order Chebyshev polynomial basis function, represents the length of the current-carrying pipeline.

[0091] Then, establish the dynamic load basis function vector, which is expressed as follows:

[0092]

[0093] wherein is the dynamic load basis function vector, corresponds to the s -order Chebyshev polynomial basis function, representing the s th element in the dynamic load basis function vector, s = 1, 2,..., n , n is the maximum order of the Chebyshev polynomial basis function participating in the load identification calculation, generally the number of test vibration response measurement points.

[0094] (4.1.3) Based on the dynamic flexibility matrix in the translational direction of the current-carrying pipeline and the dynamic load basis function vector, construct the vibration transfer matrix.

[0095] The vibration transfer matrix is expressed as follows:

[0096]

[0097] wherein Q is the vibration transfer matrix, is the length of the current-carrying pipeline.

[0098] Q The H appearing in is the same as the Q appearing in is the included in Equation (8).

[0099] is the real part of the vibration response of the x th measurement point under a unit excitation applied at the q point, is the imaginary part of the vibration response of the x th measurement point under a unit excitation applied at the q point, q = 1, 2,..., m , mis the number of measuring points arranged on the current-carrying pipeline; denotes the real part of the vibration response of the q th measuring point when loaded onto the current-carrying pipeline, denotes the imaginary part of the vibration response of the q th measuring point when loaded onto the current-carrying pipeline.

[0100] Load the dynamic model of the current-carrying pipeline with as the excitation, and calculate the vibration responses and of each measurement point, and obtain the matrix ω at frequency Q .

[0101] (4.2) Establish the objective function.

[0102] The objective function is expressed as:

[0103] (10)

[0104] In the formula, denotes the translational vibration response vector obtained from the previous steps, Q is the vibration transfer matrix obtained from the previous steps; denotes the coefficient vector to be solved; is the penalty factor, which can be obtained by calculating according to the L-curve method.

[0105] (5) Solve the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline.

[0106] Specifically, solving the dynamic load identification model of the current-carrying pipeline includes:

[0107] (5.1) Use the genetic algorithm to solve the objective function and solve for the coefficient vector A;

[0108] (5.2) Use the following formula to obtain the dynamic load applied to the current-carrying pipeline:

[0109] (11)

[0110] In the formula, is the dynamic load applied to the current-carrying pipeline, and are the real part coefficient vector and the imaginary part coefficient vector respectively.

[0111] That is, use formula (9) to establish the vibration transfer matrix Q , and introduce Q into formula (10), and combine with the measured vibration response Form an objective function for identifying the dynamic load of the current-carrying pipeline. By solving, the coefficient vector A can be obtained. Combining with Equation (11), the dynamic load applied to the global current-carrying pipeline can be obtained through the vibration responses of a small number of local degrees of freedom.

[0112] In summary, for the research object of the dynamic load of the current-carrying pipeline in Embodiment 1, a load identification method combining data and model is established: the initial dynamic simulation model of the current-carrying pipeline is updated and corrected by using the measured data based on the translational degrees of freedom to obtain a high-precision dynamic model of the current-carrying pipeline with all degrees of freedom. At the same time, a condensed dynamic stiffness matrix considering the coupling effect of translation and rotation is established based on the model condensation theory to realize the load identification of the current-carrying pipeline under the concerned degrees of freedom and avoid the error caused by the coupling of translational and rotational vibration responses to the load identification of the current-carrying pipeline. According to the dynamic stiffness matrix with all degrees of freedom of the dynamic model, a dynamic flexibility matrix in the translational direction of the pipeline is formed to ensure the accuracy and orthogonality of the dynamic model during subsequent load identification calculations, and at the same time ensure that the excitation load applied to the current-carrying pipeline can be accurately identified only relying on the vibration response in the translational direction. Relying on the established high-precision dynamic simulation model and measured vibration responses, a dynamic load identification model of the current-carrying pipeline is established. The Chebyshev orthogonal polynomial is introduced to describe the dynamic load of the current-carrying pipeline, and the dynamic load identification problem is transformed into a problem of solving the Chebyshev polynomial coefficients, so as to obtain the global load of the pipeline through local vibration responses under the condition of arranging a small number of sensors and accurately identify the dynamic load applied to the current-carrying pipeline.

[0113] Embodiment 2:

[0114] Embodiment 2 provides a system for identifying the dynamic load of a current-carrying pipeline by fusing data and model, including:

[0115] A test unit for testing and obtaining the frequency response function of a local position of the current-carrying pipeline in the pump source off state, and for testing and obtaining the translational vibration response of a local position of the current-carrying pipeline in the pump source on state;

[0116] An initial dynamic model construction unit for constructing an initial dynamic model of the current-carrying pipeline;

[0117] A model correction unit for correcting the initial dynamic model by using the frequency response function to obtain a corrected dynamic model;

[0118] An identification unit for establishing a dynamic load identification model of the current-carrying pipeline based on the corrected dynamic model and the translational vibration response, and for solving the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline.

[0119] The system for identifying dynamic loads of current-carrying pipelines by fusing data and model provided in Embodiment 2 is used to perform the steps in the method for identifying dynamic loads of current-carrying pipelines by fusing data and model as described in Embodiment 1.

[0120] Since the functions of the units in the system for identifying dynamic loads of current-carrying pipelines by fusing data and model provided in Embodiment 2 correspond to the steps in the method for identifying dynamic loads of current-carrying pipelines by fusing data and model provided in Embodiment 1, understanding can be made with reference to the description in Embodiment 1 and will not be elaborated here.

[0121] Finally, it should be noted that the above specific embodiments are only used to illustrate the technical solutions of the present invention rather than to limit them. Although the present invention has been described in detail with reference to the examples, those of ordinary skill in the art should understand that the technical solutions of the present invention can be modified or equivalently replaced without departing from the spirit and scope of the technical solutions of the present invention, and they should all be covered within the scope of the claims of the present invention.

Claims

1. A method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model, characterized in that Including the following steps: Construct an initial dynamic model of the current-carrying pipeline; obtain the frequency response function of a local position of the current-carrying pipeline in the pump source off state through testing; use the frequency response function to correct the initial dynamic model to obtain a corrected dynamic model; Obtain the translational vibration response of a local position of the current-carrying pipeline in the pump source on state through testing, where the translational vibration response is the vibration response in the translational direction; Based on the corrected dynamic model and the translational vibration response, establish a dynamic load identification model for the current-carrying pipeline; Solve the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline.

2. The method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model according to claim 1, characterized in that, Using the finite element method, construct the initial dynamic model according to the geometric parameters and material parameters of the current-carrying pipeline.

3. The method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model according to claim 1, wherein The initial dynamic model is expressed as follows: In the formula, is the dynamic stiffness matrix of the i -th order of the current-carrying pipeline, is the modal vibration mode vector of the i -th order of the current-carrying pipeline, is the natural frequency of the i -th order of the current-carrying pipeline, is the mass matrix of the current-carrying pipeline, is the damping matrix of the current-carrying pipeline, is the stiffness matrix of the current-carrying pipeline, j is the imaginary unit.

4. The method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model according to claim 1, wherein The establishment of the dynamic load identification model for the current-carrying pipeline includes: establishing a vibration transfer matrix, and establishing an objective function according to the vibration transfer matrix and the translational vibration response; Among them, the establishment of the vibration transfer matrix includes the following sub-steps: based on the corrected dynamic model, obtain the dynamic flexibility matrix of the current-carrying pipeline in the translational direction; based on the Chebyshev polynomial basis function, establish the dynamic load basis function vector of the current-carrying pipeline; based on the dynamic flexibility matrix of the current-carrying pipeline in the translational direction and the dynamic load basis function vector, construct the vibration transfer matrix.

5. The method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model according to claim 4, characterized in that, Based on the corrected dynamic model, obtain the dynamic stiffness matrix of the current-carrying pipeline with all degrees of freedom; Using elementary matrix transformation, block the dynamic stiffness matrix in two dimensions of translation and rotation; Based on the dynamically stiffened block matrix obtained by blocking, calculate the dynamic flexibility matrix of the current-carrying pipeline in the translational direction.

6. The method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model according to claim 5, characterized in that, The block representation of the dynamic stiffness matrix in two dimensions of translation and rotation is as follows: In the formula, is the dynamic stiffness matrix of the full degrees of freedom of the current-carrying pipeline, , , , are all dynamic stiffness submatrices, Ω is the translation direction identifier, and Φ is the rotation direction identifier; The dynamic flexibility matrix of the current-carrying pipeline in the translational direction is: In the formula, is the dynamic flexibility matrix in the translational direction of the current-carrying pipeline, represents the inverse matrix of 7. The method for identifying the dynamic load of a current-carrying pipeline by data and model fusion according to claim 4, characterized in that The dynamic load basis function vector is expressed as follows: wherein, is the dynamic load basis function vector, corresponding to the s -th order Chebyshev polynomial basis function, representing the s -th element in the dynamic load basis function vector, s = 1, 2, ……, n , n being the maximum order of the Chebyshev polynomial basis functions participating in the load identification calculation; The vibration transfer matrix is expressed as follows: In the formula, Q is the vibration transfer matrix, is the length of the current-carrying pipeline; is x the real part of the vibration response of the q th measurement point under a unit excitation applied at point is x the imaginary part of the vibration response of the q th measurement point under a unit excitation applied at point q = 1, 2, ……, m , m is the number of measurement points arranged on the current-carrying pipeline; represents the real part of the vibration response of the q th measurement point when is applied to the current-carrying pipeline, and represents q the imaginary part of the vibration response of the th measurement point when is applied to the current-carrying pipeline.

8. The method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model according to claim 7, wherein The objective function is expressed as: Among them, expressing the translational vibration response as a response vector to obtain a translational vibration response vector, which is expressed as follows: In the formula, represents the translational vibration response vector, and are respectively the real part and the imaginary part of the translational vibration response at the k th measurement point, where k = 1, 2, ……, m , m being the number of measurement points arranged on the current-carrying pipeline; Q is the vibration transfer matrix; represents the coefficient vector to be solved, being the penalty factor.

9. The method for identifying the dynamic load of a current-carrying pipeline by fusing data and a model according to claim 8, characterized in that The solution of the dynamic load identification model of the current-carrying pipeline includes: Using a genetic algorithm to solve the objective function to solve the coefficient vector A; Using the following formula to obtain the dynamic load applied to the current-carrying pipeline: In the formula, is the dynamic load applied to the current-carrying pipeline, and are the real part coefficient vector and the imaginary part coefficient vector respectively.

10. A dynamic load identification system for current-carrying pipelines integrating data and models, characterized in that, Including: A test unit for testing and obtaining the frequency response function of a local position of the current-carrying pipeline in the pump source off state, and for testing and obtaining the translational vibration response of a local position of the current-carrying pipeline in the pump source on state, where the translational vibration response is the vibration response in the translational direction; An initial dynamic model construction unit for constructing an initial dynamic model of the current-carrying pipeline; A model correction unit for correcting the initial dynamic model using the frequency response function to obtain a corrected dynamic model; An identification unit for establishing a dynamic load identification model for the current-carrying pipeline based on the corrected dynamic model and the translational vibration response, and for solving the dynamic load identification model of the current-carrying pipeline to obtain the dynamic load applied to the current-carrying pipeline; The system for identifying dynamic loads of current-carrying pipelines by fusing data and models is used to perform the steps in the method for identifying dynamic loads of current-carrying pipelines by fusing data and models according to any one of claims 1-9.

Images