CN111625963B - Method and system for predicting residual life of pipeline based on hydrogen diffusion dynamics - Google Patents
Method and system for predicting residual life of pipeline based on hydrogen diffusion dynamics Download PDFInfo
- Publication number
- CN111625963B CN111625963B CN202010474166.2A CN202010474166A CN111625963B CN 111625963 B CN111625963 B CN 111625963B CN 202010474166 A CN202010474166 A CN 202010474166A CN 111625963 B CN111625963 B CN 111625963B
- Authority
- CN
- China
- Prior art keywords
- pipeline
- life
- crack growth
- growth rate
- loading frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 125000004435 hydrogen atom Chemical class [H]* 0.000 title claims abstract description 56
- 239000001257 hydrogen Substances 0.000 title claims abstract description 44
- 229910052739 hydrogen Inorganic materials 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000009792 diffusion process Methods 0.000 title claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000015654 memory Effects 0.000 claims description 22
- 238000004590 computer program Methods 0.000 claims description 11
- 239000000463 material Substances 0.000 claims description 9
- 238000012216 screening Methods 0.000 claims description 7
- 230000007613 environmental effect Effects 0.000 claims description 4
- 230000000704 physical effect Effects 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 25
- 150000002431 hydrogen Chemical class 0.000 description 15
- 230000000694 effects Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 230000007547 defect Effects 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005381 potential energy Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 230000005501 phase interface Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
The present disclosure discloses a method and a system for predicting the remaining life of a pipeline based on hydrogen diffusion dynamics, comprising: acquiring basic data of a pipeline to be subjected to life prediction; inputting the pipeline basic data to be subjected to life prediction into a screened pipeline crack growth rate calculation model, and outputting the pipeline crack growth rate; and predicting the residual life of the pipeline according to the initial crack size and the pipeline crack growth rate.
Description
Technical Field
The disclosure relates to the technical field of pipeline life prediction, in particular to a method and a system for predicting residual life of a pipeline based on hydrogen diffusion dynamics.
Background
The statements in this section merely mention background art related to the present disclosure and do not necessarily constitute prior art.
Conventional hydrogen embrittlement studies often only qualitatively discuss the embrittlement mechanism of hydrogen atoms to materials, mainly HELP (Hydrogen Enhanced Local Plasticity), HEDE (Hydrogen Enhanced Decohesion) and hydrogen bubble theory (local stress improvement caused by hydrogen gas generated by hydrogen atom aggregation in defects).
HELP theory emphasizes that hydrogen atoms can promote dislocation movement at crack tips, greatly improving localized plasticity. This view, at first glance, is against the theory of brittle deformation caused by hydrogen, however, when dislocations emitted from the crack tip encounter the grain boundaries, strain energy near the grain boundaries is greatly increased, thereby making the grain boundaries fracture easier to occur, which obviously greatly promotes the formation of brittle fracture.
Another view is that, contrary to HELP, hydrogen atoms can be concentrated into the stress concentration zone, i.e., the plastic zone at the crack tip, thereby impeding the generation of dislocations, making fracture the only way to unload strain energy.
HEDE theories believe that hydrogen atoms are more trapped near grain boundaries and phase interfaces. The hydrogen atoms not only lower the surface energy, making the free surface easier to form, but also increase the strain energy near the interface, making it easier for fracture to occur at the grain boundaries.
The hydrogen atom gathering area is more likely to generate small cracks, and hydrogen atoms can combine in the cracks or holes to generate hydrogen, so that the local stress is correspondingly improved, which is the origin of hydrogen bubble theory.
The destructive action caused by the hydrogen atoms is not irreversible before the crack grows, i.e. before the free surface of the crack tip is formed. By properly raising the temperature of the material, hydrogen atoms can be exuded from the free surface, thereby effectively inhibiting the hydrogen embrittlement.
The effect of temperature on hydrogen embrittlement is not linear: the temperature is too low, and the diffusion rate of hydrogen atoms is too slow, so that the hydrogen atoms are difficult to accumulate around the defect to play a role; too high a temperature, too fast a movement of hydrogen atoms, hydrogen permeation will occur, and thus the hydrogen embrittlement is weakened.
The stress intensity and the stress change rate also have great influence on the movement of hydrogen atoms, and the hydrogen atom enrichment rate of the crack tip is shown in a formula (1):
Where F r is the axial stress related to potential energy, K B is the Boltzmann constant, D is the diffusivity of the hydrogen atom, Ω is the component volume of the hydrogen atom, v is the Poisson's ratio of iron, and K I represents the stress intensity at the crack tip.
When the stress intensity becomes large, the hydrogen atoms move to the crack tip at a high rate, and crack growth is more likely to be caused. The quantitative relationship between the hydrogen enrichment and the stress intensity provides a theoretical basis for the quantitative relationship between the stress intensity and the crack growth rate. In addition, the rate of change of stress intensity is also important, and even if the stress intensity is large, hydrogen atoms cannot be enriched for a sufficient time, and the influence of hydrogen embrittlement on crack growth is not obvious when too frequent unloading and loading operations exist.
Based on the above theory, an initial hydrogen diffusion crack propagation model was tried to quantitatively describe the effect of hydrogen atom concentration and diffusion rate on crack growth.
The rate of crack growth was obtained by calculating the rate of hydrogen atoms diffusing into the crack tip to fill the free surface at an atomic ratio of 1, and deriving the mechanical model as follows:
Wherein, Is the crack growth rate, Δa is the crack length change, t is time, and c o is the hydrogen atom concentration. However, the purely mechanical model cannot be applied to actual pipe fracture prediction due to the application conditions under which it is too idealized (ignoring the effects of plastic regions and fatigue). In practice, crack growth predictions are often solved by empirical models, such as the Paris model:
The model has obvious defects, and neglects the influence of external conditions such as loading frequency and the like on crack growth. The superposition model is shown in formula (4). As a complement to the Paris model, loading frequencies were introduced into the model, but the effect on loading frequencies was not given an accurate quantification, nor was there a clear solution to the relationship of stress corrosion (SCC) and fatigue (fatigue).
The joint factor model is shown in formula (5). The model is a model which is relatively complete in an empirical model and has high prediction accuracy. However, the model ignores the threshold value of the influence of the loading frequency on the crack growth, and cannot accurately calculate the influence of the external environment such as temperature and the concentration of hydrogen atoms on the crack growth.
In summary, the existing technology cannot accurately predict the remaining life of the pipeline, and the existing technology predicts the remaining life of the pipeline and lacks adaptability to changing environments.
Disclosure of Invention
In order to solve the defects of the prior art, the present disclosure provides a method and a system for predicting the residual life of a pipeline based on hydrogen diffusion kinetics;
in a first aspect, the present disclosure provides a method of predicting a remaining life of a pipeline based on hydrogen diffusion kinetics;
The method for predicting the residual life of the pipeline based on hydrogen diffusion dynamics comprises the following steps:
Acquiring basic data of a pipeline to be subjected to life prediction;
Inputting the pipeline basic data to be subjected to life prediction into a screened pipeline crack growth rate calculation model, and outputting the pipeline crack growth rate;
and predicting the residual life of the pipeline according to the initial crack size and the pipeline crack growth rate.
In a second aspect, the present disclosure provides a system for predicting the remaining life of a pipeline based on hydrogen diffusion kinetics;
A hydrogen diffusion dynamics based pipe remaining life prediction system comprising:
an acquisition module configured to: acquiring basic data of a pipeline to be subjected to life prediction;
An input module configured to: inputting the pipeline basic data to be subjected to life prediction into a screened pipeline crack growth rate calculation model, and outputting the pipeline crack growth rate;
A prediction module configured to: and predicting the residual life of the pipeline according to the initial crack size and the pipeline crack growth rate.
In a third aspect, the present disclosure also provides an electronic device, including: one or more processors, one or more memories, and one or more computer programs; wherein the processor is coupled to the memory, the one or more computer programs being stored in the memory, the processor executing the one or more computer programs stored in the memory when the electronic device is running, to cause the electronic device to perform the method of the first aspect.
In a fourth aspect, the present disclosure also provides a computer readable storage medium storing computer instructions which, when executed by a processor, perform the method of the first aspect.
In a fifth aspect, the present disclosure also provides a computer program (product) comprising a computer program for implementing the method of any one of the preceding aspects when run on one or more processors.
Compared with the prior art, the beneficial effects of the present disclosure are:
The present disclosure discloses a crack growth model and a program prototype based on hydrogen embrittlement theory, the method comprising: through theoretical deduction, a crack growth prediction formula based on the stress intensity change of the pipeline, the hydrogen potential energy, the hydrogen diffusivity and the pipeline working environment is established; fitting the model with the existing test data, correcting parameters and improving the model precision; and writing the theoretical model into a C language code to realize the large data processing of the pipeline transmission pressure and the prediction of the service life of the pipeline. The hydrogen embrittlement crack growth model is based on a hydrogen embrittlement crack growth model, the model theory is solid, the prediction is accurate, the unpredictability of the influence of environmental conditions such as temperature on the service life of the pipeline is solved, and the hydrogen embrittlement crack growth model has great significance for the integrity management and risk evaluation of the pipeline.
Drawings
The accompanying drawings, which are included to provide a further understanding of the disclosure, illustrate and explain the exemplary embodiments of the disclosure and together with the description serve to explain the disclosure, and do not constitute an undue limitation on the disclosure.
FIG. 1 is a flow chart of a method of a first embodiment;
FIG. 2 is a virtual loading spectrum of a first embodiment;
FIG. 3 is a graph showing peak and valley values of a first embodiment;
FIG. 4 is a schematic representation of the final crack growth rate of the first embodiment;
Detailed Description
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the present disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments in accordance with the present disclosure. As used herein, unless the context clearly indicates otherwise, the singular forms also are intended to include the plural forms, and furthermore, it is to be understood that the terms "comprises" and "comprising" and any variations thereof are intended to cover non-exclusive inclusions, such as, for example, processes, methods, systems, products or devices that comprise a series of steps or units, are not necessarily limited to those steps or units that are expressly listed, but may include other steps or units that are not expressly listed or inherent to such processes, methods, products or devices.
Embodiments of the present disclosure and features of embodiments may be combined with each other without conflict.
Example 1
As shown in fig. 1, the present embodiment provides a pipe remaining life prediction method based on hydrogen diffusion dynamics;
The method for predicting the residual life of the pipeline based on hydrogen diffusion dynamics comprises the following steps:
s101: acquiring basic data of a pipeline to be subjected to life prediction;
S102: inputting the pipeline basic data to be subjected to life prediction into a screened pipeline crack growth rate calculation model, and outputting the pipeline crack growth rate;
S103: and predicting the residual life of the pipeline according to the initial crack size and the pipeline crack growth rate.
As one or more embodiments, the screening step includes:
and comparing the actual loading frequency of the pipeline to be life-predicted with a pipeline loading frequency threshold value, and screening out a corresponding pipeline crack growth rate calculation model.
As one or more embodiments, the screened pipeline crack growth rate calculation model includes:
And constructing a corresponding pipeline crack growth rate calculation model according to the relation between the actual pipeline loading frequency and the pipeline loading frequency threshold value.
As one or more embodiments, the pipeline loading frequency threshold, the obtaining step includes:
and calculating a pipeline loading frequency threshold according to the poisson ratio of the pipeline to be predicted to be life-time, the yield strength of the pipeline to be life-time, the hydrogen atom component volume of the pipeline to be life-time, the hydrogen atom diffusivity of the pipeline to be life-time and the temperature of the pipeline to be life-time.
It should be understood that the poisson ratio of the pipeline refers to: poisson's ratio refers to the ratio of the absolute value of the positive transverse strain to the positive axial strain of a material under unidirectional tension or compression, and is a dimensionless value.
It should be understood that the yield strength of the pipe refers to: the yield strength is the yield limit of the metallic material when it undergoes a yield phenomenon, i.e. the stress against a minor plastic deformation, in Pa.
It should be understood that the volume of the hydrogen atomic component of the pipeline refers to: the individual hydrogen atoms are distributed in the volume occupied in the pipeline steel in units of m 3.
It should be understood that the pipe hydrogen atom diffusivity refers to: the diffusivity is a diffusion rate, which is a measure of the propagation velocity of hydrogen atoms, in m 2/s.
It should be understood that the pipe temperature refers to: the actual temperature of the pipeline is predicted after the service life, and the unit is K.
As one or more embodiments, calculating a pipeline loading frequency threshold according to poisson ratio of the pipeline to be predicted to be life-time, yield strength of the pipeline to be life-time, hydrogen atom component volume of the pipeline to be life-time, hydrogen atom diffusivity of the pipeline to be life-time and temperature of the pipeline to be life-time; the method comprises the following specific steps:
where v is poisson's ratio, Ω is hydrogen atom component volume, D is hydrogen atom diffusivity, K max and K min are maximum and minimum stress intensity values for a single cyclic load, R p is plastic denaturation zone size, R eq is plastic zone hydrogen supply zone size, K B is boltzmann constant, and T is temperature.
As one or more embodiments, constructing a corresponding pipeline crack growth rate calculation model according to the relation between the actual loading frequency of the pipeline and the pipeline loading frequency threshold value; the method comprises the following specific steps:
Wherein R is the ratio of K min to K max, A is an environmental factor, deltaK is the difference between K max and K min, and n is a physical property constant related to the material, and the material is obtained by fitting.
As one or more embodiments, comparing the actual loading frequency of the pipeline to be life-predicted with a pipeline loading frequency threshold value, and screening out a corresponding pipeline crack growth rate calculation model; the method comprises the following specific steps: selecting formula (7) when the actual loading frequency of the pipeline to be life-predicted is greater than the pipeline loading frequency threshold; and when the actual loading frequency of the pipeline to be life-predicted is less than or equal to the pipeline loading frequency threshold value, selecting the formula (8).
As one or more embodiments, the initial crack size calculating step includes:
Scanning to obtain detection data in the pipeline through a temperature sensor, a stress sensor and a pipe cleaner;
the in-pipe inspection data includes: running the length of the microcrack in the pipeline, the service temperature of the pipeline and the stress variation intensity variation of the pipeline;
and acquiring the initial crack size of the pipeline according to the detection data in the pipeline.
As one or more embodiments, pipeline base data for life prediction, comprising: the temperature of the pipeline to be predicted for the service life, the concentration of the hydrogen atoms of the matrix of the pipeline to be predicted for the service life, the Poisson ratio of the pipeline to be predicted for the service life and the yield strength of the pipeline to be predicted for the service life.
It should be understood that the concentration of the base hydrogen atoms of the pipe means: atomic percent of hydrogen atoms in the pipeline steel.
As one or more embodiments, the remaining life of the pipe is predicted based on the initial crack size and the pipe crack growth rate; the method comprises the following specific steps:
where a is the predicted crack length, a o is the initial crack length, N is the number of loading cycles, and da/dN is the crack growth length per cycle, which can be found by either equation (7) or (8).
The explanation of the hydrogen embrittlement theory comprises HEDE (hydrogen enhanced decohesion), HELP (hydrogen enhanced local plasticity) and HBT (hydrogen bubble theory), the core idea is that hydrogen atoms are gathered around defects such as crack tips, dislocation and the like, free surface energy is reduced, local stress and local deformation are increased, so that the steel material is more prone to failure, the relation between the gathering rate of the hydrogen atoms around the defects and the crack growth rate is quantified by a hydrogen embrittlement fracture model based on the theory, and the crack growth rate is determined; the relation between the section forming rate and the stress intensity is quantified firstly; the influence of the loading frequency has a limit above which the rate increases with increasing loading frequency, whereas below this limit the rate does not change with time; the effect of temperature and the effect of hydrogen atom concentration (i.e., pH-worthy effect) were quantified.
The crack growth is predicted by a formula which, given the change in the pipeline transport stress intensity,
Analyzing the change of the stress intensity distribution of the original data pipeline corresponding to the time: fig. 2 utilizes a virtual loading spectrogram for convenient presentation. Such regular and fluctuating stress intensity variations do not exist in actual pipeline transportation.
The effect of loading above 0.94 of the R ratio (K min/Kmax) on crack growth is negligible, and the corresponding stress intensity threshold is 0.594MPa x m 0.5.
The loading period in fig. 1 where the stress intensity range is greater than or equal to 3.3mpa x m 0.5 is determined to be a non-small period loading. For ease of illustration, no small cycle loading is provided in fig. 1.
The loading poling is broken down into a number of basic modules, each module comprising a large loading period with a ak higher than 3.3mpa x m 0.5 and a number of small loading periods with a ak distribution between 0.594 and 3.3mpa x m 0.5 before the occurrence of the next such loading period.
And recording the screened key data into a file, wherein the original data is processed in the figure 3, and the recorded time corresponds to the corresponding relation diagram of K max and K min. Since all loads Δk are greater than or equal to 3.3mpa x m 0.5, the peak and trough values shown in fig. 3 are all recorded.
(1) After the original data processing, the loading frequency f needs to be recorded.
(2) K max and Δk were recorded for each larger loading cycle, and since only the case where Δk was greater than or equal to 3.3mpa x m 0.5 was present in the initial spectrum we provided, the stress intensity values for all cycles were recorded as shown in fig. 2.
(3) The stress intensity distribution in each base load module is between 0.594 and 3.3mpa x m 0.5 small cycles (n). There is no small period in fig. 2, so fig. 3 does not show this screening condition.
(4) The temperature is provided by the pump house of the loading station.
(5) The basic hydrogen atom concentration c o in the pipeline was obtained by experiment and simulation.
Calculating crack growth rate of each basic loading module (comprising a large loading period with Δk higher than 3.3mpa×m 0.5 +n small loading periods with Δk distributed between 0.594 and 3.3mpa×m 0.5, n being 0 in this example);
There is a threshold for loading frequency variation in the pipe run, below which the crack growth rate is maximized and does not vary with loading frequency. The relation between crack growth rate and loading frequency, for loading in the same stress intensity range, the greater the loading frequency, the slower the crack growth rate. The threshold value of the loading frequency is f critical, which is a function of the corresponding loading mode, the physical properties of the material (poisson ratio v, hardness), the hydrogen atom diffusivity (component volume omega, diffusivity D) and the temperature T, as shown in a formula (6). After the processed data, the frequency is compared with f critical, and then the crack propagation formulas (7) and (8) based on a hydrogen atom diffusion model are carried in.
The calculation can be performed to obtain the crack growth rate in a single loading module. Wherein A is a factor related to temperature, hydrogen atom concentration, physical properties of the material:
If the small period acceleration condition is considered, the small period acceleration factor can be expressed as follows:
Wherein n is the number of small loading cycles with delta K distributed between 0.594 and 3.3MPa m 0.5, c o is the hydrogen atom specific concentration of the pipeline matrix, n critical and c critical are pipeline characteristic values which are 657 and 0.01568 respectively, and when the small loading cycles exist, the acceleration factor is introduced into a formula (8) to form Can be used for predicting the real crack growth rate. Since there is no small period in this example, the small period acceleration condition is negligible and the final crack growth rate is shown in fig. 4:
FIG. 4 is a graph of crack length versus number of loading modules, where there are 10 smaller modules per larger module, i.e., between the larger cycles in FIG. 1, and thus a stepwise increase in crack growth occurs. By combining theoretical research, the relation between a hydrogen atom diffusion fracture model mechanism and industrial reality is fully considered, and the theoretical model is updated and reproduced, so that the novel mechanical model can be applied to the prediction of pipeline fracture in soil solution, and the theoretical model is arranged into a program, so that the model is easier to popularize and use, and has important significance for evaluating the integrity of the pipeline. The purpose of the present disclosure is to build a pipeline crack growth model based on hydrogen embrittlement theory, and aims to quantify the dependency relationship between the service life of a pipeline and the environmental temperature and the concentration of hydrogen atoms, and to increase the reliability and applicability of prediction.
Example two
The embodiment provides a pipeline residual life prediction system based on hydrogen diffusion dynamics;
A hydrogen diffusion dynamics based pipe remaining life prediction system comprising:
an acquisition module configured to: acquiring basic data of a pipeline to be subjected to life prediction;
An input module configured to: inputting the pipeline basic data to be subjected to life prediction into a screened pipeline crack growth rate calculation model, and outputting the pipeline crack growth rate;
A prediction module configured to: and predicting the residual life of the pipeline according to the initial crack size and the pipeline crack growth rate.
Here, it should be noted that the above-mentioned obtaining module, input module and prediction module correspond to steps S101 to S103 in the first embodiment, and the above-mentioned modules are the same as examples and application scenarios implemented by the corresponding steps, but are not limited to the disclosure of the first embodiment. It should be noted that the modules described above may be implemented as part of a system in a computer system, such as a set of computer-executable instructions.
The foregoing embodiments are directed to various embodiments, and details of one embodiment may be found in the related description of another embodiment.
The proposed system may be implemented in other ways. For example, the system embodiments described above are merely illustrative, such as the division of the modules described above, are merely a logical function division, and may be implemented in other manners, such as multiple modules may be combined or integrated into another system, or some features may be omitted, or not performed.
Example III
The embodiment also provides an electronic device, including: one or more processors, one or more memories, and one or more computer programs; wherein the processor is coupled to the memory, the one or more computer programs being stored in the memory, the processor executing the one or more computer programs stored in the memory when the electronic device is running, to cause the electronic device to perform the method of the first embodiment.
It should be understood that in this embodiment, the processor may be a central processing unit CPU, and the processor may also be other general purpose processors, digital signal processors DSP, application specific integrated circuits ASIC, off-the-shelf programmable gate array FPGA or other programmable logic device, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may include read only memory and random access memory and provide instructions and data to the processor, and a portion of the memory may also include non-volatile random access memory. For example, the memory may also store information of the device type.
In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or by instructions in the form of software.
The method in the first embodiment may be directly implemented as a hardware processor executing or implemented by a combination of hardware and software modules in the processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in a memory, and the processor reads the information in the memory and, in combination with its hardware, performs the steps of the above method. To avoid repetition, a detailed description is not provided herein.
Those of ordinary skill in the art will appreciate that the elements of the various examples described in connection with the present embodiments, i.e., the algorithm steps, can be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
Example IV
The present embodiment also provides a computer-readable storage medium storing computer instructions that, when executed by a processor, perform the method of embodiment one.
The foregoing description of the preferred embodiments of the present disclosure is provided only and not intended to limit the disclosure so that various modifications and changes may be made to the present disclosure by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present disclosure should be included in the protection scope of the present disclosure.
Claims (5)
1. The method for predicting the residual life of the pipeline based on hydrogen diffusion dynamics is characterized by comprising the following steps of:
Acquiring basic data of a pipeline to be subjected to life prediction;
Inputting the pipeline basic data to be subjected to life prediction into a screened pipeline crack growth rate calculation model, and outputting the pipeline crack growth rate;
The screened pipeline crack growth rate calculation model comprises the following screening steps: comparing the actual loading frequency of the pipeline to be predicted for service life with a pipeline loading frequency threshold value, and screening out a corresponding pipeline crack growth rate calculation model; the step of obtaining the pipeline loading frequency threshold value comprises the following steps: calculating a pipeline loading frequency threshold according to the poisson ratio of the pipeline to be predicted to be life-time, the yield strength of the pipeline to be life-time, the hydrogen atom component volume of the pipeline to be life-time, the hydrogen atom diffusivity of the pipeline to be life-time and the temperature of the pipeline to be life-time; the method for calculating the pipeline loading frequency threshold according to the poisson ratio of the pipeline to be predicted for the service life, the yield strength of the pipeline to be predicted for the service life, the hydrogen atom component volume of the pipeline to be predicted for the service life, the hydrogen atom diffusivity of the pipeline to be predicted for the service life and the temperature of the pipeline to be predicted for the service life comprises the following specific steps:
；（6）
Wherein, Is Poisson's ratio,/>Is hydrogen atom component volume,/>Is hydrogen atom diffusivity,/>And/>Maximum and minimum stress intensity values for a single cyclic load,/>Size of plastic denatured area,/>Size of hydrogen-donating region for plastic region,/>Is Boltzmann constant,/>Is the temperature;
The screened pipeline crack growth rate calculation model comprises the following screening ranges: constructing a corresponding pipeline crack growth rate calculation model according to the relation between the actual loading frequency of the pipeline and the pipeline loading frequency threshold value; according to the relation between the actual loading frequency of the pipeline and the loading frequency threshold value of the pipeline, a corresponding pipeline crack growth rate calculation model is constructed, and the concrete steps comprise:
；（7）
；（8）
Wherein, For/>And/>Ratio of/(I)Is an environmental factor,/>For/>And/>Difference of/>Is a physical property constant related to the material and is obtained by fitting;
and predicting the residual life of the pipeline according to the initial crack size and the pipeline crack growth rate.
2. The method of claim 1, wherein the actual loading frequency of the pipeline to be life-predicted is compared with a pipeline loading frequency threshold value, and a corresponding pipeline crack growth rate calculation model is screened out; the method comprises the following specific steps: selecting formula (7) when the actual loading frequency of the pipeline to be life-predicted is greater than the pipeline loading frequency threshold; and when the actual loading frequency of the pipeline to be life-predicted is less than or equal to the pipeline loading frequency threshold value, selecting the formula (8).
3. A hydrogen diffusion dynamics-based pipe remaining life prediction system employing the hydrogen diffusion dynamics-based pipe remaining life prediction method according to claim 1, comprising:
an acquisition module configured to: acquiring basic data of a pipeline to be subjected to life prediction;
An input module configured to: inputting the pipeline basic data to be subjected to life prediction into a screened pipeline crack growth rate calculation model, and outputting the pipeline crack growth rate;
A prediction module configured to: and predicting the residual life of the pipeline according to the initial crack size and the pipeline crack growth rate.
4. An electronic device, comprising: one or more processors, one or more memories, and one or more computer programs; wherein the processor is coupled to the memory, the one or more computer programs being stored in the memory, which processor, when the electronic device is running, executes the one or more computer programs stored in the memory to cause the electronic device to perform the method of any of the preceding claims 1-2.
5. A computer readable storage medium storing computer instructions which, when executed by a processor, perform the method of any of claims 1-2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010474166.2A CN111625963B (en) | 2020-05-29 | 2020-05-29 | Method and system for predicting residual life of pipeline based on hydrogen diffusion dynamics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010474166.2A CN111625963B (en) | 2020-05-29 | 2020-05-29 | Method and system for predicting residual life of pipeline based on hydrogen diffusion dynamics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111625963A CN111625963A (en) | 2020-09-04 |
| CN111625963B true CN111625963B (en) | 2024-05-28 |
Family
ID=72259248
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010474166.2A Active CN111625963B (en) | 2020-05-29 | 2020-05-29 | Method and system for predicting residual life of pipeline based on hydrogen diffusion dynamics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111625963B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112666066B (en) * | 2020-12-15 | 2022-11-11 | 中国石油大学(华东) | Pipeline hydrogen embrittlement temperature threshold prediction method based on hydrogen diffusion dynamics and application |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104142272A (en) * | 2014-07-22 | 2014-11-12 | 广东电网公司电力科学研究院 | Ultra supercritical boiler super-heat and re-heater heating surface life evaluation method |
| CN106126856A (en) * | 2016-07-06 | 2016-11-16 | 西南交通大学 | It is applicable to the crack growth rate Forecasting Methodology based on low cycle fatigue property parameter of negative stress ratio |
| CN109740295A (en) * | 2019-02-27 | 2019-05-10 | 南京市特种设备安全监督检验研究院 | A kind of vibrative mechanism residual Life Calculation method with crack defect |
-
2020
- 2020-05-29 CN CN202010474166.2A patent/CN111625963B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104142272A (en) * | 2014-07-22 | 2014-11-12 | 广东电网公司电力科学研究院 | Ultra supercritical boiler super-heat and re-heater heating surface life evaluation method |
| CN106126856A (en) * | 2016-07-06 | 2016-11-16 | 西南交通大学 | It is applicable to the crack growth rate Forecasting Methodology based on low cycle fatigue property parameter of negative stress ratio |
| CN109740295A (en) * | 2019-02-27 | 2019-05-10 | 南京市特种设备安全监督检验研究院 | A kind of vibrative mechanism residual Life Calculation method with crack defect |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111625963A (en) | 2020-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tekoğlu et al. | On localization and void coalescence as a precursor to ductile fracture | |
| Farhidzadeh et al. | Monitoring crack propagation in reinforced concrete shear walls by acoustic emission | |
| Tekoglu et al. | A criterion for the onset of void coalescence under combined tension and shear | |
| Sagar et al. | A review of recent developments in parametric based acoustic emission techniques applied to concrete structures | |
| Hardin et al. | Effect of porosity on deformation, damage, and fracture of cast steel | |
| Imanian et al. | A thermodynamic entropy approach to reliability assessment with applications to corrosion fatigue | |
| Pommier et al. | Incremental model for fatigue crack growth based on a displacement partitioning hypothesis of mode I elastic–plastic displacement fields | |
| Pourmodheji et al. | Improvement of the extended finite element method for ductile crack growth | |
| Zhao et al. | Collapse pressure of randomly corroded spherical shell | |
| Marsili et al. | Seismic reliability assessment of a concrete water tank based on the Bayesian updating of the finite element model | |
| CN111625963B (en) | Method and system for predicting residual life of pipeline based on hydrogen diffusion dynamics | |
| WO2023155349A1 (en) | Method and apparatus for predicting service life of steel box girder, device, and medium | |
| Abyani et al. | A new approach for finite element based reliability evaluation of offshore corroded pipelines | |
| Maraveas | Concrete silos: Failures, design issues and repair/strengthening methods | |
| Ma et al. | Computational investigation of multi-axial damage modeling for porous sintered metals with experimental verification | |
| Xin et al. | Effects of galvanostatic and artificial chloride environment methods on the steel corrosion spatial variability and probabilistic flexural capacity of RC beams | |
| Chiu et al. | Probability-based damage assessment for reinforced concrete bridge columns considering the corrosive and seismic hazards in Taiwan | |
| Walz et al. | Probabilistic fracture mechanics assessment of flaws in turbine disks including quality assurance procedures | |
| Wee et al. | Stochastic study on effects of material and physical parameters on slow crack growth behaviors of HDPE using the crack layer theory | |
| Tian et al. | Experimental study on bond performance and damage detection of corroded reinforced concrete specimens | |
| Yang et al. | Generalized hysteretic constitutive model of marine structural steel with pitting corrosion | |
| Milenin et al. | Residual strength and reliability of corroded pipelines—Monte-Carlo approach for consideration of spatially nonuniform material properties | |
| Tian et al. | Fracture and fatigue analyses of cracked structures using the iterative method | |
| Yuan et al. | Effect of Reinforcement Corrosion Sediment Distribution Characteristics on Concrete Damage Behavior. | |
| Meng et al. | Prediction of optimal inspection time for structural fatigue life |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |