CN106295021A - A kind of Forecasting Methodology facing the strain of hydrogen heavy wall cylindrical shell elastic stress - Google Patents
A kind of Forecasting Methodology facing the strain of hydrogen heavy wall cylindrical shell elastic stress Download PDFInfo
- Publication number
- CN106295021A CN106295021A CN201610670679.4A CN201610670679A CN106295021A CN 106295021 A CN106295021 A CN 106295021A CN 201610670679 A CN201610670679 A CN 201610670679A CN 106295021 A CN106295021 A CN 106295021A
- Authority
- CN
- China
- Prior art keywords
- theta
- sigma
- cylindrical shell
- hydrogen
- strain
- 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
Links
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
- 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/06—Power analysis or power optimisation
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)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
The invention discloses a kind of Forecasting Methodology facing the strain of hydrogen heavy wall cylindrical shell elastic stress, the HELP of the method application hydrogen embrittlement is theoretical, with the elastic stress strain under hydrogen environment of the cylindrical shell directly prediction cylindrical shell of the elastic stress strain under atmospheric environment.Compared to the method for structural stress strain under existing prediction hydrogen environment, the present invention has the simple advantage of principle.Straining Forecasting Methodology compared to existing cylindrical shell elastic stress, the impact of cylindrical shell ess-strain is taken into account, engineering design is had some reference value by the present invention by hydrogen environment.
Description
Technical field
The present invention relates to face the elastic response prediction field of hydrogen carrying structure, specifically based on hydrogen embrittlement HELP theoretical prediction
The elastic stress strain of heavy wall cylindrical shell.
Background technology
Cylindrical shell is the crucial pressure restraining element in hydrogen system during hydrogen storage and defeated hydrogen, determines that heavy wall cylindrical shell is at hydrogen environment
Under the elastic response problem that must take into when being the design hydrogen-contacting equipment such as hydrogen container, hydrogenation reactor.Due to this structure of hydrogen damage
The complexity of relation, the method for existing prediction material mechanical response under hydrogen environment is to the competency profiling of research worker very
Height, simultaneously need to a large amount of programing work, it is difficult to promote in engineering.The most do not occur that being specifically designed to prediction faces hydrogen cylindrical shell stress
The method of strain.It is proposed that a kind of method that hydrogen cylindrical shell ess-strain is faced in simple prediction has engineering significance.The present invention
The hydrogen of hydrogen embrittlement is promoted, and plasticity localization theoretical (HELP) is applied to heavy wall cylindrical shell, it is proposed that a kind of with atmospheric environment cylinder bullet
Property ess-strain prediction face hydrogen cylindrical shell elastic stress strain method.
Summary of the invention
Present invention aims to the deficiencies in the prior art, it is provided that one faces the strain of hydrogen heavy wall cylindrical shell elastic stress
Forecasting Methodology, the method predicts the elasticity facing hydrogen cylindrical shell based on the cylindrical shell elastic stress strain in classical atmospheric environment
Ess-strain.
It is an object of the invention to be achieved through the following technical solutions: one faces the strain of hydrogen heavy wall cylindrical shell elastic stress
Forecasting Methodology, the method comprises the following steps:
Step 1: calculate circumferential stress σ of heavy wall cylindrical shell according to below equationθAnd radial stress σr:
Wherein a is cylindrical shell inside radius, and b is cylindrical shell outer radius, and p is inner pressuring load, and r is to need to ask for ess-strain
Position.
Step 2: solve below equation group, calculating hydrogen concentration c:
Wherein, E is the elastic modelling quantity of material, and ν is the Poisson's ratio of material, VHRepresent hydrogen partial molar volume in mother metal, VM
For the molal volume of mother metal, α is the hydrogen trap number of per unit lattice, the interstitial void binding site number of β per unit lattice, NL=
NA/VMFor the metallic atom quantity of per unit volume, NAFor Avogadro's number, NTRepresent the trap number of per unit volume.
R is ideal gas constant, and T is Kelvin's thermometric scale, c0For the hydrogen concentration in cylindrical shell during no-load, WBHydrogen trap for material combines
Energy.
Step 3: calculate axial stress σz
Step 4: calculate circumferential strain εθAnd radial strain εr
Further, described step 2 uses Newton Algorithm hydrogen concentration c, specifically includes following sub-step:
Step 201: hydrogen concentration initial value c=c is set0;
Step 202: according to current hydrogen concentration c, calculates parameters according to the following formula:
Step 203: calculate following functional value g:
Step 204: judge g value size, if | g | >=εerr, then execution step 205 is to step 206, otherwise terminates to calculate,
To hydrogen concentration c, wherein εerrFor convergence tolorence, desirable εerr=10-6;
Step 205: calculate hydrogen concentration c according to the following formula1
Wherein
Step 206: make c=c1, return step 202.
The invention have the advantages that and use the method solving Nonlinear System of Equations, according to the cylindrical shell being easy to get at air
The elastic stress strain of cylindrical shell under elastic stress strain in environment directly prediction hydrogen environment, it is not necessary to write the hydrogen loss of complexity
Hindering the finite element program of material constitutive model, application threshold is relatively low.
Accompanying drawing explanation
Fig. 1 is the objective for implementation sketch of the present invention;
Fig. 2 is the distribution on different radii of the present example calculated stress;
Fig. 3 is present example calculated strain distribution on different radii.
Detailed description of the invention
Below with the example shown in Fig. 1 and Biao 1 as objective for implementation, the invention will be further described.
Example shown in Fig. 1 is the heavy wall cylindrical shell of a two ends constraint, and its inside radius and outer radius are respectively a and b, hold
Being subject to constant internal pressure load p, under unstress state, concentration is c0Hydrogen be uniformly distributed in cylindrical shell.The present invention can predict
The cylindrical shell being in hydrogen environment stress and strain under intrinsic pressure p effect.
Material parameter that table 1 example is used and geometric parameter
The inventive method to realize process as follows:
Step 1: choose a radius r (a≤r≤b), such as, take r=a, by the inside radius a in table 1, outer radius b and
Inner pressuring load p brings below equation into and calculates circumferential stress σ of heavy wall cylindrical shellθAnd radial stress σr:
Step 2: solve below equation group, calculating hydrogen concentration c:
The most each parameter is all shown in Table 1, and step 2 specifically includes following sub-step:
Step 201: hydrogen concentration initial value c=c is set0;
Step 202: according to current hydrogen concentration c bring into following formula calculate parameters:
Step 203: parameter step 202 obtained brings following functional value into, calculating functional value g:
Step 204: judge g value size, if | g | >=εerr, then execution step 205 is to step 206, otherwise terminates to calculate,
To hydrogen concentration c, wherein convergence tolorence takes εerr=10-6;
Step 205: calculate hydrogen concentration c according to the following formula1
Wherein
Required each parameter is calculated by step 202;
Step 206: make c=c1, return step 202.
Step 3: hydrogen concentration c step 2 obtained is brought following formula into and calculated axial stress σz
Step 4: circumferential stress σ that step 1 and step 3 are obtainedθ, radial stress σrWith axial stress σzBring following formula meter into
Calculate circumferential strain εθAnd radial strain εr
Different r values is used to repeat the stress and strain distribution that step 1 to step 4 can be calculated on whole cross section,
Result as shown in Figure 2 and Figure 3, additionally can be seen that hydrogen concentration distribution in cylindrical shell is uniform.
Claims (2)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610670679.4A CN106295021B (en) | 2016-08-15 | 2016-08-15 | A kind of Forecasting Methodology for facing the strain of hydrogen heavy wall cylindrical shell elastic stress |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610670679.4A CN106295021B (en) | 2016-08-15 | 2016-08-15 | A kind of Forecasting Methodology for facing the strain of hydrogen heavy wall cylindrical shell elastic stress |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106295021A true CN106295021A (en) | 2017-01-04 |
CN106295021B CN106295021B (en) | 2018-05-29 |
Family
ID=57670685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610670679.4A Expired - Fee Related CN106295021B (en) | 2016-08-15 | 2016-08-15 | A kind of Forecasting Methodology for facing the strain of hydrogen heavy wall cylindrical shell elastic stress |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106295021B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20190137337A (en) * | 2018-06-01 | 2019-12-11 | 경희대학교 산학협력단 | A method for determining average radial stress and pressure vs. volume relationship of a compressible material |
CN110727997A (en) * | 2019-09-20 | 2020-01-24 | 湖北省工业建筑集团有限公司 | Method for calculating stability of metal cylindrical shell |
CN115238557A (en) * | 2022-07-28 | 2022-10-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | An evaluation method for hydrogen loss life of shock tube body |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002342694A (en) * | 2001-05-17 | 2002-11-29 | Araco Corp | Method and device for structural analysis of elastic material |
CN101323247A (en) * | 2008-07-11 | 2008-12-17 | 清华大学 | Concentration dilution device for hydrogen gas emitted by fuel cell buses |
CN102141466A (en) * | 2010-12-21 | 2011-08-03 | 浙江大学 | Boundary simulation device and method for thin-walled cylindrical shell structure experiment |
CN103154703A (en) * | 2010-10-05 | 2013-06-12 | 株式会社普利司通 | Method for predicting elastic response performance of rubber product, method for design, and device for predicting elastic response performance |
-
2016
- 2016-08-15 CN CN201610670679.4A patent/CN106295021B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002342694A (en) * | 2001-05-17 | 2002-11-29 | Araco Corp | Method and device for structural analysis of elastic material |
CN101323247A (en) * | 2008-07-11 | 2008-12-17 | 清华大学 | Concentration dilution device for hydrogen gas emitted by fuel cell buses |
CN103154703A (en) * | 2010-10-05 | 2013-06-12 | 株式会社普利司通 | Method for predicting elastic response performance of rubber product, method for design, and device for predicting elastic response performance |
CN102141466A (en) * | 2010-12-21 | 2011-08-03 | 浙江大学 | Boundary simulation device and method for thin-walled cylindrical shell structure experiment |
Non-Patent Citations (1)
Title |
---|
陈志平 等: "多层不等厚组合圆柱壳屈曲数值模拟分析", 《浙江大学学报(工学版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20190137337A (en) * | 2018-06-01 | 2019-12-11 | 경희대학교 산학협력단 | A method for determining average radial stress and pressure vs. volume relationship of a compressible material |
KR102085134B1 (en) | 2018-06-01 | 2020-03-05 | 경희대학교 산학협력단 | A method for determining average radial stress and pressure vs. volume relationship of a compressible material |
CN110727997A (en) * | 2019-09-20 | 2020-01-24 | 湖北省工业建筑集团有限公司 | Method for calculating stability of metal cylindrical shell |
CN115238557A (en) * | 2022-07-28 | 2022-10-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | An evaluation method for hydrogen loss life of shock tube body |
CN115238557B (en) * | 2022-07-28 | 2025-03-18 | 中国空气动力研究与发展中心超高速空气动力研究所 | A method for evaluating hydrogen damage life of shock tube |
Also Published As
Publication number | Publication date |
---|---|
CN106295021B (en) | 2018-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | A time-variant reliability analysis method based on stochastic process discretization | |
CN102567567A (en) | Finite element analysis based pipeline random-vibration fatigue life analyzing method | |
CN106295021A (en) | A kind of Forecasting Methodology facing the strain of hydrogen heavy wall cylindrical shell elastic stress | |
CN105005660A (en) | Stress calculation method for non-linear flexible member close to practical running state | |
CN103745114A (en) | Method for computing stress relaxation numerical values and resilience of titanium alloy | |
Ekh et al. | Models for cyclic ratchetting plasticity—integration and calibration | |
Moghanlou et al. | Assessment of the pitting corrosion degradation lifetime: a case study of boiler tubes | |
Wang et al. | A continuum damage mechanics-based viscoplastic model of adapted complexity for high-temperature creep–fatigue loading | |
CN114970106A (en) | Method and system for predicting radiation hardening based on microstructure | |
Thiagarajan et al. | Fracture simulation using an elasto-viscoplastic virtual internal bond model with finite elements | |
CN106372273B (en) | One kind faces hydrogen heavy wall cylindrical shell elastic limit loading prediction method | |
Ravi Kiran et al. | Fatigue-ratcheting behavior of 6 in pressurized carbon steel piping systems under seismic load: experiments and analysis | |
CN106874579A (en) | GB150-1998 analysis of Ultimate method based on yield strength and the double key points of tensile strength | |
Funke et al. | Modeling the tensile strength and crack length of wire-sawn silicon wafers | |
Zhang et al. | Identification of material parameters for high strength steel under impact loading | |
Hwang et al. | Evaluation of flow stresses of tubular materials considering anisotropic effects by hydraulic bulge tests | |
He et al. | Damage-plastic constitutive model of thin-walled circular steel tubes for space structures | |
Goyal et al. | Cyclic plasticity and fatigue-ratcheting behavior of SS304LN stainless steel material | |
Kontermann et al. | On the Evaluation and Consideration of Fracture Mechanical Notch Support Within a Creep-Fatigue Lifetime Assessment | |
Nguyen et al. | Modeling of microstructure effects on the mechanical behavior of ultrafine-grained nickels processed by severe plastic deformation by crystal plasticity finite element model | |
Liu et al. | The whole-life ratcheting behavior of internal pressurized elbow subjected to force and displacement cycling with a damage-coupled kinematic hardening model | |
CN105205342A (en) | Method of determining limit load of high temperature pressure pipeline | |
CN117497102B (en) | Prediction method and system for helium bubble evolution of metal material irradiation | |
Wijeyeratne et al. | A Comparative Study of Crystal Viscoplastic Modeling of Directionally Solidified Nickel-Base Superalloys | |
CN107766659A (en) | It is a kind of suitable for the elastic constitutive model model of rubber type of material and its application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180529 |
|
CF01 | Termination of patent right due to non-payment of annual fee |