CC BY
126
УДК 539.42
Compressive brittle fracture prediction in blunt V-notched PMMA specimens by means of the strain energy density approach
A.R. Torabi1, M.R. Ayatollahi2, M. Colussi3
1 Fracture Research Laboratory, Faculty of New Sciences and Technologies, University of Tehran, Tehran, 14395-1561, Iran 2 Fatigue and Fracture Research Laboratory, Center of Excellence in Experimental Solid Mechanics and Dynamics, School of Mechanical Engineering, Iran University of Science and Technology, Tehran, 16846-13114, Iran 3 Department of Management and Engineering, University of Padova, Vicenza, 36100, Italy
The paper aims to examine the suitability of the strain energy density criterion in predicting the fracture behavior of blunt V-notched specimens under compression load. Recent studies on local stress fields have shown that the strain energy density, averaged over a specific control volume which embraces the notch round border, could be a robust parameter in the brittle fracture assessment of several materials. A set of experimental results recently published in the literature on compressive brittle fracture of V-notched specimens made out of polymethyl methacrylate has here been considered. Finite element analyses have been performed on plane strain condition and experimental data have been summarized by means of the SED criterion. It has been shown that the proposed criterion permits a satisfactory evaluation of the fracture load of polymethyl methacrylate specimens weakened by notches having different opening angles and radii.
Keywords: brittle fracture, compression, V-notch, strain energy density, polymethyl methacrylate
DOI 10.24411/1683-805X-2018-11002
Оценка хрупкого разрушения образцов ПММА с тупым V-образным надрезом при сжатии на основе критерия плотности энергии деформации
A.R. Torabi1, M.R. Ayatollahi2, M. Colussi3
1 Тегеранский университет, Тегеран, 14395-1561, Иран 2 Иранский университет науки и технологии, Тегеран, 16846-13114, Иран 3 Падуанский университет, Виченца, 36100, Италия
В статье рассмотрено применение критерия плотности энергии деформации для оценки характера разрушения образцов с тупым V-образным надрезом в условиях сжимающей нагрузки. Согласно недавним исследованиям локальных полей напряжений, плотность энергии деформации, усредненная по объему материала вокруг закругленного конца надреза, может служить надежным параметром для оценки хрупкого разрушения ряда материалов. В настоящей работе обсуждаются опубликованные результаты экспериментов по хрупкому разрушению образцов полиметилметакрилата с V-образным надрезом при сжатии. Приведены результаты моделирования методом конечных элементов задачи для случая плоской деформации. Экспериментальные данные обобщены с помощью критерия плотности энергии деформации. Показано, что данный критерий позволяет удовлетворительно оценить разрушающую нагрузку для образцов полиметилметакрилата с надрезами, имеющими разные углы раскрытия и радиусы.
Ключевые слова: хрупкое разрушение, сжатие, V-образный надрез, плотность энергии деформации, полиметилметакрилат
1. Introduction
The two main groups of stress concentrators in components and structures are cracks and notches. A literature survey indicates that brittle fracture investigations in components containing stress raisers under compression are limited to notched members. This is probably because the clo-
© Torabi A.R., Ayatollahi M.R., Colussi M., 2018
sure phenomenon takes place in cracked members under compression and, the faces coming in contact with each other, little chance of crack propagation is allowed. Conversely, in a notched member, particularly in V- and U-shaped notches, the faces may or may not come in contact with each other, depending on the gap value and the maxi-
mum face displacement at the onset of failure, and the compressive damage may nucleate and lead to the final fracture.
Studies on notched components under tensile and shear loading conditions have been widely performed during the past two decades, see a comprehensive list of references in Refs. [1, 2]. However, it is of crusial importance to develop appropriate failure criteria for structures containing various shaped notches, such as V-, U- and O-notches, under compressive loading conditions to bridge the current gap on the topic and prevent the onset of sudden fracture.
The first use of the notch fracture mechanics principles in compressive brittle fracture prediction of notched components has probably been made in [3]. They have recently investigated the fracture behavior of graphite specimens weakened by double V-notches with end holes (VO-notches) both experimentally and theoretically by means of the average strain energy density (SED) criterion. Torabi and Ayatollahi [4] have successfully predicted the test results reported in [3] by means of the two brittle fracture criteria, known as the point-stress (PS) criterion and mean-stress (MS) criterion. The most recent works on compressive brittle fracture of notched members are those published by Ayatollahi et al. and Torabi et al. [5, 6] in which the PS and MS fracture criteria have been successfully employed for predicting the fracture loads of two new V-notched test specimens, called the V-notch stepped cottage (VSC) and the flattened V-notched semidisk (FVSD) specimens, made of polymethylmethacrylate (PMMA), for various notch angles and different notch radii.
Since the extension of the brittle fracture criteria under tension and shear loading conditions to compressive loading conditions has taken place in the most recent years, their validity should strongly be verified by means of various experimental results on different brittle materials and test specimens. The attempt made here is to examine the suitability of the average strain energy density (SED) criterion in predicting the fracture load of PMMA specimens. The criterion states that brittle fracture occurs when the strain energy density, averaged in a specific control volume, which includes a crack or notch tip, reaches a critical value dependent on the material. Good agreement has been found between numerical predictions and the experimental results, demonstrating that the proposed criterion works well not only on compressive brittle fracture of VO-notched graphite specimens, but also on that of blunt V-notched PMMA specimens.
2. Experimental results from the literature
A set of new experimental results have been recently published by Ayatollahi et al. [5] dealing with brittle fracture in blunt V-notches under compression. The details about specimens and test procedures have been reported in the following subsections.
Table 1
Some of the mechanical properties of PMMA at room temperature
Material property Value
Elastic modulus in compression, MPa 2018
Poisson's ratio 0.38
Ultimate tensile strength, MPa 70.5
Plane-strain fracture toughness, MPa m05 1.96
Compressive yield strength, MPa 96
Ultimate compressive strength, MPa 98
2.1. Material
Polymethylmethacrylate (PMMA) has been used for conducting fracture experiments at room temperature on blunt V-shaped notches under compression. The essential mechanical properties of the adopted PMMA at room temperature have been taken from the experiments performed earlier by Ayatollahi et al. [5], according to the standard ASTM D695-10, and are shown in Table 1. In Fig. 1 the true stress-strain curve of the PMMA sample tested under compression is depicted.
2.2. Specimen
The adopted specimen, named V-notch stepped cottage (VSC) specimen, is schematically shown in Fig. 2. It has firstly been proposed and utilized in Ref. [5] for conducting compressive fracture experiments on blunt V-notches. In Fig. 2, the parameters 2a, p, a, b, c, d, e, L, t and P are respectively the notch opening angle, notch radius, notch length, flattened length, step height, step width, specimen height, specimen length, specimen thickness and the com-pressive load.
Due to the bending effects, tensile stresses are unavoidably generated at some points far enough from the notch tip. By making a step on the top and bottom of the specimen and applying the compressive load to the remaining
True strain
Fig. 1. A sample true stress-strain curve of the PMMA under compression according to the results of [5]
Fig. 2. A scheme of the tested VSC specimen as in Ref. [5]
edge, the notch surroundings experience pure compressive deformations which could be significantly higher compared to tensile ones. The particular shape of the proposed speci-
men is then a consequence of the need to make sure to have compressive fracture from the notch border before tensile fracture from other locations of the specimen. Because the compressive strength of brittle materials is usually two to three times greater than the tensile one, it was necessary to adopt VSC specimens having dimensions so that the level of the compressive stresses at the notch tip was several times greater than the maximum tensile stress in the whole specimen. The parameters above introduced have been chosen as 2a = 30°, 60° and 90°; p = 0.5, 1.0 and 2.0 mm; a = 10 and 25 mm; b = 30 mm; c = 20 mm; d = 5 mm; e = 85 mm; L =110 mm and t = 6 mm. Note that these dimensions have also been chosen in order to make sure that buckling would not take place during the experiments. Tested specimens have been fabricated by laser cutting from a PMMA plate of
Table 2
The experimental fracture loads and computed SED values of the VSC specimens
a = 10 mm a = 25 mm
2a p, mm Pf, kN W, MJ/m3 2a p, mm Pf, kN W, MJ/m3
30° 0.5 -13.10 11.67 30° 0.5 -8.36 18.79
-13.31 12.05 -7.53 15.24
-17.40 20.59 -8.62 19.98
1.0 -11.42 9.13 1.0 -6.29 10.97
-10.79 8.15 -7.12 14.05
-13.27 12.33 -5.74 9.13
2.0 -12.70 11.65 2.0 -8.96 22.59
-12.74 11.73 -8.86 22.09
-10.93 8.63 -6.91 13.44
60° 0.5 -11.37 9.55 60° 0.5 -7.21 14.76
-12.24 11.07 -7.72 16.92
-12.29 11.16 -7.60 16.40
1.0 -9.84 7.22 1.0 -5.99 10.27
-9.59 6.86 -6.12 10.72
-10.15 7.68 -6.67 12.73
2.0 -11.98 10.73 2.0 -7.28 15.00
-11.81 10.43 -6.58 12.25
-17.27 22.30 -6.18 10.81
90° 0.5 -12.59 12.12 90° 0.5 -8.16 18.42
-12.76 12.45 -8.34 19.24
-10.42 8.30 -8.47 19.85
1.0 -8.91 6.03 1.0 -4.95 6.72
-9.70 7.15 -5.37 7.91
-8.68 5.72 -5.28 7.65
2.0 -10.46 8.13 2.0 -5.92 9.34
-10.20 7.73 -6.56 11.47
-8.56 5.45 -6.63 11.71
6 mm thick. The cut surfaces were accurately polished by means of appropriate abrasive papers to remove possible defects due to the cutting process. The strain rate in the tests was set to 0.000075 1/s providing quasi-static monotonic loading conditions. Three specimens have been tested for each geometry in order to check the repeatability of the test. Totally, 54 fracture experiments have been carried out.
2.3. Test results
Table 2 summarizes the experimentally recorded fracture loads. It is possible to point out that all the fracture loads are less than 18 kN and most of them are less than 13 kN. No evidence of buckling has been observed.
A sample measured load-displacement curve is shown in Fig. 3. A linear trend has been found up to a load peak, after which a sudden fall to zero occurs, confirming the predominantly linear elastic behavior of the specimens up to the final failure. For this reason, the experimental results can potentially be predicted by means of brittle fracture criteria based on linear elastic notch mechanics principles, e.g. the averaged SED criterion.
In the forthcoming sections, the SED averaged over a specified control volume which embraces the notch border is presented and adopted to predict the measured fracture loads summarized in Table 2.
3. Averaged strain energy density approach
The averaged strain energy density criterion (SED), according to [7], states that brittle failure occurs when the mean value of the strain energy density over a given control volume W, reaches the critical value Wc. In agreement with [8], named Wc the ultimate compression strength under elastic stress field conditions and E the Young's modulus of the material, the critical value of the total strain energy can be determined by the following:
(1)
W =-L a 2.
2 E
This critical value varies from material to material but it does not depend on the notch geometry and sharpness.
Under plane strain condition and tension loading, the following expression has been proposed by Lazzarin and Zambardi [7] to evaluate the control volume dimension R0 :
Ro =
(1 + v)(5 -8 v)
4n
(2)
where KIc is the material fracture toughness, v is the Poisson's ratio and at is the ultimate tensile strength of an unnotched specimen. Under compressive load, Eq. (2) cannot be used because there is no definition for fracture toughness. An empirical approach can then be a good alternative for determining R0 for notched components under compression, as described in the following section.
4. Finite element model
Linear static finite element analyses were conducted in order to numerically evaluate the averaged SED. Analyses were performed in 2D plane strain condition, by means of ANSYS® software and adopting quadratic finite elements. The mesh insensitivity of the SED approach has already been proved by Lazzarin et al. [9] and it is a consequence of the finite element method, in which the elastic strain energy is computed directly from the nodal displacements, without involving stresses and strains. Figure 4, a shows the schematic representation of the finite element model and a sample of the coarse mesh adopted.
The external load was introduced as uniformly distributed along the top flat line of the specimen and the nodes
-3 -2 -1 0
Displacement, mm
Fig. 3. A sample load-displacement curve of the VSC specimen according to the results of [5]
-250.334 -190.804 -131.274 -71.7438 -12.2138
-220.569 -161.039 -101.509 - 41.9788 17.5513
Fig. 4. Finite element model and stress contour for the case of 2a = 30°, p = 0.5 mm and a = 25 mm. Schematic representation of the finite element model (a), first principal stress contour lines (b)
0.133936 0.538935 0.943933 1.34893 1.75393 0.336435 0.741434 1.14643 1.55143 1.95643
Fig. 5. Shape of the control volume under mode I loading and SED contour lines for the case of 2a = 30°, p = 0.5 mm and a = = 25 mm. Control volume for blunt V-notches (a), sample of SED contour lines (b)
belonging to the same line were constrained so that they could move only along the loading direction. The nodes belonging to the corresponding bottom line of the specimen were completely fixed. A contour plot of the main principal stress component g1 is shown in Fig. 4, b for the case 2a = 30°, p = 0.5 mm and a = 25 mm. According to the SED approach, for rounded V-notches the critical volume is represented by a portion of circular sector. The shape of the volume is shown in Fig. 5, a while the contour lines of the SED are depicted in Fig. 5, b for the same case presented in Fig. 4.
On the notch side the control volume boundary is represented by the circular notch edge with radius p, on the other side the boundary is a circle arc with radius R0 + r0. The center of the notch radius is geometrically given while the center of the radius R0 + r0 coincides with the center of the polar coordinate system in the analytical solution for the blunt V-notch problem in which r0 = p(q -1)/q, with q = (2n - 2a)/n, according to [10].
5. Results
By considering five different values of R0, ranging from 0.5 to 1.5 mm, and by plotting the computed averaged SED as a function of R0, it has been found that assuming R0 = = 1 mm the specimens having p = 0.5 mm, 2a = 30° and p = 2 mm, 2a = 90° experience the same averaged SED value at failure. This means that at the critical load the specimens are characterized by the same SED value averaged in a control volume having radius equal to 1 mm, which is
Fig. 6. SED based summary from V-notched PMMA specimens under compression
independent of the notch shape and sharpness. This value also corresponds to the critical value that can be obtained by Eq. (1) using as compressive strength Gc = 230 MPa, according to [28], and E = 2018 MPa. For the herein considered PMMA the critical SED has been found equal to 13.1 MJ/m3. Figure 6 shows a summary of the experimental data in terms of the square root value of the SED averaged on the control volume of radius R0 =1 mm, normalised with respect to the constant value of the critical SED (13.1 MJ/m3). Indeed, this parameter is proportional to the fracture load. The aim was to investigate the range of accuracy of the SED-based fracture assessment for the PMMA specimens under compression. It has been found that about all experimental data fit in a narrow scatter band, whose limits are drawn here with an engineering judgement from 0.75 to 1.25. However, it should be pointed out that the majority of the results falls inside a band ranging from 0.8 to 1.1. It can be concluded that the scatter of the data is quite limited and substantially independent of the notch geometry and, in the present author opinion, the averaged SED criterion appears suitable for the fracture strength assessment of V-notched PMMA specimens under compression.
6. Conclusion
Brittle fracture in V-notched PMMA specimens was investigated under compression loading both experimentally and numerically. Fracture tests conducted on prismatic specimens and reported in the recent literature have been reanalyzed in this paper by means of the averaged strain energy density criterion. Being the experimental results in good agreement with the numerical assessments, it has been shown that the proposed method is suitable for predicting PMMA brittle failure under compressive loading condition and in the presence of notches having a different opening angle and sharpness.
References
1. Berto F., Lazzarin P. Recent developments in brittle and quasi-brittle failure assessment of engineering materials by means of local approaches // Mater. Sci. Eng. - 2014. - V. 75. - P. 1-48.
2. Cendon D.A., Torabi A.R., Elices M. Fracture assessment of graphite V-notched and U-notched specimens by using the cohesive crack model // Fat. Fract. Eng. Mater. Struct. - 2015. - V. 38. - No. 5. - P. 563573.
3. Berto F., Lazzarin P., Ayatollahi M.R. Brittle fracture of sharp and blunt V-notches in isostatic graphite under pure compression loading // Carbon. - 2013. - V. 63. - P. 101-116.
4. Torabi A.R., Ayatollahi M.R. Compressive brittle fracture in V-notches with end holes // Eur. J. Mech. A. Solids. - 2014. - V. 45. - P. 32-40.
5. Ayatollahi M.R., Torabi A.R., Firoozabadi M. Theoretical and experimental investigation of brittle fracture in V-notched PMMA specimens under compressive loading // Eng. Fract. Mech. - 2015. -V. 135.- P. 187-205.
6. Torabi A.R., Firoozabadi M., Ayatollahi R. Brittle fracture analysis of blunt V-notches under compression // Int. J. Solid. Struct. - 2015. -V. 67-68. - P. 219.
7. 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. - P. 275-298.
8. Beltrami E. Sulle condizioni di resistenza dei corpi elastici // Rendiconti Del Regio Istituto Lombardo. - 1885. - V. XVIII. - P. 704-714.
9. Lazzarin P., Berto F., Zappalorto M. Rapid calculations of notch stress intensity factors based on averaged strain energy density from coarse meshes: Theoretical bases and applications // Int. J. Fatig. - 2010. -V. 32. - P. 1559-1567.
10. 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. - P. 161-185.
Поступила в редакцию 23.01.2017 г.
Сведения об авторах
Ali Reza Torabi, PhD, University of Tehran, Iran, a_torabi@utac.ir
Majid R. Ayatollahi, PhD, Prof., Director, Iran University of Science and Technology, m.ayat@iust.ac.ir Marco Colussi, PhD student, University of Padova, Italy, colussim@hotmail.com
