【溫馨提示】壓縮包內(nèi)含CAD圖有下方大圖片預(yù)覽,下拉即可直觀呈現(xiàn)眼前查看、盡收眼底縱觀。打包內(nèi)容里dwg后綴的文件為CAD圖,可編輯,無(wú)水印,高清圖,壓縮包內(nèi)文檔可直接點(diǎn)開預(yù)覽,需要原稿請(qǐng)自助充值下載,所見才能所得,請(qǐng)見壓縮包內(nèi)的文件及下方預(yù)覽,請(qǐng)細(xì)心查看有疑問(wèn)可以咨詢QQ:11970985或197216396
A probabilistic method to predict fatigue crack initiation
SALIL. S. KULKARNI,? L. SUN, B. MORAN, S. KRISHNASWAMY and
J.D. ACHENBACH
Department of Mechanical Engineering, Northwestern University, Evanston, IL-60208, USA
?Author for Correspondence. (E-mail: sail@northwestern.edu)
Received 8 August 2005; accepted in revised form 6 September 2005
Abstract. A probabilistic method to predict macrocrack initiation due to fatigue damage is presented in this paper. Acoustic non-linearity is used to quantify pre-macrocrack initiation damage. This data is then used in a probabilistic analysis of fatigue damage. The probabilistic fatigue damage analy- sis consists of a suitably chosen damage evolution equation to model accumulated damage coupled with a procedure to calculate the probability of macrocrack initiation. The probability of macrocrack initiation is evaluated using the Monte Carlo Method with Importance Sampling. Numerical results for the probabilistic assessment of fatigue damage for a sample problem are presented and compared with experimental results.
Key words: Acoustic non-linearity, crack initiation, damage evolution, fatigue, probabilistic method, struc- tural health monitoring.
1. Introduction
Conventional procedures for life prediction of components subjected to fatigue are generally based on the ‘safe-life’ approach (see e.g. Suresh, 1991), coupled with Palmgren (see Palmgren, 1924) and Miner (see Miner, 1945) rules of linear cumu- lative damage. In the ‘safe-life’ approach for metal fatigue, life prediction is based on data from fatigue testing of components. Components of a structure are replaced when the probability of failure reaches a prescribed (often small) value, even though some of them may have a signi?cant remaining life. Hence it is a conservative approach with an economic penalty. To avoid this penalty, the ‘damage-tolerant’ approach is often a suitable alternative for life predictions. This approach is espe- cially useful when the rate of damage accumulated is well understood and can be monitored with a technique of quantitative non-destructive evaluation. However, in certain materials, such as high strength steels, like those used for drive trains of rotor- craft, critical damage in the form of a crack of detectable but very small length, often occurs late in the lifetime of a component. When a detectable crack has developed out of microscopic damage processes, it grows to an unacceptable length in a time that is short as compared to the total lifetime of the component. Periodic monitor- ing by a condition monitoring system can, however, signi?cantly improve the safety of the damage tolerance approach. Particularly, if pre-crack damage can be moni- tored and related to crack formation by an analytical fatigue damage procedure, sub- stantial safety and cost bene?ts can be gained. With this objective, a methodology to
Probabilistic method to predict fatigue crack initiation 17
calculate the probability of macrocrack initiation based on periodic in situ measure- ments is presented in this paper. The methodology essentially consists of two parts:
? A structural health monitoring system to monitor the fatigue process. In the pres-
ent approach, acoustic second harmonic generation (see Section 2) is used to mon- itor the damage accumulated during the fatigue process.
? A probabilistic fatigue damage procedure which takes the data from the health
monitoring system as its input and calculates the probability of macrocrack initia- tion as its output. The probabilistic damage procedure consists of a suitably cho- sen damage evolution equation to model accumulated damage (see Section 3) and a procedure to calculate the probability of macrocrack initiation (see Section 4). The procedure to calculate the probability of macrocrack initiation is capable of updating its forecast based on the latest inspection data from the monitoring sys- tem.
Using this methodology, the probability of macrocrack initiation is calculated for a sample problem taken from the literature, and the results are presented in Section 5.
2. Acoustic second harmonic generation
To quantify the damage accumulation in a component undergoing fatigue, it is ?rst necessary to relate the accumulated damage to an observable variable. The accumu- lated damage is caused by changes in the microstructure of the component which in turn introduce material non-linearity in the specimen. To characterize the material non-linearity, a single frequency ultrasonic wave generated by a transducer (gener- ator) is transmitted through the specimen and the signal received by the receiving transducer is analyzed. The material non-linearity distorts the single frequency wave and leads to generation of second and higher harmonics. As a result, the signal received at the receiver not only consists of a component at the transmitting (funda- mental) frequency but also second and higher harmonics whose amplitudes increase as the material non-linearity increases. Therefore, the material non-linearity and the underlying accumulated damage can be characterized by the ratio of the amplitude of the second harmonic A2 to the amplitude of the fundamental A1. This ratio is henceforth referred to as A2/A1 and is expected to increase as damage accumulation progresses. Note that the characterization of the acoustic non-linearity by A2/A1 (see Ogi et al., 2001) used here is different from the one usually found in the literature (see e.g. Moris et al., 1979).
Ogi et al. (2001) have used an electromagnetic acoustic transducer to measure A2/A1 in a specimen undergoing fatigue. As opposed to the procedure described ear- lier to measure the ratio A2/A1 – transmit a single (fundamental) frequency wave and measure the signal amplitude at the fundamental and the second harmonic, Ogi and his co-workers used a different but analogous procedure. They transmitted two differ- ent signals separately, one at resonance frequency fr and the other at frequency fr/2, into the specimen. For the case when the signal was transmitted at frequency fr, they measured the amplitude of the component of the received signal with frequency fr denoted A1(fr). For the case when the signal was transmitted at frequency fr/2 they measured the amplitude of the component of the received signal with frequency fr denoted by A2(fr/2). Observe that A2(fr/2) is the amplitude of the second harmonic. Note that this procedure ensures that both the measured signals are measured at
=
Figure 1. Typical evolution of the ratio A2(fr/2)/A1(fr) for 0.25 mass % C steel, Nf 56,000 (from Ogi et al., 2001).
the transducer resonance frequency fr, leading to more accurate measurements. For more details refer to Ogi et al. (2001). They observed that the ratio A2(fr/2)/A1(fr) increases nearly monotonically, and shows a distinct peak at the point of macrocrack initiation (see Figure 1).
This suggests that the state of damage in a specimen during fatigue can be tracked by measuring the ratio A2/A1 for different cycle numbers.
3. Damage model
=
=
This section presents a damage model whose evolution represents the evolution of the ratio A2(fr/2)/A1(fr) (see Section 2) up to the point of macrocrack initiation. The state of damage in a specimen at a particular cycle during fatigue is represented by a scalar damage function D(N), where N is the number of cycles. The magnitude D 0 corresponds to no damage, and D 1 corresponds to the appearance of the ?rst mac- rocrack. The following phenomenological model is assumed to represent the evolu- tion of the damage with the number of cycles
?? 1 . Oσ/2 ? rc(σˉ ) Σm 1
dD
?
dN = Nc
Oσ/2
(1 ? D)n if Oσ/2 > rc(σˉ )
0 if Oσ/2 < rc(σˉ )
(1)
ˉ
ˉ
Here, Nc is a normalizing constant, Oσ is the stress range in a cycle, rc(σ ) is the endurance limit when the mean stress in a cycle is σ , and m> 0 and n> 0 are param- eters which depend on the material and the loading conditions and are determined by correlating the evolution of the ratio A2/A1 to the evolution of D. It is assumed that rc(σˉ ) follows the Goodman relation (see Goodman, 1899), i.e.
σult
rc(σˉ ) = rc(0) .1 ? σˉ Σ
where σult is the ultimate tensile strength of the material. Assuming that Oσ and σˉ are constant during cycling and that Oσ/2 is always greater than rc(σˉ ), Equation (1)
can be solved to obtain
D(N) = 1 ?
(1 ? D0)n+1 ?
Σ
N ? N0 . Oσ/2 ? rc(σˉ ) Σm
Nc
Oσ/2
1
Σ
. (2)
n+1
(n + 1)
Here D0 is the damage present in the specimen at cycle N0.
4. Probability of macrocrack initiation
This section describes the procedure for calculation of the probability of macrocrack initiation and is based on the model described in Section 3. To begin, some of the quantities appearing in Equation (2) are taken as random variables. The probability of macrocrack initiation Pma, which is de?ned as the probability that the number of cycles to macrocrack initiation, will be less than a speci?ed number of cycles Ns, is then given by
Pma = 1 ? P r(D(Ns)< 1). (3)
From Equation (1), it follows that
D(Ns)< 1
implies
Σ(1 ? D0)n+1 ?
Ns ? N0 . Oσ/2 ? rc(σˉ ) Σm
1
Σ
(n + 1) n+1
> 0.
Nc Oσ/2
From the previous equation it is observed that D(Ns)< 1 is equivalent to
Ns < Nini (4)
. Σ
where
Nini = N0 +
Nc
n + 1
(1 ? D0)n+1
Oσ/2 m
.
Oσ/2 ? rc(σˉ )
=
Note that in a purely deterministic setting, Nini corresponds to the cycle at which a macrocrack appears. This can be seen by setting D 1 in Equation (2) and solving for N .
The probability of macrocrack initiation can therefore be calculated as
Pma = 1 ? P r(D(Ns)< 1)
= 1 ? Pr(Ns < Nini)
= Pr(Nini < Ns). (5)
To calculate Pma given by Equation (5), one ?rst de?nes a limit state surface
g = Nini ? Ns.
=
Also let X [X1 X2 ... Xk]T denote the random quantities in Equation (2) (for example, in the present paper we consider a problem with known constant stress
?? < 0 : domain of macrocrack initiation
cycle amplitude Oσ while the quantities rc(σˉ ), D0, m and n are considered ran- dom with known probability distribution. We then have X = [X1 X2 X3 X4]T = [rc(σˉ ) m n D0]T ). Observe that
g(X)
0 : limit state surface
=
∫
? > 0 : safe domain
The probability of macrocrack initiation Pma is therefore given by
Pma =
g(x)<0
fx(x)dx, (6)
=
=
where fX(x) denotes the joint probability distribution of X. To evaluate the integrals appearing in Equation (6), the random variables are ?rst mapped via a Rosenblatt transformation (see Rosenblatt, 1952) into a standard Gaussian space where the random variables denoted by U [U1 U2 ... Uk]T are independent, normally distrib- uted and have zero mean and unit standard deviation. The modi?ed Hasofer-Lind, Rackwitz-Fiessler (HL-RF) algorithm described in Kiureghian and Liu (1989) is used to obtain the point closest to the origin u? on the surface g(u) 0. In the modi- ?ed HL-RF algorithm, one adjusts the step size during each iteration to obtain a suf?cient decrease in the merit function which is based on the ?rst-order optimal- ity conditions. Monte Carlo integration with importance sampling, with the sampling density centered at u? (see Fujita and Rackwitz, 1988) is then used to calculate the integrals in Equation (6).
5. Sample problem
To demonstrate the application of the ideas presented in the paper, the probability of macrocrack initiation is calculated for a sample problem. The data, i.e. the ratio A2/A1 as a function of the number of cycles, is obtained from Ogi et al. (2001). Ogi and his co-authors have performed a rotating bending fatigue test with a four point bending con?guration on a 0.25% C (mass) steel and have obtained the variation of the ratio A2/A1 with the number of cycles (see Figure 1). The yield strength of the material is 333 MPa specimen and it is subjected to a maximum bending stress of 280 MPa. The endurance limit of the material at zero mean stress, rc(0), is assumed to be 180 MPa.
To simulate the inspection process, A2/A1 is obtained from Figure 1 at different number of cycles with the interval between the cycles getting progressively shorter. Table 1 shows the ratio A2/A1 measured as a function of number of cycles (see Figure 1).
×
It is observed that the evolution of the ratio A2/A1 is not strictly monotonic dur- ing the initial stages of fatigue. The ‘raw’ damage is obtained by normalizing the non- linearity measurements by the expected maximum value of the non-linearity which for the given problem is assumed to be 3.25 10?3 (see Figure 1). The ‘raw’ dam- age is tabulated in the third column in Table 2. Note that to be consistent with the assumed monotonic damage evolution (see Section 3), the ‘raw’ damage is corrected and tabulated as ‘corrected’ damage in the fourth column in Table 2.
Table 1. Measured values of the ratio A2/A1 during successive inspections.
j
Ninspj
(A2/A1)j × 10?3
0
0
0.90
1
11200
0.80
2
22400
0.90
3
26880
1.50
4
30800
2.00
5
33040
2.50
6
34000
3.10
Table 2. Raw and ‘corrected’ values of damage during suc- cessive inspections.
j
Ninspj
‘Raw’ damage
‘Corrected’ damage Dinspj
0
0
0.28
0.28
1
11200
0.25
0.28
2
22400
0.28
0.28
3
26880
0.46
0.46
4
30800
0.62
0.62
5
33040
0.77
0.77
6
34000
0.95
0.95
These damage values which are based on the experimental results of Ogi et al. (2001) are used for two different purposes in the paper. In the ?rst place, they are used to determine the statistics of parameters m and n in Equation (1) (see Section 3). Secondly, they are interpreted as damage measurements obtained by structural health monitoring system for input to the procedure to calculate the probability of macrocrack initiation (see Section 4).
5.1. Prior information
=
For calculating the probability of macrocrack initiation using the experimental data presented earlier, it is assumed that prior information regarding the parameters m and n is available. This is a reasonable assumption if appropriate experiments have been performed in the laboratory. In that case, the statistics for m and n can be cal- culated by ?tting the damage model to the experimental datasets using non-linear regression (see e.g. Draper and Smith 1998). Notice that in the present case, only one experimental data set (Dinspj ,j 0,... , 6) is available. To calculate the statis- tics for m and n additional datasets need to be generated. The additional datasets are generated by assuming that at inspection cycle Ninspj , the ‘corrected’ damage has a truncated normal distribution in the range (0, 1) with mean of the parent normal
distribution equal to Dinspj and standard deviation equal to 0.1, i.e.
. Σ?φ
D0 μ
σ
=
f (D0)
σ
Σ8 . 1?μ Σ ? 8
. μ?0 ΣΣ
σ
0 ? D0 ? 1
σ
= =
with μ Dinspj and σ 0.1. Here φ and 8 are the probability density and the cumu- lative distribution function of for the normal distribution. The damage model is ?tted to each of the additional datasets, i.e. m and n are calculated for each of the datasets. The statistics of m and n are then compiled. For the current problem, the simulations suggest that m and n are dependent parameters with means 2.9074 and 2.2295 respec- tively and the coef?cient of correlation equal to 0.9029. The standard deviations of m and n are found to be 0.5336 and 0.7129, respectively.
For calculating m and n using non-linear regression, the following values are used:
Oσ = 2×280 MPa, σˉ = 0 MPa, Nc = 10, 000 and rc(0) = 180 MPa. The value of Nc,
which is a normalizing constant is assumed to be of the order of the number of cycles to failure. For this particular problem, failure is observed at 56,000 cycles. Here Nc is taken to be 10,000.
5.2. Calculation of probability of macrocrack initiation
The probability of macrocrack initiation is now calculated using the data given in Table 2. The procedure essentially consists of recording the inspection data at, say, cycle Ninspj and then calculating the corresponding value of the damage, say, Dinspj . If Dinspj < 1 then the probability of macrocrack initiation is calculated for Ns > Ninspj as described in Section 4.
ˉ
For the problem under consideration, it is assumed that Oσ , σ , and Nc are ?xed
ˉ
quantities while rc(0), D0, m and n are random variables. The ?xed quantities Oσ , σ and Nc, are assumed to have the same values as those used for the non-linear regres- sion (see Section 5.1). The random variable rc(0) is assumed to have a log-normal distribution with a mean of 180 MPa and a standard deviation of 5.4 MPa. The ran- dom variable D0 at cycle Ninspj is assumed to have a truncated normal distribution in the interval (0, 1) with mean equal to Dinspj and standard deviation equal to 0.1. This distribution of D0 is consistent with the one used to generate the statistics of m and n (see Section 5.1). The random variables m and n are assumed to have a joint lognormal distribution with the statistics given in the Section 5.1. Starting from the recording of the third inspection, the probability of macrocrack initiation calculated for different cycles is given in Table 3.
It is seen from Table 3 that the probabilistic fatigue damage analysis proce- dure updates the probability of macrocrack initiation with successive inspections. The cumulative probability distribution for macrocrack initiation with inspections is shown in Figure 2.
The numbers (1,..., 6) in the ?gure indicate the cycle numbers at which inspection is recorded. Observe that immediately following an inspection, the probability of mac- rocrack initiation drops to zero. This is due to the fact that it is implicitly assumed in the probabilistic procedure that a macrocrack if present, will be de?nitely detected by the inspection procedure.
Table 3. Calculation of Pma.
Cycles Pma
Ns 3rd Insp 4th Insp 5th Insp 6th Insp
(Ninsp = 26880) (Ninsp = 30800) (Ninsp = 33040) (Ninsp = 34000)
30,000 0.09304 0.00000 0.00000 0.00000
31,000 0.15548 0.02922 0.00000 0.00000
32,000 0.22702 0.21027 0.00000 0.00000
33,000 0.30400 0.39166 0.00000 0.00000
34,000 0.38180 0.54330 0.62907 0.00000
35,000 0.45845 0.66148 0.81884 0.96343
40,000 0.74368 0.93052 0.99162 0.99989
6
5
1
2
4
3
1
0.9
0.8
0.7
Pma --- >
0.6
0.5
0.4
0.3
0.2
0.1
0
0 0.5 1 1.5 2 2.5 3 3.5 4
N --- >
x 104
Figure 2. Cumulative probability distribution for macrocrack initiation.
Note that for the fatigue problem described, the ?rst macrocrack is observed at approximately 34,160 cycles (see Ogi et al., 2001). As seen from Table 3 and Figure 2, the formation of the macrocrack is predicted quite well.
6. Summary
A methodology to predict macrocrack initiation which has the following features has been presented:
? use of acoustic non-linearity to quantify damage
? heuristic damage growth law
? material/other parameters treated as random variables
? periodic measurements to assess state of damage and update state of damage
? calculation of probability of macrocrack initiation
Reasonable results are obtained for a sample problem solved using the presented methodology.
7. Acknowledgements
This work was supported by the Federal Aviation Administration, Rotorcraft Struc- tural Integrity and Safety, under Award AGMT: 01-C-AW-NWU/A04 to the Center for Quality Engineering and Failure Prevention a