On secondary creep of anisotropic nuclear materials

On secondary creep of anisotropic nuclear materials

JOURNAL OF NUCLEAR 33 MATERIALS ON SECONDARY (1969) 52-63. CREEP OF ANISOTROPIC of Waterloo, Received Department 26 November anisotropic p...

1MB Sizes 2 Downloads 70 Views

JOURNAL

OF NUCLEAR

33

MATERIALS

ON SECONDARY

(1969)

52-63.

CREEP OF ANISOTROPIC

of Waterloo,

Received

Department

26 November

anisotropic pressure tubes may be twice

half

as high

as diametral

tubes if the uniaxial about

the third

creep rates

power

test results available in-reactor

creep

is shown

at

can

least

Using

worked

or one isotropic

and diametral

that

specimens

in accordance

to

the limited

zirconium

qualitatively

alloys

it

experimental

sont

En

disponibles,

5t la

de

moins

zirconium

with

the

sur Bchantillons des tubes accord

avec le theor&me

anisotrope

qui suppose

de puissance contrainte.

Diese

dans

de fluage en rkacteur,

sol1

an

en contrainte

la

aus

anisotropen

metralen

uniaxiale,

wenn

des

theorbme

de

moduln

etwa

20%

Bcroui.

Un

pour des alliages

dans les modules

de fluage d’environ

zirconium

peuvent

les vitesses

de fluage dans le sens diametral

der sich

Btre 2 St 1,5 fois aussi Blevbes que

designed alloy

to

materials serve

pressure

in

tubes

are the

as

form core

of

analytic *

theory

zirconium

components

Research

Assistant

on leave

creep from

in

for multiBuilding

bzw.

Spannungszustand

Reaktor-

kijnnen bei

die dia-

anisotropen

Unterschiede

betragen

und

der

die

Kriech-

einachsigen

etwa mit der dritten Anhand

qualitativ

Zusammenhiinge lassen.

ge-

zum Kriechen unter

Danach

bzw.

Zirkon

halb so gross wie bei isotropen die

variieren.

zumindest

Kriech-

Bestrahlung

Weitere

mit

der

Potenz

vorhandenen

nachgewiesen, dem

Theorem

Untersuchungen

dass be-

werden

axial states of stress. Berman and Pai 7) have formulated a simplified theorem of creep for anisotropic materials assuming that in stress space the magnitude but not the direction of the creep rate vector depends on anisotropy. Here an exact analysis of anisotropic in-reactor creep is proposed for zirconium alloys considering the work of Kachanov *), Odqvist 9), Rabotnov lo), Malinin I’), Olszak and Sawczuk

and

l-4). Recently,

of anisotropic

la

vorgeschlagen.

Ross-Ross and Hunt 5) and Ibrahim 6) studied the inreactor creep behaviour of these pressure tubes experimentally and indicated the need for an nuclear power reactors

wird die

schreiben

des tubes

fabricated

Spannung

Daten

de force

Introduction Anisotropic

avec

der

kalt-bearbeitetem

Kriechgeschwindigkeiten

20%, les vitesses des tubes

anisotrope

1.

aniso-

que pour des diffhrences

dans le sens diam&ral

fluage

de ces travanx

wie sie unter

Kriechgeschwindigkeiten sein,

de fluage

anisotropes

Zr-Legierungen

doppelt

qui indique

le

du coefficient

Diskrepanzen

vorgeschlagen.

Rohren

en reacteur

pour

Es wird ein Theorem

Rohren

trope est propose

et sur eux en

ult&ieur

mit einachsigem

Druckrohren

bedingungen

mesur&

die

par

fluage

uniaxial

propose

klgren,

constitue

zirconium

qn’au

i?tre relies entre

ou des coefficients

et sur des tubes de force d’un reacteur de

montre

expkrimentaux

une variation

Un dt?veloppement

Arbeit

von

alliages

il est

& un fluage

peuvent

condition

dans un rkacteur

ces r&ultats

soumis

de force,

en

limit&

est sugg&&

theorem

les d&accords

it la fois sur des Bchantillons

Bcroui,

qualitativement

3 de

d’essais

le fluage

d’alliage

messen wurden.

des vitesses

pour

et dans le sens diam&ral

bei Proben

corr8lation

fois

la puissance

&ultats

and pressure tubes

for further work are indicated.

d’expliquer

les

uniaxiale

coefficients

tente

Canada

8. environ

utilisant

or of the anisotropy

Ce travail

Ontario,

proportionnelles

l’effort.

geschwindigkeiten

with stress. Suggestions

*

isotropes si les vitesses de fluage en condition uniaxiale

a variation

creep

MATERIALS

Waterloo,

assuming

for anisotropic

of the power-index

of

are proportional

for both uniaxial

of cold

be correlated

rates

of stress.

results of both uniaxial proposed

creep

CO., AMSTERDAM

1968; in revised form 9 April 1969

This work attempts to explain discrepancies in the correlation of in-reactor creep rates measured on both uniaxial specimens 14) and reactor pressure tubes 5, 6) made of cold worked zirconium alloys. An in-reactor creep theorem for anisotropic zirconium alloys is proposed which indicates that for differences in creep moduli of approximately 20’36, diametral creep rates of

Engineering,

of Civil

PUBLISHING

NUCLEAR

and M. J. HOLICKY

J. SCHROEDER University

0 NORTH-HOLLAND

Research 52

Institute,

GWT,

Prague,

Czechoslovakia.

ON

is) Freudenthal

SECONDARY

and

CREEP

OF

ANISOTROPIC

53

MATERIALS

B may be calculated in view of eqs. (1) and (2)

Geiringer is). The work-

ability of the proposed analysis is demonstrated

NUCLEAR

from the condition

by an attempt to correlate diametral creep rates

that

B~jm&&m =

/%” = 8,

(5)

measured in the reactor on cold worked zirconium alloy pressure tubes 596) with uniaxial

where

creep data obtained by Fidleris 14). An extensive

B and s are good measures of creep behaviour

review

of zirconium

given 2.

of

literature

on

anisotropic

creep

is

in ref. ‘9 12).

through In the following,

Cartesian tensors are used

term. &f denotes the Kronecker delta. Considering the work *-Is) the following theorem

is

proposed

iij = /9s(“-1)Asjrlskl, &t = 0, dff represents

secondary

for

(1)

steady

creep

rates, both B and m are material constants depending on temperature T and/or fast neutron flux 9j and may be found from uniaxial creep test, the Aijkl are coefhcients of anisotropy, .Q = GIN--&ah&~ designates stress tensor aij= ajr and

the deviator

8 = (A~,xIw&, is the so-called Theorem

equivalent creep

of the

(2) or effective

(1) may be developed

a dissipative

alloys,

depends

potential

stress.

such

that

,5W = J.$ daa/ represents the complementary of the dissipative strain energy density dusij = &a~~ and consequently Baj=/lo W/bsij. Choosing W = =sm+l/(m+ 1) the power law (1) is obtained. The incompressibility condition dii = 0 may be enforced by setting Ati&

= 0.

(3)

An equivalent or effective strain rate 6 is defined by postulating 16) that gi+si3= is and it follows from eqs. (1) and (2) that the theorem (1) can be rewritten in the form i = @n.

(4)

upon

how well

in terms of test data

of this work, theorems

(1)

or (4) will be restricted to cases where the axes of anisotropy are orthogonal and approximately coincident with the principal axes x1, x2, 5s of strain rate and stress. In addition, creep rate differences in tension and compression are neglected even though theorem (1) may be generalized to account for these differences by defining All11 = A$,, if a11 is tensile, Aiiaa= if ai1 is compressive, Aliz,= A&& if = A,,, all and aa2 are tensile, Anss= A&; if ai1 is tensile and a22 is compressive, etc. Under the restrictions specified above, and in view of the symmetry of both the strain rate tensor and stress tensor, problem (1) involves six coefficients of anisotropy AIIII, A22227 A33339 A1122, 4133, which are interrelated by the three constraints shown in (3). Choosing AIIII, Az222, An22 as independent, s and 6 may be expressed for these simplified conditions as A2233

s = [A1111(~11- ~33)~

by introducing W

Whether or not

eq. (4) or eq. (1).

In the remainder

as described for example l5), i.e., the orthogonal Cartesian coordinates are designated by xl, x2, x3 and the subscripts i, i, k, 1, m, n, p, q= 1, 2, 3 are summation indices if repeated in a single

where

= Amnpp

i and s can be correlated

A creep theorem for zirconium alloys

anisotropic creep zirconium alloys :

&d~~rnnAklpq

+ A2222(S22 - ~33)~ +

+ 2A1122(S33 - Sll)(S33

- S22)]*

(5)

and d= (All1lA2222 - A2~~2$t(A,,,,~2,, + 3.

A1111i222 -

+

%41122il1i22)t

(6)

Uniaxial and biaxial in-reactor creep of cold-worked zirconium-alloys

As part of a major test programme, Fidleris 14) investigated the in-reactor creep behaviour of uniaxial specimens manufactured from coldworked zirconium alloys. The specimens were made from bar stock or machined from pressure tubes such that the longitudinal axes of the specimens were parallel to the axes of the tubes.

J. SCHROEDER

54

For

bot’h types

{ lOi%}

plane

of

specimens

normals

the

were

parallel to the longitudinal

AND 111.J. HOLICKY

prismatic

predominantly

uniaxial data 14) but using normalization proposed 5).

axis of the specimen.

Zirconium

the

linear

flux

alloy pressure tubes investigated

Ross-Ross

and Hunt 5) and Ibrahim 6) studied

by

in-reactor

creep of biaxially

exhibited a t,hickness/radius ratio 0.05 -
tubes made from cold-worked

stressed pressure zirconium

alloys

ROSS-ROSS and

and have

to be classified

of the t’ype studied 14). Results

of all t’hree

r and h denote

investigations

of creep rate

and

indicate dependence

on st’ress, temperature

and fast neutron

flux

over the range tested. [In 14, 6) test data of heat treated zirconium alloys are also report’ed but because of the scatter of results encountered in the case of heat-treated uniaxial creep tests, no correlation is attempted.] In 5) in-reactor diametral creep rates of pressure tubes tested were correlated employing an empirical creep law and an attempt was made to interrelate the uniaxial results 14) and biaxial test data of the pressure tubes 5) using an isotropic analysis as an approximation since a manageable theorem of creep for anisotropic materials was not available as remarked 5). Using the isotropic approach, diametral creep rates calculated in accordance with uniaxial test data appeared to be in agreement wit’h t’est results of tubes at hoop stresses of 11 000 psi only, while at 20 000 psi t,he calculat’ed creep rates were three to six times the rat’es indicated by the tubes tested. Uniaxinl creep rat,es 14) and diametral

thickness

complicated

Hunt 5) and as thick

Ibrahim 6) shells 18);

radius of the middle of

the

tubes

surface

respectively.

A

thick shell analysis, as for example

demonstrated 9, 19), is in this case not warranted since it may improve the results by only 5o;b I*). Pressure tubes test,ed 5) passed vertically through the reactor core a,nd contained a t,cn foot high stack of six fuel assemblies cooled by pressurized water. Each fuel assembly had a different power output which decreased the neutron flux in steps from the center towards both ends of the tubes. Due to the variation of creep rates with fast neutron flux the deformation mode of axial sections of the tubes consisted of a series of bell-shaped stel)s associated with every fuel assembly as shown in fig. 3 of ref. 5). In accordance with elast’ic theory of thin shells 2o) each portion of the pressure tube surrounding a fuel assembly may be classified as a long cylindrical shell. where the deformation in the central portion is independent of the end condition of each portioll. This assumption appears to be confirmed by

creep

the test data shown in figs. 3 and 7 of ref. 5).

rates 5, 6, were normalized in 6) for temperature and fast neutron flux using the empirical linear correction 5) based on biaxial data only.

The tube sections exposed to fast neutron flux 6, also qualify as long shells and consequently the only non-zero stresses existing in accordance with thin shell bheorem in the tube sections 5, 6, where creep rates were measured are

In addition, the tubes considered purposes by a adjustment for

hoop stress of all cold-worked 596) was reduced for correlation constant fact’or of 0.5. an biaxiality based on known

texture coefficients IT), and reasonable grouping of uniaxial and biaxial data along a convex curve was obtained by plotting normalized uniaxial and diametral creep rates versus uniaxial and adjusted stress, respectively, on a log-log graph. In this work, test data of cold-worked zirconium alloys shown 14, 5, 6)> are correlated through eqs. (1) and (4) by finding in accordance with continuum mechanics B and m from

oB= p,rlh = 2~~= constant,

(7)

assuming that the end conditions of the pressure tubes represent closed-ends. go and a, designate hoop and axial stress respectively and p is the internal pressure of the tube. The effects of elastic deformation of the tube due to temperature and internal pressure will be neglected for the creep analysis since they do not change the original shape of the tubes significantly. The stresses in the central portion

ON

of each tube recalculating

SECONDARY

section

CREEP

should

OF

ANISOTROPIC

be adjusted

by

NUCLEAR

55

MATERIALS

and

h and r for creep strains larger

than 0.05% which is the upper limit for the infinitesimal theory of small displacements. But this adjustment up to the maximum measured diametral strain of about 0.2% is insignificant in comparison assumptions limitations

with the errors due to both the made above and the experimental

indicated

(9b) Similarly,

eq. (4) may be expressed

in view

of eqs. (% (I), (7) and (5) for uniaxial biaxial cases as

and

51 6, 14). (lOa) and

4.

Correlation of in-reactor cold-worked Zircaloy-2

creep data of

As stated by Fidleris 14), uniaxial

LO= ((aA,,+A,,+A,,)t/(A,,+3A&)}~~= test data

are inconclusive with respect to the creep rate dependence on fast neutron flux +, because the range of flux tested was 0.36 to 1.16 x 1013 n/cm2esec only where 11 denotes numbers of neutrons. But the results 5) indicate an almost linear dependence of diametral strain on integrated fast neutron flux if stress and temperature are constant. Because of this linearity the variation of creep rates with fast neutron flux will be accounted for by normalizing as in ref. 5) both uniaxial and biaxial creep rates over periods of approximately constant flux, stress and temperature to a reference level of 1 x 1013 n/cm2. set and normalized creep rates will be denoted by .$i = iii/no, where n0=gj/(1013 n/cm2.sec). In view of the above simplifications,

(8) theorem

(1) or (4) may be used directly to calculate normalized creep rates corresponding to a state of stress indicated in eq. (7) using cylindrical r, 8, z coordinates, i.e., setting pairs of subscripts ll=z, 22=0, 33=r. Eq. (1) may be written for both &lo and 0~1 denoting normalized creep rate and stress respectively of uniaxial specimens with axes parallel to the longitudinal z-axis of the tubes and for 8,O the normalized hoop creep rate which is identical to the diametral creep rate of the tubes. Using eqs. (S), (7) and (5) in eq. (1) Pa)

=BO[(aAZz+Aee+Aez)tag]m=r60sm.

(lob)

Since ,80and m will be evaluated in accordance with eq. (9a) for cold-worked Zircaloy-2 using uniaxial experimental data 14)shown in table Al of the appendix, it is convenient to set AZz=l.

(11)

The data of table Al were selected as fairly reliable from a temperature range of 220 to 350 “C and a stress range of 11 000 to 40 000 psi since these are approximately the ranges at which the tubes 516) were tested. [In ref. 596~14) all measured values indicated are averages but in this approximate analysis instantaneous and average values will not be differentiated.] /?Oand m may be estimated from these limited uniaxial data and used for the pressure tubes by assuming that 1. The dependence

of normalized

creep rate

on temperature indicated by ,50 does not vary significantly with stress from 11 000 to 40 000 psi; 2. The stress dependence of creep rates measured by m is not significantly affected by temperature between 220 and 350 “C and does not vary with stress between 11 000 and 40 000 psi; 3. Neither /P nor m change significantly with fast neutron flux over the range of 0.36 to 3.1 x 1013n/cm2. sec. The last assumption is necessary, since the uniaxial tests were conducted at fluxes of 0.36

56

J. SCHROEDER

AND

M.

J.

HOLICKY

to 1.16 n/cm2. set and the tubes were tested

the estimates of ,P and m exhibited

at fluxes of 1.15 to 3.1 n/crnz. sec. The above

may be less useful outside these limits.

assumptions

had to be made in order to utilize

the limited

number

of results available.

even if the following to be qualitative suggestions

correlation

As remarked

But, and

the

concept

creep 14) is used,

of

activation

energy

/l’J= B exp ( -b/T)

B

300 “C and at fast neutron fluxes greater than 1.16 x 1013 n/cm2 -sec. Fig. 6 of ref. 5) indicates strain versus integrated neutron flux at various elevations of this tube from which average normalized creep rates d, were obtained as shown in table 2. In accordance with eqs. (9b), (12) and (ll), the coefficient

and b are constants and T is temperature (“K). Employing the method of least-squares for two independent variables 22) to the logarithmic form of eq. (9a) simplified by eq. (ll), the constants, B, b and m may be estimated using the selected data of table Al in connection with the assumptions (1) to (3) such that for cold-worked Zircaloy-2 ,P=B

A = (4 + A,, + A,JQm-1) (A,, + QA&) = 4/(B04?),

exp (-b/T)=

m=2.72;

(12)

Numerical values of /?O are tabulated in table 1 for various temperatures. /?Oand m shown in eq. (12) should be considered as estimates of cold-worked

Zircaloy-2

material

constants only, since they are based on limited experimental evidence and will change if more test results become available or if different data

A values in table 2 corrected

TABLE Values 250

i

260

270

for this tempera-

ture increase using eqs. (13) and (12) indicates that ,@J may be used over the full range of temperature and fast neutron flux shown in eq. (12). [For a consistent comparison with other tubes 270 “C will be used for U-2 Mk IV at elevation 261’-4” as the average of the inlet and outlet temperatures shown in table 1 of ref. 5) instead of the weighted average temperature of 264 “C.]

are selected from ref. 14). The numbers in eq. (12) are given to three significant figures for reasons of numerical computation only. The majority of the uniaxial data on which PO and m are based correspond to temperatures of 300 and 350 “C, stresses of 20 000 and 30 000 psi and fast neutron fluxes ranging from 0.36 to 1.16 n/ems. set as shown in table Al. Consequently

T W)

(13)

should be a constant for different combinations of pressure, temperature and fast neutron flux if the anisotropy coefficients are invariants. The above measure of anisotropy A was calculated at various elevations in accordance with eqs. (13) and (12) from the average creep rates shown in table 2, first using for each elevation 264 “C? which is the temperature at elevation 461’-4” weighted with respect to time in accordance with fig. 5 of ref. 5). Estimating an increase in temperature with elevation of 3.6 “C per foot from table 1 5) the consistency of the

= 26.0 x lo-13 exp [ - 8 390/T] ; T = 222 to 350 “C, stress= 11 000 to 40 000 psi ; pi=O.36 to 3.1 x 1013 n/cm2+sec.

the

tube in fig. 6 of ref. 5) may be used to test the validity of the estimate (12) for PO below

for

where

5)

of flux by the test data. However, detailed results available for the U-2 ;Llk IV pressure

for further work may be obtained

from the analysis. If

and Hunt

effects of temperature and stress on diametral creep rates are not as well defined as the effect

is considered

only, valid information

by Ross-Ross

in eq. (12)

1 of /?”

1

280

1 290

300

1

310

320

ON

SECONDARY

CREEP

OF

ANISOTROPIC TABLE

NUCLEAR

MATERIALS

57

2

Cold-worked Zircaloy-2 pressure tube U-2 Mk IV 5, A=(*+A,+LI,,)*(~-‘)

(ABe+*Ae,),

ue= 14 000 psi

/ Elevation

1(lO_S~i~in~h_l~)

464:--g”

11.4 9.5 8.5 7.0 6.5 6.1 6.1 6.5

463’-3” 462’-0” 461’-4” 460’~0” 458’-9” 458’-3” 456’-5”

adju:d for not adjtsted for ) temperature gradient temperature gradient 1.5 1.2 1.1 0.9 0.8 0.7 0.7 too close to end

of A calculated in accordance eq. (12) for m = 2.72 for various values

Magnitudes with

of A,,and A, are shown in table 3 and indicate changes of up to 100% for differences in A, and A, of about

20%

only.

TABLE 3 Change in A=(f+A,,+A,,)‘)“-l’(A,+Ae,) variation of ABo, A, for WL=2.72 A,, 1 1 1

/ A,, 1 1.2 0.8

1 A, -0.5 -0.3 -0.5

1 A 0.6 1.2 0.3

due to

1.0 1.0 1.0 0.9 0.9 1.0 1.0 too close to end

descriptions independent of temperature and for comparison of experiments and theorems, both sides of eqs. (9) and (10) have to be divided by 80. In order to facilitate this comparison numerical values of j30 are given in table 1. The uniaxial results shown in table Al of the appendix are plotted in fig. 1 as a,1 versus .$!J~Ousing a log-log scale and show reasonable scatter with respect to the straight line

Material type isotropic anisotropic anisotropic

(14)

Due to technical difficulties transverse normal

i.e., eq. (9a) combined with eq. (11) for values given in eq. (12). Two creep rates in fig. 1 are indicated at 45 000 psi 14) and appear to show a possible variation of m with stress. But since

strain rates of uniaxial creep specimens have not been measured and results of creep tests of uniaxial specimens with longitudinal

45 000 psi corresponds almost to first yield of the material 14) no conclusions can be drawn. Data of cold-worked Zircaloy-2 pressure

axes tangential to the circumferential direction of the tube are not available. Consequently the coefficients of anisotropy A,, and A, are unknown and direct correlation of uniaxial and biaxial test results through the power law (10) in terms of the effective strain rate and effective stress is not possible. However, in view of eq. (9) uniaxial and biaxial results from tubes should be grouped with respect to two parallel lines on a log-log graph of creep rate versus stress if the coefficients of anisotropy do not vary with tubes and m is invariant. For the following correlation, data are plotted in terms of &/~O, $//P, S/PO in order to use

tubes 5) are shown in table A2 of the appendix: the U-3 Mk V tube was not considered in this work because of low stress (about 5 000 psi) and an organic coolant. Results 6) of small cold-worked Zircaloy-2 tubes are tabulated in table A3 of the appendix and are represented with the data of table A2 in fig. 1 as $‘//?O versus c0 for comparison with uniaxial test results 14). The points associated with ref. 5) appear to be grouped with respect to a line where m is approximately one, whereas in view of eq. (9) data belonging to ref. 6) appear to be in agreement with uniaxial data if cB is greater than 20 000 psi. The value m= 1 was

58

J.

SCHROEDER

50

AND

M.

J.

HOLICKY

I I I ZIRCALOY-2

COLD-WORKED

0 *

40 rT 0 x

30

E 20b

UNIAXIAL

SPECIMEN

UNIAXIAL

SPECIMEN

A

PRESSURE

.

SMALL

.

REFERENCE

b b”

0 0

I

I

‘5.

/

Al

RX-21,

14 TUBES

TUBES

,L?= 260

TABLE

TABLE

TABLE

A2 A3

x10-‘3exp[-6390/TIoK;1

-,

/ (k

if,/ Fig.

1.

Cold-worked

represented

as

Zircaloy-2

uniaxial

or

in-reactor

hoop

stress

p”

creep versus

or data

;i

/ p”

from

uniaxial

normalized

by B”, eq.

hre’psi.e272

(psim

set cd

x IO”)

specimen

uniaxial

/

O+ 10’3.)

or

14) and tube

diametral

creep

specimens rate

5, 6)

divided

(12).

used by Ross-Ross and Hunt 5) in order to interrelate results of pressure tubes using a linear temperature correction. The disagreement with eq. (9) indicated in fig. 1 has two possible explanations : 4. The power index m varies with stress for cold-worked Zircaloy-2, or 5. The cold-worked Zircaloy-2 coefficients of anisotropy are stress dependent. Explanation (4) is feasible since only one

ends assumed in eq. (7) may not have been perfect. But since in 6) an open loop was used, i.e., perfect closed end condition existed, and the experimental data 6) show the same trend as those 5, at stresses below 20 000 psi, this explanation does not appear to be warranted. The validity of explanations (4) and (5) can be proven only by conducting more uniaxial creep tests in the react’or. The possibility that (5) is valid may be

uniaxial

investigated by considering whether or not a possible relationship between the measure of anisotropy A defined in eq. (13) is indicated

experimental result 14) exists below 18 000 psi as shown in fig. 1 and for example Kachanov 8) has indicated, that creep data are sometimes grouped along a line consisting of two segments of different slopes or along convex curves on a log-log graph of creep rate versus stress. In-reactor creep may depend on slip and/or diffusion 21) and rotation of basal planes of the hexagonal crystals under stress suggested by Hill 23) and/or dependence of magnitude and direction of diffusion on stress indicated by Garofalo 24), may account for a variation of anisotropy with stress, i.e., explanation (5). The discrepancy shown in fig. 1 could be due to a variation of stress ratio ae/az for the pressure tubes 5) with load since the condition of closed

by data 53 6) for a constant

m=2.72.

Values

for A calculated for all the tubes in accordance with eq. (12) are shown in tables A2 and A3 and are plotted versus stress in fig. 2a. There appears to be a significant dependence of A on stress and fig. 2b indicates that the overall change in A can be explained by variations in A,, or A, of reasonable magnitude which for tubes tested in the NDP reactor may be classified in the following way: A,,>l,

A,< -0.5; = -0.5;

0.8
A,,
A,>

A,= -0.5.

(15)

In view of eqs. (9b), (11) and (13) the effect

ON

SECONDARY

CREEP

OF

ANISOTROPIC

59

MATERIALS

on A of a possible variation

1 COLD-WORKED ZIRCALOY -2 . NDP PRESSURE TUBES TABLE A2 . NRU PRESSURE TUBES TABLE A2 . NRU SMALL TUBES TABLE A5 COLD-WORKED ZIRCONIW-25% ANRU PRESSURE TUBES TABLE

NUCLEAR

in Au is reduced

by a factor

of $ and may be neglected. Since (5) appears to be a possible explanation

NS A5

of the discrepancy constant

indicated

for cold-worked

in fig. 1 if m is a

Zircaloy-2

over the

stress range shown in eq. (12), it may be worthwhile to select comparison

2 011

I

IO

I

,

15

20

(0)

25

I 35

30

40

2

Fig.

2a.

eq.

(13)

Dependence on hoop

IO

12

of measure

stress for cold-worked

for anisotropy

cold-worked

14

coefficients

zirconium alloy tubes

anisotropic

strain

view of eq. (10) with respect to one straight line on a log-log graph of 9//30 versus s, if both 80 and s are adjusted for variations in A,, and A, with stress. Selecting from the curve shown in fig. 2a a value for A corresponding to the hoop stress of the tubes considered 59s), a comparison value for A,, was found from

16

of anisotropy

alloy tubes 5, 6) if m=2.72=constant. values

and (11) effective

rates do and stress s for the measured diametral creep rates and hoop stresses of the tubes 516). Because A was calculated from eq. (13) for a constant m=2.72, data from uniaxial specimens 14) and tubes 596) should be grouped in

+

08

and with

eqs. (lob)

U# (psi x IO31

06

values of A,

A, from fig. 2 and calculate in accordance

A

zirconium

b. Comparison A,,

and

5,

if m=2.72=

6)

A,

for

fig. 2b using a constant A,= - 0.5 as shown in tables A2 and A3. The effective strain rates and stresses calculated for the tubes in terms of these comparison values correlate as predicted with the uniaxial data shown in fig. 3. If different comparison

constant.

values for A,, and A,

are selected

40 30

;“I Fig. 3.

Cold-worked

represented

Zircaloy-2

as effective

in-reactor

anisotropic

p*

(psimx

IO”)

creep data from uniaxial

stress versus effective sssuming

specimen 14) and tube

anisotropic

m=2.72=constant.

strain rate divided

specimens

9

)

by PO, eq. ;12;,

60

J.

in accordance are obtained.

agreement

was

10, ref. 6) by grouping

indicated

both

a convex

curve

uniaxial

and

diametral

uniaxial

and

one

respectively.

creep

half

The hoop

Zircaloy-2

on a log-log of

graph

rates

the

hoop

of

versus stress

stresses were reduced

by one half in order to account for biaxiality in accordance with known texture coefficients assuming that creep is mainly caused by slip between grains. This reduction of the hoop stress is analogous to a correlation using an anisotropic creep theorem and the method of correlation indicated in fig. 10, ref. 6) may be simulated by setting in eq. (lob) the coefficient of $ to one and the coefficient of (T@to one half. It follows that independent of m and /30 the value of A,,= 1 and of A,,= - 1. The corresponding magnitude of A calculated for m= 2.72 from eq. (13) is 0.15 and appears to be too low even for the small tubes as shown in fig. 2b mainly because in this work a temperature dependence based on uniaxial results was used in accordance with continuum mechanics, whereas 5) a linear temperature correction based on biaxial results 5) only was employed for both uniaxial and biaxial data. COLD-WORKED ITi1

- _ _

v-l

0

ZIRCONIUM-2.5%

M.

J.

5.

in

uniaxial 14)

and biaxial 516) data of cold-worked along

AND

with fig. 2, similar correlations

Reasonable fig.

SCHROEDER

HOLICKY

Correlation of in-reactor creep data of cold-worked zirconium-2.5 wt o/o niobium Only three cold-worked

niobium

zirconium-2.5

wt, y.

uniaxial specimens were tested la), all

three at an average temperature

of 300 “C as

shown in table A4 of the appendix.

In fig. 4, the

power index m = 2.72 estimated for cold-worked Zircaloy-2

in eq. (12) appears to be in agreement

with

three

the

uniaxial

data

obtained

for

zirconium-2.5 wt o/o niobium and due to lack of experimental evidence this value of m and the variation of normalized creep rate with temperature shown in eq. (12) will be used such that as an est’imate for cold-worked zirconium-2.5 wt o/o niobium PO= 8.5 x IO-13 exp [- 8 390/T (OK)], m= 2.72,

(16)

where the coefficient 8.5 x lo-13 was estimated from eq. (14) in connection with the average of the three creep rates shown in table A4. Data of cold-worked zirconium-2.5 wt y. niobium pressure tubes obtained from 59 6) are shown in table A5 of the appendix and are plotted as &O//30versus u0 on a log-log scale in fig. 4a where a similar discrepancy with uniaxial data exists as indicated in fig. 1 for cold-worked Zircaloy-2. The measures of anisotropy A in table A5 are calculated in accordance

NB 25

25

x ‘g

t?

20

15

ii

A PRESS. TUBES TABLE A5

tf

IO

I

2 ;z,

/fl”

3

45

IO

or 2; /p” (psimx IO”)

s? z x ‘3 e YI

15

IO

I

2 8”/

3 6”

45

IO

(psimx IO”)

(b)

(a) Fig. 4.

Cold-worked zirconium-2.5 wt y. niobium in-reactor creep data from uniaxial specimen 14) and tube specimens 5. 6) represented in (a) as uniaxial and hoop stress versus normalized uniaxial and diametral creep rate divided by PO, eq. (16) and in (b) as effective anisotropic stress versus effective anisotropic

strain rate divided

by ,P, eq. (16) assuming

WL= 2.72=constant.

ON

SECONDARY

with eqs. (13) and (16) similar dependence

CREEP

OF

ANISOTROPIC

and appear to have a

effective

anisotropic

Ross-Ross, Atomic Energy of Canada

1) P. A. Limited

Report,

AECL-3126

(1968)

2, W. Evans, J. E. Le Surf and W. R. Thomas, At,omic

stress

Energy

AECL-2890

versus effective anisotropic strain rate divided by ,80 was established in a similar manner as

B. A.

of

Canada

Cheadle and W.

of Canada

Limited

Evans,

Report,

fig. 3 and indicates the same type of correlation.

Report,

AECL-1048

6.

P. A. Ross-Ross

and C. E.

Limited

Conclusions

of the material. If the power index for coldworked zirconium alloys is constant the variation of the coefhcients of anisotropy indicated by the variation of the so-called measure of anisotropy A shown in fig. 2a, may be used to draw conclusions regarding the basic inreactor creep mechanism 21) since A appears to approach a constant or at least a minimum value.

Mat. E.

The author wishes to acknowledge the valuable assistance and advice given by Mr. Ross-Ross, Fuels and Materials Division, Atomic Energy of Canada Limited, for both the formulation of the problem and the evaluation of the experimental data.

Report,

Atomic

Energy

AECL-2652

(1966)

Atomic

Energy

of Canada

(1960) L.

Hunt,

J. Nucl.

2

Ibrahim,

In-reactor

creep of zirconium

alloy tubes and its correlation with uniaxial data, ASTM

Symp.

delphia;

7)

I.

Berman

8 (1966)

8)

zirconium

Nov.,

and D.

L. M. Kachanov, tizdat,

and

hafnium

(Phila-

1968) H.

Pai,

Int.

J. Mech.

Sci.

341 The theory

Moscow,

Lending

1960)

Library

Boston

of creep (Fizma-

; Engl. transl., National

for Science

Spa, Yorkshire,

and

England

Technology,

(1967)

Q) K. G. Odqvist, Mathematical theory of creep and creep rupture

loI

(Oxford,

Y. N. Rabotnov, Progress

in

New York,

1966)

The Prager anniversary volume

Applied

Mechanics

(MacMillan,

1963) p. 307

11) N. I. Malinin, PMFT Zh. Prikl. Mat. Tekh. Fiz. 3 (1964)

16

12) W. Olszak and A. Sawczuk, Inelastic behaviour in shells (Noordhoff, Groningen, The Netherlands, 1967) 13

)

A. M. Freudenthal pedia of Physics

and H.

Geiringer,

(Springer,

V. Fidleris, J. Nucl.

15)

W. Prager, Introduction to mechanics of Continua

9

R. Hill,

(Ginn, New York, 17

)

B.

19

Mech.

A.

Cheadle

and

Techn.

4 (1966)

329

Groningen, )

D. H. S.

29

26 (1968)

51

1961)

J. Appl.

V. V. Novozhilov,

13)

Mat.

Encyclo-

1958) p. 229

14)

17 (1950)

C.

E.

64

Ells,

Electrochem.

Thin shell theory (Noordhoff,

The Netherlands,

1964)

Pai, Int. J. Mech. Sci. 9 (1967)

Timoshenko

Theory

Acknowledgment

Perryman,

26 (1968)

F.

Limited

(1967)

E. C. W.

The anisotropic creep theorem proposed in this work appears to be useful to analyse uniaxial and biaxial in-reactor creep behaviour of coldworked zirconium alloys. Using the limited results from uniaxial tests 14) it was shown how both uniaxial data and results of biaxially stressed tubes596) may be correlat,ed. In order to analyse the behaviour of zirconium alloys in the reactor throughly it is necessary to test assumptions (1) to (3) by conducting more uniaxial creep experiments which will automatically lead to conclusions regarding explanations (4) or (5) given in section 4. Furthermore, it appears to be essential to determine the coefficients of anisotropy Ace, A, and their variation with stress if the power index m is found to be invariant, in order to facilitate the correlation of these coefficients with texture

61

MATERIALS

References

on stress as the A values

of the cold-worked 7ircaloy-2 tubes for an 17~ of 2.72 as shown in fig. 2a. The log-log graph in fig. 4b indicating

NUCLEAR

Toronto,

of

and

plates

S.

and

335

Woinowsky-Krieger, shells

(McGraw-Hill,

1959)

21) G. R. Piercy, J. Nucl. Mat. 26 (1968) 18 22

)

E. F. Croxton and D. J. Cowden, Applied general statistics

23

24

)

1

(Prentice-Hall,

1955)

R. Hill,

The mathematical

(Oxford

Clarendon Press, 1964) p. 317

F.

Garofalo,

Fundamentals

rupture in metals p. 178

theory

of plasticity

of creep and creep

(MacMillan,

New York,

1965)

Appendix Al

TABLE

Data

based

on uniaxial

test results

Test

Temp.

no.

(“C)

R-6

of cold-worked

Stress (lO:‘,si)

Fast

flux

rate

table

/

0.54

3.5 I_t 107;

5.7

0.60

13.0 & loo/:,

5.8

300

18

0.58

3.0 *

20%

5.1

300

11

0.61

0.5 i

600,b

0.7

300

20

0.64

3.0 & 3076

4.1

320

20

0.68

6.0 & 20%

5.1

350

20

0.68

14.0 5 200,;

5.6

220

30

0.96

1.5 *

3076

11.5

260

30

0.93

5.0 & 4o”/b

14.2

350

30

0.94

32.0 + 10%

9.2

R-2

300

30

0.66

15.0 & 10%

19.9

R-4

300

20

0.58

3.5 & 30%

5.3

300

30

0.36

8.0 + lo”:,

20.0

Rx-2(

350

30

1.16

60.0 f

lOo,d

14.0

Rx-21

300

40

0.84

35.0 *

5076

36.5

1 and fig.

8 of ref. 5,

Zircaloy-2

TABLE Data

based Data

on results

based

Av. of inlet

Hoop

and outlet

stress

I temperature (“C)

selected

on results

I

tubes

A2

from

table

TABLE

A3

shown

in figs.

3 and

7 of ref. 6, ~

~ Diametral

A,,

comparison

ng

I(lo3psi) Table

264

IJ)

( x 101’)

20

Cold-worked

4 of ref.

i:1 IB” eqs. (8, 12)

‘z (1O-7 (:n;in)h-I)

20

Rx-1I

NDP

from

350

Rx-14

designation

Creep

~

1013 n/zmj.sec)/

selected

300

R-9

Tube

Zircaloy-2

A2

U-2 Mk III

281

~ 10.5 i 11.5

U-2 Mk IX U-2 Mk IV

281

13.4

3.1

270

~ 14.0

2.7

1.9

1.4

0.73

~

1.06

U-2 Mk IX

281

17.3

3.1

2.9

I

1.4

0.40

~

0.92

258

15.7

2.56

1.3

~

0.98

19.4

2.56

2.4

I I

0.55

258

2.6

258 258

23.9 29.2

i

2.56 2.56

2.7 3.6



2.56 2.56

1

~

~

1.15

0.54

1.1

2.7

1.8

0.9

~ 1.27 i 0.86

2.3

1.1

0.65

~

Table

258 258

I

~ 34.0

37.6

I

1

I

1.28 1.20 1.09

A3

6.1 8.8

1.4 2.9 3.9

;

6.6 9.6

I ~

0.57

0.87

0.36

0.80

0.28

0.80

0.31 0.34

0.80 0.80

ON

SECONDARY

CREEP

OF

ANISOTROPIC

NUCLEAR

63

MATERIALS

TABLE A4 Data based on uniaxial test results of cold-worked zirconium-2.5 wt ye niobium selected from table 4 of ref. 14)

R-11 R-12 R-17

~- Creep rate

Fast flux

I_

Test no.

“” (10-e (in~m)h-1)

/ / )

300 300 300

i j /

23.0 16.5 20.0

20.0 & 15% 7.0 & 10% 10.5 -& 40%

o.7

\

0.56 0.68

1 ; ;

7.7 3.4 4.2

TABLE A5 Data based on results from cold-worked zirconium-2.5 wt ye niobium pressure tubes selected from ref. 5, 6)

Tube designation

Av. of inlet and outlet temperature

EIoop stress

Diametral creep rates

Fast flux

9, (1013n/om2*sec) 1(10-7 (~~]in)h-l)l

(“C)

/ (lG:sii)

U-l Mk VIII

270 270

16.5 22.0

2.9 2.9

U-1MkV

285 285

15.5 21.0

3.1 2.0

/ /

W” eqs. (8, 12) ( X 1Olif

A eq. (13)

6.5 10.4

1.35 2.16

0.46 0.33

6.7 7.8

1.07 1.54

0.43 0.27

A,, comparison value selected from fig. 2b for A,= -0.5

1

0.96 0.80 1.01 0.80