HSNV123 - Thermo metal-worker-mechanics EDGAR

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

1/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

Organization (S):

EDF-R & D/AMA, CS IF

Manual of Validation
V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures
V7.22.123 document

HSNV123 - Thermo metal-worker-mechanics EDGAR

Summary:

This test illustrates a mechanical calculation on a material (Zircaloy) undergoing metallurgical transformations.
Concretely, initially, operator CALC_META calculates the metallurgical evolution associated with

a given thermal history. This metallurgical evolution is then provided to STAT_NON_LINE which goes

to carry out a mechanical calculation by taking of account metallurgical phases (in addition to loadings
mechanics). The material of mechanical calculation is defined with ELAS_META_FO, META_ECRO_LINE and
META_VISC_FO.

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

2/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

Problem of reference

It is about a cylindrical bar in creep.

1.1 Geometry

Appear 1.1-a: Geometry and loading of the problem of reference

It is about a cylinder H= height 1.0 m, and radius R=1.0 Mr.
The square in fat corresponds to axisymmetric modeling used to [§3].

1.2

Material properties

The properties materials are described by the following parameters:

For metal thermo calculation

(Zircaloy)
C

= 2000000 J.m

- 3

.°C

- 1

= 9999.9 W.m

- 1

.°C

- 1

Coefficients for the metallurgy:
teqd = 809 °C, K=1.135

- 2, n=2.187

CCT = 831 °C, qsr = 14614, ac = 1.58E-4
m = 4.7, tdr = 949,1°C, Ar = - 5.725, Br = 0.05

For calculation thermo metal-worker-mechanics

· Young modulus: E= 200000 AP
· Poisson's ratio: = 0.3

Definition of the elastic characteristics, dilation and elastic limits for
modeling of an undergoing material of the metallurgical transformations:

· T

ref.

= 800°C

· Thermal expansion factor average of the cold phases:

(T) = 0

· Thermal expansion factor average of the hot phase:

(T) = 0

· Temperature of definition of the expansion factor: T

= 800°C

· Choice of the metallurgical phase of reference: heat
· Deformation of the phase not of reference compared to the phase of reference to

temperature T

ref.

= 0

· Elastic limit of the cold phase 1 for a viscous behavior: F_sigm_f (T) = 0
· Elastic limit of the cold phase 2 for a viscous behavior: F_sigm_f (T) = 0
· Elastic limit of the hot phase for a viscous behavior: to see [Figure 1.2-a]

F=25 NR

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

3/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

100

150

200

250

200

400

600

800

1000

Temperature (°C)

C
_
S_

Figure 1.2.- a: Limit elastic of the hot phase for a viscous behavior

· Function used for the law of mixture on the elastic limit of multiphase material for

a viscous behavior: F

0,2

0,4

0,6

0,8

1,2

0,2

0,4

0,6

0,8

1,2

META

_
M
EL

With
NR
G
E

Appear 1.2-b: Law of mixture

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

4/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

Definition of the modules of work hardenings used in the modeling of the phenomenon
of linear isotropic work hardening of an undergoing material of the phase shifts
metallurgical:

· Slope of the traction diagram for the cold phase 1

5000

10000

15000

20000

25000

30000

200

400

600

800

1000

1200

Tem pérature (°C)

F
1_D

I
G
M_E

P
S
I

Appear 1.2-c: Traction diagram for the cold phase 1

· Slope of the traction diagram for the cold phase 2:

F (T) = 0

· Slope of the traction diagram for the hot phase:

F (T) = 0

Definition of the viscous parameters of the viscoplastic law of behavior with catch in
count metallurgy:

· Parameter of the viscoplastic law of flow, for the cold phase 1

500

1000

1500

2000

2500

200

400

600

800

1000

1200

Temperature (°C)

_
E
MT

Appear 1.2-d: Parameter

viscoplastic law of flow, for the cold phase 1

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

5/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

· Parameter of the viscoplastic law of flow, for the cold phase 2

50000

100000

150000

200000

250000

300000

500

1000

1500

Temperature (°C)

_
E
T
With

Appear 1.2-e: Parameter

viscoplastic law of flow, for the cold phase 2

· Parameter of the viscoplastic law of flow, for the hot phase:

F (T) = 0

· Parameter N of the viscoplastic law of flow, for the cold phase 1:

F (T) = 5.76

· Parameter N of the viscoplastic law of flow, for the cold phase 2:

F (T) = 2.94

· Parameter N of the viscoplastic law of flow, for the hot phase:

F (T) = 1.0

· Parameter C relating to the restoration of work hardening of viscous origin, for the cold phase

F (T) = 13.70539827

· Parameter C relating to the restoration of work hardening of viscous origin, for the cold phase

F (T) = 0

· Parameter C relating to the restoration of work hardening of viscous origin, for the hot phase

F (T) = 0

· Parameter m relating to the restoration of work hardening of viscous origin, for the phase

cold 1:

F (T) = 5.76

· Parameter m relating to the restoration of work hardening of viscous origin, for the phase

cold 2:

F (T) = 1.0

· Parameter m relating to the restoration of work hardening of viscous origin, for the phase

heat:

F (T) = 1.0

1.3

Boundary conditions and loadings

The base of the cylinder is locked according to y:
Uy = 0 on the basis of cylinder
A force of traction F=25 NR is imposed on the top of the cylinder

The temperature is imposed on all the cylinder for t=120s.
T (X, y, 120) = 800°C

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

6/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

1.4 Conditions

initial

The following variables are initialized:
T (X, y, 0) = 800°C
V1 (X, y, 0) = 1.0
V2 (X, y, 0) = 0.0
V3 (X, y, 0) = 20.

V1: proportion of the cold phase

V2: proportion of the cold phase

, mixed with the phase

V3: temperatures with the nodes

Reference solution

2.1

Results of reference

The results of reference were obtained with a former version of aster. It is about a test of
not-regression.

2.2

Uncertainty on the solution compared to the result of not-regression

Uncertainty is 10%.

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

7/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

3 Modeling

With

3.1

Characteristics of modeling

The modeling used in the case test is as follows:

Elements 2D “AXIS” (QUA8)

Appear 3.1-a: Geometry and mesh of modeling

Cutting:

2 meshs QUAD8 according to the x axis

2 meshs QUAD8 according to the y axis

Boundary conditions:

Uy=0 on D1

F=25N on D2

3.2

Characteristics of the mesh

A number of nodes: 21
A number of meshs and types: 4 QUAD8, 8 SEG3.

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

8/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

3.3 Functionalities

tested

Controls

DEFI_MATERIAU THER

RHO_CP

LAMBDA

META_ZIRC

TDEQ
K
NR
CCT
QSR_K
AC
M
TDR
AR
Br

CALC_META ETAT_INIT

COMP_INCR

META_INIT_ELNO
RELATION

“ZIRC”

DEFI_MATERIAU ELAS_META

NAKED

F_ALPHA

C_ALPHA

TEMP_DEF_ALPHA

PHAS_REFE

EPSF_EPSC_TREF

F1_S_VP
F2_S_VP
C_S_VP
S_VP_MELANGE

META_ECRO_LINE

F1_D_SIGM_EPSI
F2_D_SIGM_EPSI
C_D_SIGM_EPSI

META_VISC_FO

F1_ETA
F2_ETA
C_ETA
F1_N
F2_N
C_N
F1_C
F2_C
C_C
F1_M
F2_M
C_M

MODEL AFFE_CHAR_MECA

TEMP_CALCULEE

MODEL STAT_NON_LINE

CHAM_MATER

EXCIT

CHARGE

FONC_MULT

COMP_INCR

RELATION

“META_V_IL”

RELATION_KIT

“ZIRC”

INCREMENT

LIST_INST

NEWTON

STAMP

“TANGENT”

REAC_ITER

CONVERGENCE

RESI_GLOB_RELA

ITER_GLOB_MAXI

Code_Aster

Version

8.1

Titrate:

HSNV123 - Thermo metal-worker-mechanics EDGAR

Date:

05/09/05

Author (S):

V. CANO, S. VANDENBERGHE

Key

V7.22.123-A

Page:

9/10

Manual of Validation

V7.22 booklet: Thermomechanical nonlinear statics of the voluminal structures

HT-66/05/005/A

Results of modeling A

4.1 Values

tested

Identification Size Reference

Aster %

difference

t=120s m3 N5

EPYY

- 3.1E-2

- 2.888E-2

- 6.8%

t=120s m3 N5

SIYY

- 25.0

- 24.99

- 8.90E-5%

5 Comments

This case test of not-regression makes it possible to check the coherence of Code_Aster of a version on the other
with regard to the metallurgy.