CN1329689C - Fatigue life safety predicting method for pressure container - Google Patents
Fatigue life safety predicting method for pressure container Download PDFInfo
- Publication number
- CN1329689C CN1329689C CNB2004100677460A CN200410067746A CN1329689C CN 1329689 C CN1329689 C CN 1329689C CN B2004100677460 A CNB2004100677460 A CN B2004100677460A CN 200410067746 A CN200410067746 A CN 200410067746A CN 1329689 C CN1329689 C CN 1329689C
- Authority
- CN
- China
- Prior art keywords
- delta
- partiald
- crack
- increment
- pressure container
- 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.)
- Expired - Fee Related
Links
Abstract
The present invention discloses a fatigue life safe prediction method of pressure containers, which comprises the following steps: crack cleft area increment is obtained through inputting initial defect information and node three-dimensional address information for further obtaining crack expansion increment; the crack expansion increment and the original crack depth are added for obtaining crack depth; the crack depth is compared with tolerance size stored in computers for controlling pressure containers. The present invention carris out real-time control to pressure containers through computers and has important application value for preventing sudden destruction accidents of pressure containers in use, predicting the residual service life of pressure containers in use and ensuring the safe use of pressure containers.
Description
Technical field
The present invention relates to a kind of safety detecting method of pressurized container, particularly fatigue life safety predicting method for pressure container.
Background technique
Pressurized container is the key equipment in the development of the national economy, be again a kind of extraordinary bearing device simultaneously with explosion risk, it is bearing the pressure of high temperature, low temperature, inflammable, explosive, severe toxicity or aggressive medium, will cause irremediable catastrophic failure in case take place to destroy and leak.In the pressurized container malicious event, have greatly because the tired expansion of crackle causes, thereby can learn the propagation of crackle in advance, just can take appropriate measures to pressurized container, avoid the generation of accidents such as exploding.The fatigue life prediction of pressurized container is the emphasis problem of domestic and international breaking power educational circles research always.The method that generally adopts is a measurement method at present, but this method is subjected to the restriction of various objective condition, can't carry out on-the-spot test to inservice pressure vessel, and manpower and materials consumption is big, validity can not guarantee.A kind of numerical simulation technology of the pressurized container fatigue life prediction based on stress intensity factor width of cloth Δ K also has reported in literature at present, and has necessarily been used in actual engineering.But Δ K only is confined to the breakage problem of linear elasticity and small scale yielding, can not be suitable for for a large amount of large scale yielding that exists in the engineering reality and complete yield breakage problem, should adopt in theory comparatively perfect cyclic J integration Δ J as fracture parameter.
In the practice, people wish by computer pressurized container to be detected incessantly, make immediately closing or keeping operational order such as use, come the pilot pressure container by instruction again, guarantee being perfectly safe of container.
Summary of the invention
At the technical problem that exists in the prior art, the invention provides and a kind ofly come the pilot pressure container by computer, can prevent the sudden malicious event of inservice pressure vessel generation, predict its fatigue surplus life, guarantee the fatigue life safety predicting method for pressure container that pressure vessel safety uses.
The present invention is to realize by such technological scheme: a kind of fatigue life safety predicting method for pressure container is provided, may further comprise the steps for reaching above purpose:
1), is used to import the input step of initial imperfection information;
2), be used to import the input step of the three-dimensional address information of node;
3), be used to obtain the step of crackle area increment, adopt
Δ x wherein
n, Δ y
n, Δ z
nThe displacement of expression node n on X, Y, Z direction, h
n(s
i) the expression shape function, J (s
i) the expression Jacobian matrix, the s at node 1,2,3 places
iThe expression value is respectively-1,0 ,+1.(seeing accompanying drawing 1)
4), be used to obtain the step of elastoplasticity cyclic J integration Δ J, adopt
A wherein
nThe increment of the crackle area that the empty expansion of expression crackle causes, σ
IjThe expression stress tensor, u
iThe expression displacement vector, W represents strain energy density, δ represents the Kronecker tensor, Δ x
iThe mapping function of representing transformation of coordinates between two independent coordinate systems, V represent to be subjected in the cracks in body the empty volumetric region that influences of expanding of crackle.
5), be used to obtain the step of crack propagation increment, adopt
N wherein
i, N
I+1, a
i, a
I+1Number of times and crack length when the expression crack propagation arrives i, i+1, c, n represent material constant, m represents to calculate segments.
6), be used to obtain the step of crack depth, described crack depth is the original crack degree of depth and above-mentioned crack propagation increment sum;
7), step that above-mentioned crack depth and the tolerance limit size that is stored in computer-internal are compared, if crack depth is less than the tolerance limit size, output signal then, pressurized container keeps using; If crack depth is greater than the tolerance limit size, output signal then, closing pressure container.
The crack propagation increment is to be control signal with elastoplasticity cyclic J integration Δ J, gets by the computer conversion, thus the use of pilot pressure container.Crackle area increment can be reduced to
The present invention controls in real time pressurized container by computer and (by computer pressurized container is detected incessantly, immediately make and close or keep operational orders such as use, come the pilot pressure container by instruction again), for the generation of the sudden malicious event of prevention inservice pressure vessel, predict its fatigue surplus life, guarantee that the safe handling of pressurized container has great application value.The present invention with D elastic-plastic cyclic J integration Δ J as fracture parameter, adopt the Paris extended model, utilize finite element technique simulation tracing crackle tired expansion overall process and predict its fatigue surplus life.This is a kind of numerical simulation new technology of the pressurized container fatigue life prediction based on cyclic J integration Δ J, this technology has explicit physical meaning, intuitive is strong, calculation accuracy is high, with characteristics such as actual degree of agreement is good.
Description of drawings
Fig. 1 is the local grid floor map of the three-dimensional isoparametric element of 20 nodes;
Fig. 2 is the crack propagation schematic representation.
Embodiment
At first structure, material and the force-bearing situation according to pressurized container obtains the tolerance limit size that this pressure is easy to, and deposits this tolerance limit size in Computer Database, and be stand-by.
Qualitative by nondestructive detecting method (for example ray or ultrasound examination), quantitatively obtain defect type, position and the size of pressurized container, with this type of initial imperfection information input computer; To also import computer about the three-dimensional address information of node.
From the Paris-Erdogan extended model as can be known, can be represented by the formula for the Stable Crack Growth stage:
For the tired spreading rate of the crackle of a large amount of large scale yielding that exists in the engineering reality and complete yield breakage problem, replace Δ K with Δ J, therefore (1) formula becomes:
Wherein the representation based on the three-dimensional Δ J of tensor representation method is as follows:
In the formula: A
nBe the empty increment of expanding the crackle area that causes of crackle; σ
IjBe stress tensor; u
iBe displacement vector; W is a strain energy density; δ is the Kronecker tensor; Δ x
iIt is the mapping function of transformation of coordinates between two independent coordinate systems; V is the volumetric region that influenced by the empty expansion of crackle
The matrix form that is suitable for FEM (finite element) calculation is:
In the formula: N is the unit that influenced by the node virtual displacement; { J} is a Jacobian matrix; α is corresponding Weighting factor;
All the other cotypes (3).
Crackle area increment A
nCalculating very crucial, directly have influence on J integral Calculation precision, its concrete calculating is described below:
For 20 node isoparametric elements, because shape function has secondary, nodal displacement can cause distortion of the mesh phenomenon as shown in Figure 1.
According to finite element theory, any coordinate of crackle forward position can be expressed as among Fig. 1:
Corresponding shape function is:
In the formula: the s value at node 1,2,3 places is respectively-1,0,1.When node 3 moves to new coordinate
The time:
The coordinate of any point can be expressed as on the curve of being determined by node 1,2 and node 3 new coordinates:
So obtained the virtual displacement representation of any point on the new curve be:
Therefore obtain among Fig. 1 that the crackle area increment of dash area is between two curves:
In the formula (10): dl is the arc length increment;
Can certainly obtain the crackle area increment that is caused by node 1 or node 2 place's virtual displacement, its formula unification is expressed as:
Formula (11) is carried out Gauss integration, obtains:
In the formula (12): I is the Gauss integration point.
By above-mentioned derivation, obtained the formula of crackle area increment on 1 unit that causes by 1 node virtual displacement.It is pointed out that and when producing virtual displacement, as shown in Figure 1, can cause the variation of 2 unit grids in next door, so A by 1 cell node
nShould be the crackle area increment of 2 unit.
In the Practical Calculation process, need some disposal skill of utilization.When being positioned at node 1,3 on line centers as node 2, A
nCalculating just fairly simple.Because x
1+ x
3-2x
2=0, y
1+ y
3-2y
2=0, z
1+ z
3-2z
2=0, then Jacobi's representation J just can be reduced to:
In the formula (13): L represents wire length between the node 1,3.
Crackle area increment A correspondingly
nFormula can be reduced to:
The path of J integration should avoid the 1st layer of unit, crackle forward position to cause J integration situation bigger than normal because of splitting sharp Stress singularity as far as possible in addition, also should consider to reduce amount of calculation simultaneously, and chooses along the route that splits the 2nd layer of unit of point.
According to said method, utilize finite element technique just can calculate the cyclic J integration Δ J value of crack propagation on computers to the crackle forward position each point in each stage.
At a certain specific crackle, as shown in Figure 2, if crackle forward position curve defines with several discrete points, the change procedure of then getting 1 i on the crackle is as research object, after crack propagation, the i point becomes the i+1 point on the new leading surface of crackle, at this moment, for research point i, can be with formula (2) both sides integrating conversion:
In the formula: N
i, N
I+1, a
i, a
I+1Number of times and crack length when arriving i, i+1 for crack propagation
In the actual engineering, we can suppose that crackle forward position curve is the arbitrary shape of precognition, and the D elastic-plastic cyclic J integration Δ J value in crackle forward position can be obtained by FEM (finite element) calculation, then replaces (15) formula integration with the Eular integral equation, can get:
As long as control Δ a within the specific limits, crackle is from a
iExpand to a
I+1The time, available numerical method obtains the N that separates in the specified accuracy scope
I+1, recycling formula (16) just can calculate the whole expansion process of this point.
Concerning whole crackle forward position, suppose that the J integration in the initial forward position of crackle is tried to achieve, then obviously can use the Paris-Erdogan model at 2 for j, m:
Δa|
j=c(ΔJ)
nΔN|
j (17)
Δa|
m=c(ΔJ)
nΔN|
m (18)
In formula (17) and the formula (18), Δ N|
j=Δ N|
m, two formulas are divided by and can be got:
According to formula (19) as can be known, as long as obtain on the cyclic J integration Δ J value of each point on the crackle forward position and the crackle certain any expansion increment, just can obtain the crack propagation increment of crackle forward position each point.Can at first calculate the crackle deepest point in the specific implementation process and calculate other each point as a reference, and to define this crackle increment be Δ a
Max, just can calculate thus and the expansion situation of change of simulating crack pattern in whole fatigue process, and can obtain the fatigue surplus life of crackle thus.
Last circulation gained crack depth added the crack propagation increment of this circulation gained be the crack depth (the prediction crack depth value when this crack depth is worked for next) that this circulation can get.When if crack depth less than the tolerance limit size in the database, output signal then, pressurized container keeps using; If the crack propagation degree of depth is greater than the tolerance limit size, output signal then, closing pressure container.
The above is a specific embodiment of the present invention only, should be pointed out that for the person of ordinary skill of the art, can also make many modification and improvement, and all modification or improvement all should be considered as protection scope of the present invention.
Claims (3)
1, a kind of fatigue life safety predicting method for pressure container is characterized in that may further comprise the steps:
1), is used to import the input step of initial imperfection information
2), be used to import the input step of the three-dimensional address information of node;
3), be used to obtain the step of crackle area increment, adopt
4), be used to obtain the step of elastoplasticity cyclic J integration Δ J, adopt
5), be used to obtain the step of crack propagation increment, adopt
6), be used to obtain the step of crack depth, described crack depth is the original crack degree of depth and above-mentioned crack propagation increment sum;
7), step that above-mentioned crack depth and the tolerance limit size that is stored in computer-internal are compared, if crack depth is less than the tolerance limit size, output signal then, pressurized container keeps using; If crack depth is greater than the tolerance limit size, output signal then, closing pressure container.
2, fatigue life safety predicting method for pressure container according to claim 1 is characterized in that: described crack propagation increment is to be control signal with elastoplasticity cyclic J integration Δ J, gets by the computer conversion.
3, fatigue life safety predicting method for pressure container according to claim 2 is characterized in that: described crackle area increment can be reduced to
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100677460A CN1329689C (en) | 2004-10-31 | 2004-10-31 | Fatigue life safety predicting method for pressure container |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100677460A CN1329689C (en) | 2004-10-31 | 2004-10-31 | Fatigue life safety predicting method for pressure container |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1614294A CN1614294A (en) | 2005-05-11 |
CN1329689C true CN1329689C (en) | 2007-08-01 |
Family
ID=34765113
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100677460A Expired - Fee Related CN1329689C (en) | 2004-10-31 | 2004-10-31 | Fatigue life safety predicting method for pressure container |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1329689C (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644646B (en) * | 2009-07-07 | 2010-12-08 | 西安交通大学 | Fracture toughness measurement method based on optics |
CN103020426B (en) * | 2012-11-23 | 2015-11-18 | 北京航空航天大学 | A kind of short-cut method of rectangular slab Plate with Inclined Center Crack propagation life of fatigue prediction |
KR101517552B1 (en) * | 2013-08-22 | 2015-05-04 | 한국생산기술연구원 | Pressure vessel history management apparatus and its using method |
CN103743636B (en) * | 2014-01-16 | 2015-11-18 | 清华大学 | A kind of method predicting welding joint threshold in fatigue crack propagation |
US10288223B2 (en) * | 2015-11-20 | 2019-05-14 | Hexagon Technology As | Failure indicator supplemental vessel for primary vessel |
CN109596709B (en) * | 2018-12-19 | 2021-03-26 | 张磊 | Detection method of fixed pressure container |
CN117056686B (en) * | 2023-08-14 | 2024-02-02 | 嘉兴市安得特种设备科技有限公司 | Alarming method and system for detecting surface defects of pressure container |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4764882A (en) * | 1983-04-19 | 1988-08-16 | Kraftwerk Union Aktiengesellschaft | Method of monitoring fatigue of structural component parts, for example, in nuclear power plants |
US6491182B1 (en) * | 1994-10-13 | 2002-12-10 | Luxfer Group Limited | Treating pressure vessels |
US6594619B1 (en) * | 1999-08-02 | 2003-07-15 | Hood Technology Corporation | Apparatus and method for predicting failures of spinning disks in turbo-machinery |
-
2004
- 2004-10-31 CN CNB2004100677460A patent/CN1329689C/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4764882A (en) * | 1983-04-19 | 1988-08-16 | Kraftwerk Union Aktiengesellschaft | Method of monitoring fatigue of structural component parts, for example, in nuclear power plants |
US6491182B1 (en) * | 1994-10-13 | 2002-12-10 | Luxfer Group Limited | Treating pressure vessels |
US6594619B1 (en) * | 1999-08-02 | 2003-07-15 | Hood Technology Corporation | Apparatus and method for predicting failures of spinning disks in turbo-machinery |
Non-Patent Citations (3)
Title |
---|
压力容器及管道剩余寿命的评估方法 莫剑,化工装备技术,第25卷第5期 2004 * |
压力容器及管道剩余寿命的评估方法 莫剑,化工装备技术,第25卷第5期 2004;在役压力容器含缺陷结构腐蚀疲劳剩余寿命预测的随机分析 周昌玉,张艳丽,李强,黄文龙,压力容器,第19卷第2期 2002 * |
在役压力容器含缺陷结构腐蚀疲劳剩余寿命预测的随机分析 周昌玉,张艳丽,李强,黄文龙,压力容器,第19卷第2期 2002 * |
Also Published As
Publication number | Publication date |
---|---|
CN1614294A (en) | 2005-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhuang et al. | Extended finite element method: Tsinghua University Press computational mechanics series | |
Beer | Advanced numerical simulation methods | |
Cisilino et al. | Three-dimensional boundary element analysis of fatigue crack growth in linear and non-linear fracture problems | |
CN106952000A (en) | A kind of Karst Regional landslide disaster risk dynamic assessment method | |
Ichimura et al. | Three-dimensional nonlinear seismic ground response analysis of local site effects for estimating seismic behavior of buried pipelines | |
Wu et al. | A new method for classifying rock mass quality based on MCS and TOPSIS | |
CN1329689C (en) | Fatigue life safety predicting method for pressure container | |
Saouma et al. | Aging, Shaking, and Cracking of Infrastructures: From Mechanics to Concrete Dams and Nuclear Structures | |
Lu et al. | Analysis of FPSO dropped objects combining Monte Carlo simulation and neural network-genetic approach | |
Liu et al. | Cracking risk analysis of face slabs in concrete face rockfill dams during the operation period | |
Xie et al. | A life prediction method of mechanical structures based on the phase field method and neural network | |
Zhou et al. | Prospective forecast of sliding instability time using a precursory AE time series | |
Ahmed et al. | Homogenization of Honeycomb Core in Sandwich Structures: A Review | |
JPH08179047A (en) | Judging method of liquefaction caused by earthquake with neural network | |
Bouchet et al. | Prediction of the effects of neutron irradiation on the Charpy ductile to brittle transition curve of an A508 pressure vessel steel | |
QIU et al. | Dynamic Risk Assessment of the Subsea Tunnel Construction Process: Analytical Model. | |
Wei et al. | A novel method for identifying damage in transverse joints of arch dams from seismic responses based on the feature of local dynamic continuity interruption | |
CN103473425B (en) | Based on Discrete-time Model with Two Neurons mining induced stress effect work surface coal dilatation method of discrimination | |
Brüch et al. | Coupled poro‐mechanical tangent formulation applied to sedimentary basin modeling | |
Ichimura et al. | Comprehensive seismic response analysis for estimating the seismic behavior of buried pipelines enhanced by three-dimensional dynamic finite element analysis of ground motion and soil amplification | |
Zhu et al. | An efficient adaptive multi-mesh phase-field method for simulating rock fractures | |
Zheng et al. | Analysis of the mechanical characteristics and plasticity-hazardous zones development of natural gas pipeline under overload | |
Dascalu et al. | On a 3D micromechanical damage model | |
CN112861403B (en) | Safety evaluation method for large liquid storage tank structure | |
Korobkin et al. | Rational assessment of fluid impact loads |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |