Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 56, 4, pp. 915-925, Warsaw 2018 DOI: 10.15632/jtam-pl.56.4.915 INCLINED SURFACE CRACKS IN A GRADED HALF-PLANE SUBJECTED TO FRICTIONAL SLIDING CONTACT Onur Arslan Eskisehir Osmangazi University, Department of Mechanical Engineering, Eskisehir, Turkey e-mail: oarslan@ogu.edu.tr Fracture parameters of an inclined surface crack in a graded half-plane subjected to sli- ding frictional contact are investigated in this study. The problem is modelled via the finite elementmethod (FEM) under the plane strain assumption employing a newly developed di- splacementboundary condition.The shearmodulus of the half-plane is gradedexponentially through the thickness direction bymeans of the homogeneous finite element approach. The augmented Lagrange algorithm is selected as an iterative contact solver. The mixed mode stress intensity factors (SIFs) which are induced by the contact stresses are evaluated utili- zing theDisplacementCorrelationTechnique (DCT). The accuracy of the present procedure is ensured comparing the SIF results to those calculated by an analytical method for verti- cal surface cracks. Additional SIF curves are generated as functions of the crack inclination in order to reveal the effects of non-homogeneity, punch location, crack length and friction coefficient. The prominent conclusion of this study indicates that the crack orientation has profound effects on the behavior of SIF curves, which has not been presented in any study published so far. Hence, consideration of crack inclination is proved essential to successfully predict fracture behavior of a gradedmedium under frictional contact. Keywords: inclined surface cracks, crack/contact problems, gradedmaterials, finite element method 1. Introduction In engineering applications, keeping the structural integrity of a member is the major concern to extend service life. From this point of view, investigation of contact and fracture mechanics problems come into prominence especially for load transfer members. Brittle materials such as ceramics are well suited for load carrying surfaces due to their high wear resistance. However, because of their limited toughness they are prone to surface cracking under frictional contact loadings, which may lead to loss in fatigue strength (Hills et al., 1993). Such destructive effects of contact loadings can be partially tolerated by employing graded materials which consist of spatially varying composition of brittle and ductile materials. For instance, frictional contact and indentation related damages is shown to bemitigated upon introducingmaterial gradation (Surresh, 2001). It is revealed that non-homogeneous ceramics are more damage resistant to loadings induced by spherical intenders (Lawn, 1998). Spatial gradation of Young’s modulus is proven to prevent formation of Hertzian cone-cracks on brittle surfaces loaded by spheri- cal stamps (Jitcharoen et al., 1998). Material gradation provides an improvement in damage tolerance of ceramics which are subjected to contact and flexural loadings (Zhang, 2012). In literature, there are some computational and analytical studies analyzing behavior of cracks in elastic continua which are subjected to contact loadings. Hasebe et al. (1989) investi- gated edge crack problems in an elastic half plane under frictional contact by utilizing a rational mapping function and complex stress functions. Fracture parameters of a surface crack in a 916 O. Arslan gradedmediumwere evaluated under the action of frictional contact loadings by using the sin- gular integral equation approach (Dag, 2001; Dag and Erdogan, 2002). Embedded horizontal cracks in isotropic (El-Borgi et al., 2004) and orthotropic (El-Borgi et al., 2013) graded coatings subjected to frictional Hertzian contact assumption were examined by means of the singular integral equation approach. An edge crack in a graded coating loaded by a flat rigid stamp was investigated by means of an uncoupled solution procedure based on the singular integral equation approach (Dag et al., 2012). Sarikaya and Dag (2016) analytically calculated fracture parameters of a surface crack in an orthotropic half-plane whose surface was exposed to sliding frictional contact loadings. In all the aforementioned studies, the cracks are located parallel or perpendicular to the contact surface. However, propagation of the cracks is observed to be slanted in an elastic medium which is subjected to fretting contact (Faanes, 1995; Giner et al., 2014). Moreover, surface cracks may possibly initiate and propagate in any direction through a medium which is under non-symmetrical loading conditions, such as frictional contact loadings. Hence, effects of crack orientation on fracture parameters should be examined to clearly predict critical and sub-critical crack growth. The literature survey revealed that no computational or analytical studies related to inclined surface cracks in a gradedmedium under frictional contact has been conducted. In this study, fracture parameters of an inclined surface crack in an elastic graded medium which is subjected to frictional sliding contact are investigated. Themedium ismodelled on the plane strain assumption bymeans of a finite element analysis software ANSYS (2013). A newly developed displacement boundary condition is imposed on the model to be able to improve the simulation capability. The shear modulus of the half-plane is graded exponentially through the thickness direction utilizing the homogeneous finite element method. In the solution of the nonlinear contact problem, the augmented Lagrange algorithm is selected. The crack tip SIFs induced from the contact stresses are calculated by using DCT. Comparisons of the results to those generated by the analytical method for vertical surface cracks (Dag, 2001) validate the accuracy of the proposed finite element procedure. Extra SIF curves are depicted as functions of crack inclination for different non-homogeneity parameters, punch location, relative crack length and friction coefficient. The striking outcome of this study indicates that the change in crack inclination have remarkable influence on the mixed mode SIFs. As a result, consideration of crack inclination in crack-contact problems is proved essential to successfully predict fracture behavior of a gradedmedium under frictional contact loadings. 2. Problem definition and solution method The geometry of an inclined surface crack in a graded medium which is subjected to frictional sliding contact is depicted in Fig. 1. The surface of the gradedmedium is loaded by a rigid flat punch whose end points are located at y = b and y = c. Coulomb’s dry friction is assumed to exist between the rigid punch and themedium surface. Exponential shearmodulus gradation is introduced through the thickness direction considering Poisson’s ratio as a constant (Dag, 2001; Dag and Erdogan, 2002) µ(x)=µ0exp(γx) (2.1) where µ0 and γ represent the shear modulus at the surface (x = 0) and the non-homogeneity parameter, respectively. Thus, constitutive relations for the gradedmedium can be stated as ε=exp(−γx)Sσ (2.2) Inclined surface cracks in a graded half-plane... 917 Fig. 1. Problem configuration where S, ε and σ respectively stand for the compliance matrix, strain vector and stress vector in the global coordinate system xy. The stress and strain vectors can be written explicitly as ε= [εxx,εyy,εzz,εyz,εxz,εxy] T σ= [σxx,σyy,σzz,σyz,σxz,σxy] T (2.3) Isotropic materials have an infinite number of elastic symmetry planes and possess two independentmaterial properties. The compliance matrix in Eq. (2.2) can be expressed in terms of the shear modulus µ0 and Poisson’s ratio ν as follows (Jones, 1998) S=          C1 C2 C2 0 0 0 C2 C1 C2 0 0 0 C2 C2 C1 0 0 0 0 0 0 C1−C2 0 0 0 0 0 0 C1−C2 0 0 0 0 0 0 C1−C2          (2.4) where C1 = 1 2(1+ν)µ0 C2 =− v 2(1+ν)µ0 Constitutive relations can be restated for plane elasticity in the local coordinate system x1x2 which is defined at the crack tip      ε11 ε22 2ε12      =exp[−γ(x1+a)cosφ]    s11 s12 s16 s21 s22 s26 s61 s62 s66         σ11 σ22 σ12      (2.5) where a, φ and sij (i,j = 1,2,6) represent length of the surface crack, crack inclination angle and compliance coefficients, respectively. The compliance coefficients for the plane stress case read s11 = s22 =C1 s12 =C2 s16 = s26 =0 s66 =2(C1−C2) (2.6) The compliance coefficients for the plane strain case are defined below s11 = s22 = C21 −C 2 2 C1 s12 = C2(1−C2) C1 s16 = s26 =0 s66 =2(C1−C2) (2.7) 918 O. Arslan Analytical expressions of the openingmode (mode I) and the shearingmode (mode II) SIFs are given as follows (Perez, 2004) KI = lim x1→0+ √ 2π(x1−a)σ22(x1,0) KII = lim x1→0+ √ 2π(x1−a)σ12(x1,0) (2.8) In order to calculate SIFs through Eq. (2.8), one should determine the asymptotic stress expressions by conducting a heavy solution procedure based on the singular integral equation approach (Dag, 2001; Dag and Erdogan, 2002). On the other hand, SIFs can be calculated via DCT with a high level of accuracy within a properly constructed finite element model. DCT utilizes nodal displacements around crack tip which are computed through FEM (Kim and Paulino, 2003) KI = 1 8 √ 2π ∆a 4V1−V2 s22exp[−γacosφ] KII = 1 8 √ 2π ∆a 4U1−U2 s11exp[−γacosφ] (2.9) where U1 =u1(N1)−u1(N3) U2 =u1(N2)−u1(N4) V1 =u2(N1)−u2(N3) V2 =u2(N2)−u2(N4) (2.10) where Nj (j = 1, . . . ,4) are the nodes around the crack tip, which are depicted in Fig. 2. ∆a represents the characteristic length of the crack tip elements. u1 andu2 are the displacement components defined through x1 and x2 directions, respectively. Fig. 2. Nodes on crack tip elements,∆a/a=1/40000 The finite element model (as seen in Fig. 3) possesses boundaries on the contrary to the original half-plane problem.B,H andW represent width of the flat punch, height of themodel and width of the model, respectively. The ratios B/H and B/W are selected sufficiently small that they do not have any effect on the contact stresses and crack tip parameters. When the side-lines (lines 1 and 2) of the model are set free, they experience some displa- cements under the effect of frictional contact stresses leading to incorrect crack tip parameters. These adverse displacement effects are observed to be still active, even if the model is chosen to be very large compared to the crack/contact region. Hence, the effect of side-line displace- ments on the crack-contact region should be neutralized in order to satisfactorily simulate the half-plane problem. If line 1 and line 2 are forced to undergo the same horizontal displacements under the effect of contact stresses, the behavior of the surface crack can be retrieved from the side-line effects to a certain degree. Inclined surface cracks in a graded half-plane... 919 Fig. 3. Finite element model,B/H =1/30,B/W =1/80, φ=−30◦ The following displacement boundary condition is imposed on the side-lines of the model which increases the simulation capability of the model v ( x,− W 2 ) = v ( x, W 2 ) 0 0. In addition, normalizations of mode I (KI) andmode II (KII) SIFs at the crack tip are defined as KIn = √ a π KI P KIIn = √ a π KII P (3.1) Note that negativemode I (KI) SIF results are not valid, and the negativity reveals that the crack closure occurs. However, these negative results are still useful and applicable in considera- tion of some problems including remote loading (Dag and Erdogan, 2002). Hence, in literature, negative KI results have been presented in parametric analysis of crack/contact problem solu- tions (Dag, 2001; Dag and Erdogan, 2002; Sarikaya and Dag, 2016). Figure 6 reveals comparisons of the results to those evaluated analytically by Dag (2001) for vertical surface cracks (φ= 0◦). The comparisons are depicted for the combinations of two different normalized non-homogeneity paremeters γa and four different friction coefficients η. As can be observed in this figure, excellent agreement between the analytical (Dag, 2001) and present results is attained in each case, which indicates the validity and reliability of the present finite element solution. Fig. 6. Comparisons to the analytical results obtained by Dag (2001) for a graded isotropic half plane for ν =0.25, a/(c− b)=1.0, φ=0◦; (a)KIn, γa=0.1, (b)KIn, γa=1.0, (c)KIIn, γa=0.1, (d) KIIn, γa=1.0 922 O. Arslan Figure 7 illustrates deformed shape of the graded medium for a/(c− b) = 1.0, b/a = 1.0, φ=−30◦, η = 0.4, ν = 0.25, and γa= 1.0. Mesh densities on the crack-contact region can be seen in this figure, and the surface crack is observed to be opened under the action of frictional contact loading. Figure 8 reveals the effect of normalized non-homogeneity parameter γa on the curves of normalizedmixed-mode SIFs versus the inclination angle φ. Examining Fig. 8a, when the normalized non-homogeneity parameter γa is increased from 0.1 to 1.0, KIn decreases for φ < 0◦, however increases for φ > 0◦. Note that the crack closure (KIn = 0) occurs at higher inclination angles φ, as γa is increased from 0.1 to 1.0. Observing Fig. 8b, when γa is increased from 0.1 to 1.0,KIIn has an decreasing trend in the range−20 ◦ <φ< 60◦. However, this trend is reversed in the range −60◦ <φ<−20◦. Also note that the absolute value of KIIn is critical in design, not its sign. Fig. 7. Deformed shape of the crack-contact region, ν =0.25, a/(c− b)=1.0, b/a=1.0, φ=−30◦, η=0.4,γa=1.0 Fig. 8. Effect of normalized non-homogeneity parameter on normalizedmode I (a) andmode II (b) SIFs, ν=0.25, a/(c− b)= 1.0, b/a=1.0, η=0.4 Figure 9 demonstrates the plots of normalized mixed-mode SIFs at the crack tip versus the inclination angle φ for different coefficients of friction η. When the coefficient of friction is increased from 0.0 to 0.8, KIn elevates and KIIn decreases significantly. Note that the crack seems closed for η< 0.2. InFig. 10, the normalizedmixed-modeSIFs at the crack tip are plotted versus the inclination angle φ for different relative crack lengths a/(c− b). As seen in Fig. 10a, KIn increases remarkably when a/(c− b) is decreased from 5.0 to 0.66̄ for φ> 0. When φ< 0, KIn is not considerably affected by the change in a/(c− b). Moreover, the crack closure seems to occur at higher inclanation angles φ, as a/(c−b) is decreased from 5.0 to 0.66̄. As displayed Inclined surface cracks in a graded half-plane... 923 in Fig. 10b, when a/(c−b) is decreased from 5.0 to 0.66̄,KIIn decreases for φ> 0 and elevates for φ< 0. Fig. 9. Effect of friction on normalizedmode I (a) andmode II (b) SIFs, ν=0.25, a/(c− b)= 1.0, b/a=1.0,γa=1.0 Fig. 10. Effect of relative crack length a/(c− b) on normalizedmode I (a) andmode II (b) SIFs, ν=0.25, b/a=1.0, γa=1.0, η=0.4 Fig. 11. Variations of normalizedmode I (a) andmode II (b) SIFs versus the normalized punch location b/a for different crack inclinations, ν =0.25, a/(c− b)=1.0, γa=1.0, η=0.4 Figure 11 illustrates the plots of normalized mixed-mode SIFs at the crack tip versus the normalized punch location b/a for different crack inclanation angles φ. As seen in Fig. 11a, KIn reaches the local extremum when b/a is between 1.0 and 2.0 for 0¬ φ¬ 60. In this range 924 O. Arslan of b/a, magnitude of KIn elevates significantly when the inclination angle φ is increased from −60◦ to 60◦. It can be also infered fromFig. 11a that the crack is substantially closed (KIn < 0) for φ­ 30◦ in the range 0 0), the risk of mode I type failure may arise as the degree of gradation is increased. The risk of mode I failure seems existing especially for η > 0.2. For the surface cracks having negative inclination (φ < 0), the risks of mode II type failure can be alleviated substantially by increasing the contact length. However, for the surface cracks having positive inclination (φ> 0), the risk ofmode I type crack propagation can bemitigated upon decreasing the contact length. To conclude: effects of the problemparameters on the crack tip SIFs are shown to be altered significantly with the variation in crack orientation φ, which is the main outcome of this study. Hence, investigation of the crack tipparameters only for vertical cracks is shown tobe insufficient and, therefore, surface crack problems should be examined by taking the crack orientation into account to be able to clearly predict the critical and sub-critical surface crack growth in graded structures subjected to frictional contact loadings. A new displacement boundary condition is also proposed to improve the computational modelling of crack/contact problems. References 1. ANSYS Inc., 2013, Ansys Basic Analysis Procedures Guide – Release 15, Ansys Inc., Canonsburg 2. Arslan O., 2017, Computational contactmechanics analysis of laterally graded orthotropic half- planes,World Journal of Engineering, 14, 145-154 3. Dag S., 2001, Crack and contact problems in gradedmaterials, Ph.D. Thesis, Lehigh University, Bethlehem 4. Dag S., Apatay T., Guler M.A., Gulgec M., 2012, A surface crack in a graded coating subjected to sliding frictional contact,Engineering Fracture Mechanics, 80, 72-91 5. Dag S., Erdogan F., 2002, A surface crack in a gradedmedium loaded by a sliding rigid stamp, Engineering Fracture Mechanics, 69, 1729-1751 6. Dag S., Guler M.A., Yildirim B., Ozatag A.C., 2009, Sliding frictional contact between a rigid punch and a laterally graded medium, International Journal of Solids and Structures, 46, 4038-4053 Inclined surface cracks in a graded half-plane... 925 7. Dag S., Yildirim B. and Sarikaya D., 2007, Mixed-mode fracture analysis of orthotropic functionally gradedmaterials undermechanical and thermal loads, International Journal of Solids and Structures, 44, 7816-40 8. El-Borgi S., Keer L.,WissemB.S., 2004,An embedded crack in a functionally graded coating bonded to a homogeneous substrate under frictional Hertzian contact,Wear, 257, 760-776 9. El-Borgi S., Keer L., Wissem B.S., 2013, An embedded crack in a functionally graded ortho- tropic coating bonded to a homogeneous substrate under frictional Hertzian contact, International Journal of Solids and Structures, 50, 3898-3910 10. Faanes S., 1995, Inclined cracks in fretting fatigue,Engineering FractureMechanics, 80, 6, 71-82 11. Giner E., Sabsabi M., Rodenas J.J., Fuenmayor F.J., 2014, Direction of crack propagation in a complete contact fretting-fatigue problem, International Journal of Fatigue, 58, 172-180 12. Hasebe N., Okumura M., Nakamura T., 1989, Frictional punch and crack in plane elasticity, Journal of Engineering Mechanics, 115, 6, 1137-1149 13. Hills D.A., Nowell D., Sackfield A., 1993, Mechanics of Elastic Contacts, Butterworth- -Heinemann, London 14. Jitcharoen J., Padture N.P., Giannakopoulos A.E., Suresh S., 1998, Hertzian-Crack suppression in ceramics with elastic-modulus-graded surfaces, Journal of the American Ceramic Society, 81, 9, 2301-2308 15. Jones R.M., 1998,Mechanics of Composite Materials, Taylor & Francis, Philadelphia 16. Kim J-H., Paulino G.H., 2003, Mixed-mode J-integral formulation and implementation using graded elements for fracture analysis of nonhomogeneous orthotropicmaterials,Mechanics of Ma- terials, 35, 107-128 17. Lawn B., 1998, Indentation of ceramics with spheres: a century after Hertz, Journal of the Ame- rican Ceramic Society, 81, 8, 1977-1994 18. Mijar A.R., Arora J.S., 2004, An augmented Lagrangian optimization method for contact analysis problems, 1: formulation and algorithm, Structural Multidisciplinary Optimization, 28, 99-112 19. Perez N., 2004,Fracture Mechanics, Kluwer Academic Publishers, Boston 20. Santare M.H., Lambros J., 2000), Use of graded finite elements to model behavior of nonho- mogeneousmaterials, Journal of Applied Mechanics, 67, 819-822 21. Sarikaya D., Dag S., 2016, Surface cracking in an orthotropic medium subjected to frictional contact, International Journal of Solids and Structures, 90, 1-11 22. Surresh S., 2001, Graded materials for resistance to contact deformation and damage, Science, 292, 5526, 2447-2451 23. Zhang Y., 2012, Overview: damage resistance of graded ceramic restorativematerials, Journal of the European Ceramic Society, 32, 11, 2623-2632 Manuscript received September 17, 2017; accepted for print February 1, 2018