Microsoft Word - numero_38_art_24 A. Niesłony, Frattura ed Integrità Strutturale, 38 (2016) 177-183; DOI: 10.3221/IGF-ESIS.38.24 177 Focussed on Multiaxial Fatigue and Fracture A critical analysis of the Mises stress criterion used in frequency domain fatigue life prediction Adam Niesłony Opole University of Technology, Poland a.nieslony@po.opole.pl, http://orcid.org/0000-0002-1218-8341 ABSTRACT. Multiaxial fatigue failure criteria are formulated in time and frequency domain. The number of frequency domain criteria is rather small and the most popular one is the equivalent von Mises stress criterion. This criterion was elaborated by Preumont and Piefort on the basis of well-known von Mises stress concept, first proposed by Huber in 1907, and well accepted by the scientific community and engineers. It is important to know, that the criterion was developed to determine the yield stress and material effort under static load. Therefore the direct use of equivalent von Mises stress criterion for fatigue life prediction can lead to some incorrectness of theoretical and practical nature. In the present study four aspects were discussed: influence of the value of fatigue strength of tension and torsion, lack of parallelism of the SN curves, abnormal behaviour of the criterion under biaxial tension- compression and influence of phase shift between particular stress state components. Information contained in this article will help to prevent improper use of this criterion and contributes to its better understanding. KEYWORDS. Mises stress; Frequency domain; Multiaxial fatigue. Citation: Niesłony, A., A critical analysis of the Mises stress criterion used in frequency domain fatigue life prediction, Frattura ed Integrità Strutturale, 38 (2016) 177-183. Received: 05.05.2016 Accepted: 10.06.2016 Published: 01.10.2016 Copyright: © 2016 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 ultiaxial fatigue life criteria are used during fatigue life prediction to compute equivalent uniaxial strain or stress state. This procedure is essential for the calculation of final fatigue life based on uniaxial fatigue characteristics. There are many propositions of multiaxial fatigue failure criteria in the literature based on stress invariants, critical plane concept, integral criteria and other [1–3]. Most of them are defined in time domain and implemented for working on stress or strain time histories or on sets of amplitudes. There are only few propositions dedicated for frequency based fatigue life assessment were the frequency domain definition of the criterion is expected [4, 5]. One of them is the criterion of Equivalent von Mises Stress (EMS), which is nowadays more and more popular and well accepted by engineers [6]. In this paper advantages and possible risks of usage of this criterion are taken under discussion. Particular four aspects were presented in details which are causing significant effect on computed fatigue life:  differences between fraction of fatigue strength of pure tension and torsion and the value of square root of 3,  lack of parallelism of the S-N curves for pure tension and torsion, M A. Niesłony, Frattura ed Integrità Strutturale, 38 (2016) 177-183; DOI: 10.3221/IGF-ESIS.38.24 178  abnormal behaviour of the criterion under biaxial tension-compression loading condition in comparison to the experimental results presented in the literature,  influence of phase shift between particular stress state components on the value of equivalent stress. The aim of this paper is the clear presentation of the theoretical limitation and the area of practical application of the EMS criterion. Examples based on technical experiments are recalled which confirms the correctness of discussed limitations. Provided information are highly important for engineers which are using equivalent von Mises stress in frequency domain as well as for researchers who are working on new multiaxial fatigue failure criteria where von Mises stress is using as a part of the definition. EQUIVALENT VON MISES STRESS CRITERION IN FREQUENCY DOMAIN he main reason for creating of this document are increasingly frequent, uncritical application of the Equivalent von Mises Stress (EMS) criterion in engineering calculations in the field of fatigue assessment. Since M.T. Huber in 1904 [7] and R. von Mises in 1913 [8] publish the theoretical background of the possibility of measurement of material effort by specific work of strain EMS criterion has become the most widely used in industry and science, among others in determining equivalent uniaxial stress, yield strength limit or other parameters used for example in constitutive equations [9]. Without going into details about derivation, which are well presented in many publications, the final equation for equivalent stress EMS criterion can be written as follow:        EMS xx yy yy zz zz xx xy yz zx 2 2 2 2 2 21 6 2                      (1) or        EMS xx xx xx xx xx yy zz zz xx xy yz zx 2 2 2 2 2 2 23 3 3                     (2) using proper components of the stress tensor: xx xy xz yx yy yz zx zy zz                     σ (3) In 1994 Preumont and Piefort [10] represent a breakthrough adapt of the von Mises criterion for determining the material fatigue directly in frequency domain. They propose to calculate of the power spectral density of equivalent von Mises stress directly in frequency domain as follow  EMS MG f f( ) Trace ( ) Q G (4) However in the referenced paper [10] plane stress state was analysed presented method is well applicable also in spatial stress state with following power spectral density matrix:           xx xx xx yz yz xx yz yz G f G f f G f G f , , , ,            G      (5) defined according to vector of stress tensor components xx yy zz xy yz zx[ ]     S (6) and with von Mises coefficient matrix [11] T A. Niesłony, Frattura ed Integrità Strutturale, 38 (2016) 177-183; DOI: 10.3221/IGF-ESIS.38.24 179 M 1 0.5 0.5 0 0 0 0.5 1 0.5 0 0 0 0.5 0.5 1 0 0 0 0 0 0 3 0 0 0 0 0 0 3 0 0 0 0 0 0 3                       Q (7) After solving Eq. (4) using (5) and (7) following expression for PSD of EMS can be written       EMS xx xx yy yy zz zz xy xy yz yz zx zx xx yy yy zz zz xx G f G f G f G f G f G f G f G f G f G f , , , , , , , , , ( ) ( ) ( ) ( ) 3 ( ) 3 ( ) 3 ( ) Re ( ) Re ( ) Re ( )           (8) PSD of equivalent von Mises stress in the form of Eq. (8) can be useful while programing this criterion in low level programing languages where matrix operation (4) cannot be easily realised. LIMITATIONS ON THE USE OF EMS CRITERION IN THE FATIGUE CALCULATION Fraction of Fatigue Strengths of Tension and Torsion Let us set the loading of an abstract structure and the reference axes in such a way as to obtain only one a non-zero shear stress component tor xy. [ 0 0 0 0 0 ]S (9) In practice this can be a stress state observed at surface of round specimen under pure torsion. In such a case PSD matrix will possess only one nonzero component – the component Gxy,xy( f ). According to the Eq. (8) equivalent PSD function of stress reduce to the following form EMS xy xyG f G f,( ) 3 ( ) (10) Generally speaking PSD function of stress describes how power of a stress history is distributed over frequency. ‘The power’ should be understood as the variance of the stress history, what can be expressed as follow G f df 0 ( )    (11) where  is the variance of the stress history. Also, it is possible to calculate the expected signal amplitude for a specified small frequency range of a width of f fr f a fr G f df2 2 ( )       (12) Analysing the Eq. (10) it can be seen that according to the von Mises stress criterion for computing PSD of equivalent uniaxial stress is to multiply PSD for pure torsion times 3. Transforming this action to the stress amplitude following the relationship presented in Eq. (12) we get the equation: A. Niesłony, Frattura ed Integrità Strutturale, 38 (2016) 177-183; DOI: 10.3221/IGF-ESIS.38.24 180 af af af af 3 3 1.7321         (13) what is also valid for basic von Mises criterion in time domain, Eqs. (1) and (2). According to numerous reported experimental results such an equality is fulfilled only for few materials. Usually the ratio (13) varies between 1 and 2 for fatigue strength and fatigue limit as well. It is important to know how much influence have a deviation from the specified square root of three value (13) on calculated fatigue life. In order to present the scale of the problem fatigue life was calculated for round specimen under random, narrow-banded and Gaussian, pure torsion loading. In such a kind of loading the Probability Density Function (PDF) of amplitudes describe Rayleigh distribution [12]. On the Fig. 1a) PDF for shear stress amplitudes and equivalent tension amplitudes according EMS criterion were presented. It was also assumed that the Miner rule is applicable and constant amplitude SN curves for pure torsion and tension are known and described as follow m m a f a f N N S T N N 1 1 ,                      (14) where: Sf and Tf – fatigue limits for tension and torsion; N and N – number of cycles for knee points; m and m – slopes of the SN curves. On the Fig. 1b) two SN curve are presented which satisfy the Sf / Tf = 1.7321 condition, Eq. (13). Computed fatigue life T1 according PDF of shear stress amplitudes and SN curve for torsion are equal to T2 computed from PDF of equivalent tension amplitudes and SN curve for tension. For materials that do not meet the condition (13) computed fatigue life T1 and T2 differ significantly. Such a case is presented on Fig. 1c) for Sf / Tf = 1.5 what results in T1 / T2 = 3.16. Lack of Parallelism of the SN Curves Lack of parallelism is a special situation of the problem discussed in previous section. In this case the equality (13) cannot be fulfilled in whole range of cycles to failure. Depends on the m and m slopes of SN curves different deviation from T1 / T2 = 1 can be obtained. This effect was illustrated in Figs. 1d) 1e) 1f) where computed results obtained with the same procedure as described in previous section were presented. Abnormal Behaviour of the EMS Criterion under Biaxial Tension-Compression Biaxial tension-compression is a plane stress state where for specific reference axes only the shear component is constant and equal zero bi xx yy. [ 0 0 0 0 ] S (15) Pure biaxial tension-compression is rarely found in practice but it is used for verification of multiaxial fatigue failure criteria. Such kind of stress state is released on so called cruciform specimens through loading of two sets of arms in perpendicular direction [13]. This two loading components can be of any type, for example in-phase, out-of-phase or random with given correlation coefficient. In biaxial tension-compression the PSD of equivalent stress can be computed as follow (16) Real part of cross spectral density Re[Gxx,yy(f)] is equal 0 for fully uncorrelated loading components xx(t) and yy(t). There are some interesting results published in the literature which are showing basic biaxial fatigue behaviours of tested materials. Cláudio et al. [13] presenting results from which it appears that the EMS criterion does not fit the real behaviour of material. This can be observed on the Fig. 2 where EMS criterion for correlated data gives lower stress amplitudes and for out-of- phase loading higher stress amplitudes than expected. EMS xx xx yy yy xx yyG f G f G f G f, , ,( ) ( ) ( ) Re[ ( )]   A. Niesłony, Frattura ed Integrità Strutturale, 38 (2016) 177-183; DOI: 10.3221/IGF-ESIS.38.24 181 a) b) c) d) e) f) Figure 1: PDF of amplitudes used for computation of fatigue life T1 and T2 (a) and five sets of SN curves for torsion and uniaxial tension (b), (c), (d), (e), and (f). Influence of Phase Shift between Particular Stress State Components It is well known under fatigue community that the phase shift between particular components of the multiaxial loading are influencing the fatigue life. A special interest among scientists and researchers has a combination of tension-compression or bending and torsion, as it is commonly encountered in a responsible machine elements such as shafts. Many test results have been published in this area. All of them show a significant effect of the phase shift fatigue. Analysed criterion does not have such properties. Under a combination of tension and torsion (17) 0 100 200 300 400 500 600 700 800 0 1 2 3 4 5 6 7 x 10 −3 σ a , τ a , MPa p σ a (σ a ), p τ a (τ a ), M P a − 1 p σa (σ a ) pτa (τa) 10 2 10 4 10 6 10 8 10 2 10 3 N, cycles σ a , τ a , M P a T1 T2 = 1.0 Tension: S f = √ 3 · 200; m = 8; N0 = 2 · 106 Torsion: T f = 200; m = 8; N0 = 2 · 106 10 2 10 4 10 6 10 8 10 2 10 3 N, cycles σ a , τ a , M P a T1 T2 = 3.16 Tension: Sf = 300; m = 8; N0 = 2 · 106 Torsion: T f = 200; m = 8; N0 = 2 · 106 10 2 10 4 10 6 10 8 10 2 10 3 N, cycles σ a , τ a , M P a T1 T2 = 0.133 Tension: S f = √ 3 · 200; m = 8; N0 = 2 · 106 Torsion: T f = 200; m = 12; N0 = 2 · 106 10 2 10 4 10 6 10 8 10 2 10 3 N, cycles σ a , τ a , M P a T1 T2 = 0.421 Tension: Sf = 300; m = 8; N0 = 2 · 106 Torsion: T f = 200; m = 12; N0 = 2 · 106 10 2 10 4 10 6 10 8 10 2 10 3 N, cycles σ a , τ a , M P a T1 T2 = 0.0422 Tension: Sf = 400; m = 8; N0 = 2 · 106 Torsion: T f = 200; m = 12; N0 = 2 · 106 tt xx xy. [ 0 0 0 0 ] S A. Niesłony, Frattura ed Integrità Strutturale, 38 (2016) 177-183; DOI: 10.3221/IGF-ESIS.38.24 182 the Eq. (8) for PSD of EMS simplify to following expression (18) As we can see cross power spectrum is not present. This function is the only one which includes information about correlation (phase shift for one harmonic component) what can be treat as a proof of the omission of phase shift effect. Figure 2: SN curves uniaxial and biaxial data presented by Cláudio et al. [13]. REMARKS AND FINAL CONCLUSIONS 1. It is not recommended to use the EMS criterion in a case where the shear stresses dominate, and the ratio of fatigue limits is different from the square root of three. 2. Abnormal behaviour of EMS criterion it can be observed in comparison to experimental results under biaxial tension- compression. For example uniaxial and in-phase biaxial loading give the same equivalent amplitude according this criterion but the test results showing shortening of the fatigue life, see Fig. 2. 3. The impact of non-parallelism of fatigue characteristics on calculated life is significant and depends on loading level. 4. Correlation between normal and shear stress components are neglected. Therefore, this criterion can be used for materials that do not show sensitivity to the phase shift between these components. REFERENCES [1] Karolczuk, A., Macha, E., A review of critical plane orientations in multiaxial fatigue failure criteria of metallic materials, Int J Fract., 134 (2005) 267–304. doi:10.1007/s10704-005-1088-2. [2] Wang, Y., Susmel, L., The Modified Manson–Coffin Curve Method to estimate fatigue lifetime under complex constant and variable amplitude multiaxial fatigue loading, Int. J. Fatigue., 83 (2016) 135–149. doi:10.1016/j.ijfatigue.2015.10.005. [3] Ince, A., Glinka, G., A generalized fatigue damage parameter for multiaxial fatigue life prediction under proportional and non-proportional loadings, Int. J. Fatigue. 62 (2014) 34–41. doi:10.1016/j.ijfatigue.2013.10.007. [4] Pitoiset, X., Preumont, A., Spectral methods for multiaxial random fatigue analysis of metallic structures, Int. J. Fatigue. 22 (2000) 541–550. doi:10.1016/S0142-1123(00)00038-4. [5] Niesłony, A., Comparison of some selected multiaxial fatigue failure criteria dedicated for spectral method, J. Theor. Appl. Mech., 48 (2010) 233–254. [6] de la Fuente, E., An efficient procedure to obtain exact solutions in random vibration analysis of linear structures, Engineering Structures, 30 (2008) 2981–2990. doi:10.1016/j.engstruct.2008.04.015. [7] Huber, M.T., Specific work of strain as a measure of material effort, Archives of Mechanics, 56 (2004) 173–190. [8] v Mises, R., Mechanik der festen Körper im plastisch- deformablen Zustand, Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 1913 (1913) 582–592. EMS xx xx xy xyG f G f G f, ,( ) ( ) 3 ( )  10 4 10 5 10 6 10 7 10 8 10 2 N, cycles σ x x , a , σ E M S , a , M P a uniaxial: data uniaxial: fit biaxial δ = 0◦: data biaxial δ = 0◦: fit biaxial δ = 90◦: data biaxial δ = 90◦: fit biaxial δ = 180◦: data biaxial δ = 180◦: fit A. Niesłony, Frattura ed Integrità Strutturale, 38 (2016) 177-183; DOI: 10.3221/IGF-ESIS.38.24 183 [9] Hashiguchi, K., Cyclic plasticity models: critical reviews and assessments, in: Elastoplasticity Theory, Springer Berlin Heidelberg, (2014) 187–202. http://link.springer.com/chapter/10.1007/978-3-642-35849-4_8 (accessed January 30, 2016). [10] Preumont, A., Piefort, V., Predicting random high-cycle fatigue life with finite elements, J. Vib. Acoust., 116 (1994) 245–248. doi:10.1115/1.2930420. [11] de la Fuente, E., Von Mises stresses in random vibration of linear structures, Computers & Structures, 87 (2009) 1253– 1262. doi:10.1016/j.compstruc.2009.06.008. [12] Benasciutti, D., Tovo, R., Comparison of spectral methods for fatigue analysis of broad-band Gaussian random processes, Probabilistic Engineering Mechanics, 21 (2006) 287–299. doi:16/j.probengmech.2005.10.003. [13] Cláudio, R.A., Reis, L. , Freitas, M., Biaxial high-cycle fatigue life assessment of ductile aluminum cruciform specimens, Theor. Appl. Fract. Mech., 73 (2014) 82–90. doi:10.1016/j.tafmec.2014.08.007. << /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