[go: up one dir, main page]

CN120449534A - Pipeline structure tension identification method based on tension coefficient matrix correction - Google Patents

Pipeline structure tension identification method based on tension coefficient matrix correction

Info

Publication number
CN120449534A
CN120449534A CN202510965609.0A CN202510965609A CN120449534A CN 120449534 A CN120449534 A CN 120449534A CN 202510965609 A CN202510965609 A CN 202510965609A CN 120449534 A CN120449534 A CN 120449534A
Authority
CN
China
Prior art keywords
pipeline
tension
equation
strain
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202510965609.0A
Other languages
Chinese (zh)
Inventor
曾捷
曾诚
周煜靑
马波
张佳溪
李荣冰
姜明顺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
AVIC Chengdu Aircraft Design and Research Institute
Original Assignee
Nanjing University of Aeronautics and Astronautics
AVIC Chengdu Aircraft Design and Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics, AVIC Chengdu Aircraft Design and Research Institute filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202510965609.0A priority Critical patent/CN120449534A/en
Publication of CN120449534A publication Critical patent/CN120449534A/en
Pending legal-status Critical Current

Links

Landscapes

  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

本发明提供了基于张力系数矩阵修正的管道结构张力识别方法，能够辨识管道服役过程中的张力分布，为管道服役状态实时监测与评估提供数据支撑，从而提升管道服役安全性和稳定性。本发明方法包括：构建管道坐标系，根据管道约束端之间相对位置，确定管道张力函数参数，进而构建张力函数模型。在此基础上，推导管道受到集中力作用时的张力识别方程。引入重力约束迭代法，求解管道未受到集中力作用时的张力识别方程。同时，通过提取集成于管道的光纤传感器应变感知信息，构建张力特征修正系数矩阵，用以显著提高张力识别精度。

The present invention provides a pipeline structure tension identification method based on tension coefficient matrix correction. The method can identify the tension distribution during pipeline service and provide data support for real-time monitoring and evaluation of pipeline service status, thereby improving pipeline service safety and stability. The method of the present invention includes: constructing a pipeline coordinate system, determining the pipeline tension function parameters based on the relative positions between the pipeline constraint ends, and then constructing a tension function model. On this basis, the tension identification equation when the pipeline is subjected to a concentrated force is derived. The gravity constraint iteration method is introduced to solve the tension identification equation when the pipeline is not subjected to a concentrated force. At the same time, by extracting the strain sensing information of the optical fiber sensor integrated in the pipeline, a tension feature correction coefficient matrix is constructed to significantly improve the tension identification accuracy.

Description

Pipeline structure tension identification method based on tension coefficient matrix correction
Technical Field
The invention belongs to the technical field of pipeline structure monitoring, and particularly relates to a pipeline structure tension identification method based on tension coefficient matrix correction.
Background
Modern industry, city lifelines, and aeronautical fields (e.g., air-fuelling hoses) are highly dependent on various types of piping systems. These pipe structures often operate under complex load conditions, with pipe tension as a core safety and performance indicator, directly related to their buckling resistance, fatigue life, overall stability and reliability. However, the existing detection technology is mainly aimed at the defect of the pipeline structure, and is difficult to adapt to long-distance, continuous and real-time monitoring requirements. Therefore, the development of a high-precision and intelligent tension recognition method has the indistinct necessity, is a scientific basis for evaluating extreme environments, dynamic service and old pipeline states, optimizing maintenance/operation decisions, prolonging asset life/improving task success rate, and meets the strategic requirements of increasingly strict safety regulations and realizing pipeline integrity lean management. The research has irreplaceable importance for constructing quality safety, improving economic benefit and enhancing national defense efficiency.
In recent years, strain sensing technology based on optical fiber sensors is gradually applied to the field of structural health monitoring. The optical fiber sensor has the advantages of high precision, electromagnetic interference resistance, light weight, corrosion resistance and the like, and can monitor the strain distribution of the structure in real time. However, at present, a mature and universal tension identification method capable of effectively utilizing abundant strain data acquired by an optical fiber sensor is not available, and a mapping relation between distributed strain and tension distribution of a pipeline structure under a complex working condition is difficult to establish.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a pipeline structure tension identification method based on tension coefficient matrix correction, which comprises the following steps:
Step 1, establishing a rectangular coordinate system by taking the lower end of a pipeline as a coordinate origin, determining the length l of the pipeline, determining the height difference h 0 between two constraint ends of the pipeline, connecting an included angle theta between two ends of the pipeline, and carrying out gravity load p on the pipeline;
Step 2, deriving a corresponding tension recognition equation when the pipeline is acted by a vertical downward concentrated force according to a tension function theory;
step 3, carrying out numerical solution on key parameters of a tension identification equation which is not acted by the vertical downward concentrated force on the pipeline by adopting a gravity constraint iteration method;
And 4, sticking an optical fiber sensor on the surface of the pipeline, obtaining pipeline surface strain distribution data in the service process, constructing a tension characteristic correction coefficient matrix according to the measured strain original data matrix and the first derivative matrix, substituting the tension characteristic correction coefficient matrix into a tension identification equation, and correcting tension identification errors generated by the fact that the tension identification equation does not contain detailed material coefficient characteristics of the pipeline.
In step 1, defining the horizontal length of the pipeline AB as l, and solving the gravity load p born by the pipeline through the pipeline density.
The step 2 comprises the following steps:
Step 2.1, taking any micro-segment ds on the pipeline, setting the horizontal component force of the pipeline under the tension T as H and the vertical component force as V, and obtaining a pipeline static equilibrium equation in the service process through statics analysis:
(1),
Wherein d is a differential symbol, Y is the Y-axis coordinate of the pipeline micro-element section, and X is the X-axis coordinate of the pipeline micro-element section;
Carrying out quadratic integral solution on the formula (1), and obtaining a linear equation in a free state of the pipeline according to boundary conditions, wherein the abscissa x 1 =0 and the ordinate y 1 =0 of the starting point of the pipeline, the abscissa x 2 =l and the ordinate y 2 =l of the ending point of the pipeline:
(2),
Wherein the method comprises the steps of 、Any constant generated for the integration process is determined by equation (3):
(3),
wherein arsinh denotes an anti-hyperbolic sine function;
Step 2.2, setting the point C as the lowest point of the pipeline, deriving the formula (3) and taking the extreme point to obtain the abscissa of the lowest point C The method comprises the following steps:
(4),
substituting formula (4) into formula (2) to obtain a pipeline linear function y (x) under the condition of the known lowest point:
(5),
let the tension of the point Q (x, y) at any position of the pipeline be T, the vertical component of the tension be V, and obtain the pipeline tension recognition equation by statics analysis:
(6),
(7),
the pipeline is at the lowest point according to formula (7) The tension value at the position is minimum, from the lowest pointThe tension value gradually increases towards the two ends, and the tension value reaches the maximum at the point B at the higher end, so that the tension value of the pipeline is the maximumAnd minimum valueThe method comprises the following steps:
(8),
Wherein the method comprises the steps of Is the component force of the tension of the point B in the vertical direction;
Step 2.3, when the pipe is subjected to a vertical downward concentrated force, the pipe shape is generalized to two, one of which is applied to the center of the pipe, called pipe line type I, and the other of which is applied near both ends of the pipe, called pipe line type II.
In step 2.3, when the pipeline is in a linear form I, the combined type (2), (6) and (7) are obtained:
(9),
when the pipeline is in a linear II shape, the combined type (2), (6) and (7) are obtained:
(10),
wherein:
(11),
Wherein the method comprises the steps of The X-axis coordinates of the position for force application to be concentrated,In the form of a Y-axis coordinate,Concentrated force loads are applied to the pipeline;
Obtained by numerical analysis of (9) and (10) Will beSubstituting the values (6) and (7) to obtain the tension distribution characteristics of the pipeline when the pipeline is subjected to the vertical downward concentrated force.
The step 3 comprises the following steps:
step 3.1, according to the abscissa of the lowest point of the known pipeline Ordinate ofObtaining parameters ofThe constraint equation of (2) is:
(12),
Step 3.2, adopting a gravity constraint iteration method to carry out numerical solution and constructing an iteration function :
(13),
Step 3.3, solving parameters by adopting a gravity constraint iteration methodWhen calculating the derivative value of the iterative function:
(14),
Wherein:
(15),
And 3.4, the gravity constraint iteration method is as follows:
(16),
Wherein, the Representing the parameters obtained in the nth iteration;Is a convergence factor; Representing parameter values Functions of (a), i.e., functions M and N; the convergence condition is,The convergence accuracy is set;
And 3.5, verifying a calculation result by adopting a gravity constraint equation, and setting the vertical tension component at the pipeline A as V A and the vertical tension component at the pipeline B as V B, wherein the steps are as follows:
(17),
according to the gravity constraint equation ps=v A+VB, the calculation result needs to satisfy:
(18),
Where s represents the pipe arc length.
Step 4 comprises:
Step 4.1, arranging three strain sensing paths of the fiber Bragg grating sensors on the surface of the pipeline along the circumferential direction of 0 degree, 120 degrees and 240 degrees, wherein the actual pulling and pressing strain of the pipeline is acquired by the fiber Bragg grating sensors ;
Step 4.2, the strain obtained by each discrete measuring point is in turnThen define the measured strain matrix E as:
(19),
Wherein the method comprises the steps of Representing the measured strain of the nth discrete measurement point;
calculating and measuring a first derivative matrix of a strain matrix by adopting a central difference method :
(20),
Wherein the method comprises the steps ofThe strain value of the ith optical fiber measuring point is; The strain value is the i-1 th optical fiber measuring point strain value; I is the X-axis coordinate value of the ith optical fiber measuring point, and the value of i is 1-n;
calculating to obtain a first derivative matrix of the measured strain matrix The method comprises the following steps:
(21),
step 4.3, obtaining a tension characteristic correction coefficient matrix K according to the measured strain matrix and the first derivative of the measured strain matrix obtained by the optical fiber sensing system:
(22),
And 4.4, establishing a pipeline tension identification equation based on the tension correction coefficient matrix.
In step 4.1, the actual tension-compression strain of the pipelineThe calculation formula is as follows:
(23),
Wherein, the 、、The strain is measured by the fiber Bragg grating sensors with the angles of 0 degree, 120 degrees and 240 degrees respectively.
In step 4.1, the pipe tension identification equation based on the tension correction coefficient matrix is as follows:
(24),
(25),
wherein Y (x) represents a pipeline linear function based on a tension correction coefficient matrix; Representing a tube structure tension identification function based on a tension correction coefficient matrix.
The invention also provides an electronic device comprising a processor and a memory, the memory storing program code which, when executed by the processor, causes the processor to perform the steps of the method.
The invention also provides a storage medium storing a computer program or instructions which, when run on a computer, perform the steps of the method.
The method constructs a tension function model facing the pipeline by determining the coordinate system of the pipeline and combining the spatial position relation of two constraint ends of the pipeline. Furthermore, a gravity constraint numerical iteration method is introduced, tension distribution of the pipeline under the condition of only being subjected to dead weight is obtained through solving, and the problem of difficulty in solving the equation due to the fact that the tension function is an overrun equation is solved. Meanwhile, strain information of the optical fiber sensor integrated with the pipeline is extracted, a tension characteristic correction coefficient matrix is constructed, and tension identification accuracy is remarkably improved.
The invention aims to provide data support for monitoring and evaluating the service state of the pipeline in real time, so that the service safety and stability of the pipeline are improved.
The invention has the beneficial effects that the optical fiber sensor and the innovative tension recognition model are fused, and the invention brings remarkable improvement to pipeline monitoring. Has the following beneficial effects:
Correlating the pipeline linear function with the load born by the pipeline, deducing a tension identification equation of the pipeline under different stress conditions, and breaking through the limitation that the pipeline tension function only contains pipeline gravity parameters; the method comprises the steps of providing an iteration method for introducing gravity constraint numerical values, solving the tension distribution of the pipeline under the condition of only being subjected to dead weight, solving the problem of difficulty in solving the equation due to the fact that the tension function is an overrun equation, and meanwhile, constructing a tension characteristic correction coefficient matrix by extracting strain information of an optical fiber sensor integrated with the pipeline, so that the precision of tension identification is remarkably improved.
Drawings
The foregoing and/or other advantages of the invention will become more apparent from the following detailed description of the invention when taken in conjunction with the accompanying drawings and detailed description.
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a schematic representation of simulated constraint application made to derive one embodiment of the present invention.
FIG. 3 is a graph comparing tension profiles obtained by the method of the present invention.
FIG. 4 is a graph showing the relative errors in tension obtained by the method of the present invention.
Detailed Description
As shown in fig. 1, the embodiment provides a method for identifying the tension of a pipe structure based on the correction of a tension coefficient matrix, which comprises the following steps:
The method comprises the steps of 1, establishing a rectangular coordinate system by taking the lower end of a pipeline as a coordinate origin, determining the horizontal length l of the pipeline, determining the height difference h 0 between two constraint ends of the pipeline, connecting an included angle theta between two ends of the pipeline, and carrying out gravity load p on the pipeline;
step 2, deriving a tension recognition equation of the pipeline acted by the vertical downward concentrated force according to the tension function theory;
step 3, carrying out numerical solution on key parameters of a tension identification equation, which are not acted by vertical downward concentrated force, of the pipeline by adopting a gravity constraint iteration method;
And 4, sticking an optical fiber sensor on the surface of the pipeline, obtaining the strain distribution data of the surface of the pipeline in the service process, and constructing a tension characteristic correction coefficient matrix according to the measured strain original data matrix and the first derivative matrix thereof. Substituting the tension characteristic correction coefficient matrix into a tension identification equation to correct tension identification errors generated by not containing detailed material coefficient characteristics of the pipeline.
The step 2 comprises the following steps:
And 2.1, taking any micro-segment ds on the pipeline. Assuming that the horizontal component force of the tension T applied to the pipeline is H and the vertical component force is V, obtaining a static equilibrium equation of the pipeline in the process by statics analysis:
(1),
In the formula, d is a differential symbol, Y is the Y-axis coordinate of the pipeline micro-element section, and X is the X-axis coordinate of the pipeline micro-element section.
Carrying out quadratic integral solution on the formula (1), and obtaining a linear equation in a free state of the pipeline according to boundary conditions, wherein the abscissa x 1 =0 and the ordinate y 1 =0 of the starting point of the pipeline, the abscissa x 2 =l and the ordinate y 2 =l of the ending point of the pipeline:
(2),
Wherein the method comprises the steps of 、Any constant generated for the integration process is determined by equation (3):
(3),
wherein arsinh denotes an anti-hyperbolic sine function;
Step 2.2, setting the point C as the lowest point of the pipeline, deriving the formula (3) and taking the extreme point to obtain the abscissa of the lowest point C The method comprises the following steps:
(4),
substituting formula (4) into formula (2) to obtain a pipeline linear function y (x) under the condition of the known lowest point:
(5),
let the tension of the point Q (x, y) at any position of the pipeline be T, the vertical component of the tension be V, and obtain the pipeline tension recognition equation by statics analysis:
(6),
(7),
the pipeline is at the lowest point according to formula (7) The tension value at the position is minimum, from the lowest pointThe tension value gradually increases towards the two ends, and the tension value reaches the maximum at the point B at the higher end, so that the tension value of the pipeline is the maximumAnd minimum valueThe method comprises the following steps:
(8),
Wherein the method comprises the steps of Is the component force of the tension of the point B in the vertical direction;
Step 2.3, when the pipe is subjected to a vertical downward concentrated force, the pipe shape is generalized to two, one of which is applied to the center of the pipe, called pipe line type I, and the other of which is applied near both ends of the pipe, called pipe line type II.
In step 2.3, when the pipeline is in a linear form I, the combined type (2), (6) and (7) are obtained:
(9),
when the pipeline is in a linear II shape, the combined type (2), (6) and (7) are obtained:
(10),
wherein:
(11),
Wherein the method comprises the steps of The X-axis coordinates of the position for force application to be concentrated,In the form of a Y-axis coordinate,Concentrated force loads are applied to the pipeline;
Obtained by numerical analysis of (9) and (10) Will beSubstituting the values (6) and (7) to obtain the tension distribution characteristics of the pipeline when the pipeline is subjected to the vertical downward concentrated force.
The step 3 comprises the following steps:
step 3.1, according to the abscissa of the lowest point of the known pipeline Ordinate ofObtaining parameters ofThe constraint equation of (2) is:
(12),
Step 3.2, adopting a gravity constraint iteration method to carry out numerical solution and constructing an iteration function :
(13),
Step 3.3, solving parameters by adopting a gravity constraint iteration methodWhen calculating the derivative value of the iterative function:
(14),
Wherein:
(15),
And 3.4, the gravity constraint iteration method is as follows:
(16),
Wherein, the Representing the parameters obtained in the nth iteration;Is a convergence factor, is to ensure,Halving may be performed on the basis of the previous iteration.Representing parameter valuesFunctions of (a), i.e., functions M and N; the convergence condition is,The convergence accuracy is set;
And 3.5, verifying a calculation result by adopting a gravity constraint equation, and setting the vertical tension component at the pipeline A as V A and the vertical tension component at the pipeline B as V B, wherein the steps are as follows:
(17),
according to the gravity constraint equation ps=v A+VB, the calculation result needs to satisfy:
(18),
Where s represents the pipe arc length.
Step 4 comprises:
Step 4.1, arranging three strain sensing paths of the fiber Bragg grating sensors on the surface of the pipeline along the circumferential direction of 0 degree, 120 degrees and 240 degrees, wherein the actual pulling and pressing strain of the pipeline is acquired by the fiber Bragg grating sensors ;
Step 4.2, the strain obtained by each discrete measuring point is in turnThen define the measured strain matrix E as:
(19),
Wherein the method comprises the steps of Representing the measured strain of the nth discrete measurement point;
calculating and measuring a first derivative matrix of a strain matrix by adopting a central difference method :
(20),
Wherein the method comprises the steps ofThe strain value of the ith optical fiber measuring point is; The strain value is the i-1 th optical fiber measuring point strain value; I is the X-axis coordinate value of the ith optical fiber measuring point, and the value of i is 1-n;
calculating to obtain a first derivative matrix of the measured strain matrix The method comprises the following steps:
(21),
step 4.3, obtaining a tension characteristic correction coefficient matrix K according to the measured strain matrix and the first derivative of the measured strain matrix obtained by the optical fiber sensing system:
(22),
And 4.4, establishing a pipeline tension identification equation based on the tension correction coefficient matrix.
In step 4.1, the actual tension-compression strain of the pipelineThe calculation formula is as follows:
(23),
Wherein, the 、、The strain is measured by the fiber Bragg grating sensors with the angles of 0 degree, 120 degrees and 240 degrees respectively.
In step 4.1, the pipe tension identification equation based on the tension correction coefficient matrix is as follows:
(24),
(25),
wherein Y (x) represents a pipeline linear function based on a tension correction coefficient matrix; Representing a tube structure tension identification function based on a tension correction coefficient matrix.
In one embodiment of the present invention, a method for identifying tension of a pipe structure based on tension coefficient matrix correction is provided, including:
Firstly, deriving a tension recognition equation of the pipeline under the action of a vertical downward concentrated force according to a tension function theory, and solving the tension recognition problem when the pipeline is under an external acting force;
Secondly, carrying out numerical solution on key parameters of a tension recognition equation, which are not acted by vertical downward concentrated force, of the pipeline by adopting a gravity constraint iteration method, so as to solve the problem of tension recognition when the pipeline is only acted by gravity;
And sticking an optical fiber sensor on the surface of the pipeline, acquiring the strain data of the surface of the pipeline in the process, and constructing a tension characteristic correction coefficient matrix by using the strain original data matrix and the first derivative matrix thereof. Substituting the tension characteristic correction coefficient matrix into a tension identification equation to correct tension identification errors generated by not containing detailed material coefficient characteristics of the pipeline.
And finally, using Ansys Workbench simulation software to establish a pipeline model, applying fixed constraint to two ends of the pipeline, applying self gravity load on the pipeline as shown in figure 2, extracting pipeline simulation tension data and strain data, substituting the strain data into the method to obtain the method identification tension data, and comparing the two tension data as shown in figure 3. The tension data of the method is subjected to relative error on simulation tension data, as shown in fig. 4.
The invention provides a method for identifying the tension of a pipeline structure based on tension coefficient matrix correction, and the method and the way for realizing the technical scheme are numerous, the above description is only a preferred embodiment of the invention, and it should be pointed out that a plurality of improvements and modifications can be made to those skilled in the art without departing from the principle of the invention, and the improvements and the modifications are also considered as the protection scope of the invention. The components not explicitly described in this embodiment can be implemented by using the prior art.

Claims (10)

1.基于张力系数矩阵修正的管道结构张力识别方法，其特征在于，包括以下步骤：1. A pipeline structure tension identification method based on tension coefficient matrix correction, characterized by comprising the following steps: 步骤1，以管道较低端为坐标原点建立直角坐标系，确定管道长度l，管道两个约束端之间高度差h 0，管道两端连线夹角θ，管道所受重力载荷p；Step 1: Establish a rectangular coordinate system with the lower end of the pipeline as the origin, determine the pipeline length l , the height difference h 0 between the two constrained ends of the pipeline, the angle θ between the two ends of the pipeline, and the gravity load p on the pipeline; 步骤2，根据张力函数理论，推导得到管道受到垂直向下的集中力作用时对应的张力识别方程；Step 2: Based on the tension function theory, derive the corresponding tension identification equation when the pipeline is subjected to a vertical downward concentrated force; 步骤3，采用重力约束迭代法对管道未受到垂直向下集中力作用时的张力识别方程所属关键参数进行数值求解；Step 3: Using the gravity constraint iteration method, numerically solve the key parameters of the tension identification equation when the pipeline is not subjected to the vertical downward concentrated force; 步骤4，在管道表面粘贴光纤传感器，获取服役过程中管道表面应变分布数据，根据所测应变原始数据矩阵与一阶导数矩阵，构建张力特征修正系数矩阵；将张力特征修正系数矩阵代入张力识别方程，修正因未包含管道详细材料系数特征而产生的张力识别误差。Step 4: Attach fiber optic sensors to the pipeline surface to obtain the strain distribution data on the pipeline surface during service. Construct a tension characteristic correction coefficient matrix based on the measured strain raw data matrix and the first-order derivative matrix. Substitute the tension characteristic correction coefficient matrix into the tension identification equation to correct the tension identification error caused by not including the detailed material coefficient characteristics of the pipeline. 2.根据权利要求1所述的方法，其特征在于，步骤1中，定义管道AB的水平长度为l，通过管道密度求解管道所受重力载荷p。2. The method according to claim 1, characterized in that in step 1, the horizontal length of the pipeline AB is defined as l , and the gravity load p on the pipeline is solved by the pipeline density. 3.根据权利要求2所述的方法，其特征在于，步骤2包括：3. The method according to claim 2, wherein step 2 comprises: 步骤2.1，取管道上任意微元段ds，设定管道受到张力T的水平方向分力为H，竖直方向分力为V，则由静力学分析得到服役过程中管道静力平衡方程：In step 2.1, take any infinitesimal segment ds on the pipeline and set the horizontal component of the tension T on the pipeline to H and the vertical component to V. Then, the static equilibrium equation of the pipeline during service is obtained from static analysis: （1）， (1), 其中d为微分符号；y为管道微元段Y轴坐标；x为管道微元段X轴坐标；Where d is the differential symbol; y is the Y -axis coordinate of the pipeline element segment; x is the X- axis coordinate of the pipeline element segment; 对式(1)进行二次积分求解，并根据边界条件：管道起始点横坐标x 1=0，纵坐标y 1=0；管道终点横坐标x 2=l，纵坐标y 2=l，得到管道自由状态下的线型方程：Solve equation (1) by quadratic integration, and according to the boundary conditions: the horizontal coordinate x 1 = 0 and the vertical coordinate y 1 = 0 at the starting point of the pipeline; the horizontal coordinate x 2 = 1 and the vertical coordinate y 2 = 1 at the end point of the pipeline, the linear equation of the pipeline in the free state is obtained: （2）， (2), 其中、为积分过程产生的任意常数，由式(3)确定：in 、 is an arbitrary constant generated by the integration process, which is determined by formula (3): （3）， (3), 其中arsinh表示反双曲正弦函数；Where arsinh represents the inverse hyperbolic sine function; 步骤2.2，设C点为管道最低点，对式(3)求导并取极值点，得到最低点C处横坐标为：Step 2.2: Let point C be the lowest point of the pipeline, take the derivative of equation (3) and take the extreme point to obtain the horizontal coordinate of the lowest point C. for: （4）， (4), 将式(4)代入式(2)，得到已知最低点条件下的管道线型函数y(x)：Substituting equation (4) into equation (2), we can obtain the pipeline linear function y ( x ) under the condition of known lowest point: （5）， (5), 设管道任意位置Q(x, y)点张力为T，张力竖直分量为V，由静力学分析得到管道张力识别方程：Assume that the tension at any position Q ( x , y ) in the pipeline is T and the vertical component of the tension is V. The pipeline tension identification equation is obtained from static analysis: （6）， (6), （7）， (7), 根据式(7)，管道在最低点处的张力值最小，由最低点处向两端张力值逐渐增大，在较高端B点达到最大，因此，管道张力最大值和最小值为：According to formula (7), the pipeline is at the lowest point The tension value at the lowest point is The tension value gradually increases towards both ends and reaches the maximum at the higher end point B. Therefore, the maximum tension of the pipeline and minimum value for: （8）， (8), 其中为B点张力竖直方向分力；in is the vertical component of the tension at point B; 步骤2.3，当管道受到垂直向下的集中力作用时，管道形态归纳为两种，其中一种集中力施加在管道中心，称为管道线型Ⅰ，另一种集中力施加在靠近管道两端处，称为管道线型Ⅱ。Step 2.3, when the pipeline is subjected to a vertical downward concentrated force, the pipeline shape can be summarized into two types. One concentrated force is applied at the center of the pipeline, called pipeline line type I, and the other concentrated force is applied near the two ends of the pipeline, called pipeline line type II. 4.根据权利要求3所述的方法，其特征在于，步骤2.3中，当管道形态为线型Ⅰ时，联立式(2)、(6)、(7)得到：4. The method according to claim 3, characterized in that in step 2.3, when the pipeline shape is linear type I, the simultaneous equations (2), (6), and (7) are obtained: （9）， (9), 当管道形态为线型Ⅱ时，联立式(2)、(6)、(7)得到：When the pipeline shape is linear II, the simultaneous equations (2), (6), and (7) are: （10）， (10), 式中：Where: （11）， (11), 其中为集中力施加位置X轴坐标，为Y轴坐标，为管道所受集中力载荷；in is the X- axis coordinate of the concentrated force application position, is the Y- axis coordinate, is the concentrated load on the pipeline; 通过数值求解式(9)、(10)得到的值，将代入式(6)、(7)求得当管道受到垂直向下的集中力作用时，管道张力分布特征。By numerically solving equations (9) and (10), we can obtain The value of Substituting into equations (6) and (7) we can obtain the distribution characteristics of pipeline tension when the pipeline is subjected to a vertical downward concentrated force. 5.根据权利要求4所述的方法，其特征在于，步骤3包括：5. The method according to claim 4, wherein step 3 comprises: 步骤3.1，根据已知管道最低点横坐标，纵坐标，得到关于参数的约束方程为：Step 3.1, based on the known horizontal coordinate of the lowest point of the pipeline , vertical coordinate , get the parameters The constraint equation is: （12）， (12), 步骤3.2，采用重力约束迭代法进行数值求解，构造迭代函数：Step 3.2: Use gravity constraint iteration method to solve numerically and construct iteration function : （13）， (13), 步骤3.3，采用重力约束迭代法求解参数时，计算迭代函数的导数值：Step 3.3, use gravity constraint iteration method to solve the parameters Calculate the derivative of the iterative function when : （14）， (14), 其中：in: （15）， (15), 步骤3.4，重力约束迭代法格式为：Step 3.4, the format of gravity constraint iteration method is: （16）， (16), 其中，表示第n次迭代得到的参数；为收敛因子；表示参数值的函数，即函数M和N；收敛条件为，为设定的收敛精度；in, Indicates the parameters obtained at the nth iteration ; is the convergence factor; Indicates parameter value Functions of M and N ; the convergence condition is , is the set convergence accuracy; 步骤3.5，采用重力约束方程对计算结果进行验证，设管道A处张力竖直分量为V A ，管道B处张力竖直分量为V B ，则：In step 3.5, the calculation results are verified using the gravity constraint equation. Assume that the vertical component of the tension at pipe A is VA and the vertical component of the tension at pipe B is VB . Then: （17）， (17), 根据重力约束方程ps=V A +V B ，计算结果需满足：According to the gravity constraint equation ps = VA + VB , the calculation results must satisfy: （18）， (18), 其中s表示管道弧长。Where s represents the arc length of the pipe. 6.根据权利要求5所述的方法，其特征在于，步骤4包括：6. The method according to claim 5, wherein step 4 comprises: 步骤4.1，在管道表面沿周向0°、120°、240°布置三条光纤Bragg光栅传感器应变感知路径；由光纤Bragg光栅传感器采集的管道实际拉压应变；Step 4.1: Arrange three fiber Bragg grating sensor strain sensing paths along the circumferential direction of 0°, 120°, and 240° on the pipeline surface; the actual tensile and compressive strains of the pipeline collected by the fiber Bragg grating sensor are ; 步骤4.2，各离散测点所得应变依次为，则定义测量应变矩阵E为：In step 4.2, the strains obtained at each discrete measuring point are: , then the measurement strain matrix E is defined as: （19）， (19), 其中表示第n个离散测点的测量所得应变；in represents the strain measured at the nth discrete measuring point; 采用中心差分法计算测量应变矩阵一阶导数矩阵：The first-order derivative matrix of the measured strain matrix is calculated using the central difference method : （20）， (20), 其中为第i个光纤测点应变值；为第i-1个光纤测点应变值；为第i个光纤测点X轴坐标值；i取值1~n；in is the strain value of the i -th optical fiber measuring point; is the strain value of the i -1th optical fiber measuring point; is the X -axis coordinate value of the i -th optical fiber measurement point; i ranges from 1 to n ; 计算得到测量应变矩阵一阶导数矩阵为：Calculate the first-order derivative matrix of the measured strain matrix for: （21）， (twenty one), 步骤4.3，根据光纤传感系统所得的测量应变矩阵和测量应变矩阵一阶导数，得到张力特征修正系数矩阵K：Step 4.3: Based on the measured strain matrix and the first-order derivative of the measured strain matrix obtained by the optical fiber sensing system, the tension characteristic correction coefficient matrix K is obtained: （22）， (twenty two), 步骤4.4，建立基于张力修正系数矩阵的管道张力识别方程。Step 4.4: Establish the pipeline tension identification equation based on the tension correction coefficient matrix. 7.根据权利要求6所述的方法，其特征在于，步骤4.1中，管道实际拉压应变计算公式为：7. The method according to claim 6, characterized in that in step 4.1, the actual tensile and compressive strains of the pipeline The calculation formula is: （23）， (twenty three), 其中，、、分别为0°、120°、240°光纤Bragg光栅传感器测得应变。in, 、 、 The strain was measured by fiber Bragg grating sensors at 0°, 120°, and 240° respectively. 8.根据权利要求7所述的方法，其特征在于，步骤4.1中，所述基于张力修正系数矩阵的管道张力识别方程为：8. The method according to claim 7, characterized in that in step 4.1, the pipeline tension identification equation based on the tension correction coefficient matrix is: （24）， (twenty four), （25）， (25), 其中Y(x)表示基于张力修正系数矩阵的管道线型函数；表示基于张力修正系数矩阵的管道结构张力识别函数。Where Y ( x ) represents the pipeline linear function based on the tension correction coefficient matrix; Represents the pipeline structure tension identification function based on the tension correction coefficient matrix. 9.一种电子设备，其特征在于，包括处理器和存储器，所述存储器存储有程序代码，当所述程序代码被所述处理器执行时，使得所述处理器执行如权利要求1至8中任一项所述的方法的步骤。9. An electronic device, comprising a processor and a memory, wherein the memory stores program code, and when the program code is executed by the processor, the processor executes the steps of the method according to any one of claims 1 to 8. 10.一种存储介质，其特征在于，存储有计算机程序或指令，当所述计算机程序或指令在计算机上运行时，执行如权利要求1至8中任一项所述的方法的步骤。10. A storage medium, characterized in that a computer program or instruction is stored therein, and when the computer program or instruction is run on a computer, the steps of the method according to any one of claims 1 to 8 are executed.
CN202510965609.0A 2025-07-14 2025-07-14 Pipeline structure tension identification method based on tension coefficient matrix correction Pending CN120449534A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202510965609.0A CN120449534A (en) 2025-07-14 2025-07-14 Pipeline structure tension identification method based on tension coefficient matrix correction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202510965609.0A CN120449534A (en) 2025-07-14 2025-07-14 Pipeline structure tension identification method based on tension coefficient matrix correction

Publications (1)

Publication Number Publication Date
CN120449534A true CN120449534A (en) 2025-08-08

Family

ID=96622791

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202510965609.0A Pending CN120449534A (en) 2025-07-14 2025-07-14 Pipeline structure tension identification method based on tension coefficient matrix correction

Country Status (1)

Country Link
CN (1) CN120449534A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050216242A1 (en) * 2002-05-20 2005-09-29 Michael Flax System and method for evaluation of fluid flow in a piping system
CN104899391A (en) * 2015-06-18 2015-09-09 哈尔滨工程大学 J-shaped pipe laying work model modeling simulation method
CN112100892A (en) * 2020-09-18 2020-12-18 哈尔滨工业大学（威海） Prediction method for natural frequency of in-band flow flexible pipeline under different boundary conditions
CN116384057A (en) * 2023-02-27 2023-07-04 南京航空航天大学 Pipeline morphological inversion error correction method based on strain correction coefficient matrix optimization

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050216242A1 (en) * 2002-05-20 2005-09-29 Michael Flax System and method for evaluation of fluid flow in a piping system
CN104899391A (en) * 2015-06-18 2015-09-09 哈尔滨工程大学 J-shaped pipe laying work model modeling simulation method
CN112100892A (en) * 2020-09-18 2020-12-18 哈尔滨工业大学（威海） Prediction method for natural frequency of in-band flow flexible pipeline under different boundary conditions
CN116384057A (en) * 2023-02-27 2023-07-04 南京航空航天大学 Pipeline morphological inversion error correction method based on strain correction coefficient matrix optimization

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHENG ZENG 等: "Shape reconstruction of aerial refueling hose based on distributed fiber optic sensors and FSResNet-MLP", FIRST CONFERENCE ON DISTRIBUTED OPTICAL FIBER SENSING TECHNOLOGY AND APPLICATIONS (DOFS 2024), 20 May 2025 (2025-05-20) *
周威;丁飞鹏;徐诚侃;: "含张力柔性管道的涡激振动计算方法", 科技通报, no. 04, 30 April 2019 (2019-04-30) *

Similar Documents

Publication Publication Date Title
CN106897543B (en) Beam structure damage identification method of modal compliance curvature matrix norm
CN111691679B (en) Digital twinning-based intelligent tensioning method for prestressed steel structure
CN106156501B (en) A kind of bridge local damage quantization method based on deflection line
CN119272213B (en) A method and system for monitoring abnormality of underground cable pipelines
CN109490072B (en) Detection system for civil engineering building and detection method thereof
CN116244886B (en) Virtual-real test data matching method and system
CN101806667B (en) Progressive method for identifying damaged cable and support displacement based on cable force monitoring
CN103049608A (en) System and method of load identification based on restrained edge strain extreme value coordinates
CN117787067B (en) A tower-type derrick bearing capacity determination method, system and device
CN106528945A (en) Identification method of load position and size of plate structure based on normalization of loading angle
CN117852122A (en) Method for detecting existing stress of steel strand in PC bridge structure
CN119004919B (en) Finite Element Simulation Analysis Method for Flexible Interconnection System
CN102323075B (en) Progressive method for identifying damaged cable, slack cable and support generalized displacement on basis of angle monitoring
CN120449534A (en) Pipeline structure tension identification method based on tension coefficient matrix correction
CN101793627A (en) Health monitoring method of increment cable structure based on hybrid monitoring
CN113449401A (en) Tower crane counter-force identification method
CN101793629B (en) Progressive cable structure health monitoring method based on strain monitoring
CN102323076A (en) Strain monitoring based progressive method for recognizing damaged cable, loose cable and supporting seat generalized displacement
CN102323083A (en) Strain monitoring based method for recognizing damaged cable, loose cable and supporting seat generalized displacement
CN101819097A (en) Progressive method for identifying damaged cable and support displacement based on strain monitoring
CN115525941A (en) Model-free detection method based on statistical moment ratio
CN102288440A (en) Progressive method for recognizing damaged cable, slack cable and angular displacement of support based on space coordinate monitoring
CN102297781A (en) Progressive method for identifying damaged cable, slacking cable and bearer angular displacement based on strain monitoring
CN101813571B (en) Progressive method for identifying damaged cables and support displacement based on angle monitoring
CN119203853B (en) Water level sensor gland and diaphragm assembly method combined with simulation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination