The Life Predicting Calculations in Whole Process Realized by Calculable Materials Constants from short Crack to Long Crack Growth Process
Yangui Yu^{1, 2}
^{1}Zhejiang Guangxin New Technology Application Academy of Electromechanical and Chemical Engineering, Hangzhou, China
^{2}Wenzhou University, Wenzhou, China
Email address:
To cite this article:
Yangui Yu. The Life Predicting Calculations in Whole Process Realized by Calculable Materials Constants from short Crack to Long Crack Growth Process. International Journal of Materials Science and Applications. Vol. 4, No. 2, 2015, pp. 8395. doi: 10.11648/j.ijmsa.20150402.13
Abstract: To use the theoretical approach, to adopt the simple stress, or strainparameter method, by means of the conventional material constants, to establish some new calculation models in whole crack growth process for elasticplastic steels, which are the equations of the driving forces and the life predictions from micro to macro crack; To provide yet several crack growthratelinkingequations and life calculating expressions in whole process, for which are under different loading conditions in high cycle and low cycle fatigue. For key parameters A_{1}, A_{2} and B_{2} have proposed some new concept, and to define their physical and geometrical meanings. For the transition crack size from crack forming stage to crack growth stage, provide concrete calculation processes and methods. Thereby to realized the lifetime predicting calculations in whole process based on traditional calculable material parameters. The purpose is to try to make the modern fatiguefracture discipline depended on tests become gradually calculable subjects as the traditional material mechanics. So that will be having practical significances for saving testing manpower and funds, for promoting applying and development relevant disciplines.
Keywords: ElasticPlastic Materials, Fatigue Fracture, Crack Propagation Modeling, Low Cycle Fatigue, High Cycle Fatigue, Lifetime Prediction
1. Introduction
As everyone knows for the traditional material mechanics, that is a calculable subject, and has made valuable contributions for every industrial engineering designs and calculations. But it cannot accurately calculate the life problems for some structures when it is preexisting flaws and under repeated loading. In that it has no to contain such calculable parameters as the crack variable or as the damage variablein their calculating models. On the other hand, inside the fracture mechanics and the damage mechanics, due to there are these variables, so they can just calculate above problems. But nowadays latter these disciplines are all subjects mainly depended on fatigue, damage and fracture tests.
Author thinks, in the mechanics and the engineering fields, in which are also to exist such scientific principles of similar to genetic and clone technology in life science. Author has done some of works used the theoretical approach as above the similar principles [18]. For example, for some strength calculation models from micro to macro are provided by reference [1], for some damage growth rate calculation models from micro to macro damage growth are proposed by references [28] which are models in each stage even in whole process, under different loading conditions. Two years ago, in order to do the lifetime calculations in whole process on fatiguedamagefracture for an engineering structure, author was by means of Google Scholar to search the lifetime prediction models, as had been no found for this kind of calculation equations. After then, author continues to research this item, and bases on was provided and now is complemented on the comprehensive figure 1 of material behaviors [3]; still applies above genetic principles, to study and analyze data in references, thereby to provide some new calculable models for the new crack growth driving force and for the lifetime predictions. Try to make the fracture mechanics, step by step become calculable disciplines as the material mechanics. That way, may be having practical significances for decrease experiments, stint man powers and funds, for promoting engineering applying and developing to relevant disciplines.
2. The Life Prediction Calculations for ElasticPlastic Materials Containing PreFlaws
2.1. The Life Prediction Calculations in Short Crack Growth Process (Called the First Stage)
About life curves of short crack growth as the first stage are just described with curves 1 ( ), 2 and 3 ) in reversed direction coordinate system inside attach figure.1, where depicting for the coordinate system constituting and each curve meaning had explained in references [13].
(1) Under work tress (high cycle fatigue) condition
Under work tress condition, it is with as variable to adopt the simple stress parameter as its calculating one, and with the short crackstressfactor range to express the life prediction equation, that corresponding reversed curves 1 and 3 can be calculated by following equations
(11)
or
(12)
Where
(2)
(3)
Here the in eqn (2) is defined as short crack stress intensity factor, it is driving force of short crack growth under monotonic loading, and the crack stress intensity factor range in eqn (3) it is driving force under fatigue loading. The is defined as the comprehensive and calculable material constant. Author researches and thinks, its physical meaning of the is a concept of power, that just is a maximal increment value to give out energy in one cycle before to cause material failure. Its geometrical meaning is a maximal microtrapezium area approximating to beeline (Fig1), that is a projection of corresponding to curve 1 () or 3 on the yaxis, also is an intercept between. Its slope of microtrapezium bevel edge just is corresponding to the exponentof the formula (45).The comprehensive material constant in formulas (45) is a calculable one, it has function relation with other material constants and ,etc. theis a fatigue strength coefficient.
(4)
(5)
Here should yet explain the in eqn (4) is corresponding reversed curves 1, its mean stress ; The in eqn (5) is corresponding curves 3, its . And the correctional method for its mean stress can be corrected by reference [9].
Where
(6)
or
(7)
The in eqns (67) is defined as an effective damage history correction factor in first stage, its physical meaning is the effective crack growth rate of whole failure to cause specimen material in a cycle, its unit is. is a reduction of area. is premicrocrack value that is no effect to fatigue damage under prior cycle loading [10]. is an initial mirocrack value, is a critical fracture size before failure. is initial life in first stage,; is failure life,. By the way, here is also to adopt those material constants as "genes" in the fatigue damage subject.
So, for the eqn (1), its final expansion equation corresponded reversed to curves 1 is as below form:
(8)
And its final expansion equation corresponded reversed to curves 3should be:
(9)
Where is a transitional crack size from short crack to long crack growth process, ,is a long crack size corresponding forming macro crack.is a medial crack size between initial crack size and transitional crack size corresponding medial life .
(2) Under work stress (low cycle fatigue) condition.
Under condition, here it is still with as variable, but it should adopt the simple strain as its calculating parameter, and it should adopt the short crackstrainfactor range to express its life equation. That is corresponded to reversed direction curve in Fig1, is as following form
(10)
Here
(11)
(12)
Where theis also calculable comprehensive material constants.
(13)
If via the crack stress factor amplitude in eqn (10) to express it, due to plastic strain occur cyclic hysteresis loop effect, it should be
(14)
Where the is also a calculable comprehensive material constant:
(15)
(16)
Or
(17)
Where is a cyclic strength coefficient. is defined to be fatigue ductility exponent, , , just is a fatigue ductility exponent under low cycle fatigue, , is a strain hardening exponent. So that, its final expansion equation for (10) is as below form,
(18)
Its final expansion equation for (14) is as following form,
(19)
If to take formula (17) to replaceinto eqn. (14), its final expansion equation is as below forming
(20)
Here influence of mean stress in eqn (1920) can be ignored.
2.2. The Calculations for FatigueDamage in Long Crack Growth Process （or Called the Second Stage）
In Fig.1, the residual life curves of long crack growth in the second stage are just described with curves 1' ( ), 2' and 3' ) at reversed direction coordinate system.
(1) Under work tress condition
Here it is divided two methods: factor method and method:
1) factor method
For life prediction equation corresponded to reversed curves and should be as following
(21)
Where
(22)
(23)
As it is known the factor is just the macro crack stress intensity factor, and the is the macrocrack stress intensity factor range. Theis correction factor [11] related for long crack form and structure size. The in eqn.(21) is defined as comprehensive material constants of macrocrack, for, it is corresponding curve, and it is also calculable one as following
(24)
Or
(25)
For , the is corresponded to curve , it should be as following form
(26)
Where is mean crack stress intensity factor, is an effective crack stress intensity factor, is a critical crack stress intensity factor, which they are parameters under cyclic loading. It should be point that the physical meaning for theis also a concept of power, that just is a maximal increment value to give out energy in one cycle before failure. Its geometrical meaning is also a maximal microtrapezium area approximating to beeline (Fig1), that is a projection of corresponding to curve 1’ () or 3’ on the yaxis, also is an intercept between. Its slope of microtrapezium bevel edge just is corresponding to the exponentof the formula (24~26). Here,
(27)
Author research and think, the parameter is defined to be the virtual rate, it is an equivalent rate caused in precrack, it can take similar dimension with the "" () by reference [12]. But its unit is different, here unit of theis "". The crack size is a virtual crack size as equivalent to a precrack size, is an initial as equivalent to the initial microcrack size. is an initial life, . is a virtual life, . In references [1314], all refer to effective stress intensity factor, here to propose to take equivalent values as follow
(28)
(29)
Here the is a threshold crack stress intensity factor value.
or
(30)
(31)
(32)
(33)
It should be point that is a true fracture stress under monotonous load, and the is the fatigue strength coefficient under fatigue load.
So the effective life expanded equation corresponding reversed direction curve is following forming.
(34)
And the effective life expanded equation corresponding reversed direction curve should be
(35)
Its medial life in second stage is
(36)
Where is a transitional crack size between two stages from short crack to long crack growth process, , the is a medial size. .
2) method
Due to word stress is still ( ), in the long crack growth process, its residual life equation of corresponding reversed direction curve and in fig.1 is as following form
(37)
Where
(38)
(39)
Theis a crack tip open displacement, it is driving force of short crack growth under monotonic loading; and the is a crack tip open displacement range, it is driving force under fatigue loading. For the coefficient in eqn (39), it equal 1() under plane stress condition, under plane strain condition, . is Poisson's ratio. is an elasticity modulus. The is a comprehensive material constant in second stage, its physical meaning of theis also a concept of power. Its geometrical meaning is also a maximal microtrapezium area approximating to beeline (Fig1).is also calculable comprehensive material constant, for that is
(40)
And for that is
(41)
Where and are maximum and minimum work stress. The can be calculable effective crack size.
So its final expansion form corresponded reversed direction curve for eqn (37) is as below,
(42)
And the life equation corresponded to reversed direction curve is following
(43)
(2) Under work tress condition
Under condition, its effective life models corresponded to reversed curve in figure 1 is as below form,
(44)
Where is also a calculable comprehensive material constant,
(45)
(46)
The is defined to be ductility exponent in long crack growth process, , is a fatigue ductility exponent in second stage.
So that, the final expansion equations is derived from above mentioned eqn. (44) as follow
For ,
(47)
For , it should be
(48)
In the eqn (48), influence to mean stress usually can ignore.
Where, theis an effective crack size, it can calculate from effective crack tip opening displacement
(49)
And
(50)
Where the is critical crack tip displacement. So the in (4748) can be calculated by value by means of equations (4950).
2.3. The Life Prediction Calculations for Crack Propagation in Whole Process
Due to short crack behavior and long crack one are distinctly different, for availing to life calculation in whole process, it should take a crack size at transition point from short crack to long crack growth process, and the transition point sizecan be derived to make equal between the crack growth rate equations by two stages[3,15], for instance
(51)
Here the equation (51) is defined as the crack growth ratelinkingequation in whole process.
(1) Under work stress
For, its expanded the crack ratelinking equation for eqn (51) corresponded to positive curve is as following form
(52)
For, its expanded the crack ratelinkingequation for eqn (51) corresponded to positive curve is as following form
(53)
And the life equations in whole process corresponding to reversed direction curves and should be as below
(54)
Its expanded equation corresponding to reversed direction curves is as following form
(55)
But for , its expanded equation corresponding to reversed direction curves should be
(56)
(2) Under work stress
Under work stress, its expanded rate link equation for eqn (51) corresponding to positive curve is as following form
(57)
The life equations in whole process corresponding to reversed direction curve should be as following
(58)
And the life prediction expanded expression in whole process corresponded reversed curve , it should be
(59)
It should point that the calculations for rate and life in whole process should be according to different stress level, to select appropriate calculable equation. Here it have to explain that its meaning of the crack ratelinkingequation (5153, 57) is to make a calculable linking formula between the first stage crack rate and the second stage one, which it should be calculated by the short crack growth rate equation before the transition point size ; it should be calculated by the long crack growth rate equation after the transition point, note that it is not added together by the crack growth rates for two stages. But the life calculations for two stages can be added together. About calculation method, it can calculate by means of computer doing computing by different crack size [16].
3. Calculating Example
3.1. Contents of Example Calculations
To suppose a pressure vessel is made with elasticplastic steel 16MnR, its strength limit of material, yield limit, fatigue limit , reduction of area is , modulus of elasticity; Cyclic strength coefficient , strainhardening exponent ; Fatigue strength coefficient, fatigue strength exponent , ; Fatigue ductility coefficient , fatigue ductility exponent , . Threshold value, critical stress intensity factor . Working stress , in pressure vessel. And suppose that for long crack shape has been simplified via treatment become an equivalent throughcrack, the correction coefficientof crack shapes and sizes equal 1, i.e.. Other computing data are all included in table 1.









92.7  28.23  8.6 
 3.91  0.18  2.9  1.0  0.07 
3.2. Required Calculating Data
Try to calculate respectively as following different data and depicting their curves:
(1) To calculate the transitional point crack sizebetween two stages;
(2) To calculate the crack growth rate at transitional point ( at crack size )
(3) To calculate the life in first stage from micro crack growth to transitional point
(4) To calculate the life in second stage from transitional point to long crack size ;
(5) Calculating the whole service lifetime.
(6) Depicting their life curves in whole process.
3.3. Calculating Processes and Methods
The concrete calculation methods and processes are as follows
(1) Calculations for relevant parameters
1) Stress range and mean stress calculations:
Stress range:
Mean stress:
2) According to formulas (7), calculation for correction coefficient in first stage.
For the effective crack sizes in first stage and the second stage, it both can be calculated respectively, and can take smaller one of both, here to take same value with the second stage, . For example, according to formulas (49), Calculating effective size
Take
3）By eqn (27), to select virtual rate in second stage, here take:
(2) To calculate the crack sizeof transitional point between two stages
1) To select calculating equation of short crack growth rate
At first, calculation for comprehensive material constantby eqn (17)
Here select the crack growth rate linking equation (57), and for its growth rate to simplify as follow form,
2) By rate linking equation (57), calculating for crack growth rate in long crack growth process,
Calculation for comprehensive material constantby eqn (46)
For the long crack rate in second stage, to take brief calculations as follow form,
3) Calculation for crack sizeat transitional point
According to the equations (51) and (57), for calculation the crack size at transitional point, it can make equal between both rate expansion equations, and to simplify it as following
Result, to obtain crack sizeat the transitional point.
(3) To calculate the crack growth rate at transitional point
Here it can be seen, the growth rate at the transition point is same.
(4) Life prediction calculations in whole process
1) To select eqn (20), to calculate the life from microcrack to transitional point is as follow,
So predicting life in first stage
And for above formulas, we can derive simplified life equation corresponded to different crack size as follow form
2) To select eqn (48), to calculate the life from transitional point to is as follow,
From above formulas, we can also derive simplified life equation corresponding different crack size as follow form
Therefore, predicting lifetime in whole process is
The life data corresponded to different crack length is all included in table 2, 3 and 4. Author finds these data are calculated by the simple stress or strainparameter method are basically close as compared with another result data that are calculated by the twoparametermultiplicationmethod in damage mechanics, where it has been published recently in reference [3]. Especially, those dada in second stage are closer, and the calculating models and method is simpler than the twoparametermethod.
Data point of number  1  2  3  4  5 
Crack size (mm)  0.02  0.04  0.1  0.2  0.4 
Data of the first stage  51020408  25510204  10204082  5102041  2551020 
Data of the second stage  Invalid section 
Data point of number  5  6  7 transition point  8 
Crack size (mm)  0.5  0.6  0.789  1.0 
Data of the first stage  2040816  1700680  1293293  1020408 
Data of the second stage  4851966  2859513  1292431  650026 
Data point of number  9  10  11  12  13 
Crack size (mm)  1.5  2.0  3.0  4  5 
Data of the first stage  680272  510204  Invalid section  
Data of the second stage  200570  87085  26871  11667  6108 
(5) To depict the life curves in the whole process
By means of the data in tables 24 mentioned above have depicted the life curves for two stages and whole course are respectively in figure 2 and 3.
(A) 11 data curve in first stage obtained by singleparameter calculating method;
(B) 12data curve in second stage obtained by singleparameter calculating method;
(C) This example transition point from microcrack size 0.02mm to long crack size 5 is just at seventh point (crack size ).
(A) 11 data curve in first stage obtained by singleparameter calculating method;
(B) 12data curve in second stage obtained by singleparameter calculating method;
(C) This example transition point from microcrack size 0.02mm to long crack size 5 is just at seventh point (crack size ).
4. Discussions and Conclusions
(1) About comparison of calculating results data: The lifetime data calculated by the simple stress or strainparametermethod as compared with the twoparametermultiplicationmethod in damage mechanics [3], which is basically close, the calculating models and method is also simpler than it.
(2) About new understanding for some material constants: True material constants must show the inherent characters of materials, such as the and etc in the material mechanics; for instance theand ; and and so on in the fatigue damage mechanics; and so on in the fracture mechanics, above which could all be checked and obtained from general handbooks. But some key materials constants ,and in the fracture mechanics which are essentially functional formulas with other parameters the and etc. to have functional relations, for which can all be calculated by means of the relational expressions e.g. eqns (45), (13,1516), (2426), (4041), (4546), etc. that has to combine experiments and to verify. Therefore for this kind of material constants can be defined as the calculable comprehensive materials constants.
(3) About cognitions to the physical and geometrical meanings for key parameters and : The parameters in the first stage and in the second stage, they are all a concept of power, just are a maximal increment value paying energy in one cycle before to cause failure. Their geometrical meanings are a maximal microtrapezium area approximating to beeline.
(4) About the methods for crack propagation rate and lifetime calculations in whole process: Calculation for crack transition size it can be calculated form the crack growthratelinkingequation (5153) and (57); calculation for the short crack rate it should be calculated by the short crack rateequations before at transition point ; for the long crack rate it should be calculated via long crackrateequations after transition point . But for the lifetime calculations in whole process can be added together by life cycle number of two stages.
(5) Based on the traditional material mechanics is a calculable subject; in consideration of the conventional materials constants there are "the hereditary characters"; In view of the relatedness and the transferability between related parameters among each disciplines; And based on above viewpoints and cognitions of the (1)~(4); then for the fatigue and the fracture disciplines, if make them become calculable subjects, that will be to exist possibility.
Acknowledgments
At first, author sincerely thanks scientists David Broek, Miner, P. C. Paris, Coffin, Manson, Basquin, Y. Murakami，S. Ya. Yaliema，Morrow J D, etc, who they have be included or no included in this paper reference, for they have all made out valuable contributions for the fatiguedamagefracture subjects. Due to they hard research, make to discover the fatigue damage and crack behavioral law for materials, to form the modern fatigue, damage, and fracturemechanics; due to they work like a horse, make to develop the fatigue, damage, and fracture mechanics subjects, gain huge benefits for accident analysis, safety design and operation for which are mechanical equipments in engineering fields. Particularly should explain that author cannot have so many of discovery and provide above the calculable mathematical models and the figure 1, if have no their research results.
Author thanks sincerity the Zhejiang Guangxin New Technology Application Academy of Electromechanical and Chemical Engineering gives to support and provides research funds.
References