Modeling Reliability of Ceramics Under Transient Loads and Temperatures

1. Modeling Reliability of Ceramics Under Transient Loads and Temperatures Noel N. Nemeth Osama M. Jadaan Eric H. Baker The 26th Annual International Conference on Advanced Ceramics & Composites January 13-18, 2002, Cocoa Beach Florida E-mail: Noel.N.Nemeth@grc.nasa.gov Glenn Research Center Life Prediction Branch at Lewis Field

2. Outline • Objective • Background - CARES/Life - References to previous work • Theory - Power law & Walker law - Computationally efficient method for cyclic loading • Examples - Diesel exhaust valve - Alumina static fatigue • Conclusions

3. Objective • Develop a methodology to predict the time-dependent reliability (probability of survival) of brittle material components subjected to transient thermomechanical loading, taking into account the change in material response with time. • Transient reliability analysis

4. Fully Transient Component Life Prediction • MOTIVATION: • To be able predict brittle material component integrity over a simulated engine operating cycle • REQUIRES: • Life prediction models that account for: • - transient mechanical & temperature loads • - transient Weibull and fatigue parameters (temperature/time) • Interface codes that transfer transient analysis finite element • results into life prediction codes (CARES/Life)

5. CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures) Software For Designing With Brittle Material Structures • CARES/Life – Predicts the instantaneous and time-dependent probability of failure of advanced ceramic components under thermomechanical loading • Couples to ANSYS, ABAQUS, MARC, NASTRAN

6. CARES/Life Structure Finite Element Interface Output from FEA codes (stresses, temperatures, volumes) read and printed to Neutral Data Base Parameter Estimation Weibull and fatigue parameter estimates generated from failure data Reliability Evaluation Component reliability analysis determines “hot spots” and the risk of rupture intensity for each element

7. Some References Regarding Transient Reliability Analysis • Paluszny and Nicholls (1978) -- Discrete time steps, SCG, Weibull and fatigue parameters were constant: Paluszny, A., and Nicholls, P. F., “Predicting Time-Dependent Reliability of Ceramic Rotors,” Ceramics for High Performance Applications-II, edited by Burke, J., Lenoe, E., and Katz, N., Brook Hill, Chesnut Hill, Massachuestts, 1978. • Jakus and Ritter (1981) -- Probabilistic parameters for both applied stress (truncated Gaussian distribution) and component strength (Weibull distribution): Jakus, K., and Ritter, J, “Lifetime Prediction for Ceramics Under Random Loads,” Res Mechanica, vol. 2, pp. 39-52, 1981. • Stanley and Chau (1983) – Failure probability for non-monotonically increasing stresses (maximization procedure): Stanley, P., and Chau, F. S.; “A Probabilistic Treatment of Brittle Fracture Under Non-monotonically Increasing Stresses,” Int. J. of Frac., vol. 22, 1983, pp. 187-202. • Bruckner-Foit, A., and Ziegler (1999) – 3 Reliability formulations; no SCG, SCG governed by a power law,\ and SCG governed by a power law with a threshold: (1) Bruckner-Foit, A., and Ziegler, C., “Design Reliability and Lifetime Prediction of Ceramics,” Ceramics:Getting into the 2000’s, edited by Vincenzini, P., 1999. (2) Bruckner-Foit, A., and Ziegler, C., “Time-Dependent Reliability of Ceramic Components Subjected to High-Temperature Loading in a Corrosive Environment,” ASME paper number 99-GT-233, International Gas Turbine and Aeroengine Congress and Exhibition, Indianapolis, Indiana, 1999. • Ziegler (1998) -- SCG parameters vary with temperature/time: Ziegler, C., Bewertung der Zuverlassigkeit Keramischer Komponenten bei zeitlich veranderlichen Spannungen und bei Hochtemperaturbelastung, Ph.D. Thesis, Karlsruhe University, 1998. • Jadaan and Nemeth (2001) – Cyclic loading + Weibull and SCG parameters vary with temperature/time: (1) Jadaan, O, and Nemeth, N. N.;”Transient Reliability of Ceramic Structures.” Fatigue & Frac. Of Eng. Mater. Struct., vol. 24, pp. 475-487. (2) Nemeth, N. N., and Jadaan, O.; “Transient Reliability of Ceramic Structures For Heat Engine Applications,” Proceedings of the 5th Annual FAA/Air Force/NASA/Navy Workshop on the Application of Probabilistic Methods to Gas Turbine Engines, June 11-14, 2001, Westlake Ohio.

8. Transient Life Prediction TheoryFor Slow Crack Growth • Assumptions: • Component load and temperature history discretized into short time steps • Material properties, loads, and temperature assumed constant over each time step • Weibull and fatigue parameters allowed to vary over each time step – including Weibull modulus • Failure probability at the end of a time step and the beginning of the next time step are equal

9. Transient Life Prediction Theory -Slow Crack Growth and Cyclic Fatigue Crack Growth Laws Power Law: - Slow Crack Growth (SCG) Combined Power Law & Walker Law:SCG and Cyclic Fatigue - Denotes location and orientation

10. Transient Life Prediction Theory -Power Law General reliability formula for discrete time steps:

11. Binomial Series Approximation Used to Derive Computationally Efficient Solution For Cyclic Loading Binomial Series Expansion: When x>>y the series can be approximated as a two term expression (x + y)n xn + nxn-1 y , when x >> y

12. Computationally efficient transient reliability formula for cyclic loading - full solution

13. load time Transient Life Prediction Theory - Slow Crack Growth Modeled With Power Law Computationally efficient transient reliability formula for cyclic loading - simplified version T 2T ZT T

14. Combined Walker Law & Power Law for cyclic fatigue • Computationally • efficient version • with Z factor multiple

15. Temp m o N B Time step # 100 Time 5 230 Ieq 40 Temp 0.0021 500 1 25 9 226 100 36 100 0.021 2 1000 50 14 221 90 31 200 0.21 3 75 80 300 4 100 70 400 5 125 60 500 6 150 70 600 7 175 80 700 8 200 90 800 9 225 95 900 10 250 100 1000 Example – Tradeoff Between Accuracy and Computational Efficiency For a Cyclic Load 10 step transient uniaxial loading for a single load block – single element problem Temperature vs: material properties

16. Exact solution versus the Z approximation method for one solution increment (n = 1) • The results for one solution increment represent the least accurate but • most computationally efficient answer. Note: Load factor is 0.5 Pf = Failure Probability

17. Exact solution versus the Z approximation method for one solution increment (n = 1) • The results for one solution increment represent the least accurate but • most computationally efficient answer. Note: Load factor is 1.0 Pf = Failure Probability

18. Example of Z approximation method for various values of n. The solution increments are equally spaced (Zi = Zj= Zn). Percent error from exact solution versus number of load blocks for a failure probability prediction of 1000 cycles Increasing Computational Effort n = Number of discrete load blocks Pf = Failure Probability

19. EXAMPLE: Diesel Engine Si3N4 Exhaust Valve (ORNL/Detroit Diesel) OBJECTIVE: Contrast failure probability predictions for static loading Versus transient loading of a Diesel engine exhaust valve for the power law and a combined power & Walker law Material: Silicon Nitride NT551 Information Source:Andrews, M. A., Wereszczak, A. A., Kirkland, T. P., and Breder, K.; “Strength and Fatigue of NT551 Silicon Nitride and NT551 Diesel Exhaust Valves,” ORNL/TM1999/332. Available from the Oak Ridge National Laboratory 1999 Corum, J. M, Battiste, R. L., Gwaltney, R. C., and Luttrell, C. R.; “Design Analysis and Testing of Ceramic Exhaust Valve for Heavy Duty Diesel Engine,” ORNL/TM13253. Available from the Oak Ridge National Laboratory, 1996 DATA: MODEL: • ANSYS FEA analysis using axisymmetric elements • Combustion cycle (0.0315 sec.) discretized into 29 load steps • A 445 N (100 lb) spring pre-load applied to valve stem in • open position. 1335 N (300 lb) on valve stem on closure. • Thermal stresses superposed with mechanical stresses • Volume flaw failure assumed

20. Thermal distribution Loading and Stress Solution of Diesel Engine Exhaust Valve First principal stress at maximum applied pressure (MPa) Pressure load applied to face of a ceramic valve over the combustion cycle

21. Silicon Nitride NT551 Fast Fracture and SCG Material Properties Power Law Parameters (NT551): N and B Cyclic Fatigue Parameters: Q and A2A1 Note: Cyclic fatigue parameters are assumed values for demonstration purposes only

22. Diesel Engine Si3N4 Exhaust Valve Batdorf, SERR criterion with Griffith crack Transient and static probability of failure versus combustion cycles (1000 hrs = 1.14E+8 cycles)

23. Diesel Engine Si3N4 Exhaust Valve Transient reliability analysis with proof testing capability for combined Walker & power law

24. Diesel Engine Si3N4 Exhaust Valve Transient reliability analysis with proof testing capability Proof test: 10,000 cycles at 1.1 load level

25. EXAMPLE: Predict material reliability response of an alumina assuming time varying Weibull & Fatigue Parameters Material: Alumina Specimen: 4-pt flexure (2.2mm x 2.8mm x 50mm -- 38mm and 19mm bearing spans) Test Type: Static Fatigue Temperature: 10000 C Source: G. D. Quinn – J. Mat. Sci. – 1987 DATA: MODEL: • Single element model of specimen inner load span (2.8mm x 19mm) • with uniform uniaxial stress state (surface flaw analysis) • Loading is static (non-varying) over time • Weibull and fatigue parameters vary with the log of the time PROCEDURE: A single element CARES neutral file is constructed with discrete time steps (10 steps per decade on a log scale) spanning 8 orders of magnitude. Applied load is constant but Weibull and fatigue parameters allowed to vary with each time step.

26. EXAMPLE: Time Dependent Weibull & Fatigue Parameters G. D. Quinn, “Delayed Failure of a Commercial Vitreous Bonded Alumina”; J. of Mat. Sci., 22, 1987, pp 2309-2318. Static Fatigue Testing of Alumina (4-Point Flexure) 10000 C

27. Parameters interpolated with log of time - No extrapolation outside of range t = 1.6 sec., m = 29.4,s0= 165.8, N = 6.7, B = 2711.1 t = 31.6 sec., m = 15.8,s0= 152.7, N = 13.2, B = 9707.7 t = 1.0E+5 sec., m = 13.1,s0= 127.3, N = 36.4, B = 2276.2

28. Parameters interpolated with log of time - No extrapolation outside of range t = 1.6 sec., m = 29.4,s0= 165.8, N = 6.7, B = 2711.1 t = 31.6 sec., m = 7.4,s0= 263.3, N = 8.0, B = 2395.9 t = 316.2 sec., m = 4.5,s0= 870.1, N = 9.0, B = 10,389.0

29. Conclusions • A computationally efficient methodology for computing the transient reliability in ceramic components subjected to cyclic thermomechanical loading was developed for power law (SCG), and combined power & Walker law (SCG & cyclic fatigue). • This methodology accounts for varying stresses as well as varying Weibull and fatigue parameters with time/temperature. • FORTRAN routines have been coded for the CARES/Life (version 6.0), and examples demonstrating the program viability & capability were presented.

30. Future Plans • Goal to release CARES/Life 6.0 to engine companies for evaluation/beta testing by 9/30/02 - Continuing benchmarking activities - Continue developing GUI - Complete ANSYS and ABAQUS interfaces - User guide with example problems ?? (FY’02 - FY’03) • CARES/MEMS - Single crystal reliability - Edge recognition macro within ANSYS - Edge flaw reliability model