Microsoft Word - numero_35_art_46 R.A. Cardoso et alii, Frattura ed Integrità Strutturale, 35 (2016) 405-413; DOI: 10.3221/IGF-ESIS.35.46 405 Focussed on Crack Paths Crack path simulation for cylindrical contact under fretting conditions R.A. Cardoso, J.A. Araújo, J.L.A. Ferreira, F.C. Castro University of Brasília, Brasil rapha_2213@hotmail.com ABSTRACT. In this work different strategies to estimate crack path for cylindrical contacts under fretting conditions are carried out. The main goal is to propose and to evaluate methodologies not only to estimate the direction of crack initiation but also the subsequent propagation in its earlier stages, where the stress field is multiaxial, non-proportional and decays very fast due to the proximity with the contact interface. Such complex conditions pose a substantial challenge to the modelling of crack path. The numerical simulations are provided by a 2D Finite Element Analysis taking into account interactions between the crack faces. The results show that, under fretting conditions, models based on the critical plane method are not effective to estimate the crack initiation orientation, while models based on a so called “critical direction” applied along a critical distance provide better results. Regarding the subsequent crack propagation orientation, it was possible to see that stress intensity factor based models where one considers an infinitesimal virtual crack emerging from an original pre- existent crack are powerful mechanisms of crack orientation estimation. KEYWORDS. Fretting fatigue; Crack propagation; Critical distances. INTRODUCTION retting is the phenomenon that occur in mechanical couplings subjected to contact loads and small relative tangential displacement due to oscillatory loads. In association with remote fatigue loads, the process well known as fretting fatigue take place. Failure due to fretting fatigue usually is observed in engineering assemblies under vibration, such as riveted or bolted connections, dovetail joints in turbines and overhead conductors among, others. The fretting problem is characterized by a strong stress concentration in association with wear, which invariably leads to the nucleation of small cracks. Depending on the level of the stress gradient in such cracks they may arrest, since a threshold value at the crack tip be reached [1]. In order to conduct such analysis criteria able to precise the location, direction of initiation and the further crack path are essential. In these aspects many challenges must be faced, since the fretting problem is characterized for high levels of multiaxility and non-proportionality The aim of this work is to estimate the crack path under fretting conditions in cylindrical contacts assessing some already used methodologies and propose a new methodology to estimate the direction of crack propagation in their earlier stages. Experimental data from literature will be used to confront the theory. CRACK ORIENTATION MODELING s extensively shown in the literature, crack initiation and further propagation is usually characterized by two distinct stages, Fig. 1. The Stage I describes the earlier stage of propagation, commonly governed by Mode II while the Stage II is characterized by a Mode I dominated crack. However, according to [2], Stage I may be F A R.A. Cardoso et alii, Frattura ed Integrità Strutturale, 35 (2016) 405-413; DOI: 10.3221/IGF-ESIS.35.46 406 divided in two groups: type 1 when crack initiation is mainly influenced by the range of the shear stress Δτ, and type 2 when the range of normal stress Δσn governs the process. Figure 1: Scheme of classical stages of crack initiation and short crack propagation. In [3], where tests were conducted in low carbon steel weakened by U-notches to study the high cycle fatigue cracking behaviour, was verified that Stage I, in this kind of stress raiser, is governed by a mixed mode, since the crack profile in Stage I was not parallel to the notch bisector line. It was also noted that the irregular mixed mode dominated path was confined within a distance from the notch tip of the order of L, the material characteristic length, Eq. (1). 2 1 1           thK L (1) where Δσ-1 is the plain fatigue limit and ΔKth is the threshold value of stress intensity factor range, both for a load ratio of -1. Considering these assumptions, two multiaxial criteria based on the theory of critical distances will be applied to estimate the initial crack direction, Fig. 2. Early stage of crack propagation-The critical direction method combined with the TCD (the critical direction method) The theory of critical distances (TCD) [4] can be used in association with any fatigue criterion [5]. So based on the hypothesis that in stage I one may have cracks dominated by the maximum range of normal stress (crack type 2), the crack orientation will be defined here by a line of length 2L, inclined by an angle θσ with respect to the y axis, which maximizes θ in Eq. (2), Fig. 2a. Therefore, an important aspect that one should notice in Eq. (2) is that the maximum range of normal stress is always computed perpendicular to the line defined by angle θ, i.e, the line 2L (inclined by θ) is discretized in many material points and the maximum normal stress range is calculated for the same plane θ for all these points. This procedure provides the critical orientation θσ that a crack of length 2L would have in theory in its stages of initiation and early growth. A similar approach can be used to find the direction that minimizes the shear stress amplitude, Δτ, Eq. (3). The mechanical basis to support this hypothesis is the fact that, in these planes, less energy is wasted with friction and consequently more energy is available for crack propagation [7].         drr L L ),( 2 1 max 2 0   (2)         drr L L r ),( 2 1 min 2 0   (3) In Eqs. 2 and 3 the integrals are carried out over a constant length 2L, as defined by the so-called Line Method (LM). In a similar way one could assume the point or the area method (for 2D problems) [6], but these will not be assessed here due to space constraints. R.A. Cardoso et alii, Frattura ed Integrità Strutturale, 35 (2016) 405-413; DOI: 10.3221/IGF-ESIS.35.46 407 (a) (b) Figure 2: (a) Stress state in the vicinity of the critical zone by means of line method in cylindrical coordinates; (b) critical plane according to the MWCM. Early stage of crack propagation-The critical plane method combined with a multiaxial fatigue criterion applied in terms of the TCD (the critical plane method) In this case, one assumes that crack initiation is governed by the range of shear stress Δτ (crack type 1). Considering a critical plane approach, such as the Modified Wöhler Curve Method (MWCM), crack initiation is expected to take place on the material plane that experiences the highest shear stress amplitude. As the shear stress always has the same magnitude in two orthogonal planes, the critical plane between these two planes will be the one with the maximum normal stress along a cycle. Fig. 2b depicts the critical plane method in association with the point method (PM) to estimate the direction of crack propagation. Once again, in order of brevity only the PM will be assessed here. Crack propagation orientation under fretting conditions In fretting, the load conditions are complex and non-proportional, which makes the process to estimate the directions of propagation hard. Now a review of some relevant criteria already used to estimate the direction of crack propagation under fretting conditions will be presented. Notice that all these criteria are valid for pre-existing cracks (stage II). Classical crack path criteria, such as the maximum tangential stress (MTS) [8], the maximum strain energy density [9] or the maximum energy release rate [10] are not adequate for non-proportional multiaxial loading like in fretting conditions. Taking into account this non-proportionality [2] and [11] considered the following criteria based on [12]: (i) Crack propagates in the direction where k1(θ,t) is maximum during a cycle (ii) Crack propagates in the direction where Δk1(θ) is maximum during a cycle (iii) Crack propagates in the direction where da/dN(θ)is maximum during a cycle where k1 and k2 are the stress intensity factors in mode I and in mode II, respectively, of an infinitesimally small kinked crack emerging from the pre-existent crack with an angle θ, Fig. 3. The expression that relates k1 and k2 with the classical mode I and mode II stress intensities KI and KII is given by Eq. (4), where the angular functions Kij(θ) are reported in [11] and [13]. Figure 3: Inclusion of a small virtual crack emerging from the crack tip of a pre-existing crack. )( )( )()( )()( ),( ),( 2221 1211 2 1 tK tK KK KK tk tk II I      (4) Based on the MTS criterion, Dubourg and Lamacq [2] proposed that the direction of crack propagation may be found searching for the direction that maximizes Δσθ along a cycle, considering that if σθ < 0 then σθ = 0, once that compressive stress does not encourage crack propagation. Recently Giner et al. [7], proposed the criterion of the minimum shear stress range, which consists in finding the plane that minimizes the shear stress range at the crack tip, Fig. 4. As the shear stress R.A. Cardoso et alii, Frattura ed Integrità Strutturale, 35 (2016) 405-413; DOI: 10.3221/IGF-ESIS.35.46 408 always has the same magnitude in two orthogonal planes, the chosen plane for crack propagation direction is the one that experiences the maximum Δσn (less frictional energy is lost). Figure 4: Plane of the minimum shear stress at the crack tip. AVAILABLE DATA o evaluate the accuracy from the aforementioned methodologies experimental data available from tests under fretting conditions will be assessed. All the tests considered are from cylinder-on-plane configuration. Tests were conducted on steel alloy 35NCD16 and AISI 1034 specimens. In both cases the pads were made from different materials, steel 100C6 and chromium 52100 steel, respectively. The mechanical properties are presented on Tab. 1, and the experiments have been reported and discussed in detail in [14] and [15], hence, only essential information necessary to carry out the analyses are briefly presented here. Material E (GPa) v σY (MPa) σu (MPa) σ-1 (MPa) τ-1 (MPa) ΔKth (MPa.m1/2) AISI 1034 200 0.3 350 600 270 - 7 52100 210 0.3 1700 2000 - - - 35NCD16 200 0.3 1127 1270 590 460 - 100C6 195 0.3 1500 - - - - Table 1: Material properties. The load history applied in the tests are as follows. Firstly, a constant normal load, P, was applied to the fretting pad and held constant. Then an oscillatory shear force varying between ±Qa was applied. In the tests analysed here there is absence of a bulk load. All the tests were designed to run in partial slip regime, i.e., Qa/P> /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