Microsoft Word - numero_46_art_9 M.F.M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 84 Developments in the fracture and fatigue assessment of materials and structures Probabilistic-based analysis for damaging features of fatigue strain loadings M. F. Mod Yunoh, S. Abdullah Centre for Automotive Research, Faculty of Eng. and Built Enviroments, Universiti Kebangsaan Malaysia faridzyunus@gmail.com, shahrum@ukm.edu.my M. H. M. Saad, Z. M. Nopiah, M. Z. Nuawi, S. S. K. Singh Dept. of Mechanical and Materials, Faculty of Eng. and Built Enviroments, Universiti Kebangsaan Malaysia hanifsaad@ukm.edu.my, zmn@ukm.edu.my, mzn@ukm.edu.my, salvinder@ukm.edu.my ABSTRACT. This paper presents the behaviour of fatigue damage extraction in fatigue strain histories of automotive components using the probabilistic approach. This is a consideration for the evaluation of fatigue damage extraction in automotive components under service loading that is vital in a reliability analysis. For the purpose of research work, two strain signals data are collected from a car coil spring during a road test. The fatigue strain signals are then extracted using the wavelet transform in order to extract the high amplitude segments that contribute to the fatigue damage. At this stage, the low amplitude segments are removed because of their minimal contribution to the fatigue damage. The fatigue damage based on all extracted segments is calculated using some significant strain-life models. Subsequently, the statistics-based Weibull distribution is applied to evaluate the fatigue damage extraction. It has been found that about 70% of the probability of failure occurs in the 1.0 x 10-5 to 1.0 x 10-4 damage range for both signals, while 90% of the probability of failure occurs in the 1.0 x 10-4 to 1.0 x 10-3 damage range. Lastly, it is suggested that the fatigue damage can be determined by the Weibull distribution analysis KEYWORDS. Fatigue damage; Features extraction; Probabilistic; Wavelet; Weibull distribution. Citation: Yunoh, M. F. M., Abdullah, S., Saad, M. H. M., Nopiah, Z. M., Nuawi, M. Z., Singh, S. S. K., Determining probabilistic-based failure of damaging features for fatigue strain loadings, Frattura ed Integrità Strutturale, 46 (2018) 84-93. Received: 18.01.2018 Accepted: 20.05.2018 Published: 01.10.2018 Copyright: © 2018 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 atigue life prediction under service loading remains a challenging problem in real life, especially in the automotive industry. In an automotive application, service loadings such as stress on a car wheel are in variable amplitude loading, as reported in Sonsino [1]. This situation leads to a need for developing a new approach to predict the reliability of F http://www.gruppofrattura.it/VA/46/9.mp4 M. F. M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 85 the components that are subjected to the fatigue loading. The fatigue time histories often contain a major percentage of small amplitude cycles, and the fatigue damage for these cycles can also be small. Therefore, in many cases, the fatigue loading history is edited by removing these small amplitude cycles to produce representative and meaningful, yet economical testing, Stephens et al. [2]. Several fatigue data editing techniques have been developed for use in the time domain analysis, Abdullah et al. [3]. Some previous algorithms have been developed for eliminating low amplitude cycles in order to only retain the high amplitude cycles, El Ratal et al. [4]. In the frequency domain, the time history is a low pass filtered by the criterion that high-frequency cycles have small amplitude, which is are not damaging. The filtering method does not shorten the signal because it does not provide the time-based information, Nizwan et al. [5]. In addition, Abdullah et al. [3] have developed a method for data editing to shorten the strain signal in the time-frequency domain. In fatigue data editing, the behaviour of extraction segments also needs to be studied because it contributes many bits of information that can improve fatigue life prediction. To deal with uncertainties and variations in fatigue data, the statistical analysis, i.e. the probability analysis is the best approach that should be adopted. Zhao and Liu [6] use the approach of statistical aspects of the S-N curve by means of the Weibull distribution. This approach indicates that an appropriate distribution determination is the primary task for a rolling contact fatigue analysis. Tiryakioglu [7] uses the Weibull analysis in fatigue data and predicts the failure mechanisms due to cracks initiating from surface and interior defects. The evaluation of fatigue damage extraction for automotive components under service loading is vital in the reliability analysis. Very few analysis methods have been developed to evaluate the fatigue damage extraction. A consistent description of the probability of damage occurrences is possible only if the damage distribution function is known. In this paper, two sets of strain signals from the real component in service are used in the fatigue feature extraction. The feature extraction, i.e. fatigue damage, is analysed using statistical inferences. The objective of this study is to determine the probabilistic-based failure of damage featured in the strain signals, and at the same time, validate them through the extraction process. By considering the significant statistical tools, features extraction is combined with the Weibull distribution analysis in order to obtain a better evaluation. THEORETICAL BACKGROUND The Global Statistic he time series contains an explanation of a set of variables taken at equally spaced time intervals. A statistical analysis is normally used to determine the random signals and monitor the pattern of the analysed signals. The calculation of the root-mean-square (r.m.s.) and the kurtosis is very important in the fatigue signals in order to retain a certain number of the signal amplitude range characteristics. The r.m.s. value is the 2nd statistical moment, is used to quantify the overall energy content of the signal. For discrete data sets the r.m.s. value is defined as: (1) While kurtosis is the signal 4th statistical moment, is a global signal statistic which is highly sensitive to the spikiness of the data. For discrete data sets the kurtosis value is defined as: (2) where jx is the amplitude of signal, n is number of data and x is the mean value. The Wavelet Transform The wavelet transform (WT) is defined in the time-scale domain and is a significant tool for analysing the time-localised features of a signal. It represents a windowing technique within the variable-sized region. A wavelet transform can be classified as either a continuous wavelet transform (CWT) or a discrete wavelet transform (DWT) depending on the discretisation of the scale parameters of the analysing wavelet. The DWT based on such wavelet functions is called the orthogonal wavelet transform (OWT). Orthogonal wavelet transforms are normally applied for the compression and feature 2/1 1 21 ..          n j jx n smr     n j j xx smrn K 1 4 4)..( 1 T M.F.M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 86 selection of signals. DWT is derived from discrete CWT, and it is shown as the following expression, Purushotham et al. [8]:  /2 *0 0 0( , ) ( ) ,m mW m n x t a a t nb dt        (3) where a and j are the scale factor, both b and k are the position, and Ψ is the mother wavelet. Oh [9] has previously conducted fatigue data analysis using the wavelet transform (WT) for spike removal, denoising, and data editing. Piotrkowski et al. [10] used the Wavelet Transform application in acoustic emissions to detect damage and corrosion. Fatigue Life Assessment The Palmgren-Miner linear cumulative damaging rule is normally associated with the established strain-life fatigue models Sun et al., [11]. The fatigue damage caused by each cycle of repeated loading is calculated by reference to material life curves, such as S N or N  curves. The fatigue damage caused by multiple cycles is expressed respectively as: 1 f D N         (4) i f N D N           (5) where D is fatigue damage for one cycle and D is total fatigue damage iN is the number of cycles within a particular stress range and its mean and fN is a number of cycles. The strain-life model commonly used for the prediction of fatigue strain life is the Coffin-Manson relationship model. This model can provide a traditional prediction when there is more compressive load time history and the mean stress is zero. The following equation can define this model:     ' 2 ' 2 b cf a f fN f N E     (6) E is the material modulus of elasticity, a is the true strain amplitude, 2 fN is the number of reversals to failure, ' f is the fatigue strength coefficient, b is the fatigue strength exponent, ' f is the fatigue ductility coefficient, c is the fatigue ductility exponent, m is the mean stress, and max is the maximum stress. The inclusion of mean stress effects in the life prediction makes the process more complex. The Morrow mean stress model is given by Dowling [12]:     ' ' ' 1 2 2 b cf m a f f f f N N E               (7) where is the total strain amplitude, ' f , b, ' f and c are considered to be material properties, fN is the number of cycles to failure, and m is the mean stress. Another strain life model dealing with mean stress effects is known as the Smith- Watson-Topper (SWT) model, and its equation is written as: 2 2 ' ( 2 ) ' ' ( 2 ) f b b c a mak f f f fN N E        (8) M. F. M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 87 The Coffin-Manson relationship only considers the damaging calculation at zero mean stress. However, in real-case scenarios, some of the realistic service situations involve nonzero mean stresses. For example, in a case of the loading being predominantly compressive, particularly for wholly compressive cycles, the Morrow mean stress correction effect provides more realistic life estimates and seems to work reasonably well for steels, Ince & Glinka [13]. The Weibull Distribution In terms of the statistical analysis used in engineering, the Weibull distribution is a theoretical model that has been successfully used to model the life data. The Weibull distribution is described by the shape, scale, and threshold parameters. The Weibull distribution model is a tool to develop the probabilistic analysis because of its ability to provide reasonably accurate failure analysis and failure forecasts with extremely small samples, Sivapragash, [14]. The 2-P Weibull distribution function is shown in Eqn. 7: (9) where f (x: θ: β) represents the probability of strain-life being equal to or less than x, θ is a scale parameter, and β is a shape parameter. θ and β are estimated by observation. The Weibull distribution has been widely utilised to develop a model of extreme values, such as failure time and fracture strength, Shalabh et al. [15]. The Weibull distribution is a probabilistic analysis which has been used for the determination of static and dynamic mechanical properties of materials. This distribution has the capability to model experimental data with different characters, Li et al. [16]. METHODOLOGY n this work, the strain signal has been collected from the coil spring in the car suspension system during a road test as shown in Fig. 1 which is collected in macrostrain ( e ). An established signal, known as SAESUS, is the typical strain history for the suspension system, developed by the Society of Automotive Engineers (SAE), as reported in Oh [9]. Another signal, S1, is a set of experimental strain signals measured at 500 Hz sampling rate using a strain gauge that is positioned on the coil spring component of a car being driven on a rural road at a velocity of 50 to 60 km/h. The measurement procedure is reported in their previous work by Yunoh et al. [17]. The fatigue signal was measured at the car coil spring which subjected to the road load service. All data were recorded as strain time histories and Fig. 2 shows the set- up of fatigue data measurement during the process. The strain value was measured using strain gauge and it was connected to the specific data logger, for data acquisition purpose. Experimental parameters such as sampling frequency and type of output data were then set using the common data logger interface. (a) (b) Figure 1: The strain signals used for the work: (a) SAESUS, (b) An experimental measured data, S1. The extraction algorithm of the strain signal is developed for fatigue data editing purposes. The discrete wavelet transform is utilised to identify the high amplitude segments in the fatigue signal due to the high-energy coefficients magnitude. The magnitudes of the wavelet energy coefficients in the time-frequency domain are transposed into time histories representations to trace the location of the high amplitude segments. The respective magnitudes are obtained from the accumulation of wavelet transform magnitude distributions along the frequency band at each time interval. Thus, the energy coefficients in time representation are gained, and these signals are used to detect the presence of high amplitude events in the fatigue signal.                          xx xf exp),:( 0 20 40 60 80 100 120 -1,000 -800 -600 -400 -200 0 200 400 600 800 1,000 T ime (s) A m p li tu d e (u e) 0 20 40 60 80 100 120 -500 -400 -300 -200 -100 100 200 300 400 500 T ime (s) A m p li tu d e (u e) I M.F.M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 88 The Cut-off Level (COL) needs to be set up for the elimination process in a high amplitude event extraction process. The COL is the minimum percentage of the wavelet energy coefficients to be retained and the segments with magnitudes lower than the COL value will be removed. The retained segments are then sliced from the original signal. The sliced segment identification is performed by means of a search which identifies the point at which the signal envelope inverts from the decay behaviour. The two inversion points, one on either side of the peak value, define the temporal extent of the sliced segment. The original strain signal is then edited to remove the low amplitude cycle contained in the signal based on the time location of the wavelet transform spectrum-sliced segment. The fatigue damage for every retained segment is then calculated for the probabilistic analysis. Figure 2: A diagrammatic process flow for fatigue signal collection The fatigue damage for all segments is further analysed using the 2-P Weibull distribution to evaluate the segments’ extraction results. The advantage of using this distribution lies in its ability to explain the simple function. This approach is often used in assessing the fatigue life of the material based on its simple calculation, Glodez et al. [18]. All the fatigue damage of the retained segments is incorporated into the Weibull distribution, and as has been later suggested, based on the probability distribution, the resultant significant findings are utilised to reveal the inferences of this study. RESULTS AND DISCUSSIONS irst, the strain signals were evaluated by observation based on statistics characteristics, i.e. number of cycles counted, root-mean-square (r.m.s), kurtosis and total damage as tabulated in Tab. 1. According to the results, the SAESUS signal produced higher values of r.m.s and kurtosis. This is due to the higher amplitude segments existing in the SAESUS signal compared to the S1 signal. Higher amplitude segments contributed more energy in the oscillatory signals. The value of kurtosis for the SAESUS strain signal was found to be 4.32, which indicates that the spike and extreme values exist in the signal. This result was consistent with the r.m.s value for SAESUS, which is found to be 246.6 με, higher when compared to S1. The kurtosis value for both strain signals was not found in a Gaussian distribution because in a stationary Gaussian process, the kurtosis value is approximately 3. The fatigue damage for both signals was calculated based on three strain-life models as tabulated in Tab. 1. The values of the fatigue damage based on the three models were close to each other with minor differences. The fatigue damage values for SAESUS were found to be higher compared to S1. The existence of the higher amplitude segments in the SAESUS signal contributed more fatigue damage compared to S1. F Strain Gauge Strain gauge position Strain gauge connect to data acquisition Data transfer from data acquisition to computer Strain signal in time domain M. F. M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 89 Data N (Cycles) r.m.s (με) Kurtosis Fatigue Damage (Damage/Block) SAESUS 1253 246.6 4.28 CF Morrow SWT 1.70 x 10-3 1.70 x 10-3 1.64 x 10-3 S1 9768 122.68 2.82 CF Morrow SWT 2.97 x 10-4 2.94 x 10-4 2.85 x 10-4 Table 1: Summary of statistical characteristic of the signals The developed wavelet transform algorithm transformed the signal into a time-frequency domain to obtain the time- frequency localisation, as shown in Fig. 3. The spectrum colour intensity is proportional to the absolute energy coefficients values because it delivers the energy distribution display on time and frequency. The spectrum of SAESUS and S1 signals shows that the intensity colours of the SAESUS strain signal are more highlighted compared to S1. A large scale was indicative of low frequency, and higher amplitude segments indicated that the cycles had higher energy. This means that it can inflict higher fatigue damage. (a) (b) Figure 3: The wavelet transform spectrum of; (a) SAESUS, (b) S1. Higher amplitude segments based on wavelet transform spectrum were extracted. The location of the retained segments was shown in Fig. 4. After eliminating the low amplitude segments, the high amplitude segment should be retained to maintain the fatigue damage contained in both signals. The SAESUS strain signal contributed 21 high amplitude retained segments. Meanwhile, the S1 strain signal contributed 11 high amplitude retained segments. Based on the extraction results, the SAESUS strain signal consists of higher amplitude segments compared to S1 because it retained the higher amplitude segments after the extraction process. (a) (b) Figure 4: The location of the retained segments; (a) SAESUS, (b) S1. The features extraction, i.e. fatigue damage for every segment, was calculated to analyse and determine the probability of failure based on high amplitude segments. It is helpful to know the fatigue damage distribution based on the extracted segments for both signals. The 2-P Weibull distribution was used for the data analysis. According to Tab. 2, the values of the shape parameters, β, for both signals were less than 1.0, indicating a decreasing failure rate. According to Finkelstein, [19], if the shape parameter for Weibull distribution is less than 1.0, this distribution is always decreasing the failure rate. The value of scale parameter, θ, of fatigue damage for both signals was found to be 5.3 x 10-5 and 1.4 x 10-5. This means that when the components reach these values, the probability of failure will increase. Dat a S ca le ( 1 /H z) 0.5 1 1.5 2 2.5 x 10 4 1 27 53 79 105 131 157 183 209 235 257 50 100 150 200 250 Dat a S ca le ( 1 /H z) 1 2 3 4 5 6 x 10 4 1 27 53 79 105 131 157 183 209 235 257 50 100 150 200 250 0 20 40 60 80 100 120 -1,000 -800 -600 -400 -200 0 200 400 600 800 1,000 A m p li tu d e (u e) T ime (s) 0 20 40 60 80 100 120 -500 -400 -300 -200 -100 0 100 200 300 400 500 A m p li tu d e (u e} T ime (s) M.F.M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 90 Signals Shape parameter, β Scale parameter, θ SAESUS 0.59 5.3 x 10-5 S1 0.34 1.4 x 10-5 Table 1: Statistical data analysis based on Weibull distribution. Cumulative Distribution Function (CDF) represents the probability of the components to fail at the particular time interval. Based on the plots in Fig. 5, the fatigue damage for S1 increases more rapidly in cumulative probability values compared to the SAESUS. This finding can be suggested that SAESUS signal consists high amplitude segments that cause high fatigue damage compared to S1. Figure 5: Probability density functions for SAESUS and S1 strain signal Figure 6: The Weibull probability distribution plots of fatigue damage for SAESUS Figs. 6 and 7 show the best fit of the Weibull probability distribution plot for fatigue damage of high amplitude segments for both signals. The parameters of the distribution are estimated using the maximum likehood method. The Weibull probability distribution plot for both strain signals was estimated using 95% confidence level correction line, which includes the true population parameter. From a statistical point of view, the more test data that falls into the 95% confident level, the more reasonable the probabilistic distribution is, according to Xianmin et al [20]. SAESUS 90% 70% M. F. M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 91 Based on the plot indicated in Fig. 5 for the SAESUS strain signal, there are four data points which are out of the 95% confidence level. This is because several segments in the signal consist of very high amplitude segments. Meanwhile, all points in S1 strain signal were distributed at the 95% confidence level. Figs. 5 and 6 also show that 70% of the probability of failure occurred from 1.0 x 10-5 to 1.0 x 10-4 fatigue damage for both signals. Around 90% of the probability of failure occurred from 1.0 x 10-4 to 1.0 x 10-3 fatigue damage for both signals. Figure 7: The Weibull probability distribution plots of fatigue damage for S1. Figs. 8 and 9 show the best fit of the Weibull probability distribution plot for fatigue life of high amplitude segments for both signals. All of data for the SAESUS and S1 strain signal were distributed at the 95% confidence level. Fig. 7 show that 70% of the probability of failure for SAESUS occurred at 5.2 x 103 cycles. Meanwhile the probability of failure for SAESUS occurred at 6.9 x 103 cycles. Around 90% of the probability of failure occurred at 1.6 x 105 to 3.5 x 105 fatigue life for both signals. Thus, the value of fatigue damage and life based on extraction segments approaching the failure can be determined roughly by the Weibull analysis. Figure 8: The Weibull probability distribution plots of fatigue damage for SAESUS. 70% 90% 70% 90% M.F.M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 92 Figure 9: The Weibull probability distribution plots of fatigue damage for S1. CONCLUSION o summarise the statistics and strain life analysis, the SAESUS strain signal contributes more damaging segments compared to S1. This situation is due to the high amplitude segments obtained in the signal. The fatigue damage based on high amplitude segments fits well to a Weibull distribution. Based on the results, it shows that 70% of the probability of failure occurred between 1.0 x 10-5 and 1.0 x 10-4 fatigue damage for both signals. Around 90% of the probability of failure occurred between 1.0 x 10-4 to 1.0 x 10-3 fatigue damage for both signals. Thus, the value of fatigue damage based on extraction segments that are approaching the failure can be determined using the Weibull analysis with the 95% confident level, which is a reasonable probabilistic distribution. This approach presents an alternative technique to predict the fatigue damage probability function for automotive components. However, further works and analyses need to be performed for the purpose of validation, and also to produce high data accuracy by means of life assessment. REFERENCES [1] Sonsino, C. M. (2006). Fatigue testing under variable amplitude loading, International Journal of Fatigue, 29(6), pp. 1080-1089. [2] Stephens, R. I., Dindinger, P. M., and Gunger, J. E. (1997). Fatigue damage editing for accelerated durability testing using strain range and SWT parameter criteria, International Journal of Fatigue, 19(8-9), pp. 599-606. [3] Abdullah, S., Nizwan, C. K. E., Yunoh, M. F. M., Nuawi, M. Z., and Nopiah, Z. M. N. (2013). Fatigue features extraction of road load time data using the S-Transform, International Journal of Automotive Technology, 14(5), pp. 805-815. [4] El-Ratal, W., Bennebach, M., Lin, X. and Plaskitt, R. Fatigue life modelling and accelerated test for components under variable amplitude loads, Proc. of Symposium on Fatigue Testing and Analysis under Variable Amplitude Loading Conditions, Tours, France, 2002. [5] Nizwan, C. K. E., Abdullah, S., Nuawi M. Z. and Lamin, F. (2007). A study of fatigue data editing using frequency spectrum filtering technique, Proceeding of the World Engineering Congress (WEC), Penang, Malaysia, 2007, pp. 372- 378. [6] Zhao, Y. X. and Liu, H. B. (2014). Weibull modelling of the probabilistic S-N curves for rolling contact, International Journal of Fatigue, 66, pp. 47-54. [7] Tiryakioglu, M. (2015). Weibull analysis of mechanical data for casting II: Weibull mixtures and their interpretation, Metallurgical and Materials Transactions A, 46A, pp. 270-280. T 70% 90% M. F. M. Yunoh et alii, Frattura ed Integrità Strutturale, 46 (2018) 84-93; DOI: 10.3221/IGF-ESIS.46.09 93 [8] Purushotham, V., Narayanan, S. and Prasad, S. A. N. (2005). Multi-fault diagnosis of rolling bearing elements using Wavelet analysis and hidden Markov model-based fault recognition, NDT & E International, 38, pp. 654-664. [9] Oh, C. S. (2001). Application of Wavelet transform in fatigue history editing, International Journal of Fatigue, 23(3), pp. 241-250. [10] Piotrkowski, R.,Castro, E and Gallego, A. (2009). Wavelet power, entropy and bispectrum applied to AE signals for damage identification and evaluation of corroded galvanised steel, Mechanical System and Signal Processing, 23(2), pp. 432-445. [11] Sun, Q., Dui, H. N. and Fan, X. L. (2014). A statistically consistent fatigue damage model based on Miner’s rule, International Journal of Fatigue, 69, pp. 16-21. [12] Dowling, N. E. (2013). Mechanical Behavior of Materials, 4th ed.: Pearson Education Limited. [13] Ince, A. and Glinka, G. A modification of Morrow and Smith-Watson-Topper mean stress correction models, Fatigue Fracture Engineering Materials and Structure, 34, pp. 854-867. [14] Sivapragash, M., Lakshminarayanan, P. R., Karthikeyen, R., Raghukanda, K. and Hanumantha, M. (2008). Fatigue life prediction of ZE41A magnesium alloy using Weibull distribution, Material & Design, 29 pp.1549–53. [15] Shalabh G, Dheeraj S.S, Asok R. Statistical pattern analysis of ultrasonic signals for fatigue damage detection in mechanical structures, NDT & E International, 41, pp. 491-500. [16] Li W, Sakai T, Wang P. (2011). Statistical analysis of fatigue crack growth behavior for grade B cast steel. Material and Design, 32, pp. 1262-1272. [17] Yunoh, M. F. M., Abdullah, S., Saad, M. H. M., Nopiah, Z. M. and Nuawi, M. Z. (2014). Fatigue time history analysis for determining the strain signal behavior, International Journal Vehicle System Modelling and Testing, 9 (3/4), pp.363- 371. [18] Glodez, S., Sori, M., and Kramberger. J. (2013). Prediction of micro-crack initiation in high strength steels using Weibull distribution. Engineering Fracture Mechanics, 108, pp. 263-274. [19] Finkelstein, M. (2009). Understanding the shape of the mixture failure rate (with engineering and demographic applications), Applied Stochastic Models in Business and Industry, 25(6): 643-663. [20] Xianmin, C., Qin, S., Hongna, D. and Junling, F. (2016). A statistically self-cnsistent fatigue damage accumulation model including load sequence effects under spectrum loading, Frattura ed Integrita Strutturale, 38, pp. 319-330. << /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