Microsoft Word - numero_41_art_4.docx M. Kurek et alii, Frattura ed Integrità Strutturale, 41 (2017) 24-30; DOI: 10.3221/IGF-ESIS.41.04 24 Focused on Multiaxial Fatigue Determination of the critical plane orientation depending on the fatigue curves for bending and torsion Marta Kurek, Tadeusz Łagoda University of Technology - Opole, Poland ma.kurek@po.opole.pl, t.lagoda@po.opole.pl ABSTRACT. This paper contains a proposition setting out a new way of determining an orientation angle of the critical plane. The applied method involves an analysis of fatigue life results’ scatter of several selected materials on the basis of a new criterion of multiaxial fatigue. The expression for the equivalent stress was derived on the basis of pure bending and pure torsion. The paper also contains a verification of the effectiveness of the proposed model by analyzing scatter bands and the value of the ratio of normal to shear stresses. KEYWORDS. Critical plane; Fatigue life; Equivalent stress; Scatter analysys. Citation: Kurek, M., Łagoda, T., Determination of the critical plane orientation depending on the fatigue curves for bending and torsion Frattura ed Integrità Strutturale, 41 (2017) 24-30. Received: 28.02.2017 Accepted: 15.04.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 literature contains a number of propositions including models for fatigue life assessment [1-4]. The work aimed at finding an algorithm for predicting fatigue life accompanying random loads was undertaken by Macha in 1979 [5], who applied the notion of stress for this purpose. Over the successive years, the propositions were improved and modified, and each of them accounted for one or a few aspects of the fatigue life. With the time passing, the strain model was established [6], followed by the one based on energy [7], which applied the parameter of strain energy density parameter acting on the critical plane. Despite the large number of propositions, there is still a need to modify the existing models. The justification for the adopted foundations is based on the scatter of the gained calculation results of the fatigue [8]. They form an important factor used for verifying the effectiveness of the proposed way of assessing fatigue strength. The objective of this paper is to present a new model applicable for assessment of fatigue life on the basis of a concept of the critical plane accompanied by the analysis of the scatter of the results. T M. Kurek et alii, Frattura ed Integrità Strutturale, 41 (2017) 24-30; DOI: 10.3221/IGF-ESIS.41.04 25 THE ALGORITHM OF FATIGUE LIFE EVALUATION he proposed algorithm of predicting fatigue life was formulated on the basis of a concept by Carpintieri [9]. The principal way of modifying the concept by Carpintieri involves a new relation that is established for the angle defining the critical plane orientation, as it plays a critical role in the assessment of fatigue strength. Calculations were done using multiaxial fatigue criterion, based on the critical plane concept. The coefficients B and K in formulas for equivalent stress are calculated according to typical fatigue limits for pure bending and pure torsion. The general form of the equivalent stress [10] according to the criteria on critical plane is expressed as: ( ) ( ) ( )t B t Keq s t     (1) Authors use criterion (1) where weighing factors B and K [11] using pure bending and pure torsion can be defined as: 2 2 2 sin(90 2 ) cos sin 2 sin(90 2 ) cos(90 2 ) 2 cos o o o B B             (2) 2 sin 2 2 2 222 cos B af K B af         (3) 2 af B af    (4) ση(t) is the normal stress and τηs(t) is the shear stress, both acting on the critical plane: 2( ) ( )cos ( )sin 2xx xyt t t      (5) 1 ( ) ( )sin 2 ( )cos 2 2 s xx xyt t t       (6) where     (7) being αη the angle defined by the direction of the maximum normal stress if the above damage accumulation method is applied and β is 2 ( )1 3 1 3 (4 ) 22.5 22.5 ( ) 2 ( ) 2 a fi fi a fi N ctg B N N                        (8) The proposition (8) was derived on the basis of the variability of calculated fatigue strength depending on angle β for several selected materials [11-13]. We can note here that formula (4) is based on the ratio of the fatigue life boundaries. In the proposition (8), this coefficient is relative to the number of cycles, and can be stated the following form ( ) ( )2 ( ) Na fi B N fi Na fi    (9) T M. Kurek et alii, Frattura ed Integrità Strutturale, 41 (2017) 24-30; DOI: 10.3221/IGF-ESIS.41.04 26 A graphical interpretation of formula (8) is presented in Fig. 1, whereas the broken line marks the proposition developed by Carpintieri (10). 22 2 2 3 1 3 1 1 45 ( ) 1 45 2 2 ( ) fi fi B N B B N                             (10) Figure 1: The relationship between the angle β and the parameter B2 according to different models. The final step is the calculation of the fatigue strength. For constant amplitude cyclic loading, the fatigue strength is evaluated by using Basquin’s fatigue characteristics ( A and m ) in compliance with the relevant ASTM standard [14]. The formula for strength calculation under cyclic loading is expressed as follows: ,lg10 eq acal A m N   (11) where ,eq a is the amplitude of the equivalent stress related to the critical plane (Eq.(1)) , A, m - coefficients of regression equation for oscillatory bending. FATIGUE STRENGTH SCATTER ANALYSYS he proposition in Eq. (8) presented in Fig. 1 was developed on the basis of the results of a study involving several materials and requires verification. An analysis of the scatter of the results of fatigue life depending on the proposed angle and value of the ratio of the normal to shear stress was undertaken for this purpose. The analysis applied the values of the scatter bands for the particular materials just as it was calculated in [8] for angle β in the range <0°, 45°>. Fig. 2 presents the methodology followed during the analysis of the proposed expression (8) with regard to the relation of angle β and scatter T. The values of angles β were calculated so that the scatter T is in the range . On this basis, βopt ϵ <βmin, βmax>. Fig. 3 illustrates relationship between scatter values T and the angle β for D30 [15]. T M. Kurek et alii, Frattura ed Integrità Strutturale, 41 (2017) 24-30; DOI: 10.3221/IGF-ESIS.41.04 27 Figure 2: The methodology of presented calculations. Figure 3: Relationship between scatter values T and the angle β for D30. In addition, calculations were performed with regard to the relation of the ratio of normal stress to shear one depending on angle β. The study took into consideration the level of error of 10%. Accordingly, it was adopted that the ratio is established with the probability level of ± 5% 0.95 1.05a a a a a a         . One of the characteristics of fatigue tests is that the results demonstrate a considerable degree of scatter. This paper is based on an assumption that the error is at a level of 10%, which is sufficient during calculations of fatigue life with regard to loads. The results that were gained, were subsequently plotted onto the diagram representing the relation of angle β depending on the ratio of normal to shear stresses (Fig. 4). M. Kurek et alii, Frattura ed Integrità Strutturale, 41 (2017) 24-30; DOI: 10.3221/IGF-ESIS.41.04 28 Figure 4: Exemplary analysis of material D30. The hatched field is the area for which the results of the fatigue life are the same or differ insignificantly. Such analysis was performed for all test materials. Fig. 5 shows the distribution of the scatter results for GTS45 material [16]. Figure 5: Relationship between scatter values T and the angle β for GTS45. The subsequent analysis for exemplary materials is presented in Figs. 6-7. From the analysis of Fig. 7, we can note down that the model also applies to materials, which are not covered by the proposed function (e.g. D30). In the future, it would be necessary to perform a more comprehensive static analysis involving the determination of fatigue characteristics with regard to pure bending and pure torsion for the respective materials. M. Kurek et alii, Frattura ed Integrità Strutturale, 41 (2017) 24-30; DOI: 10.3221/IGF-ESIS.41.04 29 Figure 6: Exemplary analysis of selected materials. Figure 7: Exemplary analysis of selected materials. CONCUSIONS n the basis of the analysis of the experiments as well as calculations and simulations, we can draw the following conclusions: 1. The proposed model has been developed on the basis of relations for materials in elastic-brittle state (Galileo’s hypothesis) and in elastic-plastic state (HMH hypothesis). These states are nevertheless not limiting like in case of the Carpinteri’s model. 2. The proposed model accounts for a wider range of the ratio of normal to shear stresses than the Carpintieri model. 3. The presented analysis is indispensable for the verification of the foundations and effectiveness of the new model. 4. Taking into consideration the scatter of the determined β angle and also the ratio of the fatigue limits for bending and torsion in may be noticed, that most of the materials coincide with the proposed model. 5. In the future more materials should be analyzed in order to continue the verification of the proposed model or to look for an evet better suiting model for experimental tests. O M. Kurek et alii, Frattura ed Integrità Strutturale, 41 (2017) 24-30; DOI: 10.3221/IGF-ESIS.41.04 30 ACKNOWLEDGMENTS roject financed by the National Science Centre. Decision number: 2016/21/D/ST8/02007 REFERENCES [1] Carpinteri, A., Fortese, G., Ronchei, C., Scorza, D., Spagnoli, A., Vantadori S., Fatigue life evaluation of metallic structures under multiaxial random loading, International Journal of Fatigue, 90 (2016) 191-199. DOI: 10.1016/j.ijfatigue.2016.05.007. [2] Skibicki, D., Pejkowski, Ł., Low-cycle multiaxial fatigue behaviour and fatigue life prediction for CuZn37 brass using the stress-strain models, International Journal of Fatigue, 102 (2017) 18–36. [3] Wanga, Y., Susmel, L., The Modified Manson–Coffin Curve Method to estimate fatigue lifetime under complex constant and variable amplitude multiaxial fatigue loading, International Journal of Fatigue, 83 (2016) 135–149. [4] Karolczuk, A., Kluger, K., Łagoda, T., A correction in the algorithm of fatigue life calculation based on the critical plane approach, International Journal of Fatigue, 83 (2016) 174-183. [5] Macha, E., Modele matematyczne trwałości zmęczeniowej materiałów w warunkach losowego złożonego stanu naprężenia, Seria: Monografie, Wrocław, 13 (1979) 100 (in Polish). [6] Macha, E., Generalization of strain criteria of multiaxial cyclic fatigue to random loadings, Fortschr. – Ber. VDI, Reiche 18, Nr 52, VDI-Verlag, Dusseldorf, (1988) 102. [7] Łagoda, T., Macha, E., Generalization of energy multiaxial cyclic fatigue criteria to random loadings, Multiaxial Fatigue and Deformation: Testing and Prediction, ASTM STP 1387, American Society for Testing and Materials, West Conshohocken, PA, (2000) 173. [8] Walat, K., Łagoda, T., Lifetime of semi ductile materials through the critical plane approach, International Journal of Fatigue, 67 (2014) 73-77. [9] Carpinteri, A., Spagnoli, A., Vantadori, S., Multiaxial fatigue assessment using a simplified critical plane-based criterion, International Journal of Fatigue, 33(8) (2011) 969-976. DOI: 10.1016/j.ijfatigue.2011.01.004. [10] Macha, E., Generalization of fatigue fracture criteria for multiaxial sinusoidal loadings in the range of random loadings, in: Biaxial and Multiaxial Fatigue, EGF 3, (Edited by M.W. Brown and K.J. Miller), Mechanical Engineering Publications, London, (1989) 425–436. [11] Carpinteri, A., Kurek, M., Łagoda, T., Vantadori, S., Estimation of fatigue life under multiaxial loading by varying the critical plane orientation, International Journal of Fatigue, (2016). DOI: 10.1016/j.ijfatigue.2016.10.028. [12] Mrzygłód, M., Kurek, M., Łagoda, T., The application of the criteria of multiaxial fatigue in the critical plane for the topology optimization, Application of the fictitious radius for the calculation of the fatigue life using the elastic-plastic body model, AIP Conference Proceedings, 1780, Fatigue Failure and Fracture Mechanics XXVI, (2016) 030003-1-6. [13] Kurek, M., Łagoda, T., Including of ratio of fatigue limits from bending and torsion for estimation fatigue life under cyclic loading, Procedia Materials Science, 12 (2016) 30-35. [14] ASTM E 739–91, Standard practice for statistical analysis of linearized stress–life (S–N) and strain life fatigue data, in: Annual Book of ASTM Standards, Philadelphia, 03.01 (1999) 614–628. [15] Nishihara, T., Kawamoto, M., The Strength of Metals under Combined Alternating Bending and Twisting, Memoirs of the College of Engineering, Kyoto Imperial University, Japan, (1941). [16] Muller, A., Zum Festigkeitsverhalten von mehrachsig stochastisch beanspruchten Gußeisen mit Kugelgraphit und Tempergu. Fraunhofer - Institut fur Betriebsfestigkeit, Darmstadt, (1994). P << /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