Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
1/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
Organization (S):
EDF-R & D/AMA
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
Document: V6.04.505
SSNV505 - Contact of 2 beams into large
displacements
Summary:
This test represents a calculation of contact without friction between two beams in great displacements subjected to
a specific displacement. Initially the beams are not in contact, it occurs for a level of
loading given and then the two beams slip one on the other.
The analyzed result is the normal reaction of contact according to displacement. Results obtained with
various modelings are compared between them.
·
Modeling a: the structure is modelized with elements of COQUE_3D associated a mesh
QUAD9, with the FORCED method.
·
Modeling b: the structure is modelized in 3D plane deformations using elements HEXA8,
with the FORCED method.
·
Modeling C: the structure is modelized in 2D plane deformations using elements QUAD4,
with the FORCED method.
·
Modeling D: the structure is modelized in 3D plane deformations using elements HEXA8,
with the method CONTINUES.
·
Modeling E: the structure is modelized in 2D plane deformations using elements QUAD4,
with the method CONTINUES.
·
Modeling F: the structure is modelized in 3D using elements of beam POU_D_E associated with
meshs SEG2, with the methods LAGRANGIAN and CONTINUOUS. This modeling is distinguished
others, in what it is carried out in small transformations, while all the others are
carried out in great transformations.
·
Modeling G: the structure is modelized in 3D using elements of beam POU_D_E associated with
meshs SEG2, with the LAGRANGIAN method with friction. This modeling validates
possibility of taking into account the radius of the section of beam like initial play.
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
2/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
1
Problem of reference
1.1 Geometry
1.2
Properties of material
E = 6.8948x10
9
AP
Young modulus
= 0.3333
Poisson's ratio
µ = 0.
Coefficient of friction
1.3
Boundary conditions and loadings
·
C.L. : Sections at embedded points A and B
·
Loading: vertical displacement of the point C: D = - 790 mm
1.4 Conditions
initial
Without object.
2
Reference solution
2.1
Method of calculation used for the reference solution
One does not have reference solution, one makes only tests of not-regression.
Uy
2032
2032
508
914.4
X
Y
Z
Y
63.5
10.16
Section
1320.8
With
B
C
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
3/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
3 Modeling
With
3.1
Characteristics of modeling
Boundary conditions:
This modeling is carried out in COQUE_3D. In fact, the finite elements used are QUAD9.
When one modelizes contact between this type of elements, one binds the nodes mediums to the nodes
nodes. If the user does not take guard there, that can easily cause the appearance of null pivots (with
cause redundant conditions of Dirichlet). Also one tried to specify the conditions well with
following limits:
·
sides LAA3 and LBB1: DX=DY=DZ=DRX=DRY=DRZ=0
·
Nodes nodes C and C1: DY=-790 mm
·
All the nodes nodes of the structure: DZ=0, DRX=0
Conditions of contact:
·
surface main: group meshs POUTRE1 minus the group of meshs whose edge is
embedded according to LBB1,
·
surface slave: group meshs POUTRE2 minus the group of meshs whose edge is
embedded according to LAA3.
3.2
Characteristics of the mesh
A number of nodes:
261
A number of meshs:
158
SEG3 107
QUAD8 51
Y
Z
X
LBB1
LAA3
C1
C
Poutre1
Poutre2
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
4/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
3.3 Functionalities
tested
Controls Key word
factor Key word
MODI_MAILLAGE MODI_MAILLE
ORIE_NORM_COQUE
QUAD8_9
VECT_NORM
AFFE_MODELE
AFFE
MODELING = “COQUE_3D”
DEFI_MATERIAU ELAS
AFFE_CARA_ELEM HULL
THICK
COEF_RIGI_DRZ=0.001
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
NOEUD_BOUCLE
METHOD = “FORCED”
FRICTION = “WITHOUT”
DIST_1 = 31.75
DIST_2 = 31.75
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
DEFORMATION = ' GREEN_GR'
4
Results of modeling A
4.1 Values
tested
One tests the reaction to the displacement imposed on the hull. To obtain it, one calculates the reactions to
embeddings in LBB1 and LAA3 (one cannot simply measure the reaction to the nodes C and C1
because it contains also the forces due to the connections of the nodes mediums to the nodes nodes).
Identification Moments Reference
Aster
Reaction 0.2
Not-regression 4.37897E+02
Reaction 0.4
Not-regression 1.65012E+03
Reaction 0.6
Not-regression 3.35705E+03
Reaction 0.8
Not-regression 3.01163E+03
Reaction 1.
Not-regression
3.14692E+03
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
5/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
5 Modeling
B
5.1
Characteristics of modeling
Boundary conditions:
·
Surfaces SA and SB: DX=DY=DZ= 0
·
Lines LCC1: DY=-790 mm
·
All the nodes nodes of the structure: DZ=0
Conditions of contact:
·
surface main: group meshs SMAI
·
surface slave: group meshs SESC
5.2
Characteristics of the mesh
A number of nodes:
662
A number of meshs:
1247
POI1 1
SEG2 351
QUAD4 655
HEXA8 240
Y
Z
X
SB
SA
Poutre1
Poutre2
LCC1
SMAI
SESC
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
6/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
5.3 Functionalities
tested
Controls Key word
factor Key word
MODI_MAILLAGE ORIE_PEAU_3D
AFFE_MODELE
AFFE
MODELING = “3D”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
NOEUD_BOUCLE
METHOD = “FORCED”
FRICTION = “WITHOUT”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
DEFORMATION = “GREEN”
6
Results of modeling B
6.1 Values
tested
One tests the reaction to the displacement imposed on the POUTRE2. To obtain it, the reactions are calculated
with embeddings in SA and SB.
Identification Moments Reference
Aster
Reaction 0.2
Not-regression
4.36995E+02
Reaction 0.4
Not-regression
1.66834E+03
Reaction 0.6
Not-regression
3.26454E+03
Reaction 0.8
Not-regression
2.73716E+03
Reaction 1.
Not-regression
3.09524E+03
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
7/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
7 Modeling
C
7.1
Characteristics of modeling
Modeling is made in 2D plane deformations to find blocking in DZ imposed on the model
3D.
Boundary conditions:
·
sides LA3A and LB3B: DX=DY= 0
·
Nodes C: DY=-790 mm
Conditions of contact:
·
surface main: group of mesh SMAI
·
surface slave: group of mesh SESC
7.2
Characteristics of the mesh
A number of meshs:
415
SEG2 175
QUAD4 240
7.3 Functionalities
tested
Controls Key word
factor Key word
MODI_MAILLAGE ORIE_PEAU_2D
AFFE_MODELE
AFFE
MODELING = “2D”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
NOEUD_BOUCLE
METHOD = “FORCED”
FRICTION = “WITHOUT”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
DEFORMATION = “GREEN”
LA3A
LB3B
C
Poutre1
Poutre2
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
8/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
8
Results of modeling C
8.1 Values
tested
One tests the reaction to the displacement imposed on the hull. To obtain it, one calculates the reactions to
embeddings in LB3B and LA3A
Identification Moments Reference
Aster
Reaction 0.2
Not-regression
4.35668E+01
Reaction 0.4
Not-regression
1.69598E+02
Reaction 0.6
Not-regression
3.23491E+02
Reaction 0.8
Not-regression
2.67162E+02
Reaction 1.
Not-regression
3.09623E+02
8.2 Notice
To obtain in 2D results comparable with the COQUE_3D results and 3D, it is necessary to multiply them
preceding reactions by the width of the beam, is 10.16 Misters.
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
9/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
9 Modeling
D
9.1
Characteristics of modeling
Boundary conditions:
·
Surfaces SA and SB: DX=DY=DZ= 0
·
Line LCC1: DY=-790 mm
·
All the nodes nodes of the structure: DZ=0
Conditions of contact:
·
surface main: group meshs SMAI
·
surface slave: group meshs SESC
9.2
Characteristics of the mesh
A number of nodes:
662
A number of meshs:
1247
POI1 1
SEG2 351
QUAD4 655
HEXA8 240
Y
Z
X
SB
SA
Poutre1
Poutre2
LCC1
SMAI
SESC
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
10/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
9.3 Functionalities
tested
Controls Key word
factor Key word
MODI_MAILLAGE ORIE_PEAU_3D
AFFE_MODELE
AFFE
MODELING = “3D”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
METHOD = “CONTINUES”
FRICTION = “WITHOUT”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
DEFORMATION = “GREEN”
10 Results of modeling D
10.1 Values
tested
One tests the reaction to the displacement imposed on the POUTRE2. To obtain it, the reactions are calculated
with embeddings in SA and SB. The percentage of difference indicates the difference between this modeling
and equivalent modeling with the method
STRESS
.
Identification Moments Reference
Aster
% difference
Reaction 0.2
Not-regression 4.36995E+02
4.25E-05
Reaction 0.4
Not-regression 1.668034E+03 - 0.018
Reaction 0.6
Not-regression 3.26453E+03
- 3.13E-04
Reaction 0.8
Not-regression 2.74197E+03 0.176
Reaction 1.
Not-regression 3.09252E+03 - 0.088
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
11/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
11 Modeling
E
11.1 Characteristics of modeling
Modeling is made in 2D plane deformations to find blocking in DZ imposed on the model
3D.
Boundary conditions:
·
sides LA3A and LB3B: DX=DY= 0
·
Nodes C: DY=-790 mm
Conditions of contact:
·
surface main: group of mesh SMAI
·
surface slave: group of mesh SESC
11.2 Characteristics of the mesh
A number of meshs:
415
SEG2 175
QUAD4 240
11.3 Functionalities
tested
Controls Key word
factor Key word
MODI_MAILLAGE ORIE_PEAU_2D
AFFE_MODELE
AFFE
MODELING = “2D”
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
PAIRING = “MAIT_ESCL”
METHOD = “CONTINUES”
FRICTION = “WITHOUT”
STAT_NON_LINE
COMP_ELAS
RELATION = “ELAS”
DEFORMATION = “GREEN”
LA3A
LB3B
C
Poutre1
Poutre2
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
12/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
12 Results of modeling E
12.1 Values
tested
One tests the reaction to the displacement imposed on the hull. To obtain it, one calculates the reactions to
embeddings in LB3B and LA3A. the percentage of difference indicates the difference between this modeling
and equivalent modeling with the method
STRESS
.
Identification Moments Reference
Aster
% difference
Reaction 0.2
Not-regression 4.35668E+01 6.9E-07
Reaction 0.4
Not-regression 1.69562E+02 - 0.021
Reaction 0.6
Not-regression 3.23491E+02 9.1E-05
Reaction 0.8
Not-regression 2.67527E+02 0.14
Reaction 1.
Not-regression 3.09623E+02 1.7E-04
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
13/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
13 Modeling
F
13.1 Characteristics of modeling
One carries out here a modeling using elements beam in 3D. Deformations being plane in
the plan (DX, DY), one imposes DZ = 0 on the model 3D. The goal of this case-test is to compare
deformation of the beams with a formulation in small rotations with that obtained with the large ones
rotations. The first model is of course abusive (forgery) compared to the second (true), but allows
to illustrate the difference of the results obtained in one or the other case. The true motivation of this case
test is however of exhiber an example of validation of the contact between beams with the method
continue for the preliminary stage of the small transformations. The test will be enriched soon
of a modeling in great rotations, more delicate to implement.
Boundary conditions:
·
Nodes A and b: DX=DY= 0
·
Nodes C: DY=-790 mm
Conditions of contact:
·
surface main: group of mesh POU1
·
surface slave: group of mesh POU2
13.2 Characteristics of the mesh
A number of meshs:
4
SEG2 4
13.3 Functionalities
tested
Controls Key word
factor Key word
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
METHOD = “CONTINUES”
VECT_Y = (0,0, - 1)
DIST_MAIT=31.75
DIST_ESCL=31.75
STAT_NON_LINE
COMP_INCR
RELATION = “ELAS”
A1
With
B1
B
C
POU2
POU1
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
14/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
14 Results of modeling F
14.1 Values
tested
One tests the reaction to the displacement imposed on the beam. To obtain it, one calculates the reactions to
embeddings in A and B. the percentage of difference indicates the difference between this modeling and
equivalent modeling with the method
LAGRANGIAN
.
On the curve of bottom, one represented the force at the point B (embedding of the higher beam) in
function of the evolution of the loading. One compares the methods Lagrangian and continues with
solution obtained in great rotations. It appears that the two methods give results
almost identical, but that the latter differ notably from those obtained into large
rotations. This point is rather logical, since within the framework of small rotations where calculations were
carried out, one neglects the terms of deformation of the second command, which as one notes it are not
not negligible in great rotations.
Identification Moments Reference
Aster
% difference
Reaction 0.2
Not-regression 3.07483D+02
8.47E-11
Reaction 0.4
Not-regression 1.25148D+03
- 5.07E-09
Reaction 0.6
Not-regression 3.02391D+03
- 2.23E-07
Reaction 0.8
Not-regression 4.82957D+03
- 7.10E-07
Reaction 1.
Not-regression 6.65614D+03
- 1.12E-06
SSNV505F: NODAL REACTIONS
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
15/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
15 Modeling
G
15.1 Characteristics of modeling
One carries out here a modeling using elements beam in 3D. Deformations being plane in
the plan (DX, DY), one imposes DZ = 0 on the model 3D. The goal of this case-test is to validate the catch in
count real section of the beam, that the user informed by the keyword “BEAM” in
AFFE_CARA_ELEM.
Boundary conditions:
·
Nodes A and b: DX=DY= 0
·
Nodes C: DY=-790 mm
Conditions of contact:
·
surface main: group of mesh POU1
·
surface slave: group of mesh POU2
Characteristics of the beam:
· tubular section of radius 31.75 mm and thickness 1mm
15.2 Characteristics of the mesh
A number of meshs:
4
SEG2 4
15.3 Functionalities
tested
Controls Key word
factor Key word
DEFI_MATERIAU ELAS
AFFE_CHAR_MECA
CONTACT
METHOD = “LAGRANGIAN”
VECT_Y = (0,0, - 1)
DIST_MAIT=31.75
DIST_POUTRE = “YES”
STAT_NON_LINE
COMP_INCR
RELATION = “ELAS”
A1
With
B1
B
C
POU2
POU1
Code_Aster
®
Version
8.3
Titrate:
SSNV505 - Contact of two beams in great displacements
Date:
04/05/06
Author (S):
Mr. KHAM, Mr. ABBAS
Key
:
V6.04.505-B
Page:
16/16
Manual of Validation
V6.04 booklet: Nonlinear statics of the voluminal structures
HT-62/06/005/A
16 Results of modeling G
16.1 Values
tested
One tests the reaction to the displacement imposed on the beam. One compares the case where the section of
beam entered via a constant play DIST_ESCL and the case where one takes the real section via
keyword DIST_POUTRE.
Identification Moments
Aster (DIST_ESCL)
Aster
(DIST_POUT)
% difference
Reaction 0.16
2.45986E+02
2.45986E+02
0
Reaction 0.4
1.25149D+03
1.25149D+03
0
17 Summary of the results
The graph below presents the evolution of the component DY of the force of reaction to
displacement imposed according to this last.
One notices a very good agreement between various modelings up to 500 mm then the curve
COQUE_3D separates from the 2D and of the 3D before meeting with 700 Misters This variation is normal: it
appears when the end of beam 2 is orthogonal with beam 1.