The Computing of Intersectant Relations for Its Strength Problem on Damage and Fracture toMaterials with Short and Long Crack

Adopt two types of damage variables, a andD, and make the bidirectional combined coordinate system and the bidirectional curves in the whole process, describing their damage evolutive behaviors on fatigue damage-fracture to elastic-plastic steels; communicate the cross-referencing between their computing models and between describing curves at each stage among varied disciplines; bring forward a viewpoint about the driving force of material damage, that is the damage stress factor at crack forming stage, Provide the computation expressions and computing methods of the strength problem of materials with short crack and long crack at each stage; reveal the geometrical and the physical meaning of force triangle and its edge vector at each stage, Provide the conversion methods between the variables, the equations, the material constants and the dimensional units; Indicate the physical and the geometrical meanings for some key parameters. This will be having practical significance for promoting developing, and applying each discipline.


Introduction
As is now well known, we adopt the crack size a as a variable in the fracture mechanics to describe crack growth process undergoing damage for a material and we adopt a damage variable D in the damage mechanics to describe an evolutive process undergoing damage for one.Either the sign a or the sign D are all virtually damage variables so, we could also adopt the damage D to describe the evolutive law of a structure material with crack.References [1][2][3] had made out such research.In each discipline we have all in-house features and its advantages.If we can communicate and convert corresponding relations for that among the damage variables, and the equations, the material constants, the dimensional units which describe the material behavior for varied discipline and provide some conversion methods, thus we are also able to adopt the same variables D 1 and D 2 or the variable a 1 and a 2 to compute the strength and the life at each stage or even in overall process for structures and materials undergoing fatigue damage [4,5].And the conventional materials and damage mechanics are made by inheritance and development and the modern one is made by all better combination and application.Based on this aim THE authors adopt the mathematical derivation and computational analyses used with computer by long-range research, educe a series of the computation expressions and the computing methods.Thus, this may be having practical significance for promoting, developing, and applying of some disciplines.

Bidirectional Combined Coordinate System and Bidirectional Curves in the Whole Process
In some of branch disciplines on fatigue-damage fracture, for finding their correlations among variables, curves, equations, and material constants of describing material behaviors at each stage and for connecting their relations to each other, we must put up analysis and developments for a number of problems that are above mentioned.Here, it is by means of bidirectional combined coordinate system that is, adopted in Figure 1 [6] we express the damage evolving process of material behavior at each stage and in the whole course, which consists of six

Relations among Force Triangles and Relations among Edge Vectors of Force Triangle Itself at Each Stage
Due to extent of material undergone damage at each stage is different and due to the variance in the material behaviors after undergoing damage, the rigidity of material is also different to ensue from change, so the force triangles consist of edge vectors at different stages which are also varied relationships among force triangles, and relationships among the edge vectors themselves at each stage, their geometrical and physical meanings are all compiled in Table 1.We can make out from the force triangle in each area from Figure 1 find the mathematic model of driving force is different in each discipline.The driving force seen from force triangle between the axes O I and OI is the damage stress intensity factor under ultra-high cycle fatigue.We should also point out that this is the calculation domain in micro-damage mechanics or in conventional material mechanics where the ΔG uh and ΔG uh are, respectively, the microcrock stress intensity factor or the microdamage stress intensity factor under ultra-high cycle fatigue to correspond a crystal grain or microdamage defect a μ (D μ ) which is originated from interior materials.And the mathematic models of driving force between the axes O 1I and O 3III are the stress intensity factor range ΔH 1 of shortcrack or the damage stress intensity factor range ΔH 1 . or that is calculation domain of the microfracture mechanics or damage mechanics.And the mathematic models of driving force between the axes O 3III and O 4IV are the stress intensity factor range ΔK 2 of longcrack or the damage stress intensity factor range ΔK 2 at the second stage or that is calculation domain of the macro-fracture mechanics.
Relations among force triangles can be expressly made out in Figure 1; the relations among edge vectors of force triangle iteself and their geometrical and physical meaning at each stage are all compiled in Table 1.

Computing the Damage Strength for Elastic-Plastic Material under High-Cycle Fatigue
Based upon above-mentioned viewpoints we bring forward the computation expressions and calculation methods of strength for elastic-plastic steels which are the material  [7,8].Murakami [9] provides the computation expression as follows: And research and analysis of the driving force of force triangle in Figure 1 and according to above explanation and the present author provide again a computational model to describe stress-strain about the tip of short crack that is the stress intensity factor at the first stage as following [10][11][12]: where H mac = Critical stress intensity factor of short crack that is corresponding to macro-crack size at about the threshold level ΔK th and more than one, that is artificially definite according to the size of structure member.y 1 is a correction coefficient concerned with shape and size, for example, of a structure member.And it became the damage stress intensity factor at the same stage when we adopt damage variable D 1 to describe [13] where H mac : Critical damage stress intensity factor that is equivalent to the H mac .D mac : Critical damage level that is equivalent to a mac ; if a mac = 0.7 mm, then D mac = 0.7.
Conversion method between variables a and D and between the dimensional units defines 1 mm equivalent to one damage unit (nondimensional value), 1 m equivalent to 1000 damage unit.If critical a mac = 0.7 mm, then critical damage value of the equivalent is D mac = 0.7 damage unit.And define that a 0 < a 1 ≤ a mac ; D 0 < D 1 ≤ D mac .The D 0 is an initial damage value corresponding to micro-crack size a 0 , and it is advised to take the average value of 10 crystal sizes or to take the maximum size of crystal size.It should also be pointed out that the stresses are all local ones in (2)-( 5) and ( 7)- (10).

Computing of Strength at the Second Stage.
The computational model to describe stress strain about the tip of long crack at the second stage had been provided by famous scientists Broek and Hellan as follows [14,15]: 11) where K Ic (= K 2c ) = Critical stress intensity factor of long crack y 2 is a correctional coefficient concerned with the shape of crack and shape and size of structure member.And it becomes the damage stress intensity factor at the same stage if we adopt damage variable D 2 to describe [16,17] K 2c is a critical damage stress intensity factor equivalent to the K Ic (K 2c ).D 2c is a critical damage value correspondent to the critical crack size a 2c .
Conversion method between variables a 2 and D 2 and between the dimensional units at the second stage is the same as the first stage.It should be pointed out that the dimensional unit of the damage stress intensity factor H and the K 2 at each stage become identical with the unit of stress because the dimensional unit between the variables a and D is converted.
It should be pointed out that both the values between the stress intensity factor of short-crack and the stress intensity of long crack are different under the condition of the same crack size (e.g., a 1 = a 2 = 0.1 mm) because the mechanisms of material damage and because their mathematical models and dimensional units of both driving forces are all different.But both growth rates at the turning point from short crack to long crack should be accordant or near.

Computing Example
A pressure vessel adopts the steel 16 MnR to make its strength limit of material σ b = 545 MPa, yield limit σ y = 349 MPa, the strain-hardening exponent n = 0.136, fatigue strength exponent in short crack growth m 1 = 11.478,threshold level ΔK th = 6.87 MPa Try respectively, to compute the stress intensity factor H and the damage stress intensity factor H at short crack size a 0 = 100 μm undergone after damage for a crystal grain and the stress intensity factor K 2 = K I and the damage stress intensity factor K 2 at long crack size a 2 = 2 mm.Computing approach is as follows.

Computing the Stress Intensity Factor H and the Damage
Stress Intensity Factor H for Short Crack

Computing the Stress Intensity Factor H for Short Crack.
According to (7), select the computing parameter: a 1 = 100 μm = 1.0 −4 m; the local stress is 840 MPa at focal point of stress undergone damage; If we take y 1 = 1.1, This results in So this short crack is growth.

Computing of the Damage Stress Intensity
Factor H for Short Crack.According to (9), when the a 1 = 100 μm = 1.0 −4 m, equivalent damage value of short crack is the D 1 = 100 = 1.0 −4 (damage unit).So ( This results in So this short crack is also growth.(11), select a parameter, take the correction coefficient y 2 = 1.05, a 2 = 2 mm, working stress σ = 280 MPa;

Computing the
This results in So the pressure vessel is safe.

Computing the Damage Stress Intensity Factor K 2 for
Long Crack.According to equation ( 13), take the computing parameter y 2 = 1.05; when take a 2 = 2 mm, the equivalent damage value of long crack is This results in So the pressure vessel is also safe.

Summarization
The bidirectional combined coordinate system, the bidirectional curves in the whole process and their force triangles at each stage are important scientific method and tool to communicate the cross-referencing which is to describe the evolutive process of a material behavior undergoing of the fatigue damage in each of the disciplines, which are able to make available communications and conversions for those complicated correlations among some variables, some equations, curves, and dimensional units, which are clearely able to explain the geometrical and physical meanings for the key parameters.Thus, it is also able to adopt same variables D 1 and D 2 or the variables a 1 and a 2 to compute the strength and the lift at each stage or even in overall process for structures and materials undergoing a fatigue damage.
And the conventional material and damage mechanics are able of making inheritance and development, and the modern one are all able to make better combination and application.Thus, that may be having practical significance for promoting developing, and applying some disciplines.

Conclusions
(1) About the Problem of the Relation between the Mathematical Model and the Material Behavior.Under identical loading, when the structure material is undergoing fatigue damage, the differences between the mathematical models to describe the material behavior are due to the degree of damage undergone at varied stages that makes the stiffness of the material change, which find that expressions in the curves cbaBA 1 A 2 in evolutive process are turned to take place at the points b, B, and A 1 .The slopes of the curves at each stage also are brought to change.The exponents in the equation also became from the m 1 = −1/b 1 at crack forming stage to the m 2 = −1/b 2 at crack growth stage.The driving forces became from the H 1 (H 1 ) to the K 2 (K 2 ).The compositive material constants became from the A 1 to the A 2 .
(2) About the Problem of Driving Force.It well known that the driving force of long crack growth in macro fracture mechanics is the stress intensity factor K I but the driving force of microdamage at crack forming stage is defined by the author as the damage stress intensity factor H 1 , and the driving force in microfracture mechanics is defined as the stress intensity factor H 1 of short-crack growth.
(3) About Problem of Dimensional Units.The stress intensity factor H 1 and the damage stress intensity factor H 1 of short crack at crack forming stage are all to describe the stress strain about the crack tip of micro-crack, the unit of the H 1 is MPa • m 1 √ m, and the unit of the H 1 is the MPa m 1 D 1 .Both the units are different, but the unit of the damage stress intensity factor H 1 is the same as the unit of the stress, and the H 1 and the H 1 are all essentially the stress intensity factor and a relation of equivalents.And the stress intensity factor K 2 and the damage stress intensity factor K 2 of long crack at crack growth stage are all to describe the stress strain about the crack tip of long crack, the unit of the K 2 is MPa • √ m, and the unit of K 2 is the MPa • D 2 .Both the units are also different, but the unit of the damage stress intensity factor K 2 is also the same as the unit of the stress, and the K 2 and the K 2 are all essentially the stress intensity factor and a relation between equivalents.In addition, both values and dimensional units between the stress intensity factor of short crack and the stress intensity factor of long crack are also different under the condition of same crack size.

Nomenclature
(1) D = damage variable in the whole process, D μ , D 1 = microdamage variable equivalent to micro-crack a μ or short crack a 1 at the crack forming stage (first stage), and D 2 = damage variable at the crack growth stage (the second stage).
(2) a μ , a 1 = micro-, short-crack size at the crack forming stage (variable of the first stage), a 2 = macro-, longcrack size at the crack growth stage (variable of the second stage), a 10 = initial size of micro-crack forming (ordaining value), crack size a th = corresponding to the threshold level ΔK th , original size a 20 = a mac of macro-crack forming stage (ordaining value), and a 2c = critical size of long crack.
(4) ΔG or ΔG = micro-crack stress intensity factor range or microdamage stress intensity factor range corresponding to microcrack size a μ or microdamage D μ under ultra-high cycle fatigue, H 1 = stress intensity factor of short crack, H mac = critical stress intensity factor of short crack, H 1 = damage stress intensity factor of short crack and ΔH 1 or ΔH 1 or ΔH 1 /2 = stress intensity factor range or stress intensity factor amplitude relative to short-crack a 1 ; ΔH 1 /2 = damage stress intensity factor range or damage stress intensity factor amplitude relative to damage variable D 1 .
ΔI or ΔI/2 = damage strain factor range or damage strain factor amplitude relative to short crack a 1 .(6) dD/dN = damage evolutive rate, dD 1 /dN 1 = damage evolutive rate at the crack forming stage, dD 2 /dN 2 = damage evolutive rate at the macro-crack growth stage, da/dN = crack growth rate, da 1 /dN 1 = short crack growth rate at the crack forming stage, and da 2 /dN 2 = its rate at the macro-crack growth stage.
(7) N oi = life of correspondance to medial damage variable D oi or short-crack medial size a oi at the first stage, and N oj = life of correspondance to medial damage variable D oj or long-crack medial size a oj at the second stage.
(8) K 2 = K I = stress intensity factor of long crack, K 2 = damage stress intensity factor of long crack, J-integral of long crack; and crack tip opening displacement of long crack, K m = mean stress intensity factor, ΔK/2 = stress intensity factor amplitude of correspondance to macro-crack a 2 ,ΔJ/2 = J-integral amplitude corresponding to macro-crack a 2 , Δδ t /2 = crack tip opening displacement amplitude corresponding to macro-crack a 2 .
(9) K 1c = critical stress intensity factor corresponding to macro-crack critical a 2c , K eff = effective stress intensity factor to be applicable in Paris's equation, J c = critical J-integral value corresponding to macrocrack critical a 2c , and δ c = critical crack tip opening displacement corresponding to macro-crack critical a 2c .

√m, 1 √m
critical stress intensity factorK 1c = 97.3MPa √ m = K 2c = K Ic = 97.3(MPa √ m), critical damage stress intensity factor K 2c = 97.3MPa D 2 equivalent to the K Ic (K 2c), the critical stress intensity factor H mac = 350 m and the damage stress intensity factor H mac = 350 (MPa m 1 D 1 ) of short crack corresponding to the threshold level ΔK th , and the mean sizes a 0 = 20 ∼ 30 μm of crystal grains; its working stress is 280 MPa of pressure vessel and local stress is 840 MPa at focal point of stress.

( 5 ) a 1 =
fatigue strength exponent under ultra-high cycle fatigue, b 1 = fatigue strength exponent under high cycle fatigue, c 1 = fatigue ductility exponent under low cycle fatigue, n 1 = fatigue strength exponent in micro-crack growth rate equation under ultra-high cycle fatigue, n 1 = −1/a 1 ; m 1 = fatigue strength exponent in short-crack growth rate equation under high-cycle fatigue, m 1 = −1/b 1 , m 1 = fatigue ductility exponent in short-crack growth rate equation under low cycle fatigue, m 1 = −1/c 1 , b 2 = fatigue strength exponent of the macro-crack growth stage under high-cycle fatigue, c 2 = the fatigue ductility exponent at the macro-crack growth stage under low-cycle fatigue, m 2 = the fatigue strength exponent in crack growth rate equation under high-cycle fatigue, m 2 = −1/b 2 , m 2 = the fatigue ductility exponent in crack growth rate equation under low-cycle fatigue, m 2 = −1/c 2 .
abscissa axes O I , O I , O 1I , O 2II , O 3III , and O 4IV and two bidirectional ordinate axis O 1 O 4 and O 1 O 4 .Between the axes O I and O 1I , it is the calculation domain of the conventional material mechanics; between the axes O I and O I , it is the calculation domain of the current ultrahigh cycle fatigue; among the axes O I , O 1I , and O 2II , it is the calculation domain of the damage mechanics and the microfracture mechanics; between the axes O 3III and O 4IV , it is calculation domain of the macrofracture mechanics; Between the axes O 2II and O 3III , it is all applicable calculation domain for the microfracture mechanics and macrofracture mechanics.The upward direction along the ordinate axis is presented as damage evolving rate dD/dN or crack growth rate da/dN (that can also carve up the damage evolving rate dD 1 /dN 1 or short crack growth rate da 1 /dN 1 at crack forming stage or the damage evolving rate dD 2 /dN 2 or long crack growth rate da Consequently, the distance O O 3 is shown as the region relating to life N mac from grains size to microcrack initiation until macrocrack forming; the distance O O 4 is shown as the region relating to the lifelong life 2N from micro-crack initiation until fracture of structure material.The coordinate system combined from upward axis O O 4 and abscissa axes O I , O 1I , and O 2II is presented to be the relationship between the damage evolving rate dD 1 /dN 1 (or the short crack growth rate da 1 /dN 1 ) and the damage stress factor amplitude ΔH/2 (or damage strain factor amplitude ΔI/2) at uhf .The ABA 1 shows the varying regularities of elastic material behaviors as under high cycle loading at macro-crack forming stage; positive direction ABA 1 shows the relation between dD 1 /dN 1 (or da 1 /dN 1 ) and ΔH/2; inverted A 1 BA shows the relation between the ΔH/2 and 2N.The curve CBC 1 shows the varying regularities of plastic material behaviors, as is under low-cycle loading at macrocrack forming stage; positive direction CBC 1 shows the relation between da 1 /dN 1 and ΔI/2; inverted C 1 BC shows the relation between the Δε p /2-2N.And the curve A 1 A 2 at crack growth stage is showed as under high cycle loading: positive direction A 1 A 2 ; shows da 2 /dN 2 -ΔK/2 (ΔJ/2); inverted A 2 A 1 , shows the relation between the ΔK/2 and ΔJ/2-2N.The C 1 C 2 shows the positive direction relation between the da 2 /dN 2 and Δδ t /2 under low-cycle loading, inverted C 2 C 2 /dN 2 at crack growth stage), and the downward direction is presented as each stage life 2N.The distance O O between axis O I and O I is shown as the region of the nominal stress S or remote stress σ o ; the distance O O 2 between axis O I and O 2II is shown as the region from uncrack to microcrack initiation; the distance O 2 O 3 between axes O 2II and O 3III is shown as the region relative to life N oi mic-mac from microcrack growth to macrocrack forming.
da 1 /dN IFigure1: Bidirectional combined coordinate system and bidirectional curves in the whole process.

Table 1 :
Relations among the edge vectors of the force triangle iteself and their geometrical and physical meanings at each stage.
4.1.Computing of Strength for Material with Crack 4.1.1.Computing of Strength at the First Stage.About computing the material strength with short crack, many researchers had presented varied computing models and made out valuable contributions Stress Intensity Factor K 2 = K I and the Damage Stress Intensity Factor K 2 for Long Crack 5.2.1.Computing the Stress Intensity Factor K 2 for Long Crack.According to