19.01.2014 Views

Structure and Mechanical Properties of Oxidation Resistant Alloys

Structure and Mechanical Properties of Oxidation Resistant Alloys

Structure and Mechanical Properties of Oxidation Resistant Alloys

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Structure and Mechanical

Properties of Oxidation

Resistant Alloys

Kevin Hemker

Johns Hopkins

hemker@jhu.edu


“The development of

ever more efficient gas

turbines has always been paced by the

results of research and development in

the concurrent fields of design and

materials technology.”

- G.W Goward

Surface and Coatings Tech, 108-109 (1998) 73-79


“Improved structural design

and cooling technology

applied to higher strengthat-temperature

alloys cast

by increasingly complex

methods, and coated with

steadily improved coating

systems, have lead to

remarkably efficient

turbine engines.”

- G.W Goward

Top coat

TGO

Bond coat

Substrate

A modern TBS

Surface and Coatings Tech, 108-109 (1998) 73-79


Drivers for bond coat design

Improve oxidation resistance and

slow TGO formation.

Increase bond coat strength.

Alloy to avoid phase transformations

and CTE mismatch with substrate.


Paradigm for materials development

Performance

- reduced weight

- increased temperature

Microstructure

- crystal structure

- grain size

- defects

- phases

Processing

- vapor deposition

- electrodeposition

-cast

- wrought

Properties

- modulus, strength, CTE

- temperature resistance


Outline

A brief history of oxidation resistant

alloys and coatings

Diffusion aluminide bond coats

• Microstructure

Properties

• Factors influencing ratcheting

Overlay bond coats

• Microstructure

Properties

• Factors influencing delamination


Delhi Iron Pillar

The famous iron pillar in Delhi is a

metallurgical wonder. This huge wrought iron

pillar, 24 feet in height 16.4 inches in

diameter at the bottom, and 6 1/2 tons in

weight has stood exposed to tropical sun and

rain for 1,600 years without showing any

visible signs of rusting or corrosion.

High concentrations of P catalyze the

formation Why? of δ-FeOOH (iron oxyhydroxide),

which is amorphous in nature and forms as

an adherent compact layer next to the

metal-scale How would interface, you Upon protect formation, the

corrosion resistance is significantly

enhanced this because iron pillor? δ-FeOOH forms a barrier

between the rust and the metal.”

- Balasubramaniam (IIT)

http://www.iupac.org/publications/ci/2005/2706/3_veazey.html#1


Paths to oxidation resistance

Alloying

• With select oxide formers

• Self healing

• Requires considerable quantities of solute

• Must not degrade other properties (e.g. creep)

Coatings

• Must be applied

• Must retard oxidation

• Must be adherent

• Cannot degrade properties of the substrate


Select oxide formers

Preferential oxidation

• Must form oxide or

protective layer before

base metal oxide

is formed

Protective

• Must form a dense, slow

growing, protective film

• E.g. Cr 2 O 3 , Al 2 O 3 , SiO 2

E for Oxide Formation (kJ/mol)


Alloying of steels

Low carbon steels

• Cheap, easily formed, and mechanically strong, but

rust at low T and oxidizes at high T.

Stainless steels

• Add significant quantities of Cr (~18%)

which forms a very protective oxide film

• Cuts down the rate of attack by 100 times

@ 900°C

• Affects microstructure and mechanical properties

• Other elements (e.g. Al, Si, P) also beneficial


Influence of Cr in superalloys

Metal loss (a la Ashby)

• Ni -> NiO

• Ni loss ~0.1mm in 600h at 0.7T m ~950°C

• Cr -> Cr 2 O 3

• Cr loss ~0.1mm in 1600h at 0.7T m ~1,231°C

• Cr loss ~0.1mm in 1,000,000h at 950°C

• Ni20% Cr -> Cr 2 O 3

• Ni20%Cr loss ~0.1mm in 6,000h at 950°C

– Foreign elements increase diffusion rates

• Commercial superalloys limited to


Alloying of Ni-base superalloys

Modern superalloys have been engineered

over decades and contain ~10 elements.

Composition of Single-crystalline René N5

Element: Ni Co Cr Ta Al W Re Mo Hf Ti C

Wt% 63.38 7.33 7.03 6.42 6.05 5.13 3.05 1.40 0.15 0.01 0.05

At% 64.92 7.48 8.13 2.13 13.48 1.68 0.99 0.88 0.05 0.01 0.25

The composition of most elements is

restricted by processing windows and

mechanical property requirements.


Diffusion aluminide coatings

Pack cementation aluminizing

• Parts packed in Al powder, sal ammoniac (NH 4 Cl),

graphite and alumina filler and heated to 450°C

for 2h.

• First used for Fe wire or ribbon heating elements,

Cu condenser tubes, etc.

First patent

• Van Aller, Allison, Hawkins (GE) 1911

• Way to ‘calorize’ metals to render then ‘inoxidizable’

• Attributed to selective formation of alumina scales


Schematic of the pack cementation

process


Pack cementation reactions

Pack Components

Source (Cr, Al, Si or their alloys)

Activator (NaCl, NH 4

Cl or other halide)

Inert Filler (often alumina)

Reactions between Source and Activator (Aluminizing)

3NaCl(g) + Al(l) = AlCl 3

(g) + 3Na(g)

2AlCl 3

(g) + Al(l) = 3AlCl 2

(g)

AlCl 2

(g) + Al(l) = 2AlCl(g)

Deposition on Substrate (Aluminizing)

Disproportionation 2AlCl(g) = Al + AlCl 2

(g)

3AlCl 2

(g) = Al + 2AlCl 3

(g)

Displacement

AlCl 2

(g) + Ni = Al + NiCl 2

(g)

Direct Reduction AlCl 3

(g) = Al + 3/2Cl 2

(g)

And, for activators which contain hydrogen e.g. NH 4

Cl,

Hydrogen Reduction AlCl 3

(g) + 3/2H 2

(g) = Al + 3HCl(g)

(The underlined symbols refer to species in the solid substrate.)


Pack cementation coatings

High-activity diffusion coating

on a Ni-base superalloy

- Al diffuses in -

Low-activity diffusion coating

on a Ni-base superalloy

- Ni diffuses out -


Schematic of CVD processing


Aero history of diffusion aluminides

1950’s

• Allison and Curtiss Wright – hot dipping of Ni blades

1960’s

• P&W – aluminizing of Ni blades by slurry processing

1970’s

• Most vane and blade coatings applied by “pack

cementation” and more recently by CVD

1990’s

• Recognized as useful TBC bond coats

• Forms an adherent α-Al 2 O 3 TGO


Making a diffusion aluminide

coating

NiAl

30-40 µm

Ni substrate

• Ni and Al diffuse together and form an

intermetallic coating


Ni-Al binary phase diagram:


Ni-Al binary phase diagram:

γ

FCC


Ni-Al binary phase diagram:

γ’

L1 2


Ni-Al binary phase diagram:

β

B2


Commercial Pt aluminide coatings

Pt

(Pt,Ni)Al + X i

30-40 µm

Superalloy substrate

• Bond coat chemistry and microstructure effected

by interdiffusion with substrate … reactive elements

appear to be good!!


TEM observations of diffusion

alumined coating

As-deposited

Thermal cycled

SADP

SADP

{110} T

100 nm

Micro-beam DP

B2

L1 0 Martensite


Compositional changes:

Concentration (at. %)

BC IDZ

60

40

(a) As-TBC'd

Ni

20

0

Al

Pt

60 (b) 28% FCT life

Ni

40

20

Al

0

Pt

60 (c) 100% FCT life

Ni

40

20

Al

0

Pt

0 20 40 60 80 100

Distance from surface of bond coat (µm)

β

β

M M +γ’ +

M +γ’ +

γ’+γ

γ’

γ’+γ

As-TBC’d

28% FCT life

γ’+γ 100% FCT life

Sample Ni Al Pt Cr Co Ta Re W

0% 36.72 43.50 12.10 3.9 3.78 0 0 0

28% 46.01 35.92 8.09 4.45 5.11 0.28 0.049 0.057

100% 51.32 34.44 4.08 4.21 5.17 0.57 0.08 0.13


Ni diffuses outward


Bond coat evolution …

Martensite transformation has been

reported in this composition range.

For example:

Smileck, et al, , 1973, Metall. . Trans.

Khadkikar, et al, , 1993, Metall. . Trans.


Failure by ratcheting and spallation

Key observations:

Isothermal exposure is not

the same as thermal cycling.

TGO thickens, lengthens and

roughens as a result of

thermal cycling.

Bond coat must creep to allow

TGO to ratchet.

Top coat cannot creep;

it cracks instead

Mumm, Evans et al. (2001)


Advanced model of TGO ratcheting:

Balint-Hutchinson model

A special purpose multilayer code that takes input from

micromechanical parameters and predicts the progressive

development of undulations of the TGO layer upon thermal cycling.

Key microstructural parameters:

E sub , α sub , E bc , α bc , YS bc , (creep rate) bc , t bc , ε mart , M s , A s ,

α tgo , E tgo , (growth rate) tgo , (creep rate) tgo , YS tgo , t tgo ,

(shape) tgo , E tbc , α tbc , (creep rate) tbc , YS tbc , …

Note: property measurements and microstructural

observations of bond coat are crucial to model

development.

Balint and Hutchinson, Acta 2003, JMPS 2005


Measuring the mechanical properties

of diffusion aluminide bond coats


Bond coat microsample preparation:

Bond coat

Substrate

600 µm slice is scalped from 1 inch diameter button using

a wire EDM.

A lapping machine is employed to grind the superalloy side

of the slice to get a uniform thickness of ~150µm.

500 µm

A sinking EDM is then used to punch microsamples from

each slice.

Faces of the microsample are polished to a mirror finish

and a final thickness of 30~50 µm with TEM ‘dimpler’.

100µm

Pt lines are deposited on both sides of the microsample

using a Focused Ion Beam (FIB).


Microsample testing

• Microsamples:

3.5 mm x 1.2 mm x 25 µm

(overall)

250 µm x 500 µm x 25 µm

(in gage)

Fringe

Detector

Laser

Fringe

Detector

Linear Air

Bearing

Load Cell

Piezoelectric

actuator

Grips

Manual

turnbuckle

• Temperature range of 25-1200°C can be

achieved by resistive heating.

• ISDG provides a strain resolution of

better than 5 µstrain.


Measuring bond coat CTE (T):

Engineering Strain (%)

1.5

1

0.5

y = -0.15515 + 0.0015501x R 2 = 0.98857

As-bond-coated

0

400 500 600 700 800 900

Temperature ( o C)

σ = ~2 MPa

• CTE of as-coated

bond coat has been

determined to be

15.5 ppm o C -1

between 400 o C and

850 o C.

•CTE measurement

above 850 o C was

hindered by creep.


Optical CTE measurements

Rene N5


xx

( T )

CTE =

dT

=

d

dT




∂u

∂x




Optical CTE measurements

?


Measuring bond coat E(T):

400

200

350

Stress (MPa)

300

250

200

150

100

50

0

0.2

T = 25 o C

x

E

Strain (%)

T = 650 o C

Young's Modulus (GPa)

150

100

50

BC-1, as-coated

BC-2, as-coated

BC-1, cycled

NiAl literature values

0

0 200 400 600 800 1000

Temperature ( o C)


Microsample tests of Rene N5 :

1200

Engineering Stress (MPa)

1000

800

600

400

200

0

E = 129 GPa

E = 126 GPa

E = 126 GPa

T = 25 o C

0 0.01 0.02 0.03 0.04

Engineering Strain

• Modulus of [100] oriented

single-crystalline Rene N5

microsamples measured to be

126 ± 3 GPa (8 measurements).

• GE reported modulus to be

E = 129 GPa.


Texture effects on E:

Young's Modulus

300

250

200

150

100

The Elastic Constants of NiAl

(Ref. Rusovic et al, 1977)

C 11 C 12 C 44

199 GPa 137 GPa 116 GPa

50

0 50 100 150

In-plane Angle (Degree)

(001) plane

(011) plane

(111) plane

The Average E in Certain Texture of NiAl

E(GPa) (001)

plane

(011)

plane

(111)

plane

Voigt 125 184 181

Reuss 117 159 181

• Each angle corresponds to one certain

orientation in that plane.

• The measured value of E is comparable to

that of NiAl with a [001] texture.

• X-ray texture determination did not

evidence any obvious out-of-plane texture

of bond coats.


Measuring bond coat σ yield (T):

300

700

Engineering Stress (MPa)

250

200

150

100

50

E

T = 650 o C

T = 890 o C

T = 980 o C

Yielding Strength (MPa)

600

500

400

300

200

100

As-BC'd -1

As-BC'd -2

TC100

0

0.2

Engineering Strain(%)

0

600 700 800 900 1000 1100 1200

Temperature ( o C)


Effect of thermal cycling

600

BC1,as-coated

Yield Strength (MPa)

500

400

300

200

100

BC1, cycled

BC2, as-coated

BC2, cycled

BC3, as-coated

BC3, cycled

Deng Pan et al., Acta 2003

0

600 700 800 900 1000 1100 1200

Temperature (C)


Effect of thermal cycling

600

BC1,as-coated

500

BC1, cycled

BC2, as-coated

Yield Strength (MPa)

L1 0

Martensite

400

300

200

100

BC2, cycled

BC3, as-coated

BC3, cycled

100 nm

B2

Deng Pan et al., Acta 2003

0

600 700 800 900 1000 1100 1200

Temperature (C)


High T deformation


Stress relaxation behavior:

Stress (MPa)

300

250

200

150

BC-1, as-coated

400 o C -experimental

650 o C -experimental

900 o C -experimental

1150 o C -experimental

100

50

0

0 200 400 600 800 1000

Time (s)

Strain rate*exp(Q/RT) (1/s)

.

ε pl =A(σ/E sample ) n exp(-Q/RT)

10 4

10 10 10 3

650 o C

10 82

n = 4

10 1

10 6 n = 4

900 o C

10 0

10 4

10 -1 1150 o C

Q = 250kJ/mol 125kJ/mol

10 -2 2

10 -4 10 -3 -3 10 -2 -2

•Empirical fit suggests that Q creep is only half of Q D for NiAl and n = 4.

σ/E


Stress relaxation behavior

300

250

as-deposited

Stress (MPa)

200

150

100

400 o C -experimental

650 o C -experimental

900 o C -experimental

1150 o C -experimental

Power law

50

Empirical fit:

. .

0

0 200 400 600 800 1000

Time (s)

ε pl =-σ/E machine = 1.2 x 10 14 (σ/E sample

) 4 exp(-125(kJ/mol)/RT)


Characterizing the

microstructural changes in

diffusion aluminide bond coats


In-situ TEM observations:

25°C -> 900°C

L1 0

B2

Transformation from martensitic L1 0 to B2 observed

at approximately 800°C.


n situ heating X-ray diffraction:

Reversible transformation

Crystallographic relationship

c 25°C, cooled from 1150°C

between B 2 and L1 0

b

1150°C

Intensity

25°C

(100)

(110) B2

L1o

(110) L12

(111) L12

(111) L1o

(110) B2

(200) L1o

(200) L12

(200) L12

(111) B2

(201) L1o

(200) B2

(220) L1o

(210) B2

(202) L12

(221) L1o

(211) B2

a

L1 B 2

0

(011)B2//(11-1)L1 0

(311) L1o

30 35 40 45 50 55 60 65 70 75 80 85 90

2θ (degree)


Transformation to B2:

Intensity (counts)

7500

7000

6500

6000

5500

5000

4500

4000

3500

3000

2500

2000

1500

1000

500

(111)L10

(110)B2

0

-500

41.0 41.5 42.0 42.5 43.0 43.5 44.0 44.5 45.0 45.5 46.0 46.5 47.0 47.5 48.0 48.5

2θ (degree)

450 ºC

500 ºC

550 ºC

600 ºC

625 ºC

650 ºC

675 ºC

685 ºC

700 ºC

(200)L10

Volume Ratio(B2/ M)

60

50

40

30

20

10

0

Volume ratio between B2 and

Martensite calculated from

the intensities of (110)B2 and

(111)L10

450 500 550 600 650 700

Temperature (ºC)

As

The A s temperature is approximately 620ºC.


Martensite transformation temperatures:

Heat Flow Endo Down (mW)

80

60

40

20

0

-20

-40

DTA (Differential Thermal Analyzer)

Peak T of B2 to M

Cooling at ~25 ° C/min

Peak T of M to B2

Heating at 40C/min

°

0 200 400 600 800 1000 1200

‣The M to B2

transformation occurs

at 600-700ºC

‣The B2 to M

transformation occurs

at 550-650 ºC

‣M s is approximately

40-50ºC less than A s

Temperature(C)


Thermal and transformation strains:

2.2

2.0

1.8

Thermal strain of B2

1.6

1.4

Transformation

M s

Strain (%)

1.2

1.0

0.8

0.6

strain

cooling

A s

heating

0.4

0.2

Thermal strain of L1 0

(° C)

0.0

0 200 400 600 800 1000

Temperature


Effect of martensitic transformation

Balint and Hutchinson, JMPS 2005


Strategies for bond coat design

Improve oxidation resistance and

slow TGO formation.

Increase bond coat strength.

Alloy to avoid martensitic and

other phase transformations.


Overlay coatings for superalloys

Developed to allow for greater

flexibility in bond coat compositions.

• Addition of reactive elements known to be

beneficial

• Peg formation, S gettering, increased TGO adhesion ???

• Can tailor the phases present in the bond coat

• β + γ, β + γ ’, γ + γ ’, etc.

Requires deposition

• Initially PVD

• More commonly now LPPS


Schematic of an EBPVD process


EBPVD coatings

CoCrAlY coating as-deposited

CoCrAlY coating after peening

and heat treatment


Schematic of a plasma deposition

system


NiCoCrAlY argon-shrouded

plasma sprayed coating

NiCoCrAlY coating after peening and heat treatment


Aero history of overlay coatings

FeCrAlY

• Alloy merged from GE nuclear programs (1964)

• P&W applied (EBPVD) as a coating (1970)

• Limited at high T by formation of a brittle NiAl layer

CoCrAlY (1972)

• Useful oxidation resistance but too brittle

NiCrAlY (1973)

• Limited hot corrosion resistance but ductile

NiCoCrAlY (1975; 1986)

• Good compromise; 40+ patented MCrAlY’s


Microstructural observations of

Pratt & Whitney’s “two-phase”

NiCoCrAlY(Hf,Si) bond coat


Effect of thermal cycling

As-received 100 hrs @ 1100 o C

Near-TGO γ

100

80

60

100 hrs @ T i

Average across entire thickness

γ

Vf [%]

40

20

β

0

0 20 40 900 1000 1100 1200

Temperature [C]

“two-phase bond coat”

(β-NiAl + γ-Ni)


TEM observations of the bond coat

1750

1500

β-phase

Ni Ka

1250

Counts

1000

750

Co L/ Ni L

500

Al Ka

250

0

Co Ka

Cu Ka

Cr Ka

Ni Kb

Cr Kb

Cu Kb

0 1 2 3 4 5 6 7 8 9 10

Energy (keV)

1750

1500

γ-phase

Ni Ka

1 µm

Counts

1250

1000

750

Cr Ka

Co Ka

As–received NiCoCrAlY

500

250

0

Co L/ Ni L

Cr L

C K

Al K

Si K

Cr Kb

0 2 4 6 8 10

Energy (keV)

Cu Ka

Co Kb

Cu Kb


TEM observations of the γ grains:

As-received

BF

1 10

DF

110

100 nm

200 nm

ELTB

BF

500 nm

DF

1 10

110

γ grains filled with

very fine γ’

precipitates

modest coarsening

suggests the γ’

dissolves at T.


EFTEM of γ′ in the γ grains:

200nm

Elastic image Thickness map Cr L (574 eV) Co L (779 eV)

The γ’ precipitates are

rich in Ni and Al while

the matrix is rich in Cr

and Co.

Ni L (854 eV)

Al L (73 eV)


TEM observations of the β grains

As-received

ELTB

100 hrs @ 1100 o C

750 nm

g 110

g 110

250 nm

-Mottled-

250 nm

-Tweed-

-Lath-

[001] [001]


Do β-grains contain martensite?

As-received

M s


Microstructural evolution of NiCoCrAlY

Beta Gamma

As-received

Thermal cycled


Precipitates found in an ELTB

8000

6000

Cr

Counts

4000

?

2000

Energy (keV)

Co Ka

Ni Ka

Cr

Cr L Co Kb

Co Kb

Si Ka

Ni

0

0 1 2 3 4 5 6 7 8 9 10

500 nm

Cr-rich intergranular precipitates

widely dispersed in the ELTB

specimen but not observed in the

as–received samples.

Composition is 55.5 at% Cr, 32.1 %

Co, 11.0 % Ni and 1.4 % Si.


Crystal structure of Cr-rich precipitate

110/-1-10 110/-1-10

10 o (10 o ) 10 o (7 o )

B = [-110] B = [-221]

B = [-111]

30 o (31 o )

001/00-1

413/-4-1-3

B = [-140]

The precipitate is a (CrCoNi)

σ-phase (P4 2 /mnm space group).

Tetraganol unit cell with lattice parameters

a = 8.37 Å and c = 4.30 Å.


Elemental maps of the σ-precipitates

Mendis

Backscattered image

Al map

Cr map

Tryon

Backscattered image Cr map Si map


Potential impact of σ-phase …

Notoriously brittle but bond coat is not a load

carrying structure.

Appears to trap Si, which would mitigate its effect.

Is observed at the TGO interface and may affect

adhesion.

Bond coat

TBC

TGO

Elastic image

Cr L (574 eV)

Note: only observed in ELTB;

not as-received or after 100hrs at 1100C!


Predicting phase stability…

Achar, Munoz, Arroyo, Singheiser, Quadakkers, Modeling of phase

equilibria in MCrAlY coating systems, Surface and Coating Tech (2004)


TGO and PEG

Chemistry

3500

3000

2500

O/ Cr

Al

EDX of TGO peg

columnar TGO

Counts

2000

1500

1000

500

0

Y

Co

Co/ Ni

Ni

Hf

Cr

Hf/ Cu

S

Cr

Ni Hf

0 2 4 6 8 10

Energy (kV)

1 µm

4000

3500

Al

Counts

3000

2500

2000

1500

O

EDX of columnar TGO

peg

1000

500

0

Y

0 2 4 6 8 10

Energy (kV)

The peg contains a large number of elements

including sulfur, although the neighboring columnar

TGO is nearly pure Al 2 O 3 .


STEM shows S in pegs but not

the TGO

Bond coat

Al

O

Ni

Co

S

Hf

Y

Cr


Mechanical behavior of P&W’s

NiCoCrAlY(Hf,Si) bond coat


RT strength of NiCoCrAlY

1600

Tensile test of NiCoCrAlY at RT

Eng. Stress [MPa]

1400

1200

1000

800

600

400

E=155 GPa

as received

thermally cycled

E=130 GPa

1.4 GPa

with RT

ductility!

200

E=130 GPa

0

0.0 0.5 1.0 1.5 2.0 2.5

True Strain [%]


Effect of thermal cycling

Diffusion aluminde

NiCoCrAlY

- strength increases - strength decreases

Yield Strength (MPa)

600

500

400

300

200

100

Deng Pan et al., Acta 2003

BC1,as-coated

BC1, cycled

BC2, as-coated

BC2, cycled

BC3, as-coated

BC3, cycled

Eng. Stress [MPa]

1600

1400

1200

1000

800

600

400

200

E=155 GPa

E=130 GPa

as received

thermally cycled

E=130 GPa

0

600 700 800 900 1000 1100 1200

Temperature (C)

0

0.0 0.5 1.0 1.5 2.0 2.5

True Strain [%]


High temperature deformation:

Strain [%]

12

10

8

6

4

2

0

250

800

Stress [MPa]

200

150

100

50

0

400 450 500 550 600 650 700

700

600

500

400

300

Temperature [C]

T im e [s ]


High temperature deformation:

Eng. Stress [MPa]

Constant strainrate Constant strainrate

225

Cross head fixed Cross head fixed

200

175

150

125

100

75

3.0

2.5

2.0

1.5

1.0

True Strain [%]

50

25

0

0.5

0.0

0 50 100 150 200 250

Time [s]

500ºC


Comparison of high T strength:

500

Flow stress [MPa]

400

300

200

100

MCrAlY

(Ni,Pt)Al corr

(Ni,Pt)Al

0

400 600 800 1000 1200

Temperature [C]


Failure of NiCoCrAlY TBC systems

e.g. NiCoCrAlY bond coat systems

Salient observations:

No TGO ratcheting.

Cracks occur below the

TGO.

Crack growth is mode-

II and effected by

interfacial energy and

plasticity in the

adjacent layers.

Cracks must be initiated

before they can grow.

Xu, Evans et al., Acta 2004


Proposed failure mechanism map:

No !!

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!