Microsoft Word - numero_59_art_27_3359.docx S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 405 The critical influence of some “tiny” geometrical details on the stress field in a Brazilian Disc with a central notch of finite width and length Stavros K. Kourkoulis, Christos F. Markides, Ermioni D. Pasiou National Technical University of Athens, School of Applied Mathematical and Physical Sciences, Department of Mechanics, Laboratory for Testing and Materials, Theocaris Building, Zografou Campus, 157 73, Athens, Greece stakkour@central.ntua.gr, http://orcid.org/0000-0003-3246-9308 markidih@mail.ntua.gr, http://orcid.org/0000-0001-6547-3616 epasiou@teemail.gr, http://orcid.org/0000-0003-1580-3415 Andronikos Loukidis, Dimos Triantis University of West Attica, Department of Electrical & Electronics Engineering, Electronic Devices and Materials Laboratory, 250 Thivon Avenue, 122 44, Athens, Greece a.loukidis@uniwa.gr, http://orcid.org/0000-0002-5764-5366 triantis@uniwa.gr, http://orcid.org/0000-0003-4219-8687 ABSTRACT. The role of some geometrical characteristics of the notches ma- chined in circular discs, in order to determine the mode-I fracture toughness of brittle materials, is discussed. The study is implemented both analytically and numerically. For the analytic study advantage is taken of a recently intro- duced solution for the stress- and displacement-fields developed in a finite disc with a central notch of finite width and length and rounded corners. The variation of the stresses along strategic loci and the deformation of the peri- meter of the notch obtained analytically are used for the calibration/validation of a flexible numerical model, which is then used for a parametric investiga- tion of the role of geometrical features of the notched disc (thickness of the disc, length and width of the notch, radius of the rounded corners of the notch). It is concluded that the role of the width of the notch is of critical im- portance. Both the analytic and the numerical studies indicate definitely that ignoring the accurate geometric shape of the notch leads to erroneous results concerning the actual stress field around the crown of the notch. Therefore, it is possible that misleading values of the fracture toughness of the material of the disc may be obtained. KEYWORDS. Brazilian disc; Central notch; Analytical solution; Numerical ana- lysis; Finite elements. Citation: Kourkoulis, S.K., Markides, Ch.F., Pasiou, E.D., Loukidis, A., Triantis, D., The critical influence of some “tiny” geometrical details on the stress field in a Brazilian Disc with a central notch of finite width and length, Frattura ed Integrità Strutturale, 59 (2022) 405- 422. Received: 11.11.2021 Accepted: 21.11.2021 Published: 01.01.2022 Copyright: © 2022 This is an open access ar- ticle under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and re- production in any medium, provided the orig- inal author and source are credited. https://youtu.be/p8BWx0005AU S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 406 INTRODUCTION ode-I fracture toughness, KIC, is a critical mechanical parameter for structural engineering applications because it quantifies the resistance of materials against the initiation (and propagation) of pre-existing cracks due to ten- sile loads acting normally to the crack axis. Both ASTM and ISRM have long ago standardized the procedure for the laboratory determination of KIC [1,2]. Especially for brittle, rock-like materials, three standardized tests are usually adopted: The “Short Rod”, the “Chevron Bend” and the “Cracked Chevron Notched Brazilian Disc (CCNBD)” ones. The latter is, perhaps, the one most widely used worldwide. However, and in spite of its wide application, some issues con- cerning the validity of its outcomes are still under study [3,4] and some critical questions are still open [5]. Among them special attention is paid to: (i) the reasonability of the formulae available for the determination of the respective Stress Intensity Factor (SIF), (ii) the exact shape of the cracks machined and (iii) the actual boundary conditions prevailing along the disc-jaw contact arc. Currently, the formulae adopted for the determination of the SIF in the CCNBD test are obtained from the respective “Cracked Straight Through Brazilian Disc” (CSTBD) configuration, i.e., considering a disc with a central “mathematical” crack (i.e., a linear discontinuity with zero distance between its lips and singular tips) loaded by a pair of diametral forces. It is obvious, however, that in practical applications the cracks machined are rectangular notches rather than mathematical cracks [6] and the disc is loaded by a complicated combination of radial (and shear) stresses rather than by a pair of point forces [7,8]. Things become even more complicated taking into account that for obvious practical restrictions the “cracks” machined cannot be “short” enough (with respect to the radius of the disc) and, therefore, the common assumption that the boundary conditions do not influence the stress field in the immediate vicinity of the crack becomes questionable. In this direction, an attempt is described in this study to address some of the aforementioned issues. The main target is to quantify the influence of some “tiny” geometrical details on the stress field (and, therefore on the respective SIF or, equiva- lently, on the value of KIC) in a Brazilian Disc with a central notch of finite width and length and rounded corners. In addition, in the analytic approach, the loading scheme comprises a parabolic distribution of radial stresses, acting along two finite arcs. The specific loading scheme is proven to approach extremely closely the stress field developed during the compression of the disc between the curved jaws suggested by ISRM [2] (Fig.1a). The geometrical details considered in this study include: (i) the length of the notch with respect to the radius of the disc, (ii) the width of the notch, (iii) the radius of the corners of the notch, and, finally, (iv) the thickness of the disc. As a first step the problem is confronted analytically by means of the complex-potentials technique [9]. Analytical (but lengthy) full-field formulae, in series form, are obtained for the components of the stress- and displacement-fields. Taking into account the difficulties to handle these expressions for practical purposes and exhaustive parametric investigations, a numerical model is developed, calibrated and validated with the aid of the analytic solution. The numerical model is then used for the parametric study of the role of the geometrical characteristics (length, width and corner-radius) of the notch. Moreover, the numerical model is used to enlighten the influence of the third dimension (i.e., that along the disc’s thick- ness), a parameter that is not taken into account by the analytic solution, which is based on two-dimensional elasticity. The results of the study indicate clearly that ignoring the exact geometric shape of the notch may lead to erroneous results concerning the stress field in the immediate vicinity of the crowns of the notch and therefore to unacceptable results concerning the actual fracture toughness of the disc’s material. (a) (b) Figure 1: (a) Schematic representation of the notched disc squeezed between the curved jaws suggested by ISRM [2]; (b) The mapping function provided by Savin [11], which was used for the analytic solution of the problem. Pframe Pframe RO ISRM’s jaw L w α2 α1=eikπ α3 α3 θ ζ=ρeiθ η ξ Ο A2 A1 A3 A3 y xΟ  z=reiϕ ∞ z=ω(ζ) L/2 w/2 γ RO 1 M S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 407 ANALYTICAL SOLUTION rom the theoretical point of view, the problem of the finite circular disc with a rectangular notch (of finite length and width) with rounded corners, squeezed between the curved jaws suggested by the respective ISRM standard (Fig.1a), is here confronted as a first fundamental problem of the plane theory of linear elasticity. In this context, and assuming, in addition, that the material of the disc is homogeneous and behaves isotropically, Muskhelishvili’s com- plex potentials technique is adopted [9]. The general solution of the problem was obtained, relatively recently, by Markides and Kourkoulis [4,10] by taking advantage of the following mapping function, provided by Savin [11] (Fig.1b):                          1 2 1 1 1 2 2 2 4 4 2 2 1 2 3 4 5 5 3 3 6 6 4 4 2 2 5 6 2 1 2 , 1 1 2 5 4 2 ; ; ; ; 2 24 80 896 7 5 2 21 14 5 4 ; 2304 11264 ( , 0 0, , ,i k ik ik c b z R R c c c c c c c e k for w e e                                                                                          3 4, )ik ike e      (1) The above function maps the region outside the area A1A2A3A4 (see Fig.1b) to the region outside γ in the mathematical plane. In other words, one introduces Eq.(1) into the already available solution [4, 10] of the circular ring of outer radius RO. Then, demanding that the rectangular notch A1A2A3A4 must be free from stresses, and assuming, also, that the notch has too little affection on the boundary of the disc, the solution of the problem is sought in the form:            ;* *o oz z z z z z         (2) where  o z and  o z are given by the following expressions:                     2 2 1 2 2 14 1 4 13 1 4 1 4 34 42 0 2 1 2 2 1 4 1 2 2 14 33 1 3 4 1 4 34 42 0 2 3 4 1 4 1 4 3 4 1 3 3 4 1 4 3 4 3 n nn nn nc n n o n o n n nnn nc n B BP B BB z b z z B z z z z z n n n n z B B BP B BB z z b z z z z z n n n                                                            4 1 1 4 1 n n z n            (3) providing the solution of the circular ring, as it is analytically described by Kourkoulis et al. [12]. In the above expressions φ * (z) and ψ * (z) are analytic functions perturbing the ring’s solution due to the introduction of the notch. Alternatively, instead of the ring the intact disc could have been used as the solution basis. In this direction, it is assumed that the inner radius of the ring tends to zero. Combining Eqs.(1)-(3) and introducing them in the boundary condition for zero stresses on the rectangular notch:           0 s s s s s          (4) (where s is the point ζ on γ) one obtains (after some relatively lengthy algebra) the quantities φ * (z) and ψ * (z) solving the problem [10]. The configuration of this alternative problem is shown in Fig.2, where, in accordance with the solution of the circular ring, the disc with the notch is to be subjected to a parabolic pressure along two symmetric arcs of its peri- phery. In the same figure the introduction of the curvilinear coordinate system ρ=const., θ=const., at any point of the disc with the notch is shown, as it is dictated by the mapping function described in Eq.(1). F S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 408 Figure 2: The notched disc under a parabolic distribution of radial stresses along two symmetric finite arcs of its perimeter (left) and the curvilinear coordinate system used in the analytic solution (right). Then, stresses and displacements at any point of the disc are obtained using the well-known expressions [9]:                                     2 2 , 4 2 ( ) , , ie i u iv                                                                  (5) Using these formulae, the deformed configuration of the notch, as well as the stress-field components are determined. The deformed shape of the notch and the variation of the stresses along characteristic loci of the notched disc will be considered in juxtaposition to the respective numerical results that are obtained according to the analysis of the following section. As a typical example, highlighting the capabilities of the above-described analytic approach, the transverse stress σθ (nor- malized over the amplitude, Pc, of the parabolic distribution of radial stresses exerted on the disc) is plotted in Fig.3, along a part of the perimeter of the notch, which corresponds to the [0o, 60o] interval of the angle θ in the ζ-plane and the θ=const. hyperbolas in the z-plane (Fig.3a). Two cases are studied, one corresponding to a relatively wide notch, i.e., with w=5 mm (Fig.3b) and one corresponding to a relatively narrow one, i.e., with w=0.5 mm (Fig.3c). The crucial role of the parameter w is obvious by simply comparing Fig.3b and Fig.3c. Indeed, for the wider notch the normalized values of the σθ stress component vary in the (-0.80Pc, 0.85Pc) interval. On the other hand, for the narrower notch the respective values of the σθ stress component vary in the (-0.50Pc, 5.20Pc) interval. The maximum value of σθ in the case of the narrow notch is more than six times higher compared to that of the wider notch. It is obvious that ignoring the role of w (by assuming, for ex- ample, that the notch resembles a mathematical crack, for which w tends to zero, an assumption that is impossible to be realized for practical reasons) definitely undermines the validity of the outcomes of the standardized techniques employed for the determination of the fracture toughness, KIC. NUMERICAL STUDY he formulae provided by the analytic solution have the advantage of being full-field (i.e., they provide the stress- and displacement-field components at any point of the notched disc) and of closed form, however they are quite lengthy. Therefore, employing them in exhaustive parametric analyses is rather tedious. Moreover, the fact that before applying these formulae, convergence issues of the series expansions have to be carefully considered, renders their T S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 409 (a) (b) (c) Figure 3: (a) The locus along which the stresses are plotted (red line); (b, c) The variation of the transverse stress σθ (normalized over the amplitude, Pc, of the parabolic distribution of radial stresses exerted on the disc) along the part of the perimeter of the notch shown in Fig.3a, for w=5 mm (b) and w=0.5 mm (c). practical use inflexible and time consuming. In this direction, it was thought that one could take advantage of the analytic formulae in order to validate a numerical model, which, in turn, could be used for thorough parametric investigation of the geometric factors. In this context, a three-dimensional numerical model is designed using the finite elements method and employing the commercially available software ANSYS-19, as it is described in next sections. Reference numerical model A disc of diameter equal to 100 mm and thickness equal to 10 mm, with a central notch of 50 mm length and 5 mm width, was chosen as the reference model. The radius of curvature at the notch’s corners was set equal to 0.33 mm. The disc was assumed to be squeezed between the curved loading platens suggested by the ISRM standard. The whole complex consisting of the disc and the two loading platens was modelled (Fig.4), in an attempt to approach the actual loading scheme developed along the disc-loading platen interface (without any additional assumption for this issue). The disc was considered to be made of plexiglass (modulus of elasticity and Poisson’s ratio equal to 3.2 GPa and 0.36, respectively) while the mechanical properties of steel (modulus of elasticity and Poisson’s ratio equal to 210 GPa and 0.30, respectively) were assigned to the loading platens. The coefficient of friction between the plates and the disc was set equal to 0.01, in order to approach the assumptions adopted in the analytic solution, where the role of friction was ignored. The volumes of the model were meshed using the SOLID185 element while the TARGE170 and CONTA173 elements were used for the two material interfaces between the platens and the disc. The final meshing (Fig.5) was decided based on a thorough convergence analysis carried out as a first step of the numerical study. -2.0 0.0 2.0 4.0 6.0 0 20 40 60 T he σ θ st re ss ( no rm al iz ed o ve r P c) a lo ng pa rt o f th e no tc h’ s pe ri m et er ( 0o ≤ θ ≤6 0o ) Angle θ in the ζ-plane [deg], corresponding to the θ=constant hyperbola in the z-plane -1.0 -0.5 0.0 0.5 1.0 0 20 40 60 T he σ θ st re ss ( no rm al iz ed o ve r P c) a lo ng pa rt o f th e no tc h’ s pe ri m et er ( 0o ≤ θ ≤6 0o ) Angle θ in the ζ-plane [deg], corresponding to the θ=constant hyperbola in the z-plane Angle θ in the ζ-plane [deg], corresponding to the θ=constant parabola in the z-plane Angle θ in the ζ-plane [deg], corresponding to the θ=constant parabola in the z-plane T he n or m al iz ed σ θ st re ss a lo ng p ar t o f th e pe ri m et er o f th e no tc h w ith 0o ≤ θ ≤6 0o T he n or m al iz ed σ θ st re ss a lo ng p ar t o f th e pe ri m et er o f th e no tc h w ith 0o ≤ θ ≤6 0o S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 410 The nodes of the base of the lower platen were rigidly clamped while the nodes of the lateral surfaces of both platens were not permitted to move along the horizontal x-axis. Finally, a vertical downwards displacement equal to 2.275 mm was imposed on the nodes of the upper surface of the upper platen (Fig.6). The specific displacement value corresponds to the 20 kN load applied in the analytic solution used for the validation of the numerical model. Validating the numerical model taking advantage of the analytical solution Due to the double symmetry of the model, the results obtained by the numerical model and those provided by the analytic solution will be presented only for one quarter of the disc-loading platens complex. More specifically, in Fig.7a the initial shape of one quarter of the perimeter of the notch (black line) is drawn, in juxtaposition to its deformed shape, as it was obtained by the analytical (red line) and the numerical (blue line) solutions. The difference between them does not exceed 4% at any point of the locus considered, which is obviously well accepted. Even, this small difference may be explained by the fact that the series expressions of the analytic solution used for the calculation of the displacements are too complicated imposing limitations on the number of terms that can be used. Figure 4: The reference numerical model. Figure 5: The meshing of the reference model and the material interfaces. S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 411 Figure 6: The loading and the boundary conditions. Figure 7: (a) Undeformed and deformed geometries of one quarter of the perimeter of the notch; normal stresses developed along the vertical y-axis (b) and the horizontal x-axis (c), as they were obtained by both the analytical and the numerical approaches. -80 -60 -40 -20 0 20 2.5 12.0 21.5 31.0 40.5 50.0 σ n or m al [ M P a] x [mm] σxx-central σyy-central σxx-front σyy-front σxx-analytical σyy-analytical 25 30 35 40 45 50 -150 -75 0 75 150 y [m m ] σnormal [MPa] σxx-central σyy-central σxx-front σyy-front σxx-analytical σyy-analytical 0 10 20 30 0 1 2 3 y [m m ] x [mm] Initial geometry Analytical Numerical (a) (b) (c) S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 412 Concerning now the stress field components, it is here recalled that the numerical model deals with the problem in three dimensions (contrary to the analytic solution which is by default two dimensional). In this context, the normal stresses de- veloped both at the central cross section of the disc as well as these developed within the plane of its front face are here taken into account. The comparison between numerical and analytical data is seen in Fig.7b for the normal stresses which are developed along the vertical y-axis (i.e., the line of symmetry along the longitudinal axis of the notch) and in Fig.7c for the normal stresses developed along the horizontal x-axis (i.e., the line of symmetry normal to the longitudinal axis of the notch). It is clearly seen from these figures that non-ignorable differences between the data of the analytic and numerical solutions appear only as one approaches the immediate vicinity of the disc-loading platens contact arc. On the contrary, in the immediate vicinity of the crown of the notch (i.e., in the area of utmost importance for the purposes of the present study) the mutual agreement is excellent. It is thus concluded that the numerical model can be safely used to study either more complicated loading schemes or for exhaustive parametric analyses of the role of various geometrical characteristics. Parametric investigation After its validation, the reference model was properly modified to design six additional models, in order to quantify the influence of the geometrical characteristics of the notch on the stress field developed in the disc. The modified geometries of the discs are presented in Fig.8 (the loading platens are not shown in this figure, although the whole disc-loading platens complex is again simulated, as it was the case for the reference model). The radius of curvature of the corners of the notch, and the length and width of the notch were the characteristics considered and the values assigned to each one of them are shown in Table 1. Moreover, the role of the third dimension, namely that along the disc’s thickness was, also, studied. (b) (c) (d) Figure 8: The numerical models (only the discs are presented here) of the parametric analysis. (a) The reference disc, (b) enlarged views of one of the corners of the notch (parameter r); (c) the discs for the parameter L and (d) the discs for the parameter w. Parameter Model Radius of the corners of the notch, r [mm] Length of the notch, L [mm] Width of the notch, w [mm] Reference model 0.33 50 5 r=0.66 mm 0.66 50 5 r=1.00 mm 1.00 50 5 L=30 mm 0.33 30 5 L=75 mm 0.33 75 5 w=1 mm 0.33 50 1 w=3 mm 0.33 50 3 Table 1: The parameters considered during the numerical analysis and the values assigned to each model. r=0.66 mm L=30 mm w=3 mm r=1.00 mm L=75 mm w=1 mm Reference (a) S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 413 RESULTS The stress field for the reference model and the role of the third dimension he distribution of the equivalent stress developed in the whole complex is presented in Fig.9a. As it was expected, the most severely stressed areas are those in the immediate vicinity of the rounded corners of the notch. In ad- dition, it is seen that the equivalent stress attains lower values as one is directed from the face of the disc towards its central plane. Plotting the equivalent stress along one quarter of the notch, the previously mentioned comment is more clearly seen (Fig.9b). The first point of the plots (position at 0 mm, Point M) corresponds to the midpoint of the width of the notch and the last one (position at 28 mm, Point K) corresponds to the midpoint of the notch’s length (clockwise di- rection). For clarity, the locus considered was divided into three segments (embedded sketch in Fig.9(b2)), i.e., the horizontal segment or in other words the one along the x-axis (continuous line), the corner of the notch (dashed line) and the vertical one or in other words the segment which is parallel to the y-axis (dotted line). The same line code was used for the stresses, for which, in addition, blue colour indicates the stresses developed at the central cross section of the disc while the red one the stresses at its front face. It is seen from Fig.9b that the maximum equivalent stress at the front surface of the disc is about 15% larger than the respective value at the central cross section. It is, also, interesting to note that the equivalent stress at the front surface of the disc becomes significantly larger from the equivalent stress at the central cross section only at the immediate vicinity of the corner of the notch. (a) (b1) (b2) Figure 9: Reference model. (a) The overall distribution of the equivalent stress in [Pa] over the disc-loading platens complex; (b1) its variation along one quarter of the perimeter of the notch and (b2) a detailed view around the corner of the notch. 60 100 140 180 0 7 14 21 28 σ e qv [ M P a] Along the notch [mm] central front 60 100 140 180 0 2 4 6 σ e qv [ M P a] Along the notch [mm] central front Upper view M K T S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 414 Contrary to the equivalent stress, the maximum principal stress, σ1, seems to attain larger values as one is directed from the face of the disc towards its center, as it is seen from Fig.10a. A more conceivable interpretation of the differences of the maximum principal stresses, σ1 over the plane of the front face (red colour) and over that of the central cross section (black colour) is depicted in Fig.10b. In this figure one quarter of the notch was again isolated, due to the double symmetry of the model, and the stress developed at each point of the notch is drawn (for clarity reasons, the numerical values of σ1, attained at some strategic points of the locus are added in the figure). As it was expected, the maximum principal stress along the length of the notch (free boundary) equals zero. Based on this figure, it can be stated that, contrary to the equiva- lent stress, the maximum principal stress is significantly larger at the central cross section (black colour) along the major part of the perimeter of the notch. Only along the notch’s rounded corners the two plots tend to coincide, or in other words the role of the third dimension becomes negligible. (a) (b) Figure 10: Reference model. (a) The distribution of the maximum principal stress in [Pa] and (b) its variation along one quarter of the notch over the plane of the front face (red colour) and over that of the central cross section (black colour). S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 415 The effect of the radius of the corners of the notch The influence of the radius of curvature at the corners of the notch on the stress field is depicted in Fig.11, in which the equivalent stress at the central cross section along the notch’s perimeter is plotted for three values of r, i.e., r=0.33 mm, 0.66 mm and 1.00 mm. It is seen that the equivalent stress developed is rather marginally influenced by the changes of r. Indeed, even around the notch’s crown, the differences detected are of the order of 5% (Fig.11a). Regarding the maximum principal stress, σ1, plotted in Fig.11b, the maximum difference is detected at the beginning of the rounded corner (point N in Fig.11c) and it is of the order of 15%. (a) (b) Figure 11: The role of the radius of the corners of the notch: The variation of (a) the equivalent stress and an enlarged view of the vicinity around the corner of the notch and (b) the maximum principal stress along one quarter of the notch. The effect of the width of the notch Concerning the role of the width of the notch, it can be said that the qualitative distribution of the equivalent stress, σeqv, in the body of the disc is not drastically different for the two additional cases studied here (i.e., w=1 mm and w=3 mm), in comparison with the reference model with w=5 mm. Indeed, the distributions of σeqv are quite similar to that presented in Fig.9a for the reference model. From a quantitative point of view, however, things are drastically different, especially if atten- 60 85 110 135 160 0 7 14 21 28 σ e qv [ M P a] Along the notch [mm] r=0.33 mm r=0.66 mm r=1.00 mm 60 85 110 135 160 0.0 1.5 3.0 4.5 σ e qv [ M P a] Along the notch [mm] r=0.33 mm r=0.66 mm r=1.00 mm 0 45 90 135 180 0 1 2 3 σ 1 [ M Pa ] Along the notch [mm] r=0.33 mm r=0.66 mm r=1.00 mm N S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 416 tion is focused on the area around the crown of the notch, namely, the locus most important for the determination of KIC, taking into account that crack is expected to start propagating from a point of the notch’s perimeter located in this area. In an attempt to gain better insight on the above-mentioned differences, the distribution of the equivalent stress along the upper half of the y-axis (i.e., the axis parallel to the symmetry axis of the externally applied stress distribution) at the central cross section of the disc is plotted in Fig.12a. It is seen from this figure that approaching the crown of the notch the equivalent stress is about 60% higher for the case of the narrow notch. On the contrary, as one is directed far from the notch’s crown (i.e., moving towards the disc-loading platen contact arc), the equivalent stress is not seriously affected by the differences of the value of w. As far as it concerns the second strategic locus, i.e., that along the perimeter of the notch, it can be seen from Fig.12b, that the values of the equivalent stress are almost insensitive to changes of w, along the segment of the notch that is parallel to the loading axis (see the dotted lines in Fig.12b). However, approaching the mid-point of the notch’s crown the equivalent stress for the narrow notch exceeds dramatically that of the wider notches (see the continuous lines in Fig.12b). (a) (b) Figure 12: The role of the width of the notch: The variation of (a) the equivalent stress along the upper half of the vertical loading y- axis and (b) the maximum principal stress along one quarter of the notch with a detailed view around its corner. 25 30 35 40 45 50 0 50 100 150 200 y [m m ] σeqv [MPa] w=5 mm w=3 mm w=1 mm 50 90 130 170 210 0 7 14 21 28 σ e qv [ M P a] Along the notch [mm] w=5 mm w=3 mm w=1 mm 50 90 130 170 210 0 2 4 6 σ e qv [ M P a] Along the notch [mm] w=5 mm w=3 mm w=1 mm S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 417 (c) (d) Figure 12 (continued): (c) The distribution of the maximum principal stress for the models with w=1 mm and w=3 mm, (d) the variation of the maximum principal stress along one quarter of the notch. Quite similar conclusions are drawn for the maximum principal stress, the distribution of which is plotted in Fig.12c for the cases with w=1 mm and w=3 mm (the distribution for reference case, i.e., that with w=5 mm has been already pre- 0 50 100 150 200 250 0 2 4 6 σ 1 [ M P a] Along the notch [mm] w=5 mm w=3 mm w=1 mm S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 418 sented in Fig.10a). The values of the maximum principal stress, σ1, along one quarter of the notch’s perimeter (again at the central cross section of the disc) are plotted in Fig.12d. Differences of the order of 70% are observed between the model with w=1 mm and that with w=5 mm, indicating once again the paramount influence of the width of the notch on the stress field developed in the vicinity of the crown of the notch. The effect of the length of the notch The last geometrical parameter considered in this study is the length of the notch with respect to the disc’s radius. Again, the distributions of the equivalent stress in the three models (i.e., the models with L=30 mm, 50 mm and 75 mm) are quali- tatively similar to each other (Fig.13a), however the differences along some characteristic paths exhibit significant quantita- tive differences. Indeed, the larger the length of the notch, or in other words the closer the notch’s crown to the boundary of the disc, the larger the equivalent stress developed in the disc, as it can be clearly seen in Fig.13b, where the variation of the equivalent stress along the vertical y-axis of the central cross section is plotted. It is concluded from this figure that the equivalent stress developed at the midpoint of the width of the notch in case L=75 mm is about 60% higher compared to the respective stress developed at the same point when L=30 mm (Fig.13b). The difference is even larger, approximately 75%, at the point where the straight segment of the width meets the starting point of the curved part of the notch’ crown, as it can be seen in Fig.13c. Differences of the same order are observed, also, for the maximum principal stress, σ1. This is clearly seen in Fig.14a where the variation of the maximum principal stress along one quarter of the perimeter of the notch, at the central cross section of the disc, is shown. The difference is about 60% at the midpoint of the width of the notch and it reaches 100% at the point where the straight segment of the width meets the starting point of the curved part of the notch’ crown. Even larger differences, ranging from about 25% to more than 100% (depending on the specific point of the locus) can be detected along the perimeter of the disc (Fig.14b). DISCUSSION AND CONCLUDING REMARKS he role of some geometrical features of the notches machined in circular discs in order to experimentally determine the fracture toughness of brittle rock-like materials was studied, both analytically and numerically. The motive of the study is the difference between the configuration of the theoretical model adopted for the determination of the Stress Intensity Factor (i.e., that of the Crack Straight Through Brazilian Disc - CSTBD) and the actual configuration obtained in the laboratory, when a disc made of rock-like materials is mechanically notched. Indeed, while the theoretical model considers a relatively short “mathematical” crack (i.e., a discontinuity of short length with respect to the disc’s diam- eter, zero distance between its lips and singular tip) the discontinuities mechanically machined are notches of finite width, without singular tip and of length well comparable to the radius of the disc. Based on the results of a recently introduced analytical solution [4,10] a numerical model was designed and validated. The calibration procedure was successful given that the results of the analytic solution and those of the numerical model are in very good agreement ignoring minor discrepancies in the immediate vicinity of the platens-disc interface (attributed to the very complex nature of the analytic solution, which renders convergence of the respective series expressions at specific loci satisfactory only when a large number of terms is used, which is inconvenient from the practical point of view). The validated model was used for a thorough parametric investigation. The investigation highlighted the role of the third dimension (i.e., along the disc’s thickness), which is usually considered of minor importance in two-dimensional studies. Differences between the stress levels at the front face of the disc and at the central vertical section approaching even 20% were detected. The role of the radius of the corners of the notch was proved to be of rather secondary importance. Variations of r of the order of 200% resulted to changes of the equivalent stress lower than about 15%. On the contrary, the role of the length of the notch is proven to be quite catalytic. Indeed, increasing L from 30 to 50 and 75 mm, results to changes of the equivalent stress from 100 to 140 and to 180 MPa, respectively. The parameter, the influence of which on the magnitude of the stress field in the immediate vicinity of the rounded corners of the notch was found to be of paramount importance is the width, w, of the artificially machined notch: Increasing w from 1 to 3 and 5 mm results to changes of the equivalent stress from 140 to 160 and to 200 MPa, respectively. What is to be stressed out, is that a realistic overview of the stress field developed when a notched disc is compressed between the loading platens suggested by either ASTM or ISRM is not singular and the concept of stress intensity does not properly reflect experimental reality. In fact, it is the concept of stress concentration that is to be considered since the stress field at the crown of the notch is intensified by the presence of the notch-shaped discontinuity but by no means it approaches that of a singular field predicted by the theoretical model of the CSTBD approach. T S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 419 (a) (b) (c) Figure 13: (a) The distribution of the equivalent stress in the disc for the three models with L=30, 50, 75 mm; The variation of the equivalent stress along (b) the vertical y-axis and (c) along one quarter of the perimeter of the notch with a detailed view around its corner. To make the above comment clearer, advantage is here taken of the analytic solution described in previous sections. The stress concentration, defined here as the ratio of the equivalent stress developed over the amplitude of the externally applied 15 20 25 30 35 40 45 50 0 50 100 150 y [m m ] σeqv [MPa] L=50 mm L=30 mm L=75 mm 60 100 140 180 0 7 14 21 28 σ e qv [ M P a] Along the notch [mm] L=50 mm L=30 mm L=75 mm 60 100 140 180 0 1 2 3 4 σ e qv [ M P a] Along the notch [mm] L=50 mm L=30 mm L=75 mm S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 420 (a) (b) Figure 14: The variation of the maximum principal stress along one quarter of the notch with a detailed view around its corner (a); and the variation of the equivalent stress around one quarter of the perimeter of the disc (b). applied load, is plotted in Fig.15, versus the width of a notch machined in a disc of radius 50 mm with a notch of length equal to 50 mm, i.e., equal to the radius of the disc. Two strategic points are considered, i.e., point M on the longitudinal axis of symmetry of the notch and point C, the mid-point of the rounded corner of the notch. It is definitely seen that criticality is achieved at point M where the stress concentration is considerably exceeding that at point C, and therefore crack is indeed expected to start propagating along the symmetry axis of the configuration (in full accordance with a series of experimental results [13]). It is seen from Fig.15 that even for widths of the notch smaller than 0.5 mm (a value that is well below the limits of any practical method used to machine the notch) the maximum equivalent stress developed is about 4 times higher compared to the amplitude of the externally applied parabolic stress distribution, and by no means it could be considered as tending to infinite values. Recapitulating, it can be definitely stated that the stress field developed in the disc is strongly influenced by the geometric characteristics of the notch especially in the vicinity of its crowns. Obviously, ignoring the exact geometric shape of the notch may lead to erroneous results concerning the actual severity of the stress field in this critical region and therefore it is quite possible that misleading results about the actual fracture toughness of the disc’s material will be obtained. 0 50 100 150 200 0 7 14 21 28 σ 1 [ M Pa ] Along the notch [mm] L=50 mm L=30 mm L=75 mm 0 50 100 150 200 0 1 2 3 4 σ 1 [ M Pa ] Along the notch [mm] L=50 mm L=30 mm L=75 mm 0 40 80 120 0 30 60 90 σ e qv [ M P a] Along the perimeter [deg] L=50 mm L=30 mm L=75 mm S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 421 Figure 15: The stress concentration (i.e., the ratio of the equivalent stress developed over the amplitude of the externally applied load) versus the width of a notch in a disc of radius 50 mm with a notch of length equal to 50 mm. ACKNOWLEDGEMENTS his research is co-financed by Greece and the European Union (European Social Fund - ESF) through the Opera- tional Programme «Human Resources Development, Education and Lifelong Learning 2014-2020» in the context of the project “Critical assessment and revision of the standard regarding the determination of the Mode I fracture toughness, KIC, of brittle building materials using a combination of innovative experimental techniques, numerical simula- tions and analytical methods” (MIS 5049181). REFERENCES [1] ASTM (2014). E399-12e3 Standard test method for linear-elastic plane-strain fracture toughness KIC of metallic materials. ASTM volume 03.01: Metals-Mechanical Testing; Elevated and Low Temperature Tests; Metallography. [2] ISRM (1995). Suggested methods for determining Mode-I fracture toughness using CCNBD specimens, Int. J. Rock Mech. Min., 32(1), pp. 57–64. [3] Kourkoulis, S.K. and Markidis, Ch.F. (2014). Fracture toughness determined by the centrally cracked Brazilian disc test: Some critical issues in the light of an alternative analytic solution, ASTM Materials Performance & Character- ization, 3(3), pp. 45–86. DOI: 10.1520/MPC20130056. [4] Markides, Ch.F. and Kourkoulis, S.K. (2016). “Mathematical” cracks versus artificial slits: Implications in the determination of fracture toughness, Rock Mech. Rock Eng., 49(3), pp. 707–729. DOI: 10.1007/s00603-015-0794-y. [5] Fowell, R.J., Xu, C. and Dowd, P.A. (2006). An update on the fracture toughness testing methods related to cracked CCNBD specimen, Pure Appl. Geophys., 163, pp. 1046–1057. DOI: 10.1007/s00024-006-0057-7. [6] Kaklis, K., Mavrigiannakis, S., Saltas, V., Vallianatos, F. and Agioutantis, Z. (2017). Using acoustic emissions to enhance fracture toughness calculations for CCNBD marble specimens. Frattura ed Integrità Strutturale, 11(40), pp. 1–17. DOI: 10.3221/IGF-ESIS.40.01. [7] Kourkoulis, S.K., Markides, Ch.F. and Chatzistergos, P.E. (2013). The standardized Brazilian disc test as a contact problem, Int. J. Rock Mech. Min., 57, pp. 132–141. DOI: 10.1016/j.ijrmms.2012.07.016. [8] Markides, Ch.F. and Kourkoulis, S.K. (2013). Naturally accepted boundary conditions for the Brazilian disc test and the corresponding stress field, Rock Mech. Rock Eng., 46(5), pp. 959–980. DOI: 10.1007/s00603-012-0351-x. [9] Muskhelishvili, N.I. (1963). Some Basic Problems of the Mathematical Theory of Elasticity, Noordhoff, Groningen, The Netherlands. [10] Markides, Ch.F. and Kourkoulis, S.K. (2015). The displacement field in a finite circular disc with a central rectangular slit, Procedia Engineering, 109, pp. 257–267. DOI: 10.1016/j.proeng.2015.06.231. -1.5 0.0 1.5 3.0 4.5 0.0 2.5 5.0 7.5 10.0 St re ss co nc en tra tio n Width of slit, b [mm] Midpoint Corner Width of the notch, w [mm] Μ C T S.K. Kourkoulis et alii, Frattura ed Integrità Strutturale, 59 (2022) 405-422; DOI: 10.3221/IGF-ESIS.59.27 422 [11] Savin, G.N. (1970). Stress Distribution Around Holes, NASA, TT F-607, Washington DC. [12] Kourkoulis, S.K., Markides, Ch.F. and Pasiou, E.D. (2015). A combined analytic and experimental study of the displace- ment field in a circular ring, Meccanica, 50(2), pp. 493–515. DOI: 10.1007/s11012-013-9846-0. [13] Kourkoulis, S.K., Pasiou, E.D., Loukidis, A., Stavrakas, I. and Triantis, D. (2021). The determination of mode-I Frac- ture Toughness (by means of the Brazilian-Disc configuration) in the light of experimental evidence provided by novel sensing techniques. SyMBoL final conference, Online, 2-3 September. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice