# Analytical solutions for crack initiation angle of mixed mode crack in solid material/Kietojo kuno misraus pavidalo plysio pradzios kampo analitinis nustatymas.

1. IntroductionThere are a lot of cracks in solids; a crack will propagate under biaxial loading conditions, forming wing cracks [1]. As the crack slips under shear stress, the slip itself causes wing cracks to grow at the crack tip. The study of mixed mode crack initiation and propagation is very important for understanding the mechanical characteristics of solids [2, 3]. The initiation and propagation problem of mixed mode cracks is of great interest in the field of fracture mechanics. To accurately predict and simulate crack propagation paths, a valid criterion needs to be defined for the crack initiation angle. The crack propagation path is influenced by the value of the crack initiation angle [4]. Because the actual loading condition at the crack tip in solid material is very complex, the crack initiation angle and propagation under mixed-mode loading conditions cannot be simplified to either pure mode I or pure mode II analysis.

Erdogan and Sih [5] were the first to study the problem of mixed mode crack initiation angles. They presented the maximum tangential stress (MTS) criterion based on the tangential stress for crack initiation under mixed mode stress conditions. Since then, several criteria have been proposed for the mixed mode crack initiation problem, such as strain energy, stress triaxiality and strain fields [6-8]. These fracture criteria suppose that a plastic core zone exists at the crack tip, where in simplified form most of the plastic core region is a circle with a constant radius and some non-uniform regions are defined based on the mechanical characteristics of the materials. For the MTS criterion [5] and the minimum strain energy density criterion [9], a circle plastic region is assumed to exist at the crack tip. However, the maximum dilatational strain energy criterion (T-criterion) [7] and the modified MTS criterion [10] use a variable plastic region. In fact, the shape of plastic core region is influenced by the microstructure of the materials, because most solid material have nonlinear, discontinuous, anisotropic mechanical characteristics. The shape of the plastic core region computed by an empirical formula cannot reflect the real situation, and numerical methods need to be adopted to simulate the plastic region at the crack tip [11-15].

The main objective of this study is to present a method for mixed mode crack initiation based on a failure criterion. Key issues include: the stress state at the crack tip, computation of stress intensity factors based on crack tip stress and the crack initiation angle. The remainder of the paper is organized as follows. In Section 2, a computational method for determining the stress intensity factors at the crack tip is established based on the crack tip elastic stresses. In Section 3, the minimum [J.sub.2] (second invariant variable of the stress deviation) yield function is used to define the variable radius of the plastic core region and crack initiation angle. In Section 4, the crack initiation angle of mixed mode cracks under different loading conditions are simulated, and compared with the simulated results for the MTS and T criteria. Conclusions are presented in Section 5.

2. Stress intensity factors at the crack tip

In this section, the mixed mode cracks are divided into two types: tension cracks and compression-shear cracks. A computational method for determining the stress intensity factors at the crack tip is established based on the crack tip elastic stresses.

2.1. Stress components at the crack tip

The existing research data show that the stress state at a mixed mode crack tip is very complex, depending on the different loading conditions. Mixed mode cracks are divided into two types, combined with the external loading conditions: tension cracks and compression-shear cracks. The crack propagation paths are assumed to be wing cracks. Fig. 1 shows the mixed mode cracks under biaxial loading conditions.

As shown in Fig. 1, a it is assumed that for a crack with an inclination [phi] under biaxial tensile loading conditions, there is no friction effect at the crack surface. As shown in Fig. 1, b a crack with an inclination [phi] under biaxial compressive loading condition and the friction effect at the crack surface is considered in this study. The crack length is 2c, the crack initiation angle is [theta], the axial stress is [[sigma].sub.1] and the confining stress is [[sigma].sub.2]. r is the radius of one element far from the crack tip, and [theta] is the crack initiation angle, the stress components are described by [[sigma].sub.r] and [[sigma].sub.[theta]]. If the external loading is a tensile stress, it has a positive value; and if the external loading is a compressive stress, it has a negative value.

For a crack under a mixed mode condition, simplifying assumptions are used to analyze the stress state at the crack tip. The linear elastic stress field at the crack tip under biaxial loading condition can be expressed as follows (plane stress problem) [3]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where [[sigma].sub.r], [[sigma].sub.[theta]], and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are the radial, circumferential and shear stresses; [K.sub.I] and [K.sub.II] are the stress intensity factors of mode I and mode II cracks.

For a crack in the center of an infinite plane, the crack length is 2c and the crack inclination is p. The stress intensity factors at the crack tip can be represented in a unified form as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)

where a is the external stress, which can be a tensile or compressive stress, and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a function of the different loading conditions and inclinations [phi] of the crack.

The stress intensity factors of mode I and mode II cracks for the plane stress problem are [6]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)

Combining Eqs. (1)-(3), the stress components can be rewritten as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are the function of crack inclination and stress intensity factors.

The stress state at the crack tip is very complex and a plastic core region exists at the crack tip where crack initiation and propagation takes place. In the crack initiation and propagation process, the stress intensity factors are key fracture mechanical parameters. The crack is propagated step by step and results in variation of stress intensity factors at the crack tip.

2.2. Stress intensity factor at the tension crack tip

For the mixed mode crack under tensile stress, as shown in Fig. 1, a the tensile loading conditions can be divided into three types: uniaxial tension; biaxial tension and tensile-compressive stress.

When the crack is under the uniaxial tension condition, there is only axial tensile stress [[sigma].sub.1], and the stress intensity factors can be computed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

where [K.sub.I] is the stress intensity factor of a mode I crack, [K.sub.II] is the stress intensity factor of a mode II crack, [phi] is the inclination of the crack and c is the half length of the crack.

When the crack is under biaxial tension, the axial stress is [[sigma].sub.1], the confining stress is [[sigma].sub.2], and the relationship between them is assumed to be:

[[sigma].sub.2] = [lambda][[sigma].sub.1] (6)

where [lambda] is the ratio of the confining and axial stresses with a range from 0 to 1.

The stress intensity factors under biaxial tension can be computed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (7)

When the crack is under tensile-compressive stress, it is assumed that the axial stress is a tensile stress, while the confining stress is a compressive stress. The relationship between the axial and confining stresses is assumed to be:

[[sigma].sub.2] = -[lambda][[sigma].sub.1], (8)

where [lambda] is the ratio of the confining and axial stresses also with a range from 0 to 1.

The stress intensity factors under tensile-compressive stress can be computed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9)

For the stress intensity factors at the crack tip, we defined a scale factor ^ between mode I and model II stress intensity factors as:

[eta] = [K.sub.I]/[K.sub.II] (10)

Then the scale factor under different external loading condition can be computed as:

[eta] = cos[phi] / sin[phi] (Uniaxial tension) (11, a)

[eta] = ([cos.sup.2] [phi] + [lambda][sin.sup.2] [phi])/[(1 - [lambda]) sin[phi] cos[phi]] (Biaxial tension) (11, b)

[eta] = ([cos.sup.2] [phi] + [lambda][sin.sup.2] [phi])/[(1 - [lambda]) sin[phi] cos[phi]] (Tensile-compressive stress) (11, c)

The scale factor for different external loading condition will be used for the crack initiation analysis process. It has an important impact on the crack propagation path.

2.3. Stress intensity factor at compression-shear crack tip

According to the mixed mode crack under compressive stress, as shown in Fig. 1, b the compressive loading conditions can be divided into two types: uniaxial compression and biaxial compression. Here, the friction effect at the crack surface is considered for the crack initiation and propagation problem for a compression-shear crack.

When the crack is under uniaxial compression, there is only an axial compressive stress [[sigma].sub.1], and the normal and shear stresses at the crack surface are computed as follows [16]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (12)

where [sigma] is the normal stress at the crack surface and t is the shear stress at the crack surface. The effective shear stress is influenced by the mechanical characteristics of the crack surface, can be computed as:

[[tau].sup.eff] = [[sigma].sub.1] sin [phi] con [phi] + [[tau].sub.c] - [mu][[sigma].sub.1] [cos.sup.2] [phi], (13)

where [[tau].sup.c] is the cohesion of the crack surface and f is its friction coefficient.

Here, the cohesion of the crack surface tc is assumed equal to 0 and is not considered for the shear stress at the crack surface, and the effective shear stress is then:

[[tau].sup.eff] = [[sigma].sub.1] sin [phi] con [phi] - [mu][[sigma].sub.1] [cos.sup.2] [phi]. (14)

The stress intensity factors can be computed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

When the crack is under biaxial compression, the axial compressive stress is [[sigma].sub.1] and the confining stress is [[sigma].sub.2], and the normal and shear stresses at the crack surface are computed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

Let the cohesion of the crack surface also be 0, the effective shear stress at the crack surface is:

[[tau].sup.eff] = ([[sigma].sub.1] - [[sigma].sub.2]) sin [phi] con [phi] - [mu] ([[sigma].sub.1][cos.sup.2][phi] + [[sigma].sub.2][sin.sup.2][phi]). (17)

The ratio of confining and axial stresses is introduced so that:

[[sigma].sub.2] = [lambda][[sigma].sub.1], (18)

where [lambda] is the ratio between the confining and axial stresses with the range from 0 to 1.

The stress intensity factors can be computed as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

For a mixed mode crack under compression, the friction effect is considered in this study, and the crack initiation and propagation process is influenced by the shear strength of the crack surface. The crack propagation is limited by the shear strength of the crack surface where the crack closes under compression, but opens under tension. The crack propagates more easily under tension than under compression.

3. Analytical solution for crack initiation of mixed mode crack

For a mixed mode crack under biaxial loading, there is a plastic region around the crack tip and crack propagation will occur with the minimum elastic fracture energy, as shown in Fig. 2. For the plane stress problem, the minimum [J.sub.2] (second invariant variable of the stress deviation) yield function is used to define the variable radius of the plastic core region.

In the plastic core region, crack propagation will occur when the fracture energy is at a minimum. For the crack initiation angle, it is assumed that the minimum angle which will allow the plastic core radius to form is when the minimum [J.sub.2] (second invariant variable of the stress deviation) yield function is used to define the crack initiation angle.

The crack initiation will happen when [J.sub.2] is in the minimum direction, because [J.sub.2] can reflect the real stress state at the crack tip, and the volume strain energy density is larger than the shape change energy density. The crack initiation angle is determined by:

[partial derivative][J.sub.2]/[partial derivative][theta] = 0, (20)

where the second invariant variable of the stress deviation can be computed as follow,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (21)

Moreover, the second invariant variable of the stress deviation [J.sub.2] should be at the minimum value:

[[partial derivative].sup.2][J.sub.2]/[partial derivative][[theta].sub.2] 0 and [J.sub.2min] = [J.sub.2] ([[theta].sub.0]). (22)

Combining Eqs. (20) and (21), we then obtain:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)

The crack initiation angle [[theta].sub.0] can be solved by Eq. (23), and the computed results for the crack initiation angle should also satisfy Eq. (22).

If we know the stress intensity factors at the tip of mixed mode crack, the crack initiation angle can be determined based on the failure criterion of solid material. The presented method can compute the crack initiation angle for the simple mode I or mode II crack problem.

When the crack is a mode I fracture, the stress intensity factor of [K.sub.II] is equal to 0, and the Eq. (23) can be simplified as:

3[K.sup.2.sub.I]sin(2[theta]) - 3[K.sup.2.sub.I] sin(2[theta]) = 0. (24)

The solution results of Eq. (32) are:

[[theta].sup.0] = 0 and [[theta].sup.0] = [+ or -] arccos 1/3. (25)

The crack initiation angle [[theta].sub.0] should also satisfy Eq. (22), then:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (26)

So the crack initiation angle for the mode I fracture is:

[[theta].sub.0] = 0 (27)

The analysis result is the same for several classic fracture criterions, such as MTS, T and other criteria.

When the crack is a mode II fracture, the stress intensity factor of [K.sub.I] is equal to 0, and Eq. (23) can then be simplified as:

2[K.sup.2.sub.II]sin[theta] - 9[K.sup.2.sub.II] sin(2[theta]) = 0. (28)

The solution results for Eq. (36) are:

[[theta].sub.0] = 0 and [[theta].sub.0] = [+ or -] arccos 1/2. (29)

The crack initiation angle [[theta].sub.0] should also satisfy Eq. (22), then:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (30)

So the crack initiation angle for the mode II fracture is:

[[theta].sub.0] = [+ or -] 83.602 6[degrees]. (31)

The crack initiation angle of the mode II fracture is 70.53[degrees] computed using the maximum circumferential stress fracture criterion, and the crack initiation angle of the mode II fracture is -70.53[degrees] computed using the maximum tensile stress fracture criterion II. The simulation results computed by this method are larger than for other classic fracture criteria, but the computing error is very small.

4. Simulated results of crack initiation angle

In this section, the crack initiation angle of a mixed mode crack under different loading conditions is simulated, and compared with the simulated results for the MTS and T criteria. The crack initiation angle influenced by the loading ratio between the confining and axial stresses, and the friction effect at the crack surface under compression, are analyzed. The crack propagation characteristics under different loading conditions are presented in the end.

4.1. Tension crack

According to the mixed mode crack under tensile stress, for the plane stress problem, the scale factor n between mode I and mode II stress intensity factors is introduced in Eq. (23), giving:

12[eta]cos(2[theta]) + 3([[eta].sup.2] - 3)sin(2[theta]) -

-4[eta]cos[theta] - 2([[eta].sup.2] - 1)sin[theta] = 0. (32)

For the mixed mode crack under uniaxial tensile loading, the loading ratio [lambda] is 0 and Eq. (11, a) is substituted into Eq. (32), and a function for the crack inclination and crack initiation angle can then be obtained. Table 1 shows the computed results for the crack initiation angle for different crack inclinations under uniaxial tensile loading.

Fig. 3 shows the crack initiation angle for a mixed mode crack under uniaxial tensile loading using different fracture criteria.

As shown in Fig. 3, the computed results for the crack initiation angle under uniaxial tensile loading by the presented criterion is less than the T-criterion and larger than the MTS-criterion. The crack initiation angle increases with increasing crack inclination. The relationship between the crack initiation angle and crack inclination is a nonlinear function, and can be simplified into a linear relationship as:

[[theta].sub.0] + [phi] [congruent to] 0[degrees], (33)

where [phi] is the crack inclination and [[theta].sub.0] is the crack initiation angle.

From Fig. 1, a we observe that for the mixed mode crack under uniaxial tensile loading, the crack initiation angle is perpendicular to the axial loading direction.

For the mixed mode crack under biaxial tensile loading, the loading ratio [lambda] is given as 0.2, 0.5 and 0.8, and substituting Eq. (11, b) into Eq. (32), we then obtain a function for the crack inclination and crack initiation angle. Fig. 4, a shows the compute results for the crack initiation angle for different crack inclinations under biaxial tensile loading when the loading ratio is 0.5. Fig. 4, b shows the

As shown in Fig. 4, a for the mixed mode crack under the loading ratio of 0.5 for biaxial tensile stress, the crack initiation angle first increases and then decreases with the increasing of crack inclination, and the error between the presented criterion and the MTS and T criteria is very small. The crack initiation angle is influenced by the external loading pattern. As shown in Fig. 4, b the crack initiation angle decreases with increasing loading ratio, and the crack propagation is restrained by the lateral tensile stress. The crack is easy to propagate under biaxial tensile loading, while with increasing lateral tensile stress, the crack propagation is progressively more difficult. If the lateral tensile stress is equal to the axial tensile stress, the crack ceases to propagate.

For the mixed mode crack under tensile-compressive loading, the loading ratio [lambda] is set at 0.2, 0.5 and 0.8, and substituting Eq. (11, c) into Eq. (32), we then obtain a function for the crack inclination and the crack initiation angle. Fig. 5, a shows the compute results for the crack initiation angle for different crack inclinations under tensile-compressive loading when the loading ratio is 0.5. Fig. 5, b shows the compute results for the crack initiation angle for different crack inclinations under different tensile-compressive loading conditions.

As shown in Fig. 5, a for the mixed mode crack under a loading ratio of 0.5 for the tensile-compressive stress condition, the crack initiation angle increases with increasing crack inclination. The error between the presented criterion and the MTS and T criteria is very small. As shown in Fig. 5, b the crack initiation angle increases with increasing loading ratio, and the crack propagation is accelerated by the lateral compressive stress. With increasing lateral tensile compression, the crack propagation accelerates.

4.2. Compression-shear crack

For the mixed mode crack under a compressive loading condition, the fracture type is a compressive-shear crack. A scale factor [eta] between mode I and mode II stress intensity factor is defined for the compressive-shear crack:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (34)

When the compressive-shear cracks under uniaxial compression, the stress intensity factors at the crack tip are:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (35)

When the compressive-shear cracks under biaxial compression, the stress intensity factors at the crack tip are:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (36)

For the mixed mode crack under uniaxial compression, substituting Eq. (35) into Eq. (32), we then obtain a function for the crack inclination and crack initiation angle. Fig. 6, a shows the compute results for the crack initiation angle for different crack inclinations under uniaxial compression when the friction coefficient of the crack surface is 0.2. Fig. 6, b shows the compute results for the crack initiation angle for different crack inclinations under different friction coefficients for the crack surface.

As shown in Fig. 6, a there is a great difference between different fracture criteria for mixed mode cracks, and the crack initiation angle decreases with increasing crack inclination, following a nonlinear relationship. As shown in Fig. 6, b the crack initiation angle is obviously influenced by the crack friction coefficient. The crack surface friction prevents the crack propagating, and it increases as the crack friction coefficient increases while the other parameters remain constant.

The relationship between the crack initiation angle and crack inclination can be described by a linear function:

[[theta].sub.0] + [phi] [congruent to] 180[degrees]. (37)

For the mixed mode crack under biaxial compression, substituting Eq. (36) into Eq. (32), we then obtain a function for the crack inclination and the crack initiation angle. Fig. 7, a shows the computed results for the crack initiation angle for different crack inclinations under biaxial compression when the friction coefficient of the crack surface is 0.2, and the loading ratio is also 0.2. Fig. 7, b shows the computed results for the crack initiation angle for different crack inclinations under different loading ratios, [lambda] of 0.2, 0.5 and 0.8.

As shown in Fig. 7, a for the mixed mode crack under a loading ration of 0.2 for the biaxial compression condition, the crack initiation angle first decreases and then increases with the increasing crack inclination. The crack initiation angle is influenced by the external loading pattern. As shown in Fig. 7, b, the crack initiation angle decreases with increasing loading ratio, and the crack propagation is restrained by the lateral compressive stress. The crack propagation is limited by the lateral stress under biaxial compression.

When the crack inclination is less than about 70[degrees], the relationship between the crack initiation angle and crack inclination can be described by a linear function:

[[theta].sub.0] + [phi] [congruent to] 90[degrees]. (38)

The crack initiation angle of a mixed mode crack is influenced by the external loading pattern, crack inclination, and the crack friction effect.

5. Conclusions

Cracks propagate forming wing cracks under biaxial loading conditions, and the study of mixed mode crack initiation and propagation is very important for understanding the mechanical characteristics of solids. The stress state at the crack tip is very complex. A plastic core region exists at the crack tip, and the crack initiation and propagation will occur in this region. In the crack initiation and propagation process, the stress intensity factors are key fracture mechanical parameters. A computational method for determining the stress intensity factors at the crack tip is established based on the crack tip elastic stresses. The stress intensity factors at the tension crack tip are computed with three tensile loading conditions: a. uniaxial tension; b. biaxial tension; c. tensile-compressive stress. The stress intensity factors at the compressive-shear crack tip are computed with two tensile loading conditions: a. uniaxial compression; b. biaxial compression. The friction effect at the crack surface is considered for the crack initiation problem for the compression-shear crack.

A plastic core region exists at the crack tip, the minimum [J.sub.2] (second invariant variable of the stress deviation) yield function is used to define the variable radius of the plastic core region. The crack initiation angle of a mixed mode crack under different loading condition was simulated. The computed results show that the crack initiation angle under uniaxial tensile loading using the presented criterion is less than for the T criterion and larger than for the MTS criterion. The crack initiation angle is influenced by the loading ratio between the confining and axial stresses, and the friction effect at the crack surface under compression.

Acknowledgements

This work is supported by the National Natural Science Foundation of China (No. 41102194 and No. 51209156) and the China Postdoctoral Science Foundation (No. 2012T50785).

Received April 25, 2012

Accepted September 05, 2013

References

[1.] Sih, G.C. 2008. Modeling of multiscale fatigue crack growth: nano/micro and micro/macro transitions, in: B. Farahmand (Ed.),Virtual Testing and Predictive Modeling: Fatigue and Fracture Mechanics Allowables, Springer, 187-219.

[2.] Lawn, B.R.; Wilshaw. T.R. 1975. Fracture of brittle solids, London: Cambridge University Press.

[3.] Sih, G.C. 2009. Crack tip mechanics based on progressive damage of arrow: hierarchy of singularities and multiscale segments, Journal of Theoretical and Applied Fracture Mechanics 51(1): 11-32. http://dx.doi.org/10.1016/j.tafmec.2009.01.007.

[4.] Nuismer, R.J. 1975. An energy release rate criterion for mixed mode fracture, International Journal of Fracture Mechanics 11(2): 245-250. http://dx.doi.org/10.1007/BF00038891.

[5.] Erdogan, F.; Sih, G.C. 1963. On the crack extension in plates under plane loading and transverse shear, Journal of Basic Engineering 85(4): 519-527. http://dx.doi.org/10.1115/1.3656897.

[6.] Ziliukas, A.; Surantas, A.; Ziogas, G. 2010. Strength and fracture criteria application in stress concentrators areas, Mechanika 3(83): 17-20.

[7.] Theocaris, P.S.; Kardomateas, G.A.; Andrianopoulos, N.P. 1982. Experimental study of the T-criterion in ductile fracture, Engineering Fracture Mechanics 17(5): 439-447. http://dx.doi.org/10.1016/0013-7944(83)90040-1.

[8.] Kong, X.M.; Schluter, N.; Dahl, W. 1995. Effect of triaxial stress on mixed-mode fracture, Engineering Fracture Mechanics 52(2): 379-388. http://dx.doi.org/10.1016/0013-7944(94)00228-A.

[9.] Zhou, J.W.; Xu, W.Y.; Yang, X.G. 2010. A mi crocrack damage model for brittle rocks under uniaxial compression, Mechanics Research Communications 37(4): 399-405. http://dx.doi.org/10.1016/j.mechrescom.2010.05.001.

[10.] Khan, S.M.A.; Khraisheh, M.K 2000. Analysis of mixed mode crack initiation angles under various loading conditions, Engineering Fracture Mechanics 67(5): 397419. http://dx.doi.org/10.1016/S0013-7944(00)00068-0.

[11.] Shariati, M.; Sedighi, M.; Saemi, J.; Eipakchi, H.R.; Allahbakhsh, H.R. 2010. Numerical and experimental investigation on ultimate strength of cracked cylindrical shells subjected to combined loading, Mechanika 4(84): 12-19.

[12.] Sih, G.C.; Tang, X.S. 2004. Dual scaling damage model associated with weak singularity for macroscopic crack possessing a micro/mesoscopic notch tip, Journal of Theoretical and Applied Fracture Mechanics 42(1): 1-24. http://dx.doi.org/10.1016/j.tafmec.2004.06.001.

[13.] Marie, S.; Chapuliot, S. 2002. 2D crack growth simulation with an energetic approach, Nuclear Engineering and Design 212(1-3): 31-40. http://dx.doi.org/10.1016/S0029-5493(01)00475-7.

[14.] Khan, S.M.A.; Khraisheh, M.K. 2004. A new criterion for mixed mode fracture initiation based on the crack tip plastic core region, International Journal of Plasticity 20(1): 55-84. http://dx.doi.org/10.1016/S0749-6419(03)00011-1.

[15.] Zhang, Y.; Ji, H.G.; Ye, J.H.; Li, S.Y. 2012. Analytical solutions of the tensile strength and preponderant crack angle for the I-II mixed crack in brittle material, Mechanika 18(1): 22-28.

[16.] Zhou, J.W.; Yang, X.G.; Li, H.T.; Zhou, H.W. 2011. A plasticity criterion for mixed mode fracture initiation with crack surface frictional characteristics under compression, The Open Mechanics Journal 5: 1-6. http://dx.doi.org/10.2174/1874158401105010001.

Jiawen Zhou, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China, E-mail: jwzhou@scu.edu.cn

Xingguo Yang, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China, E-mail: 89022251@163.com

Huige Xing, College of Architecture and Environment, Sichuan University, Chengdu 610065, China, E-mail: hgxing@scu.edu.cn

http://dx.doi.org/ 10.5755/j01.mech.19.5.5530

Table 1 Compute results of the crack initiation angle for different crack inclinations under uniaxial tensile loading No. 1 2 3 4 Crack inclination, [degrees] 0.000 10.000 20.000 30.000 Crack initiation angle, [degrees] 0.000 -18.486 -31.798 -41.490 [K.sub.I]/[K.sub.II] none 5.6713 2.7475 1.7321 No. 5 6 7 Crack inclination, [degrees] 40.000 50.000 60.000 Crack initiation angle, [degrees] -49.470 -59.660 -63.486 [K.sub.I]/[K.sub.II] 1.1918 0.8391 0.5774 No. 8 9 10 Crack inclination, [degrees] 70.000 80.000 90.000 Crack initiation angle, [degrees] -69.839 -76.857 -83.621 [K.sub.I]/[K.sub.II] 0.3740 0.1763 0.0000

Printer friendly Cite/link Email Feedback | |

Author: | Zhou, Jiawen; Yang, Xingguo; Xing, Huige |
---|---|

Publication: | Mechanika |

Article Type: | Report |

Geographic Code: | 9CHIN |

Date: | Sep 1, 2013 |

Words: | 4819 |

Previous Article: | Steel abrasive wear forecasting by wearing surfaces microgeometric parameters/Plienu abrazyvinio dilimo prognozavimas pagal dylancio pavirsiaus... |

Next Article: | A new non-linear continuum damage mechanics model for fatigue life prediction under variable loading/Naujas netiesinio koontiniumo pazeidimo... |

Topics: |