Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
1/10
Organization (S): EDF/IMA/MN
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
Document: R3.06.03
Calculation of the constraints to the nodes by local smoothing
Summary:
One presents a local method of calculation of constraints at the nodes starting from the constraints at the points of
GAUSS. It is used in options SIGM_ELNO_DEPL and SIEF_ELNO_ELGA of the command
CALC_ELEM [U4.61.02].
This method is summarized to calculate the constraints at the tops of an element by multiplying the constraints
at the points of GAUSS by a matrix of smoothing, constant for each type of element.
For the isoparametric elements of degree 2, the constraints with the nodes mediums are obtained by
average of the values of the constraints at the 2 tops of the edge.
This method of smoothing has two advantages:
· the nodal constraints obtained have a command of precision moreover than by direct calculation with the nodes,
· the method is inexpensive in time CPU.
This method was generalized:
· with calculations of the deformations (option EPSI_ELNO_DEPL) and variables intern (option
VARI_ELNO_ELGA) with the nodes in mechanics,
· with the calculation of flows (option FLUX_ELNO_TEMP) to the nodes in thermics.
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
2/10
1 Preliminaries
This method is based on the observation [bib1] which it exists points where the calculation of the constraints, with
to leave displacements in a primal formulation in displacements, is more precise.
In the case of isoparametric finite elements of command 2 (SEG3 in 1D, QUAD8 and QUAD9 in 2D, HEXA20
in 3D), one shows that points of GAUSS of the formula of quadrature to 2n points (N: dimension of
space) are such as one can hope, without that being formally shown, for the calculation of
the same command of precision as for the calculation of the field of displacement U.
The idea of the method is to calculate for each element the constraints! with the nodes from
K at the points of GAUSS, these last being calculated on each element by the formula:
NR 0
K
K
K
= dB U = D iB Ui
i=1
where:
D is the matrix of elasticity,
Bk is the matrix connecting the deformations to displacements at the point of GAUSS K,
U are nodal displacements.
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
3/10
2
Local method of minimization by least squares
Generally, one wishes to approximate, within the meaning of least squares, a function (X, y)
by a polynomial function:
! (X, y) = xi yi has
ij
i=0, p
i=0, Q
The problem amounts finding the coefficients aij which minimize the functional calculus:
= (
2
-
!) dx Dy
The values of the function are known here only at the points of Gauss: K = (xk, yk)
The minimum will be reached if and only if:
I = 0,…, p
=
0
has
J = 0,…, Q
ij
Within the framework of the finite element method in displacement, the function of smoothing is written:
N
! (X, y) = NR (X, y
I
)!I
i=1
where:
Ni is related to form associated with node I,
!
I is the value of the sought constraint to node I,
N the number of nodes retained for smoothing.
One must thus solve the system:
= 0 I =,
1…, N
éq 2-1
!I
One can choose between two methods of local smoothing: continuous smoothing or discrete smoothing.
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
4/10
3
Methods of local smoothing (ref. [bib2] and [bib3])
3.1
Continuous local smoothing
2
2
This type of smoothing led to solve the system [éq 2-1] with = (-!) =
NR
.
E
(- I!I)
E
What leads to M E
Fe
! =
with:
npg
M E = NR NR dx Dy = NR
ij
I J
() Nj (K) (det J)
E
I
K
K
K
K =1
npg
and
F E = NR dx Dy =
NR
I
I
I (K
) K (det J)
K
K
E
K =1
where
K are the points of GAUSS in the element of reference
(det J)
the jacobien of the geometrical transformation enters the element of reference and
K
the element running at the point K.
K: the weight associated with the point K
K: the constraint at the point K
Nor (K
): the value of the function of form in the element of reference to the point K
^
^
^
3
2
^
1
^
4
3
3
2
2
4
4
1
1
direct calculation of the constraints
smoothed constraints
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
5/10
Note:
· If spaces of interpolation of and of! are the same ones, one has obviously =! .
In practice, one retains for space of! a space smaller than that where is defined
by the finite element.
· One sees the link between the approximation at the points of GAUSS of where thus converges
better and this process of smoothing whose justification is on the contrary continuous.
· The way in which is calculated at the points of GAUSS does not intervene. Generalization
with the nonlinear problems is thus obvious, although it cannot concern same
justification.
This method is however not adopted because it requires a resolution of system linear for
each calculation of!
.
3.2
Discrete local smoothing
In this case, the functional calculus is replaced by the summation:
npg
npg
N
2
~
2
= ((K) -!(K) =
(K) - Ni!I (K)
K =1
K =1
i=1
^
^
3
3
^
3
III
4
^
2
4
IV
II
2
4
2
^
1
I
1
1
constraints at the points of GAUSS
smoothed constraints
~
The system to be solved is written there still:
=
0 are:
!I
npg
npg
NR () NR ()! = NR () () I
I
K
J
K
J
I
K
K
K =1 J
K =1
maybe in matric form:
M {!node} = P {GAUSS}
The matrices M and P being then independent of the element running E.
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
6/10
They can thus be calculated once and for all on the element of reference.
Note:
· This method is more economic than the preceding one and gives results
comparable [bib2],
· There still, the manner thus K is calculated in each point of GAUSS is indifferent
(since the number of points of GAUSS used for the calculation of and that of!
is the same one). One will be able to thus use this method into nonlinear (option
SIEF_ELNO_ELGA).
4
Application of the method to calculation of the constraints to
nodes for various elements
The local smoothing adopted in Code_Aster is the discrete local smoothing [§2.2], which makes it possible to avoid it
calculation of integrals on the element.
On all the elements of continuous medium 2D and 3D, one chose a space of smoothing resting on
functions of form relating to the nodes of the element.
The method thus makes it possible to obtain the constraints at the tops. In the case of elements of command 2,
one calculates the constraints with the nodes mediums by taking the average value of the two nodes
“framing” the node medium considered.
One gives hereafter the matrices of passage allowing to calculate the constraints with the nodes
nodes starting from the constraints at the points of GAUSS. These matrices can be square or
rectangular. Indeed, matrices of passage Mr. P
1 is calculated once and for all with
the initialization of each type of finite element (in AFFE_MODELE). Two types of matrices exist:
·
matrices
Square M-1 P, which are to be used when the number of points of GAUSS used
for the calculation of the constraints at the points of GAUSS K is identical to the number of nodes
nodes,
·
matrices
Rectangular M-1 P, which are to be used when the number of points of GAUSS
K is different (in general higher) nodes nodes.
4.1
Square matrices of passage
These matrices are used by all the elements for option SIGM_ELNO_DEPL. The option calculates in
first constraints in a number of points of GAUSS equal to the number of nodes. Then them
matrices M-1 P (given afterwards) are used to calculate the constraints with the nodes. These matrices
are also used for option SIEF_ELNO_ELGA, in the elements for which the number
points of GAUSS of the calculation of SIEF_ELGA (in STAT_NON_LINE) is equal to the number of
nodes. They are the elements:
· in 2D: QUAD4, TRIA6, under-integrated QUAD8,
· in 3D: TETRA4, PENTA6, HEXA8, PYRAM5 and under-integrated HEXA20.
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
7/10
4.1.1 Square matrices of passage for the elements 2D
4.1.1.1 Triangles
5 - 1 -
1
1
1
Mr. P =
- 1
5 -
1
3 - 1 - 1 5
4.1.1.2 Quadrangles
3
1
3
1
1+
-
1
-
2
2
2
2
1
3
1
3
-
1+
-
1
1
2
2
2
2
Mr. P =
3
1
3
1
1
-
1+
-
2
2
2
2
1
3
1
3
-
1
-
1+
2
2
2
2
4.1.2 Square matrices of passage for the elements 3D
4.1.2.1 Tetrahedrons
has
has
- 1 has
has
has has -
1
has
has
Mr. P
1
1
=
- B has - 1 has
has
has
has
has
has
has
has -
1
5 - 5
5 + 3 5
has =
B =
20
20
4.1.2.2 Pentahedrons
-
1 - - 1 1 -
-
1 - 1 - - 1
-
- 1 1 - 1 -
M-1 P =
1 - - 1 1 -
-
1 - 1 - - 1
-
- 1 1 - 1 -
-
3 +1
= 2
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
8/10
4.1.2.3 Hexahedrons
B B C B C C D has
5 + 3 3
has =
4
B C C D has B B C
(1+ 3)
C D B C B a.c. B
B =
4
B a.c. B C D B C
Mr. P
1
=
B has C B C B D C
(3 -) 1
C =
C B D C B has C B
4
D C C B C B B has
(5-3 3)
C B B has D C C B
D =
4
7
6
VIII
VI
8
5
VII
IV
3
V
II
Z
2
III
y
X
I
4
1
Appear 4.1.2.3-a: Numérotation of the points of GAUSS
on the hexahedron with 8 nodes
4.2
Matrices of passage Mr. P
1 rectangular
Into nonlinear for certain types of elements (TRIA3, QUAD8 and QUAD9 in 2D, TETRA10, PENTA15 and
HEXA20 in 3D), the internal constraints and variables at the points of GAUSS are calculated on one
family of points of richer GAUSS (9 points for the quadrangles, 15 points for the tetrahedrons, 21
points for the pentahedrons, 27 points for the hexahedrons).
Discrete local smoothing is then carried out starting from these fields and made transport with the nodes
to intervene of the matrices different from the preceding ones. They are not square any more, because of dimension
(Nb nodes, Nb points of GAUSS). Matrices of passage Mr. P
1 is not calculated
explicitly, in particular M is reversed by Aster.
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
9/10
In the particular case of the triangle with 3 nodes, the fields are supposed to be constant by element (one
only point of GAUSS) and:
1
1
Mr. P =
1
1
For example, the calculation carried out by option SIEF_ELNO_ELGA is then the following:
If the constraints were calculated (in STAT_NON_LINE for example) on a family
having a number of points of GAUSS higher than the number of nodes (for the elements
Nb nodes Nb pts Gauss
announced above). Mr. P
1 is then rectangular, and!
- 1
K
I =
(MR. P).
ik
i=1
K =1
If not, if the number of points of GAUSS equal to the number of nodes, (Mr. P
1) is then square.
One calculates!
-
K
I = (MR. P
1) [§4.1].
ik
5
Other options of calculation using the same method
The method described previously is used in Code_Aster to calculate the deformations, them
internal variables and flows with the nodes.
The list of modelings supporting these options is given below, with the numbers of
routines of elementary terms YOU corresponding.
The produced fields are cham_elem with the nodes.
5.1 Phenomenon
:
“MECANIQUE”
MODELISATION
SIGM_ELNO_DEPL
EPSI_ELNO_DEPL
SIEF_ELNO_ELGA
VARI_ELNO_ELGA
AXIS
TE0086
TE0087
TE0098
AXIS_SI
TE0086
TE0087
TE0098
C_PLAN
TE0086
TE0087
TE0098
D_PLAN
TE0086
TE0087
TE0098
D_PLAN_SI
TE0086
TE0087
TE0098
AXIS_FOURIER
TE0116
TE0114
3D
TE0023
TE0025
TE0020
3d_SI
TE0023
TE0025
TE0020
AXIS_META
TE0352
TE0087
TE0098
3d_META
TE0357
TE0025
TE0020
AXIS_INCO
TE0448
TE0447
not disp
PLAN_INCO
TE0448
TE0447
not disp
3d_INCO
TE0454
TE0453
not disp
COQUE_AXIS
TE0230
TE0229
not disp
COQUE_C_PLAN
TE0230
TE0229
not disp
COQUE_D_PLAN
TE0230
TE0229
not disp
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Code_Aster ®
Version
3
Titrate:
Calculation of the constraints to the nodes by local smoothing
Date:
23/01/97
Author (S):
X. DESROCHES
Key:
R3.06.03-B
Page:
10/10
5.2 Phenomenon
:
“THERMIQUE”
MODELISATION
FLUX_ELNO_TEMP
META_ELNO_TEMP
META_INIT_ELNO
AXIS
TE0069
TE0067
TE0320
PLAN
TE0069
TE0067
TE0320
AXIS_FOURIER
TE0265
not disp
not disp
3D
TE0063
TE0064
TE0321
COQUE
not disp
not disp
not disp
6
Other methods of smoothing of constraints
There are two other methods of smoothing, relating only to the constraints, used by
estimators of Zhu-Zienkiewicz version 1 and 2 [R4.10.01 §3].
The stress fields to the produced nodes are then cham_no.
The corresponding options of calculation are accessible by command CALC_ELEM [U4.61.02].
7 Bibliography
[1]
BARLOW J. - Optimal stress hirings in finite element models - International Journal for
Numerical Methods in Engineering Vol.10 p 243 - 251 (1976).
[2]
HINTON E., CAMPBELL JJ. - Total Room and smoothing off discontinuous finite element
functions using has least public gardens method - International Journal for Numerical Methods in
Engineering Vol.8 p 461 - 480 (1974).
[3]
HINTON E., SCOTT F.C., RICKETTS R.E. - Local least public gardens stress smoothing for
parabolic isoparametric elements - Int. J. for Num. Meth. in Eng. Vol 9 p 235 - 256 (1975)
Handbook of Référence
R3.06 booklet: Machine elements and thermal for the continuous mediums
HI-75/97/004/A
Outline document