SElogo.gif (5955 bytes)

 

"Probabilistics" 3

There is more to Monte Carlo simulation than replacing constants with probability densities.

Monte Carlo Modeling Specifics:

After fitting individual Paris equations to each of the 68 specimens, the mean and standard deviation for the individual Paris parameters, intercept, C, and slope, n, were computed. The well-known Paris model for fatigue crack growth is given in equation 1.

dadn_dk.gif (519 bytes)  (1)

where da/dN is the crack growth rate, in mm per cycle, and DK is the applied stress intensity factor, in MPam, given by equation 2.

   dk.gif (683 bytes)  (2)

Here, Ds is the testing stress range, smax - smin, a is the crack length, and f( ) is a function of the specimen (or component) geometry and crack length. Of course, when equation 1 is plotted on a log-log grid this is a straight line with intercept C and slope n.

Assuming for the sake of simplicity that there was no variation in the starting crack size, the final crack size, or the test stress, the calculated cyclic lifetime can be computed from the individual Paris fits using equation 3.

eqns.gif (1943 bytes)

eqn3.gif (2924 bytes)            (3)

In practice this integration is usually carried out numerically.

To conduct the usual MC simulation Ni is computed from h(Ci , ni) where h( ) is equation 3, and i ranges from 1 to say 1000 (or 10 000).

Many MC practitioners then calculate a mean and standard deviation for N, or log10(N), report the results and stop there, since there is nothing against which to compare the distribution of computed values for Ni.

In this case, however, Virkler's data show the observed distribution of actual specimen lives and thus provide a direct comparison for these calculations.

continued ...          

SElogo.gif (5955 bytes)

[ Home ] [  Probs 1  Probs 2  SEbullet_2.gif (863 bytes)  Probs 4  Probs 5  Probs 6  Download ]

Mail to Charles.Annis@StatisticalEngineering.com
Copyright � 1998-2008 Charles Annis, P.E.
Last modified: June 08, 2014