Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
1/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Organization (S):
EDF-R & D/AMA
Manual of Reference
R3.08 booklet: Machine elements with average fiber
R3.08.06 document
Finite elements of right pipe and curve
with ovalization, swelling and roll
in elastoplasticity
Summary:
This document presents the modeling of a finite element of pipe usable in calculations of pipings in
elasticity or in plasticity. The pipes, curves or rights, can be relatively thick (thickness report/ratio on
radius of the transverse section up to 0.2) and are subjected to various combined loadings - internal pressure,
cross-bendings and anti-plane, torsion, extension - and can have a nonlinear behavior.
This linear element combines at the same time properties of hulls and beams. The average fiber of the pipe
comprise like a beam and the surface of the pipe like a hull. The element carried out is an element of
right pipe or curve in small rotations and deformations, with an elastoplastic behavior in
plane stresses.
Three modelings, corresponding to three various types of elements, are available:
·
TUYAU_3M
, which takes into account 3 modes of Fourier to the maximum, and which can rest on
meshs with 3 nodes or 4 nodes.
·
TUYAU_6M
, which takes as a count up to 6 modes of Fourier, and is pressed on meshs with 3 nodes.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
2/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Count
matters
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
3/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
4/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
1 Introduction
There is an important bibliography on the modeling of pipings and many elements
stop pipes right and bent are available in the great codes of finite elements.
syntheses were already realized [bib1], [bib5], [bib6], in the past that one supplemented in
incorporating the last developments known in the field [bib11]. Important effects with
to take into account are swelling due to the internal pressure and the ovalization of the transverse sections
by combined bendings plane and anti-plane. One places oneself on the assumption of small rotations and
deformations within the framework of this document.
It is about a linear element with 3 or 4 nodes, of curved or right beam type with local plasticity
taking into account ovalization, the roll and swelling. The kinematics of beam is
enriched by a kinematics of hull for the description of the behavior of the transverse sections.
This kinematics is discretized in M modes of Fourier of which the number M must at the same time be sufficient
to obtain good results in plasticity and not too large to limit the calculating time.
literature encourages us to use M=6 [bib9], [bib13] in plasticity. In elasticity, for thick pipes, one
can be satisfied with M=2 or 3.
2
Various theories of hulls and beams for
finite elements of right or bent pipes
One presents in this chapter the elements of kinematics in three-dimensional curvilinear geometry,
as their restrictions within the framework of the models of beam and hull. Indeed, to build
the finite element of piping enriched which answers the schedule of conditions defined in introduction, one exploits
a technique of decomposition of three-dimensional kinematics. The kinematics of hull y
bring the description of ovalization, swelling and roll, while kinematics
of beam described there the general movement of the line of piping.
The various theories of hulls and beams used for each element translate them
assumptions chosen a priori on the type of deformations and behaviors.
2.1
The pipe in theory of beam
2.1.1 Case of a bent pipe
A first approaches relatively simple come down to consider the elbow represented below
like a beam digs circular section. The beam is obtained by rotation of angle
circular section around OZ. A point of the beam is identified by its distance R compared to the axis
beam and by the two angles
,
where
is the longitudinal angle with OY indicated above and
the trigonometrical angle with OZ measured on the circular section.
·
·
R
Z
y
X
Z
y
X
1
2
U
Z
U
y
U
X
·
Z
O
Z
Y
X
·
·
R
Z
y
X
1
2
·
E
R
E
E
Z
O
Z
Y
X
()
R
Appear 2.1.1-a: Geometry and kinematics of the elbow in theory of beam
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
5/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
In the curvilinear frame of reference (
R
,
,
), relations between displacements
U
points of the elbow of position
)
,
(
)
(
R
y
0
E
E
OM
R
R
R
+
-
=
=
and deformations of
Green-Lagrange are given by the following tensor in the natural base (
R
,
,
):
(
) (
)
0
0
0
0
R
R
U
R
U
R
.
.
2
-
+
+
=
F
,
(
) {
}
,
,
,
R
.
The vectors units in the directions (
R
,
,
) are:
0
0
0
R
R
E
R
E
R
E
B
With
R
1
,
1
,
=
=
=
where
0
0
0
0
R
R
R
R
.
.
=
=
B
With
and
.
If one expresses the position of a point of the elbow in the local toric orthonormée base (
E
E
E
R
,
,
)
by (
y
y
y
R
,
,
) there are the following relations:
y
y
y
R
0
0
0
R
R
E
R
E
R
E
=
=
=
,
,
.
The form of the tensor of the deformations of Green-Lagrange in this base is then:
(
) (
)
y
y
y
y
0
0
0
0
R
R
U
R
U
R
.
.
2
-
+
+
=
.
Relations of passage between the expression of the deformations of Green-Lagrange in the system of
co-ordinates curvilinear and in the local toric base previously definite are:
B
F
With
F
AB
F
B
F
With
F
F
R
R
R
R
rr
rr
=
=
=
=
=
=
2
2
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
6/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
The use of this base is particularly interesting because the relations of behavior in
base toric orthonormée are simple of use. For elbow Ci above, if it is considered that them
deformations remain small, one obtains then [bib4] after linearization of the deformations of
Green-Lagrange:
R
U
R
U
U
B
R
U
R
U
U
With
With
AB
U
U
With
U
B
R
U
U
B
R
U
With
AB
U
U
With
R
U
R
R
R
R
R
R
R
R
R
rr
+
-
=
=
+
-
=
=
-
+
=
=
+
=
+
+
=
=
1
2
1
2
1
1
2
1
1
with:
.
,
sin
sin
,
,
sin
R
R
R
R
R
R
B
R
R
With
=
+
=
=
+
=
The expressions of the deformations established above are written then:
R
U
U
U
R
R
U
U
U
R
R
U
R
U
U
R
R
U
U
R
U
U
U
R
R
R
U
R
R
R
R
R
R
R
R
R
rr
+
-
=
=
+
-
+
=
=
+
-
+
=
=
+
=
+
+
+
=
=
)
(
1
2
)
sin
(
sin
1
2
1
)
cos
(
sin
1
2
)
(
1
)
sin
cos
(
sin
1
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
7/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Displacement
R
U
,
U
U,
of a point of the elbow in the toric base associated the transverse section
of observation can easily express itself according to displacements and rotations associated with
center transverse section. Indeed, if one notes
3
2
1
,
,
U
U
U
displacement in the curvilinear base
local
() () () ()
(
)
Z
y
X
,
,
,
O
associated the transverse section as indicated on [Figure 2.1.1-a]
there are the following, valid relations within the framework of the kinematics of the beams of Timoshenko
[R3.08.01]:
sin
)
(
)
(
)
,
,
(
cos
)
(
)
(
)
,
,
(
cos
)
(
sin
)
(
)
(
)
,
,
(
3
2
1
R
U
R
U
R
U
R
U
R
R
U
R
U
X
Z
X
y
y
Z
X
-
=
+
=
-
+
=
where
Z
y
X
U
U
U
,
,
is the displacement of translation of the section and
X
,
y
,
Z
the rotation of its center
O
. The expression of the components of displacement in the local toric orthonormée base
(
E
E
E
R
,
,
) is obtained by change of reference mark:
]
sin
)
(
cos
)
(
[
]
sin
)
,
,
(
cos
)
,
,
(
[
)
,
,
(
)
(
cos
)
(
sin
)
(
cos
)
,
,
(
sin
)
,
,
(
)
,
,
(
cos
)
(
sin
)
(
)
(
)
,
,
(
)
,
,
(
2
3
2
3
1
y
Z
R
X
y
Z
y
Z
X
U
U
R
U
R
U
R
U
R
U
U
R
U
R
U
R
U
R
R
U
R
U
R
U
+
-
=
+
-
=
-
-
=
-
=
-
+
=
=
The introduction of this field of displacement into the expression of the linearized deformations us
allows to obtain the expression of the three-dimensional deformations associated the kinematics of
beam:
0
2
)
cos
sin
(
)
sin
cos
sin
cos
sin
sin
(
sin
1
2
)
sin
cos
(
)
cos
sin
cos
sin
cos
cos
(
sin
1
2
0
)
cos
sin
cos
(
sin
1
0
2
,
,
2
,
,
,
,
,
,
=
-
+
-
+
-
-
-
+
=
+
+
-
+
-
+
-
-
+
=
=
-
+
-
-
+
=
=
R
y
Z
Z
y
Z
y
X
R
y
Z
Z
y
X
Z
y
X
y
Z
X
y
X
rr
R
R
U
U
U
R
R
R
R
R
U
U
U
R
R
R
R
R
U
U
R
R
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
8/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
2.1.2 Case of the right pipe
The expressions of the deformations established above also apply to the case of the right pipe, where one
replace
by
S
where
S
is the curvilinear X-coordinate along average fiber of the pipe, with:
.
,
0
/
1
,
,
1
R
R
R
R
B
With
=
=
=
=
The expressions given for the elbow are written then for the right pipe:
R
U
U
U
R
R
U
X
U
U
R
X
U
U
U
R
X
U
R
U
R
R
R
X
R
X-ray
X-ray
X
X
X
R
X
xx
R
rr
+
-
=
=
+
=
=
+
=
=
+
=
=
=
)
(
1
2
2
1
2
)
(
1
Z
y
X
Z
y
X
1
2
U
y
·
Z
O
Z
Y
X
·
X
U
Z
U
Z
y
Z
X
U
R
U
U
2
U
3
With
Transverse section: sight of 1 towards 2
U
X
=u
1
O
Appear 2.1.2-a: Geometry and kinematics of a right pipe in theory of beam
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
9/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
As previously, displacement
R
U
,
U
U
X
,
of a point of the pipe in the associated toric base
with the transverse section of observation can easily express itself according to displacements and
rotations associated with the center with the transverse section. Indeed, if one notes
3
2
1
,
,
U
U
U
displacement
in the local curvilinear base
(
)
Z
y
X,
,
,
O
associated the transverse section as indicated on the figure
below there are the following relations:
sin
)
(
)
(
)
,
,
(
cos
)
(
)
(
)
,
,
(
cos
)
(
sin
)
(
)
(
)
,
,
(
3
2
1
R
X
X
U
X
R
U
R
X
X
U
X
R
U
R
X
R
X
X
U
X
R
U
X
Z
X
y
y
Z
X
-
=
+
=
-
+
=
and:
]
sin
)
(
cos
)
(
[
]
sin
)
,
,
(
cos
)
,
,
(
[
)
,
,
(
)
(
cos
)
(
sin
)
(
cos
)
,
,
(
sin
)
,
,
(
)
,
,
(
cos
)
(
sin
)
(
)
(
)
,
,
(
)
,
,
(
2
3
2
3
1
X
U
X
U
X
R
U
X
R
U
X
R
U
R
X
X
U
X
U
X
R
U
X
R
U
X
R
U
R
X
R
X
X
U
X
R
U
X
R
U
y
Z
R
X
y
Z
y
Z
X
X
+
-
=
+
-
=
-
-
=
-
=
-
+
=
=
The introduction of this field of displacement into the expression of the deformations given below
us allows to obtain the expression of the deformations associated with kinematics with beam:
0
2
cos
)
(
sin
)
(
2
cos
)
(
sin
)
(
2
0
cos
sin
0
,
,
,
,
,
,
,
,
=
+
-
-
=
-
+
+
+
-
=
=
-
+
=
=
R
X
Z
y
X
y
Z
X-ray
X
y
Z
X
Z
y
X
X
X
X
y
X
Z
X
X
xx
rr
U
U
U
U
R
R
R
U
2.1.3 Remarks
The fact that
R
rr
and
,
are simultaneously null shows that the kinematics of beam cannot
to represent the deformations of the transverse sections to average fiber of the pipe. Indeed, them
transverse sections are actuated by a rigid movement of body, which prohibits to modelize it
roll, swelling and ovalization.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
10/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
2.2
The pipe in linearized theory of hull
2.2.1 Case
General
The bent pipe is regarded as a thin hull of revolution (portion of core). Surface
average is obtained by rotation of angle
of a circle of radius has whose center is at a distance
R of the axis of revolution OZ. One indicates by H the thickness of the elbow. One imposes on this thickness
to remain constant like with the section of the elbow being perfectly circular. A point on surface
average is characterized by the two angles
,
and its position
2
/
2
/
H
H
+
-
compared to
surface average, where
is the longitudinal, variable angle between 0 and
, and
the angle measured on
transverse section.
R: Radius of curvature
R: Radius of the cross section
H: Thickness of the elbow
: Longitudinal angle
: Angle of cross section
U: Axial displacement of average surface
v: Orthoradial displacement of average surface
W: Radial displacement of average surface
R
R
H
U
v
W
O
Z
: Rotation of average surface compared to E
: Rotation of average surface compared to E
y
X
O
Z
E
R
E
E
E
y
X
Appear 2.2.1-a: Geometry and kinematics of the elbow in theory of hull
One places oneself first of all within the framework of the linearized theory of the hulls with shearing
transverse such as it was described for example in Washizu [bib14]. This choice had already was made
for the linear elements of hulls [R3.07.01]. It limits our study to the framework of small
deformations. Moreover, great rotations of average surface are not taken into account.
Displacements and rotations are thus defined compared to the initial geometry of the elbow. If them
displacements of the points of average surface in the three directions
axial,
orthoradiale and
radial are noted
v
U,
and
W
those of any point of the elbow
are written in the following way:
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
W
U
v
U
U
U
=
-
=
+
=
where
and
are rotations compared to the vectors
E
and
E
respectively.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
11/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
The deformations in any point are thus given by [bib14]:
(
)
(
)
(
)
(
)
R
E
R
E
R
R
E
R
E
R
E
/
1
2
2
/
1
2
2
/
1
/
1
2
2
2
/
1
/
1
+
=
=
+
=
=
+
+
+
=
=
+
+
=
+
+
=
with:
.
,
sin
sin
,
,
sin
has
R
has
R
R
has
B
has
R
With
=
+
=
=
+
=
where
E
E,
and
E
are the membrane deformations of average surface,
,
,
deformations of bending of average surface and
E
E,
transverse distortions.
deformations of average surface are connected to displacements of average surface in
replacing the field of displacement of the preceding paragraph by that given above. One finds
then:
R
v
W
B
E
R
U
W
With
E
With
AB
U
v
With
R
U
B
R
With
With
AB
B
B
With
AB
With
With
AB
U
v
With
U
B
E
R
W
v
B
E
R
W
With
AB
v
U
With
E
-
+
-
=
-
+
=
-
+
+
-
-
=
-
=
-
=
-
+
=
+
=
+
+
=
1
2
1
2
)]
1
(
1
1
1
[
1
1
2
1
1
1
1
2
1
1
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
12/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
That is to say still:
)
(
1
)
sin
(
sin
1
)]
cos
(
1
sin
1
1
sin
sin
[
)
cos
(
sin
1
1
1
)
cos
(
sin
1
1
)
cos
(
sin
1
)
(
1
)
sin
cos
(
sin
1
v
W
has
U
W
has
R
U
v
has
has
R
U
has
has
R
has
R
has
has
has
R
U
has
U
v
has
R
W
v
has
E
W
v
U
has
R
E
-
+
-
=
-
+
+
=
-
+
+
+
+
+
+
-
=
-
=
-
+
=
+
-
+
=
+
=
+
+
+
=
In this theory there are thus five unknown factors; 3 displacements
v
U,
and
W
like two
rotations
,
. If the assumption of Coils-Kirchhoff is applied (thin tube) shearings
transverses are null and there are nothing any more but 3 displacements
v
U,
and
W
since:
)
(
1
)
sin
(
sin
1
v
W
has
U
W
has
R
-
=
-
+
-
=
2.2.2 Case of the right pipe
If one applies these equations to the case of the right pipe with:
.
,
0
/
1
,
,
1
has
R
R
has
B
With
=
=
=
=
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
13/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
One finds the more usual expression for this kind of geometry:
)
(
1
2
2
]
1
[
1
2
1
1
2
)
(
1
v
W
has
E
X
W
E
X
v
has
X
has
has
X
U
has
X
v
E
W
v
has
E
X
U
E
X
X
X
X
X
xx
X
xx
-
+
-
=
+
=
+
-
=
-
=
=
+
=
+
=
=
In this theory there are thus five unknown factors; 3 displacements
v
U,
and
W
like two
rotations
,
X
. If the assumption of Coils-Kirchhoff is applied (thin tube) shearings
transverses are null and there are nothing any more but 3 displacements
v
U,
and
W
since:
)
(
1
v
W
has
X
W
X
-
=
-
=
2.2.3 Notice
One can directly introduce the kinematics of hull into the field of deformation 3D. In it
case one a:
E
E
E
E
E
2
2
2
2
2
2
2
=
=
=
=
+
=
=
+
=
+
=
where expressions
E
E,
and
E
for the membrane deformations,
,
,
for
deformations of bending and
E
E,
for the transverse distortions are given by the expression
following in the general case:
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
14/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
)
(
1
2
)
sin
(
sin
1
sin
sin
2
)
cos
(
sin
1
1
2
1
)
cos
(
sin
1
)
cos
(
sin
1
1
2
)
(
1
)
sin
cos
(
sin
1
v
W
R
R
has
E
U
W
R
R
R
R
has
R
E
R
R
R
R
R
R
U
v
R
R
U
R
E
W
v
R
E
W
v
U
R
R
E
-
+
-
=
-
+
+
+
+
=
+
+
-
=
-
=
-
+
=
-
+
+
=
+
=
+
+
+
=
One notices that with command 1 in
the two ways of proceeding give identical results. It is
the definition of the deformation of membrane or bending which changes. In the first case it is
independent of the position in the thickness and is calculated for the average radius of the section
transverse of the pipe, whereas it depends on it in the case on the approach 3D. The term between hook
in the expression of [§2.2.1]. represent a coupling between the bending and the membrane which appears
when one expresses
sin
R
R
+
and
R
according to
sin
has
R
+
and
has
. In the continuation of our analysis
we will use this expression 3D degenerated of the kinematics of hull.
If moreover we use the assumption of Coils-Kirchhoff for transverse shearings,
0
=
=
E
E
one finds well the following expressions of rotations:
)
(
1
)
sin
(
sin
1
v
W
has
U
W
has
R
-
=
-
+
-
=
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
15/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
and:
)
sin
(
1
)
cos
sin
(
)
sin
(
1
]
)
sin
(
1
)
sin
(
1
[
]
)
sin
(
cos
)
sin
) (
sin
(
cos
) [
sin
(
2
)
(
1
)]
(
cos
)
sin
(
sin
1
[
sin
1
2
2
2
2
2
2
has
R
R
U
U
R
R
has
v
has
R
R
R
R
has
W
has
R
R
has
has
R
R
R
U
W
v
W
rear
v
W
has
U
W
has
R
R
R
+
+
+
+
+
+
+
+
-
+
+
+
+
-
=
-
-
=
-
-
-
+
-
+
=
expressions for
E
E
E
and
,
remainder unchanged.
One can easily extend this remark to the case of the right pipe.
2.3
Analyze right and bent pipes
In conclusion of the two preceding analyzes one can modelize the pipe like an element of
beam whose section is a thin hull. This interpretation is made in the majority of the codes
([bib2], [bib8], [bib9], [bib10], [bib12], etc…). In the absence of roll of the transverse sections
(i.e the transverse sections remain plane) the axial displacement of beam gives the new position
transverse section and displacements of ovalization (it is enough to take then u=0 in
mean equations of hulls) make it possible to know how the aforementioned becomes deformed. Total deflection
is obtained like superposition of the deformations of beam and the deformations of ovalization.
field of displacement which one represents on the figure below writes:
S
p
U
U
U
+
=
.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
16/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
In the first field of displacement the image of the transverse section is a transverse section
identical obtained by translation and rotation of the first. In the second field of displacement,
transverse section is deformed.
M
M
bending-torsion of a right beam
In theory of the beams-Euler
In theory of the hulls
Cross
Transverse section
Transverse section
ovalization
Cross
roll
U
v
W
Appear 2.3-a: Decomposition of displacement in fields of beam and hull
Modeling finite element must thus give an account of two different mechanical answers: that
beam and that of the hull for ovalization, swelling and the roll. These three
last modelings utilize degrees of freedom which are not nodal (decomposition
in Fourier series for example).
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
17/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
3
Mixed elements hull-beam for the right pipes and
curves
3.1 Kinematics
One breaks up the field of displacement into a macroscopic part of “beam” and a part
additional local of “hull”.
V
is the useful space of the fields of displacements
three-dimensional definite on an unspecified section of pipe.
For the beam part, as in [R3.03.03], one introduces space
T
fields associated with a torque
(defined by two vectors):
(
)
()
{
}
GM
v
v
T
+
=
=
T
T
M
V
that
such
,
/
For the fields of displacement of
T
,
T
is the translation of the section (or point
G
),
infinitesimal rotation and fields
v
are displacements preserving the section
S
plane and not
deformation there (One uses still the assumptions of NAVIER-BERNOULLI).
T
is a vectorial subspace of finished size equalizes to 6. It has additional orthogonal
for the scalar product on
V
:
=
=
T
W
W
v
v
T
0
.
/
S
V
.
Any field
U
of
V
breaks up then in a single way all in all of an element of
T
and of one
element of
T
:
+
=
T
U
T
U
U
U
U
S
p
S
p
,
.
One postulates then for displacements of surface of the pipe defined in [§2.2] the decomposition in
following Fourier series who check the preceding principle of orthogonality with displacements of
beam until command 3 in the thickness of the pipe:
=
=
M
m
im
m
X
U
X
U
2
cos
)
(
)
,
(
=
+
M
m
O
m
m
X
U
2
sin
)
(
=
+
=
M
m
inm
in
m
X
v
X
W
X
v
2
1
sin
)
(
sin
)
(
)
,
(
............
=
+
-
M
m
O
Nm
one
m
X
v
X
W
2
1
cos
)
(
cos
)
(
)
(
)
,
(
X
W
X
W
O
N
=
(uniform radial expansion)
.........
=
+
M
m
inm
m
X
W
1
cos
)
(
.........
=
+
M
m
O
Nm
m
X
W
1
sin
)
(
where
X
is the curvilinear X-coordinate along the elbow or right pipe, indifferently, and
M
the number of
modes of Fourier. Rotations
)
,
(
X
X
and
)
,
(
X
result from
() ()
()
,
,
,
,
X
W
X
v
X
U
and
by
the relations of Coils-Kirchhoff [§2.2.1].
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
18/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Note:
One can note that in the decomposition of
()
()
,
,
X
W
X
v
and
terms in
sin
cos
and
are not completely independent because of orthogonality with displacements of
beam. This makes it possible moreover to avoid the movements of rigid body, because if
O
N
in
O
N
in
W
W
v
v
1
1
1
1
,
,
and
are independent, one can find a solution nonnull giving
null deformations. In addition in the expression of
()
,
X
U
the absence of the terms is noted
in
sin
cos
and
already present in the beam part.
If one neglects the variation of metric with the thickness of the pipe the conditions of orthogonality
rigorous between displacements of beam and those of the surface of the pipe are satisfied.
In the contrary case, to satisfy this condition rigorously one would need one
development in Fourier series of rotations
)
,
(
X
X
and
)
,
(
X
starting with
command 2. This is incompatible with the assumptions of Love_Kirchhoff for these rotations.
3.2
Law of behavior
The behavior of the new element is a behavior 3D in plane stresses, because it
total behavior of the structure is that of a thin hull. It results from it that
0
=
and the law of
behavior is written in a general way in the following way:
=
+
+
+
+
+
=
X
X
xx
sx
p
xr
S
p
R
sx
p
X
S
p
sxx
p
xx
X
X
xx
C
C
In our case one will neglect transverse shearings for the hull part of our field of
displacement. It thus results from it that
0
=
=
S
sx
. As in addition [§2.1.2] it was shown that
0
=
p
R
it results from it that
0
=
. For an elastic behavior one has as follows:
-
-
-
=
X
X
xx
X
X
xx
v
v
v
E
2
1
0
0
0
0
2
1
0
0
0
0
1
0
0
1
1
2
and
-
-
-
=
2
1
0
0
0
0
2
1
0
0
0
0
1
0
0
1
1
2
v
v
v
E
C
.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
19/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
3.3
Work of deformation
The general expression of the work of deformation 3D for the element of elbow with the type of
above mentioned behavior is worth:
FD
W
L
X
X
X
X
H
H
xx
xx
def
+
+
+
=
-
0
2
0
2
/
2
/
)
(
where
L
is the curvilinear X-coordinate which is worth
)
sin
(
R
R
L
+
=
for an elbow where
is the traversed angle
to describe the elbow. In the case of an elbow, one has thus
D
rd
D
R
R
FD
)
sin
(
+
=
and for one
right pipe
D
dxrd
FD
=
where
. is the position in the thickness of the elbow which varies between - H/2 and
+h/2. In the continuation, in order to reduce the notations, one will employ the second expression.
3.4
Energy interns elastic elbow
In the case of an elastic behavior, energy interns elastic elbow is expressed way
following:
FD
G
E
X
X
xx
xx
L
H
H
))
(
)
2
(
1
(
2
1
2
2
2
2
0
2
0
2
2
/
2
/
+
+
+
+
-
=
-
int
This energy can be broken up into part of energy of beam, part of energy for
surface pipe and terms of coupling of the type
-
2
/
2
/0
2
0
.
H
H
L
sxx
p
xx
FD
.
3.5
Work of the forces and couples external
With the decomposition of displacements stated at the head of paragraph, the work of the forces
being exerted on the pipe expresses itself in the following way:
S
p
S
C
S
S
S
v
p
C
p
S
p
v
S
P
C
S
P
S
S
P
v
U
F
U
F
U
F
U
F
U
F
U
F
U
U
F
U
U
F
U
U
F
ext.
ext.
H
H
L
L H
H
H
H
L
L H
H
H
H
L
L H
H
ext.
W
W
D
rd
dx
D
H
has
FD
D
rd
dx
D
H
has
FD
D
rd
dx
D
H
has
FD
W
+
=
+
±
+
+
+
±
+
=
+
+
±
+
+
+
=
+
-
+
-
+
-
+
-
+
-
+
-
.
)
2
/
(
.
.
.
)
2
/
(
.
.
)
.(
)
2
/
(
)
.(
)
.(
2
/
2
/
2
0
0
2
0
0
2
/
2
/
2
0
2
/
2
/
2
0
0
2
0
0
2
/
2
/
2
0
2
/
2
/
2
0
0
2
0
0
2
/
2
/
2
0
by simple linear decomposition, where
C
S
v
F
F
F
,
,
are the voluminal, surface efforts and of contour
being exerted on the pipe, respectively. One determines as follows:
D
rd
dx
D
H
has
FD
W
H
H
L
L H
H
ext.
p
C
p
S
p
v
p
U
F
U
F
U
F
.
)
2
/
(
.
.
2
/
2
/
2
0
0
2
0
0
2
/
2
/
2
0
+
-
+
-
+
±
+
=
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
20/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
and:
D
rd
dx
D
H
has
FD
W
H
H
L
L H
H
ext.
S
C
S
S
S
v
S
U
F
U
F
U
F
.
)
2
/
(
.
.
2
/
2
/
2
0
0
2
0
0
2
/
2
/
2
0
+
-
+
-
+
±
+
=
The work of the external forces can thus be separate in two contributions distinct from same
forces, on the kinematics of beam and its additional.
3.6
Principle of virtual work
It is written in the following way:
int
ext.
ext.
ext.
W
W
W
W
=
+
=
S
p
with:
FD
W
X
X
X
X
xx
xx
L
H
H
int
)
(
0
2
0
2
/
2
/
+
+
+
=
-
3.6.1 External efforts part and couples for the beam part
The discretization of the principle of virtual work for the external efforts gives:
L
Z
Z
y
y
X
X
Z
Z
y
y
X
X
L
Z
Z
y
y
X
X
Z
Z
y
y
X
X
ext.
U
U
U
dx
m
m
m
U
F
U
F
U
F
W
,
0
0
]
[
)
(
µ
µ
µ
+
+
+
+
+
+
+
+
+
+
+
=
p
Z
y
X
F
F
F
,
,
: linear forces acting according to
X
,
y
and
Z
passing by the center of gravity of
transverse sections:
D
H
has
D
rd
F
H
H
I
)
2
/
(
.
.
2
0
2
/
2
/
2
0
±
+
=
+
-
I
S
I
v
E
F
E
F
where
Z
y
X
E
E
E
,
,
are the vectors of the base
curvilinear local.
Z
y
X
m
m
m
,
,
: linear couples acting around the axes
X
,
y
and
Z
:
D
H
has
D
rd
m
H
H
I
)
2
/
(
.)
(
.)
X
(
2
0
2
/
2
/
2
0
±
+
=
+
-
I
v
I
v
E
F
R
E
F
R
X
where
Z
y
X
E
E
E
,
,
are the vectors of the base
curvilinear local.
Z
y
X
,
,
: concentrated forces acting according to
X
,
y
and
Z
passing by the center of gravity of
transverse sections:
D
rd
H
H
I
+
-
=
2
/
2
/
2
0
.
I
C
E
F
where
Z
y
X
E
E
E
,
,
are the vectors of the local curvilinear base.
Z
y
X
µ
µ
µ
,
,
: moments concentrated around the axes
X
,
y
and
Z
:
+
-
=
2
/
2
/
2
0
.)
X
(
H
H
I
D
rd
µ
I
C
E
F
R
where
Z
y
X
E
E
E
,
,
are the vectors of the local curvilinear base.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
21/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
3.6.2 External efforts part and couples for the hull part
It is supposed that the external efforts applied to the elbow are independent the thickness of
bend:
D
W
v
U
D
D
M
M
W
F
v
F
U
F
W
L
R
X
L
X
X
R
X
ext.
+
+
+
+
+
+
+
+
+
=
2
0
,
0
0
2
0
]
[
)
(
S
where:
R
X
F
F
F
,
,
: surface forces acting according to
X
,
and
R
:
)
2
/
(
.
.
2
/
2
/
H
has
rd
F
H
H
I
±
+
=
+
-
I
S
I
v
E
F
E
F
where
R
X
E
E
E
,
,
are the vectors of the local toric base.
M
M
X
,
: surface couples acting around
X
and
:
)
2
/
(
).
2
/
(
).
(
2
/
2
/
H
has
H
D
R
M
H
H
I
±
±
+
=
+
-
I
S
R
I
v
R
E
F
E
E
F
E
X
X
where
R
X
E
E
E
,
,
are the vectors of the base
toric local.
R
X
,
,
: linear forces acting according to
X
,
and
R
:
+
-
=
2
/
2
/
.
H
H
I
rd
I
C
E
F
where
R
X
E
E
E
,
,
are the vectors of the local toric base.
,
X
: linear couples acting around
X
and
:
+
-
=
2
/
2
/
).
(
H
H
I
rd
I
C
R
E
F
E X
where
R
X
E
E
E
,
,
are the vectors of the local toric base.
Note:
When the external forces applied are independent of
external work on
kinematics of hull is null except that of the compressive forces corresponding to the forces
according to
R
E
. It is also noticed that the expressions of the moments linear and concentrated by
report/ratio with R are null. One finds although exerted moment ago
perpendicular to the plan of the hull.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
22/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
3.7 Efforts
generalized
If
S
is the surface of the transverse section
S
pipe, one poses:
dS
NR
S
xx
=
: normal effort in the center of gravity of the transverse section.
dS
dS
T
X
S
S
X
xy
y
)
cos
(sin
+
-
=
=
and
dS
dS
T
X
S
S
X
xz
Z
)
cos
(sin
-
=
=
sharp efforts according to
y
and
Z
.
dS
dS
Z
y
M
S
X
xy
xz
S
X
-
=
-
=
)
(
: torque around
X
.
dS
R
dS
Z
M
S
xx
S
xx
y
-
=
=
cos
: bending moment around
y
.
dS
R
dS
y
M
S
xx
S
xx
Z
=
-
=
sin
: bending moment around
Z
.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
23/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
4
Numerical discretization of the variational formulations
4.1 Discretization of the fields of displacement and deformation for
the beam part
In a point of average fiber, the field of displacement of beam is in the curvilinear reference mark
room defined in [§2.1]:
=
Z
y
X
Z
y
X
p
U
U
U
U
This field can be discretized in the following way:
()
=
+
+
=
NR
K
K
K
Z
K
K
y
K
kx
K
U
U
U
H
1
]
[
Z
y
X
U
and
()
=
+
+
=
NR
K
K
kz
K
ky
K
kx
K
H
1
]
[
Z
y
X
It should be noted that the nodal values are given in the local reference marks attached to the nodes and
that
U
and
must be expressed in the local reference mark associated the current point.
4.1.1 Beam
curve
One obtains then:
()
=
+
+
=
NR
K
K
kz
K
K
y
K
kx
K
K
y
K
K
X
K
Z
y
X
U
U
U
U
U
H
U
U
U
1
)
.
(
)
.
(
)
.
(
)
.
(
Z
y
y
y
X
X
y
X
X
and
()
=
+
+
=
NR
K
K
kz
K
ky
K
kx
K
ky
K
kx
K
Z
y
X
H
1
)
.
(
)
.
(
)
.
(
)
.
(
Z
y
y
y
X
X
y
X
X
According to the kinematics of beam presented higher to [§2.1]:
0
2
)
cos
sin
(
)
sin
cos
sin
cos
sin
sin
(
sin
1
2
)
sin
cos
(
)
cos
sin
cos
sin
cos
cos
(
sin
1
2
0
)
cos
sin
cos
(
sin
1
0
2
,
,
2
,
,
,
,
,
,
=
-
+
-
+
-
-
-
+
=
+
+
-
+
-
+
-
-
+
=
=
-
+
-
-
+
=
=
R
y
Z
Z
y
Z
y
X
R
y
Z
Z
y
X
Z
y
X
y
Z
X
y
X
rr
R
R
U
U
U
R
R
R
R
R
U
U
U
R
R
R
R
R
U
U
R
R
Knowing that
X
y
y
X
=
-
=
,
,
and
with moreover
K
K
K
K
K
K
K
K
S
C
=
-
=
-
=
=
-
=
=
)
sin (
.
.
)
cos (
.
.
y
X
X
y
y
y
X
X
and
.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
24/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
That implies for the field of deformation:
sin
))
cos (
)
sin (
(
cos
]
sin
))
cos (
)
sin (
(
cos
sin
cos
))
sin (
)
cos (
(
sin
))
cos (
)
sin (
(
sin
))
sin (
)
cos (
(
sin
[
sin
1
cos
))
cos (
)
sin (
(
sin
]
cos
sin
))
cos (
)
sin (
(
cos
))
cos (
)
sin (
(
))
sin (
)
cos (
(
sin
))
sin (
)
cos (
(
cos
))
cos (
)
sin (
(
cos
))
sin (
)
cos (
(
cos
[
sin
1
0
]
sin
))
cos (
)
cos (
(
cos
))
cos (
)
sin (
(
cos
))
sin (
)
cos (
(
cos
))
cos (
)
sin (
(
))
cos (
)
sin (
(
))
sin (
)
cos (
(
[
sin
1
2
1
2
1
1
kz
K
K
ky
K
kx
K
kz
K
K
ky
K
kx
K
kz
K
K
K
y
K
kx
K
K
K
y
K
kx
K
K
K
y
K
kx
K
NR
K
R
kz
K
K
ky
K
kx
K
kz
K
K
ky
K
kx
K
K
ky
K
kx
K
K
ky
K
kx
K
kz
K
K
K
y
K
K
X
K
K
K
y
K
kx
K
K
K
y
K
kx
K
NR
K
kz
K
K
ky
K
kx
K
K
ky
K
kx
K
K
ky
K
kx
K
K
K
y
K
K
X
K
K
K
y
K
kx
K
K
K
y
K
NR
K
kx
K
H
H
H
R
H
R
U
H
U
U
H
U
U
H
U
U
H
R
R
H
H
H
R
H
R
H
R
H
R
U
H
U
U
H
U
U
H
U
U
H
R
R
H
R
H
R
H
R
H
R
U
U
H
U
U
H
U
U
H
R
R
+
-
+
-
-
-
-
+
-
+
-
-
-
-
-
-
+
-
-
-
-
-
-
+
=
+
-
+
-
+
-
-
+
-
+
-
-
-
-
-
-
-
-
-
+
-
-
-
-
-
+
-
-
-
-
-
-
+
=
=
+
-
-
-
-
-
+
-
-
-
-
-
-
-
+
-
-
-
-
-
-
+
-
-
-
+
=
=
=
=
Maybe in matric form:
=
=
NR
K
P
K
P
K
P
P
P
1
U
B
where
=
kz
ky
kx
kz
K
y
kx
p
K
U
U
U
U
is the field of displacement to the node
K
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
25/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
and
B
K
P
K
K
K K
K K
K
K
K K
K
K
K
K
K K
K
K K
K
K
K
H C
R R
H S
R R
H S
R R
H C
R R
R H S
R R
R H C
R R
R H C
R R
R H S
R R
R H
R R
H S
R R
H C
R R
H
=
+
-
+
-
+
-
+
-
+
-
+
-
+
+
+
+
-
+
-
+
-
sin
sin
sin
sin
cos
sin
cos
sin
cos
sin
cos
sin
sin
sin
cos
sin
cos
sin
2
2
0
2
2
0
0
0
0
0
0
2
C
R R
H S
R R
H
R R
H S R
R
R R
R H C
R R
H C R
R
R R
R H S
R R
RH
R R
H S
R R
H C
R R
H C
R R
H S
R R
K
K K
K
K K
K
K
K
K
K K
K
K K
K
K
K
K
K K
cos
sin
cos
sin
sin
sin
(sin
)
sin
sin
(sin
)
sin
sin
cos
sin
sin
sin
sin
sin
sin
sin
sin
+
+
+
+
+
+
-
+
+
+
+
+
+
-
+
-
+
-
+
+
+
2
2
2
2
2
sin
cos
sin
cos
sin
cos
sin
sin
sin
-
+
-
+
-
+
+
H
R R
RH S
R R
RH C
R R
RH
R R
K
K K
K
K
K
The matrix of passage of the deformations to the field of displacement is written as follows:
(
)
P
P
P
NR
B
B
B
L
1
=
4.1.2 Beam
straight line
()
=
=
NR
K
kz
K
y
kx
K
Z
y
X
U
U
U
X
H
U
U
U
1
and
()
=
=
NR
K
kz
ky
kx
K
Z
y
X
X
H
1
According to the kinematics of beam presented higher [§2.1]:
0
2
cos
)
(
sin
)
(
2
cos
)
(
sin
)
(
2
0
cos
sin
0
,
,
,
,
,
,
,
,
=
+
-
-
=
-
+
+
+
-
=
=
-
+
=
=
R
X
Z
y
X
y
Z
X-ray
X
y
Z
X
Z
y
X
X
X
X
y
X
Z
X
X
xx
rr
U
U
U
U
R
R
R
U
that implies for the field of deformation:
()
()
()
()
()
()
()
()
()
()
(
)
=
=
=
+
-
-
-
=
+
+
-
+
-
=
=
+
-
=
NR
K
kz
K
ky
K
kz
K
K
y
K
X-ray
NR
K
kz
K
ky
K
kx
K
kz
K
K
y
K
X
NR
K
kz
K
ky
K
K
X
K
xx
H
H
U
H
U
H
H
H
R
H
U
H
U
H
R
H
R
H
U
H
1
'
'
1
'
'
'
1
'
'
'
sin
cos
cos
sin
cos
sin
sin
cos
0
sin
cos
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
26/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Maybe in matric form:
=
=
NR
K
P
K
P
K
P
X
P
X
P
P
X
1
U
B
where
=
kz
ky
kx
kz
K
y
kx
p
K
U
U
U
U
is the field of displacement to the node
K
and:
-
-
-
-
-
-
=
K
K
K
K
K
K
K
K
K
K
K
K
P
H
H
H
H
H
H
H
R
H
H
H
R
H
R
H
K
)
sin (
)
cos (
0
)
cos (
)
sin (
0
)
cos (
)
sin (
)
sin (
)
cos (
0
0
0
0
0
0
0
)
sin (
)
cos (
0
0
0
B
The matrix of passage of the deformations to the field of displacement is written as follows:
(
)
P
P
P
NR
B
B
B
L
1
=
4.2 Discretization of the fields of displacement and deformation for
the additional part
One discretizes the field of displacement for the surface of the pipe in the form:
()
=
=
NR
K
S
K
K
S
X
H
1
U
U
with:
=
O
O
I
om
om
om
im
im
im
S
W
W
W
W
v
U
W
v
U
1
1
U
and
M
m
W
W
W
W
v
U
W
v
U
O
K
O
K
ik
om
K
om
K
om
K
I m
K
I m
K
I m
K
S
K
,
2
1
1
=
=
U
.
One has as follows:
4
4
4
4
4
4
4
4
4
4
4
3
4
4
4
4
4
4
4
4
4
4
4
2
1
M
m
m
m
m
m
m
m
X
W
X
v
X
U
,
2
1
)
sin (
)
cos (
)
sin (
0
0
)
cos (
0
0
0
)
cos (
)
sin (
0
)
cos (
0
0
)
sin (
0
0
0
0
0
0
)
sin (
0
0
)
cos (
)
,
(
)
,
(
)
,
(
=
-
=
S
U
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
27/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
if indices
m
of
S
K
U
are ordered in the following way:
=
=
=
=
=
=
=
=
=
=
=
=
=
O
K
O
K
ik
O
M
m
K
O
M
m
K
O
M
m
K
I
M
m
K
I
M
m
K
I
M
m
K
om
K
om
K
om
K
I m
K
I m
K
I m
K
S
K
W
W
W
W
v
U
W
v
U
W
v
U
W
v
U
1
1
2
2
2
2
2
2
M
U
The kinematics of hull presented higher to [§2.2] is:
0
2
0
2
2
2
=
=
=
=
+
=
+
=
+
=
E
E
E
E
E
X
4.2.1 Bend
With:
)
cos
(
sin
1
1
2
)
(
1
)
sin
cos
(
sin
1
U
v
R
R
U
R
E
W
v
R
E
W
v
U
R
R
E
-
+
+
=
+
=
+
+
+
=
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
28/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
and:
)
sin
(
1
)
cos
sin
(
)
sin
(
1
]
)
sin
(
1
)
sin
(
1
[
]
)
sin
(
cos
)
sin
) (
sin
(
cos
) [
sin
(
2
)
(
1
)]
(
cos
)
sin
(
sin
1
[
sin
1
2
2
2
2
2
2
has
R
R
U
U
R
R
has
v
has
R
R
R
R
has
W
has
R
R
has
has
R
R
R
U
W
v
W
rear
v
W
has
U
W
has
R
R
R
+
+
+
+
+
+
+
+
-
+
+
+
+
-
=
-
-
=
-
-
-
+
-
+
=
allows to break up the field of deformation of hull on the modes of Fourier in the way
following:
=
=
NR
K
S
K
S
K
sx
sx
S
sxx
1
U
B
with:
(
)
sg
K
so M
m
K
som
K
if
M
m
K
sim
K
S
K
B
B
B
B
B
B
=
=
=
=
=
L
L
2
2
where
B
km
if
K
K
K
K
K
K
K
K
H
m
R R
R has
H
m
R R
has
H
m
R R
H
m
R R
R has
H m
m
R R has
m
R H
has
m
R H
m
has
m
m
R H
m
H
=
+
+
+
+
+
+
-
+
+
+
+
+
+
-
-
cos
sin (
sin
(
sin))
sin
cos
sin
cos
sin
sin
cos
(
sin) (
sin)
sin
cos
(
sin)
cos
cos
sin
1
1
0
1
1
1
2
K
K
K
K
K
K
K
m
R R
R has
H
m
R R
aH
m
R R has
H
m
m
m
R R has
H
has
m
R R
mH
m
R R has
mH
m
R R has
R has
m
cos
cos
sin
cos sin
(
sin) [
cos
(
sin)
cos
(
sin)]
(cos
cos
sin
sin)
(
sin)
sin
sin
sin
(
sin)
sin
(
sin)
cos
sin [
cos
+
-
+
+
+
+
+
-
+
+
+
+
+
+
+
+
1
H
R R
has
m H
R R has
K
K
(
sin)
cos
(
sin)]
+
+
+
0
0
0
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
29/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
B
km
so
K
K
K
K
K
K
K
K
H
m
R R
R has
H
m
R R
has
H
m
R R
H
m
R R
R has
H m
m
R R has
m
R H
has
m
R H
m
has
m
m
R H
m
H
=
+
+
+
+
+
+
-
+
+
-
+
-
+
+
-
sin
sin [
sin
(
sin)]
cos
sin
cos
sin
sin
sin
sin
(
sin) (
sin)
cos
cos
(
sin)
sin
sin
cos
1
1
0
1
1
1
2
K
K
K
K
K
K
K
K
m
R R
R has
H
m
R R
aH
m
R R has
H
m
m
m
R R has
H
has
m
R R
mH
m
R R has
mH
m
R R has
R has
H
m
sin
cos
sin
cos sin
sin [
sin
sin
sin
(
sin)]
(cos
sin
sin
cos)
(
sin)
cos
sin
cos
(
sin)
cos
(
sin)
cos
sin [
sin
+
-
+
+
+
+
+
+
+
+
+
-
+
+
+
+
+
1
R R
aH
m
R R has
K
+
+
+
sin
sin
(
sin)]
0
0
0
and
B
K
sg
K
K
K
K
K
K
K
K
K
H
R R
H
R R
R has
H
R R has
H
R R
H
R R
R has
H
R R has
H
R R
H
R has
R R
R
H has
=
+
-
+
+
+
+
-
+
-
+
+
-
+
+
-
+
+
+
2
2
2
2 1
2
2
2
cos sin
sin
cos
(
sin) (
sin)
cos sin
(
sin)
(sin
cos
)
sin
sin
(
sin) (
sin)
cos
(
sin)
sin
(
sin)
(
sin) (
sin)
cos
sin
sin [
sin
(
sin)]
cos
sin [
cos
sin
cos
(
sin)]
cos [
sin
(
sin)]
sin
sin [
cos
sin
cos
(
2 1
1 2
1 2
R
H has
H
R
H
has
R R
R R has
H
R has
R R
has
R R has
H
has
R R
R R has
H
R has
R R
has
R R has
K
K
K
K
K
K
+
+
+
+
+
+
+
+
+
+
-
+
+
+
+
+
+
+
+
+
sin)]
cos
(
sin) (
sin)
cos
(
sin)
+
+
+
+
H
R R
R has
H has
R R has
K
K
2
0
0
0
4.2.2 Pipe
right
With:
X
R
has
X
U
R
X
v
E
W
v
R
E
X
U
E
X
X
X
xx
X
xx
-
=
-
=
=
+
=
+
=
=
1
2
1
1
2
)
(
1
)
(
1
v
W
has
X
W
X
-
=
-
=
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
30/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
the field of deformation of hull breaks up on the modes of Fourier in the following way:
(
)
(
)
(
)
(
)
(
)
(
)
(
)
(
)
(
)
(
)
0
sin
cos
)
(
cos
sin
cos
sin
)
(
sin
cos
)
(
1
cos
sin
cos
sin
)
(
sin
cos
)
(
cos
sin
cos
sin
)
(
sin
cos
)
(
1
cos
sin
cos
sin
)
(
1
sin
cos
)
(
sin
cos
)
(
1
1
2
1
2
1
1
1
2
1
2
1
1
1
1
2
2
1
2
1
1
1
1
1
2
1
1
1
1
2
2
1
2
=
+
+
+
-
+
+
-
+
+
+
+
+
-
=
+
+
-
+
+
-
+
+
+
+
+
+
-
=
+
+
-
+
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
X
NR
N
M
m
O
Nm
inm
O
N
N
NR
N
M
m
onm
inm
O
N
in
N
NR
N
M
m
O
Nm
inm
N
NR
N
M
m
O
Nm
inm
O
N
in
N
X
NR
N
M
m
O
Nm
inm
O
N
N
NR
N
M
m
O
Nm
inm
O
N
in
N
NR
N
M
m
O
Nm
inm
O
N
N
NR
N
M
m
O
Nm
inm
O
N
in
N
NR
N
M
m
O
Nm
inm
one
N
NR
N
M
m
O
Nm
inm
N
xx
W
m
W
m
W
X
H
X
R
R
v
m
v
m
W
W
X
H
X
has
U
m
U
m
X
H
R
v
m
v
m
W
W
X
H
X
W
m
W
m
W
X
H
R
has
v
m
v
m
W
W
X
H
R
has
W
m
W
m
W
X
H
R
v
m
v
m
W
W
X
H
R
W
m
W
m
W
X
H
X
U
m
U
m
X
H
X
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
31/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
That is to say still:
(
)
(
)
(
)
(
)
(
)
(
)
(
)
(
)
(
)
(
)
0
cos
sin
)
(
cos
sin
cos
sin
)
(
cos
sin
)
(
1
cos
sin
cos
sin
)
(
sin
cos
)
(
sin
cos
sin
cos
)
(
sin
cos
)
(
1
sin
cos
sin
cos
)
(
1
sin
cos
)
(
sin
cos
)
(
1
1
1
2
1
1
1
2
1
2
1
1
1
1
2
2
1
2
1
1
1
1
1
2
1
1
1
1
1
2
=
+
-
+
-
+
-
+
+
-
+
+
-
=
+
+
-
+
+
+
+
+
+
-
+
+
=
+
+
-
+
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
X
NR
N
M
m
O
Nm
inm
N
NR
N
M
m
O
Nm
inm
O
N
in
N
NR
N
M
m
O
Nm
inm
N
NR
N
M
m
onm
inm
O
N
in
N
X
NR
N
M
m
O
Nm
inm
N
NR
N
M
m
O
Nm
inm
one
in
N
NR
N
M
m
O
Nm
inm
O
N
N
NR
N
M
m
O
Nm
inm
O
N
in
N
NR
N
M
m
O
Nm
inm
O
N
N
NR
N
M
m
O
Nm
inm
N
xx
W
m
m
W
m
m
X
H
has
R
v
m
v
m
W
W
X
H
has
U
m
m
U
m
m
X
H
R
v
m
v
m
W
W
X
H
W
m
m
W
m
m
X
H
R
has
v
m
m
v
m
m
W
W
X
H
R
has
W
m
W
m
W
X
H
R
v
m
m
v
m
m
W
W
X
H
R
W
m
W
m
W
X
H
U
m
U
m
X
H
This gives in matric form:
=
=
NR
K
S
K
S
K
sx
sx
S
sxx
1
U
B
with:
(
)
sg
K
so M
m
K
som
K
if
M
m
K
sim
K
S
K
B
B
B
B
B
B
=
=
=
=
=
L
L
2
2
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
32/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
where:
+
+
-
+
+
-
=
0
0
0
)
sin (
1
)
sin (
)
sin (
1
)
cos (
1
1
)
cos (
0
)
cos (
0
)
cos (
2
m
H
m
has
R
has
m
H
m
H
R
m
has
m
m
H
R
has
m
H
R
m
m
H
m
H
K
K
K
K
K
K
K
if
km
B
,
+
-
+
+
+
-
-
=
0
0
0
)
cos (
1
)
cos (
)
cos (
1
)
sin (
1
1
)
sin (
0
)
sin (
0
)
sin (
2
m
H
m
has
R
has
m
H
m
H
R
m
has
m
m
H
R
has
m
H
R
m
m
H
m
H
K
K
K
K
K
K
K
so
km
B
and
+
+
-
+
+
+
+
-
-
-
=
0
0
0
0
)
cos (
2
1
)
sin (
2
1
)
sin (
1
2
)
cos (
1
2
)
sin (
)
cos (
K
K
K
K
K
K
K
K
sg
K
H
R
has
H
R
has
R
H
H
has
R
H
has
R
H
H
H
B
4.3
Discretization of the field of total deflection
=
=
=
=
=
+
=
+
=
NR
K
K
K
NR
K
S
K
S
K
NR
K
P
K
P
K
sx
sx
S
sxx
P
X
P
X
P
P
xx
X
X
xx
1
1
1
B U
U
B
U
B
U
B
with
(
)
NR
K
S
K
P
K
,
1
=
=
B
B
B
and
NR
K
S
K
P
K
,
1
=
=
U
U
U
4.4
Stamp rigidity
The variational formulation of the work of deformation is:
(
)
dxd
rd
W
L
H
H
X
X
xx
X
X
xx
def
-
=
0
2
0
2
/
2
/
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
33/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
that is to say still:
(
)
dxd
rd
W
L
X
X
xx
X
X
xx
H
H
def
=
-
0
2
0
2
/
2
/
C
=
=
=
=
-
=
=
-
=
=
-
NR
T
NR
T
L
NR
T
T
NR
T
H
H
L
NR
K
K
K
NR
K
T
K
T
K
H
H
L
NR
K
K
K
T
NR
K
K
K
H
H
def
dxd
rd
dxd
rd
dxd
rd
W
U
U
K
U
U
U
U
B
C
B
U
U
U
B
C
B
U
U
B
C
U
B
M
L
M
L
1
1
0
2
0
1
1
2
/
2
/
0
2
0
1
1
2
/
2
/
0
2
0
1
1
2
/
2
/
The principle of virtual work is written then
U
F
KU
U
T
=
where
K
is the matrix of rigidity which
is worth:
{
}
dxd
rd
L
H
H
-
=
0
2
0
2
/
2
/
CB
B
K
T
Note:
One makes no assumption on the law of behavior. This expression is thus in
private individual validates the nonlinear behaviors in the case of (plasticity).
4.5
Stamp of mass
The terms of the matrix of mass are obtained after discretization of the variational formulation
following of the noncentrifugal terms of inertia:
D
rdxd
W
L
H
H
ac
farmhouse
v
U.
0
2
0
2
/
2
/
+
-
=
&&
with
(
)
(
)
(
)
(
)
(
)
(
)
=
R
X
W
R
X
v
R
X
U
R
X
U
R
X
U
R
X
U
,
,
,
,
,
,
,
,
,
,
,
,
3
2
1
U
.
The notations used are those of [§2.1]:
3
2
1
,
U
U
U
and
are displacements of beam in a point
section and
W
v
U
and
,
are displacements of average fiber of this section in this same
not.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
34/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
The discretization gives then:
M
m
W
W
W
W
v
U
W
v
U
U
U
U
H
NR
K
ok
O
K
ik
O
km
O
km
O
km
ikm
ikm
ikm
K
Z
ky
K
X
kz
K
y
kx
K
K
,
2
1
1
1
=
=
=
NR
U
where matrices
NR
K
have as an expression:
NR
X X there X
X y
y y
X there there y
X X
y X
X X
y X
K
K
K
K
K
K
K
K
K
K
K
R
R
R
R
R
R
R
m
m
m
m
=
-
-
-
-
-
.
.
cos (
. )
cos (
. )
sin
.
.
cos (
. )
cos (
. )
sin (. )
sin (
. )
cos (
)
sin (
)
sin (
)
cos (
)
sin ()
cos ()
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
L
L
L
L
L
L
cos (
)
sin (
) cos ()
sin ()
,
m
m
m
M
0
0
1
2
=
1
2
444444444444
3
444444444444
The matrix of mass has then as an expression:
D
rdxd
L
H
H
T
+
-
=
0
2
0
2
/
2
/
NR
NR
M
.
with
(
)
NR
K
K
K
H
,
1
=
=
NR
NR
.
Note:
In the case of the right pipe, one has
1
.
.
=
=
y
y
X
X
K
K
and
0
.
.
=
=
X
y
y
X
K
K
.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
35/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
4.6
Functions of form
One chooses at least quadratic functions of form for the part beam (displacements and
rotations) in order to avoid the phenomena of numerical blocking [bib3]. This choice implies the use of one
finite element with three or four nodes. In the case of an element with 3 nodes, the functions of form are
quadratic, and for an element with 4 nodes, the functions of form will be cubic. For the part
additional, one chooses to take the same functions of form as for the beam part.
The quadratic functions of form (element with 3 nodes) are as follows:
()
()
()
-
-
=
-
=
-
-
=
1
4
1
2
1
1
2
3
2
1
L
X
L
X
X
H
L
X
L
X
X
H
L
X
L
X
X
H
The cubic functions of form (element with 4 nodes) are related to Lagrange of command 3:
() (
) (
) (
)
(
) (
) (
)
() (
) (
) (
)
(
) (
) (
)
() (
) (
) (
)
(
) (
) (
)
() (
) (
) (
)
(
) (
) (
)
1
1
3
4
2
4
1
4
3
2
1
4
4
3
2
3
1
3
4
2
1
3
4
2
3
2
1
2
4
3
1
2
4
1
3
1
2
1
4
3
2
1
-
-
-
-
-
-
-
=
-
-
-
-
-
-
=
-
-
-
-
-
-
=
-
-
-
-
-
-
=
X
H
X
H
X
H
X
H
4.7 Integration
numerical
Integration is done by the method of Gauss along average fiber, the method of Simpson
in the thickness and on the circumference. For the integration of Gauss, one uses 3 points of integration
for the elements with 3 nodes, as for the elements with 4 nodes (those under-are thus integrated).
Integration in the thickness is an integration by layers of which the number could be fixed
later on by the user. For each layer one takes 3 points of Simpson, the 2 points
ends being common with the close layers. Thus for N layers one uses 2n+1 points.
a many sectors for integration on the circumference, could also be fixed later on by
the user. Currently the numbers of layers and sectors are fixed at their maximum value: 3
layers (7 points) and 16 sectors (33 points), which gives 693 points of integration on the whole.
The integration of Simpson amounts calculating the sum of the values of the function at the points
integrations (ends and medium of each layer or sector) affected of the weights given
by the table below.
Cordonnées of the points
Weight
-
3 5
/
= - 0,77459 66692 41483
5/9=0,55555 55555 55556
0
8/9=0,88888 88888 88889
3 5
/
=0,77459 66692 41483
5/9=0,55555 55555 55556
Integration of Gauss on average fiber
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
36/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
1
3
4
3
2
3
4
3
2
3
.....
2
3
4
3
2
3
4
3
1
3
Weight of the points of integration for the method of Simpson
Thus for a function
)
,
,
(
X
F
on a right geometry one a:
[
]
=
+
=
+
=
-
-
-
=
=
NPG
K
NR
N
NR
m
NR
NR
NPG
NR
m
N
K
SECT
NECK
H
H
L
H
H
NECK
SECT
NECK
SECT
NECK
X
F
R
W
W
W
NR
NR
H
L
D
D
X
D
R
X
F
L
D
D
dx
R
X
F
1
1
2
1
1
2
1
1
1
2
0
2
/
2
/
0
2
0
2
/
2
/
)
,
,
~
(
~
2
2
2
2
~
)
,
,
~
(
~
2
)
,
,
(
m
N
K
W
W
W
being respectively weights of integration over the length, the circumference and
in the thickness, ordered as the two tables show it above.
4.8
Discretization of external work
The variational formulation of external energy for the beam part is written
:
L
Z
Z
y
y
X
X
Z
Z
y
y
X
X
L
Z
Z
y
y
X
X
Z
Z
y
y
X
X
ext.
U
U
U
dx
m
m
m
U
F
U
F
U
F
W
,
0
0
]
[
)
(
µ
µ
µ
+
+
+
+
+
+
+
+
+
+
+
=
p
and for the additional part she is written by taking into account only the loading of
pressure:
L
R
L
R
ext.
W
dx
W
F
W
,
0
0
]
[
+
=
S
By taking account of the discretization of displacements, one can write:
W
F H X
U
U
F H X
U
U
F H X U
m H X
m H X
m H X
dx
H X U
ext.
X
K
xk
y
K
y
K
xk
y
K
Z
K
zk
L
K
NR
X
K
xk
yk
y
K
xk
yk
Z
K
zk
X
K
xk
p
K
K
K
K
K
K
K
K
X X
y X
X y
y y
X X
y X
X y
y y
=
+
+
+
+
+
+
+
+
+
+
+
=
(
() (.
.
)
() (.
.
)
()
() (
.
.
)
() (
.
.
)
()
)
[
()
0
1
(
)
(
)
µ
µ
µ
µ
y
K
y
K
Z
K
zk
X
K
xk
y
K
yk
Z
K
zk
L
K
L
X
y
X
K
L
K
X
y
L
y
K
L
Z
K
X
K
H X U
H X U
H X
H X
H X
H X F
F
dx
H X
H X F
F
dx
H X
F H X dx
H X
()
()
()
()
()
]
()
(
. )
(. )
[
()]
()
(
. )
(
. )
[
()]
()
[
(
,
,
,
+
+
+
+
=
+
+
+
+
+
0
0
0
0
0
X X
X y
y X
y y
K
K
K
K
(
)
(
)
)]
()
(
. )
(. )
[
()]
()
(
. )
(
. )
[
()]
()
[
()]
,
,
,
,
0
0
1
0
0
0
0
0
0
1
L
L
K
NR
K
X
y
L
X
K
L
K
X
y
L
y
K
L
Z
K
L
Z
K
L
K
p
K
NR
kP
K
p
P
p
H X m
m
dx
H X
H X m
m
dx
H X
m H X dx
H X
+
+
+
+
+
=
=
=
=
X X
X y
y X
y y
U
F U
F U
K
K
K
K
µ
µ
µ
and
W
F H X W dx
H X W
F H X dx
H X
ext.
K
NR
R
K
KB
L
R
K
KB
L
R
K
L
R
K
L
K
NR
K
S
ks
K
S
K
NR
S
S
S
U
F U
F U
=
+
=
+
=
=
=
=
=
1
0
0
0
0
1
1
0 0 0 0 0 0
()
[
()
]
()
[
()]
,
,
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
37/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Note:
For the extreme nodes of the elbow one has
1
.
.
=
=
y
y
X
X
K
K
and
0
.
.
=
=
X
y
y
X
K
K
. In the case of it
right pipe, one has
1
.
.
=
=
y
y
X
X
K
K
and
0
.
.
=
=
X
y
y
X
K
K
for all the element.
5
Geometrical characteristics of the pipe section
One has in this chapter some useful results to characterize the element pipe and which are
calculated by the option of calculation
MASS_INER
of Code_Aster. In the continuation the index
D
indicate them
results for the right pipe and the index
C
for the curved pipe.
·
Volume:
ah
R
D
D
D
has
has
R
D
D
dx
R
V
lah
D
D
dx
R
V
H
H
L
H
H
C
L
H
H
D
=
+
+
+
=
=
=
=
-
-
-
2
)
(
]
sin
)
(
[
2
0
2
0
2
/
2
/
0
2
0
2
/
2
/
0
2
0
2
/
2
/
.
·
Center of gravity: The position of this last is calculated starting from the point medium to both
extreme nodes of the pipe section, in the reference mark associated with the node interns element
(cf [§2.1.1]). In this reference mark, the co-ordinates of the center of gravity are:
=
0
0
0
Gd
Gd
Gd
Z
y
X
and
-
+
+
-
=
0
2
cos
])
4
[
2
1
1
(
2
sin
2
0
2
2
2
H
has
R
R
Z
y
X
Gc
Gc
Gc
·
Stamp inertia: It is relatively easy to calculate the matrix of inertia in the center of
curvature of the elbow
O
in the reference mark defined above. To have his expression one uses
then the fact that:
2
2
)
(
)
(
)
(
)
(
)
(
)
(
mb
O
I
G
I
O
I
G
I
mb
O
I
G
I
zz
zz
yy
yy
xx
xx
-
=
=
-
=
where B is the distance between the center of gravity and the center of curvature which is worth:
])
4
[
2
1
1
(
2
sin
2
2
2
2
H
has
R
R
B
+
+
=
.
In the case of a right pipe, the concept of center of curvature does not have a direction.
O
and
G
are confused
with the node interns element and the medium of the segment joining the two nodes node.
If one notes
With
the surface of the transverse section and
I
its inertia compared to the center of the section one can
to write:
)
12
/
2
/
(
)
(
)
12
/
2
/
(
)
(
)
(
2
2
Al
I
L
O
I
Al
I
L
O
I
Li
O
I
D
zz
D
yy
D
xx
+
=
+
=
=
and
)
2
3
(
)
(
])
4
sin
2
1
] [
2
3
[
2
(
)
(
])
4
sin
2
1
] [
2
3
[
2
(
)
(
2
2
2
I
AR
R
O
I
I
AR
I
R
O
I
I
AR
I
R
O
I
czz
C
yy
cxx
+
=
-
+
+
=
+
+
+
=
.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
38/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
6 Connection
pipe-pipe
In order to be able to represent a line of piping correctly where the elbows are not
coplanar, it is necessary to choose a common origin of
. Thus for two elbows belonging to two
perpendicular plans between them, it is necessary to be able to take into account the fact that displacements in
the plan of the first elbow are equal to displacements except plan of the second in the cross-section of
connection.
·
·
·
·
· ·
1
2
Generating line
Appear 6-a: Representation of two noncoplanar elbows connected by a right pipe
In [bib12], this common origin is defined by a generating line continues along
piping as indicated above. This generator intersects each transverse section in one
not. The angle enters
Z
defined on [Figure 2.1.1-a] and line passing by the center of the section
transverse and this point are worth
.
6.1
Construction of a particular generator
For a transverse section end of the line of piping, one defines a vector origin
1
Z
unit in the plan of this section. The intersection enters the direction of this vector and surface
average of the elbow determines the trace of the generator on this section. One calls
1
1
1
Z
y
X
,
,
direct trihedron associated this section where
1
X
is the unit vector perpendicular to the section
transverse built with [Figure 2.1.1-a]. For the whole of the other transverse sections, the trihedron
K
K
K
Z
y
X
,
,
is obtained either by rotation of the trihedron
1
K
1
K
1
K
Z
y
X
-
-
-
,
,
in the case of bent parts,
maybe by translation of the trihedron
1
K
1
K
1
K
Z
y
X
-
-
-
,
,
for the right parts of piping. The intersection
between the transverse section and the line resulting from the center of this section directed by
K
Z
is the trace
of a generator represented below.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
39/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Z
X
Y
O
Z
x2
Z
y2
x1
Z
y1
Z
With
B
generator
1
2
Appear 6.1-a: Representation of the generator of reference
The origin of
commune with all the elements is defined compared to the trace of this generator
on the transverse section. The angle enters the trace of the generator and the current position on the section
transverse is then called
.
6.2
Connection from one element to another
The kinematics of [§3.1] is given in the plan of the elbow. The aforementioned is determined by the arc of circle
generated by the axis of the elbow. The origin of the angles is the normal in the plan chosen as with [§2.1].
To define the origin starting from a generator makes it possible to raise the problems of continuity of
displacements of an element to another. Indeed if one postulates that relative displacements of
transverse sections are of the type
=
+
M
p
COp
IP
p
U
p
U
1
sin
cos
where
is the angle with the trace of
generator on the transverse section, the continuity of displacements is automatically assured
from one element to another.
One notes
Z
the vector perpendicular to the plan of the elbow corresponding at the selected origin of the angles
up to now. It is noticed that the vectors
Z
and
K
Z
are in the plan of the section
K
.
is the angle
defined compared to
Z
. If one introduces
the angle counted starting from the trace of the generator on
transverse section (thus compared to
K
Z
) one with the following relation:
K
-
=
where
)
,
(
K
K
Z
Z
=
angle enters
Z
and
K
Z
in the plan of the transverse section. Thus displacements are from now on
type
=
-
+
-
M
p
K
COp
K
IP
p
U
p
U
1
)
(
sin
)
(
cos
. It should be noted that for an elbow given the angle
K
whatever the selected transverse section is identical. It is at the time of the passage from one elbow to another
that the value of
K
change.
Note:
When piping consists of right elements colinéaires, one chooses arbitrarily
0
=
.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
40/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
6.3
Numerical establishment
The line of piping is with a grid by right or curved elements to order. The first element
indicate the beginning of the line of piping. One determines for this element the associated trihedron
1
1
1
Z
y
X
,
,
. If this element is right, one chooses
0
1
=
, if not one calculates
1
as indicated in
preceding paragraph. If the first element is right the trihedron associated with the first section
transverse of the second element
2
2
2
Z
y
X
,
,
is obtained by translation of
1
1
1
Z
y
X
,
,
. If the first
element is curved, the associated trihedron
2
2
2
Z
y
X
,
,
is obtained by rotation of
1
1
1
Z
y
X
,
,
in the plan
elbow. In this case
0
2
=
if the second element is right and
)
,
(
2
2
2
Z
Z
=
if the second
element is curved where
2
Z
is built like
Z
[Figure 2.1.1-a]. The continuation of construction
results easily by recursion from the preceding diagram.
7
Connections hull-pipe and 3D-pipe
7.1 Step
followed
One adopts here a step similar to the cases 3D-beam [R3.03.03], and hull-beam [R3.03.06]: it
acts to characterize the connection between a node end of an element pipe and a group of meshs of
edge of elements of hulls or 3D. This makes it possible to net part of piping (for example one
bend) in hulls or elements 3D, and the remainder in right pipes.
Appear 7.1-a: Connection between a right mesh COQUE_3D and pipes [HI75-98/001]
Thanks to the kinematics introduced into the element pipe, the connections hull - pipe and 3D - pipe must
to allow to net in elements of hulls or 3D only the elbow, without right parts,
since the damping of ovalization (and the roll) is taken into account in the element
pipe.
The connection results in relations kinematics between the degrees of freedom of the nodes of
S
(which
represent the section of connection, modelized by elements of edge of hull or 3D), and it
node NR of pipe.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
41/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
So that the connection is effective, it is necessary [R3.03.03] that it checks the following properties:
1) to be able to transmit efforts of beam to the mesh hulls or 3D, and to be able to transmit
all degrees of freedom of the element pipe (or duaux efforts of those),
2) not to generate in the elements of hulls or 3D of secondary stresses,
3) not to support the static relations kinematics or conditions ones by report/ratio
with the others,
4) to admit unspecified behaviors and to function in dynamics.
The linear relations will have the same form as in the case hull-beam, with equations
additional specific to the degrees of freedom of the pipe.
One already introduced with [§3.1] space
T
fields associated with a torque (defined by two vectors):
(
)
()
{
}
GM
v
v
T
+
=
=
T
T
M
V
that
such
,
/
where for the fields of displacement of
T
,
T
is the translation of the section (or point
G
),
infinitesimal rotation and fields
v
are displacements preserving the section
S
plane and not
deformation there (One uses still the assumptions of NAVIER-BERNOULLI).
The displacement of the pipe is worth then:
+
=
T
U
T
U
U
U
U
S
p
S
p
T
,
where:
=
=
T
W
W
v
v
T
0
.
/
S
V
The step consists in breaking up the field of displacement of hull
C
U
or the field of
three-dimensional displacement
D
3
U
in three fields:
U
U
U
U
p
+
+
=
S
C
·
a field of displacement following a kinematics of beam
p
U
(torque),
·
a field of local displacement of the section according to the kinematics of pipe
S
U
(series of
Fourier) defined in [§3.1],
·
an additional field
U
orthogonal with the two first within the meaning of the scalar product.
Note:
When the decomposition in Fourier series of [§3.1] is infinite one has
0
=
U
.
To translate the equation above into linear relations, it is shown that the integrals should be calculated
following, for the hull (or the 3D) and the pipe:
·
average displacement:
S
dS
C
U
·
average rotation:
S
dS
C
U
GM
·
average swelling:
S
dS
C
U
GM
GM.
·
modes of Fourier
:
S
S
dS
p
dS
p
sin
,
cos
C
C
U
U
,
S
dS
p
cos
GM
GM
U
C
,
S
dS
p
sin
GM
GM
U
C
,
S
dS
p
cos
.
C
U
GM
GM
,
S
dS
p
sin
.
C
U
GM
GM
. For
modes of Fourier one will choose the relations simplest to exploit since some are
redundant (see remark of [§6.7]).
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
42/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Note:
One passes easily from the analytical expressions of the connection 3D - pipe to those of the connection
hull - pipe. It is enough to substitute
D
3
U
with
C
U
in the whole of the integrals of
connection proposed above. One will thus speak again of this connection only to [§7.8] treating
numerical establishment.
7.2
Kinematics of the pipe.
In the curvilinear base (O, X, y, Z) associated the transverse section of [§2.1] one notes displacements
3
2
1
,
U
U
U
and
where:
sin
)
,
(
cos
)
,
(
sin
)
,
(
sin
)
(
)
(
)
,
,
(
cos
)
,
(
sin
)
,
(
cos
)
,
(
cos
)
(
)
(
)
,
,
(
)
,
(
)
,
(
sin
)
(
cos
)
(
)
(
)
,
,
(
3
2
1
X
X
W
X
v
R
X
X
U
X
R
U
X
X
W
X
v
R
X
X
U
X
R
U
X
X
U
R
X
R
X
X
U
X
R
U
X
Z
X
y
Z
y
X
-
-
+
-
=
+
-
-
+
=
+
+
+
-
=
.
Once discretized this expression becomes:
U R X
H X
U
H X
U
H X
R
H X
R
H X
R
U X
X
U R X
H X
U
H X
U
H X
K
K
NR
xk
K
y
K
K
yk
K
yk
K
zk
K
xk
K
K
NR
y
K
K
1
1
2
1
(,)
() (. )
() (
. )
() (
. )
cos
() (
. )
cos
()
sin
(,)
(,)
(,)
() (
. )
() (
. )
() (
=
+
-
-
+
+
+
=
+
+
=
=
X X
y X
X y
y y
X y
y y
X
K
K
K
K
K
K
K
K
K
K
X
y X
X X
y X
. )
cos
() (
. )
cos
(,) cos
(,) sin
(,) cos
(,)
()
() (
. )
sin
() (
. )
sin
(,) sin
(,) cos
(,) sin
xk
K
yk
K
K
NR
zk
K
xk
K
yk
R
H X
R
v X
W X
X
U R X
H X U
H X
R
H X
R
v X
W X
X
+
-
-
+
=
-
-
+
-
-
=
3
1
where
() () () ()
()
,
,
,
,
,
,
,
,
X
X
X
W
X
v
X
U
and
are discretized as with [§3.1].
Displacement in a node
K
of X-coordinate
K
X
end of pipe is written then:
U R X
U
R
R
U
m
U
m
X
U R X
U
R
v
m
v
m
W
m
W
m
W
W
W
K
xk
yk
zk
km
I
m
M
km
O
K
K
y
K
xk
km
I
m
M
km
O
km
I
m
M
km
O
K
I
KB
KB
1
2
2
2
2
1
1
2
2
(,
,)
cos
sin
(
cos
sin
)
(,)
(,
,)
cos
cos
(
sin
cos
) sin
(
cos
sin
)
sin
cos
sin
(
=
-
+
+
+
+
=
+
-
+
-
+
-
+
-
+
=
=
=
X
U R X
U
R
v
m
v
m
W
m
W
m
W
W
W
X
K
K
zk
xk
km
I
m
M
km
O
km
I
m
M
km
O
K
I
KB
KB
K
,) cos
(,
,)
sin
sin
(
sin
cos
) cos
(
cos
sin
)
cos
sin
cos
(,) sin
3
2
2
1
1
2
2
=
-
+
+
-
+
-
-
-
-
=
=
For the pipe the kinetic vector moment
)
(M
U
GM
and swelling
)
(
.
M
U
GM
GM
have for
respective expressions:
-
-
+
+
-
=
sin
)
,
,
(
cos
)
,
,
(
)
,
(
)
(
cos
)
(
sin
)
(
)
(
1
1
2
X
R
Ru
X
R
Ru
X
rv
X
R
R
X
U
R
X
U
M
X
y
Z
U
GM
,
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
43/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
and:
)
,
(
sin
)
(
cos
)
(
)
(
.
X
W
X
U
X
U
M
y
Z
+
-
-
=
U
GM
GM
.
Note:
The first component of the field of displacement utilizes
()
,
X
U
in an isolated way. It
in goes in the same way for the first component of the vector rotation with respect to
()
,
X
v
and of
swelling with respect to
()
,
X
W
. This remark will be used with [§5.6] to bind the modes of
Fourier with the degrees of freedom of the edges of hull.
7.3
Kinematics of hull
The kinematics of hull of Coils-Kirchhoff or of Naghdi-Mindlin is written in the thickness:
()
()
()
(
)
3
.y
Q
Q
M
N
U
U
C
C
C
+
=
·
)
(
Q
C
U
constitute the vector displacement of average surface in
Q
,
·
)
(
Q
C
constitute the vector rotation in
Q
normal according to directions' T
1
and T
2
plan
tangent with
Q
.
E
1
G
X
3
X
1
H
y
2
Q
H
Q
M
E
2
E
3
y
3
T
1
=e
1
y
3
N
Section of hull:
S L I
= ×
L
: line of the points
Q
on the average layer
I
H H
= -
2 2
,
.
interval describing the thickness
This displacement and this rotation are calculated in the total reference mark. It is possible by change
of reference mark to have their expressions in the curvilinear base (O, X, y, Z) of [§2.1] associated the section
transverse of the junction enters the hull and the pipe.
For each node, the program calculates the coefficients of the 9+6 (M-1) linear relations which connect:
·
6 ddl of the node of beam P of the pipe,
·
2+3x2 (M-1) ddl of Fourier of the pipe,
·
the ddl of swelling of the pipe,
·
with the ddl of all the nodes of the list of the meshs of the edge of hull.
These linear relations will be dualisées, like all the linear relations resulting for example from
key word
LIAISON_DDL
of
AFFE_CHAR_MECA
. They are built as for the connection 3D-beam has
to leave the assembly of elementary terms.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
44/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
7.4
Calculation of average displacement on the section S
It is a question of calculating the integral
S
dS
C
U
, where
C
U
of hull (comprising 6 ddl is the displacement by
node),
S
is the edge of hull.
Average displacement on the section
S
is written:
()
()
()
(
)
+
=
-
S
L
L
H
H
DLL
Dy
y
Q
DLL
Q
H
dS
M
2
/
2
/
3
3
N
U
U
C
C
C
that is to say
()
()
=
S
L
DLL
Q
H
dS
M
C
C
U
U
.
In addition one also has for the pipe part:
()
()
()
=
=
+
=
kz
K
y
kx
S
K
Z
K
y
K
X
S
S
S
p
U
U
U
S
dS
U
U
U
dS
M
M
dS
M
]
[
U
U
U
C
.
One establishes as well as the average displacement of the section of the pipe to the node
K
is equal to
displacement of beam of the node
K
. One can thus linearly bind the degrees of freedom of beam of
translation with the node
K
with the average of the degrees of freedom of displacement of the edge of the hull.
One neglects in this expression the variation of metric in the thickness of the hull.
7.5
Calculation of the average rotation of the section S
()
()
(
)
()
() ()
(
)
()
() ()
(
)
()
()
()
() ()
(
)
-
-
-
-
+
+
+
=
+
+
=
L
L
H
H
H
H
L
L
H
H
S
L
H
H
DLL
Dy
y
Q
Q
Q
DLL
Dy
y
Q
Q
Dy
y
DLL
Q
Q
DLL
Q
H
Dy
DLL
y
Q
Q
Q
Q
y
dS
M
2
2
3
2
3
2
/
2
/
3
3
2
/
2
/
3
3
2
/
2
/
3
3
3
.
N
N
U
N
N
GQ
U
GQ
N
U
N
GQ
U
GM
C
C
C
C
C
C
C
that is to say
()
()
()
DLL
Q
H
DLL
Q
H
dS
M
L
S
L
+
=
C
C
C
U
GQ
U
GM
12
.
3
.
In addition one also has for the pipe part:
()
()
()
GM U
GM
U
U
C
=
+
=
=
M dS
M
M ds
R
R
R
dS
I
p
S
S
S
X
K
y
K
Z
K
S
X
K
y
K
Z
K
[
]
cos
sin
2
2
2
2
2
.
where I is the tensor of inertia of the beam. One establishes as well as the average rotation of the section of the pipe
with the node
K
is equal to the rotation of beam to the node
K
. One can thus bind the degrees linearly of
freedom of beam of rotation to the node
K
with the degrees of freedom of rotation of the edge of the hull.
One neglects in this expression the variation of metric in the thickness of the hull.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
45/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
7.6
Calculation of the average swelling of the section S
It is a question of calculating the integral
=
S
S
dS
dS
GM
GM
U
N
U
C
C
.
.
, where
GQ
GQ
GM
GM
N
=
=
is the normal with
surface average hull.
Average displacement on the section
S
is written:
()
()
()
(
)
()
=
+
=
-
L
S
L
L
H
H
DLL
Q
H
DLL
Dy
y
Q
DLL
Q
H
dS
M
N
U
N
N
N
U
N
U
C
C
C
C
.
.
.
.
2
/
2
/
3
3
.
In addition one also has for the pipe part:
()
()
()
dS
W
dS
M
M
dS
M
S
ok
S
S
S
p
=
+
=
]
.[
.
U
U
GM
GM
U
GM
GM
C
One establishes as well as the average swelling of the section of the pipe to the node
K
is equal to the degree of
freedom of swelling of the pipe to the node
K
. One can thus linearly bind the degree of freedom of
swelling of pipe to the node
K
with the degrees of freedom of displacement of the edge of the hull.
One neglects in this expression the variation of metric in the thickness of the hull.
7.7
Calculation of the coefficients of Fourier on the section S
It is a question of calculating the six integrals
S
dS
p
cos
C
U
,
S
dS
p
sin
C
U
,
S
dS
p
cos
GM
GM
U
C
,
S
dS
p
sin
GM
GM
U
C
,
S
dS
p
cos
.
C
U
GM
GM
and
S
dS
p
sin
.
C
U
GM
GM
, where
C
U
is it
displacement of hull (comprising 6 ddl by node),
S
is the edge of hull.
One with the following relation for displacements on the section
S
:
()
()
()
(
)
+
=
-
S
L
L
H
H
DLL
p
Dy
y
Q
DLL
p
Q
H
dS
p
M
cos
cos
cos
2
/
2
/
3
3
N
U
U
C
C
C
that is to say
()
()
=
S
L
DLL
p
Q
H
dS
p
M
cos
cos
C
C
U
U
and
()
()
=
S
L
DLL
p
Q
H
dS
p
M
sin
sin
C
C
U
U
.
In addition one also has for the pipe part:
()
dS
p
X
R
U
X
R
U
X
R
U
dS
p
M
S
K
K
K
S
cos
)
,
,
(
)
,
,
(
)
,
,
(
cos
3
2
1
=
C
U
.
The first component of this relation then enables us to connect the coefficient linearly of
Fourier
ikp
U
with the components of displacements of the edge of hull in the following way:
-
=
=
1
p
cos
1
p
cos
cos
)
,
,
(
cos
)
(
2
2
1
1
if
if
dS
p
U
=
dS
R
dS
p
X
R
U
DLL
p
Q
U
H
S
ikp
S
ky
S
K
L
C
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
46/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
In the same way
()
dS
p
X
R
U
X
R
U
X
R
U
dS
p
M
S
K
K
K
S
sin
)
,
,
(
)
,
,
(
)
,
,
(
sin
3
2
1
=
C
U
from where one deduces from it that:
=
=
1
p
sin
1
p
sin
sin
)
,
,
(
sin
)
(
2
2
1
1
if
if
dS
p
U
=
dS
R
dS
p
X
R
U
DLL
p
Q
U
H
S
O
kp
S
kz
S
K
L
C
One with the following relation for rotations on the section
S
:
()
()
()
DLL
p
Dy
y
Q
DLL
p
Q
H
dS
p
M
L
H
H
S
L
+
-
+
=
cos
)
(
)
(
cos
cos
.
3
2
/
2
/
3
N
N
GQ
GQ
U
GM
GM
U
C
C
C
.
The first component of this relation then enables us to connect the coefficient linearly of
Fourier
O
kp
v
with the components of displacements and rotations of the edge of hull in the manner
following:
()
+
-
=
1
p
)
(
cos
1
p
]
cos
cos
[
cos
]
[
2
2
1
2
1
if
if
S
O
kp
S
O
K
K
y
L
dS
p
rv
=
dS
rw
Ru
DLL
p
Q
H
GQ
GQ
U
C
In the same way one a:
()
+
=
1
p
)
(
sin
1
p
]
sin
sin
[
sin
]
[
2
2
1
2
1
if
if
S
ikp
S
ik
kz
L
dS
p
rv
=
dS
rw
Ru
DLL
p
Q
H
GQ
GQ
U
C
One with the following relation for swelling on the section
S
:
()
()
=
S
L
DLL
p
Q
H
dS
p
M
cos
.
cos
.
.
C
C
U
GQ
GQ
U
GM
GM
.
This relation enables us to connect the coefficient of Fourier linearly
ikp
W
with the components of
displacements of the edge of hull in the following way:
()
+
-
=
1
p
)
(
cos
1
p
]
cos
cos
[
cos
.
2
2
1
2
if
if
S
ikp
S
ik
kz
L
dS
p
W
=
dS
W
U
DLL
p
Q
H
C
U
GQ
GQ
.
In the same way, one a:
()
+
-
=
1
p
)
(
sin
1
p
sin
sin
sin
.
2
2
1
2
if
if
S
O
kp
S
O
K
K
y
L
dS
p
W
=
dS
W
U
DLL
p
Q
H
C
U
GQ
GQ
One uses for all these relations the fact that
=
S
dS
Q
p
.
if
Q
p
0
cos
cos
One neglects in this expression the variation of metric in the thickness of the hull.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
47/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Note:
One can note that some of the relations established in this paragraph for p=1 are
redundant with those established in the paragraphs [§7.4] and [§7.5]. On the six relations
established starting from the calculation of the integral forms
U
C
cos
dS
S
,
U
C
sin
dS
S
,
U
GM
GM
C
cos
dS
S
,
U
GM
GM
C
sin
dS
S
,
GM
GM U
C
.
cos
dS
S
and
GM
GM U
C
.
sin
dS
S
, only
two among the four last are linearly independent of the others. Thus both
first were already established with [the §7.4] and of the combinations of the four last give again
those of [§7.5].
7.8
Establishment of the method
The calculation of the coefficients of the linear relations is done in three times:
·
calculation of elementary quantities on the elements of the list of the meshs of edges of hulls
(mesh of type SEG2):
- surface
=
elt
elt
elt
elt
Z
y
X
;
;
;
1
.
-
summation of these quantities on (
S
) from where the calculation of:
-
S
With
=
- position
of
G
.
·
knowing
G
, elementary calculation on the elements of the list of the meshs of edges of
hulls of:
{
}
element
of
form
of
functions
:
where
=
=
Ni
Z
y
X
zNi
yNi
xNi
Ni
elt
elt
elt
elt
,
,
;
;
;
GM
It should be noticed that in the case of the connection hull - pipe, the integrals on the elements of edge
are to be multiplied by the thickness of the hull:
=
L
I
elt
I
NR
H
NR
where
L
represent average fiber of
the element of edge of hull. Moreover, one adds the additional term:
L
I
NR
H
12
3
.
·
“assembly” of the terms calculated above to obtain of each node of the section
of connection, coefficients of the terms of the linear relations,
·
connection between the modes of Fourier and displacements of hull as shown at the beginning
[§7].
More precisely:
·
for the connection hull - pipe, one carries out elementary calculations on all the elements of
edge of the section of connection
S
type:
()
()
()
=
=
=
=
S
elements
Nb
N
Z
y
X
R
X
N
N
D
U
U
U
m
D
m
U
U
U
D
m
1
2
0
2
0
2
1
cos
1
cos
1
cos
1
C
C
C
C
C
C
C
cm
P
U
U
and
()
()
()
=
=
=
=
S
elements
Nb
N
Z
y
X
R
X
N
N
D
U
U
U
m
D
m
U
U
U
D
m
1
2
0
2
0
2
1
sin
1
sin
1
sin
1
C
C
C
C
C
C
C
Sm
P
U
U
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
48/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
if
1
>
m
where
P
is the matrix of passage of the local reference mark of the element to the total reference mark and
1
position orthoradiale of the element. By expressing displacement according to the degrees of freedom
nodal:
()
()
=
=
=
=
S
elements
Nb
N
nodes
Nb
N
N
Z
N
y
N
X
N
I
m
O
m
I
m
N
N
D
U
U
U
NR
m
W
v
U
1
_
1
2
1
cos
1
P
U
cm
and
()
()
=
=
=
=
S
elements
Nb
N
nodes
Nb
N
N
Z
N
y
N
X
N
O
m
I
m
O
m
N
N
D
U
U
U
NR
m
W
v
U
1
_
1
2
1
sin
1
P
U
Sm
where them
NR
are related to form of the element, one obtains for each calculation twice
9 coefficients with the nodes of the element running of
S
:
() () ()
() () ()
=
=
=
=
=
L
N
ij
N
ij
N
ij
N
Z
N
y
N
X
S
elements
nodes
Nb
N
I
m
O
m
I
m
dx
X
NR
X
P
MX
R
D
NR
P
m
has
U
U
U
has
has
has
has
has
has
has
has
has
W
v
U
0
_
1
33
32
31
23
22
21
13
12
11
cos
1
cos
1
.
2
1
cm
U
and an equivalent expression for
Sm
U
where
I
is the length of the element of edge of hull.
·
for the connection 3D - pipe, one carries out elementary calculations on all the elements of
edge of the section of connection
S
type:
()
()
()
=
=
=
=
=
S
elements
Nb
N
N
D
Z
D
y
D
X
nodes
Nb
N
N
H
H
S
R
X
S
N
N
N
N
D
rd
U
U
U
NR
m
S
D
rd
m
U
U
U
S
dS
m
S
1
3
3
3
_
1
2
1
2
1
)
,
(
cos
2
cos
2
cos
2
P
U
U
3D
3D
3D
3D
cm
and
()
()
()
=
=
=
=
=
S
elements
Nb
N
N
D
Z
D
y
D
X
nodes
Nb
N
N
H
H
S
R
X
S
N
N
N
N
D
rd
U
U
U
NR
m
S
D
rd
m
U
U
U
S
dS
m
C
S
1
3
3
3
_
1
2
1
2
1
)
,
(
sin
2
sin
2
sin
2
P
U
U
3D
3D
3D
3D
Sm
if
1
>
m
where
P
is the matrix of passage of the local reference mark of the element to the total reference mark,
1
the position
orthoradiale of the element,
1
H
its radial position and them
NR
are related to form of the element.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
49/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
8
Establishment of the element PIPE in Code_Aster
8.1 Description
This new element (of name
METUSEG3
) is pressed on a mesh
SEG3
or
SEG4
curvilinear. It supposes
that the section of the pipe is circular. Contrary to the elements
POU_D_E
,
POU_D_T
, [R3.08.01] this
element is not <<exact>> with the nodes for loadings or torques concentrated with
ends. It is thus necessary to net with several elements to obtain correct results.
8.2
Introduced use and developments
The element is used in the following way:
AFFE_MODELE (MODELING = “TUYAU_3M'…)
The meshs with 4 nodes are generated starting from the meshs with 3 nodes using:
MALL =CREA_MAILLAGE (MAILLAGE=MAIL,
MODI_MAILLE=_F (OPTION = “SEG3_4”, ALL = “YES”)
)
One calls upon the routine
INI090
for the functions of form, their derivative and their derivative
seconds (for the hull part) at the points of Gauss, as well as the corresponding weights.
The characteristics of the section are defined in
AFFE_CARA_ELEM
AFFE_CARA_ELEM (BEAM = _F (SECTION = “CIRCLE”,
CARA = (“R” “EP”),
VALE = (.......),),
ORIENTATION=_F (GROUP_NO=D, CARA=' GENE_TUYAU', VALE= (X Y Z),),
TUYAU_NCOU = ' A NUMBER OF COUCHES', TUYAU_NSEC = ' A NUMBER OF SECTEURS',),
)
R
and
EP
represent, as for the elements of beams conventional, respectively the radius
external and the thickness of the section. One also defines on one of the nodes end of the line of
piping the vector whose projection on the transverse section is the origin of the angles for
decomposition in Fourier series. This vector should not be colinéaire with the average line of the elbow
with the node end considered.
One also defines in this level the number of layers and angular sectors to use for integration
numerical.
AFFE_CHAR_MECA (DDL_IMPO = _F (
DX =., DY =., DZ =., DRX =., DRY =., DRZ =.,
DDL of beam
UI2 =., VI2 =., WI2 =., UO2 =., VO2 =., WO2 =.,
DDL related to mode 2
UI3 =., VI3 =., WI3 =., UO3 =., VO3 =., WO3 =.,
DDL related to mode 3
WO =., WI1 =., WO1 =.,
DDL of swelling and mode 1 on W
FORCE_NODALE = _F (FX =., FY =., FZ =., MX =., MY =., MZ =. )
They are the conventional forces of beam, which work only on displacements of beam.
Addition of a key word in
AFFE_CHAR_MECA (FORCE_TUYAU = _F (CLOSE =. ))
for the calculation of
work of the internal pressure.
The pressure works on the DDL of swelling
WO
, one calculates then:
O
K
L
NR
K
K
L
NR
K
O
K
K
L
O
near
W
dx
D
R
p
H
dx
D
R
W
H
p
dx
D
R
W
p
W
=
=
=
=
=
0
2
0
1
0
2
0
1
0
2
0
int
int
int
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
50/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
8.3
Calculation in linear elasticity
The matrix of rigidity and the matrix of mass (respectively options
RIGI_MECA
and
MASS_MECA
)
are integrated numerically in
TE0582
. Calculation holds account owing to the fact that the terms
corresponding to the DDL of beam are expressed classically in total reference mark, and that DDL of
Fourier are in the local reference mark with the element. In the case or the element does not belong to any elbow, it
identify local is defined by the generator and the directing vector carried by average fiber of the element
as indicated on [Figure 8.3-a]. If the element belongs to an elbow, the local reference mark
is defined starting from the plan of the elbow as mentioned with [§2.1].
Z
X
Y
O
Z
X
Z
y
generator
Appear 8.3-a: Identifies local for a right pipe
Elementary calculations (
CALC_ELEM
) currently available correspond to the options:
·
EPSI_ELGA_DEPL
and
SIEF_ELGA_DEPL
who provide the strains and the stresses
at the points of integration in the local reference mark of the element. Calculation is carried out in
TE0584
,
and currently gives the values to the 693 points of integration (for an element to 3 modes of
Fourier). These fields are called fields at “under-points” of integration. One stores these
values in the following way:
-
for each point of Gauss in the length, (n=1, 3)
-
for each point of integration in the thickness, (n=1, 2N
NECK
+1=7)
-
for each point of integration on the circumference, (n=1, 2N
SECT
+1=33)
-
6 components of strain or stresses:
EPXX
EPYY
EPZZ
EPXY
EPXZ
EPYZ
or
SIXX
SIYY
SIZZ
SIXY
SIXZ
SIYZ
where X indicates the direction given by the two nodes nodes of
the element, Y represents the angle
describing the circumference and Z represents
the radius.
EPZZ
and
EPYZ
agent with
R
rr
,
in the case of them
deformations and
SIZZ
and
SIYZ
agent with
R
rr
,
in the case of them
stresses are taken equal to zero.
·
EFGE_ELNO_DEPL
: who gives the efforts generalize of beam conventional:
NR, VY, VZ,
MT, MFY, MFZ
. These efforts are given in the local curvilinear reference mark of the element. This
option is calculated in
TE0585.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
51/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
·
SIGM_ELNO_TUYO
, and
VARI_ELNO_TUYO
allow, starting from the fields
SIEF_ELGA
or
VARI_ELGA
, and of the data of an angular position and a position in the thickness, of
to calculate the stresses or variables intern with the nodes of the elements, with this position. One
then obtains a field without under-points, i.e. a conventional field, usable for
operations of postprocessing (for example the statement of a component of stresses length
one line of piping).
·
Options
EQUI_ELGA_SIGM
and
EQUI_ELGA_EPSI
allow the calculation of the invariants, in
each point of integration (fields at “under-points”)
·
The option
VALE_NCOU_MAXI
allows to extract, of each of the 3 linear points of Gauss
of an element, the values maximum and minimum of a component of a field.
·
In version 7, one can extract from other fields in a point of the section:
EPSI_ELNO_TUYO
for the deformations,
SIEQ_ELNO_TUYO
for the equivalent stresses,
EPEQ_ELNO_TUYO
for the equivalent deformations.
Finally it
TE0585
calculate also the option
FORC_NODA
for the operator
CALC_NO
.
8.4
Plastic design
The matrix of tangent rigidity (options
RIGI_MECA_TANG
and
FULL_MECA
) as well as projection
plastic (options
FULL_MECA
and
RAPH_MECA
) are integrated numerically in
TE0586
. One makes
call to the option of calculation
STAT_NON_LINE
. All laws of plane stresses available in
Code_Aster can be used: if they are not integrated directly, it is always possible
to use a law of behavior formulated in plane deformation, and to treat the assumption of
plane stresses using the method of Borst:
STAT_NON_LINE (…
COMP_INCR = _F (RELATION = ' “, ALGO_C_PLAN=' DEBORST”),
.....)
Elementary calculations (
CALC_ELEM
) currently available correspond to the options:
·
SIEF_ELNO_ELGA
who allows to obtain the efforts generalized by element with the nodes in
identify beam. This option is calculated in
TE0587.
·
VARI_ELNO_ELGA
who calculates the field of internal variables by element with the nodes for
all layers and all sectors, in the local reference mark of the element. This option is
calculated in
TE0587.
8.5
Test: SSLL106A
It is about a right pipe of directing vector (4, 3, 0) fixed in its end O and which is with a grid with 18
elements
PIPE
.
O
B
L
L = 5
Y
X
Z O
4
3
X
y
Z
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
52/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
The pipe is subjected to different types of load:
·
a tensile load,
·
2 sharp efforts,
·
2 bending moments,
·
1 torque,
·
an internal pressure.
Displacements at the point B, the strain and the stresses in certain points are calculated
of integration of the section containing B, as well as the first clean modes.
This makes it possible to test the DDL of beam, the DDL of swelling and modes 1 of the development in
Fourier series.
9 Conclusion
The finite elements of elbow which we describe here are usable for calculations of piping in
elasticity or in plasticity. Pipings can be subjected to various combined loadings -
internal pressure, cross-bendings and anti-plane, torsion, extension.
For the moment, the element carried out is a linear element of beam type, right or curve, to three
nodes, in small rotations and deformations, with a local elastoplastic behavior in
plane stresses. It makes it possible to take into account ovalization, the roll and swelling. It
combine the properties of hulls and beams. The kinematics of beam for the axis of the elbow is
increased by a kinematics of hull, of type Coils-Kirchhoff without transverse shearing, for
description of the behavior of the transverse sections. This last kinematics is discretized in
M modes of Fourier, of which the number M, which the literature encourages us to choose equal to 6 [bib8], [bib13],
must at the same time be sufficient to obtain good results in plasticity and not too large to limit it
calculating time. In elasticity, for relatively thick pipings (the thickness report/ratio on
radius of the transverse section higher than 0.1), one can be satisfied with M=2 or 3.
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
53/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
10 Bibliography
[1]
ABAQUS: “Theory manual” Hibbit, Karlsson and Sorensen Inc. (1984) chapter 4.4.1
[2]
K.J. BATHE and A.C. ALMEIDA: “A simple and effective pipe elbow element-linear analysis”,
Newspaper off Applied Mechanics, Vol.47 (1980) pp93-100.
[3]
J.L. BATOZ, G. DHATT: “Modeling of the structures by finite elements - Hulls”, Vol.3,
Hermès, Paris, 1992.
[4]
W.B. BICKFORD, B.T. STROM: “Vibration off planes curved beams”, Journal off Sound and
Vibration (1975) 39 (2), pp.135-146.
[5]
J.T. BOYLE and J. SPENCE: “Inelastic analysis methods for piping systems”, Nuclear
Engineering and Design 57 (1980) 369-390.
[6]
J.T. BOYLE and J. SPENCE: “A state off the art review off inelastic (static & dynamic) piping
analysis methods, with particular application to LMFBR ", Glasgow, Scotland, the U.K. (1983).
[7]
P. GUIHOT: “Accidental mechanical Behavior of the Primary education Circuit”, Notes EDF
HP-52/96/012/A.
[8]
H.D. HIBBIT: “Special structural elements for piping analysis”, Proc. Conf., Presses Vessels
and Piping: Analysis and Computers, Miami, ASME, 1974.
[9]
H.D. HIBBIT & E.K. LEUNG: “Year approach to detailed inelastic analysis off thin-walled
pipelines ", Special ASME Pub., “Nonlinear Finite Element Analysis off Shells” (1981).
[10]
A. KANARACHOS, R.N. KOUTSIDES: “A new Approach off Shell displacements in A
Beam-type pipe element “, Proc. 8th Int. Conf. SMIRT, paper B9/2, Brussels (1985).
[11]
P. MASSIN, A. BEN HAJ YEDDER: “Bibliography for the development of an element
bend with internal pressure and local plasticity ", Note HI-74/97/026, EDF-DER 1998.
[12]
A. MILLARD: “Year enriched beam finite element for accurate piping analysis”, Proc. 14th Int.
Conf. SMIRT, Lyon (1997).
[13]
H. OHTSUBO, O. WATANABÉ: “Stress analysis off pipe bends by boxing ring elements”, Trans.
ASME, Journal off Presses Vessel technology, Vol.100 (1978) 112-122.
[14]
K. WASHIZU: “Variational Methods in Elasticity & Plasticity”, 3rd ED. Pergamon (1981).
[15]
J. PELLET: “Connection 3D-beam”, Reference material of Code_Aster [R3.03.03].
[16]
J.M. PROIX
: “Connection hull-beam”, Reference material of Code_Aster
[R3.03.06].
[17]
F. VOLDOIRE, C. SEVIN: “Axisymmetric thermoelastic Hulls and 1D”, Documentation
of Reference of Code_Aster [R3.07.01].
[18]
J.M. PROIX, P. MIALON, m.t. BOURDEIX: “Exact Elements of right and curved beams”,
Reference material of Code_Aster [R3.08.01].
Code_Aster
®
Version
6.4
Titrate:
Finite elements of right pipe and curve
Date:
12/12/03
Author (S):
P. MASSIN, J.M. PROIX, A. BEN HAJ YEDDER
Key
:
R3.08.06-B
Page
:
54/54
Manual of Reference
R3.08 booklet: Machine elements with average fiber
HT-66/03/005/A
Intentionally white left page.