Microsoft Word - numero_41_art_61.docx V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 491 Analysis of longitudinal cracked two-dimensional functionally graded beams exhibiting material non-linearity Victor Rizov Department of Technical Mechanics, University of Architecture, Civil Engineering and Geodesy, 1 Chr. Smirnensky blvd., 1046 – Sofia, Bulgaria. v_rizov_fhe@uacg.bg ABSTRACT. An analytical study of longitudinal fracture in two-dimensional functionally graded cantilever beam configurations is carried-out with taking into account the non-linear behavior of material. A longitudinal crack is located arbitrary along the beam cross-section height. The material is functionally graded along the width as well as along the height of beam. The external loading consists of a bending moment applied at the free end of lower crack arm. Fracture is studied in terms of the strain energy release rate by considering the beam complementary strain energy. The solution derived is verified by analyzing the longitudinal crack with the help of the J-integral. The distribution of J-integral value along the crack front is studied. The effects of crack location, material gradients and non-linear behavior of material on the fracture are elucidated. The analysis reveals that the material non-linearity has to be taken into account in fracture mechanics based safety design of structural members and components made of two-dimensional functionally graded materials. KEYWORDS. Functionally graded beam; Crack; Material non-linearity; Analytical solution. Citation: Rizov, V., Analysis of longitudinal cracked two-dimensional functionally graded beams exhibiting material non-linearity, Frattura ed Integrità Strutturale, 41 (2017) 491-503. Received: 23.04.2017 Accepted: 28.05.2017 Published: 01.07.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 he functionally graded materials are new inhomogeneous materials produced by mixing of two constituent materials in different proportions in order to satisfy the requirements of material properties in different parts of a structural member [1-12]. Since the composition of constituent materials varies continuously along one or more spatial coordinates during manufacturing without sudden changes of material properties, the interfaces and stress concentrations are largely avoided in contrast to laminated composite materials [13-15]. Besides, the performance of functionally graded structural members and components to the external load can be optimized by tailoring the variation of material properties in a desired way during manufacturing. Therefore, the functionally graded materials have been widely applied as advanced structural materials in aeronautics, electronics, nuclear energy, optics, engineering and biomedicine for T V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 492 the last three decades [16]. Fracture behavior is crucial for structural integrity and functionality of members and components made of functionally graded materials. Therefore, considerable attention has been paid to fracture in these novel materials by the international research community [17-22]. The present paper is concerned with a longitudinal fracture analysis of a two-dimensional functionally graded cantilever beam configuration with taking into account the non-linear behavior of material. Fracture is analyzed in terms of the strain energy release rate. The material is functionally graded along the height as well as along width of the beam. The mechanical behavior of material is described by a power-law stress-strain relation. An additional fracture analysis is carried-out by applying the J-integral approach and the results obtained are compared with the strain energy release rate. The effects of material gradients, crack location and material non-linearity on the fracture behavior are elucidated. ANALYSIS OF THE STRAIN ENERGY RELEASE RATE he present paper analyses the longitudinal fracture in the two-dimensional functionally graded cantilever beam configuration shown schematically in Fig. 1. Figure 1: The geometry and loading of the two-dimensional functionally graded cantilever beam. A longitudinal crack of length, a, is located arbitrary along the beam cross-section height (it should be mentioned that the present paper is motivated also by the fact that functionally graded materials can be built up layer by layer [16] which is a premise for appearance of longitudinal cracks between layers). The lower and upper crack arms have different thicknesses denoted by 1h and 2h , respectively. The beam is loaded by one bending moment, yM , applied at the free end of lower crack arm (Fig. 1). Thus, the upper crack arm is free of stresses. The beam has a rectangular cross-section of width, b, and height, 2h. The beam is clamped in section, D. The beam non-linear mechanical behavior is described by the following power-law stress-strain relation [23]: mB  (1) where σ is the longitudinal normal stress, ε is the longitudinal strain, B and m are material properties. Zircona-titanium is one of the functionally graded materials which exhibit non-linear mechanical behavior. Data for mechanical properties of zircona-titanium can be found in [24]. The material is functionally graded along the width as well as along the height of the beam (Fig. 1). The distribution of material property, B, in the beam cross-section is expressed as a function of 3y and 3z by the following quadratic law: T V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 493 2 2 3 3 0 1 22 2 4 y z B B B B b h    , (2) where 3 2 2 b b y   , 3h z h   . (3) In (2), 0B is the value of B in the centre of the beam cross-section, 1B and 2B are material properties which govern the material gradients along axes, 3y and 3z , respectively (Fig. 1). Figure 2: Distribution of the material property, B, in the beam cross-section. It should be mentioned that Eq. (2) is suitable for two-dimensional functionally graded beam configurations which have high mechanical properties at the lower, upper and lateral surfaces of the beam as shown in Fig. 2. The fracture behavior is studied in terms of the strain energy release rate. For this purpose, an elementary increase, da , of the crack length is assumed (the external load is kept constant). The strain energy release rate, G, is expressed via the changes of external work, extdW , and strain energy, dU , as extdW dUG dA   (4) where the increase of crack area is dA bda (5) The change of external work is written as * extdW dU dU  (6) where *dU is the change of complementary strain energy. By combining of (4), (5) and (6), the strain energy release rate can be expressed as *dU G bda  (7) It should be noted that the present analysis holds for non-linear elastic behavior of material. The analysis is applicable also for elastic-plastic behavior of material if the beam considered undergoes active deformation, i.e. if the external loading increases only [25]. It should also be mention that the present study is based on the small strains assumption. V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 494 Figure 3: The free end of the lower crack arm (. 1 1n n . is the neutral axis). The beam complementary strain energy, *U , is written as   22 2 2 * 01 2 2 2 2 1 * 0 * 1 1 dzdyualdzdyuaU h h b b U h h b b L                             (8) where * 0Lu and * 0Uu are the complementary strain energy densities in the lower crack arm and the un-cracked beam portion  lxa  3 , respectively. Axes, 1y and 1z , are shown in Fig. 3, 2y and 2z are the centroidal axes of the cross- section of the un-cracked beam portion. In principle, the complementary strain energy density is equal to the area OQR that supplements the area OPQ, enclosed by the stress-strain curve, to a rectangle (Fig. 4). Therefore, the complementary strain energy density in the lower crack arm is written as LL uu 0 * 0   (9) where the strain energy density, Lu0 , is equal to the area OPQ    0 0 )( du L (10) From (1), (9) and (10), the complementary strain energy density can be obtained as 1 1 * 0    m m Bu m L  (11) The strain distribution is analyzed by applying the Bernoulli’s hypothesis for plane sections since the span to height ratio of the beam under consideration is large. Concerning the application of the Bernoulli’s hypothesis in the present study, it should also be noted that since the beam is loaded in pure bending (Fig. 1) the only non-zero strain is the longitudinal strain, ε. Therefore, according to the small strain compatibility equations, ε is distributed linearly along the height of the beam cross-section. Thus, ε in the lower crack arm is written as V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 495  1 1 11nz z   (12) where 1 1nz is the coordinate of the neutral axis, 1 1n n , (the neutral axis shifts from the centroid since the material is functionally graded (Fig. 3)), 1 is the curvature of the lower crack arm. Figure 4: Schematic of a non-linear stress-strain curve ( 0u and * 0u are the strain energy density and the complementary strain energy density, respectively). The following equilibrium equations of the lower crack arm cross-section (Fig. 3) are used in order to determine 1 1nz and 1 : 1 2 2 1 1 1 1 2 2 h b h b N dy dz                 (13) 1 2 2 1 1 11 1 2 2 h b y h b M z dy dz                 (14) where 1N and 1yM are the axial force and the bending moment, respectively. It is obvious that 1 0N  , 1y yM M (15) By using (2), the distribution of material property, B, in the lower crack arm cross-section is written as  22 11 0 1 22 2 4 r zy B B B B b h     (16) where 1 2 2 b b y   , 1 1 1/ 2 / 2h z h   , 1 2 h r h  . (17) By substituting of (1), (12) and (16) in (13) and (14), one derives V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 496       2 1 10 1 1 1 2 1 1 1 221 3 1 1 m m m m mB b B b B brN m m h m                       12 2 1 12 1 11 2 1 22 2 1 2 1 m nm m m m zB br m mh                   (18) 3 3 2 2 12 1 1 1 2 1 22 21 3 2 f q f q f q f q m nq q q qzB b q f q f qh                                  2 1 1 1 2 f q f q n q qz f q               2 2 21 1 2 1 0 1 1 221 1 3 2 m m m m m y B b B br M B b mh                   3 3 1 1 11 2 1 1 2 1 22 2 1 1 3 f q f q m n q qm m z B br q m f qh                          2 2 2 1 11 1 1 2 1 2 2 2 f q f q f q f q n nq q q qz z f q f q                           (19) 4 4 3 3 12 1 1 1 2 1 22 31 4 3 f q f q f q f q m nq q q qzB b q f q f qh                                 2 22 3 1 11 1 1 2 1 2 3 2 f q f q f q f q n nq q q qz z f q f q                            where 1 1 1 1/ 2 nh z   , 2 1 1 1/ 2 nh z    and /f q m ( f and q are positive integers). In (19), 1N and 1yM are determined by (15). It is obvious that at m=1 the non-linear stress-strain relation (1) transforms into the Hooke’s law, assuming that B=E (here E is the modulus of elasticity). This means that at m=1 Eq. (19) should transform into the formula for curvature of linear-elastic beam. Indeed, at m=1, 0B E and 1 2 0B B  Eq. (19) yields 11 3 1 12 yM Ebh   (20) which is exact match of the formula for curvature of a linear-elastic homogeneous beam of width, b, and height, 1h . Eqs. (18), and (19) should be solved with respect to 1 1nz and 1 by using the MatLab computer program. By substituting of (12) and (16) in (11) the distribution of complementary strain energy density in the lower crack arm is written as     1 122 1 1 111* 1 0 0 1 22 2 4 1 m m n L m z zr zy u B B B mb h             (21) V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 497 Formula (21) can also be applied to obtain the distribution of complimentary strain energy density, *0Uu , in the un-cracked beam portion. For this purpose, r , 1y , 1z , 1 1nz and 1 have to be replaced with 0 , 2y , 2z , 2 2nz and 2 , respectively ( 2 2nz is the neutral axis coordinate of the cross-section of un-cracked beam portion, 2 is the curvature of un-cracked beam portion). The quantities, 2 2nz and 2 , can be determined from Eqs. (18) and (19). For this purpose, r , 1 1nz , 1h and 1 have to be replaced with 0 , 2 2nz , 2h and 2 , respectively. Finally, by substituting of *0Lu , * 0Uu and (8) in (7), one derives         21 1 10 1 21 1 221 1 3 1 1 m m mu u u u u B b B b B brm G b m m m h m                       12 2 1 12 11 2 1 22 2 1 2 1 nm m m mu u u u u u zB br m mh                  3 3 2 2 12 1 1 2 1 22 21 3 2 f q f q f q f qu u u u u u u u nq q q qu u u u u u u u u zB bq f q f qh                                 2 1 1 1 2 f qf q u uu u n qq u u u z f q               (22)       1 1 10 12 1 2 1 1 3 1 m m mu u u u u u B b B bm b m m m                  3 3 2 2 22 2 1 2 1 22 21 3 2 f q f q f q f qu u u u u u u u nq q q qu u u u u u u u u u u u u zB bq f q f qh                                 2 2 2 1 2 f qf q u uu u n qq u u u u u z f q               where 1um m  , /u u uf q m ( uf and uq are positive integers), 1 2 2u nh z   and 2 2 2u nh z    . Formula (22) calculates the strain energy release rate in the beam configuration shown in Fig. 1 when the beam mechanical behavior and the material gradient are described by formulae (1) and (2), respectively. It should be noted that at m=1, 0B E , 1 2 0B B  and 1h h formula (22) yields 2 2 3 21 4 yM G Eb h  (23) which is exact match of the expression for the strain energy release rate when the beam considered is linear-elastic and homogeneous [26]. In order to verify (22), the fracture is analyzed also by using the J-integral written as [27] 0 cos x y u v J u p p ds x x               (24) V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 498 where Γ is a contour of integration going from the lower crack face to the upper crack face in the counter clockwise direction, 0u is the strain energy density, α is the angle between the outwards normal vector to the contour of integration and the crack direction, xp and yp are the components of stress vector, u and v are the components of displacement vector with respect to the crack tip coordinate system xy (x is directed along the crack), ds is a differential element along the contour. The J-integral is solved by using an integration contour, Γ, that coincides with the beam contour (Fig. 1). It is obvious that the J-integral value is non-zero only in segments 1 and 2 of the integration contour. Therefore, the J-integral solution is written as 1 2 J J J   (25) where 1 J and 2J are the J-integral values in segments 1 and 2 , respectively ( 1 and 2 coincide with the free end of lower crack arm and the clamping, respectively). The components of J-integral in segment, 1 , of the integration contour are written as (Fig. 3) xp   mB  , 0yp  (26) 1ds dz , cos 1   (27) where 1z varies in the interval 1 1[ / 2, / 2]h h . Partial derivative, /u x  , in (24) is found as  1 1 11n u z z x        (28) where 1 1nz and 1 are obtained from Eqs. (18) and (19). The strain energy density is calculated by substituting of (1) in (10) 1 0 1 mB u m     (29) By substituting of (1), (12), (16), (26), (27), (28) and (29) in (24) and integrating in boundaries from 1 / 2h to 1 / 2h , one derives         2 21 1 10 1 1 21 1 22 21 4 1 1 1 1 m m mu u u u u B y B B rm J m m b m h m                          12 2 1 12 11 2 1 22 2 1 2 1 nm m m mu u u u u u zB r m mh                  3 3 2 2 12 1 1 2 1 22 21 3 2 f q f q f q f qu u u u u u u u nq q q qu u u u u u u u u zB q f q f qh                                 (30) 2 1 1 1 2 f qf q u uu u n qq u u u z f q               V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 499 where 1y varies in the interval  / 2, / 2b b . The J-integral solution in segment, 2 , of the integration contour (Fig. 1) can be determined by (30). For this purpose, r , 1 1n z , 1 , 2 and 1 have to be replaced with 0 , 2 2nz , 1u , 2u and 2 , respectively. Besides, the sign of (30) must be set to „minus” because the integration contour is directed upwards in segment, 2 . The J-integral final solution is found by substituting of 1 J and 2J in (25)         2 21 1 10 1 1 21 1 22 2 4 1 1 1 1 m m mu u u u u B y B B rm J m m b m h m                       12 2 1 12 11 2 1 22 2 1 2 1 nm m m mu u u u u u zB r m mh                  3 3 2 2 12 1 1 2 1 22 21 3 2 f q f q f q f qu u u u u u u u nq q q qu u u u u u u u u zB q f q f qh                                 2 1 1 1 2 f qf q u uu u n qq u u u z f q               (31)       21 1 10 1 12 1 22 4 1 1 1 m m mu u u u u u B y Bm m m b m                  3 3 2 2 22 2 1 2 1 22 21 3 2 f q f q f q f qu u u u u u u u nq q q qu u u u u u u u u u u u u zB q f q f qh                                 2 2 2 1 2 f qf q u uu u n qq u u u u u z f q               Formula (31) describes the distribution of the J-integral value along the crack front. The average value of the J-integral along the crack front is written as 2 1 2 1 b AV b J J dy b    (32) It should be noted that the J-integral solution derived by substituting of (31) in (32) is exact match of the strain energy release rate (22). This fact verifies the fracture analysis developed in the present paper. RESULTS he distribution of J-integral value along the crack front is analyzed. For this purpose, calculations are performed by using formula (31). It is assumed that b=0.020 m, h=0.004 m, m=0.7, f=7, q=10, 1.7um  , 17uf  , 10uq  and 20yM  Nm. The J-integral value is presented in non-dimensional form by using the formula,  0/NJ J B b . The material gradient along the beam width is characterized by 1 0/B B ratio. Fig. 5 shows the distribution of the J-integral value in non- T V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 500 dimensional form along the crack front for two 1 0/B B ratios at 1 / 2 0.75h h and 2 0/ 2B B  . Only the right-hand half of the crack front is shown in Fig. 5 since the distribution is symmetrical with respect to the crack front centre. The horizontal axis is defined such that 12 / 0.0y b  is in the crack front centre. Thus, 12 / 1.0y b  is in the right-hand lateral surface of the beam. One can observe in Fig. 5 that for 1 0/ 1.2B B  , the J-integral value is maximum in the crack front centre and decreases towards the right-hand lateral surface (this is due to the fact that for 1 0/ 1.2B B  the material property, B, is minimum in the crack front centre and increases towards the beam lateral surfaces). Fig. 5 shows that for 1 0/ 0.4B B   , the J-integral value increases from its minimum in the crack front centre towards the right-hand lateral surface of the beam. Figure 5: The distribution of the J-integral value in non-dimensional form along the crack front (curve 1 - at 1 0/ 1.2B B  , curve 2 – at 1 0/ 0.4B B   ). The effects of the material gradient and the crack location on the strain energy release rate are evaluated too. For this purpose, calculations of the strain energy release rate are carried-out by using formula (22). The strain energy release rate is presented in non-dimensional form by using the formula,  0/NG G B b . Figure 6: The strain energy release rate in non-dimensional form plotted against 1 0/B B ratio at three 1 / 2h h ratios. V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 501 The crack location along the height of the beam cross-section is characterized by 1 / 2h h ratio. In the calculations, 0B is kept constant. Therefore, 1B is varied in order to generate various 1 0/B B ratios. The strain energy release rate in non- dimensional form is plotted against 1 0/B B ratio at three 1 / 2h h ratios and 2 0/ 2.0B B  in Fig. 6. The curves in Fig. 6 indicate that the strain energy release rate decreases with increasing of 1 / 2h h ratio. This finding is attributed to the increase of the lower crack arm stiffness. It can be observed also in Fig. 6 that the strain energy release rate decreases with increasing of 1 0/B B ratio. The influence of non-linear behavior of the material on the fracture is also analyzed. For this purpose, the strain energy release rate in non-dimensional form is presented as a function of 2 0/B B ratio at 1 / 2 0.5h h  and 1 0/ 0.5B B  in Fig. 7. It can be observed that the strain energy release rate decreases with increasing of 2 0/B B ratio (Fig. 7). The strain energy release rate obtained assuming linear-elastic behavior of the two-dimensional functionally graded material is also presented in Fig. 7 for comparison with the non-linear solution (the linear-elastic solution is derived by substituting of m=1 in formula (22)). The curves in Fig. 7 indicate that material non-linearity leads to increase of the strain energy release rate. Figure 7: The strain energy release rate in non-dimensional form presented as a function of 2 0/B B ratio (curve 1 – at non-linear elastic behavior of material, curve 2 – at linear-elastic behavior of material). CONCLUSIONS he longitudinal fracture behavior of two-dimensional functionally graded cantilever beam that exhibits material non-linearity is investigated analytically. The beam is loaded by a bending moment applied at the free end of the lower crack arm. The fracture is studied in terms of the strain energy release rate. The beam mechanical behavior is described by using a power-law stress-strain relation. The material property, B, varies continuously in the beam cross- section according to a quadratic law. The solution derived is applicable for longitudinal crack located arbitrary along the height of the beam cross-section. In order to verify the solution, the fracture is analyzed also by applying the J-integral approach. The distribution of the J- integral value along the crack front is investigated. The physical meanings of the dependence of the J-integral value on material gradient are discussed. The effects of material gradients and crack location along the beam height on the fracture behavior are also analyzed. The analysis revealed that the strain energy release rate decreases with increasing of the lower crack arm thickness. It is found also that the strain energy release rate decreases with increasing of 1 0/B B and 2 0/B B ratios. This finding is attributed to the increase of the beam stiffness. The influence of material non-linearity on the fracture is studied too. It is found that the material non-linearity leads to increase of the strain energy release rate. Therefore, non-linear behavior of the material has to be taken into account in fracture mechanics based safety design of structural members and components made of two-dimensional functionally graded materials. The results obtained in the T V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 502 present study indicate that the longitudinal fracture in two-dimensional functionally graded beams with material non- linearity can be efficiently regulated in their design stage by employing appropriate material gradients. Besides for the beam shown in Fig. 1, the analysis developed in the present paper can be applied to determine the strain energy release rate for longitudinal cracks in two-dimensional functionally graded non-linear elastic beam configurations loaded in pure bending (for instance, the double cantilever beam loaded with uneven bending moments, the four-point bending beam when the crack tip is located in middle beam portion which is loaded in pure bending). It should be noted that the approach developed in the present paper is applicable also for beam configurations loaded by a vertical load when the shear stresses can be neglected, i.e. when the span to height ratio of the beam is large. ACKNOWLEDGMENTS he present study was supported financially by the Research and Design Centre (CNIP) of the UACEG, Sofia (Contract BN – 198/2017). REFERENCES [1] Butcher, R.J., Rousseau, C.E., Tippur, H.V., A functionally graded particulate composite: Measurements and Failure Analysis, Acta Matererialia, 47 (1999) 259-268. [2] Gasik, M.M., Functionally graded materials: bulk processing techniques, International Journal of Materials and Product Technology, 39 (1995) 20-29. [3] Hirai, T., Chen, L., Recent and prospective development of functionally graded materials in Japan, Material Science Forum, 308-311 (1999) 509-514. [4] Mortensen, A., Suresh, S., Functionally graded metals and metal-ceramic composites: Part 1 Processing, International Materials Review, 40 (1995) 239-265. [5] Nemat-Allal, M.M., Ata, M.H., Bayoumi, M.R., Khair-Eldeen, W., Powder metallurgical fabrication and microstructural investigations of Aluminum/Steel functionally graded material, Materials Sciences and Applications, 2 (2011) 1708-1718. [6] Neubrand, A., Rödel, J., Gradient materials: An overview of a novel concept, Zeit f Met, 88 (1997) 358-371. [7] Suresh, S., Mortensen, A., Fundamentals of functionally graded materials, IOM Communications Ltd, London (1998). [8] Tokova, L., Yasinskyy, A., Ma, C.-C, Effect of the layer inhomogeneity on the distribution of stresses and displacements in an elastic multilayer cylinder, Acta Mechanica, (2016) DOI: 10.1007/s00707-015-1519-8, 1 – 13. [9] Tokovyy, Y., Ma, C.-C, Three-dimensional temperature and thermal stress analysis of an inhomogeneous layer, Journal of Thermal Stresses, 36 (2013) 790 – 808, DOI: 10.1080/01495739.2013.787853. [10] Tokovyy, Y., Ma, C.-C, Axisymmetric stresses in an elastic radially inhomogeneous cylinder under length-varying loadings, ASME Journal of Applied Mechanics, 83 (2016), DOI: 10.1115/1.4034459. [11] Uslu Uysal, M., Kremzer, M., Buckling behaviour of short cylindrical functionally gradient polymeric materials, Acta Physica Polonica A, 127 (2015) 1355-1357, DOI:10.12693/APhysPolA.127.1355. [12] Uslu Uysal, M., Buckling behaviours of functionally graded polymeric thin-walled hemispherical shells, Steel and Composite Structures, An International Journal, 21 (2016) 849-862. [13] Szekrenyes, A., Fracture analysis in the modified split-cantilever beam using the classical theories of strength of materials, Journal of Physics: Conference Series, 240 (2010) 012030. [14] Szekrenyes, A., Vicente, W.M., Interlaminar fracture analysis in the GII-GIII plane using prestressed transparent composite beams, Composites Part A: Applied Science and Manufacturing, 43 (2012) 95-103. [15] Szekrenyes, A., Semi-layerwise analysis of laminated plates with nonsingular delamination - The theorem of autocontinuity, Applied Mathematical Modelling, 40 (2016) 1344 – 1371. [16] Bohidar, S.K., Sharma, R., Mishra, P.R., Functionally graded materials: A critical review, International Journal of Research, 1 (2014) 289-301. [17] Erdogan, F., Fracture mechanics of functionally graded materials, Computational Engineering, 5 (1995) 753-770. [18] Paulino, G.C., Fracture in functionally graded materials, Engineering Fracture Mechanics, 69 (2002) 1519-1530. T V. Rizov, Frattura ed Integrità Strutturale, 41 (2017) 491-503; DOI: 10.3221/IGF-ESIS.41.61 503 [19] Shi-Dong Pan, Ji-Cai Feng, Zhen-Gong Zhou, Wu-Lin-Zhi, Four parallel non-symmetric Mode –III cracks with different lengths in a functionally graded material plane, Strength, Fracture and Complexity: an International Journal, 5 (2009) 143-166. [20] Tilbrook, M.T., Moon, R.J., Hoffman, M., Crack propagation in graded composites, Composite Science and Technology, 65 (2005) 201-220. [21] Upadhyay, A.K., Simha, K.R.Y., Equivalent homogeneous variable depth beams for cracked FGM beams; compliance approach, International Journal of Fracture, 144 (2007) 209-213. [22] Uslu Uysal, M., Güven, U., A bonded plate having orthotropic inclusion in adhesive layer under in-plane shear loading, The Journal of Adhesion, 92 (2016) 214-235, DOI:10.1080/00218464.2015.1019064. [23] Petrov, V.V., Non-linear incremental structural mechanics, Infra-Injeneria, M. (2014). [24] Tsukamoto, H., Mechanical properties of zirconia – titanium composites, International Journal of Materials Science and Applications, 3 (2014) 260-267. [25] Lubliner, J., Plasticity theory (Revised edition), University of California, Berkeley, CA (2006). [26] Hutchinson, W., Suo, Z., Mixed mode cracking in layered materials, Advances in Applied Mechanic, 64 (1992) 804– 810. [27] Broek, D., Elementary engineering fracture mechanics, Springer (1986). << /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