Document not found! Please try again

Using OOF to Model Mechanical Behavior of Thermal ... - ctcms - NIST

2 downloads 0 Views 4MB Size Report
Andrew C. E. Reid & Andrew R. Roosen,. CTCMS, MSEL, NIST. • Andrew J. Allen, Jan Ilavsky, Jay S. Wallace,. Ceramics Div., NIST .... Modulus (G. Pa). 0. 20. 40.
Using OOF to Model Mechanical Behavior of Thermal Barrier Coatings Mark R. Locatelli and Edwin R. Fuller, Jr. National Institute of Standards and Technology Gaithersburg , MD 20899-8521, U.S.A. Symposium on Advanced Ceramic Coatings Intl. Conf. on Advanced Ceramics & Composites Cocoa Beach, FL - January 15, 2002

Acknowledgments • Stephen A. Langer, Information Tech. Lab., NIST • Andrew C. E. Reid & Andrew R. Roosen, CTCMS, MSEL, NIST • Andrew J. Allen, Jan Ilavsky, Jay S. Wallace, Ceramics Div., NIST • Scott Terry & Carlos Levi, Univ. of Calif. at Santa Barbara

Technical Issues for TBC’s • Correlate properties with microstructure Ø Ø Ø

to shorten materials development cycle to improve materials & processing to enable more reliable design

• Increase thermal protection • Increase life • Increase reliability,

i.e., predict life, or coating spallation APPROACH: Develop computational tools for elucidating influences of stochastic microstructural features (e.g., porosity) on physical properties; and provide insights into mechanisms that lead to TBC spallation via predictive micro-mechanical models of reliability.

PPM2OOF Tool

Micrograph Binary Image

§ § § §

Mesh

Convert micrograph to “.ppm” (portable pixel map) file Select & identify phases to create binary image Assign constitutive physical properties to each phase Mesh in PPM2OOF via “Simple Mesh” or “Adaptive Mesh” – multiple algorithms that allow elements to adapt to the microstructure

OOF Tool Virtual Experiments: constrained cooling

Visualize & Quantify: normal residual stresses

dT

Perform virtual experiments on finite-element mesh: § To determine effective macroscopic properties § To elucidate parametric influences § To visualize microstructural physics

Young’s Modulus of EB-PVD TBC’s Intercolumnar Gap Porosity: essentially constant through TBC thickness via area fraction measurements

5 µm

t = 60 µm

Ts = 1100 °C 8 RPM 0.9 µm/min

t = 40 µm 66 mm

t = 20 µm

TBC

0.75 ° t = 1-2 µm

FeCrAlY t = 0 µm Courtesy of Scott Terry & Carlos Levi, University of California, Santa Barbara

EB-PVD TBC Cross Sections: 11716

11708

E÷÷ = 126 GPa

E÷÷ = 173 GPa E^ = 156 GPa 1 mm

E^ = 59 GPa 1 mm

Near interface microstructure

Surface microstructure

1.3 m m from TGO/TBC interface

60 m m from TGO/TBC interface

Courtesy of Scott Terry & Carlos Levi, University of California, Santa Barbara

Porosity: key key to to Porosity: TBC Performance Performance TBC Intercolumnar î Compliance

Intracolumnar î Insulation

200 nm

Different scales of shadowing

10 µm

100 nm

TEM courtesy of E. Sommer and M. Rühle MPI-Stuttgart

Courtesy of Scott Terry & Carlos Levi, University of California, Santa Barbara

Young’s Modulus of EB-PVD TBC’s Young’s Modulus (GPa)

Modulus is anisotropic and position dependent E÷÷

E^ Ts = 1100 °C

Height above TGO (mm)

Effect of of Temperature Temperature on on Surface Surface Effect Morphology and and Texture Texture Morphology Rotating Substrates

10 µm

Ts = 1100 °C Ts = 1000 °C

8 RPM, ~0.9 µm/min (Flux ~ 3.2 µm/min)

Ts = 900 °C Courtesy of Scott Terry & Carlos Levi, University of California, Santa Barbara

Young’s Modulus of EB-PVD TBC’s Young’s Modulus (GPa)

Modulus is anisotropic and position dependent E÷÷

˜ ¢  £

Ts = 900 °C Ts = 1100 °C

E^

Height above TGO (mm)

Section View of a ZrO2 – 8 wt% Y2O3 Plasma Sprayed Thermal Barrier Coating Free Standing Monolith more than 50 plasma spray parameters! ZrO2 – MgO

50 mm

AM100

Fabricated at Thermal Spray Laboratory, SUNY, Stony Brook (Jan Ilavsky)

Young’s Modulus versus Porosity 220

ZrO 2 - 8 w/o Y 2O3 (sintered) 2.381 mm f WC Indenter 4.8 N Load (Evaluated at 4 N)

200

Literature Value

Modulus [GPa]

180 160

Sintered Samples

140 120 annealed, 1 h

100

T

80

cross section

60 40

Amdry 142: 115 mm spray distance

20

intralamellar microcracks

plan section

0 80

85

90

95

100

Density [% rth] J. S. Wallace and J. Ilavsky, J. of Thermal Spray Tech., 7, [4], 521-526 (1998).

Calculating Average Elastic Properties of a Representative Region

85 mm

Effective Elastic Young’s Modulus Calculated From Microstructural Finite Elements Porosity Section Moduli (%) (GPa) Ex Calc.:

12 ± 1

Expt.:

11

Porosity (%)

Ey

39 ± 10 13 ± 7

6

28

11

Plan Moduli (GPa) Ex

Ey

83

90 19

ZrO2 – 8 wt% Y2O3: E = 214 GPa n = 0.310

How to Generate Effective Meshes?

Original Image

Low Res. Mesh Ex = 136 GPa Ey = 100 GPa Nodes = 10,894

High Res. Mesh Ex = 83 GPa Ey = 40 GPa Nodes = 43,887

Modulus versus Mesh Resolution (

160

)

AMLH (Ex) AMLH (Ey) AMTTA (Ex) AMTTA (Ey) SXLD (Ex) SXLD (Ey) SXTTC (Ex) SXTTC (Ey)

140

Modulus (GPa)

120

100

80

60

40

20

0 0

2 10

4

4 10

4

6 10

4

8 10

4

Number of Nodes

1 10

5

1.2 10

5

Modulus versus Mesh Resolution (

160

140

If the structure is truly random, then errors in simulation might vary with the spacing between nodes:



1 N

Modulus (GPa)

How accurate should mesh be?

120

100

)

AMLH (Ex) AMLH (Ey) AMTTA (Ex) AMTTA (Ey) SXLD (Ex) SXLD (Ey) SXTTC (Ex) SXTTC (Ey)

80

60

40

????

20

0 0

0.002

0.004

0.006 1/Root(N)

0.008

0.01

0.012

Critical Issues § Does the mesh resolution capture the

essential features that affect behavior? E.g., Does the mesh capture the essence of the fine cracks?

§ If the mesh resolution is changed, does the simulated behavior also change?

§ This there asymptotic behavior? § Can these techniques be validated in a general manner?

Modeling Mechanical Behavior of TBC’s SUMMARY: § Microstructure-based, finite-element simulations

§ § §

provide a new paradigm for property measurements of complex materials, such as, TBC’s. Sample preparation & image analysis are critical for obtaining accurate, quantitative measures of behavior. Mesh resolution can have significant influences on determined properties. Finite-element simulations help to elucidate the influences of stochastic microstructural features (e.g., porosity) on the elastic behavior of complex TBC microstructures.