Code - Aster: Constraints, Efforts, Forces and Deformations
Code - Aster: Constraints, Efforts, Forces and Deformations
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 1/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
Summary:
This document defines the sizes characterizing the constraints, the forces and the deformations inside a
structure in a calculation by finite elements in displacement and how that is translated in Code_Aster. The
expression of these sizes is given for the finite elements of mechanics: continuous medium 2D or 3D, hulls and
beams.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 2/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
1 Statics
1.1 Constraints
The postulate of Cauchy is that the efforts of contacts exerted in a point by part of a continuous
medium on another depends only on the normal on the surface in this point delimiting the parts.
In accordance with this postulate, one calls vector forced, for the nonmicropolar mediums, F n the
vector which characterizes the forces of contact exerted through an element of surface DS of normal
n on part of a continuous medium [bib1].
It is shown [bib3], then, that the dependence in a point fixed of F compared to the normal n is linear
and that there exists a tensor which one calls tensor of the constraints such as:
F n = n
1.2 Effort
With regard to the structures of beams or hulls, contrary to the case of the continuous medium, it
should be noted that:
• only normal directions n cuts according to tangent space with the variety are possible,
• the characteristic sizes are obtained by integration in the section or the thickness of the sizes
defined for the continuous mediums.
[ ]
F
M
=k .
D
R
One calls effort, the end cells F , M in P , geometrical centre of inertia of the cross-section ,
torque resulting from the forces of contact exerted on the section [bib2].
F =∫ F ds N
Mp =∫ PM∧F ds N . m
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 3/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
P
The force F breaks up into a normal effort N and of the efforts cutting-edges T in the plan of the
section while moment M being exerted at the point P breaks up into one bending moment and
torque at the point P .
For the beams whose cross-section is not regarded as rigid these end cells are not sufficient: for
example, for the beams taking of account the warping of the sections one is brought to consider an
additional size of effort due to warping (bimoment).
The multifibre beams (with local behavior 1D, connecting constraints to deformations, in a certain
number of points of the section) and the pipes (local behavior in plane constraints) provide at the
same time the end cells of the beams but also a stress field for each fibre or sector.
That is to say a point P of a surface medium S of thickness h , that is to say an element length dl
on S , that is to say n the normal directing the hull in this point.
Are the end cells F , M in this point of a torque resulting from the forces exerted through an
element of surface dS =h dl of normal n tangent with S on part of S .
With the preceding notations:
h/ 2
F P = ∫ F dh N
−h/ 2
h/ 2
M P = ∫ PM∧F dh N . m
−h/ 2
n
dl
h
In the same way that for the continuous mediums, one shows that there exist two symmetrical tensors
N and M and a vector Q , defined in the tangent plan in S , such as:
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 4/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
F= N
T =Q⋅
M=n∧ M
Under these conditions is a reference mark whose third component is carried by n , one has (
, =1 ou 2 ):
h/ 2
N =N = ∫ dh
−h/ 2
h/ 2
M =M = ∫ x 3 dh
−h/ 2
h/ 2
Q = ∫ 3 dh
−h/ 2
u x =∑ qi i x
i
w u =∑ q i F i
i
Note:
• The concept of node here is very general and wants to say, in fact, carrying degree of
freedom (that it is of Lagrange or Hermit).
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 5/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
• The concept of displacement is also very general and includes the concept of generalized
displacement including of the translations and rotations.
• for the continuous-current fields on all the field, one uses the values with the nodes
(CHAM_NO of Code_Aster)
u x =∑ ui i x
i
one speaks then about displacements to the nodes, constraints to the nodes or efforts to the
nodes,
Note:
The effort or stress fields are generally calculated at the points of Gauss, if they
continuously are represented it is only at ends of visualization.
• for the other fields, one uses the values in certain points characteristic of the elements (points
of Gauss or nodes).
One speaks then about constraints by elements to the nodes or efforts by elements to the
nodes, or about constraints at the points of Gauss or efforts at the points of Gauss.
NR VY VZ MT MFY MFZ
• for the beams with warping, it is necessary to add bimoment (necessarily in reference mark
“user” with fibre):
BX
• fields of efforts of hull (necessarily in reference mark “user” on the surface):
Moreover, it is sometimes convenient to be able to directly exploit the fields of efforts of beam
and discrete in the total reference mark:
FX FY FZ MX MY MZ
It is also interesting to represent the components of a stress field on the elements of beams or
hulls in the reference mark “user”. For that, one will use the same components as in total
reference mark, although confusion is possible.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 6/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
1.5.3 DEPL_R
Since in Code_Aster,
• a field can be attached only to only one size,
• that the methods finite elements mixed (mixing unknown of standard displacement and
unknown factors of nodal the forces type) are not excluded,
• that the dualisation of the boundary conditions results in having for unknown a comprising
vector of the variables of Lagrange which are nodal forces with the direction where it higher
was specified,
• that it is necessary to be able to carry out any type of linear combination on the nodal forces,
• that the classification of the unknown factors must be the same one as that of the second
members,
nodal forces (dual within the meaning of energy W nodal displacements) have necessarily the same
components as displacements namely:
DX DY DZ DRX DRY MARTINI DRZ
more, for the beams with warping, the degree of freedom associated with bimoment: GRX.
It is the field representative of the state of constraint and making it possible to continue calculations
(geometrical rigidity, nodal forces, etc). It is expressed at the points of Gauss (and is possibly under-
points for the elements of structure). The prefix of this field is SIEF, because according to the
elements, it contains constraints or efforts.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 7/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
F ,M
RAPH_MECA in
SIEF_ELGA into nonlinear
FULL_MECA in reference reference mark
mark “user” “user” *
(*) for the elements of plate and hull, the reference mark “user” is that definite starting from the data of
the user (keyword ANGL_REP or VECTOR in AFFE_CARA_ELEM/COQUE).
It is the field representative of the state of constraint at the points of Gauss (or possibly at the under-
points for the elements of structure). The prefix of this field is SIGM because this field contains only
constraints. It is an extraction of the constraints contained in the field SIEF_ELGA.
They are fields representative of the state of constraint at ends of exploitation (impression or
postprocessing of visualization) to the nodes by element (or possibly at the under-points for the
elements of structure) and to the nodes of the element. According to the elements, they contain
constraints or efforts.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 8/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
by arithmetic
mean with the
nodes of the
F ,M in
* SIEF_NOEU in reference reference mark
quantities to
the nodes by mark “user” “user” *
element
Note:
1) In this case, confusion is possible between the components in reference mark user and those
in total reference mark which bear the same name.
2) The 6 components delivered in the local reference marks by the beams and the hulls contain
possibly worthless terms according to the models used. For the most standard models:
• three worthless terms for the beams,
• two worthless terms for the hulls.
Thus, the stress field will be complete and, especially, it could be enriched each time
modeling requires it (beam with shearing, hull with pinching, etc…)
They are fields containing the efforts on the elements of beam or hull at ends exploitation (impression
or postprocessing of visualization) at the points Gauss, the nodes by elements and the nodes.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 9/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
starting from a
field of F ,M (N , M ,V )
EFGE_ELNO EFGE_ELNO displacement non-available in reference in reference
in linear mark “user” mark “user”
elasticity
by integration
of the
F ,M (N , M ,V )
EFGE_ELNO EFGE_ELNO non-available in reference in reference
constraints into
non-linear mark “user” mark “user”
by arithmetic
mean with the
nodes of the
F ,M (N , M ,V )
* EFGE_NOEU non-available in reference in reference
quantities to
the nodes by mark “user” mark “user”
element
Nodal forces (dual within the meaning of energy W nodal displacements) have the same
components as displacements namely:
DX DY DZ DRX DRY MARTINI DRZ
The option REAC_NODA of the operator CALC_CHAMP carry out a call to FORC_NODA and withdrawn:
• the loading in statics,
• the loading, the inertias and viscous in dynamics (in the facts, the viscous contribution in
dynamics is currently neglected in CALC_CHAMP).
For the solid elements, them FORC_NODA in general have the dimension of a force. It is about a field
on the nodes of the grid where the value in a node is obtained starting from the constraints calculated
on the convergent elements with this node, thus their values thus vary when the grid changes. In the
absence of loading distributed, balance imposes their nullity in an interior node, while they correspond
to the reaction on the supports where a kinematic relation is imposed (case of an imposed
displacement).
In the case of hulls, components DX, DY and DZ they give FORC_NODA (of dimension of a force) in the
total reference mark of the grid. These components are built with the normal efforts and cutting-edges
in the hull. Components DRX, DRY MARTINI and DRZ they give FORC_NODA (of one moment
dimension) in the total reference mark of the grid, built with the bending moments in the hull.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 10/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
2 Kinematics
2.1 Deformations
2.1.1 Continuous medium
In this case, displacements of the structure are represented by a field of vector u with three
components in general.
The deformation (on the assumption of the small disturbances) is defined by the tensor of deformation
by (option EPSI_ELGA and EPSI_ELNO):
1
ij u = u u j , i
2 i,j
One can want to calculate the “mechanical” deformation, i.e. by cutting off thermal dilations (options
EPME_ELGA and EPME_ELNO):
m 1
ij u = u u j ,i −th
2 i,j
In the case of great displacements, the deformations of Green-Lagrange are (options EPSG_ELGA and
EPSG_ELNO ):
1
E ij u = u u j ,i u k ,i u k , j
2 i, j
To which one can want to cut off the thermal deformations (options EPMG_ELGA and EPMG_ELNO):
m 1
E ij u = u u j , iuk , i uk , j −th
2 i, j
For the assumption of the plane deformations (D_PLAN), it is important to note that the condition of
plane deformation is written on the deformation total :
ε33 ( u ) =0
If the total deflection is written as the sum of a purely mechanical deformation and deformation
coming from the variables of order (like the thermal deformation of dilation):
ε33 ( u ) =εm33 +εth33=0
Then it comes naturally that in plane deformations, the deformation mechanics except plan is not
worthless in the presence of a thermal dilation.
m
ε33 ( u ) ≠0
In the theories of traditional beams, each point P beam represents a cross-section. They are thus
the end cells of the torque T s , s of displacement of the presumedly rigid cross-section which
characterizes the displacement of the point P with the curvilinear X-coordinate s . T is the
translation of the centre of inertia of the section, s the vector rotation of the section in this point.
The application of the theorem of virtual work (cf [bib2]) naturally led to define as deformation the
torque , derived from T s , s compared to the curvilinear X-coordinate s :
dT
= ∧
ds
d
=
ds
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 11/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
P(s)
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 12/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
= L T
= t K
Note:
For modelings of beam with taking into account of warping, kinematics is more complicated to
describe, but they lead however to concepts close to those presented above.
• space derivations use too complicated mathematical concepts for the framework of this
document, [R3.07.04],
• the hulls are very often modelled by elements of assembled plates.
In this case, in fact only the material normals are supposed to be rigid. The displacement of these
normals is thus represented by the end cells of a torque T , . T is the translation of the point
located on the average layer, the vector rotation of the normal in this point.
It is clear that the normal component of is worthless (in the case of nonmicropolar mediums). One
introduces, the vector I in the tangent plan defined by:
I=∧n
n
e2
e1
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 13/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
Maybe, decomposition:
T =w nuT
uT is tangent displacement,
w is the arrow.
In the same way that for the beams, the application of the theorem of virtual work (cf [bib2]) led to
define as deformation the unit formed by the tensors E and K and the vector , all these sizes
being defined in the tangent plan by:
1
E = u u
2 , ,
1
K = l l
2 , ,
= l w ,
Note:
There still, there is no universal convention and the disparity of conventions is even larger
than for the tensors of efforts.
= E x3 K
3 =
33 =0
Sizes DEPL_R and DEPL_C as components the degrees of freedom of modeling by finite elements
have and thus do not have necessarily only the components of the fields of displacement which are:
DX DY DZ
For the hulls, we need the three components of the vector of rotation, because the equation with the
finite elements can be expressed only in one total Cartesian reference mark.
2.2.2 EPSI_R
Size EPSI_R represent the structural deformations, therefore it must have, at least, the components:
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 14/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
• fields of deformations of beam (in reference mark “user” with the beam):
They are fields containing the deformations at the points of Gauss and possibly at the under-points
elements.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 15/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
They are fields containing the deformations whatever the modeling at ends of exploitation (impression
or postprocessing of visualization) to the nodes and possibly at the under-points of the elements.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : Contraintes, efforts, forces et déformations Date : 21/07/2017 Page : 16/16
Responsable : FAYOLLE Sébastien Clé : U2.01.05 Révision :
ae95609ea945
by arithmetic
mean with the
nodes of the
* EPMG_NOEU
quantities to Em non-available non-available
the nodes by
element
They are fields containing the deformations generalized on the elements of beam or hull at ends of
exploitation (impression or postprocessing of visualization) at the points of Gauss or the nodes of the
structure.
3 Bibliography
1) F. SIDOROFF: Course of mechanics of the solids Volume 1 E.C.L.
3) C. TRUESDELL, W. NOLL: Encyclopedia of Physics volume III/3 - The non-linear Field Theories
of Mech anics Springer-Verlag, 1965.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2019 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)