Microsoft Word - numero_40_art_10 Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 108 Focussed on Recent advances in “Experimental Mechanics of Materials” in Greece The multi-layered ring under parabolic distribution of radial stresses combined with uniform internal and external pressure Christos F. Markides, Ermioni D. Pasiou, Stavros K. Kourkoulis National Technical University of Athens, School of Applied Mathematical and Physical Sciences, Department of Mechanics, 5 Heroes of Polytechnion Avenue, Theocaris Bld., Zografou Campus, 157 73 Athens, Greece ABSTRACT. A recently introduced solution for the stress- and displacement- fields, developed in a multi-layered circular ring, composed of a finite number of linearly elastic concentric layers, subjected to a parabolic distribution of radial stresses, is here extended to encompass a more general loading scheme, closer to actual conditions. The loading scheme includes, besides the parabolic radial stresses, a combination of uniform pressures acting along the outer- and inner- most boundaries of the layered ring. The analytic solution of the problem is achieved by adopting Savin’s pioneering approach for an infinite plate with a hole strengthened by rings. Taking advantage of the results provided by the analytic solution, a numerical model, simulating the configuration of a three- layered ring (quite commonly encountered in practical applications) is validated. The numerical model is then used for a parametric analysis enlightening some crucial aspects of the overall response of the ring. KEYWORDS. Multi-layered ring; Parabolic pressure; Stress and displacement fields; Uniform pressure; Complex potentials. Citation: Markides Ch.F., Pasiou E.D., Kour- koulis S.K., The multi-layered ring under para- bolic distribution of radial stresses combined with uniform internal and external pressure, Frattura ed Integrità Strutturale, 40 (2017) 108-128. Received: 12.02.2017 Accepted: 12.03.2017 Published: 01.04.2017 Copyright: © 2017 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION n a recently published paper [1] the stress- and displacement-fields developed in a circular multi-layered ring were con- sidered both analytically and numerically. The ring was composed by a finite number of concentric layers made of different linearly elastic materials and it was loaded by a parabolic distribution of radial stresses acting along two finite arcs of the outer boundary of the outer-most layer. In general, the configuration of multi-layered rings is common in quite a few engineering applications of increased interest, ranging from Biomechanics (human aorta) to Fluid Mechanics (insulated pipes). Such composite rings are usually loaded either by uniformly distributed internal or external pressure or by a combination of them. In the study by Markides et al. [1] a more complicated loading scheme was considered consisting of radial stresses acting along two finite arcs of the external boundary of the outer-most layer (antisymmetric with respect to the geometric centre of the ring) distributed according to a parabolic law. Such a loading scheme appears in case a circular ring is compressed between either plane or curved jaws (as it is for example the jaws of the device introduced by the International Society for Rock Mechanics (ISRM) [2] for the standardized implementation of the Brazilian-disc test). The solution by Markides et al. [1] is here generalized to confront a more general loading scheme including, besides the parabolic radial stresses, a combination of uniform pressures acting along the outer- and inner-most boundaries of the I Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 109 layered ring. For the analytic solution of the problem, the procedure proposed by Savin [3] for an infinite plate with a hole strengthened by rings is adopted as it was done also in ref. [1]. The specific procedure is based on the complex potentials technique introduced by Muskhelishvili [4]. The analytic solution achieved provides full-field closed-form expressions for both the stresses and the displacements. However, the specific expressions, in spite of their “elegance” and their numerous advantages (related mainly to their generality and the fact that they are expandable to rings that are made up of any number of concentric layers), are somehow “lengthy” and cumbersome. As a result, detailed parametric analyses of the (quite a few) factors influencing the final outcomes become rather difficult. To overcome the above mentioned difficulty, the problem is revisited, also, numerically using the Finite Element method and the commercially available software ANSYS. For the validation of the numerical model, the data obtained from the analytic solution for the case of a ring made up of three concentric layers are used. Then the validated numerical model is used for a parametric study of the role of some critical parameters (i.e., the elastic modulus and thickness of the inter- mediate layer and also the length of the arcs loaded by the parabolic distribution of radial stresses) on the overall stress- and displacement-fields developed in the ring. THEORETICAL CONSIDERATIONS The problem and the basic assumptions onsider a multi-layered hollow cylinder of length w consisting of n concentric hollow cylinders of different thicknesses, perfectly joined together without any gaps at all. The cylinders are assumed to be made of homogeneous, isotropic and linearly elastic materials. The multi-layered cylinder as a whole is in equilibrium under the simultaneous action of three different kinds of loading: inner and outer uniform pressure all along its inner and outer lateral surfaces and a parabolic pressure acting along two finite parts of its outer periphery, antisymmetric with respect to the section’s center. All three types of loading act in the body’s normal cross-section and remain constant along its length. Ignoring its weight, stresses and displacements are to be determined at any point of the multi-layered cylinder. Clearly, the as above described configuration corresponds to a 1st fundamental problem of plane linear elasticity for the body’s cross- section, i.e., the multi-layered circular ring. In this context, Muskhelishvili’s method of complex potentials [4] will be employed for the analytic solution. Mathematical formulation Under the above assumptions, the multi-layered ring is considered lying in the z=x+iy=reiθ complex plane, Fig.1. The origin of the Cartesian reference coincides with the center of the ring. The n constituent concentric rings are numbered in such a way so that they are encountered in the order 1, 2,…, n as one moves from the origin towards the outer perimeter of the multi-layered ring. The arbitrary ring is denoted by the index j (1≤j≤n) and its boundaries are Lj and Lj+1 corresponding to the radii r=Rj and r=Rj+1 (Rj+1> Rj), respectively. In general, it holds that (R2–R1)≠(R3–R2)≠…(Rn+1–Rn). Adjacent constituent rings are perfectly joined together along their common interfaces. The inner boundary of the innermost ring, L1, for r=R1, is subjected to a uniform pressure of magnitude pI>0. The outer boundary of the outermost ring, Ln+1, for r=Rn+1, is under the simultaneous action of a uniform pressure of magnitude pE>0 exerted all over Ln+1, and a parabolic pressure of magni- tude p(θ)>0 acting along two finite arcs of Ln+1, antisymmetric with respect to the ring’s center. Each one of these arcs has a length equal to 2ωoRn+1, where ωo corresponds to the half loaded rim. Particularly, p(θ) is here considered equal to:                           2 2 2 1 2 sin 1 sin sin , max sin 2 2 cos2 o c o o c n o o o F p P P p R w (1) where ϕo is the arbitrary angle formed by the axis of symmetry of the parabolic pressure and x-axis (measured from x-axis in the anticlockwise direction) and F (F>0) is the resultant force due to p(θ). Obviously, for the as above described loading conditions, the multi-layered ring as a whole and each constituent j-ring separately, are in equilibrium and for that configuration the stress- and displacement-fields are to be determined for each j-ring. According to Muskhelishvili, the latter can be implemented by obtaining on each j-ring two analytic functions of the complex variable z, the complex potentials φj(z) and ψj(z), in terms of which stresses and displacements are expressed as [4]: C Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 110       ( ) ( ) 4 ( )j j r j z (2) ii e                  ( )( ) 2( ) ( ) ( ) ( )jj r j j j jr z z z z z (3) ie i               ( )( ) ( ) ( ) ( ) 2 jj r j j j j j u u z z z z (4) Figure 1: The configuration of the problem and definition of symbols. In Eqs.(2-4) prime and double prime denote first and second derivative, bar denotes the complex conjugate and is the real part; in addition, μj denotes the shear modulus and κj is the Muskhelishvili’s constant, equal to 3–4νj or (3–νj)/(1+νj), for plane strain or plane stress respectively, with νj indicating the Poisson’s ratio. Regarding φj(z) and ψj(z), their analytic character on each j-circular ring permits representing them by Laurent series as:          ( ) ( )( ) , ( ) , 1, 2,...,j jk k j jk kz a z z b z j n (5) Therefore the problem reduces to the determination of the coefficients ( )j ka and ( )j kb , a task that will be achieved by ful- filling the boundary conditions of the problem. Before implementing this step, however, it is to be mentioned that it is not possible to obtain directly all ( )j ka and ( )j kb ; namely, for a 1st fundamental problem φj(z) is determined apart from an arbitrary complex constant and an arbitrary real constant C appearing in a term of the form iCz. In this context, ( ) 0 ja and the ima- ginary part,   ( ) 1 ja , of the coefficient ( ) 1 ja remain arbitrary. In the particular problem studied, both are set equal to zero for convenience. What is more, ψj(z) is obtained apart from an arbitrary complex constant so ( ) 0 jb remains arbitrary. Actually, it can be seen that for the present problem only one out of all ( ) 0 jb can be arbitrarily fixed while the remaining ones are pE p(θ) p(θ) x y Ring n (En, νn) R1 Rj Rn+1Rj+1 RnRing j (Ej, νj) Ring 1 (E1, ν1) pΙ Ln Lj Lj+1 Ln+1 L1 ϕο O z=reiθ θ r 2ωο 2ωο Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 111 subjected to determination. Here it is assumed that (1) 0 0.b It is only under these simplifying assumptions that the remaining coefficients ( )j ka and ( )j kb can be determined with the aid of the boundary conditions of the problem. These conditions are classified into two categories: i. The boundary values of stresses on the multi-layered ring’s boundaries as a whole, reading as:        11 , 0r I rp , for  1r R and   0 2 (6)              11 , 0nn r E rp p , for  1nr R and                      o o o o o o o o (7)          11 , 0nn r E rp , for  1nr R and                        2o o o o o o o o (8) ii. The values of stresses and displacements along the interfaces Lj+1, j=1,…,n–1, between the j- and (j+1)-ring, due to their firmly joining, reading as:                11 , j jj j r r r r , for  1jr R and   0 2 (9)             11 , j jj j r ru u u u , for  1jr R and   0 2 (10) Combining Eqs.(3-10) and comparing terms of eikθ of the same k-order, in the expressions resulting, the following equations are obtained involving ( )j ka and ( )j kb : Due to the condition of uniform pressure on L1:                                             (1) (1)2 11 1 (1) (1)3 1 12 2 (1) (1) (1)2 2 1 13 1 1 (1) (1) (1)2 2 4 1 1 13 1 3 (1) ( 1) (1) (1)1 3 1 11 ( 2) 2 (1) ( 1) (1) ( 1) (1)1 1 1 1( 2) 2 0 3 0 0 0 ( 2) 0 I kk k k k k k kk k k k R b p R a R b R a R a b R a R a R b kR a R a R b R a k R a R b    , 4,5,...k (11) where  (1) 1 is the real part of (1) 1a . Due to the condition of uniform and parabolic pressure on Ln+1, (after the factor p(θ), appearing in the first of Eqs.(1) has been expanded in Fourier series):                           ( ) ( )2 11 1 2 ( ) ( )3 1 12 2 2 sin 2 2 2 2sin 0 n n c o o n E o o n n n n P R b p R a R b Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 112 i i                                                        ( ) ( ) ( ) 22 2 1 13 1 1 2 ( ) ( ) ( ) 22 2 4 1 1 13 1 3 2 sin 2 cos21 3 sin 2 sin 2 e 22sin sin 2 cos21 sin 2 sin 2 e 3 22sin ( 2) n n n c o o o n n o o o o n n n c o o o n n n o o o o n P R a R a b P R a R a R b k kR       ii                                            ( ) ( 1) ( ) ( )1 3 1 11 ( 2) 2 2 ( 1) 2 ( ) (1 1 1 sin( 1) sin( 1)1 1 12sin ( 1)cos2 sin( 1) 2sin 2 cos( 1) 1 e e , 4,5, ... 4 1 ( 2) n k n nk k c o o nnk k k o kko o o o o nk n nk P k k a R a R b k k k k k k k k R a k R      ii                                       1) ( ) ( 1) ( ) 1( 2) 2 ( 1) 2 sin( 1) sin( 1)1 1 12sin ( 1)cos2 sin( 1) 2sin 2 cos( 1) 1 e e , 4,5, ... 4 1 k n k n c o o nk k o kko o o o o P k k a R b k k k k k k k (12) where  ( ) 1 n is the real part of ( ) 1 na . Due to the condition for equal action-reaction at facing points of the j- and (j+1)-ring, because of perfect matching of the layers:                                                ( ) ( ) ( 1) ( 1)2 2 1 11 1 1 1 ( ) ( ) ( 1) ( 1)3 3 1 1 1 12 2 2 2 ( ) ( ) ( ) ( 1) ( 1) ( 1)2 2 2 2 1 1 1 13 1 1 3 1 1 ( ) ( ) (2 2 4 1 1 13 1 3 2 2 0 0 3 3 0 j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j R b R b R a R b R a R b R a R a b R a R a b R a R a R b                                                 ) ( 1) ( 1) ( 1)2 2 4 1 1 13 1 3 ( ) ( 1) ( ) ( )1 3 1 11 ( 2) 2 ( 1) ( 1) ( 1) ( 1)1 3 1 11 ( 2) 2 ( ) ( 1)1 1 1 ( 2 0 1, 2,..., 1 0 ( 2) j j j j j j j k j jk k j jjk k k j k j jk k j jjk k k j kk j jk k R a R a R b j n kR a R a R b kR a R a R b R a k R a                                  ( ) ( 1) ( ) 1) ( 1) ( 1) ( 1) ( 1) ( 1)1 1 1 1( 2) , 4,5, ... ( 2) 0 j k j j k j k j k jk j j jk k k k R b R a k R a R b (13) where  ( ) 1 j is the real part of ( ) 1 ja . Due to the condition for the same displacements of facing points of the j- and (j+1)-ring, because of perfect matching of the layers: Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 113                                                                (1) (2) (2)2 22 2 22 0 2 1 1 ( ) ( ) ( 1) ( 1)2 2 1 10 2 0 2 1 ( ) ( ) ( 1) ( 1)1 1 1 1 1 1 11 1 1 1 1 ( ) ( )2 2 1 1 12 2 2 2 0 2 2 0, 2,3, ..., 1 1 1 0 j j j j j j j j j j j j j j j j j j j j j j j j j j j j R a b R a b R a b R a j n R R b R R b R a R b R                                                           ( 1) ( 1)2 2 1 12 2 1 ( ) ( ) ( ) ( 1) ( 1) ( 1)3 1 3 1 1 1 1 1 1 1 13 1 1 3 1 1 1 ( ) ( ) ( ) ( 1) ( 1)3 1 3 3 1 3 1 1 1 1 1 1 13 1 3 3 1 0 3 3 0 j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j a R b R a R a R b R a R a R b R a R a R b R a R a R                                                            ( 1) 3 1 ( ) ( 2) ( ) ( )2 1 11 ( 2) 2 ( 1) ( 2) ( 1) ( 1)2 1 1 11 ( 2) 2 1 ( ) ( 2) ( ) ( ) 1 11 ( 2) ( 1 1 1 0 0 2 j j j k j jk k j j jjk k k j j k j jk k j j jjk k k j j k j jk k j j jjk k k j jk j j k b kR a R a R b kR a R a R b R a k R a R b R a                                                   ) ( 2) ( 1) ( 1) 11 ( 2) , 4,5, ... 1, 2,..., 1 2 0k j jk jj k k k j n k R a R b (14) Properly rearrangement of Eqs.(11-14), yields the following linear systems of equations for the determination of ( )j ka and ( )j kb , for any specific number n of constituent rings:                                                                     (1) (1)2 11 1 ( ) ( )2 11 1 2 ( ) ( ) ( 1) ( 1)2 2 1 11 1 1 1 1 ( ) ( ) ( 1) ( 1)1 1 1 1 1 1 11 1 1 1 2 2 sin 2 2 2 2sin 2 2 0 , 1 1 0 I n n c o o n E o o j j j j j j j j j j j j j j j j j j R b p P R b p R b R b j R R b R R b   1, 2,..., 1n (15) Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 114                                                  (1) ( 2) (2)2 22 2 22 0 2 1 1 ( ) ( ) ( 1) ( 1)2 2 1 10 2 0 2 (1) (1)3 1 12 2 ( ) ( )3 1 12 2 ( ) ( ) ( 1) ( 1)3 3 1 1 1 12 2 2 2 1 2 1 2 2 0 2 2 0, 2,3,..., 1 0 0 0 j j j j j j j j n n n n j j j j j j j j j j j j R a b R a b R a b R a j n R a R b R a R b R a R b R a R b R                        ( ) ( ) ( 1) ( 1)2 2 2 1 1 1 12 2 2 2 , 1, 2,..., 1 0j j j j j j j j j n a R b R a R b (16) (where, obviously, all coefficients appearing in Eqs.(16) are zero) i                                                 (1) (1) (1)2 2 1 13 1 1 (1) (1) (1)2 2 4 1 1 13 1 3 ( ) ( ) ( ) 22 2 1 13 1 1 2 ( ) ( ) ( )2 2 4 1 1 13 1 3 3 0 0 sin 2 cos21 3 sin 2 sin 2 e 22sin sin 2 3 n n n c o o o n n o o o o n n n c n n n o R a R a b R a R a R b P R a R a b P R a R a R b i                                                     2 2 ( ) ( ) ( ) ( 1) ( 1) ( 1)2 2 2 2 1 1 1 13 1 1 3 1 1 ( ) ( ) ( ) ( 1) ( 1) ( 1)2 2 4 2 2 4 1 1 1 1 1 13 1 3 3 1 3 sin 2 cos21 sin 2 e 22sin 3 3 0 o o o o o o j j j j j j j j j j j j j j j j j j j j j j R a R a b R a R a b R a R a R b R a R a R b                                                                       1 ( ) ( ) ( ) ( 1) ( 1) ( 1)3 1 3 1 1 1 1 1 1 1 13 1 1 3 1 1 1 ( ) ( ) ( ) ( 1) ( 1) ( 1)3 1 3 3 1 3 1 1 1 1 1 1 13 1 3 3 1 3 0 3 3 0 0 j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j j R a R a R b R a R a R b R a R a R b R a R a R b     , 1, 2, ..., 1j n (17) Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 115                                              (1) ( 1) (1) (1)1 3 1 11 ( 2) 2 (1) ( 1) (1) ( 1) (1)1 1 1 1( 2) ( ) ( 1) ( ) ( )1 3 1 11 ( 2) 2 2 0 ( 2) 0, 4,5,... sin( 1) sin( 1)1 ( 2) 1 12sin kk k k k k k kk k k k n k n nk k c o o n nnk k k o kR a R a R b R a k R a R b k P k k k kR a R a R b k k     ii                                                ( 1) 2 ( ) ( 1) ( ) ( 1) ( )1 1 1 1( 2) 2 ( 1)cos2 sin( 1) 2sin 2 cos( 1) 1 e e 4 1 sin( 1) sin( 1)1 ( 2) 1 12sin ( 1)cos2 s kko o o o o n k n k nk c o o n n nk k k o o k k k k P k k k R a k R a R b k k k     ii              ( 1) 2 in( 1) 2sin 2 cos( 1) 1 e e 4 1 kko o o ok k k                                                ( ) ( 1) ( ) ( )1 3 1 11 ( 2) 2 ( 1) ( 1) ( 1) ( 1)1 3 1 11 ( 2) 2 ( ) ( 1) ( ) ( 1) ( )1 1 1 1( 2) ( 1) ( 1)1 1 1 ( 2) 0 ( 2) ( 2) j k j jk k j jjk k k j k j jk k j jjk k k j k j k jk j j jk k k j kk j jk k kR a R a R b kR a R a R b R a k R a R b R a k R a                                                            ( 1) ( 1) ( 1) 1 1 ( ) ( 2) ( ) ( )2 1 11 ( 2) 2 ( 1) ( 2) ( 1) ( 1)2 1 1 11 ( 2) 2 1 ( ) ( 2) ( ) ( ) 1 11 ( 2) 0 0 2 j k j j k j j k j jk k j j jjk k k j j k j jk k j j jjk k k j j k j jk k j j jjk k k j R b kR a R a R b kR a R a R b R a k R a R b                                            ( 1) ( 2) ( 1) ( 1) 1 1 11 ( 2) , 1, 2, ..., 1 2 0 4,5,... j k j jk k j j jjk k k j n R a k R a R b k (18) For the 13n–1 unknowns appearing in Eqs.(15-18), i.e., for    (1) ( 2) ( ) 1 1 1, , ..., n ,  (1) ( 2) ( ) 2 2 2, , ..., na a a ,  (1) ( 2) ( ) 3 3 3, , ..., na a a ,     (1) ( 2) ( ) 1 1 1, , ..., na a a , (1) ( 2) ( ), , ..., n k k ka a a ,             (1) (2) ( ) 2 2 2, , ..., n k k ka a a ,  ( 2) (3) ( ) 0 0 0, , ..., nb b b ,  (1) ( 2) ( ) 1 1 1, , ..., nb b b ,     (1) ( 2) ( ) 1 1 1, , ..., nb b b ,     (1) ( 2) ( ) 2 2 2, , ..., nb b b ,     (1) ( 2) ( ) 3 3 3, , ..., nb b b ,     (1) ( 2) ( ) 2 2 2, , ..., n k k kb b b and     (1) ( 2) ( ), , ..., n k k kb b b , there, indeed, exists a number of 13n–1 equations, as expected and thus the problem has a solution. In other words, Eqs.(15), (16), (17) and (18) comprise 2n, 3n–1, 4n and 4n equations, respectively, i.e., a total sum of 13n–1. Thus, ( )j ka and ( )j kb are indeed obtainable from Eqs.(15-18), and in this context φj(z) and ψj(z) can be recapitulated as follows:                  ( ) ( ) ( ) ( ) ( )3 1 ( 2) 1 3 1 ( 2) 5,7,9,... 5,7,9,... ( ) , 1, 2,...,j j j j jk k j k k k k z z a z a z a z a z j n (19) Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 116                  ( ) ( ) ( ) ( ) ( )2 1 3 1 1 32 5,7,9,... 5,7,9,... ( ) , 1, 2,...,j j j j jk k j k k k k z b z b z b z b z b z j n (20) where ‘even’ terms are missing since are all zero. Before concluding it should be mentioned that: (a) by setting pI=0 in the first of Eqs.(15), the present formulae [Eqs.(15-20)] provide the solution in case internal uniform pressure is zero and the multi-layered ring is only under the action of external uniform pressure pE and external parabolic pressure p(θ) (along two finite arcs of its outer periphery, antisymmetric with respect to x-axis), (b) by setting pE=0 in the second of Eqs.(15), the present formulae [Eqs.(15-20)] provide the solution in case external uniform pressure is zero and the multi-layered ring is only under the action of internal uniform pressure pI and external parabolic pressure p(θ), (c) by setting pI=0 in the first of Eqs.(15) and pE=0 in the second of Eqs.(15), the present formulae [Eqs.(15-20)] provide the solution in the case the multi-layered ring is only under the action of external parabolic pressure p(θ), (d) by setting Pc=0 in Eqs.(15, 17 and 18), the present formulae [Eqs.(15-20)] provide the solution in the case external parabolic pressure is zero and the multi-layered ring is only under the action of internal uniform pressure pI and external uniform pressure pE, (e) by setting Pc=0 in Eqs.(15, 17 and 18) and pI=0 in the first of Eqs.(15), the present formulae [Eqs.(15-20)] provide the solution in the case internal uniform and external parabolic pressure are zero and the multi-layered ring is only under the action of externally imposed pressure pE and (f) by setting Pc=0 in Eqs.(15, 17 and 18) and pE=0 in the second of Eqs.(15), the present formulae [Eqs.(15-20)] provide the solution in the case external uniform and parabolic pressure are zero and the multi-layered ring is only under the action of internal uniform pressure pI. Particularly, in all last three cases (d, e, f), i.e., when only uniform pressure over one or both the two boundaries (internal and external) of the multi-layered ring is present, the complex potentials are significantly simplified to:        ( ) ( ) 1 1( ) , ( ) , 1, 2,..., j j j j b z z z j n z (21) where  ( ) 1 j and  ( ) 1 jb are now just provided by only Eqs.(15), for Pc = 0, and for: pI ≠ 0, pE ≠ 0 in case (d), pI = 0, pE ≠ 0 in case (e) and pI ≠ 0, pE = 0 in case (f). Figure 2: Distribution of stresses along the diameter corresponding to the symmetry axis of the parabolic load, for three loading schemes. Internal pressure=External pressure=10 MPa, Pframe=5 kN, ωo=10o, Pc=43.1 MPa, n=3. Composite ring’s thickness w=0.01 m. Inner ring (1): R1≤r≤ R2, R1=0.025 m, R2=0.030 m, E1=30 GPa, ν1=0.25. Intermediate ring (2): R2≤r≤ R3, R2=0.030 m, R3=0.034 m, E2=10 GPa, ν2=0.30. Outer ring (3): R3≤r≤ R4, R3=0.034 m, R4=0.050 m, E3=3 GPa, ν2=0.35). -100 0 100 200 0.025 0.030 0.035 0.040 0.045 0.050 P ol ar s tr es se s [M P a] R1 ≤ r ≤ R4 [m] Parabolic & internal pressure Parabolic & external pressure All three kinds of loading σθ σr Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 117 In order to demonstrate the capabilities of the solution introduced, the distribution of the radial and tangential stresses along the diameter corresponding to the symmetry axis of the parabolic load, is plotted in Fig.2, for three combinations of the constituent loading schemes, i.e., parabolic distribution of radial stresses and uniform internal pressure, parabolic distribution of radial stresses and uniform external pressure, and finally all three schemes simultaneously. The discontinuity of the distribution of transverse stress at the interfaces of the rings is the main characteristic of the plots. As it was expected, this discontinuity is stronger for the case of simultaneous action of the parabolic distribution and external uniform pressure. NUMERICAL ANALYSIS - APPLICATION FOR n=3 CONSTITUENT RINGS s a next step, a parametric numerical analysis was decided to be carried out, in order to explore the influence of various factors on the displacement- and stress-fields developed in a multi-layered ring. The ring with three layers under plane strain conditions was decided to be considered, due to its increased practical interest. The principal advantage of the numerical scheme is that it models the “multi-layered ring - jaw” complex (as an integrated elastic system) rather than the isolated ring. Before the main numerical study, the model was properly validated as it is described in next section. Construction and validation of the reference numerical model For the proper validation of the numerical model advantage was, initially, taken of the analytic solution for the most complex loading case, i.e., the partial parabolic pressure. For symmetry reasons (materials and geometry) only a quarter of the ring- jaw complex (0ο≤θ≤90ο) was constructed with the aid of ANSYS software. The geometry of the reference model is shown in Fig.3. For the radii of the layers it holds that Rj=25, 30, 34, 50 mm, for 1≤j≤4. All three interfaces (i.e., between the con- stituent rings and between the jaw and the outer ring) were taken into account. The model was meshed using the PLANE182 element. Contact elements (CONTA172 and TARGE169) were used to model the interfaces. The interfaces between the successive rings were considered as mutually bonded without separation. A very small (tending to zero) value was assigned to the coefficient of friction at the interface between the jaw and the ex- ternal surface of the outer ring, in order to minimize (eliminate) friction stresses. This assumption was dictated by the fact that the analytic solution, described in the previous section, considers an isolated ring under exclusively radial stresses without any friction stresses at all. The analysis performed was linearly elastic. The moduli of elasticity and the Poisson’s ratios assigned to each one of the rings were Ej=30, 10, 3 GPa and νj=0.25, 0.30, 0.35, respectively, for 1≤j≤3. The jaw was considered to be made of steel (Ejaw=210 GPa and νjaw=0.30). Figure 3: The reference numerical model and detailed views of the finally adopted mesh. Due to the model’s symmetry the following boundary conditions were adopted: i. The motion along the y-axis (loading axis) of all the nodes of the lower line of the ring was restricted. A Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 118 ii. All the nodes of the left (vertical) line of the ring as well as the nodes of both sides of the jaw were restricted of shifting horizontally (i.e., normally to the loading axis). iii. A uniform vertical displacement on the nodes of the upper line of the jaw was imposed. The optimum element size for the meshing of the model was chosen based on a proper convergence analysis. The variation of the polar stresses developed along “Path 0”, i.e., the outer periphery of the outermost ring, and especially around the position where the major differences are observed is plotted in Fig.4a. The embedded graph corresponds to the stresses along the whole length of “Path 0”. In addition, the total stress applied on point A (see Fig.3) is plotted in Fig.4b. Both graphs indicate that convergence is satisfactory for about 20000 elements. In order to obtain more accurate results (and given that the CPU-time consumption is not prohibitive), a number of about 38000 elements was chosen for the final meshing. Figure 4: (a) The variation of the polar stresses developed along “Path 0” (the embedded graph corresponds to the whole length of “Path 0”) and (b) the total stress applied on point A (see Fig.3) for various numbers of elements. Before proceeding to the parametric study, the numerical model was used to check the assumption of the analytic solution concerning the parabolic distribution of radial stresses along the loaded arc of the ring. This was considered absolutely necessary because the specific assumption is based on the solution of a different problem, i.e., that of the “compact disc - jaw” contact problem [5] and therefore adopting a similar distribution for the “ring - jaw” problem appears perhaps somehow arbitrary. In this context, the vertical displacement of the upper node (point A in Fig.3) of the ring (which is, in fact, the only point of the ring in contact with the jaw before the onset of loading) was determined according to the analytic Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 119 solution described in previous section, assuming now that ωo is not arbitrarily prescribed (as it was considered in Eqs.(1)), but it is rather provided, together with Pc, by the relevant compact disc-jaw contact problem [5].Then a uniform displacement was applied at the steel jaw (simulating in fact the load exerted by the loading frame during the laboratory experiments with the device suggested by ISRM for the standardized Brazilian-disc test), the magnitude of which was determined by demanding that it should result to a vertical displacement of point A equal to the respective one provided by the analytic solution. The results of the numerical analysis for the distribution of pressure along the common contact arc, developed at the ring-jaw interface are plotted in Fig.5, in conjunction with the parabolic pressure applied on the ring during the theoretic- al solution. It is seen from this figure that the two distributions are in relatively good agreement, despite some differences concerning their amplitude (equal to about 140 MPa for the analytic solution and 122 MPa for the numerical one) and also their flatness (the numerically obtained distribution appears somehow flatter). These differences should be expected, given that the parabolic distribution is an approximation of the actual “cyclic” distribution [5, 6] (developed during the contact of two elastic cylindrical bodies), which was adopted in an attempt to achieve closed-form analytic expressions for the dis- tribution of stresses a displacements all over the disc [7, 8]. (It is here recalled, that for the actual “cyclic” distribution analytic expressions for the stress field cannot be straightforwardly obtained [5]). Figure 5: The applied pressure applied on the ring based on both the analytical and the numerical analysis. Besides the above mentioned slight differences, it is quite astonishing that the analytic and the numerical approaches provide almost identical results for the “total” pressure applied on the ring indicating the static equivalency of the two solutions. Even more astonishing is the fact that the length of the contact arcs obtained by the two solutions (numerical and analytic) is again almost identical: The contact arc 2ωο determined from the numerical model equals to 24.67o while the respective one of the analytical solution is equal to 24.59o. The same excellent agreement was observed, also, for the variation of displacements and stresses along strategic paths of the model, as it can be seen, for example, in Fig.6, in which analytical versus numerical results concerning the vertical displacement (Fig.6a) and the polar stresses (Fig.6b) are plotted for “Path 1”, i.e., along the diameter corresponding to the axis of loading symmetry. The differences observed are quite marginal. For example, for the displacements, the difference does not exceed 4% (Fig.6a) while for the polar stresses the differences are almost tending to zero (Fig.6b). As an additional validation step, results of the numerical analysis were compared against the respective ones analytically obtained, also, for some combinations of the three loading schemes considered (parabolic pressure, external uniform pres- sure and internal uniform pressure). For this specific step the jaw was removed and the multi-layered ring was modeled as an isolated elastic body (for convenience). The following loading schemes were imposed on the ring: Case 1: Pressure was applied on the nodes of the outer perimeter of the outermost ring along a contact length 2ωο which was assumed equal to 20o. The parabolic pressure variation corresponds to a statically equivalent force equal to 5 kN applied by the loading frame. Case 2: Internal pressure equal to 10 MPa was applied on all nodes along the inner perimeter of the three-layered ring. Case 3: External pressure equal to 10 MPa was applied on all nodes along the outer perimeter of the three-layered ring. 0 40 80 120 160 76 78 80 82 84 86 88 90 σ r [M P a] θ [degrees] Analytical Numerical Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 120 Case 4: Combination of Cases 1 and 2. Case 5: Combination of Cases 1 and 3. Case 6: All three types of loading, i.e., combination of Cases 1, 2, 3. The normal polar stresses (both radial and transverse) obtained by the numerical model for all six loading cases mentioned above are presented in Fig.7 in juxtaposition to the respective analytical results. The plots are realized, as in Fig.6, along the diameter corresponding to the symmetry axis of the loading distribution. It is clear that the correlation between analytic and numerical results is quite satisfactory, for both stresses and for all loading schemes without any exception. Figure 6: Analytical vs. numerical results for the (a) vertical displacement and (b) the polar stresses developed along “Path 1”. Recapitulating the data of the validation procedure, it can be definitely concluded that the numerical model approaches the analytic solution in an excellent manner. Therefore, it is safe for the model to be used for a thorough parametric ana- lysis of the role of the (quite a few) factors, which are expected to crucially affect the overall mechanical response of the multi-layered ring, in case it is subjected to various combinations of the three constituent loading schemes considered (internal and external pressure and parabolic variation of radial stresses). Moreover, the results of the present section, besides validating the numerical model, indicated that the actual stress distribution along the contact arc at the ring-jaw interface is Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 121 rather insensitive to whether the disc squeezed between the ISRM jaws is compact or hollow, i.e., ring (either it is homo- geneous or it consists of successive layers with different mechanical properties), at least for low to moderate values of the externally applied load and provided relatively thick-walled rings are considered. Figure 7: Analytical vs. numerical results for (a) simple and (b) complex loading schemes. Blue lines correspond to radial stress while red ones correspond to transverse stress. Parametric numerical analysis The as above exhaustively validated numerical model is now used for a parametric analysis. In this direction, a series of new models was constructed for the most complicated loading constituent, i.e., the parabolic pressure. Three parameters were taken into account in the present study, the influence of which is expected to be more crucial. These parameters are: -100 0 100 200 0.025 0.030 0.035 0.040 0.045 0.050 P ol ar s tr es se s [M P a] R1 ≤ r ≤ R4 [m] Parabolic & internal-analytical Parabolic & external-analytical All three-analytical Parabolic & internal-numerical Parabolic & external-numerical All three-numerical -80 0 80 160 0.025 0.030 0.035 0.040 0.045 0.050 P ol ar s tr es se s [M P a] R1 ≤ r ≤ R4 [m] Parabolic-analytical Internal-analytical External-analytical Parabolic-numerical Internal-numerical External-numerical (a) (b) Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 122 a. The modulus of elasticity of the intermediate ring, E2. b. The thickness of the intermediate ring, t2=R3-R2, keeping the ring’s outer diameter constant, i.e., increasing t2 results to decrease of the thickness of the outermost ring, t3=R4-R3. c. The contact length, or equivalently the length of the loaded arc, 2ωo, assuming that the overall externally imposed load is kept constant. In this case the jaw was not taken into account. Concerning the distribution of the pressure applied along each contact arc, it was determined by taking advantage of Eqs.(1). The numerical values assigned to each one of the three parameters studied are summarized in Tab. 1. In addition, the variation of the applied pressure for the values considered for the length of the half contact arc ωo is plotted in Fig.8. As it is expected, by decreasing the length of the loaded arc (while keeping the overall load constant) the amplitude of the distri- bution increases, leading to increased stress concentration in the immediate vicinity of the loaded arc, which in turn increases the possibility for premature failure in this region of the ring. This aspect should be carefully taken into account in case the external layer of the ring is made of relatively brittle materials. Parameter Values E2 [GPa] 5.0 7.5 12.5 15.0 17.5 20.0 t2 [mm] 1 2 3 5 6 ωo [deg] 1 3 5 7 9 11 Table 1: The numerical values assigned to each parameter. Figure 8: The variation of the pressure applied on the ring for various ωo-values. RESULTS AND DISCUSSION The influence of the elastic modulus of the intermediate ring, E2 he displacement along the loading axis (i.e., y-axis) developed all over the ring’s area is presented in Fig.9a for the reference model. Its variation along y-axis for various E2 values is exhibited in Fig.9b. It is interesting to observe that, although the values of E2 vary within a relatively wide range (from -50% to +100% of the respective value of the reference model), the displacements along y-axis are more or less insensitive to this variation. The maximum difference does not exceed 3% and it is observed for E2=20 GPa. On the contrary, the stresses developed in the ring are much more sensitive to the changes of E2. To enlighten this point (and considering that the polar reference system is much more convenient for the specific application), the variation of the transverse stress (σθ) along the loaded diameter is plotted in Fig.10 for all E2 values studied. It is seen that the influence of E2 is very strong, especially at the interface points between the constituent rings. As it is perhaps expected, the “stress jump” 0 450 900 1350 1800 76 78 80 82 84 86 88 90 σ r [M Pa ] θ [degrees] ω=12.3 deg ω=11.0 deg ω=9.0 deg ω=7.0 deg ω=5.0 deg ω=3.0 deg ω=1.0 deg ωο= ωο= ωο= ωο= ωο= ωο= ωο= T Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 123 is amplified as the difference of the modulus of elasticity of ring 2 from the elastic modulus of the adjacent ring increases. Therefore, the variation of σθ from ring 1 (innermost) to ring 2 (intermediate) is smoother when E2=20 GPa while from ring 2 to ring 3 (outermost) the variation is smoother when E2=5 GPa. It is also noted that σθ is negative all along rings 2, 3 (along the specific path) for all E2 values studied, while it changes sign somewhere along ring 1. The variation of the polar stresses developed along the interfaces is quite interesting. For the interface between rings 1, 2 and more specifically along the outer periphery of ring 1 (Fig.11a), σθ is affected more severely. When the value of E2 is between the elastic moduli of the two surrounding rings, the variation of σθ is again very smooth with small variations. For E2<10 GPa the transverse stress becomes larger as θ increases while the opposite is true for E2>10 GPa. It is also seen from the same figure that the values of σrθ developed along the interface are of the same order of magnitude with σθ. Figure 9: (a) The distribution of the uy component of the displacement field in [m]. (b) The variation of uy along “Path 1” for various values of the elastic modulus of the intermediate ring (E2). The colours of the lines of the graph are in agreement to the colour code of the embedded figure. Figure 10: The transverse stress along “Path 1” for various values of the elastic modules, E2, of the intermediate ring. The colours of the lines of the graph are in agreement to the colour code of the embedded figure. Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 124 As the modulus of elasticity of ring 2 increases, σrθ also increases, attaining its maximum value when θ~65o. The radial stress also attains its maximum value when θ is equal to about 65o, although a slight “move” to a lower value of θ is observed as E2 increases. It is also to be noted that the changes of σr are not significant all along the outer periphery of ring 1. The variation of stresses along the inner periphery of ring 3 (interface between rings 2 and 3) is plotted in Fig.11b. It is clear that the variation of both radial and shear stresses is qualitatively similar to the respective one along the outer peri- phery of ring 1 but their magnitudes are lower. The radial stress attains its maximum value when θ is equal to about 60o and the differences are mainly detected when θ varies from about 75o to 90o. Concerning the shear stress, it attains its maximum value when θ~70o. In addition, it seems that the exact position of the maximum shear stress “moves” slightly to larger θ values with increasing E2. On the other hand, the variation of the transverse stress is monotonous all along the specific path, contrary to what it was observed along the outer periphery of ring 1. Figure 11: (a) The variation of the polar stresses along “Path 2” (outer periphery of the innermost ring) for various E2-values. (b) The variation of the polar stresses along “Path 3” (inner periphery of the outermost ring) for various E2-values. The influence of the thickness of the intermediate ring, t2 The effect of the thickness of the intermediate ring on the variation of stresses along the outer periphery of ring 1 (Fig.12a) is very similar to the one presented in Fig.11a (where the parameter studied was the modulus of elasticity of ring 2) both qualitatively and quantitatively. For t2=4 mm (reference model) the variation of σθ is again smoother. For t2<4 mm σθ in- creases with increasing θ while for t2>4 mm σθ increases with decreasing θ. The shear stress, on the other hand, increases continuously with increasing t2 and its maximum value is again detected around θ~65o independently of the values of t2 considered here. Increase of t2 leads also to increase of the radial stress developed. Again its maximum value seems to “move” slightly towards larger θ-value as t2 increases. Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 125 The variation of stresses along the inner periphery of ring 3 is exhibited in Fig.12b and it is similar to the respective one in Fig.11b. However, and despite this qualitative similarity, the thickness of the intermediate ring affects more the stress field along the inner periphery of the outermost ring 3. The differences are generally larger, especially around θ=0o and θ=90o. This is true for both radial and transverse stresses. On the contrary, the shear stress attains its maximum value at θ~65o when t2=1 mm while as t2 increases the maximum shear stress appears at slightly larger values of θ. Figure 12: The variation of the polar stresses developed along (a) “Path 2” - outer periphery of innermost ring and (b) “Path 3” - inner periphery of outermost ring for various values assigned to the thickness of the intermediate ring (t2). Another interesting point is that the contact length between the outermost ring and the jaw, as well as the pressure applied along it, are only slightly influenced by the thickness of the intermediate ring, as it can be seen in Fig.13. The lower t2 is, the smaller the contact length is (i.e., for t2=1 mm, 2ωo=23.34o while for t2=6 mm, 2ωo=24.67o) and the lower the maximum pressure applied on the ring is (keeping the externally applied displacement on the jaw constant). The role of the contact length, 2ωo Keeping the “total” pressure applied constant and at the same time narrowing the contact length ωo, results to interesting changes of the polar stresses developed along the interfaces. It is recalled here that the specific parameter is studied for an isolated ring rather than the ring-jaw complex. Therefore the contact length for the reference model is here set equal to that obtained during the validation procedure, which was equal to 2ωο=24.60ο. Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 126 Figure 13: The variation of the total pressure applied on the ring for various t2-values. Figure 14: The variation of the polar stresses developed along (a) “Path 2” - outer periphery of innermost ring and (b) “Path 3” - inner periphery of outermost ring for various values assigned to the half contact length (ωo). 0.00E+00 3.50E+07 7.00E+07 1.05E+08 1.40E+08 76 78 80 82 84 86 88 90 T ot al p re ss ur e [P a] 0o ≤ θ ≤ 90o 4 mm 1 mm 2 mm 3 mm 5 mm 6 mm Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 127 Indeed, along the outer periphery of ring 1 (Fig.14a) the transverse stress σθ seems to be almost independent of the contact length for ωo≥5o. The changes are noticeable only for ωo <5o and the respective variations are almost similar to the aforementioned ones (attaining of course larger values). The same is more or less true concerning the shear stress. As the contact length becomes narrower, the shear stress increases (the changes are larger as ωo decreases) attaining its maximum value at about θ~70o. Similar conclusions are drawn, also, for the variation of the radial stresses. The variation of the transverse stress σθ along the inner periphery of ring 3 is exhibited in Fig.14b. As ωo decreases the changes of σθ are more pronounced and the maximum difference recorded is equal to about 15% for ωo=3o (in comparison to the reference model for which ωo=12.33o). In case ωo=1o, almost half of the periphery is affected and the difference now reaches about 30%. The radial stress is also strongly affected, reaching 50% larger value at θ=90o when ωo=1° compared to the respective stress when ωo=12.33o. Concerning the shear stress, it is again larger as ωo decreases attaining its maximum value around θ=70o. CONCLUSIONS n analytic solution for the mechanical response of a multi-layered ring under a combined loading scheme, including both uniform internal and external pressure and also parabolic distribution of radial stresses along two finite arcs of its external periphery was introduced. The solution was obtained in series form with the aid of Savin’s pioneering approach for the infinite plate with a hole strengthened by rings. The main advantage of the analytic procedure described is that it is capable of providing the full-field stress- and displacement-fields for any combination of the three constituent loading schemes considered. As an application the three-layered ring was studied due to its increased practical importance. In order for a closed form full-field solution to be achieved, some critical assumptions were adopted, besides the familiar ones related to the linearity, homogeneity and isotropy of the ring’s constituent materials. These assumptions are related to the ideal bonding between the ring’s successive layers, the absence of any kind of friction and the distribution of radial stresses along the loaded arcs, which was considered of parabolic nature, in accordance with the respective distribution for an intact disc compressed between the jaws suggested by ISRM for the implementation of the Brazilian-disc test. The analytic solution was used to validate a flexible numerical model, in the direction of implementing thorough parametric studies. The advantage of the numerical model was the fact that it considered the ring and the jaw as a system of elastic bodies in equilibrium, approaching experimental reality in a much more accurate manner. The model was exhaustively validated taking advantage of the analytic results. During the validation process it was concluded that the actual stress distribution along the contact arc at the ring-jaw interface is more or less independent of some crucial geometric character- istics. In fact, it was proven that the specific distribution in case the body compressed between the ISRM jaws is intact (disc) does not exhibit significant differences compared to the respective distribution for a body weakened by the presence of a circular hole (irrespectively of whether it is a homogeneous ring or a ring consisting of successive layers with different mechanical properties). In general, the assumptions adopted during the formulation of the analytic problem for the loading scheme of the isolated ring were in excellent agreement with the conditions along the ring-jaw contact arc obtained by the numerical study of the ring-jaw complex. The parametric study revealed quite a few interesting aspects concerning the influence of three critical parameters on the stress- and displacement-fields developed in the multi-layered ring. Generally speaking, this influence is of complicated, non- monotonous nature, rendering further parametric analysis with simultaneous variation of these parameters unavoidable. REFERENCES [1] Markides, Ch.F, Pasiou, E.D., Kourkoulis, S.K., The multi-layered ring under parabolic pressure, Engineering Trans- actions, 64(4) (2016) 433-440. [2] ISRM, Suggested methods for determining tensile strength of rock materials, International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 15(3) (1978) 99–103. [3] Savin, G.N., Stress concentration around holes, Pergamon Press, Oxford (1961). [4] Muskhelishvili, N.I., Some basic problems of the mathematical theory of elasticity, 4th edition, P. Noordhoff Groningen, The Netherlands (1963). [5] Kourkoulis, S.K., Markides, Ch.F., Chatzistergos, P.E., The standardized Brazilian disc test as a contact problem, International Journal of Rock Mechanics and Mining Sciences, 57 (2012) 132-141. A Ch. F. Markideset alii, Frattura ed Integrità Strutturale, 40 (2017) 108-128; DOI: 10.3221/IGF-ESIS.40.10 128 [6] Timoshenko, S.P., Goodier, J.N., Theory of elasticity, McGraw-Hill, New York (1970). [7] Markides, Ch.F., Kourkoulis, S.K., The stress field in a standardized Brazilian disc: The influence of the loading type acting on the actual contact length, Rock Mechanics and Rock Engineering, 45(2) (2012) 145-158. [8] Kourkoulis, S.K., Markides, Ch.F., Chatzistergos, P.E., The Brazilian disc under parabolically varying load: Theoretical and experimental study of the displacement field, International Journal of Solids and Structures, 49(7-8) (2012) 959-972. << /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