УДК 539.42
Оценка разрушения полиметилметакрилата на примере образцов в форме диска с острым надрезом при нагружении смешанного типа I + III
M.R.M. Aliha1, F. Berto23, A. Bahmani1, Sh. Akhondi1, A. Barnoush3
1 Иранский университет науки и технологии, Тегеран, 16846-13114, Иран 2 Падуанский университет, Виченца, 36100, Италия 3 Норвежский университет естественных и технических наук, Тронхейм, 7491, Норвегия
В статье экспериментально и теоретически исследовано поведение акрилового стекла Perspex (полиметилметакрилата) при нагружении смешанного типа I/III с использованием новой упрощенной схемы лабораторных испытаний. В качестве образца использовали круглый диск с V-образным надрезом вдоль диаметра, нагружаемый по схеме трехточечного изгиба. Получены критические значения коэффициента интенсивности напряжений для надреза (KjV, Кщ) для ряда комбинаций режимов I и III при изменении угла наклона надреза относительно нагружающих цилиндров. Показано, что вязкость разрушения в вершинах надреза в условиях чистого нагружения или преимущественно растягивающей нагрузки выше соответствующего значения при нагружении типа III или в условиях преимущественно скручивающей нагрузки. Экспериментальные результаты с высокой точностью описываются критерием локальной плотности энергии деформации.
Ключевые слова: полиметилметакрилат, нагружение смешанного типа I/III, сопротивление распространению надреза, критерий плотности энергии деформации
Fracture assessment of polymethyl methacrylate using sharp notched disc bend specimens under mixed mode I + III loading
M.R.M. Aliha1, F. Berto23, A. Bahmani1, Sh. Akhondi1, and A. Barnoush3
1 Welding and Joining Research Center, School of Industrial Engineering, Iran University of Science and Technology, Narmak, Tehran, 16846-13114, Iran 2 Department of Management and Engineering, University of Padova, Vicenza, 36100, Italy
3 Department of Engineering Design and Materials, Norwegian University of Science and Technology, Trondheim, 7491, Norway
Mixed mode I/III behavior of Perspex (polymethyl methacrylate (PMMA)) is studied experimentally and theoretically in this research using a new and simple laboratory test configuration. The specimen is a circular disc containing a sharp V-notch along the diameter that is loaded by the conventional three-point bend fixture. The critical values of notch stress intensity factors (KjV and Кщ) were obtained for the whole combinations of modes I and III simply by changing the notch inclination angle relative to the loading rollers. The value of notch fracture toughness under pure or dominantly tension loads was greater than its corresponding value under mode III or dominantly torsion loads. The experimental results were also predicted very well by employing the local strain energy density (SED) criterion.
Keywords: PMMA, mixed mode I/III loading, notch fracture toughness, SED criterion
1. Introduction
Integrity assessment of engineering components and structures in the presence of high stress concentration sources such as notches or sharp cracks, is an interesting subject and of practical importance. Different aspects associated with fracture of cracked or notched components have been extensively investigated in the past by means of
laboratory test specimens under controlled testing conditions and by means of model materials. Polymethyl methacrylate (PMMA) is a well-known model material for investigating brittle fracture phenomenon. Optical transparency of PMMA makes ease study of fracture initiation and its propagation trajectory under any desired complex loading situation. Fracture test specimen in any arbitrary shapes
© Aliha M.R.M., Berto F., Bahmani A., Akhondi Sh., Barnoush A., 2016
can be easily manufactured from a PMMA sheet by laser cutting or water jet machines. Fracturing load of such material with laboratory scale samples is relatively low and hence fracture toughness experiments on PMMA can be conducted even by means of conventional and low capacity testing machines. Consequently, a large number of researchers have used this material in their experimental fracture studies. Based on the principles of fracture mechanics a cracked or notched body can experience mode I (opening or tensile load), mode II (in-plane sliding or shear deformation), mode III (out-of-plane sliding or tearing deformation) or any combinations of these three basic modes. Under each loading condition, study of the onset of fracture initiation, the direction of fracture growth and the path of brittle fracture are interesting subjects and hence PMMA has been widely used in the past for investigating the mentioned problems. Using different test specimens such as (center cracked tensile plate [1], single edge notched bend beam [2, 3], diametral compressed Brazilian disc [4], semicircular bend [5-7], triangular specimen subjected to three-point bending [8], diagonally loaded square plate [9], compact tension-shear specimen [10] and etc.), mode I, mode II and mixed mode I/II fracture behavior of PMMA has been studied experimentally and extensive data are available in the literature for K Ic and K IIc of cracked PMMA samples. Furthermore, there are also several research papers dealing with the fracture behavior of notched PMMA specimens under different loading modes. For example, mode I and mixed mode I/II notch fracture toughness of PMMA has been determined using some test configurations such as edge notched three point bend beam and compact tension specimens and by considering different U- and V-notch geometries [11-14]. However, in comparison with the mixed mode I/II loading case, there are relatively fewer experimental results presented for mode III or mixed mode
I/III fracture behavior of PMMA [15, 16]. In other words, only a limited number of test specimens are available for mixed mode I/III studies. However, the available test specimens have two major limitations: (i) many of them are not able to introduce full combinations of modes I and III [17, 18], (ii) some of them have complex geometries and require difficult preparation process and complex testing rigs or fixtures [19, 20]. In order to overcome these two drawbacks, in this paper a simple and suitable disc type bend test configuration is proposed for mixed mode I/III fracture toughness studies. The applicability of the suggested notched specimen is then examined by performing some fracture tests on PMMA. It is shown that the suggested specimen can easily introduce complete range of modes, including mixed ones, from pure mode I to pure mode III. Furthermore, it is demonstrated that the notch fracture toughness test data obtained from the tested PMMA can be predicted by means of the strain energy density (SED) criterion [21, 22].
2. Suggested mixed mode I/III fracture specimen
As shown in Fig. 1, the proposed specimen is a circular disc of diameter D and height B which contains a sharp edge notch of angle 2a and depth a along one side of disc. The specimen is subjected to three-point bend loading with loading span of 2S, in which the bottom loading supports are symmetric with respect to the top roller. Hence, the specimen is so called edge notched disc bend (ENDB). The state of notch flank deformation in the ENDB specimen can be controlled by the angle P, defined in Fig. 1. When P is zero (i.e. symmetric ENDB shown in Figs. 1, a, c), the specimen would be subjected to pure mode I due to symmetry of both geometry and loading relative to the notch plane. By increasing P, the specimen experiences mixed mode opening and out-of-plane sliding deformation.
s —I
Mid-section z = D/2
\ N1 D
i "1
z = D
\>-j • S---
0
V
Fig. 1. Geometry and loading condition of the ENDB specimen containing a sharp notch of angle 2a for introducing pure mode I, pure mode III and mixed mode I/III conditions
Recently Aliha and coworkers [23, 24], studied numerically this specimen containing a sharp crack and computed its stress intensity factors for different a/B, D/2S and P values. They showed that the mode II component has negligible value relative to KT and KIIT in this specimen and the full combinations of mode I/III mixities can be achieved by changing the crack inclination angle. However, in practice the artificial notches created by means of razor or thin saw blades in laboratory fracture study of brittle materials like PMMA are not real sharp cracks and hence in order to provide more realistic estimations for fracture behavior of such test specimens and investigating the possible effects of notch type, it is better to consider the created sluts as sharp notch instead of real crack. Hence in the following, by assuming a sharp notch in the ENDB specimen, the notch stress intensity factors (NSIFs) are computed for this specimen using finite element analyses and for different mixed mode I/III loading conditions.
modes I, II and III notch stress intensity factors and can be determined from the following equations:
K V = lim[(2 nrl~X t) g00 ], (6)
r
kV = lim[(2 nrl~X "k0 ], (7)
r ^0
K V = lim[(2nr1-X Zr ]. (8)
r ^0
The values of Xt, Xn and Xm can also be determined by solving the below equations [30]:
sin (2Xy) + Xsin (2y) = 0 (mode I), (9)
sin(2Xy) - Xsin(2y) = 0 (mode II), (10)
cos(Xy) = 0 (mode III). (11)
Corresponding values of X t, X it and X m for small notch opening angles (up to 2a = 15°) were determined as 0.5 from Eqs. (9) to (11), implying that singularity character of sharp notches is the same as cracks. Thus as a good engineering assumption such notches can be considered as a real crack.
3. NSIF calculation for the ENDB specimen
Based on the William's expression, the stress field for a sharp notch in a polar coordinate (as defined in Fig. 1) can be written as following equations extracted from [25-28]:
CTrr = frr (0) + (0)'
r i r
K v k v
a00 = ^i-XT f00(0) + /-XT g 00(0),
°r0 = fr0 (0) + -rk gr0 (0),
r 1 r II
KIII
,1-Xttt
sin(X tii 0),
a z0 =
-^cos( XIII0),
(1) (2)
(3)
(4)
(5)
in which Xt, Xit and Xiit are the eigenvalues that depend on the notch angle 2a and f, g are angular functions given in [29]. The parameters kV (i = I, II and III) are the
Fig. 2. Finite element model of the ENDB specimen and a zoomed view of notch tip
0.16
8 30
ISI I
E
10
b b
-10
V 0
\3
\
\
\
\2 s
0.16
».00 0.04
0.08
r, mm
0.16
Fig. 3. Variations of stress components along the notch line and for different through thickness locations (z = D/ 2 (1), 3 D/ 4 (2), D (3))
Meanwhile, g00 , tr0 and tzr in Eqs. (6) to (8) are the polar stress components on the notch bisector line (i.e. along 0 = 0°). Therefore, in order to compute the NSIFs for the investigated ENDB specimen, it is necessary to obtain these three stress components along the nodes located at 0 = 0° line. By performing a static finite element analysis for each loading angle P, the required stress components were obtained numerically. Following geometrical parameters were considered for the finite element modeling of the sharp notched = ENDB specimen: D = 70 mm, B = 20 mm, a = = 8 mm, 2S = 66.5 mm and 2a = 5°. The material properties of PMMA was also used in finite element models as v = 0.3 and E = 2950 MPa. Figure 2 shows the finite element model of the ENDB specimen and a zoomed view of the sharp notch tip. The model was created by means of 65000 solid C3D20 elements in the ABAQUS code, in which very fine elements were employed along the notch front.
The variations g00 , Gr0 and Gzr stress components along the notch bisector have been presented in Fig. 3 for typical P angle of 53o and for different crack front locations defined by z in Fig. 1. The curves in these graphs have been normalized to the corresponding value of stress component obtained at the mid-section (i.e. z = D/2) and at r = 0. It is seen from this figure that the notch tip stress field depends noticeably on the front location and in general a complex 3D stress state is seen at the tip of notch. But as one moves along the ligament, most of the stress components tends towards zero. Figure 4 shows the results of modes I, II and III NSIFs along the crack front computed from Eqs. (6) to (8) for different crack inclination angles. The presented results have been normalized by dividing the NSIFs to the KIV at 0° (i.e. the value of KI component at the mid-section and P = 0 °). According to this figure, in general, three-
dimensional effects are observed for the distribution of NSIFs along the notch front. By ignoring the unreliable results for the corner edges of the specimen due to edge effects (i.e. typically for z < 7 mm and z > 63 mm), it is seen that the mode I and mode III NSIFs increase gradually by moving towards the mid-section but mode II component (i.e. KV) tends to zero. Thus, at the mid-section of specimen (i.e. z = 35 mm), both KIV and K^ components become maximum but KV is zero for any crack inclination angle P. In other words, although both modes II and III shear components are always coupled in a cracked or notched specimen [31-33], but the value of mode II component becomes negligible in the mid-section zone of the ENDB specimen. Therefore, the ENDB specimen at the onset of fracture would experience mixed mode I/III loading in this critical point (which corresponds to the point of fracture initiation). When P is zero, it is seen from Fig. 4, a that both KV and K^ components are zero but KV is nonzero that indicates the ENDB specimen is subjected to pure mode I loading. Conversely, at a specific angle P (for instance P = 65 ° in the case of ENDB of this research) which depends on the specimen's geometry and notch length, KI becomes zero at the mid-section but K^ is nonzero (Fig. 4, d). This situation corresponds to pure mode III loading. Any inclination angle between zero and PIII (i.e. pure mode III inclination angle) may introduce an intermediate mode I/III mixity. As a conclusion, the ENDB specimen can be used for producing the whole mode mixities from pure mode I to pure mode III notch fracture modes simply by setting the notch inclination angle P in appropriate value relative to the vertical direction. In the next sections, the practical applicability of the suggested notched ENDB specimen is investigated via some fracture experiments on PMMA and also theoretical studies.
fe
Xfl £
1.0
0.6
0.2
o
& L5
Xfl
* 0.5
<D
I-0.5 o
S-1.5
a
1^ Mid 2---------- ■section
20
40
z, mm
60
80
, ^
----
Midi section i
20
40
z, mm
60
80
Fig. 4. Variations of normalized notch stress intensity factors (KV, KV and KjjI) for the investigated notched ENDB specimen thorough the crack front. KjKIm (1), Kn/KIm (2) and Km/KIm (3) at 0°
Fig. 5. Test setup (a), typical load-displacement (b) and fracture trajectories of the ENDB specimen made of PMMA (c-e): pure mode I (c), mixed mode I/III (d), fracture surface (e)
4. Notch fracture toughness experiments on PMMA
In order to conduct the experiments, several disc specimens with radius of 35 mm and thickness of 20 mm were cut from a PMMA sheet using laser cutting machine. Then an edge notch of 8 mm depth were cut in one side of each disc by means of a very narrow rotary saw blade. The loading span (i.e. 2S) was considered equal to 64.75 mm. The manufactured ENDB specimens were then placed inside a three-point bend fixture as shown typically in Fig. 5, a with different inclination angles P (in the range of zero and 65°) and were loaded monotonically until the fracture of specimens. All the tested ENDB samples had linear and elastic behavior and showed different trajectories depending on the mode I and III mixity. Figures 5, b, c show typical load-displacement curves and fracture trajectories observed for
the tested ENDB specimens made of PMMA under mixed mode I and III. It is seen from Fig. 5, c that the pure mode I fracture path is along the plane of initial notch, but by increasing the contribution of mode III component, fracture extends along out of the plane of initial notch and its deviation from the initial notch plane becomes more by moving toward pure mode III.
The notch fracture toughness experiments were conducted for six P angle of 0°, 27°, 44°, 53°, 60° and 65° and for each angle three replicates of ENDB specimen made of PMMA were tested. By applying the critical fracture load (i.e. the peak load) of each specimen to the finite element models of the ENDB specimen, the corresponding values of polar stresses in front of the notch was obtained from the finite element analyses. Then, the critical values of NSIFs
Table 1
Average of mixed mode I/III notch fracture toughness results for the PMMA tested with ENDB specimen
P
27°
44°
53°
60°
65°
Average of P, N
2784.50
3512.75
6267.00
6856.00
8636.75
9181.25
K V,
MPa • m
12
1.92
2.00
1.73
1.03
0.50
0.00
K
II :
MPa • m
12
0
KIII,
MPa • m1/2
0.00
0.55
1.37
1.31
1.50
1.40
at the onset of fracture were obtained from Eqs. (6) to (8). Table 1 presents the corresponding values of critical notch stress intensity factors for different notch inclination angles in the investigated ENDB specimens. This table obviously reveals that the ENDB specimen is able to introduce full combinations of mode I/III mixities from pure mode I to pure mode III in PMMA when P changes from zero to 65°. Mode II component was also nearly zero in all mode mixities. By increasing P, the value of KTV decreases and conversely K^ increases. However, the value of critical NSIF in this specimen under pure mode I is noticeably greater than its corresponding value under pure mode III loading case. In the next section the experimental results are predicted theoretically using an energy-based fracture criterion.
5. Strain energy density averaged over a control volume. The fracture criterion
With the aim to assess the fracture load in notched PMMA components, an appropriate fracture criterion is required which has to be based on the mechanical behavior of material around the notch tip. In this section, a criterion proposed by Lazzarin and co-authors [21, 22] based on the strain energy density is briefly described.
The averaged strain energy density criterion as presented in Refs. [22, 34] states that brittle failure occurs when the mean value of the strain energy density over a given con-
trol volume is equal to a critical value Wc. This critical value varies from material to material but it does not depend on the notch geometry and sharpness. The control volume is considered to be dependent on the ultimate tensile strength at and the fracture toughness KIc in the case of brittle or quasi-brittle materials subjected to static loads.
The method based on the averaged SED was formalized and applied first to sharp (zero radius) V-notches under mode I and mixed mode I/II loading [21] and later extended to blunt U- and V-notches [22, 35, 36]. When dealing with cracks, the control volume is a circle of radius Rc centered at the crack tip (Fig. 6, a). Under plane strain conditions, the radius Rc can be evaluated according to the following expression:
RIc =
(1 + v)(5 - 8v)
4n
K
(12)
where KIc is the fracture toughness, v is the Poisson's ratio, and a t is the ultimate tensile stress of a plane specimen.
For a sharp V-notch, the critical volume becomes a circular sector of radius Rc centred at the notch tip (Fig. 6, b). When only failure data from open V-notches are available, Rc can be determined on the basis of some relationships reported in [22], where KIc is substituted by the critical value of the notch stress intensity factors as determined at failure from sharp V-notches.
Dealing here with sharp notches under torsion loading, the control radius R3c can be estimated by means of the following equation [31]:
\1/(1-X 3)
R3c =
e3
K
IIIc
1 + v T
(13)
where KIIIc is the mode III critical notch stress intensity factor and tt is the ultimate torsion strength of the un-notched material. Moreover, e3 is the parameter that quantifies the influence of all stresses and strains over the control volume, and (1 -X3) is the degree of singularity of the linear elastic stress fields [29], which depends on the notch opening angle. Values of e3 and X3 are 0.4138 and 0.5 for the crack or very sharp notch cases (2a ~ 0°).
o
x (n- 2a)/(2n- 2a)
6. SED approach in fracture analysis of the tested PMMA specimens
The fracture criterion described in the previous section is employed here to estimate the fracture loads obtained from the experiments conducted on the PMMA specimens.
As originally thought for pure modes of loading the averaged strain energy density criterion states that failure occurs when the mean value of the strain energy density over a control volume W reaches a critical value Wc, which depends on the material but not on the notch geometry.
Under tension loads, this critical value can be determined from the ultimate tensile strength at according to Beltrami's expression for the unnotched material: _2
Wic = (14)
2 E
By using the values of at = 70 MPa and E = 2950 MPa, the critical SED for the tested PMMA is W1c = = 0.83 MJ/m3.
Under torsion loads, this critical value can be determined from the ultimate shear strength Tt according to Beltrami's
expression for the unnotched material:
w'c=<i5)
By using the values of Tt = 43 MPa and G = 1100 MPa, the critical SED for the tested graphite is W3c = 0.84 MJ/m3.
In parallel, the control volume definition via the control radius Rc needs the knowledge of the mode I and mode III critical notch stress intensity factor KV and KVIc and the Poisson's ratio v, see Eqs. (12) and (13). For the considered material KV and KVIc have been obtained from specimens weakened by sharp V-notches with an opening angle 2a < 5° and a notch radius less than 0.1 mm. A pre-crack was also generated with a razor blade at the notch tip. The resulting values are KV = 1.92 MPa-m05 and KVIc = = 1.40 MPa-m05 which provide the control radii R1c = = 0.183 mm and R3c = 0.324 mm, under pure tension and pure torsion, respectively.
The approach proposed here is a reminiscent of the work by Gough and Pollard [32] who proposed a stress-based expression able to summarize together the results obtained from bending and torsion. The criterion was extended in
Table 2
Experimental results. Overview of the data by using the proposed SED method
M W1, M j/m3 W3, M j/m3 WJ Wl0 + W3/W3c (Wi/Wic + W3IW30)12
0.00 0.83 0.00 1.00 1.00
0.00 0.89 0.00 1.07 1.03
0.00 0.94 0.00 1.14 1.07
0.00 0.91 0.00 1.09 1.05
0.19 0.66 0.12 0.94 0.97
0.19 0.73 0.13 1.04 1.02
0.19 0.81 0.15 1.15 1.07
0.19 0.71 0.13 1.01 1.00
0.48 0.33 0.57 1.07 1.04
0.48 0.35 0.62 1.15 1.07
0.48 0.38 0.66 1.23 1.11
0.48 0.36 0.63 1.18 1.09
0.51 0.37 0.75 1.33 1.16
0.51 0.39 0.80 1.41 1.19
0.51 0.41 0.84 1.49 1.22
0.51 0.38 0.77 1.37 1.17
0.78 0.06 0.93 1.17 1.08
0.78 0.06 0.93 1.17 1.08
0.78 0.06 0.96 1.20 1.10
0.78 0.06 0.95 1.19 1.09
1.00 0.00 0.92 1.08 1.04
1.00 0.00 0.94 1.10 1.05
1.00 0.00 0.97 1.14 1.07
1.00 0.00 0.95 1.12 1.06
terms of the local SED to V-notches under fatigue loading in the presence of combined tension and torsion [33] and recently under static loading [37].
In agreement with [33, 37] and extending the method to the static case, the following elliptic expression:
W_+W=1
W1C W3C
(16)
is obtained. In Eqs. (14) and (15) W1c and W3c are the critical values of SED under pure tension and pure torsion. For the considered graphite, W1c = 0.83 MJ/m3 and W3c = = 0. 84 MJ/m3. Each specimen reaches its critical energy when the sum of the weighted contributions of mode I and mode III is equal to 1, which represents the complete damage of the specimen.
The values of W1 and W3 have, instead, to be calculated as a function of the geometry and of the applied mode mixity ratio, where KI and KIII represent the mode I and mode III SIF ranges, R1c and R3c are the radii of the control volume related to mode I and mode III loadings, while e1 and e3 are equal to 0.11861 and 0.4138, respectively.
The detailed calculations employing this second criterion are reported in Table 2. The square root of the left-hand side term of Eq. (16), which is, in fact, proportional to the critical load, is given in the last column of this table.
A synthesis in terms of the square root value of the considered parameter, that is the sum of the weighted energy contributions related to mode I and mode III loading, is shown in Fig. 7 as a function of the mode mixity parameter M = 21 n arctg( Km/ KI). Many of the results are inside a scatter band ranging from 0.9 to 1.1 with only few exceptions.
The fracture models proposed in this paper can be used for predicting the onset of brittle fracture in notched PMMA components which are subjected to a combination of tension and torsion loading.
7. Conclusions
Pure mode I, pure mode III and different intermediate mode I/III mixities fracture toughness values of commercial PMMA material were obtained experimentally using sharp notched ENDB specimens.
Fig. 7. Synthesis based on SED of the PMMA results from combined tension and torsion tests
The proposed ENDB specimen in this research was capable of introducing the complete mixed mode I/III cases simply by changing the crack inclination angle relative to three-point bend rollers.
Although the fracture load of ENDB specimens made of PMMA was increased by moving from mode I to mode III loading condition, the average mode I fracture toughness of tested specimen was about 30% greater than K IIIc value.
Mixed mode I/III PMMA notch fracture toughness data were predicted successfully using the SED fracture criterion.
The SED approach has been found very suitable to summarize all the data in a narrow scatterband.
Acknowledgments
The authors would like to thanks Mr. Farhad Zivari (Machining workshop, School of Industrial Engineering at Iran University of Science and Technology) for his helps in manufacturing the test specimens.
References
1. Erdogan F., Sih G. C. On the crack extension in plates under plane loading and transverse shear // J. Basic. Eng. - 1963. - V. 85. - No. 4. -P. 519-525.
2. Maccagno T.M., Knott J.F. The fracture behaviour of PMMA in mixed modes I and II // Eng. Fract. Mech. - 1989. - V. 34. - No. 1. - P. 6586.
3. Aliha M.R.M., Ayatollahi M.R. Geometry effects on fracture behaviour
of polymethyl methacrylate // Mater. Sci. Eng. A. - 2010. - V. 527. -No. 3. - P. 526-530.
4. Awaji H., Kato T., Honda S., Nishikawa T. Criterion for combined mode I-II brittle fracture // J. Ceramic Soc. Japan. - 1999. - V. 107. -No. 10. - P. 918-924.
5. Ayatollahi M.R., Aliha M.R.M., Hassani M.M. Mixed mode brittle fracture in PMMA—an experimental study using SCB specimens // Mater. Sci. Eng. A. - 2006. - V. 417. - No. 1. - P. 348-356.
6. Ayatollahi M.R., Aliha M.R.M., Saghafi H. An improved semi-circular bend specimen for investigating mixed mode brittle fracture // Eng. Fract. Mech. -2011. - V. 78. - No. 1. - P. 110-123.
7. Saghafi H., Zucchelli A., Minak G. Evaluating fracture behavior of brittle polymeric materials using an IASCB specimen // Polym. Test. -2013. - V. 32. - No. 1. - P. 133-140.
8. Aliha M.R.M., Bahmani A., Akhondi S. Mixed mode fracture toughness testing of PMMA with different three-point bend type specimens // Eur. J. Mech. A/Solid. - 2016. - V. 58. - P. 148-162.
9. Ayatollahi M.R., Aliha M.R.M. Analysis of a new specimen for mixed mode fracture tests on brittle materials // Eng. Fract. Mech. - 2009. -V. 76. - No. 11. - P. 1563-1573.
10. Richard H.A. Some Theoretical and Experimental Aspects of Mixed Mode Fractures // 6 Int. Conf. on Fracture. - New Delhi, India, 1984.
11. Gómez F.J., Elices M. Fracture of components with V-shaped notches // Eng. Fract. Mech. - 2003. - V. 70. - No.14. - P. 1913-1927.
12. GómezF.J., Elices M. A fracture criterion for sharp V-notched samples // Int. J. Fract. - 2003. - V. 123. - No. 3-4. - P. 163-175.
13. Gómez F.J., Elices M., Planas J. The cohesive crack concept: Application to PMMA at -60o C // Eng. Fract. Mech. - 2005. - V. 72. -No. 8.- P. 1268-1285.
14. Gómez F.J., Elices M.., Berto F., Lazzarin P. Local strain energy to assess the static failure of U-notches in plates under mixed mode loading // Int. J. Fract. - 2007. - V. 145. - No. 1. - P. 29-45.
15. Susmel L., Taylor D. The theory of critical distances to predict static strength of notched brittle components subjected to mixed-mode loading // Eng. Fract. Mech. - 2008. - V. 75. - No. 3. - P. 534-550.
16. Berto F., Cendon D.A., Lazzarin P., Elices M. Fracture behaviour of notched round bars made of PMMA subjected to torsion at -60 °C // Eng. Fract. Mech. - 2013. - V. 102. - P. 271-287.
17. Ahmadi-Moghadam B., Taheri F. An effective means for evaluating mixed-mode I/III stress intensity factors using single-edge notch beam specimen // J. Strain Anal. Eng. Des. - 2013. - V 48. - No. 4. -P. 245-257.
18. Liu S., Chao Y.J., Zhu X. Tensile-shear transition in mixed mode I/III fracture // Int. J. Solids Struct. - 2004. - V. 41. - No. 22. - P. 61476172.
19. Ayatollahi M.R., Saboori B. A new fixture for fracture tests under mixed mode I/III loading // Eur. J. Mech. A/Solids. - 2015. - V. 51. -P. 67-76.
20. RichardH.A., Schramm B., Schirmeisen N.H. Cracks on mixed mode loading—Theories, experiments, simulations // Int. J. Fatigue. -2014.- V. 62. - P. 93-103.
21. Lazzarin P., Berto F. Some expressions for the strain energy in a finite volume surrounding the root of blunt V-notches // Int. J. Fract. -2005. - V. 135. - No. 1-4. - P. 161-185.
22. Lazzarin P., Zambardi R. A finite-volume-energy based approach to predict the static and fatigue behavior of components with sharp V-shaped notches // Int. J. Fract. - 2001. - V. 112. - No. 3. - P. 275-298.
23. Aliha M.R.M., Bahmani A., Akhondi S. Numerical analysis of a new mixed mode I/III fracture test specimen // Eng. Fract. Mech. - 2015. -V. 134. - P. 95-110.
24. Aliha M.R.M., Bahmani A., Akhondi S. Determination of mode III fracture toughness for different materials using a new designed test configuration // Mater. Des. - 2015. - V. 86. - P. 863-871.
25. Williams M.L. Stress singularities resulting from various boundary conditions // J. Appl. Mech. - 1952. - V. 19. - No. 4. - P. 526-528.
26. Zappalorto M., Lazzarin P., Berto F. Elastic notch stress intensity factors for sharply V-notched rounded bars under torsion // Eng. Fract. Mech. - 2009. - V. 76. - No. 3. - P. 439-453.
27. Zappalorto M., Lazzarin P., Filippi S. Stress field equations for U-and blunt V-shaped notches in axisymmetric shafts under torsion // Int. J. Fract. - 2010. - V. 164. - No. 2. - P. 253-269.
Поступила в редакцию 29.02.2016 г.
Сведения об авторах
Mohammad Reza Mohammad Aliha, Dr. PhD, Assist. Prof., IUST, Iran, m_aliha@yahoo.com, mrm_aliha@iustac.ir
Filippo Berto, Prof. PhD, Prof., University of Padova, Italy, berto@gest.unipd.it, Filippo.berto@unipd.it
Aram Bahmani, PhD Stud., IUST, Iran, bahmaniaram1@gmail.com
Shima Akhondi, PhD Stud., IUST, Iran, akhondi1994@gmail.com
Afrooz Barnoush, Prof. PhD, Prof., NTNU, Norway, afrooz.barnoush@ntnu.no
28. Zappalorto M., Lazzarin P., Yates J.R. Elastic stress distributions for hyperbolic and parabolic notches in round shafts under torsion and uniform antiplane shear loadings // Int. J. Solids Struct. - 2008. -V. 45. - No. 18. - P. 4879-4901.
29. Qian J., Hasebe N. Property of eigenvalues and eigenfunctions for an interface V-notch in antiplane elasticity // Eng. Fract. Mech. -
1997. - V. 56. - No. 6. - P. 729-734.
30. Berto F., Lazzarin P., Kotousov A., Harding S. Out-of-plane singular stress fields in V-notched plates and welded lap joints induced by inplane shear load conditions // Fatigue Fract. Eng. Mater. Struct. -2011. - V. 34. - No. 4. - P. 291-304.
31. Berto F., Lazzarin P. Fatigue strength of structural components under multi-axial loading in terms of local energy density averaged on a control volume // Int. J. Fatigue. - 2011. - V 33. - No. 8. - P. 10551065.
32. Gough H.J., Pollard H.V. Properties of some materials for cast crankshafts, with special reference to combined stresses // Proc. Inst. Auto. Eng. - 1963. - V. 31. - No. 1. - P. 821-893.
33. Lazzarin P., Sonsino C.M., Zambardi R. A notch stress intensity approach to assess the multiaxial fatigue strength of welded tube-to-flange joints subjected to combined loadings // Fatigue Fract. Eng. Mater. Struct. - 2004. - V. 27. - No. 2. - P. 127-140.
34. Lazzarin P., Berto F., Elices M., Gomez J. Brittle failures from U-and V-notches in mode I and mixed, I + II, mode: A synthesis based on the strain energy density averaged on finite-size volumes // Fatigue Fract. Eng. Mater. Struct. - 2009. - V. 32. - No. 8. - P. 671-684.
35. Berto F., Lazzarin P. A review of the volume-based strain energy density approach applied to V-notches and welded structures // Theor. Appl. Fract. Mech. - 2009. - V. 52. - No. 3. - P. 183-194.
36. Berto F., Lazzarin P. Recent developments in brittle and quasi-brittle failure assessment of engineering materials by means of local approaches // Mater. Sci. Eng. R. Rep. - 2014. - V. 75. - P. 1-48.
37. Berto F., Campagnolo A., Ayatollahi M.R. Brittle fracture of rounded V-notches in isostatic graphite under static multiaxial loading // Phys. Mesomech. - 2015. - V. 18. - No. 4. - P. 283-297.