Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
1/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
Organization (S):
EDF/MTI/MN, DeltaCAD
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
Document: V6.02.501
SSNL501 - Beam fixed 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.
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
2/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
1
Problem of reference
1.1 Geometry
With
B
D
C
Z
y
X
H
p
L
L = 1.5 m
B = 0.15m
H = 0.025 m
B
, X
, Z
, y
E
F
1.2
Properties of material
The properties of material constituting the beam are:
E
= 2. 10
11
AP
Young modulus
= 0.3
Poisson's ratio
The material follows a perfectly plastic law of elastic behavior:
E
= 2.35 10
8
AP
Yield stress
E
= 1.175 10
- 3
Limiting elastic strain
1.3
Boundary conditions and loadings
-
C.L. : Dimensioned embedded AB and CD
-
Imposed displacement following Z in E (x=L/2):
DZ
E
= 6.60910
- 3
m (
DZ
Q L
I.E.(internal excitation)
E
E
=
4
384
)
DZ (E) varies from 0 to 30 DZ
E
1.4 Conditions
initial
Without object
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
3/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
2
Reference solution
2.1
Method of calculation used for the reference solution
With
D
Q = p B
L
With
D
M
p
M
p
L/2
Q
(A)
(b)
(c)
E
E
With
E
The ruin of the beam appears when there are plastic ball 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)
M
M
Q
L L
With
p
L
=
-
=
2
2 4
0
() ()
=>
Q
M
L
L
p
= 16
2
where:
Q
L
represent the limiting pressure
M
p
represent the plastic moment (
M
bh
p
E
=
2
4
)
The appearance of the first plastic point on fiber external of the beam has place at points A and D, them
other fibers being in elastic mode. The pressure limits elastic is of
Q
bh
L
E
E
= 2
2
2
.
2.2
Results of reference
Pressure limits Q
L
= 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
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
4/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
3 Modeling
With
3.1
Characteristics of modeling
Modeling COQUE_C_PLAN
Not A (2.1; 0.7) (node n1)
E
With
y
X
Z
78.5°
Boundary conditions: Not a: U = v =
Z
=
0
Conditions of symmetry: Not E
:
Z
=
0
the beam being tilted with 78.5°, the value of displacement
imposed is then: DX
E
=6.609 10
- 3
Mr. sin (78.5°)
3.2
Characteristics of the mesh
A number of nodes: 21
A number of meshs and type: 10 SEG3
3.3 Functionalities
tested
Controls Key word
factor
Key word
AFFE_MODELE AFFE
MODELING
= ' COQUE_C_PLAN'
DEFI_MATERIAU ECRO_LINE
D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
THICK
A_CIS =0.833333
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR
PILOTING
VMIS_ISOT_LINE
COQUE_NCOU
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”
4
Results of modeling A
4.1 Values
tested
DX E DX
E
()
Identification Moments Reference
Aster %
difference
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.
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
5/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
4.3
Influence refinement of the mesh
The reference being taken compared to the solution beam with 20 elements, one observes one
improvement of results COQUE_C_PLAN when the mesh is refined.
4.4 Parameters
of execution
Version:
NEW 5.04.19
Machine:
SGI-Origin2000 R12000
Overall dimension memory: 16 megabytes
Time CPU To use: 16.26 seconds
displacement (m)
E
T
With
_
P
I
L
OT
AGE
EDF
Electricity
from France
Mechanical department and Numerical Models
Eta according to displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999
COQUE1D mall 20
COQUE1D mall 10
BEAM mall 20
INFLUENCE MESH
ON SOLUTION COQUE1D
COMPARISON WITH
THE SOLUTION BEAM
0.0
0.2
0.4
0.6
0.8
1.0
0
2
4
6
8
10x10-2
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
6/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
5 Modeling
B
5.1
Characteristics of modeling
Modeling DKT (TRIA3)
OE = OA =
L/2 2
AB//EFF//axis Z
y
X
Z
With
E
F
B
O
Boundary conditions (reference mark gobal)
- side AB: U = v =w =
X
=
y
=
Z
=
0
Conditions of symmetry
- side EFF: U =
0 (local reference mark X
1
y
1
Z
1
)
- side EFF:
Z
=
0 (total reference mark)
X
L
y
L
Z
L
the beam being tilted with 45°, the value of displacement
imposed is then: DX
E
=6.609 10
- 3
Mr. sin (45°)
5.2
Characteristics of the mesh
A number of nodes: 43
A number of meshs and type: 20 TRIA3
5.3 Functionalities
tested
Controls Key word
factor
Key word
AFFE_MODELE AFFE
MODELING
= ' DKT'
DEFI_MATERIAU ECRO_LINE
D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
THICK
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR
PILOTING
VMIS_ISOT_LINE
COQUE_NCOU
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”
6
Results of modeling B
6.1 Values
tested
DX E DX
E
()
Identification Moments Reference
Aster %
difference
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.
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
7/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
6.3
Influence refinement of the mesh
displacement (m)
E
T
With
_
P
I
L
OT
AGE
EDF
Electricity
from France
Mechanical department and Numerical Models
Eta according to displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999
DKT 10
DKT 20
DKT 50
DKT 100
POU_D_T 50
INFLUENCE MESH
ON MODELING
DKT TRIA3
COMPARISON WITH
THE SOLUTION BEAM
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0
2
4
6
8
10
12
14x10-2
The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of results DKT TRIA3 when the mesh is refined.
6.4 Parameters
of execution
Version:
NEW 5.04.19
Machine:
SGI-Origin2000 R12000
Overall dimension memory: 16 megabytes
Time CPU To use: 20.48 seconds
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
8/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
7 Modeling
C
7.1
Characteristics of modeling
Modeling DKT (QUAD4)
OE = OA =
L/2 2
AB//EFF//axis Z
y
X
Z
With
E
F
B
O
Boundary conditions (reference mark gobal)
- side AB: U = v =w =
X
=
y
=
Z
=
0
Conditions of symmetry
- side EFF: U =
0 (local reference mark X
1
y
1
Z
1
)
- side EFF:
Z
=
0 (total reference mark)
X
L
y
L
Z
L
the beam being tilted with 45°, the value of displacement
imposed is then: DX
E
=6.609 10
- 3
Mr. sin (45°)
7.2
Characteristics of the mesh
A number of nodes: 43
A number of meshs and type: 10 QUAD4
7.3 Functionalities
tested
Controls Key word
factor
Key word
AFFE_MODELE AFFE
MODELING
= ' DKT'
DEFI_MATERIAU ECRO_LINE D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
THICK
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR
PILOTING
VMIS_ISOT_LINE
COQUE_NCOU
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”
8
Results of modeling C
8.1 Values
tested
DX E DX
E
()
Identification Moments Reference
Aster %
difference
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.
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
9/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
8.3
Influence refinement of the mesh
displacement
E
T
With
_
P
I
L
OT
AGE
EDF
Electricity
from France
Mechanical department and Numerical Models
Eta according to displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999
DKT 10 QUAD4
DKT 20 QUAD4
DKT 50 QUAD4
DKT 100 QUAD4
BEAM POU_D_T 50
INFLUENCE MESH
ON MODELING
DKT (QUAD4)
COMPARISON WITH
THE SOLUTION BEAM
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0
2
4
6
8
10
12
14x10-2
The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of results DKT QUAD4 when the mesh is refined.
8.4 Parameters
of execution
Version:
NEW 5.04.19
Machine:
SGI-Origin2000 R12000
Overall dimension memory: 16 megabytes
Time CPU To use: 27.55 seconds
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
10/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
9 Modeling
D
9.1
Characteristics of modeling
Modeling COQUE_3D (TRIA7)
OE = OA =
L/2 2
AB//EFF//axis Z
y
X
Z
With
E
F
B
O
Boundary conditions (reference mark gobal)
- side AB: U = v =w =
X
=
y
=
Z
=
0
Conditions of symmetry
- side EFF: U =
0 (local reference mark X
1
y
1
Z
1
)
- side EFF:
Z
=
0 (total reference mark)
X
L
y
L
Z
L
the beam being tilted with 45°, the value of displacement
imposed is then: DX
E
=6.609 10
- 3
Mr. sin (45°)
9.2
Characteristics of the mesh
A number of nodes: 83
A number of meshs and type: 20 TRIA7
9.3 Functionalities
tested
Controls Key word
factor
Key word
AFFE_MODELE AFFE
MODELING
= ' COQUE_3D'
DEFI_MATERIAU ECRO_LINE D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
THICK
A_CIS =0
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR
PILOTING
VMIS_ISOT_LINE
COQUE_NCOU
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”
10 Results of modeling D
10.1 Values
tested
DX E DX
E
()
Identification Moments Reference
Aster %
difference
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.
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
11/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
10.3 Influence refinement of the mesh
displacement (m)
E
T
With
_
P
I
L
OT
AGE
EDF
Electricity
from France
Mechanical department and Numerical Models
Eta according to displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999
10 COQUE_3D TRIA7
20 COQUE_3D TRIA7
50 COQUE_3D TRIA7
100 COQUE_3D TRIA7
BEAM POU_D_T 50
INFLUENCE MESH
ON MODELING
COQUE_3D TRIA7
COMPARISON WITH
THE SOLUTION BEAM
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0
2
4
6
8
10
12
14x10-2
The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of the results hulls 3D TRIA7 when the mesh is refined.
10.4 Parameters
of execution
Version:
NEW 5.04.19
Machine:
SGI-Origin2000 R12000
Overall dimension memory: 16 megabytes
Time CPU To use: 72.54 seconds
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
12/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
11 Modeling
E
11.1 Characteristics of modeling
Modeling COQUE_3D (QUAD9)
OE = OA =
L/2 2
AB//EFF//axis Z
y
X
Z
With
E
F
B
O
Boundary conditions (reference mark gobal)
- side AB: U = v =w =
X
=
y
=
Z
=
0
Conditions of symmetry
- side EFF: U =
0 (local reference mark X
1
y
1
Z
1
)
- side EFF:
Z
=
0 (total reference mark)
X
L
y
L
Z
L
the beam being tilted with 45°, the value of displacement
imposed is then: DX
E
=6.609 10
- 3
Mr. sin (45°)
11.2 Characteristics of the mesh
A number of nodes: 54
A number of meshs and type: 10 QUAD9
11.3 Functionalities
tested
Controls Key word
factor
Key word
AFFE_MODELE AFFE
MODELING
= ' COQUE_3D'
DEFI_MATERIAU ECRO_LINE D_SIGM_EPSI
SY
AFFE_CARA_ELEM HULL
THICK
A_CIS = 10
5
AFFE_CHAR_MECA FORCE_COQUE NEAR
STAT_NON_LINE COMP_INCR
PILOTING
VMIS_ISOT_LINE
COQUE_NCOU
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”
12 Results of modeling E
12.1 Values
tested
DX E DX
E
()
Identification Moments Reference
Aster %
difference
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.
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
13/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
12.3 Influence refinement of the mesh
displacement (m)
E
T
With
_
P
I
L
OT
AGE
EDF
Electricity
from France
Mechanical department and Numerical Models
Eta according to displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999
COQUE_3D 10 QUAD9
COQUE_3D 20 QUAD9
COQUE_3D 50 QUAD9
BEAM POU_D_T 50
COQUE_3D 100 QUAD9
INFLUENCE MESH
ON MODELING
COQUE_3D (QUAD 9)
COMPARISON WITH
SOLUTION BEAM
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0
2
4
6
8
10
12
14x10-2
The reference being taken compared to the solution beam with 50 elements, one observes one
improvement of the results hulls 3D QUAD9 when the mesh is refined.
12.4 Parameters
of execution
Version:
NEW 5.04.19
Machine:
SGI-Origin2000 R12000
Overall dimension memory: 16 megabytes
Time CPU To use: 74.93 seconds
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
14/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
13 Modeling
F
13.1 characteristics of modeling
78.5°
Modeling POU_D_T
A= 2.1 0.7
0.
AE =0.75 m
X
Z
E
With
O
Boundary conditions (reference mark gobal)
- node a: U = v =w = Q
X
= Q
y
= Q
Z
= 0
Conditions of symmetry
- node E: U = 0, Q
Z
= 0
y
the beam being tilted with 78.5°, the value of displacement
imposed is then: DX
E
=6.609 10
- 3
Mr. sin (78.5°)
13.2 characteristics of the mesh
A number of nodes: 21
A number of meshs and type: 10 SEG2
13.3 Functionalities
tested
Controls Key word
factor Key word
AFFE_MODELE AFFE
MODELING
= ' POU-D_T'
DEFI_MATERIAU ECRO_LINE
VMIS_POUTRE
D_SIGM_EPSI
SY
AFFE_CARA_ELEM BEAM
SECTION
CARA
AFFE_CHAR_MECA FORCE_POUTRE
FY
STAT_NON_LINE COMP_INCR
PILOTING
VMIS_POU_LINE
TYPE: “DDL_IMPO”
TEST_RESU
RESU
PARA: “ETA_PILOTAGE”
14 Results of modeling F
14.1 Values
tested
DX E DX
E
()
Identification Moments Reference
Aster %
difference
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.
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
15/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
14.2 Influence refinement of the mesh
displacement (m)
E
T
With
_
P
I
L
OT
AGE
EDF
Electricity
from France
Mechanical department and Numerical Models
Eta according to displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999
POU_D_T 10
POU_D_T 20
POU_D_T 100
POU_D_T 50
INFLUENCE MESH
ON MODELING
BEAM POU_D_T
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0
2
4
6
8
10x10-2
the solution beam improves appreciably when the mesh is refined.
14.3 Parameters
of execution
Version:
NEW 5.04.19
Machine:
SGI-Origin2000 R12000
Overall dimension memory: 16 megabytes
Time CPU To use: 74.93 seconds
Code_Aster
®
Version
5.0
Titrate:
SSNL501 - Fixed beam subjected to a uniform pressure
Date:
19/10/01
Author (S):
P. MASSIN, F. LEBOUVIER
Key
:
V6.02.501-A
Page:
16/16
Manual of Validation
V6.02 booklet: Nonlinear statics of the linear structures
HI-75/01/010/A
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 piloting.
displacement (m)
P
With
R
With
M
ETR
E ETA
EDF
Electricity
from France
Mechanical department and Numerical Models
Eta according to displacement
agraf 12/04/2001 (c) EDF/DER 1992-1999
Modeling A
Modeling B
Modeling C
Modeling D
Modeling E
Modeling F
COMPARISON OF THE SOLUTIONS
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0
2
4
6
8
10
12
14
16x10-2
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 mesh used, the results obtained are
satisfactory. The errors observed are for modelings:
·
COQUE_C_PLAN: 6.9% (
With),
·
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 mesh 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.