CN111625963A - Pipeline residual life prediction method and system based on hydrogen diffusion dynamics - Google Patents

Pipeline residual life prediction method and system based on hydrogen diffusion dynamics Download PDF

Info

Publication number
CN111625963A
CN111625963A CN202010474166.2A CN202010474166A CN111625963A CN 111625963 A CN111625963 A CN 111625963A CN 202010474166 A CN202010474166 A CN 202010474166A CN 111625963 A CN111625963 A CN 111625963A
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.)
Granted
Application number
CN202010474166.2A
Other languages
Chinese (zh)
Other versions
CN111625963B (en
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.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202010474166.2A priority Critical patent/CN111625963B/en
Publication of CN111625963A publication Critical patent/CN111625963A/en
Application granted granted Critical
Publication of CN111625963B publication Critical patent/CN111625963B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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 kinetics, comprising: acquiring basic data of a pipeline to be subjected to life prediction; inputting the basic data of the pipeline with the service life to be predicted into the 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

Pipeline residual life prediction method and system based on hydrogen diffusion dynamics
Technical Field
The disclosure relates to the technical field of pipeline life prediction, in particular to a method and a system for predicting the residual life of a pipeline based on hydrogen diffusion dynamics.
Background
The statements in this section merely provide background information related to the present disclosure and may not constitute prior art.
The traditional hydrogen embrittlement research can only qualitatively discuss the embrittlement mechanism of hydrogen atoms on materials, mainly including theory of HELP (hydrogen Enhanced Local property), HEDE (hydrogen Enhanced dehydrogenation) and hydrogen bubble resonator (Local stress improvement caused by hydrogen generated by hydrogen atoms gathering in defects).
The HELP theory emphasizes that hydrogen atoms can promote dislocation movement at the crack tip, greatly enhancing local plasticity. This view is at first glance contradictory to the theory that hydrogen causes brittle deformation, however, when dislocations emitted from the crack tip meet the grain boundary, the strain energy near the grain boundary is greatly increased, so that intergranular fracture is more likely to occur, which apparently greatly promotes the formation of brittle fracture.
Another aspect is that, in contrast to HELP, hydrogen atoms tend to concentrate in the stress concentration zone, i.e., the plastic zone at the crack tip, thereby hindering the generation of dislocations, making fracture the only way to unload the strain energy.
HEDE theorizes that hydrogen atoms are more trapped near grain boundaries and phase boundaries. Hydrogen atoms not only reduce the surface energy, making the free surface more easily formed, but also increase the strain energy near the interface, making the fracture more likely to occur at the grain boundaries.
The hydrogen atom accumulation region is more prone to small cracks, and hydrogen atoms can combine in cracks or holes to generate hydrogen gas, in which case the local stress is correspondingly increased, which is the origin of the hydrogen bubble theory.
The destructive effect caused by hydrogen atoms is not irreversible before the crack growth, i.e. before the crack tip free surface is formed. By appropriately raising the material temperature, hydrogen atoms can be made to exude from the free surface, thereby effectively suppressing 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 gather around the defects to play a role; too high a temperature causes hydrogen atoms to move too fast, and hydrogen permeation occurs, thereby weakening hydrogen embrittlement.
The stress intensity and the stress change rate also greatly affect the movement of hydrogen atoms, and the hydrogen atom enrichment rate of the crack tip is shown in the formula (1):
Figure BDA0002515307920000021
wherein, FrIs the axial stress, k, related to potential energyBIs Boltzmann's constant, D is the diffusivity of a hydrogen atom, Ω is the component volume of a hydrogen atom, v is the Poisson's ratio of iron, KIRepresenting the stress intensity at the crack tip.
When the stress intensity becomes large, the hydrogen atoms move to the crack tip at a high speed, and the crack growth is more easily caused. The quantitative relation of hydrogen enrichment and stress intensity provides a theoretical basis for the quantitative relation of the stress intensity and the crack growth rate. In addition, the rate of change of the stress intensity is also important, and when too frequent unloading and loading operations exist, 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.
Based on the above theory, an initial hydrogen diffusion crack propagation model, attempts 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 at which hydrogen atoms diffuse to the crack tip to fill the free surface at an atomic ratio of 1, and deriving the mechanical model as follows:
Figure BDA0002515307920000031
wherein the content of the first and second substances,
Figure BDA0002515307920000032
is the crack growth rate and Δ a is the crack lengthChange, t is time, coIs the hydrogen atom concentration. However, due to the highly idealized application conditions of the purely mechanical model (neglecting the effects of plasticity zones and fatigue), the model cannot be applied to actual pipe break prediction. In practice, crack growth prediction is often solved by empirical models, such as Paris model:
Figure BDA0002515307920000033
the model has obvious defects, and influences of external conditions such as loading frequency and the like on crack growth are ignored. The overlay model is shown in equation (4). In addition to the Paris model, loading frequency is introduced into the model, but the effect on the loading frequency is not accurately quantified, and no explicit solution is provided for the relationship between stress corrosion (SCC) and fatigue (fatigue).
Figure BDA0002515307920000034
The joint factor model is shown in equation (5). The model is a relatively complete model with high prediction accuracy in the empirical model. However, the model neglects that there is a threshold for the influence of the loading frequency on the crack growth, and the influence of the external environment, such as temperature, and the concentration of hydrogen atoms on the crack growth cannot be accurately calculated.
Figure BDA0002515307920000035
In conclusion, the prior art cannot accurately predict the residual service life of the pipeline, and the prediction of the residual service life of the pipeline in the prior art is lack of adaptability to the changing environment.
Disclosure of Invention
In order to solve the deficiencies of the prior art, the present disclosure provides a method and system for predicting the remaining life of a pipeline based on hydrogen diffusion kinetics;
in a first aspect, the present disclosure provides a method for predicting 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 basic data of the pipeline with the service life to be predicted into the 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 remaining life of a pipeline based on hydrogen diffusion kinetics;
a system for predicting remaining life of a pipeline based on hydrogen diffusion kinetics comprises:
an acquisition module configured to: acquiring basic data of a pipeline to be subjected to life prediction;
an input module configured to: inputting the basic data of the pipeline with the service life to be predicted into the 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 a processor is connected to the memory, the one or more computer programs are stored in the memory, and when the electronic device is running, the processor executes the one or more computer programs stored in the memory, so as to make the electronic device execute the method according to the first aspect.
In a fourth aspect, the present disclosure also provides a computer-readable storage medium for 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 first aspects when run on one or more processors.
Compared with the prior art, the beneficial effect of this disclosure is:
the present disclosure discloses a crack growth model and a procedure prototype based on hydrogen embrittlement theory, the method comprises: establishing a crack growth prediction formula based on the stress intensity change of the pipeline, the hydrogen potential energy, the hydrogen diffusivity and the working environment of the pipeline through theoretical derivation; fitting the model with the existing test data, correcting parameters and improving the precision of the model; the theoretical model is written into C language code, so that the large data processing of the pipe transmission pressure and the prediction of the service life of the pipeline are realized. The method 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 method has great significance on the integrity management and the risk evaluation of the pipeline.
Drawings
The accompanying drawings, which are included to provide a further understanding of the disclosure, illustrate embodiments of the disclosure and together with the description serve to explain the disclosure and are not to limit the disclosure.
FIG. 1 is a flow chart of a method of the first embodiment;
FIG. 2 is a diagram of a virtual load spectrogram of the first embodiment;
FIG. 3 is a schematic illustration of the peak and valley values of the first embodiment;
FIG. 4 is a graph illustrating the final crack propagation 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 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 example embodiments according to the present disclosure. As used herein, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise, and it should be understood that the terms "comprises" and "comprising", and any variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
The embodiments and features of the embodiments in the present disclosure may be combined with each other without conflict.
Example one
As shown in fig. 1, the present embodiment provides a method for predicting the remaining life of a pipe based on hydrogen diffusion kinetics;
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 basic data of the pipeline with the service life to be predicted into the 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.
In one or more embodiments, the screened pipeline crack growth rate calculation model includes:
and comparing the actual loading frequency of the pipeline to be subjected to life prediction with a pipeline loading frequency threshold value, and screening out a corresponding pipeline crack growth rate calculation model.
As one or more embodiments, the selected pipeline crack growth rate calculation model includes:
and constructing a corresponding pipeline crack growth rate calculation model according to the relation between the actual loading frequency of the pipeline and the threshold value of the loading frequency of the pipeline.
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 subjected to life prediction, the yield strength of the pipeline to be subjected to life prediction, the volume of the hydrogen atom component of the pipeline to be subjected to life prediction, the hydrogen atom diffusivity of the pipeline to be subjected to life prediction and the temperature of the pipeline to be subjected to life prediction.
It is understood that the poisson's ratio of the conduit means: the poisson ratio is the ratio of the absolute value of transverse positive strain and axial positive 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 means: the yield strength is the yield limit at which the metal material yields, i.e. the stress resisting a slight amount of plastic deformation, in Pa.
It should be understood that the pipeline hydrogen atom component volume refers to: volume of single hydrogen atom distributed in pipeline steel, unit m3。
It should be understood that the hydrogen atom diffusivity of the pipeline refers to: diffusivity is a diffusion rate, a measure of the speed of propagation of a hydrogen atom, and is expressed in m2/s。
It should be understood that the pipe temperature refers to: and predicting the actual temperature of the pipeline to be subjected to life prediction, wherein the unit is K.
As one or more embodiments, calculating a pipeline loading frequency threshold according to the poisson ratio of the pipeline to be subjected to life prediction, the yield strength of the pipeline to be subjected to life prediction, the volume of hydrogen atom components of the pipeline to be subjected to life prediction, the hydrogen atom diffusivity of the pipeline to be subjected to life prediction and the temperature of the pipeline to be subjected to life prediction; the method comprises the following specific steps:
Figure BDA0002515307920000071
wherein v is Poisson's ratio, Ω is hydrogen atom component volume, D is hydrogen atom diffusivity, KmaxAnd KminMaximum and minimum stress intensity values for a single cyclic load, rpFor the size of the plastic deformation zone, ReqSize of hydrogen-donating region, k, for plastic regionBBoltzmann constant, 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 threshold value of the loading frequency of the pipeline; the method comprises the following specific steps:
Figure BDA0002515307920000081
Figure BDA0002515307920000082
wherein R is KminAnd KmaxA is an environmental factor and Δ K is KmaxAnd KminN is a material-related property constant, and needs to be obtained by fitting.
As one or more embodiments, the actual loading frequency of the pipeline to be subjected to life prediction 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: when the actual loading frequency of the pipeline to be subjected to life prediction is larger than the pipeline loading frequency threshold value, selecting a formula (7); and when the actual loading frequency of the pipeline to be subjected to life prediction 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, the calculating step comprises:
acquiring detection data in the pipeline through scanning of a temperature sensor, a stress sensor and a pipe cleaner;
the in-pipeline detection data comprises: the length of the micro-crack inside the operation pipeline, the service temperature of the pipeline and the change strength of the stress of the pipeline are changed;
and acquiring the initial crack size of the pipeline according to the detection data in the pipeline.
As one or more embodiments, the pipeline basic data to be life-predicted includes: the temperature of the pipeline to be subjected to life prediction, the concentration of matrix hydrogen atoms of the pipeline to be subjected to life prediction, the Poisson's ratio of the pipeline to be subjected to life prediction and the yield strength of the pipeline to be subjected to life prediction.
It should be understood that the substrate hydrogen atom concentration of the tube refers to: atomic number percentage of hydrogen atoms in the pipeline steel.
As one or more embodiments, predicting the remaining life of the pipeline based on the initial crack size and the pipeline crack growth rate; the method comprises the following specific steps:
Figure BDA0002515307920000091
wherein a is the predicted crack length, aoFor the initial crack length, N is the number of loading cycles, and da/dN is the crack growth length per cycle, which can be determined from equation (7) or (8).
The explanation of the hydrogen embrittlement theory comprises HEDE (hydrogen enhanced condensation), HELP (hydrogen enhanced local plasticity) and HBT (hydrogen bubble the), the core idea is that hydrogen atoms gather around defects such as crack tips, dislocations and the like, the free surface energy is reduced, the local stress and the local deformation are increased, the steel material is easy to fail, and the relationship between the gathering rate of the hydrogen atoms around the defects and the crack growth rate is quantified on the basis of the hydrogen embrittlement fracture model of the theory to determine the crack growth rate; the relationship between the formation rate of the cross section and the stress intensity is firstly quantified; the influence of the loading frequency has a limit value, when the limit value is higher, the speed is increased along with the increase of the loading frequency, and when the limit value is lower, the speed is not changed along with the time; the effect of temperature and the effect of hydrogen atom concentration (i.e., the effect of pH value) were quantified.
The crack growth is predicted by a formula, and under the condition that the stress intensity of the pipeline is known to change,
analyzing the change of the original data pipeline stress intensity distribution corresponding to time: fig. 2 utilizes a virtual loading spectrogram for convenience of demonstration. There are no such large changes in stress intensity that are so regular and fluctuate in actual pipe transport.
R ratio (K)min/Kmax) The effect of loading above 0.94 on crack growth is negligible, corresponding to a stress intensity threshold of 0.594MPa m0.5。
The stress intensity range in figure 1 is greater than or equal to 3.3MPa m0.5Load cycle ofIs determined to be a non-small cycle load. For convenience of illustration, the small cycle loading is not set in fig. 1.
Decomposing the loaded wave-front into a plurality of basic modules, each module comprising a delta-K higher than 3.3MPa m0.5And a number of deltak distributions before the next such loading period occurs in the range of 0.594 to 3.3MPa m0.5Small load cycles.
Recording the screened key data into a file, processing the original data, recording the time and K in the graph 3maxAnd KminThe corresponding relationship diagram of (1). All loadings Δ K are greater than or equal to 3.3MPa m0.5The peak and trough values shown in fig. 3 are all recorded.
(1) After raw data processing, the loading frequency f needs to be recorded.
(2) Recording K for each larger load cyclemaxAnd Δ K, since only Δ K greater than or equal to 3.3MPa m is present in the initial spectrum we provide0.5So that the stress intensity values for all cycles are recorded as shown in fig. 2.
(3) Stress intensity distribution in each basic loading module is 0.594 to 3.3MPa m0.5The number (n) of small periods. There are no small periods in fig. 2, so fig. 3 does not reflect this screening condition.
(4) The temperature is provided by the pump house of the loading station.
(5) Base hydrogen atom concentration c in a pipeoObtained by experiments and simulations.
Calculating each basic load module (comprising a Δ K higher than 3.3MPa m)0.5The large loading period + n delta K is distributed between 0.594 and 3.3MPa m0.5A small loading period of (a), n is 0 in this example);
the loading frequency change in the pipeline transportation has a threshold value, and when the loading frequency is lower than the threshold value, the crack growth rate is maximized and does not change along with the loading frequency. The relationship between the crack propagation rate and the loading frequency, and for the loading in the same stress intensity range, under the condition that the frequency is greater than the threshold value, the larger the loading frequency is, the crack growth rate isThe slower. Threshold value of loading frequency is fcriticalThe function of the 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 is shown in the formula (6). Processed data, frequency and fcriticalAfter comparison, the crack propagation equations (7) and (8) based on the hydrogen atom diffusion model were entered.
Calculations were performed to obtain the crack growth rate in a single loaded module. Wherein A is a factor related to temperature, hydrogen atom concentration, material physical properties:
Figure BDA0002515307920000111
if a small period acceleration condition is considered, the small period acceleration factor can be expressed as follows:
Figure BDA0002515307920000112
wherein n is Δ K distributed between 0.594 to 3.3MPa m0.5Small number of loading cycles of coIs the hydrogen atom ratio concentration, n, of the pipe matrixcriticalAnd ccriticalFor the pipeline characteristic values, 657 and 0.01568, respectively, this acceleration factor is introduced into equation (8) in the form of a small loading period when it exists
Figure BDA0002515307920000113
Can be used to predict the true crack growth rate. Since there are no small cycles in this example, the small cycle acceleration conditions are negligible and the final crack propagation rate is shown in FIG. 4:
FIG. 4 crack length versus number of modules loaded, where there are 10 smaller modules per larger module, i.e., between the larger cycles in FIG. 1, and thus a stepwise crossover increase in crack growth occurs. By combining theoretical research, the connection between the mechanism of the hydrogen atom diffusion fracture model and the industrial reality is fully considered, the theoretical model is upgraded and reconstructed, the new mechanical model can be applied to pipeline fracture prediction in soil solution, the theoretical model is arranged into a program, the model is easier to popularize and use, and the method has important significance for pipeline integrity evaluation. The purpose of the disclosure is to establish a pipeline crack growth model based on a hydrogen embrittlement theory, and aim to quantify the dependence of the service life of a pipeline, the ambient temperature and the hydrogen atom concentration and increase the reliability and the applicability of prediction.
Example two
The embodiment provides a system for predicting the residual life of a pipeline based on hydrogen diffusion dynamics;
a system for predicting remaining life of a pipeline based on hydrogen diffusion kinetics comprises:
an acquisition module configured to: acquiring basic data of a pipeline to be subjected to life prediction;
an input module configured to: inputting the basic data of the pipeline with the service life to be predicted into the 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.
It should be noted here that the above-mentioned obtaining module, the input module and the prediction module correspond to steps S101 to S103 in the first embodiment, and the above-mentioned modules are the same as the 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 as part of a system may be implemented in a computer system such as a set of computer-executable instructions.
In the foregoing embodiments, the descriptions of the embodiments have different emphasis, and for parts that are not described in detail in a certain embodiment, reference may be made to related descriptions of other embodiments.
The proposed system can be implemented in other ways. For example, the above-described system embodiments are merely illustrative, and for example, the division of the above-described modules is merely a logical functional division, and in actual implementation, there may be other divisions, for example, multiple modules may be combined or integrated into another system, or some features may be omitted, or not executed.
EXAMPLE III
The present embodiment also provides an electronic device, including: one or more processors, one or more memories, and one or more computer programs; wherein, a processor is connected with the memory, the one or more computer programs are stored in the memory, and when the electronic device runs, the processor executes the one or more computer programs stored in the memory, so as to make the electronic device execute the method according to 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 arrays FPGA or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, and so on. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory may include both read-only memory and random access memory, and may 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 device type information.
In implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or instructions in the form of software.
The method in the first embodiment may be directly implemented by a hardware processor, or may be implemented by a combination of hardware and software modules in the processor. The software modules may be located in ram, flash, rom, prom, or eprom, registers, among other storage media as is well known in the art. The storage medium is located in a memory, and a processor reads information in the memory and completes the steps of the method in combination with hardware of the processor. To avoid repetition, it is not described in detail here.
Those of ordinary skill in the art will appreciate that the various illustrative elements, i.e., algorithm steps, described in connection with the embodiments disclosed herein may 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 implementation. 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 four
The present embodiments also provide a computer-readable storage medium for storing computer instructions, which when executed by a processor, perform the method of the first embodiment.
The above description is only a preferred embodiment of the present disclosure and is not intended to limit the present disclosure, and various modifications and changes may be made to the present disclosure by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present disclosure should be included in the protection scope of the present disclosure.

Claims (10)

1. The method for predicting the residual service life of the pipeline based on hydrogen diffusion dynamics is characterized by comprising the following steps:
acquiring basic data of a pipeline to be subjected to life prediction;
inputting the basic data of the pipeline with the service life to be predicted into the 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.
2. The method of claim 1, wherein the screened pipeline crack growth rate calculation model, the screening step comprises:
and comparing the actual loading frequency of the pipeline to be subjected to life prediction with a pipeline loading frequency threshold value, and screening out a corresponding pipeline crack growth rate calculation model.
3. The method of claim 1, wherein the selected pipeline crack growth rate calculation model is selected over a range comprising:
and constructing a corresponding pipeline crack growth rate calculation model according to the relation between the actual loading frequency of the pipeline and the threshold value of the loading frequency of the pipeline.
4. The method of claim 1, wherein the pipeline is loaded with a frequency threshold, the obtaining step comprising:
and calculating a pipeline loading frequency threshold according to the Poisson ratio of the pipeline to be subjected to life prediction, the yield strength of the pipeline to be subjected to life prediction, the volume of the hydrogen atom component of the pipeline to be subjected to life prediction, the hydrogen atom diffusivity of the pipeline to be subjected to life prediction and the temperature of the pipeline to be subjected to life prediction.
5. The method as claimed in claim 4, wherein the pipeline loading frequency threshold is calculated according to the Poisson's ratio of the pipeline to be life-predicted, the yield strength of the pipeline to be life-predicted, the volume of the hydrogen atom component of the pipeline to be life-predicted, the hydrogen atom diffusivity of the pipeline to be life-predicted and the temperature of the pipeline to be life-predicted; the method comprises the following specific steps:
Figure FDA0002515307910000021
wherein v is Poisson's ratio, Ω is hydrogen atom component volume, D is hydrogen atom diffusivity, KmaxAnd KminMaximum and minimum stress intensity values for a single cyclic load, rpFor the size of the plastic deformation zone, ReqSize of hydrogen-donating region, k, for plastic regionBBoltzmann constant, T is temperature.
6. The method of claim 3, wherein the corresponding pipeline crack growth rate calculation model is constructed based on a relationship between an actual pipeline loading frequency and a pipeline loading frequency threshold; the method comprises the following specific steps:
Figure FDA0002515307910000022
Figure FDA0002515307910000023
wherein R is KminAnd KmaxA is an environmental factor and Δ K is KmaxAnd KminN is a material-related property constant, and is obtained by fitting.
7. The method as claimed in claim 6, 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: when the actual loading frequency of the pipeline to be subjected to life prediction is larger than the pipeline loading frequency threshold value, selecting a formula (7); and when the actual loading frequency of the pipeline to be subjected to life prediction is less than or equal to the pipeline loading frequency threshold value, selecting the formula (8).
8. The system for predicting the residual life of the pipeline based on hydrogen diffusion dynamics is characterized by comprising the following steps:
an acquisition module configured to: acquiring basic data of a pipeline to be subjected to life prediction;
an input module configured to: inputting the basic data of the pipeline with the service life to be predicted into the 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.
9. An electronic device, comprising: one or more processors, one or more memories, and one or more computer programs; wherein a processor is connected 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 any of the preceding claims 1-7.
10. A computer-readable storage medium storing computer instructions which, when executed by a processor, perform the method of any one of claims 1 to 7.
CN202010474166.2A 2020-05-29 2020-05-29 Method and system for predicting residual life of pipeline based on hydrogen diffusion dynamics Active CN111625963B (en)

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 true CN111625963A (en) 2020-09-04
CN111625963B 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112666066A (en) * 2020-12-15 2021-04-16 中国石油大学(华东) Pipeline hydrogen embrittlement temperature threshold prediction method based on hydrogen diffusion dynamics and application

Citations (3)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112666066A (en) * 2020-12-15 2021-04-16 中国石油大学(华东) Pipeline hydrogen embrittlement temperature threshold prediction method based on hydrogen diffusion dynamics and application

Also Published As

Publication number Publication date
CN111625963B (en) 2024-05-28

Similar Documents

Publication Publication Date Title
Ai et al. Probabilistic modeling of fatigue life distribution and size effect of components with random defects
CN106886663B (en) Method and device for predicting bending fatigue life of gear
Susmel et al. An elasto-plastic reformulation of the theory of critical distances to estimate lifetime of notched components failing in the low/medium-cycle fatigue regime
CN107273649A (en) A kind of Forecasting Methodology of fragile material failure probability under high-temerature creep state
Zhao et al. Collapse pressure of randomly corroded spherical shell
Recupero et al. Experimental tests on corroded prestressed concrete beams subjected to transverse load
Zhangabay et al. Analysis of a stress-strain state of a cylindrical tank wall vertical field joint zone
CN111625963A (en) Pipeline residual life prediction method and system based on hydrogen diffusion dynamics
Maraveas Concrete silos: Failures, design issues and repair/strengthening methods
Yu et al. Constitutive analysis of pressure-insensitive metals under axisymmetric tensile loading: A stress triaxiality-dependent plasticity damage model
Scales et al. Accurate estimation of yield strength and ultimate tensile strength through instrumented indentation testing and chemical composition testing
Ma et al. Computational investigation of multi-axial damage modeling for porous sintered metals with experimental verification
Zhao et al. Effect of pressure sampling methods on pipeline integrity analysis
Zhao et al. Fatigue Properties Estimation and Life Prediction for Steels under Axial, Torsional, and In‐Phase Loading
Chattopadhyay et al. Development of new correlations for improved integrity assessment of pipes and pipe bends
CN111797456B (en) Prediction method for mechanical property degradation rule of post-rust steel bar
Tian et al. Experimental study on bond performance and damage detection of corroded reinforced concrete specimens
Karlén et al. A study of the volume effect and scatter at the fatigue limit–experiments and computations for a new specimen with separated notches
Zheng et al. Dynamics of necking and fracture in ductile porous materials
Andreikiv et al. Computational models of fatigue cracks growth in metallic materials under the action of force and physicochemical factors
Mahmud et al. On the need to adopt strain-based probabilistic approach in predicting fatigue life
Wang et al. Effect of Concrete Age and Creep on the Behavior of Concrete‐Filled Steel Tube Columns
Savkin et al. Effect of Overload on the Near-Threshold Fatigue Crack Growth Rate in a 2024-T3 Aluminum Alloy: II. Fatigue Crack Growth Simulation for Calculating the Fatigue Life under Alternating Loading
Lei et al. A new fatigue damage model for pavement concrete beams bearing multi-level bending loads
Pandey et al. Fatigue crack growth model for constant amplitude loading

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