Microsoft Word - V604301a. Doc.

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

1/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Organization (S)

: EDF/IMA/MN, SAMTECH SA

Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
Document: V6.04.301

SSNV301 - Cylindrical ring in subjected rotation
with a thermal shock and an internal pressure -
Von Mises (isotropic Work hardening)

Summary:

This test of quasi-static mechanics non-linear 3D relates to the modeling of a piece of ring
cylindrical and the resumption of a field of deformation checks after a transitory linear thermal calculation (shock
thermics: application of a heat flux inside the cylindrical ring). Then, one takes calculation and one again
check the transfer of the field of thermal deformation of origin and the use of a following loading (force
centrifugal). Then, one takes again calculation once again by applying a pressure interns (change of
boundary conditions). One uses the criterion of Von Mises with isotropic work hardening defined by a function which
depends on the temperature.

Comparison with the numerical results of a modeling with the SAMCEF software V7.0 (module THERNL for
calculation of the field of temperature transitory and taken again with module MECANL for mechanical calculation
non-linear geometrical (“Green”) and material).

Modeling with elements 3D of the meca_hexa20 type and axisymmetric elements of meaxtr6 type. One
use a rather fine mesh in order to check the robustness and the performances of

Code_Aster.

In two modelings A and B, the pressure applied is not following. Results obtained by
The SAMCEF software are provided with and without following pressure.

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

2/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Problem of reference

1.1 Geometry

Roll with:

Interior radius

= 0.04 m

External radius

= 0.08 m

Height

L = 0.04 m

1.2

Material properties

The material is homogeneous isotropic, thermoelastic linear. Material mechanical properties
depend on the temperature. The properties vary linearly between the 2 extreme temperatures:
0°C and 100°C.

Law of behavior to linear isotropic work hardening:

E = 195000 MPa

= 0.3

= 181 MPa

= 2857.504 MPa

Temperature = 0

E = 175500 MPa

= 0.3

= 145 MPa

= 1979.695 MPa

Temperature = 100

Thermal properties and density:

7860

452

6 1

msK

kgK

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

3/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

1.3

Boundary conditions and loadings

The thermal shock is modelized by a heat flux of 300000 W/m

S inside the cylinder.

structure has an initial temperature of 0°C.

At the same time a centrifugal force is applied: on an interval of 15 S, the rotational speed is
increased up to 2400 rad/S, then kept constant.

(rad/S)

2400

T (S)

After 15 S, an internal pressure is applied: the pressure increases linearly up to a value
from 5.5 MPa which is reached to 20 S, then the pressure goes down to zero.

5.5

T (S)

p (Mpa)

Thermal calculation and mechanical calculation are carried out with 50 pitches of identical times.

A quarter of the cylinder is modelized by fixing the transverse components of displacements at both
cross sections of the cylinder. Moreover, the degree of translation according to the length of the cylinder is
fixed on the interior radius of the with dimensions inferior of the cylinder.

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

4/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Reference solution

2.1

Method of calculation used for the reference solution

Comparison with the numerical results of a modeling voluminal with the SAMCEF software V7.0 with 50
no time (module THERNL for the calculation of the field of temperature transitory and taken again with
modulate MECANL for non-linear mechanical calculation).
An implicit scheme is used for the integration of the equation of heat.

Non-linear mechanical calculation was carried out with the SAMCEF software by using a modeling
axisymmetric with 5 pitches of time (geometrical nonlinear modeling (“Green”)) and one
voluminal modeling (a quarter of cylinder) with 5 pitches of time (modeling nonlinear
geometrical (“Green”)).

2.2

Results of reference

Evolution of the field of temperature according to the thickness on with dimensions lower of the cylinder

Axisymmetric calculation with 50 pitches of time

Temp = 0.5 dryness

Temp = 25 dryness

1.25355E+01

8.64267E+01

5.71233E+00

7.65695E+01

2.44526E+00

6.77355E+01

1.12189E+00

5.98610E+01

4.83644E-01

5.28476E+01

2.22443E-01

4.66462E+01

9.62036E-02

4.11801E+01

4.43444E-02

3.64106E+01

1.92310E-02

3.22765E+01

8.87973E-03

2.87468E+01

3.85960E-03

2.57723E+01

1.78524E-03

2.33283E+01

7.78650E-04

2.13761E+01

3.63114E-04

1.98963E+01

1.64448E-04

1.88596E+01

8.80577E-05

1.82514E+01

6.54904E-05

1.80507E+01

Voluminal calculation with 50 pitches of time

Temp = 0.5 dryness

Temp = 25 dryness

1.26658E+01

8.64175E+01

5.64556E+00

7.65753E+01

2.50690E+00

6.77332E+01

1.10448E+00

5.98571E+01

5.06396E-01

5.28510E+01

2.23330E-01

4.66453E+01

9.98984E-02

4.11819E+01

4.41255E-02

3.64095E+01

2.02976E-02

3.22781E+01

8.97177E-03

2.87460E+01

4.05352E-03

2.57736E+01

1.79338E-03

2.33275E+01

8.25938E-04

2.13766E+01

3.68531E-04

1.98979E+01

1.74367E-04

1.88647E+01

8.91805E-05

1.82450E+01

6.98548E-05

1.80607E+01

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

5/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Evolution of the longitudinal stress

With

0.5 L

0.25 (R

- R

)

Case where the pressure is nonfollowing

Nonlinear voluminal calculation

geometrical with 5 pitches of time

T (S)

Not A

Not B

15.0

1.21692E8

1.67285E8

20.0

1.39804E8

1.65266E8

25.0

1.28967E8

1.57980E8

Nonlinear axisymmetric calculation

geometrical with 5 pitches of time

T (S)

Not A

Not B

15.0

1.20472E8

1.66917E8

20.0

1.39164E8

1.70412E8

25.0

1.28421E8

1.63168E8

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

6/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Case where the pressure is following

Nonlinear voluminal calculation

geometrical with 5 pitches of time

T (S)

Not A

Not B

15.0

1.21692E8

1.67285E8

20.0

1.39894E8

1.65313E8

25.0

1.29032E8

1.58007E8

Nonlinear axisymmetric calculation

geometrical with 5 pitches of time

T (S)

Not A

Not B

15.0

1.20472E8

1.65934E8

20.0

1.39252E8

1.64075E8

25.0

1.28485E8

1.56886E8

2.3

Uncertainty on the solution

Uncertainty lower than 1% for thermal calculation, uncertainty lower than 0.5% for calculation
mechanics.

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

7/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

3 Modeling

With

3.1

Characteristics of modeling A

A quarter of cylinder with a regular mesh.

3.2

Characteristics of the mesh

A number of nodes:

4037

A number of meshs and types:

800 MECA_HEXA20 (8 meshs on
the thickness, 10 meshs on
length and 10 meshs on the quarter
of circle)

3.3 Functionalities

tested

Controls

Keys

THER_LINE

PARM_THETA

[U4.33.01]

DEF_MATERIAU

ECRO_LINE_FO

D_SIGM_EPSI

[U4.23.01]

STAT_NON_LINE

COMP_INCR

DEFORMATION

PETIT_REAC

[U4.32.01]

RELATION

VMIS_ISOT_LINE

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

8/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Results of modeling A

4.1 Values

tested

Identification

Moment

Reference

Aster

% difference

Temperature (N21)

0.5s

1.266580E1

1.266584E1

Temperature (N201)

0.5s

5.645560E0

5.645554E0

Temperature (N202)

0.5s

2.506900E0

2.506894E0

Temperature (N203)

0.5s

1.104480E0

1.104478E0

Temperature (N204)

0.5s

5.063960E1

5.063949E1

Temperature (N205)

0.5s

2.233300E1

2.233298E1

Temperature (N206)

0.5s

9.989840E2

9.989826E2

Temperature (N207)

0.5s

4.412550E2

4.412536E2

Temperature (N208)

0.5s

2.029760E2

2.029754E2

Temperature (N209)

0.5s

8.971770E3

8.971745E3

Temperature (N210)

0.5s

4.053520E3

4.053507E3

Temperature (N211)

0.5s

1.793380E3

1.793372E3

Temperature (N212)

0.5s

8.259380E4

8.259340E3

Temperature (N213)

0.5s

3.685310E4

3.685292E4

Temperature (N214)

0.5s

1.743670E4

1.743665E4

Temperature (N215)

0.5s

8.918050E5

8.918007E5

Temperature (N28)

0.5s

6.985480E5

6.985445E5

Temperature (N21)

25s

8.641750E1

8.641753E1

Temperature (N201)

25s

7.657530E1

7.657529E1

Temperature (N202)

25s

6.773320E1

6.773321E1

Temperature (N203)

25s

5.985710E1

5.985709E1

Temperature (N204)

25s

5.285100E1

5.285104E1

Temperature (N205)

25s

4.664530E1

4.664533E1

Temperature (N206)

25s

4.118190E1

4.118188E1

Temperature (N207)

25s

3.640950E1

Temperature (N208)

25s

3.227810E1

3.227812E1

Temperature (N209)

25s

2.874600E1

2.874597E1

Temperature (N210)

25s

2.577360E1

2.577356E1

Temperature (N211)

25s

2.332750E1

2.332747E1

Temperature (N212)

25s

2.137660E1

2.137655E1

Temperature (N213)

25s

1.989790E1

1.989793E1

Temperature (N214)

25s

1.886470E1

1.886472E1

Temperature (N215)

25s

1.824500E1

Temperature (N28)

25s

1.806070E1

1.806074E1

SIYY (M473, N2931)

15s

1.216920E8

1.212172E8

0.390

SIYY (M74, N204)

15s

1.672850E8

1.672363E8

0.029

SIYY (M473, N2931)

20s

1.398040E8

1.396296E8

0.390

SIYY (M74, N204)

20s

1.652660E8

1.654056E8

0.029

SIYY (M473, N2931)

25s

1.289670E8

1.287516E8

0.390

SIYY (M74, N204)

25s

1.579800E8

1.580668E8

0.029

4.2 Parameters

of execution

Version: 3.09
Machine:

System:

Overall dimension memory:

16 megawords

Time CPU To use: 4971.07 seconds

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

9/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

5 Modeling

5.1

Characteristics of modeling B

5.2

Characteristics of the mesh

A number of nodes:

357

A number of meshs and types: 160 MEAXTR6

(8 meshs on the thickness, 10 meshs over the length)

5.3 Functionalities

tested

Controls

Keys

THER_LINE

PARM_THETA

[U4.33.01]

DEF_MATERIAU

ECRO_LINE_FO

D_SIGM_EPSI

[U4.23.01]

STAT_NON_LINE

COMP_INCR

DEFORMATION

PETIT_REAC

[U4.32.01]

RELATION

VMIS_ISOT_LINE

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

10/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Results of modeling B

6.1 Values

tested

Identification

Moment

Reference

Aster

% difference

Temperature (N1)

0.5s

1.25355E1

1.253935E1

0.031

Temperature (N161)

0.5s

5.712330E0

5.710409E0

0.034

Temperature (N162)

0.5s

2.445260E0

2.446690E0

0.058

Temperature (N163)

0.5s

1.121890E0

1.121793E0

0.009

Temperature (N164)

0.5s

4.836440E1

4.839514E1

0.064

Temperature (N165)

0.5s

2.224430E1

2.224511E1

0.004

Temperature (N166)

0.5s

9.620360E2

9.626885E2

0.068

Temperature (N167)

0.5s

4.434440E2

4.434991E2

0.012

Temperature (N168)

0.5s

1.923100E2

1.924445E2

0.070

Temperature (N169)

0.5s

8.879760E3

8.881443E3

0.019

Temperature (N170)

0.5s

3.859600E3

3.862363E3

0.072

Temperature (N171)

0.5s

1.785240E3

1.785670E3

0.024

Temperature (N172)

0.5s

7.786500E4

7.792119E3

0.072

Temperature (N173)

0.5s

3.631140E4

3.632162E4

0.028

Temperature (N174)

0.5s

1.644480E4

1.645681E4

0.073

Temperature (N175)

0.5s

8.805770E5

8.808440E5

0.030

Temperature (N14)

0.5s

6.549040E5

6.553566E5

0.069

Temperature (N1)

25s

8.642670E1

8.642786E1

0.001

Temperature (N161)

25s

7.656950E1

7.656903E1

0.01

Temperature (N162)

25s

6.773550E1

6.773627E1

0.001

Temperature (N163)

25s

5.986100E1

5.986069E1

0.001

Temperature (N164)

25s

5.284760E1

5.284816E1

0.001

Temperature (N165)

25s

4.664620E1

4.664591E1

0.001

Temperature (N166)

25s

4.118010E1

4.118056E1

0.001

Temperature (N167)

25s

3.641060E1

3.641040E1

0.001

Temperature (N168)

25s

3.227650E1

3.227691E1

0.001

Temperature (N169)

25s

2.874680E1

2.874659E1

0.001

Temperature (N170)

25s

2.577230E1

2.577260E1

0.001

Temperature (N171)

25s

2.332830E1

2.332814E1

0.001

Temperature (N172)

25s

2.137660E1

2.137655E1

0.001

Temperature (N173)

25s

1.989790E1

1.989793E1

Temperature (N174)

25s

1.886470E1

1.886472E1

0.001

Temperature (N175)

25s

1.824500E1

Temperature (N14)

25s

1.806070E1

1.806074E1

0.001

SIYY (M90, N222)

15s

1.204720E8

1.200544E8

0.347

SIYY (M101, N164)

15s

1.659340E8

1.661848E8

0.151

SIYY (M90, N222)

20s

1.391640E8

1.383875E8

0.558

SIYY (M101, N164)

20s

1.640310E8

1.643171E8

0.174

SIYY (M90, N222)

25s

1.284210E8

1.275833E8

0.652

SIYY (M101, N204)

25s

1.568610E8

1.570781E8

0.138

6.2 Parameters

of execution

Version: 3.09
Machine:

System:

Overall dimension memory:

16 megawords

Time CPU To use: 58.37 seconds

Code_Aster

Version

4.0

Titrate:

Thermal SSNV301 Shock and mechanical loading

Date:

01/12/98

Author (S)

E. LORENTZ, L. LAMMERANT, H. DRION

Key:

V6.04.301-A

Page:

11/12

Manual of Validation

V6.04 booklet: Nonlinear statics of the voluminal structures

HI-75/98/040 - Ind A

Summary of the results

It was checked that one obtains the same results with or without recovery, with or without amendment of
boundary conditions.

It should be noted that the nodal efforts that one provides to Aster are divided by 2

compared to those which

one provides to the SAMCEF software into axisymmetric.