CN115238557A - Method for evaluating hydrogen loss life of shock tube body - Google Patents
Method for evaluating hydrogen loss life of shock tube body Download PDFInfo
- Publication number
- CN115238557A CN115238557A CN202210911522.1A CN202210911522A CN115238557A CN 115238557 A CN115238557 A CN 115238557A CN 202210911522 A CN202210911522 A CN 202210911522A CN 115238557 A CN115238557 A CN 115238557A
- Authority
- CN
- China
- Prior art keywords
- tube body
- current
- pipe body
- hydrogen
- hydrogen concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 title claims abstract description 60
- 239000001257 hydrogen Substances 0.000 title claims abstract description 60
- 229910052739 hydrogen Inorganic materials 0.000 title claims abstract description 60
- 230000035939 shock Effects 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 12
- 239000000463 material Substances 0.000 claims description 9
- 150000002431 hydrogen Chemical class 0.000 claims description 8
- 238000011156 evaluation Methods 0.000 abstract description 10
- 238000009792 diffusion process Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
The invention belongs to the technical field of hypersonic test equipment, and discloses an online evaluation method for hydrogen loss life of a shock tube body. According to the method for the online evaluation of the hydrogen loss life of the shock wave tube body, the radial strain value and the radial hydrogen concentration distribution of the current tube body are calculated online by collecting the service working condition parameters of the tube body; obtaining a breaking strain value under the current tube body hydrogen concentration by utilizing the influence of the hydrogen concentration fitted by the test data on the elongation; obtaining a danger coefficient according to the ratio of the current strain value to the fracture strain value; and predicting the service life of the pipe body through the historical risk coefficient. According to the method for evaluating the hydrogen loss life of the shock tube body on line, the strength of the tube body is evaluated through the hydrogen concentration, and the hydrogen loss degree of the tube body is evaluated on line; evaluating the hydrogen loss life of the pipe body on line through the time fitting relation of the historical risk coefficients; the economic efficiency and the timeliness are good, and the method is suitable for engineering popularization.
Description
Technical Field
The invention belongs to the technical field of hypersonic test equipment, and particularly relates to an online evaluation method for hydrogen loss life of a shock tube body.
Background
The shock tunnel is a pulse type tunnel which utilizes shock waves to compress test gas and generates hypersonic test airflow through the steady expansion of a spray pipe. The hydrogen driving mode is the main driving mode of the shock tunnel, and stronger shock waves can be generated by using hydrogen driving to obtain higher test parameters. The heating of the hydrogen gas also allows for greater driving ability.
When the shock wave tube body is actually used, the shock wave tube body directly contacts hydrogen, hydrogen damage such as hydrogen embrittlement and hydrogen corrosion can be generated on the material, the service life of the tube body is adversely affected, but the shock wave tube body cannot be directly judged on the appearance, and great hidden danger is generated on the service safety; in addition, the shock tube body has high pressure and large temperature range under the service condition, and the damage of hydrogen leakage is huge, so that the service life of the tube body cannot be judged, and the tube body has a large safety risk.
At present, the development of an online evaluation method for hydrogen loss life of a shock tube body for a hydrogen-driven shock tunnel is urgently needed.
Disclosure of Invention
The invention aims to solve the technical problem of providing an online evaluation method for the hydrogen loss life of a shock tube body.
The invention discloses an online evaluation method for hydrogen loss life of a shock tube body, which comprises the following steps:
the method comprises the following steps: the outer layer temperature, the inner layer temperature and the inner layer pressure of the pipe body are collected by a sensor arranged on the pipe body;
step two: establishing a one-dimensional axisymmetric finite element model;
step three: introducing the collected outer layer temperature, inner layer temperature and inner layer pressure of the pipe body into a one-dimensional axisymmetric finite element model, and calculating the radial stress distribution of the current pipe body to obtain a radial strain value of the current pipe body;
step four: calculating the current saturated hydrogen concentration of the inner layer of the tube body according to the acquired temperature of the inner layer of the tube body and the pressure of the inner layer of the tube body by using the Sieverts' law;
step five: introducing the current inner layer saturated hydrogen concentration of the tube body in a one-dimensional axisymmetric finite element model, and calculating the radial hydrogen concentration distribution of the current tube body according to the historical hydrogen concentration of the tube body;
step six: introducing test data of materials with different hydrogen contents, fitting a radial elongation threshold of the pipe body under the current hydrogen concentration distribution condition, and further calculating a fracture strain value under the current hydrogen concentration distribution condition;
step seven: calculating the ratio of the radial strain value of the current pipe body to the fracture strain value of the current pipe body to obtain a current risk coefficient;
step eight: and (5) predicting the service life of the pipe body under the current risk coefficient by utilizing the historical risk coefficient fitting.
The method for evaluating the hydrogen loss life of the shock tube body on line has the following advantages:
a. only a small amount of conventional parameters can be collected to realize the evaluation of the hydrogen loss life of the pipe body;
b. the hydrogen damage degree of the pipe body material can be evaluated on line, and the pipe body material does not need to be sampled and detected.
c. The service life of the pipe body can be predicted on line, and the maintenance management of the pipe body is guided.
The on-line evaluation method for the hydrogen loss life of the shock tube body is established on the basis of the following four parts: 1. collecting service working condition data of a pipe body; 2. establishing a calculation model of pipe strain distribution and hydrogen concentration distribution; 3. obtaining test data of the influence of the hydrogen concentration on the elongation of the material; 4. and obtaining a historical risk coefficient.
According to the method for the on-line evaluation of the hydrogen loss life of the shock wave tube body, the radial strain value of the current tube body and the radial hydrogen concentration distribution of the current tube body are calculated on line by acquiring the service working condition parameters of the tube body; obtaining a current fracture strain value of the tube body under the radial hydrogen concentration distribution of the current tube body by utilizing the influence of the hydrogen concentration fitted by the test data on the elongation; obtaining a current danger coefficient according to the ratio of the radial strain value of the current pipe body to the fracture strain value of the current pipe body; and finally, predicting the service life of the pipe body under the current risk coefficient through the historical risk coefficient.
In short, the method for evaluating the hydrogen damage life of the shock tube body on line evaluates the strength of the tube body through the hydrogen concentration, and evaluates the hydrogen damage degree of the tube body on line; evaluating the hydrogen loss life of the pipe body on line through the time fitting relation of the historical risk coefficients; the economy and the timeliness are good, and the method is suitable for engineering popularization.
Drawings
FIG. 1 is a flow chart of the method for online evaluation of hydrogen loss life of a shock tube body according to the present invention.
Detailed Description
The present invention is described in detail below with reference to the drawings and examples.
Example 1
As shown in fig. 1, the method for online evaluating the hydrogen loss life of the shock tube body in the embodiment includes the following steps:
the method comprises the following steps: the temperature T of the outer layer of the pipe body is acquired by a sensor arranged on the pipe body 0 Inner layer temperature T of pipe body i And the inner layer pressure P of the pipe body;
step two: establishing a one-dimensional axisymmetric finite element model;
step three: introducing the collected outer layer temperature T of the pipe body in a one-dimensional axisymmetric finite element model 0 Inner layer temperature T of pipe body i And calculating the radial stress distribution of the current pipe body to obtain the radial strain value epsilon of the current pipe body r ；
The specific process is that the temperature T of the outer layer of the pipe body is used 0 And inner layer temperature T of pipe body i As a boundary condition, by the equation for heat conduction
Performing discrete solution to obtain the thermal strain value epsilon of the current pipe body th Is shown as
Wherein,
in order to output heat flux, lambda is a heat conduction coefficient, alpha is a linear expansion coefficient, r is the radius of the pipe body, and T is temperature;
the outer layer pressure is 0, the inner layer pressure P of the pipe body is taken as a boundary condition, and the radial strain value epsilon of the current pipe body is obtained by solving an equation set of force balance, deformation coordination and a physical equation r And hoop strain value ε θ (ii) a The system of equations for force balance, deformation coordination and physical equations is:
wherein σ r 、σ θ Respectively radial stress, hoop stress,. Epsilon r 、ε θ Respectively radial strain and hoop strain, sigma and epsilon are equivalent stress and equivalent strain, and K, n is a material constant.
Step four: by using Sieverts' law, the collected temperature T of the inner layer of the tube body i And calculating the current inner saturated hydrogen concentration C of the pipe body according to the inner pressure P of the pipe body H ；
C H The calculation formula is as follows:
a is a constant, R is a gas constant, Δ H is the heat of solution;
step five: introducing the saturated hydrogen concentration C of the inner layer of the current tube body into a one-dimensional axisymmetric finite element model H And calculating the radial hydrogen concentration distribution C of the current tube body according to the historical hydrogen concentration of the tube body H (r)；
The specific method comprises the following steps: solving a diffusion equation by taking the radial hydrogen concentration distribution of the current pipe body calculated at the previous acquisition moment as an initial condition
Wherein the boundary condition of the inner layer is the saturated hydrogen concentration C of the inner layer of the current tube body at the current moment H (ii) a Where C is the hydrogen concentration, τ is the time, and D is the diffusion coefficient.
Step six: introducing test data of materials with different hydrogen contents, fitting the radial elongation threshold of the tube body under the current hydrogen concentration distribution condition, and further calculating the fracture strain value epsilon under the current hydrogen concentration distribution condition f ；
Introducing the test data of the influence of the hydrogen concentration on the elongation of the material into the radial hydrogen concentration distribution of the current tube body, and calculating the fracture strain value epsilon of the current tube body f ＝ln(1+A * (C) Wherein A) is * (C) The radial elongation of the current tube at the corresponding hydrogen concentration obtained by fitting;
step seven: respectively calculating the radial strain value epsilon of the current pipe body r Strain at break epsilon of current tube body f Obtaining the current risk coefficient;
ζ r ＝(ε r +ε th )/ε f
therein, ζ r Is the current radial risk factor;
step eight: predicting the service life of the pipe body under the current risk coefficient by utilizing the historical risk coefficient fitting;
and defining a danger threshold value through the fitting relation between the historical danger coefficient and the time, fitting the time required by the current danger coefficient to reach the threshold value, and predicting the service life of the pipe body.
Although the embodiments of the present invention have been disclosed above, it is not limited to the applications listed in the description and the embodiments, but it can be applied to various fields suitable for the present invention. Additional modifications and refinements will readily occur to those skilled in the art without departing from the principles of the present invention, and the present invention is not limited to the specific details and illustrations shown and described herein.
Claims (1)
1. The method for evaluating the hydrogen loss life of the shock tube body on line is characterized by comprising the following steps of:
the method comprises the following steps: the outer layer temperature, the inner layer temperature and the inner layer pressure of the pipe body are collected by a sensor arranged on the pipe body;
step two: establishing a one-dimensional axisymmetric finite element model;
step three: introducing the collected outer layer temperature, inner layer temperature and inner layer pressure of the pipe body into a one-dimensional axisymmetric finite element model, and calculating the radial stress distribution of the current pipe body to obtain the radial strain value of the current pipe body;
step four: calculating the current saturated hydrogen concentration of the inner layer of the tube body according to the acquired temperature of the inner layer of the tube body and the pressure of the inner layer of the tube body by using the Sieverts' law;
step five: introducing the current inner layer saturated hydrogen concentration of the tube body in a one-dimensional axisymmetric finite element model, and calculating the radial hydrogen concentration distribution of the current tube body according to the historical hydrogen concentration of the tube body;
step six: introducing test data of materials with different hydrogen contents, fitting the radial elongation threshold of the tube body under the current hydrogen concentration distribution condition, and further calculating a breaking strain value under the current hydrogen concentration distribution condition;
step seven: calculating the ratio of the radial strain value of the current pipe body to the fracture strain value of the current pipe body to obtain a current risk coefficient;
step eight: and (5) predicting the service life of the pipe body under the current risk coefficient by utilizing the historical risk coefficient fitting.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210911522.1A CN115238557A (en) | 2022-07-28 | 2022-07-28 | Method for evaluating hydrogen loss life of shock tube body |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210911522.1A CN115238557A (en) | 2022-07-28 | 2022-07-28 | Method for evaluating hydrogen loss life of shock tube body |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115238557A true CN115238557A (en) | 2022-10-25 |
Family
ID=83677289
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210911522.1A Pending CN115238557A (en) | 2022-07-28 | 2022-07-28 | Method for evaluating hydrogen loss life of shock tube body |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115238557A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116380367A (en) * | 2023-06-06 | 2023-07-04 | 中国空气动力研究与发展中心超高速空气动力研究所 | Hydrogen leakage monitoring device and monitoring method for shock tube of high-pressure hydrogen driver |
-
2022
- 2022-07-28 CN CN202210911522.1A patent/CN115238557A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116380367A (en) * | 2023-06-06 | 2023-07-04 | 中国空气动力研究与发展中心超高速空气动力研究所 | Hydrogen leakage monitoring device and monitoring method for shock tube of high-pressure hydrogen driver |
| CN116380367B (en) * | 2023-06-06 | 2023-08-01 | 中国空气动力研究与发展中心超高速空气动力研究所 | Hydrogen leakage monitoring device and monitoring method for shock tube of high-pressure hydrogen driver |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115238557A (en) | Method for evaluating hydrogen loss life of shock tube body | |
| CN106371480B (en) | A kind of catalytic oxidation treatment ammonia-contaminated gas electricity consumption heating and temperature control system and method | |
| WO2023087443A1 (en) | Digital twin-based pressure vessel safety evaluation and risk warning method | |
| CN102539325A (en) | Method for monitoring structural damage of coke tower based on strain monitoring | |
| CN111090957B (en) | Dangerous point stress-strain calculation method for high-temperature structure | |
| CN102495094A (en) | Online calculation method of enthalpy-entropy values of water and water vapor | |
| CN106680372B (en) | Correction method for evaluating grain size of metal material by stress influence ultrasonic wave | |
| CN101915786B (en) | System and method for monitoring combustible gas | |
| CN111312414B (en) | Fatigue evaluation method of nuclear power plant fatigue monitoring and life evaluation system | |
| CN103207567B (en) | Low-false-alarm-rate improved principal component analysis process monitoring method and system | |
| Fukumoto et al. | Investigation of High-Temperature Oxidation in Steam for Ni–Al Alloys Using the Combination of a Hydrogen Sensor and an Oxygen Pump–Sensor | |
| CN112818579A (en) | Multi-parameter coupling-based efficient detection method for thermal stress field of equipment | |
| EP0165675B2 (en) | Apparatus for measuring thermal stress of pressure-tight tube | |
| CN116949129A (en) | Nucleic acid annealing process | |
| CN111536430B (en) | Method for three-dimensionally positioning pipeline leakage monitoring based on humidity | |
| CN110705135A (en) | High-cycle fatigue life estimation method for mixing three-way pipeline | |
| Moorthy et al. | Diagnostics of direct CT-PT contact of the coolant channels of PHWRs | |
| CN105606159B (en) | Heat source tower chemical heat pump with life-span real-time estimate function | |
| Yang et al. | On the Interrupted Creep Test Under Low Stress Levels for the Power Law Parameter Extraction | |
| JPS57157138A (en) | Foretelling method and device for device life | |
| Pu et al. | The refinement of dynamic deformation equation for AISI1215 steel using the Zener-Hollomon parameter | |
| CN213517992U (en) | Conversion section safety control system for porous ammonium nitrate production | |
| CN117272793A (en) | Real-time monitoring system for service life of high-temperature part of boiler | |
| CN114279607B (en) | Cable joint interface pressure monitoring method and device based on acoustic elastic effect | |
| Wahlen et al. | Computer controlled measurement and analytical modelling of flow stresses during hot deformation of the cooper alloy CuZn42Mn2 |
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 |