Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
1/16

Organization (S): EDF/MTI/MN, DeltaCAD

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
Document: V6.02.501

SSNL501 - Poutre embedded at the two ends
subjected to a uniform pressure

Summary:

This test represents a quasi-static calculation of a fixed beam subjected to a uniform pressure,
composed of a perfectly plastic elastic material. This test makes it possible to validate modelings elements
stop following:

·
COQUE_C_PLAN (SEG3),
·
DKT (TRIA3 and QUAD4),
·
COQUE_3D (TRIA7 and QUAD9),
·
POU_D_T (SEG2).

The limiting pressure is compared with an analytical reference solution.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
2/16

1
Problem of reference

1.1 Geometry
y
Z
p
B

F
, y
C
, Z
H
B
X
L = 1.5 m
With

E
D
B = 0.15m
, X
L
H = 0.025 m


1.2
Properties of material

The properties of material constituting the beam are:

E = 2. 1011 Pa
Young modulus
= 0.3
Poisson's ratio

The material follows a perfectly plastic law of elastic behavior:

= 2.35 108 Pa
Yield stress
E
= 1.175 10-3
Limiting elastic strain
E

1.3
Boundary conditions and loadings

-
C.L. : Dimensioned embedded AB and CD

-
Imposed displacement following Z in E (x=L/2):

Q L4
DZ
E
=
E = 6.60910-3 m (DZ
)
E
384EI

DZ (E) varies from 0 to 30 DZe

1.4 Conditions
initial

Without object
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
3/16

2
Reference solution

2.1
Method of calculation used for the reference solution

Q = p B
E
With
D
L
(A)
Q
E
With
D
M
With
E
M
p
p
L/2
(b)
(c)

The ruin of the beam appears when there are plastic knee joints at points A, D and E (figure b).
The static balance of the left half of the beam, makes it possible to determine the pressure limits (figure c)

L L

16M
M = 2 M - Q () () = 0
p
=> Q =

With
p
L
2 4
L
L2
where:
qL represents the limiting pressure
bh2
M
=
p
represent the plastic moment (M p
E
4)

The appearance of the first plastic point on fiber external of the beam has place at points A and D, them
bh2
other fibers being in elastic mode. The pressure limits elastic is of Q = 2
E
E
.
L2

2.2
Results of reference

Pressure limits qL= 39 166.67 NR/m

2.3
Uncertainties on the solution

Analytical solution

2.4 References
bibliographical

[1]
WILLIAM A. NASH: Theory and problems off Strength off material, Schaum' S outline series,
2/ED, McGRAW-HILL
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
4/16

3 Modeling
With

3.1
Characteristics of modeling

E
y
Modeling COQUE_C_PLAN
Not A (2.1; 0.7) (node n1)
78.5°
With
Boundary conditions: Not a: U = v = z= 0
Conditions of symmetry: Not E: z= 0
Z
X
the beam being tilted with 78.5°, the value of displacement
imposed is then: DXe=6.609 10-3 Mr. sin (78.5°)


3.2
Characteristics of the grid

A number of nodes: 21
A number of meshs and type: 10 SEG3

3.3 Functionalities
tested

Commands Key word
factor
Key word
AFFE_MODELE AFFE
MODELISATION
= ' COQUE_C_PLAN'
DEFI_MATERIAU ECRO_LINE
D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
EPAIS
A_CIS =0.833333
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR
VMIS_ISOT_LINE

COQUE_NCOU
PILOTAGE
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”

4
Results of modeling A

4.1 Values
tested

DX (E) DX Identification Moments Reference
Aster %
difference
E
10
ETA_PILOTAGE
10.1.0 1.0699
6.9
12
ETA_PILOTAGE
12.1.0 1.0699
6.9

4.2 Remarks

The width for modeling COQUE_C_PLAN is imposed on 1 in Code_Aster. In
consequence, we multiplied by 1/B the loading to take account of the real width of
beam.
In this analysis one uses to find the solution, a technique of resolution of the type
imposed displacement (“DDL_IMPO”). This method provides for each value of displacement
imposed, a multiplying coefficient of the loading (“ETA_PILOTAGE”). The value of the loading
imposed in “AFFE_CHAR_MECA” is equal to the limiting pressure, consequently the value of
reference of parameter “ETA_PILOTAGE” is equal to 1.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
5/16

4.3
Influence refinement of the grid

EDF
Department Mecanique and Modeles Numeriques
Electricity
Eta according to displacement
from France
INFLUENCE GRID
ON SOLUTION COQUE1D
COMPARISON WITH
1.0
THE SOLUTION BEAM
0.8
AGE
0.6
OT
I
L
P
_
With
T
COQUE1D mall 20
E
COQUE1D mall 10
POUTRE mall 20
0.4
0.2
0.0
0
2
4
6
8
10x10-2
displacement (m)
agraf 12/04/2001 (c) EDF/DER 1992-1999

The reference being taken compared to the solution beam with 20 elements, one observes one
improvement of results COQUE_C_PLAN when the grid is refined.

4.4 Parameters
of execution

Version:
NEW 5.04.19

Machine:
SGI-Origin2000 R12000

Obstruction memory: 16 megabytes
Time CPU To use: 16.26 seconds
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
6/16

5 Modeling
B

5.1
Characteristics of modeling

Z
zl
Modeling DKT (TRIA3)
yl
F
OE = OA = L/2 2
xl
AB//EFF//axis Z
O
Boundary conditions (reference mark gobal)
y
E
- side AB: U = v =w =
B
x= y= z= 0
Conditions of symmetry
- side EFF: U = 0 (local reference mark x1 y1 z1)
- side EFF:
X
z= 0 (total reference mark)
With
the beam being tilted with 45°, the value of displacement
imposed is then: DXe=6.609 10-3 Mr. sin (45°)

5.2
Characteristics of the grid

A number of nodes: 43
A number of meshs and type: 20 TRIA3

5.3 Functionalities
tested

Commands Key word
factor
Key word
AFFE_MODELE AFFE
MODELISATION
= ' DKT'
DEFI_MATERIAU ECRO_LINE
D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
EPAIS
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR
VMIS_ISOT_LINE

COQUE_NCOU
PILOTAGE
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”

6
Results of modeling B

6.1 Values
tested

DX (E) DX Identification Moments Reference
Aster %
difference
E
5
ETA_PILOTAGE
5.1.0 1.11133 11.11
25
ETA_PILOTAGE
15 1.0
1.142 11.42

6.2 Remarks

In this analysis, one uses to find the solution, a technique of resolution of the type
imposed displacement (“DDL_IMPO”). This method provides for each value of displacement
imposed, a multiplying coefficient of the loading (“ETA_PILOTAGE”). The value of the loading
imposed in “AFFE_CHAR_MECA” is equal to the limiting pressure, consequently the value of
reference of parameter “ETA_PILOTAGE” is equal to 1.
Calculations were stopped when the value of parameter ETA_PILOTAGE was stabilized.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
7/16

6.3
Influence refinement of the grid

EDF
Department Mecanique and Modeles Numeriques
Electricity
Eta according to displacement
from France
1.2
INFLUENCE GRID
ON MODELING
DKT TRIA3
COMPARISON WITH
1.0
THE SOLUTION BEAM
0.8
AGE
OT
I
L
P
0.6
_
With
T
DKT 10
E
DKT 20
DKT 50
DKT 100
POU_D_T 50
0.4
0.2
0.0
0
2
4
6
8
10
12
14x10-2
displacement (m)
agraf 12/04/2001 (c) EDF/DER 1992-1999


The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of results DKT TRIA3 when the grid is refined.

6.4 Parameters
of execution

Version:
NEW 5.04.19

Machine:
SGI-Origin2000 R12000

Obstruction memory: 16 megabytes
Time CPU To use: 20.48 seconds
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
8/16

7 Modeling
C

7.1
Characteristics of modeling

Z
Modeling DKT (QUAD4)
zl
yl
OE = OA = L/2 2
F
AB//EFF//axis Z
xl
Boundary conditions (reference mark gobal)
O
- side AB: U = v =w =
y
x= y= z= 0
E
B
Conditions of symmetry
- side EFF: U = 0 (local reference mark x1 y1 z1)
- side EFF: z= 0 (total reference mark)
X
With
the beam being tilted with 45°, the value of displacement
imposed is then: DXe=6.609 10-3 Mr. sin (45°)

7.2
Characteristics of the grid

A number of nodes: 43
A number of meshs and type: 10 QUAD4

7.3 Functionalities
tested

Commands Key word
factor
Key word
AFFE_MODELE AFFE
MODELISATION
= ' DKT'
DEFI_MATERIAU ECRO_LINE D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
EPAIS
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR VMIS_ISOT_LINE

COQUE_NCOU
PILOTAGE
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”

8
Results of modeling C

8.1 Values
tested

DX (E) DX Identification Moments Reference
Aster %
difference
E
5
ETA_PILOTAGE
5.1.0 1.0837 8.37
25
ETA_PILOTAGE
25.1.0 1.0998 9.98

8.2 Remarks

In this analysis, one uses to find the solution, a technique of resolution of the type
imposed displacement (“DDL_IMPO”). This method provides for each value of displacement
imposed, a multiplying coefficient of the loading (“ETA_PILOTAGE”). The value of the loading
imposed in “AFFE_CHAR_MECA” is equal to the limiting pressure, consequently the value of
reference of parameter “ETA_PILOTAGE” is equal to 1.
Calculations were stopped when the value of parameter ETA_PILOTAGE was stabilized.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
9/16

8.3
Influence refinement of the grid

EDF
Department Mecanique and Modeles Numeriques
Electricity
Eta according to displacement
from France
1.2
INFLUENCE GRID
ON MODELING
DKT (QUAD4)
COMPARISON WITH
1.0
THE SOLUTION BEAM
0.8
AGE
OT
I
L
P
0.6
_
With
T
DKT 10 QUAD4
E
DKT 20 QUAD4
DKT 50 QUAD4
DKT 100 QUAD4
BEAM POU_D_T 50
0.4
0.2
0.0
0
2
4
6
8
10
12
14x10-2
displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999


The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of results DKT QUAD4 when the grid is refined.

8.4 Parameters
of execution

Version:
NEW 5.04.19

Machine:
SGI-Origin2000 R12000

Obstruction memory: 16 megabytes
Time CPU To use: 27.55 seconds
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
10/16

9 Modeling
D

9.1
Characteristics of modeling

Z
Modeling COQUE_3D (TRIA7)
zl
yl
OE = OA = L/2 2
F
AB//EFF//axis Z
xl
Boundary conditions (reference mark gobal)
O
- side AB: U = v =w =
y
x= y= z= 0
E
B
Conditions of symmetry
- side EFF: U = 0 (local reference mark x1 y1z1)
- side EFF: z= 0 (total reference mark)
X
With
the beam being tilted with 45°, the value of displacement
imposed is then: DXe=6.609 10-3 Mr. sin (45°)


9.2
Characteristics of the grid

A number of nodes: 83
A number of meshs and type: 20 TRIA7

9.3 Functionalities
tested

Commands Key word
factor
Key word
AFFE_MODELE AFFE
MODELISATION
= ' COQUE_3D'
DEFI_MATERIAU ECRO_LINE D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
EPAIS
A_CIS =0
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR VMIS_ISOT_LINE

COQUE_NCOU
PILOTAGE
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”

10 Results of modeling D
10.1 Values
tested

DX (E) DX Identification Moments Reference
Aster %
difference
E
5
ETA_PILOTAGE
5.1.0 1.1143 11.43
15
ETA_PILOTAGE
15.1.0 1.1682 16.82

10.2 Remarks

In this analysis, one uses to find the solution, a technique of resolution of the type
imposed displacement (“DDL_IMPO”). This method provides for each value of displacement
imposed, a multiplying coefficient of the loading (“ETA_PILOTAGE”). The value of the loading
imposed in “AFFE_CHAR_MECA” is equal to the limiting pressure, consequently the value of
reference of parameter “ETA_PILOTAGE” is equal to 1.
Calculations were stopped when the value of parameter ETA_PILOTAGE was stabilized.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
11/16

10.3 Influence refinement of the grid

EDF
Department Mecanique and Modeles Numeriques
Electricity
Eta according to displacement
from France
1.2
INFLUENCE GRID
ON MODELING
COQUE_3D TRIA7
COMPARISON WITH
1.0
THE SOLUTION BEAM
0.8
AGE
OT
I
L
P
0.6
_
With
T
10 COQUE_3D TRIA7
E
20 COQUE_3D TRIA7
50 COQUE_3D TRIA7
100 COQUE_3D TRIA7
BEAM POU_D_T 50
0.4
0.2
0.0
0
2
4
6
8
10
12
14x10-2
displacement (m)
agraf 12/04/2001 (c) EDF/DER 1992-1999


The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of the results hulls 3D TRIA7 when the grid is refined.

10.4 Parameters
of execution

Version:
NEW 5.04.19

Machine:
SGI-Origin2000 R12000

Obstruction memory: 16 megabytes
Time CPU To use: 72.54 seconds

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
12/16

11 Modeling
E

11.1 Characteristics of modeling

Z
zl
Modeling COQUE_3D (QUAD9)
yl
F
OE = OA = L/2 2
xl
AB//EFF//axis Z
O
Boundary conditions (reference mark gobal)
y
E
- side AB: U = v =w =
B
x= y= z= 0
Conditions of symmetry
- side EFF: U = 0 (local reference mark x1 y1 z1)
- side EFF:
X
z= 0 (total reference mark)
With
the beam being tilted with 45°, the value of displacement
imposed is then: DXe=6.609 10-3 Mr. sin (45°)

11.2 Characteristics of the grid

A number of nodes: 54
A number of meshs and type: 10 QUAD9

11.3 Functionalities
tested

Commands Key word
factor
Key word
AFFE_MODELE AFFE
MODELISATION
= ' COQUE_3D'
DEFI_MATERIAU ECRO_LINE D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
EPAIS
A_CIS = 105
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR VMIS_ISOT_LINE

COQUE_NCOU
PILOTAGE
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”

12 Results of modeling E

12.1 Values
tested

DX (E) DX Identification Moments Reference
Aster %
difference
E
5
ETA_PILOTAGE
5.1.0 1.0978 9.78
25
ETA_PILOTAGE
25.1.0 1.1085 10.85

12.2 Remarks

In this analysis, one uses to find the solution, a technique of resolution of the type
imposed displacement (“DDL_IMPO”). This method provides for each value of displacement
imposed, a multiplying coefficient of the loading (“ETA_PILOTAGE”). The value of the loading
imposed in “AFFE_CHAR_MECA” is equal to the limiting pressure, consequently the value of
reference of parameter “ETA_PILOTAGE” is equal to 1.
Calculations were stopped when the value of parameter ETA_PILOTAGE was stabilized.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
13/16

12.3 Influence refinement of the grid

EDF
Department Mecanique and Modeles Numeriques
Electricity
Eta according to displacement
from France
1.2
INFLUENCE GRID
ON MODELING
COQUE_3D (QUAD 9)
1.0
COMPARISON WITH
SOLUTION BEAM
0.8
AGE
OT
I
L
P
0.6
_
With
T
COQUE_3D 10 QUAD9
E
COQUE_3D 20 QUAD9
COQUE_3D 50 QUAD9
BEAM POU_D_T 50
COQUE_3D 100 QUAD9
0.4
0.2
0.0
0
2
4
6
8
10
12
14x10-2
displacement (m)
agraf 12/04/2001 (c) EDF/DER 1992-1999


The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of the results hulls 3D QUAD9 when the grid is refined.

12.4 Parameters
of execution

Version:
NEW 5.04.19

Machine:
SGI-Origin2000 R12000

Obstruction memory: 16 megabytes
Time CPU To use: 74.93 seconds
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
14/16

13 Modeling
F

13.1 characteristics of modeling

Modeling POU_D_T
E
y
A= 2.1 0.7
0.
AE =0.75 m
78.5°
With
O
Boundary conditions (reference mark gobal)
X
- node a: U = v =w = qx= qy= qz= 0
Conditions of symmetry
- node E: U = 0, Q
Z
z= 0
the beam being tilted with 78.5°, the value of displacement
imposed is then: DXe=6.609 10-3 Mr. sin (78.5°)

13.2 characteristics of the grid

A number of nodes: 21
A number of meshs and type: 10 SEG2

13.3 Functionalities
tested

Commands Key word
factor Word-key
AFFE_MODELE AFFE
MODELISATION
= ' POU-D_T'
DEFI_MATERIAU ECRO_LINE
D_SIGM_EPSI

SY
VMIS_POUTRE

AFFE_CARA_ELEM BEAM
SECTION
CARA
AFFE_CHAR_MECA FORCE_POUTRE
FY
STAT_NON_LINE COMP_INCR
VMIS_POU_LINE
PILOTAGE
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”

14 Results of modeling F

14.1 Values
tested

DX (E) DX Identification Moments Reference
Aster %
difference
E
5
ETA_PILOTAGE
5.1.0 1.1083 10.83
25
ETA_PILOTAGE
25.1.0 1.1107 11.07

In this analysis, one uses to find the solution, a technique of resolution of the type
imposed displacement (“DDL_IMPO”). This method provides for each value of displacement
imposed, a multiplying coefficient of the loading (“ETA_PILOTAGE”). The value of the loading
imposed in “AFFE_CHAR_MECA” is equal to the limiting pressure, consequently the value of
reference of parameter “ETA_PILOTAGE” is equal to 1.
Calculations were stopped when the value of parameter ETA_PILOTAGE was stabilized.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
15/16

14.2 Influence refinement of the grid

EDF
Department Mecanique and Modeles Numeriques
Electricity
Eta according to displacement
from France
1.2
INFLUENCE GRID
ON MODELING
BEAM POU_D_T
1.0
0.8
AGE
OT
I
L
P
0.6
_
With
T
POU_D_T 10
E
POU_D_T 20
POU_D_T 100
POU_D_T 50
0.4
0.2
0.0
0
2
4
6
8
10x10-2
displacement (m)
agraf 12/04/2001 (c) EDF/DER 1992-1999


the solution beam improves appreciably when the grid is refined.

14.3 Parameters
of execution

Version:
NEW 5.04.19

Machine:
SGI-Origin2000 R12000

Obstruction memory: 16 megabytes
Time CPU To use: 74.93 seconds
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Code_Aster ®
Version
5.0
Titrate:
SSNL501 - Embedded Poutre subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER Key
:
V6.02.501-A Page:
16/16

15 Summary of the results

On the figure below, we present the evolution of normal displacement at the center of the plate in
function of the parameter of control.

EDF
Department Mecanique and Modeles Numeriques
Electricity
Eta according to displacement
from France
1.2
COMPARISON OF THE SOLUTIONS
1.0
0.8
E ETA
ETR
M
0.6
With
R
Modeling A
With
P
Modeling B
Modeling C
Modeling D
Modeling E
0.4
Modeling F
0.2
0.0
0
2
4
6
8
10
12
14
16x10-2
displacement (m)
agraf 12/04/2001 (c) EDF/DER 1992-1999


It is observed that:

·
modelings comprising of the quadrangles give better results, by report/ratio
with the meshs triangles.

Calculations were stopped when the value of parameter ETA_PILOTAGE was stabilized, or
when calculation was not possible any more. Taking into account the grids used, the results obtained are
satisfactory. The errors observed are for modelings:

·
COQUE_C_PLAN: 6.9% (A),
·
DKT:11.4% for mesh TRIA3 (B) and 9.9% for mesh QUAD4 (C),
·
COQUE_3D: 16% for mesh TRIA7 (D) and 10.8% for mesh QUAD9 (E),
·
POUT_D_T:11% (F).

But one notes that with a finer grid at the ends and the center of the plate, place or
plasticization appears, it is possible to minimize the error compared to the reference solution.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A

Outline document