Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
1/8

Organization (S): EDF-R & D/AMA

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

SSNL122 - Poutre Multi-Fibers cantilever subjected
with an effort

Summary:

This test relates to the validation of beam multifibre with a modeling in POU_D_TGM.

This problem makes it possible to test:

·
linear finite elements of beams type with a modeling in POU_D_TGM,
·
the taking into account of the orientation,
·
the calculation of the SIEF_ELGA and the SIEF_ELNO_ELGA.

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
2/8

1
Problem of reference

1.1 Geometry

Z
X
B
Y
C
Lo
With


Length of the bar: 3m
Section of the bar:
Embedding in A
height: 0.04m
Forces out of B
width: 0.02m

1.2
Properties of material

Material with elastoplastic behavior with a linear work hardening, for the linear element:

Elasticity:
·
Young modulus E = 2.1E+11 Pa
Plasticity:
·
slope of the traction diagram in the plastic range D_SIGM_EPSI = 1.0E+08 Pa
·
yield stress SY = 400.0E+06 Pa

1.3
Boundary conditions and loadings

At the point A, perfect embedding (blocking of displacements and rotations),
·
blocking of the ddl: dx, Dy, dz, DRX, DRY, DRZ.
Loading at point b: F = (X
F, Fy, Z
F).

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
3/8

2
Reference solution

2.1
Sizes and results of reference

Arrow out of B following Z
Arrow out of B following Y
Z
F L
. 3

Fy L
. 3
Z =

=

E
.
3
I
.
y
y
E
.
3
I
. Z

Constraint in a point C of co-ordinates (vy, vz) of the section of the beam

NR
M
M = +F Lo
.
y
M Z
=
+
v
.
-
v
.
Z
y
Z
y
with

S
I
I
y
Z
M = - F Lo
y
Z.

from where:
NR
ELo
.
3
=
-
(.y. vy+z.vz)
éq
2.1-1
S
L3

2.2 Reference
bibliographical

[1]
“M7-01-72 Project. Elastoplastic behavior of the beams. New approach.” Note
HM77/01/140/A.
Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
4/8

3 Modeling
With

3.1
Characteristics of modeling and the grid

Linear element: POU_D_TGM.

E00200
E00201
E00202
E00203
E00204
E00100
E00101
E00102
S002: 2m
S001: 1m


Mechanical characteristics of the section (homogeneous units with [m]):

WITH IY
IZ AY AZ JX JG
8.0e-04 2.666667e-08 1.066667e-07 1.191790e+0
1.172840e+0 7.093682e-08 1.438125e-12
0
0

Loading at the point B.


Fx Fy Fz
Moment 1
80.000N
­ 150N
­ 200N
Moment 2
80.000N
­ 280N
­ 400N

At moment 1 the section remains elastic, at moment 2 the section is partially plasticized.

3.2
Classification of fibers

The principal reference mark of inertia (Y0, Z0) of the beam must be turned of 90° so that the strongest inertia
“IZ” is along the axis “Y” of the total reference mark. The goal is to test key word ORIENTATION of
order AFFE_CARA_ELEM.

In the local reference mark of the beam
In the total reference mark of study


Y0 Y

35
36
49
52
Z0
48
32
31
30
29
53
34
35
41
48
4
32
40
49
39
3
31
Y0
38
2
6
30
Z
Z0
37
1
5
29
52
33
60
53
36

6
5
41
4
3
2
1
60
34
33
40
39
38
37


Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
5/8

Several “types” of fibers are used:

·
fibers whose numbers go from 1 to 32. They are affected of a nonnull surface in
command file,
·
fibers whose numbers are 33, 34, 35, 36. They are located at the 4 corners of the section. These
fibers, in the command file, are affected of a section equal to zero. They
are useful that with postprocessing,
·
fibers whose numbers go from 37 to 60. They are located on the edge external of the section.
These fibers, in the command file, are affected of a section equal to zero. They
are used that for postprocessing.

3.3 Functionalities
tested

Commands


DEFI_MODELE POU_D_TGM

AFFE_CARA_ELEM BEAM
AFFE_FIBER
ORIENTATION
IMPR_RESU
SIEF_ELGA
SIEF_ELNO_ELGA

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
6/8

4
Results of modeling A

4.1
Sizes tested and results

The sizes tested and analyzed are them:
· SIEF_ELGA, at the first point of Gauss of the E00200 element. It is the point of Gauss it
more close to embedding, Lo = 2.95491933m.
· SIEF_ELNO_ELGA, with the 2 nodes of the E00200 element. The G00202 node is the point
of embedding (Lo = 3.0m), the I00200 node is in Lo = 2.6m.

4.1.1 Behavior
rubber band

Constraints at the point of Gauss: SIEF_ELGA
The constraints calculated by the equation [éq 2.1-1] and by Code_Aster are given in both
tables below (values in MPa). The provision of the tables takes again the diagram of provision
fibers in the total reference mark. The most requested fiber is the n°36 with a constraint of 390 MPa

Constraints calculated in

Constraints calculated in
fibers by the equation [éq 2.1-1]
fibers by Code_Aster.

35
80.168.257
346.390 35
80
168
257.346.390
21
66.154.243
331.376 21
66
154
243.331.376
- 7
37.126.215
303.348 - 7
37
126
215.303.348
- 35
9 98 187
275.320 - 35
9
98
187.275.320
- 63 - 19 70 158
247.291 - 63
- 19
70
158.247.291
- 91 - 47 42 130
219.263 - 91
- 47
42
130.219.263
- 120 - 75 13 102
191.235 - 120
- 75
13
102.191.235
- 148 - 103 - 15 74
163.207 - 148 - 103
- 15
74.163.207
- 176 - 131 - 43 46
134.179 - 176 - 131
- 43
46.134.179
- 190 - 146 - 57 32
120.165 - 190 - 146
- 57
32.120.165

The relative error between two calculations is given in the table below. Whatever the fiber,
it remains lower than 0.1%.

- 0.022% - 0.013% - 0.009% - 0.008% - 0.007% - 0.007%
- 0.026% - 0.012% - 0.009% - 0.008% - 0.007% - 0.007%
0.014% - 0.010% - 0.007% - 0.006% - 0.006% - 0.006%
- 0.010% 0.010% - 0.004% - 0.005% - 0.005% - 0.005%
- 0.013% - 0.029% 0.000% - 0.003% - 0.004% - 0.004%
- 0.014% - 0.021% 0.012% 0.000% - 0.002% - 0.003%
- 0.014% - 0.019% 0.069% 0.004% - 0.001% - 0.002%
- 0.015% - 0.018% - 0.094% 0.012% 0.002% 0.000%
- 0.015% - 0.018% - 0.043% 0.029% 0.006% 0.003%
- 0.015% - 0.018% - 0.036% 0.049% 0.009% 0.005%

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
7/8

Constraints with the nodes: SIEF_ELNO_ELGA
The constraints calculated by the equation [éq 2.1-1] and by Code_Aster are given in the table
below (values in MPa)

Analytical node Code_Aster
Relative error
G00202
SIMAX = 394.307
SIMAX = 394.276
- 0.008%

SIMIN = - 194.307
SIMIN = - 194.276
- 0.016%
I00200
SIMAX = 355.066
SIMAX = 355.057
- 0.002%

SIMIN = - 155.066
SIMIN = - 155.057
- 0.006%

4.1.2 Behavior
plastic

Constraints at the point of Gauss: SIEF_ELGA
The table below gives the values of the constraints, after plasticization partial of the section,
obtained with Code_Aster. The behavior of material is elastoplastic “almost perfect”,
slope of work hardening is weak. The maximum constraint, which is beyond the yield stress, remains
thus very close to the elastic threshold of 400MPa.

8
79
223
366
400
400
- 16
55
199
342
400
400
- 64
7
151
294
400
400
- 112
- 41
103
246
390
400
- 160
- 89
55
198
342
400
- 208 - 137
7
150
294
365
- 256 - 185
- 41
102
246
317
- 305 - 233
- 89
54
198
269
- 353 - 281 - 137
6
150
221
- 377 - 305 - 161
- 18
126
197

This calculation is carried out for the test of nonregression of Code_Aster.

Constraints with the nodes: SIEF_ELNO_ELGA
The table below gives the values of the constraints to the nodes, after plasticization partial of
section, obtained with Code_Aster. These values are interpolated by Code_Aster from
SIEF_ELGA. The behavior of material is elastoplastic “almost perfect”, the slope
of work hardening is weak. The maximum constraint, which is beyond the yield stress, thus remains very
neighbor of the elastic threshold of 400MPa.

Node
Code_Aster
G00202
SIMAX = 400.091

SIMIN = - 385.441
I00200
SIMAX = 400.053

SIMIN = - 306.905

This calculation is carried out for the test of nonregression of Code_Aster.

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Code_Aster ®
Version
7.2
Titrate:
SSNL122 - Poutre Multi-Fibers cantilever subjected to an effort
Date:
01/03/04
Author (S):
J.L. Key FLEJOU
:
V6.02.122-B Page:
8/8

5
Summary of the results

This case test shows the correct operation of a modeling of the behavior of the beams by one
approach multifibre. The postprocessing of the beams does not require any more, with this modeling, of
to calculate the SIGM_ELNO_SIEF or the SIPO_ELNO_SIEF.

Handbook of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HT-66/04/005/A

Outline document