Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
1/54
Organization (S): EDF/MTI/MN
Handbook of Référence
R3.07 booklet: Machine elements on average surface
R3.07.03 document
Elements of plate DKT, DST, DKQ, DSQ and Q4g
Summary:
These elements of plate are intended for calculations in small deformations and small displacements of structures
thin curves or plane. In fact plane elements do not take into account the geometrical curvature
structures, contrary to the elements of hull which are curved: it results from it from the parasitic inflections which
can be reduced by using more elements in order to be able to approach the geometries correctly
curves. The formulation is thus simplified by it and numbers it degrees of freedom reduced. These elements are
considered as being among most precise for the calculation of displacements and the modal analysis.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
2/54
Contents
1 Introduction ............................................................................................................................................ 4
2 Formulation ............................................................................................................................................ 5
2.1 Geometry of the elements plates [bib1] .......................................................................................... 5
2.2 Theory of the plates ........................................................................................................................ 6
2.2.1 Kinematics ........................................................................................................................... 6
2.2.2 Law of behavior ............................................................................................................. 7
2.2.3 Taking into account of transverse shearing [bib2] ................................................................. 8
2.2.3.1 The theory known as of Hencky .......................................................................................... 8
2.2.3.2 The theory known as of Reissner (DST, DSQ and Q4g) ........................................................ 9
2.2.3.3 Equivalence of the approaches Hencky-Love-Kirchhoff and Reissner .............................. 9
2.2.3.4 Remarks ................................................................................................................ 9
3 Principle of virtual work ................................................................................................................ 10
3.1 Work of deformation ................................................................................................................... 10
3.1.1 Expression of the resulting efforts ......................................................................................... 10
3.1.2 Relation efforts resulting-deformations .............................................................................. 10
3.1.3 Energy interns elastic of plate ..................................................................................... 11
3.1.4 Remarks ........................................................................................................................... 12
3.2 Work of the forces and couples external ........................................................................................ 12
3.3 Principle of virtual work ............................................................................................................... 13
3.3.1 Kinematics of Hencky ....................................................................................................... 13
3.3.2 Kinematics of Coils-Kirchhoff ............................................................................................ 14
3.3.3 Principal boundary conditions met [bib1] ........................................................... 16
4 numerical Discretization of the variational formulation resulting from the principle of virtual work ............. 17
4.1 Introduction .................................................................................................................................... 17
4.2 Discretization of the field of displacement ..................................................................................... 18
4.2.1 Approach Q4g ...................................................................................................................... 19
4.2.2 Approach DKT, DKQ, DST, DSQ ......................................................................................... 20
4.3 Discretization of the field of deformation ....................................................................................... 21
4.3.1 Discretization of the membrane field of deformation: ..................................................... 22
4.3.2 Discretization of the transverse distortion ............................................................................. 22
4.3.2.1 For the elements Q4g ............................................................................................ 22
4.3.2.2 For the elements of the type DKT, DST ....................................................................... 24
4.3.3 Discretization of the field of deformation of inflection: ........................................................... 27
4.3.3.1 For the elements Q4g ............................................................................................ 27
4.3.3.2 For the elements of the type DKT, DST: .................................................................... 28
4.4 Stamp rigidity .......................................................................................................................... 30
4.4.1 Stamp elementary rigidity for the Q4g elements ........................................................ 30
4.4.2 Stamp elementary rigidity for elements DKT, DKQ .............................................. 31
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
3/54
4.4.3 Stamp elementary rigidity for elements DST, DSQ .............................................. 31
4.4.4 Assembly of the elementary matrices ............................................................................... 32
4.4.4.1 Degrees of freedom ...................................................................................................... 32
4.4.4.2 Fictitious rotations ....................................................................................................... 33
4.5 Stamp of mass ........................................................................................................................... 33
4.5.1 Stamp of traditional elementary mass .............................................................................. 33
4.5.1.1 Element Q4g ............................................................................................................. 33
4.5.1.2 Elements of the type DKT, DST ..................................................................................... 34
4.5.2 Stamp of improved elementary mass ............................................................................. 34
4.5.2.1 Elements of type DKT .............................................................................................. 36
4.5.2.2 Elements of the type DST .............................................................................................. 36
4.5.2.3 Elements of the type Q4g .............................................................................................. 38
4.5.2.4 Notice ................................................................................................................. 38
4.5.3 Assembly of the elementary matrices of mass .............................................................. 38
4.5.4 Stamp of lumpée diagonal mass .................................................................................... 38
4.6 Numerical integration for elasticity ........................................................................................... 40
4.7 Numerical integration for plasticity ......................................................................................... 40
4.8 Discretization of external work ................................................................................................... 41
4.9 Taking into account of the thermal loadings ............................................................................. 43
4.9.1 Thermo elasticity of the plates .............................................................................................. 43
4.9.2 Thermomechanical chaining ................................................................................................ 45
4.9.3 Case-test ................................................................................................................................. 46
5 Establishment of the elements of plate in Code_Aster .................................................................. 47
5.1 Description: ................................................................................................................................... 47
5.2 Introduced use and developments:....................................................................................... 47
5.3 Calculation in linear elasticity:........................................................................................................... 48
5.4 Plastic design ......................................................................................................................... 49
6 Conclusion ............................................................................................................................................ 50
7 Bibliography ......................................................................................................................................... 51
Appendix 1 Plaques orthotropic ................................................................................................................ 52
Appendix 2 Facteurs of transverse correction of shearing for orthotropic plates or stratifiées53
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
4/54
1 Introduction
The elements of hulls and plates are particularly used to model structures
thin where the relationship between dimensions (characteristic thickness/length) is with more than 1/10.
They thus intervene particularly in fields like the civil engineering, the interns of core
REP, vibratory analysis ..... One limits oneself to the framework of small displacements and the small deformations.
Contrary to the elements of hull, the plane elements of plate do not make it possible to take in
hope the geometrical curvature of the structure to be represented and induce parasitic inflections. It is
thus necessary to use a great number of these elements in order to approach correctly
geometry of the structure, and this, more especially as it is curved. On the other hand, one gains in simplicity
of formulation and the number of degrees of freedom is reduced. In addition, the formulations “Discrete
Shear " (DST, DSQ and Q4g) or “Discrete Kirchhoff” (DKT and DKQ) of kinematics, with or without
transverse distortion respectively, allow good results in terms of displacements and
of modal analysis.
The way in which these elements in Code_Aster like certain receipts are established
of use are given to [§5] present note.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
5/54
2 Formulation
2.1
Geometry of the elements plates [bib1]
For the elements of plate one defines a surface of reference, or surfaces average, planes (plane X y
for example) and a thickness H (X, y). This thickness must be small compared to the different one
dimensions (extensions, radii of curvature) of the structure to be modelled. [Figure 2.1-a] below
illustrate our matter.
Solid 3D
Z
H
Y
B
X
L
R1
R2
Thickness H < L, B, R1, R2
Z
y
H
N
X

B
L
Plate
Appear 2.1-a
One attaches to average surface a reference mark orthonormé local Oxyz associated with the tangent plan with
structure different from total reference mark OXYZ. The position of the points of the plate is given by
Cartesian co-ordinates (X, y) of average surface and rise Z compared to this surface.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
6/54
2.2
Theory of the plates
These elements are based on the theory of the plates in small displacements and small deformations according to
which:
2.2.1 Kinematics
The cross-sections which are the sections perpendicular to average surface remain right;
material points located on a normal at not deformed average surface remain on a line
in the deformed configuration. It results from this approach that the fields of displacement vary
linearly in the thickness of the plate. If one indicates by U, v, W displacements of a point
Q (X, y, Z) according to X, y and Z, one has the kinematics of Hencky-Mindlin thus:
U (X, y, Z) U (X, y)
(X, y
y
)
U (X, y)
(X, y
X
X
)








U (X, y, Z) = v (X, y) + Z (X, y) = v (X, y) + Z (X, y
y
X
y
)

U (X, y, Z) W (X, y)

0

W (X, y

Z
)




0

where U, v, W are displacements of average surface and X and there rotations of this surface by
report/ratio with two axes X and y respectively. One prefers to introduce two rotations
X (X, y) = y (X, y), y (X, y) =
- X (X, y).
The three-dimensional deformations in any point, with kinematics introduced previously, are
thus given by:
xx = exx + zxx
yy = eyy + Z yy
2xy = xy = 2exy + 2zxy
2xz = X
2yz = y
where E, E
and E
xx
yy
xy is the membrane deformations of average surface, X and there them
deformations associated with transverse shearings, and xx, yy, xy deformations of inflection
average surface, which is written:
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
7/54
U
exx = X

v
eyy = y
v

U
2exy =
+
X

y
X

xx = X
y
yy = y

X

y

2 xy =
+
y
X

W
X = X
+ X
W

y = y + y
Note:
In the theories of plate the introduction of X and makes it possible it to symmetrize the formulations of
deformations and will see we it thereafter the equilibrium equations. In the theories of hull
one uses X and y rather and the associated couples M X and M y compared to X and Y.
2.2.2 Law of behavior
The behavior of the plates is a behavior 3D in “plane constraints”. The constraint
transversal zz is null bus regarded as negligible compared to the other components of
tensor of the constraints (assumption of the plane constraints). The most general law of behavior
is written then as follows:
xx
xx
E

0
xx
xx









yy


yy
E

0
yy
yy




2nd, 2

and = 0
xy
= C (,) xy = This + zC + C




with E =
xy
=
.


xy

xz
X
0
0
X








yz
y
0
0
y
where
(
C,) is the matrix of local tangent rigidity combining forced plane and distortion
transverse and represents the whole of the internal variables when the behavior is nonlinear.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
8/54
For behaviors where the distortions are uncoupled from the deformations of membrane and from
inflection, (
C,) is put in the form:
H
0
C =

0 H


where
(
H,) is a matrix 3x3 and (
H,) a matrix 2x2. One will remain within the framework of this
assumption.
For an isotropic homogeneous linear behavior elastic, one has as follows:
1 v
0
0
0



v 1
0
0
0
1 -


v

E
0 0
0
0
C =

2

1 - v2
K 1
(- v)

0 0
0
0


2
K 1
(-
v)

0 0
0
0
2


where K is factor of transverse correction of shearing whose significance is given to
following paragraph.
Note:
One does not describe the variation thickness nor that of the transverse deformation zz that one
can however calculate by using the preceding assumption of plane constraints. In addition
no restriction is made on the type of behavior that one can represent.
2.2.3 Taking into account of transverse shearing [bib2]
The taking into account of transverse shearing depends on factors of correction determined a priori by
energy equivalences with models 3D, so that rigidity in shearing
transverse of the model of plate is nearest possible to that defined by the theory of elasticity
three-dimensional. Two theories including the deformation due to the shearing action exist and are
presented in [bib2].
2.2.3.1 The theory known as of Hencky
This theory as that of Coils-Kirchhoff which results from this immediately rests on
kinematics presented to [§2.2.1]. The relation of behavior is usual and the factor of
correction of shearing is worth k=1.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
9/54
Note:
The model of Coils-Kirchhoff (DKT and DKQ): When one does not take into account the distortions
transverses X and y in the theory of Hencky, the model obtained is that of Coils-Kirchhoff.
Two rotations of average surface are then related to displacements of surface
average by the following relation:
W
X = - X
W
y = - y
2.2.3.2 The theory known as of Reissner (DST, DSQ and Q4g)
The second theory, known as of Reissner, is developed starting from the constraints. Variation of
stress membrane (xx, yy and xy) is supposed to be linear in the thickness like in the case of
the theory of Hencky where that results from the linearity of the variation of the deformations of membrane with
the thickness. However, whereas one supposes, in the theory of Hencky, the constant distortion in
the shear thickness and thus stresses, which violates the boundary conditions xz=yz=0 on
faces higher and lower of the plate because of law of behavior stated than the §2.2.2., one
use within the framework of the theory of Reissner the equilibrium equations to deduce the variation from it from
shear stresses in the thickness of the plate, by observing the conditions in particular
of balance on the faces higher and lower of plate. Energy interns model afterwards obtained
resolution of the equilibrium equations in 3D, for inflection only, with the variation of
plane constraints according to Z, reveals, for an elastic material, a relation between the efforts
resulting and rotations and the arrow averages. It is in this relation that the factor appears of
correction of shearing of k=5/6 instead of 1 in the relation which binds the sharp effort to the distortion
for a homogeneous and isotropic plate. Determination of the factors of correction of shearing
for orthotropic plates or laminated plates is left in appendix.
2.2.3.3 Equivalence of the approaches Hencky-Love-Kirchhoff and Reissner
If one assimilates the slopes of average surface X, y with the averages of the slopes in the thickness of
the plate and the arrow W with the average arrow, the only difference between the theory of Hencky and that of
Reissner is the coefficient of transverse correction of shearing of 5/6 instead of 1. This difference
is due to the fact that the starting assumptions are of different nature and especially that the variables
chosen are not the same ones. Indeed, the arrow on average surface is not equal to
average of the arrows on the thickness of the plate. It is thus normal that relations of
behavior which utilizes different variables are not identical.
The fact of having to solve on the level finite elements of the problems in displacements rather than of
problems in constraints by interpolation of displacements leads us to use the approach
equivalent in displacements of the problem of Reissner formulated in constraints.
2.2.3.4 Remarks
Because of preceding equivalence one presents here only the model in displacement for all them
elements. In the facts elements DKT and DKQ are based on the theory of Hencky-Love-Kirchhoff and
the elements DST, DSQ and Q4 are based on the theory of Reissner.
The determination of the factors of correction rests within the framework of another theory, that of Mindlin,
on equivalences of Eigen frequency associated the mode of vibration by transverse shearing.
One obtains k=2/12 then, value very close to 5/6 for the DST elements, DSQ and Q4 in the case
isotropic.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
10/54
Within the framework of plasticity the problem of the choice of the coefficient of correction of shearing
transverse arises because the equivalent approach in displacements of the problem of Reissner formulated in
constraints utilizes the non-linearity of the behavior. One cannot thus deduce some, like
it is the case for elastic materials a value of the coefficient of correction of shearing
transverse. Plasticity is thus not developed for these elements.
3
Principle of virtual work
3.1
Work of deformation
The general expression of the work of deformation 3D for a plate is worth:
H/2
W
(+ + + + D
) V
def = xx xx
yy yy
xy xy
X xz
y yz
S - H/2
where S is average surface and the position in the thickness of the plate varies between ­ H/2 and +h/2.
3.1.1 Expression of the resulting efforts
By adopting the kinematics of [§2.2.1], one identifies the work of the interior efforts:
W
E NR
+ E NR + 2nd NR + M + M + 2 M + T + T) dS
def = (xx xx
yy
yy
xy
xy
xx
xx
yy
yy
xy
xy
X X
y y
S
where:
NR

M

/2

/2

xx +h xx
xx +h xx
T
H/2
X
+
xz
NR = NR yy = yy dz; M = Myy = yyzdz; T = =

dz.
T



y
yz
- H/2


- H/2

NR
- H/2
xy
xy
Mxy
xy
NR, NR, NR
xx
yy
xy is the efforts resulting from membrane (in NR/m);
M, M, M
xx
yy
xy is the efforts resulting from inflection or moments (in NR);
T, T
X
y, are the efforts resulting from shearing or sharp efforts (in NR/m);
3.1.2 Relation efforts resulting-deformations
The expression of the work of deformation is also written:
H/2
H/2
W
[(
C,) D
] V
[eCe Z

EC.
zCe Z
2 C C D
] V
def =
=
+
+
+
+


S - H/2
S - H/2
where (
C,) is the local matrix of behavior.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
11/54
By using the expression obtained for Wdef in the preceding paragraph one finds the relation following
between the resulting efforts and the deformations:
NR = H E + H
m
MF
+h/2
+h/2
+h/2
M = H E + H
2
MF
F with H
=
Hdz, H
=
Hzdz, H =
H
m

MF

F
Z dz
T = H
- H/2
- H/2
- H/2
ct
G11
0
and: Hct =

0
G22
and:
E

xx
xx
X
E = eyy, = yy, =




y
2exy
2xy
The matrices H, H and H
m
F
ct is the matrices of rigidity out of membrane, inflection and shearing
transverse, respectively. The Hmf matrix is a matrix of rigidity of coupling between
membrane and the inflection.
For an isotropic homogeneous elastic behavior of plate these matrices have as an expression:




1 v
0
3
1 v
0
Eh
Eh
kEh 1

0
H
=


v 1
0
, H =


v 1
0
, H
m
F
ct =

, and
1 - 2
v
1 -
v
12 1
(- 2
v)
1 -
v
2 1
(+ v) 0

1

0 0
0 0
2


2
Hmf = 0 bus there is material symmetry compared to the z=0 plan.
For an orthotropic material, the behavior is given in appendix.
3.1.3 Energy interns elastic of plate
Taking into account the preceding remarks, energy interns elastic plate is expressed more
usually for this kind of geometry in the following way:
1
int =
[E (H E
m + H
)
MF
+ (H E
MF + H)
F
+ H

] dS.
2
ct
S
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
12/54
3.1.4 Remarks
Relations binding H, H, H
with H and H
with H
m
F
MF
ct
are valid whatever the law of
elastic behavior, with anelastic deformations (thermoelasticity, plasticity,….).
For a plate made up of NR orthotropic layers in elasticity, matrices
H, H, H
and H
m
F
MF
ct is written:
NR
NR
NR
NR
1
H = H H, H
3
3
I
= H H, H
I I
= (zi+1 - Z) H, H
I
= H H
m
I
MF
I
F
I
ct
I
I
3

i=1
i=1
i=1
i=1
1
where: H = Z +1 - Z, = (Z +1 + Z
I
I
I
I
I
I) and H, H
and
2
I
I
the matrices H represent
H for the layer
I.
The homogenization for multi-layer hulls can lead to matrices of rigidity of
membrane and of inflection nonproportional of the type:
C
C
0
D
D
0
1111
1122

1111
1122

G11
0
H = C
C
1122
2222
0, H = D
D
1122
2222
0, H
m
F
ct =

0
G22

0
0
C1212

0
0
D1212
for which one cannot find equivalent values of the Young modulus and thickness
allowing to find the traditional expressions of rigidity, cf [bib7].
3.2
Work of the forces and couples external
The work of the forces and couples being exerted on the plate is expressed in the following way:
+h/2
+h/2
W
F. D
U V
F. D
U.S.
F. D
U zds
ext. =
v
+
S
+


C
S - H/2
S
C - H/2
where F, F, F
v
S
C are the voluminal, surface efforts and of contour being exerted on the plate,
respectively. C is the part of the contour of the plate on which the efforts of contour FC are
applied. With the kinematics of [§2.2.1], one determines as follows:
W
(F U + F v + F W + C + c) dS + (U + v + W + +) ds
ext. =
X
y
Z
X X
y y
X
y
Z
X X
y y
S
C
= (F U + F v + F W + C - c) dS + (U + v + W + -) ds
X
y
Z
y X
X y
X
y
Z
y X
X y
S
C
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
13/54
· where are present on the plate:
F, F, F
X
y
Z: surface forces acting according to X, y and Z
+h/2
F =
. dz
I
+
F E
F .e
v
I
S
I where ex and E are there the basic vectors of the tangent plan and E Z them
- H/2
normal vector.
C, C
X
y: surface couples acting around axes X and Y.
+h/2
H
C =
ze F E dz
I
+ ± E
(
).
(
F) .e
Z
v
I
Z
S
I where E, E, E are the basic vectors
2
X
y
Z
- H/2
previously definite.
· and where are present on the contour of the plate:
X, y, Z: linear forces acting according to X, y and Z
+h/2
I =
. dz
F E
C
I
where ex, ey, ez are the basic vectors previously definite.
- H/2
X, y: linear couples around axes X and Y.
+h/2
I =
zez F E dz
(
).
C
I
where E, E, E
X
y
Z are the basic vectors previously definite.
- H/2
Note:
The moments compared to Z are null.
3.3
Principle of virtual work
It is written in the following way: W
= W
ext.
def for all virtual displacements and rotations
acceptable.
3.3.1 Kinematics of Hencky
With this kinematics, it results after integration by parts of work of deformation them
equilibrium equations static of the following plates:
NR
+ NR
+ F
xx, X
xy, y
X = 0,
· For the efforts: NR
+ NR
+ F
yy, y
xy, X
y = 0,
T
+ T
+ F
X, X
y, y
Z = 0.
M
+ M
- T + C
xx, X
xy, y
X
y = 0,
· For the couples: M
+ M
- T - C
yy, y
xy, X
y
X = 0.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
14/54
as well as the boundary conditions following on contour C of S:
NR N + NR N
xx X
xy y = X,
U = U,
NR N + NR N
yy y
xy X = y,
v = v,
T N + T N
X X
y y = Z
,
or W = W,
M N + M N


xx X
xy y = y,
X = y,
M N + M N
yy y
xy X = - X.
y = - X.
where nx and ny are the cosine directors of the normal with C directed towards the outside of the plate.
The physical interpretation of these efforts (NR, T and M) starting from the preceding equations is given
below:
Z
Ty
Myy
P
P
Tx
Nyy
Myx
y
Nyx
Nxy
Mxx
X
Nxx
Mxy
Appear 3.3.1-a: Efforts resulting for an element from plate
Note:
NR, NR
xx
yy represents the tractive efforts and NR xy plane shearing. M xx and M yy
the couples of inflection and M xy represent the torque. Tx and Ty, are the efforts of
transverse shearing.
3.3.2 Kinematics of Coils-Kirchhoff
One recalls that within the framework of this kinematics, one with the following relation binding the derivative of
W
X = - X
arrow with rotations:
W. After a double integration by parts of the work of deformation,
y = - y
one obtains the following equilibrium equations static:
NR
+ NR
+ F
xx, X
xy, y
X = 0,
· For the efforts of membrane: NR
+ NR
+ F
yy, y
xy, X
y = 0,
· For the transverse shearing and bending stresses:
M
+ 2 M
+ M
+ F + C
- C
xx, xx
xy, xy
yy, yy
Z
y, X
X, y = 0,
M
+ M
- T + C
xx, X
xy, y
X
y = 0,
M
+ M
- T - C
yy, y
xy, X
y
X = 0.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
15/54
as well as the boundary conditions on contour C and with the angular points O of contour C of S:
NR N + NR N
xx X
xy y = X
,
U = U,
NR N + NR N
yy y
xy X = y,
v = v,
T + M
N
NS, S = Z
- N, S,
or W = W,
Mn = S,
= - W
N
N = S.
,
M (O+) - M (O) = - [(O+) - (O
NS
NS
N
N
-)].
T = T N + T N
N
X X
y y,
with M
= M N2 + 2 M N N + M N2
N
xx X
xy X y
yy y,
.
M
= - M N N + M (N2 - N2) + M N N
NS
xx X y
xy
X
y
yy X Y.
Surface average S
S
Contour
Million (+)
S Million (-)
Discontinuity
Appear 3.3.2-a: Condition in the limits with angular points for an element of plate
Note:
The kinematics of Coils-Kirchhoff implies that on the contour of the plate the shearing force
transverse is related to the torque. It is noted that the command of the equilibrium equations of
inflection is higher than with the kinematics of Hencky. Thus, to choose the kinematics of
Coil-Kirchhoff come down to increase the degree of the functions of interpolation because a regularity is needed
larger for the terms of arrow compared to the terms of membrane because of presence
from derived seconds of the arrow in the expression of the work of the deformations. No element
of plate of Code_Aster does not use this kinematics. One can thus have differences between
results obtained with the elements of Code_Aster and the analytical results obtained in
using the kinematics of Coils-Kirchhoff for structures with angular contours.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
16/54
3.3.3 Principal boundary conditions met [bib1]
S
Surface average S
y
nx
Contour
X
ny
N
S, S
N, N
Appear 3.3.3-a: Condition in the limits for an element of plate
The boundary conditions frequently met are gathered in the table which follows. They
are given for the kinematics of Hencky in the reference mark defined by S and the normal external with
plate:
Embedding
Simple support
Free edge
Symmetry compared to
Antisymetry by
an axis S
report/ratio with an axis S
U = 0,
one = 0,
one = 0,
custom = 0,
v = 0,
W = 0,

W = 0,
S = 0.
W = 0,
N
= 0.
N
= 0.
S
= 0,
N
= 0.

=
=
S = 0,
S
0,
S
0,
N = 0,

= 0,
= 0,
S = 0.
N
Z
S = 0.
Z = 0,
N = 0.
S = 0,
N = 0
U = one + vn; U = - one + vn
N
X
y
S
y
X,
= N + N; = - N + N
N
X X
y y
S
X y
y X,
= N2 +
2
N N + N2
N
X X
xy X y
y y,
with: = - N N + (N2 - N2) + N N
S
X X y
xy
X
y
y X y,
= N2 + 2 N N + N2
N
X X
xy X y
y y,
= - N N + (N2 - N2) + N N
S
X X y
xy
X
y
y X Y.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
17/54
S =
-,
Note: one has
N.
N = S.
4
Numerical discretization of the variational formulation
exit of the principle of virtual work
4.1 Introduction
By exploiting the law of behavior, the virtual work of the interior efforts is written (with Hmf =0
until [§4.4], which does not remove anything with the general information following results, but allows to reduce them
notations):
Wint = eH E + H + H
m
F
ct dS
(
)
S
U





, X


X, X

W, X + X
with: E =
v, y, =
y, y
, =
.




W, y + y
U, y + v, X
X, y + y, X
It results from it that the elements of plate are elements with five degrees of freedom per node. These
degrees of freedom are displacements in the plan of the element U and v, except plan W and both
rotations X and Y.
The elements DKT and DST are triangular isoparametric elements. Elements DKQ, DSQ
and Q4 are quadrilateral isoparametric elements. They are represented below:
4
3
3
1
y
2
1
2
X
Appear 4.1-a: Eléments realities
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
18/54
The elements of reference are presented below:


(0,1)
(- 1,1)
(1,1)
3
4
3
1
2
(0,0)
(1,0)


1
2
(1, - 1)
(- 1, - 1)
Appear 4.1-b: Eléments of reference triangle and quadrangle
One defines the reduced reference mark of the element as the reference mark (,) of the element of reference. The local reference mark
element, in its plan (X, y) is defined by the user. The X1 direction of this local reference mark is
projection of a direction of reference D in the field of the element. This direction of reference D is
chosen by the user who defines it by two nautical angles in the total reference mark. The normal NR with
plan of the element (12 13 for a triangle numbered 123 and 12 14 for a numbered quadrangle
1234) fix the second direction. The vector product of the two vectors previously definite
Y1=N X1 makes it possible to define the local trihedron in which will be expressed the generalized efforts
representing the state of stresses. The user will have to take care that the selected reference axis
do not find parallel with the normal of certain elements of plate. By defect, direction of
reference D is axis X of the total reference mark of definition of the grid.
The essential difference between elements DKT, DKQ on the one hand and DST, DSQ, Q4 on the other hand comes
owing to the fact that for the first the transverse distortion is null is still = 0. The difference between Q4
and elements DST and DSQ comes from a choice different of interpolation for the representation of
transverse shearing.
4.2
Discretization of the field of displacement
If one discretizes the fields of displacement in the usual way for elements
isoparametric i.e.:
NR
NR
NR
NR
NR
U = NR U, v = NR v, W = NR W, = NR, = NR
I I
I I
I I
X
I xi
y
I yi,
i=1
i=1
i=1
i=1
i=1
and that one introduces this discretization into the variational formulation of [§4.1] it a blocking results from it
in transverse shearing analyzed in [bib1] which returns the solution in inflection controlled by the effects of
transverse shearing, and not by the inflection, when the thickness of the plate becomes small by report/ratio
with its characteristic dimension.
To cure this disadvantage the variational form presented in introduction is slightly
modified so that:
W
1
int =
eH E
m + Hf + Hct dS =
eH E
m + Hf + TH
T dS


-
(
)
(
)
ct
S
S
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
19/54
where are deformations of substitution checking = in a weak way (integral on the sides of
J
the element) and such as T = Hct. One checks thus that on the sides ij of the element (S -)
S ds =

0
I
with S = W, S + S.
Two approaches are then possible; in the first, that of the element Q4, one uses
bilinear discretization of the fields of displacement and the fact that is constant on the sides of
the element. The relations on the sides ij then make it possible to express the values of on the sides in
function of the degrees of freedom of inflection. In the second approach, which is that of the elements of the type
DKT and DST, one uses the weak formulation of the preceding paragraph which makes it possible to bind the inflection to
shearing forces to deduce the interpolation from it from the terms of inflection.
4.2.1 Approach
Q4g
It rests on the linear discretization of the fields of displacement presented above:
NR
NR
NR
NR
NR
U = NR U, v = NR v, W = NR W, = NR, = NR
I I
I I
I I
X
I xi
y
I yi,
i=1
i=1
i=1
i=1
i=1
where the functions Nor are given below.
Nor (i=1, N)

Q4
i=1 with 4
4
3
1
N1 (,) =
1
(-) 1
(-)
4
1
NR (,) =
1
(+) 1
(-)

2
4
1
NR 3 (,) =
1
(+) 1
(+)
4
1
2
1
NR 4 (,) =
1
(-) 1
(+)
4
Functions Nor for the Q4 elements
Note:
1
One also notes Nor (,) = (1+ I
) (1+ I
) with (,
1,
2,
3) = (-, -) and
4
4
111 1
(,
1,
2,
3) = (-, -,)
4
1.111.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
20/54
4.2.2 Approach DKT, DKQ, DST, DSQ
Like T = M, + M
and T
,
= M, + M
X
xx X
xy y
y
yy y
xy, X and M = Hf one deduces from it that is defined
according to the derived seconds of X and y via two equilibrium equations internal
and of the law of behavior in inflection. The discretization retained for X and y, such as S is
quadratic on the sides and N linear, then utilizes of the quadratic functions of forms
incomplete in the form:
NR
2 NR
NR
2 NR
X = Nk xk + xk
P K, y = Nk yk + yk
P K with P = PC and P = P S
xk
K K
yk
K K
K =1
K = NR +1
K =1
K = NR +1
where Ck and Sk are the cosine and directing sines on the side ij to which belongs the node K defined by:
/
C = X/L = (X - X)/L; S = y/L = (y - y)/L; L = (x2 + y2 1 2
K
ji
K
J
I
K
K
ij
K
J
I
K
K
ji
ji)
.
Note:
To introduce the preceding discretization amounts adding like degrees of freedom to the element of
rotations K in the middle of the sides K of the element. Indeed, rotations S and N such as:
C
S
S
X
=

S - C
N
y
are quadratic for S and linear for N with:
S = (1 - S) if + ssj + 4s (1 - S
) K; N = (1 - S) nor + snj where 0 S = S/Lk 1.
1
1
One observes thus that: sk = S (S =) = (if + sj) + K.
2
2
J
It is the relation (S -)
S ds =

0 with S = W, S + S which will make it possible to eliminate the degrees from
I
freedom additional and to express them according to displacements and of nodal rotations.
sk
sj
K
nj
if
nor
1/2 (si+sj)
I
K
J
S
I
J
S
Variation of S
Variation of N
Appear 4.2.2-a: Variations of S and N
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
21/54
Nor (i=1, N)
Pi (i=n+1,2n)
DKT, DST
i=1 with 3
i=4 with 6

NR (,) = = 1 - -
P (,) = 4
3
1
4
NR (,) =
P (,) = 4
5
2
5
6
NR 3 (,) =
P6 (,) = 4
1
4
2

DKQ, DSQ
i=1 with 4
i=5 with 8

1
1
4
7
3
NR
2
1 (,) =
1
(-) 1
(-)
P (,) =
1
(-) 1
(-)
4
5
2
1
1
NR (,) =
1
(+) 1
(-)
P
2
(,) =
1
(-) 1
(+)
8
6
2
6

4
2
1
1
NR
2
3 (,) =
1
(+) 1
(+)
P (,) =
1
(-) 1
(+)
4
7
2
1
5
2
1
1
NR
2
4 (,) =
1
(-) 1
(+)
P (,) =
1
(-) 1
(+)
4
8
2
Functions Ni and pi for elements DKT, DST, DKQ, DSQ
4.3
Discretization of the field of deformation
The matrix jacobienne J (,) is:
NR
NR


NR X
NR y
X
y
I, I

,
,
,

I
I
J
J
i=1
i=1
11
12
J =
=
=
NR
NR

.
X
y
,
,



J
J
NR X
NR y
21
22
I, I


I, I
i=1
i=1


Moreover:



X

J
J
11
12
1
1 J22
- J

12
-
=
J

front
EC. J =
= J =

or J = det J = J J
- J J


J
J

11 22
12 21
21
22
J - J
J
21
11
y

It is pointed out that the field of displacement is discretized by:
W
W

0

U
NR

the U.K.
NR
K
2 NR





= NR,
K (
)
and = NR
(,) + [P
(,)
], the term enters
v
X
K
xk
xk
K
=1
vk
K

k=1


K = NR +1



P
y
yk
yk
(,)
hooks being present for the elements of the type DKT, DST, but not for the Q4 elements.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
22/54
4.3.1 Discretization of the membrane field of deformation:
NR
NR
E
= U, = NR (,) the U.K.
,
= (J NR
11
, + J
NR
) the U.K.
xx
X
K X
K
12
K,
,

K =1
K =1
NR
NR
E
= v, = NR (,) vk
,
= (J NR
21
, + J
NR
) v K
yy
y
K y
K
22
K,
,

K =1
K =1
NR
2nd
= U, + v, =
NR
(,) the U.K.
,
+ NR (,) vk
xy
X
y
K y
K, X
K =1
NR
= (J NR
K
K
21
, + J
NR
K
22
K,) U + (J NR
11
K + J NR
)
,
12
K,
v

K =1
Maybe in matric form:
E
xx
NR


U
E
K
yy = B

U
= is the membrane field of displacement to the node K


mk
K where the U.K.
v
2nd
K 1
K

=
xy
and:
J NR
11
K, + J NR
0


12
K,

Bmk =
0
J NR
21
K, + J
NR
22
K,


J NR
21
K, + J
NR
J NR
22
K,
11
K, + J
NR
12
K,
u1


v
1
The matrix of passage of the membrane deformations to the field of displacement U


m =!
U
NR

v
NR
in the plan of the element is written as follows: B [3×2N = (B
B
m
]
m1
mn).
4.3.2 Discretization of the transverse distortion
4.3.2.1 For the Q4g elements
One linearly discretizes the constant field by side so that:
1 -
1
12
+ 34




+



= = 2
2
.

1 -
1
23
-



41
+

2
2

Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
23/54
+1
(- (
W +))
D = 0;
,
- 1
By using the relations then:
,
+1
(- (
W +))
,
D = 0
- 1
1
ij = (wj - I
W + I +


J);
2
it is established that:
for (ij) = (12,34) and (kp) = (23,41).
1
kp
=
(wp - wk + p + K);
2
By deferring the two results above in the expression of, one establishes that:
w1



1



1


= = B U where U =! and B = (B, B

1 “
)



NR



wN

NR




NR
NR K, K NR K,
0

with B K =
.
NR K,
0
K Nk,
w1


x1


y1
I
J
J
11
12 xi


Like
=

one deduces from it that = B U
=! and

F where U F
I
J
J
21
22 yi
W
NR


xN



yN
NR K, K Nk J
,
11
K Nk J
,
12
B = (B, B


)
1 “
NR with B K =
.
NR
K, K NR K J
,
21
K Nk J
,
22
X J
J
11
12
Finally:
=

= B U with B
= jB


=
C
F
C [2×3N]
.
y
J
J
21
22
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
24/54
4.3.2.2 For the elements of the type DKT, DST
With regard to the transverse distortions one deduces from
T = M, + M
and T
,
= M, + M
X
xx X
xy y
y
yy y
xy, X with M = H F that T = H F, xx where:
T, xx = (X, xx X, yy X, xy y, xx y, yy y, xy) and
H
H
2H
H
H
H
11
33
13
13
23
12 + H33
H F =


where H
H
H
H
ij is the terms (I, J) of
13
23
12 + H
H
H
2H
33
33
22
23

Hf.
NR
2 NR
NR
NR

= NR (,) + P (,) = NR (,) + (J2 P
2
2
11
+ J J P
11 12
+ J P
)
,
X, xx
K, xx xk
xk, xx K
K, xx xk
xk,
xk
,
12 xk
,
K
K =1
K = NR +1
K =1
K = NR +1
NR
2 NR
NR
NR

= NR (,) + P (,) = NR (,) +
2
2
(J P

, + 2 J
J P, + J P
)
,
X, yy
K, yy xk
xk, yy K
K, yy xk
21 xk
21 22 xk
22 xk,
K
K =1
K = NR +1
K =1
K = NR +1
NR
2 NR



,
= NR (,)
,
+ P (,)
X xy
K xy
xk
xk, xy
K
K =1
K = NR +1
NR
NR
= NR (,)

,
+ (J J P, + [J J + J J] P, + J J P),
K xy
xk
11 21 xk
11 22
12 21
xk
11 21 xk,
K
K =1
K = NR +1
NR
2 NR
NR
NR


(,) + P (,) = NR (,) + (J2 P
2
2
11
+ J J P
11 12
+ J P
)


12

,
=
NR
y xx
K, xx
yk
yk, xx
K
K, xx
yk
yk,
yk,
yk,
K
K =1
K = NR +1
K =1
K = NR +1
NR
2 NR
NR
NR

= NR (,) + P (,) = NR (,) + (J2 P
2
2
J J P
+ J P
),
21
+
y, yy
K, yy
yk
yk, yy
K
K, yy
yk
yk,
21 22
yk,
22
yk,
K
K =1
K = NR +1
K =1
K = NR +1
NR
2 NR

= NR (,) + P (,)
y, xy
K, xy
yk
yk, xy
K
K =1
K = NR +1
NR
NR
= NR (,) + (J J P + [J J + J J] P + J J P)
K, xy
yk
11 21 yk,
11 22
12 21
yk,
11 21 yk,
K
K =1
K = NR +1
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
25/54
that is to say still in matric form that:
X, xx


X, yy

T = H X, xy =
F y, xx



y, yy

y, xy
0
j2 NR
2 J J NR
j2 NR
0

11
K
,
+ 11 12 K
,
+ 12 K,


0
j2 NR
2 J J NR
j2 NR
0

21
K
,
+ 21 22 K, + 22 K
,
W
NR 0
J J NR
[J J
J J] NR
J J NR
0

K


11 21
K
,
+ 11 22 + 12 21 K
,
+
H
11 21
K
,




F
= 0
0
j2 NR
2 J J NR
+ j2 NR
xk
K 1
11
K
,
+ 11 12 K,
12
K,



0
0
j2 NR
+ 2 J J NR
+ j2 NR
yk
21
K,
21 22
K,
22
K,



0
0
J J NR
+ [J J + J J] NR
+ J J NR

11 21
K,
11 22
12 21
K,
11 21
K,

C (J2 P
+ 2 J J P + j2 P)

K
11 K,
11 12 K,
12 K
,


C (J2 P
2
2
)

K
21 K,
+ J J P
21 22 K,
+


J P
22 K
,

2 NR
C (J J P
[
]
)
K
11 21 K,
+ J J
11 22 + J
J P
12 21
K,
+ J J P
+ H


11 21 K
,

=
F K
S (J2 P
2
2
)
K
11 K,
+ J J P
11 12 K,
+ J P
K = NR +1


12 K
,


S (J2 P
2
2
)

K
21 K,
+ J J P
21 22 K,
+


J P
22 K
,



S (J J P
[
+ J J] P
)
K, +

J J P
K
11 21 K,
+

J J
11 22
12 21
11 21 K,
C P
K
K,


C P
W
K K,
NR
K
2


NR
C P
NR
2 NR
H
K
K,
P

H T

+ H T T H P U
H T T



B U
B



F
2
ck
K =
F
F
F
+ F 2
= C F +
F
2


H
P U
F
F K
xk +



K =
S P
F
F K

F K

C
K =1
K = NR 1


+
K K,

K

=1
K = NR +1
yk


S P
K
K,

S P
K
K,

2
2
J
J
2 J J

T
11
12
11 12


2
0
where T
(T
T
2
2
=
C NR +
)
(
)
1 “c2 NR and T2 =

= J
J
2 J J

with T
.
0
T
2
21
22
21 22
2



j11 j21 j12 j22 j11 j22 + 12
J j21
J
We use then the relation (S -)
S ds =

0 with S = W, S + S for each side ij of
I
the element which makes it possible to obtain K since she is still written:
L
2
W - W
K
+
(C + S + C + S) + L = L
J
I
K xi
K yi
K xj
K yj
K K
K sk where:
2
3

- 1
-
= (C
S) = (C
S) H T = (C
S) H 1 [B U + B
]
sk
K
K
K
K
ct
K
K
ct

C

F
C
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
26/54
The relation above is still written in matric form: WITH = A U

W
F
L
0
0
NR 1
LN 1CN 1 LN 1SN 1
2
+

+
+
+
+
with: With
=
0
0
#
-
H - 1B



!
!
ct

3
C

0
0
L2N L2 C
NR
2
L
NR
2 S
N2 NR
and:
- 2 L C
L
S
2
L
C
L
S
0
0
0

NR +1 NR +1
NR +1 NR +1
NR +1 NR +1
NR +1 NR +1


0
0
0
- 2 L C
L S
2
L C
L S

K +1 K +1
K +1 K +1
K +1 K +1
K +1 K +1
0
0
0
0
0
0
- 2 L
C
L
S

2 NR - 1 2 NR - 1
2 NR - 1 2 NR - 1


1 2
L C
L S
0
0
0
0




With
2 NR
2 NR
2 NR
2 NR
= -
W
2
0
0
0





0
0
0



2
L
C
L

S

2 NR - 1 2 NR - 1
2 NR - 1 2 NR - 1



“- 2 L C
L S


2 NR
2 NR
2 NR
2 NR
L C
L
S

NR +1
NR +1
NR +1 NR +1

1
+
!
!
H B
ct
C

L C
L S


2 NR
2 NR
2 NR
2 NR
Thus =

front
EC.
1
=
-
WITH U
With
WITH A
F
W, which implies T = [B + B A
U
]
C
C

F.
Note:
For the DST elements, this expression is simplified a little since B
= 0 because of linearity of
C
functions of form NR (k=1,2,3).
K
This expression is simpler for elements DKT and DKQ since they are without transverse distortion,
1 0

0


i.e. = 0, which implies A =

0 #
0 and

0 0 1
- 2/LN
C
+1
NR
S
+1
NR
2/L
+1
NR
C
+1
NR
S
+1
NR +
0
0
0
1



0
0
0
- 2/L

K
C
+1
K
S
+1
K
2/L
+1
K
C
+1
K
S
+1
K +1


0
0
0
0
0
0
- 2/L

2 NR
C
- 1
2 NR
S
- 1
2 NR - 1



3
2/L
C
S
0
0
0
0




With
2 NR
2 NR
2 NR
W = - 4
0
0
0





0
0
0



2/L


2 NR
C
- 1
2 NR
S
- 1
2 NR - 1





“- 2/L2N
C2N
S2N

It is also noticed that for elements DKT the expression of the efforts sharp is calculated from
balance and not starting from the behavior (on the basis of the behavior one would find a value zero
sharp efforts what would not make it possible to carry out balance!). It results from it according to the §3.1.1 from
transverse stress shear nonnull in the thickness of the plate that one is in formulation
DKT or DST.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
27/54
4.3.3 Discretization of the field of deformation of inflection:
4.3.3.1 For the Q4g elements
The relation binding the deformations of inflection to the field of displacement of inflection is written:
NR
NR
xx = X, X = j11x, + j12x, = j11 Nk, xk + j12 Nk,



xk,
K =1
K =1
NR
NR
yy = y, y = j21y, + j22y, = j21 Nk, yk + j22 Nk,



yk,
K =1
K =1
NR
NR

2 xy = y, X + X, y = j11y, + j12y, + j21x, + j22x, = J NR
21
K
, xk + j22 NR K, xk
K 1
=
K 1
=
NR
NR
+ 11
J Nk, yk + 12
J Nk, yk.
K 1
=
K 1
=
That is to say still in matric form:



xx
W
NR


K


yy = B fk U fk where U fk = xk represents the field of displacement of inflection to the node K,





K
2
1
=
xy
yk
with:
0 J NR
0

11
K, + J
NR
12
K

,

B fk = 0
0
J NR
21
K, + J
NR
22
K.
0 J NR


21
K, + J
NR
22
K
J NR
,
11
K, + J
NR
12
K,
w1


x1

y
1
The matrix of passage of the field of displacement of inflection U F =! with the deformations of
W
NR


xN



yN
inflection is written then: B
= B, B
FN]
(FfN)
3 3
1
.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
28/54
4.3.3.2 For the elements of the type DKT, DST:
The relation binding the deformations of inflection to the field of displacement of inflection is written:
NR
2 NR
NR
2 NR
xx = X, X = j11,
12
,
(
11
,
,)
(
12
,

X + J X = J
NR K xk +
xk
P K + J
NR K xk +
xk
P, K),
K =1
K = NR +1
K =1
K = NR +1
NR
2 NR
NR
2 NR


,
21
,
22
,
(
21
,
,)
(
22
yy = there y = J y + J y = J
NR K yk +
yk
P K + J
NR K, yk +
Pyk, K),
K =1
K = NR +1
K =1
K = NR +1

2 xy = y, X + X, y = j11y + J
,
12 y
+ J
,
21x
+ J
,
22 X, =
NR
2 NR
NR
2 NR
NR
2 NR
J (
21 NR
+
+ 22
+
+ 11
+
K, xk
xk
P, K) J (
NR K, xk
xk
P, K) J (
NR K, yk
yk
P, K)
K =1
K = NR +1
K =1
K = NR +1
K =1
K = NR +1
NR
2 NR
+ J (
12 NR
+ P
K,
yk
yk, K).


K =1
K = NR +1
For elements DKT, DKQ:
In matric form the preceding relation is also written by introducing the relation = A U
F:

B
B

xx
j11
X + j12
X



yy =
j21B
y + j22B
y
U F = B F 3×3N U
[
]
F where




2 xy
J
11B
y + j12By + j21B
X + j22B
X
w1


x1

y
1
U F =! represent the field of displacement in inflection for the element with:
W
NR


xN



yN
6 NR
P +1C
,
NR
6P
+1
2 NR C
,
2 NR
3
B
2
2

X = (
-
, N1, - (NR
P +1C
,
NR +1 + P2 NR C
,
2 NR),
4 LN
4L
+1
2 NR
4
3
- (NR
P +1C
,
NR
S
+1 NR +1 + P2N C
,
2 NR S2 NR),
,
4

6 NR
P +K C
,
NR +k
6 NR
P +k-1C
,
NR +k 1
3
-
-, NR
2
2
K
, -
(NR
P +K C
+ P
C
),
4 L
, NR +k
NR +k - 1
,
NR +k - 1
NR +k
4 LN+k-
4
1
3
- (PN+k C
S
,
+ P
NR +K NR +K
NR +K - 1 C
S
,
NR +K - 1 NR +K - 1),
4

(K = 2., NR))
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
29/54
6 NR
P +1C
,
NR
6P
+1
2 NR C2
3
B
,
NR
2
2

X
= (
-
, N1, - (NR
P +1C
,
NR +1 + P2 NR C
,
2 NR),
4LN
4L
+1
2 NR
4
3
- (NR
P +1C
,
NR
S
+1 NR +1 + P2N C
,
2 NR S2 NR),
,
4

6 NR
P +K C
,
NR +k
6 NR
P +k-1C
,
NR +k 1
3
-
-, NR
2
2
K
, -
(NR
P +K C
+ P
C
),
4L
, NR +k
NR +k - 1
,
NR +k - 1
NR +k
4LN+k-
4
1
3
- (PN+k C
S
,
+ P
NR +K NR +K
NR +K - 1 C
S
,
NR +K - 1 NR +K - 1),
4

(K = 2., NR))
6 NR
P +1S
,
NR
6P
+1
2 NR S
,
2 NR
3
B y = (
-
, - (NR
P +1C
,
NR
S
+1 NR +1 + P2 NR C
,
2 NR S2 NR),
4 LN
4 L
+1
2 NR
4
3
NR
2
2
1
, -
(NR
P +1 S
,
NR +1 + P2 NR S
,
2 NR),
,
4

6 NR
P +K S
,
NR +k
6 NR
P +k-1S
,
NR +k 1
3
-
-, - (NR
P +K C
,
NR
S
+ P
C
S
),
4 L
+k NR +k
NR +K - 1, NR +K - 1 NR +K - 1
NR +k
4LN+k-
4
1
3
NR
2
2
K, -
(NR
P +K S
, NR +K + NR
P +k
S
- 1, NR +k-1),
4

(K = 2., NR))
6 NR
P +1S
,
NR
6P
+1
2 NR S2
3
B
,
NR
y = (
-
, - (NR
P +1C
,
NR
S
+1 NR +1 + P2 NR C
,
2 NR S2 NR),
4 LN
4 L
+1
2 NR
4
3
NR
2
2
1
, -
(NR
P +1S
,
NR +1 + P2 NR S
,
2 NR),
,
4

6 NR
P +K S
,
NR +k
6 NR
P +k-1S
,
NR +k 1
3
-
-, - (NR
P +K C
,
NR
S
+ P
C
S
),
4 L
+k NR +k
NR +K - 1, NR +K - 1 NR +K - 1
NR +k
4LN+k-
4
1
3
NR
2
2
K, -
(NR
P +K S
, NR +K + NR
P +k
S
- 1, NR +k-1),
4

(K = 2., NR))
For elements DST, DSQ:
The relation binding the deformations of inflection to the field of displacement in inflection is also written under
matric form:



xx
W
NR
2 NR


K


yy = B fk U fk + B fkU fk where U F K = xk and U F K
= K represent it





K =
K = NR
2
1
1
+
xy
yk
field of displacement of inflection to the node K, so that:
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
30/54
0 J NR
0

11
K, + J
NR

12
K

,

B F K = 0
0
J NR

21
K, + J
NR

22
K
, and


0 J NR

21
K, + J
NR

22
K
J NR

,
11
K, + J
NR

12
K
,

J

11 xk
P, + J

12 xk
P

,

B F K =
J

21 yk
P, + J

22 yk
P,
.


J
11 yk
P, + J

12 yk
P, + J

21 xk
P, + J

22 xk
P,
w1

x1

y
1
The matrix of passage of the field of displacement of inflection U
= (U
F
F,)
with U F =! and
W
NR


xN

yN

1
=! with the deformations of inflection is written then:

NR
B
= (B, B
, B
,
, B
)





= (B
, B
FNR]
F
F NR
F (NR +)
F
NR
FNR]
FNR])
3 4
1
1
2
3 3
3
.
4.4
Stamp rigidity
The principle of virtual work is written in the following way: W
= W
ext.
int is still in elasticity
UT K U = F U in matric form where K is the matrix of rigidity coming from the assembly
in the total reference mark of the whole of the elementary matrices of rigidity.
4.4.1 Stamp elementary rigidity for the Q4g elements

W eint =
[E (
Hme + Hmf) +
(Hmf E + H F) +
Hct D
] S =
E
T T
T
T
T
T
T
T
(
UmBmHmBmUm +UmBmHmf B F U F +U F B F Hmf BmUm +U F B F H F B F U F
E
+ T T
U B H B U) dS
F
C
ct
C
F
=
T
T
U (B H B dS
T
T
T
T
m m m m) Um + U (B H B dS
F F
F
F
) U F + U (B H B dS) U
m m MF
F
F
E
E
E
+ UT (BT H B dS) U
F F MF m
m
E
+ UT (BT H B dS) U
T
T
T
T
T
F C ct C
F = the U.K.U
m
m
m + the U.K.U
F
F
F + the U.K.
U
m
MF
F + the U.K.
U
F
Fm
m + the U.K.U
F
C
F
E
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
31/54
with K
= K T
MF
Fm.
Um
This is still written: W E
int = (Um, U F) K

U
where
F
K
K
[
M2 NR ×2 NR]
MF [2 NR ×3N]

K


[5N ×5N] =
is the matrix of rigidity of the element.
K T
K

+ K
MF [3N ×2 NR]
F [3N ×3N]
[
C 3N ×3N]
4.4.2 Stamp elementary rigidity for elements DKT, DKQ
Since the relation = 0 is satisfied, one can write:

W eint =
E (
Hme + Hmf) +
(Hmf E + H F) dS =
E
T T
T
T
T
T
T
T
(
UmBmHmBmUm +UmBmHmf B F U F +U F B F Hmf BmUm +U B H B U) dS
F
F
F
F
F
=
E
T
T
U (B H B dS
T
T
T
T
m m m m) Um + U (B H B dS
F F
F
F
) U F + U (B H B dS
m m MF
F
) U F
E
E
E
+ UT
T
T
T
T
T
F (B F Hmf B mdS) U m
U
mKmUm
U
F K F U F
U
mKmf U F
U
F K fmU

=
+
+
+
m
E
with K
= K T
MF
Fm.
Um
K
K
[
M2 NR ×2 NR]
MF [2 NR ×3N]
This is still written: W E


int = (Um, U F) K

=
U
where K [5N×5N]
T
F
K
K
MF [3N×2N]
F [3N ×3N]
is the matrix of rigidity of the element.
4.4.3 Stamp elementary rigidity for elements DST, DSQ

W E =
E (H E + H
) +
(H E + H) +
-
TH
T
1
=

int
m
MF


MF
F

dS
ct
E
T T
T
T
T
T
T
T
(U B H B U + U B H B U + U B H B U +

U B H B U
m
m
m
m
m
m m MF F F F F MF m m F F F F F
E
+
T
T
- 1
T
T
- 1
T
T
U
B H
B U + U B H
B
+
-
B H
B
1
U + TBT
- 1
) dS =
C H
B
ct
C
F
C
ct
C
F

F

C
ct
C

C ct C F
UT (BTH B dS) U
T
T
T
T
T
T
+
(
dS)
+
(
dS)
+
(
dS




)
m
m
m
m
m
U
B H B
U
F
F
F
F
F
U
B H B
U
m
m
MF
F
F
U
B H B
U
F
F
MF
m
m
E
E
E
E
+
T
T
- 1
T
T
- 1
T
U
T
(
dS)
+
(
dS
-
) +


1
(
+
B dS U
T BT H -
)
(
B
1
dS


) =
F
B fH B
ct

U
C
F

U F B fH B
ct
C


B H
C
ct
C
F
C
ct
C
E
E
E
E
UTK U + UT K U + UTK U + UT K U + UT K U + UT K
T
+ K U
T
+ K
m
m
m
F
F
F
m
MF
F
F
Fm
m
F

F
F


F
C
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
32/54
It is also known that U
= (U
F
F,)
from where it results that:
K
= BT H B
F 11
F F F dS

;
K
K
S
F 11
F 12
K


T
F =
with: K
= B H B dS;
K T
K


F 12
F
F
F
F 12
22
S
K
= BT H B
F 22
F F F dS

.
S
K
= BT H B
MF 11
m MF F dS

;
K
= (K
K
S
MF
MF 11
MF 12) with: K
= BT H B
MF 12
m MF F dS

.
S
K
= K T
Fm
MF.
Using the fact that = WITH one deduces from it that:
F
Wint = T
T
T
T
UmK mUm + U fK F U F + UmKmf
U F + U fK Fm
Um where:
K
T
T
T
T
(
)
(
)
(
)
F = K F 11 + K + A
K

F 22 + K
With
C
+ K F 12 + K
With

+ A
K

F 12 + K.
Kmf = Kmf 11 + K
With
MF 12

Um
This is still written: W E
int = (Um, U F) K

U
where
F
K
K
[
M2 NR ×2 NR]
m
F [2N×3N]
K


[5N ×5N] =
is the elementary matrix of rigidity for an element
K T

K
MF [3N ×2 NR]
F [3N ×3N]
of plate.
4.4.4 Assembly of the elementary matrices
The principle of virtual work for the whole of the elements is written:
nbelem
W = We = T
int
int
U KU where U is the whole of the degrees of freedom of the discretized structure and
e=1
K comes from the assembly of the elementary matrices.
4.4.4.1 Degrees of freedom
The process of assembly of the elementary matrices implies that all the degrees of freedom are
expressed in the total reference mark. In the total reference mark, the degrees of freedom are three displacements
compared to the three axes of the total Cartesian reference mark and three rotations compared to these three axes.
One thus uses matrices of passage of the local reference mark to the total reference mark for each element. However one
saw previously that the degrees of freedom of the elements of plate are two displacements
in the plan of the plate, displacement except plan and two rotations. These rotations not being
exactly rotations compared to the axes of the plate since
X (X, y) = y (X, y), y (X, y) =
- X (X, y) it is necessary to hold account of it with the level of the assembly for
to reveal the good degrees of freedom xi, yi.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
33/54
4.4.4.2 Rotations
fictitious
Rotation compared to the normal with the plate is regarded as not being a degree of freedom.
To ensure compatibility between the passage of the local mode the total mode, a degree is thus added
of additional freedom room of rotation to the plate which is that corresponding to rotation by
report/ratio with the normal in the plan of the element. This implies an expansion of the blocks of dimension (5,5)
matrix of local rigidity in blocks of dimension (6,6) by adding a line and a column
correspondent with this rotation. These additional lines and these columns are a priori null. One
then carry out the passage of the matrix of local rigidity extended to the matrix of total rigidity.
In the preceding transformation, one was satisfied to add rotations compared to the normals
in the plan of the elements without modifying the deformation energy. The contribution to the energy brought by
these additional degrees of freedom is indeed null and no rigidity is associated for them.
The matrix of total rigidity thus obtained presents the risk however to be noninvertible. For
to avoid this nuisance it is allowed to allot a small rigidity to these additional degrees of freedom
on the level of the matrix of widened local rigidity. Practically, one chooses it between 10­6 and 10­3 time it
diagonal minor term of the matrix of rigidity of local inflection. The user can choose it
multiplicative coefficient COEF_RIGI_DRZ itself in AFFE_CARA_ELEM; by defect it is worth 10­5.
4.5
Stamp of mass
The terms of the matrix of mass are obtained after discretization of the variational formulation
following:
+h/2
W ac
U U
$
$
farmhouse =
$ dzdS = (U

m $$ U + $v$ v + $ W W) + (U
MF $ $ + v
X
$$ y + X U + y v)
- H/2 S
S
+ $
$
F (X X + y) dS
y
+h/2
+h/2
+h/2
with
2
m =
dz,

MF =
zdz, and

F =
Z dz

.
- H/2
- H/2
- H/2
Note:
If the plate is homogeneous or symmetrical compared to z=0 then MF =0. One considers in
continuation of the talk that it is always the case.
4.5.1 Stamp of traditional elementary mass
4.5.1.1 Element
Q4g
The discretization of displacement for this isoparametric element is:
the U.K.


NR
vk
U =
W
NR K K =1,…, NR
K
K =1
xk

yk
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
34/54
The matrix of mass, in the base where the degrees of freedom are gathered according to the directions of
translation and of rotation, has then as an expression:
M
0
0
M
m
MF
0


0
M
0
0
M

m
MF
M 0
0
M
=
0
0

m
T

M
0
0
M
MF
F
0
0
MT
0
0
M

MF
F
with: M
= NTNdS, M
= NT dS and
NR
M = NT NR
m
m
MF
MF
F
F
dS and NR = (NR
NR
1 “K).
S
S
S
4.5.1.2 Elements of the type DKT, DST
W
W

0

NR
K
2 NR





Like = NR
(,) + P
X
K
xk
xk
(,)
K where = WITH F one deduces from it that:
k=1

k=N+1



P
y
yk
yk
(,)
W
NR
(,)
0
0
W
NR
K
K




= NR

(,) NR

(,) NR
X
kxw
kxx
kxy
(,) xk.
k=1



NR



(,) NR

(,) NR
y
kyw
kyx
kyy
(,) yk
The membrane part of the elementary matrix of mass is the same one as for Q4 with k=3 with the place
from 4 in NR. The inflection part is composed of the blocks kp (kth line and pième column) following:
NR
/

kxw NR pxw + NR kyw NR pyw + m NR K NR p
F Nkxw NR pxx + Nkyw NR pyx Nkxw NR pxy + Nkyw NR pyy


F
NR kxx NR pxw + Nkyx NR pyw
NR kxx NR pxx + Nkyx NR pyx
NR kxx NR pxy + Nkyx NR pyy


NR

kxy NR pxw + NR kyy NR pyw
NR kxy NR pxx + Nkyy NR pyx
NR kxy NR pxy + Nkyy NR pyy
4.5.2 Stamp of improved elementary mass
As the arrow of a flexbeam can be represented with difficulty by an approximation
linear, one can enrich the functions by form for the terms of inflection. This approach is used
in Code_Aster for the elements of type DKT, DST and Q4 where functions of form used
in the calculation of the matrix of mass of inflection are of command 3. The interpolation for W is written as follows:
NR
W = NR
(,) W
(-)
1
+1
+ NR
(,) W
(-)
1
+2
, + NR
(,) W
K
NR
K
K
NR
K
(K -)
1 NR +3
, K

K =1
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
35/54
where the functions of form are given for the triangle and the quadrangle in the following table:
DKT, DST
DKQ, DSQ, Q4


3
4
7
3
5
6
1
4
2

8
6

1
5
2
Interpolation
= 1 - -
i=1 with 12
for W
i=1 with 9
1
NR
2
2
(,) =
1
(- 1
) (-) (2 - - - -)
NR
2
3
1
8
1 (,) = 3 - 2 + 2
1
NR
2
2
2 (,) = +/2
N2 (,) =
1
(- 1
) (- 1
) (-)
8
NR
2
3 (,) = +/2
1
NR
2
(,) =
1
(- 1
) (- 1
) (-)
NR
2
3
3
4 (,) = 3
- 2 + 2
8
NR
2
1
2
2
5 (,) = (- 1 +) -
NR 4 (,) =
1
(+ 1
) (-) (2 - - + -)
8
NR
2
6 (,) = +/2
1
NR
2
(,) = -
1
(+ 1
) (- 1
) (-)
NR
2
3
5
7 (,) = 3 - 2 + 2
8
NR
2
1
2
8 (,) = +/2
NR 6 (,) =
1
(+) 1
(-) (1 -)
8
NR
2
9 (,) = (- 1 +) -
1
NR (,) = (1+) (1+) (2
2
2
- - + +)
7
8
1
NR (,) = - (1+) (1+) (1
2
-)
8
8
1
NR (,) = - (1+) (1+) (1
2
-)
9
8
1
NR (,) = (1 -) (1+) (2
2
2
- - - +)
10
8
1
NR (,) = (1 -) (1+) (1
2
-)
11
8
1
NR (,)
12
= (1) (1) (1
2
-
+
-)
8
Functions of interpolation for the arrow of elements DKT, DST, DKQ, DSQ and Q4G,
in dynamics and modal.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
36/54
4.5.2.1 Elements of type DKT
One knows that in the approximation of one Coils-Kirchhoff has X = - W, X and y = - W, there of any point of
the element.
Because of discretization stated above one a:
NR
W = NR
(,) W
(-)
1
+1
+ (J NR
11
(-)
1
+ (,)
2
+ J NR
(,))W
21
(-)
1
+3
,
+ (J NR
K
NR
K
K
NR
K
NR
xk
12
(K -)
1 NR + (,)
2
K =1
+ J NR
(,))W
22
(K -)
1 NR +3
, yk
W, K

J
J
11
12 W, xk
since:
=


W
.
, K

J
J
21
22 W
, yk
This is still written:
NR
W = NR
(,) W
(-) +
+ NR (-) + (,) + NR
K
NR
K
K
NR
xk
(K
-) N+ (,)
1
1
1
2
1
3
yk
K =1
NR (K) 1N+ (,)
1
= NR (K) 1N+ (,)
1
where: NR (K) 1N+ (,)
2
= - J NR
11
(K -)
1 NR + (,)
2
- J NR
21
(K -)
1 NR + (,)
3
.
NR (K) 1N+ (,)
3
= - J NR
12
(K -)
1 NR + (,)
2
- J NR
22
(K -)
1 NR + (,)
3
By not taking account of the effects of inertia, the matrix of mass has the following form thus:
M
0
0
m

M = 0
Mm
0 where M = NR NR
F
m dS.
0
0
M

S
F
4.5.2.2 Elements of the DST type
It is known that for these elements one has X = X - W, X and y = y - W, there where the distortion is
constant on the element.
Like:
NR
W = NR
(,) W
(-)
1
+1
+ (J NR
11
(-)
1
+ (,)
2
+ J NR
(,))W
21
(-)
1
+3
,
+ (J NR
K
NR
K
K
NR
K
NR
xk
12
(K -)
1 NR + (,)
2
K =1
+ J NR
(,))W
22
(K -)
1 NR +3
, yk
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
37/54
one can also write:
NR
W = NR
(,) W
(-)
1
+1
+ NR (-) 1 + (,)
2
+ NR
K
NR
K
K
NR
xk
(K
-) 1N+ (,)
3
yk
K =1
+ (J11 + J) NR
12
(-) 1 + (,)
2
+ (J21 + J) NR
X
y
K
NR
X
22 y (K -)
1 NR + (,)
3
NR (K) 1N + (,)
1
= NR (K) 1N+ (,)
1
where: NR (K -) 1N + (,)
2
= - J NR
11 (K -)
1 NR + (,)
2
- J NR
21 (K -)
1 NR + (,)
3
,
NR (K) 1N + (,)
3
= - J NR
12 (K -)
1 NR + (,)
2
- J NR
22 (K -)
1 NR + (,)
3
NR
NR
(,)
1
1
= NR
(K -) NR +
(K -)
1 NR
(,)
+1
K =1
NR
NR
(,)
1
2
= NR
(K -) NR +
(K -)
1 NR
(,)
+2
K =1
NR
NR
(,)
1
3
=
NR
(K -) NR +
(K) 1N (,)
+3
K =1
w1
w1



x1


x1



y1


y



1
X
and
=
-
H 1 [B + B A
]! = T! .

ct
C
C

W
y
W
W
NR
NR




xN
xN



yN
yN
One obtains the interpolation for W then:
NR
W = NR
(,) W
(-) +
+ NR (-) + (,) + NR
K
NR
K
K
NR
xk
(K
-) N+ (,)
1
1
1
2
1
3
yk
K =1
NR (K) 1N+ (,)
1
= NR (K) 1N+ (,)
1
+
(J T (,
1 (K
11 W

-)
1 NR +)
1 + J T (2, (K
12
W

-)
1 NR +)
1) NR
(J) 1N+ (,)
2
+
(J T (,
1 (K
21 W
-)
1 NR +)
1 + J T (2, (K
22 W
-)
1 NR +)
1) NR
(J) 1N+ (,)
3
NR (K) 1N+ (,)
2
= NR (K) 1N+ (,)
2
+
where:
(J T (,
1 (K
11 W

-)
1 NR + 2) + J T (2, (K
12
W

-)
1 NR + 2)) NR
(J) 1N+ (,)
2
+
(J T (,
1 (K
21 W
-)
1 NR + 2) + J T (2, (K
22
W

-)
1 NR + 2)) NR
(J) 1N+ (,)
3
NR (K) 1N+ (,)
3
= NR (K) 1N+ (,)
3
+
(J T
1
(, (K - 1) NR + 3) + J T (2, (K - 1) NR + 3)) NR
11 W

12 W
(J 1
-) NR +2
(,) +
(J T
1
(, (K - 1) NR + 3) + J T (2, (K - 1) NR + 3)) NR
21 W
22
W

(J 1
-) NR +3
(,)
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
38/54
By not taking account of the effects of inertia, the matrix of mass has the following form thus:
M
0
0
m

M = 0
M m
0 where M =
NR dS


.


F
m
0
0
M

S
F
4.5.2.3 Elements of the Q4g type
One proceeds in the same way that for the elements of the DST type but with:
w1

x1


y1
X


= B! where B

C
C is the matrix established with [§4.3.2.1].
y
W
NR


xN

yN
4.5.2.4 Notice
One neglects in the form of the elementary matrix of mass the terms of inertia of
rotation
$
+ $
F (X X there y dS
)
because the latter are negligible [bib3] compared to the others. In
S
effect a multiplicative H2 factor/12 the dregs with the other terms and they becomes negligible for one
thickness report/ratio over characteristic length lower than 1/20.
4.5.3 Assembly of the elementary matrices of mass
The assembly of the matrices of mass follows same logic as that of the matrices of rigidity.
degrees of freedom are the same ones and one finds the processing specific to normal rotations to
plan of the plate. For modal calculations utilizing at the same time the calculation of the matrix of rigidity
and that of the matrix of mass, it is necessary to take a rigidity or a mass on the degree of rotation
normal in the plan of the plate from 103 to 106 times smaller than the diagonal minor term of the matrix
of rigidity or mass for the terms of inflection. That makes it possible to inhibit the modes being able to appear
on the additional degree of freedom of rotation around the normal in the plan of the plate. By
defect, one takes a rigidity or a mass on the degree of normal rotation in the plan of plate 105
time smaller than the diagonal minor term of the matrix of rigidity or mass for the terms of
inflection
4.5.4 Stamp of lumpée diagonal mass
The use of a matrix of lumpée diagonal mass has two advantages: it is simpler with
to implement numerically and it converges better. However the results are worse
that with the traditional diagram for which the error is minimal [bib5].
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
39/54
The matrix of mass is made diagonal, for the degrees of translation, according to a technique of
lumping near to that developed by Hinton [bib6]. The advantage of this method is of always
to produce positive coefficients of mass, contrary to other methods like the summation
terms in column [bib5]. One uses the diagonal coefficients of the elementary matrix of mass
in the total reference mark by making sure that the total mass is well represented in each of the three
directions of translation:
dS
m
m = M
S
xi
X
xi; X = NR
Mxj
j=1
dS
m
m = M
S
yi
y
yi;; y = NR
Myj
j=1
dS
m
m = M
S
zi
Z
zi; Z = NR
M
zj
j=1
where Mxi, Myi, Mzi are the diagonal terms of the elementary matrix of mass in the total reference mark,
MT
0
in each of the three directions of translation. The matrix of mass obtained: m =

0
m
m
0
0
X


with m = 0
m
T
y
0 and m = 0 are diagonal.

0
0
m
Z
Note:
The creation of a matrix of diagonal mass in the directions of rotation with a technique
analog with that of the terms of translation (coefficients I are then the coefficients
previously defined in three directions X, y and Z) is inappropriate for the modal analysis and
the results obtained are for the moment better with a matrix of reduced diagonal mass
with the only degrees of translation. The difficulty of the use of this matrix arises then in
dynamics, when it is necessary to reverse the matrix of mass. One thus considers for the moment of
to make diagonal only the terms of translation and to preserve a matrix of mass
m supplements for the terms of rotation.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
40/54
4.6
Numerical integration for elasticity
For the triangular elements DKT, DST where [H] is constant, the matrices of rigidity are obtained
exactly with three points of integration of Hammer since [B] is linear in.
Cordonnées of the points
Weight

1/6
1 = 1/6; 1 = 1/6

1/6
2 = 2/3; 2 = 1/6

1/6
3 = 1/6; 3 = 2/3
1 -
1
N
y (,) D D
=

iy (I, I)
0 0
i=1
Numerical formulas of integration on a triangle (Hammer)
For the elements quadrangle an integration of Gauss 2x2 is used.
Cordonnées of the points
Weight

1
1 = 1/
3;1 = 1/3

1
2 = 1/
3;2 = - 1/3

1
3 = - 1/
3;3 = 1/3

1
3 = - 1/
3;3 = - 1/3
1 -
1
N
y (,) D D
=

iy (I, I)
0 0
i=1
Numerical formulas of integration on a quadrangle (Gauss)
4.7
Numerical integration for plasticity
Integration on the surface of the element is supplemented by an integration on the thickness of
behavior since:
+h/2
+h/2
+h/2
H =
Hdz H
,
=
Hzdz, H =
H2
m

MF

F
Z dz

where H is the plastic matrix of behavior
- H/2
- H/2
- H/2
local.
The initial thickness is divided into NR identical layers thicknesses. There are three points of integration by
sleep. The points of integration are located in higher skin of layer, in the middle of the layer and
in lower skin of layer. For NR layers, the number of points of integration is of 2N+1. One
advise to use from 3 to 5 layers in the thickness for a number of points of integration being worth 7, 9
and 11 respectively.
For rigidity, one calculates for each layer, in plane constraints, the contribution to the matrices of
rigidity of membrane, inflection and coupling membrane-inflection. These contributions are added and
assemblies to obtain the matrix of total tangent rigidity.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
41/54
For each layer, one calculates the state of the constraints (xx, yy, xy) and the whole of the variables
interns, in the middle of the layer and in skins higher and lower of layer, from
local plastic behavior and of the local field of deformation (xx, yy, xy). The positioning of
points of integration enables us to have the rightest estimates, because not extrapolated, in skins
lower and higher of layer, where it is known that the constraints are likely to be maximum.
Cordonnées of the points
Weight

1/3
1 = - 1

4/3
2 = 0

1/3
3 = +1
1
N
y () D =

iy (I, I)
- 1
i=1
Formulate numerical integration for a layer in the thickness
Note:
One already mentioned with [§2.2.3] that the value of the coefficient of correction in shearing
transverse for the elements DST, DSQ and Q4 was obtained by identification of energies
complémenatires elastic after resolution of balance 3D. This method is not usable any more
in elastoplasticity and the choice of the coefficient of correction in transverse shearing is posed then.
Plasticity is thus not developed for these elements.
4.8
Discretization of external work
The variational formulation of work external for the elements of plate is written
:

W
= (F U + F v + F W + m + m) dS + (X
U + y
v + Z
W + X
µ X +


y
µ) ds
ext.
X
y
Z
X
X
y
y
y
S
C
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
42/54
By taking account of a linear discretization of displacements, one can write for an element:
NR

W E = (F NR
ext.
X K (,) U + F NR
K
y
K (,)

v + F NR
K
Z
K (,)


wk
K =1S
+ m NR
X
K (,)
+ m NR
xk
y
K (,)
) dS
yk
+
(NR
X K (,)
the U.K. + NR
y
K (,)
vk + NR
Z
K (,)
wk
C
+ µ NR
X
K (,
) xk + µ NR
y
K (,
)) ds
yk
NR
=

F NR


(,) dS + NR
y K (,)
X
K (,) dS + NR
X K (,) ds F NR
y
K
ds
K =1 S
C
S
C
F NR
Z K (,) dS + NR
Z K (,) ds m NR
X K (,) dS + µ NR
X K (,) ds
S
C
S
C
m NR
E
y K (,) dS + µ NR
y K (,) ds) the U.K.
S
C
NR
E
E
E
E
= K
F The U.K. = F U
K =1
The variational formulation of the work of the efforts external for the unit of the elements is written
then:
Nb elem
W = We = FU = T T
ext.
ext.
U F where U is the whole of the degrees of freedom of the structure
e=1
discretized and F comes from the assembly of the vectors forces elementary.
As for the matrices of rigidity, the process of assembly of the vectors forces elementary
imply that all the degrees of freedom are expressed in the total reference mark. In the total reference mark,
the degrees of freedom are three displacements compared to the three axes of the total Cartesian reference mark
and three rotations by report/ratio these three axes. One thus uses matrices of passage of the local reference mark
with the total reference mark for each element.
Note:
The external efforts can also be defined in the reference mark user. One then is used
stamp passage of the reference mark user towards the local reference mark of the element to have the expression of
these efforts in the local reference mark of the element and to deduce the vector from it forces elementary room
corresponding. For the assembly one passes then from the local reference mark of the element to the total reference mark.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
43/54
4.9
Taking into account of the thermal loadings
4.9.1 Thermo elasticity of the plates
The temperature is represented by the model of thermics to three fields according to [R3.11.01]:
T (X, x3) = T m (
X). P1 (x3) + T S (X). P2 (x3) + Ti



(X) .P3 (x3),
with: P (X
- H/2, +h/2
J
3): three polynomials of LAGRANGE in the thickness: ]
[:
2
X
X
P (X) = 1 - (2x/H); P
3
3
2 (x3) =
(1+2x/H
3
); P3 (x3) = - (1-2x/H
1
3
3
3
);
H
H
From the representation of the temperature above, one obtains:
· the average temperature in the thickness:
1 +h/2
1
T (X) =
T
m
S
I

,
=
4
3
3
+
+

;
- H/2 (X
X

) dx
T
X
T X
T X
H
6 (
() () ()
· the average variation in temperature in the thickness:
% (
12 + 2
)
H/
T X
=
T
;
2 - 2 (X, X X dx
3
3
3 = T S X
- T I X
H/

)
() ()
H
Thus the temperature can be written in the following way:
~
T (X, x3) = T (X) + T % (X). X/H + T
3
(X, X




3) such as:
H/2 ~
H2
T

(
/
~
X, X

3 = 0
= 0

.
- 2
);
X T X X
H
- H/
3 (
,
/

3
2
)
If the temperature is indeed closely connected in the thickness one has, ~
T = 0.
Code_Aster treats three thermo situations elastic different, where characteristics
thermo rubber bands E, depend only on the average temperature T in the thickness:
· the case where the material is thermo homogeneous isotropic rubber band in the thickness;
· the case where the plate models an orthotropic grid (concrete reinforcing steels);
· the case where the behavior of the plate is deduced from an elastic thermo homogenization,
cf [bib7].
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
44/54
For the elements of plate in thermo elasticity, the heating effects are taken into account by
the intermediary of generalized efforts, out of membrane and inflection. Thus, in the case of a plate
homogeneous, knowing the dilation coefficient, the generalized thermal efforts are defined in
to leave the plane constraints in the thickness by:
+h/2
+h/2
NR ther =
C
ther

dx3 =
.C
ref.

-

/2

- H

- H/
(T T) dx


3
2
+h/2
+h/2
M ther =
X C
ther
3

dx3 =
.x C
ref.

3
-

/2

- H

- H/
(T T) dx


3
2
V ther

= 0
Maybe in the elastic thermo case isotropic homogeneous in the thickness:
Eh
NR ther = .C
.
-
H

(T Tréf) =.
.
2 (T - T ref.)
;



1 -
H2
Eh2
M ther = .C
.
T

% =.
.
%
;
=.
0
12
(
12 1 - 2) T
V ther


The thermal constraints of origin withdrawn from the usual mechanical constraints are calculated
in three positions (sup., moy. and inf.) in the thickness:
ther


. E
=
2 (
ref.
T - T
+ T %. X/H
3
)

1 -

In the case deduced from the elastic thermo homogenization, cf [bib7], the thermal efforts generalized
are defined by the general relation, starting from the “correctors” of membrane, those of
inflection, and that of dilation udil, like averages on representative elementary volume
(cell Z):
~
NR ther =
C


.
-
+ %
.
/
3
+
, 3
+
. U
;

(T Tréf T (X) Z H T (X X

)
dil

ij (
) Cijkl kl ()
Z
Z
~
M ther =
Z .C


.
3
. -
+ %
.
/
3
+
, 3
+
. U
;

(T Tréf T (X) Z H T (X X

)
dil

ij (
) Cijkl kl ()
Z
Z
V ther

= 0
In this case when one limits oneself to the orthotropic situations without coupling inflection-membrane, one
~
neglect the role of T (X, X

3) on the corrector udil, and it is thus found that the thermal efforts which
appear to the second member have as an expression:
NR ther =. H m

.
-

;
(T
T ref.)
M ther =. H F
.T
%
; V ther

= 0



Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
45/54
One cannot however go up with the complete three-dimensional constraints: it would be necessary of
to know the “correctors” within the basic cell having been used for the determination of
coefficients of homogenized behavior.
In the elastoplastic thermo situations, or for the hulls (elements of the family
COQUE_3D), it is necessary to evaluate the three-dimensional constraints, of which constraints
thermics, in each point of integration in the thickness.
Note:
To go up with the complete three-dimensional constraints is not immediate for the hulls
multi-layer (laminated) because it is necessary to know layer by layer the state of stress; in elasticity,
this one results from the state of deformation and the behavior on the level of each layer.
4.9.2 Chaining
thermomechanics
For the resolution of chained thermomechanical problems, one must use for thermal calculation
finite elements of thermal hull [R3.11.01] whose field of temperature is recovered like
input datum of Code_Aster for mechanical calculation. It is necessary thus that there is compatibility between
thermal field given by the thermal hulls and that recovered by the mechanical plates. It
the last is defined by the knowledge of 3 fields TEMP_SUP, TEMP and TEMP_INF given in skins
lower, medium and higher of hull.
The table below indicates compatibilities between the elements of plate and the elements of hull
thermics:
Modeling
Net
Finite element
to use
Net
Finite element
Modeling
THERMIQUE
with
MECANIQUE
COQUE
QUAD4 THCOQU4
//////////////
QUAD4
MEDKQU4
DKT
///////////////
MEDSQU4
DST
/////////////
MEQ4QU4
Q4
COQUE
TRIA3
THCOTR3
///////////////
TRIA3
MEDKTR3
DKT
/////////////
MEDSTR3
DST
Note:
The nodes of the thermal elements of hulls and mechanical plates must
to correspond. The grids will be identical.
The elements of thermal hulls surface are treated like plane elements by
projection of the initial geometry on the level defined by the first 3 nodes.
The chaining with definite multi-layer materials via command DEFI_COQU_MULT [U4.23.03]
is not available for the moment.
The thermomechanical chaining is also possible if one knows by experimental measurements
variation of the field of temperature in the thickness of the structure or certain parts of
structure. In this case one works with a card of temperature defined a priori; the field of
temperature is not given any more by three values TEMP_INF, TEMP and TEMP_SUP of thermal calculation
obtained by EVOL_THER. It can be much richer and contain an arbitrary number of points
of discretization in the thickness of the hull. Operator DEFI_NAPPE allows to create such profiles
temperatures starting from the data provided by the user. These profiles are affected by
order AFFE_CARTE (cf the case-test hsns100b). It will be noted that it is not necessary for calculation
mechanics that the number of points of integration in the thickness is equal to the number of points of
discretization of the field of temperature in the thickness. The field of temperature is
automatically interpolated at the points of integration in the thickness of the elements of plates or of
hulls.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
46/54
4.9.3 Case-test
The case-tests for the thermomechanical chaining enters of the thermal elements of hulls and of
elements of plate are the hpla100e (elements DKT) and hpla100f (elements DKQ). It is about a cylinder
heavy thermoelastic hollow in uniform rotation [V7.01.100] subjected to a phenomenon of dilation
thermics where the fields of temperature are calculated with THER_LINEAIRE by a calculation
stationary.
Z
IH
Re
Interior radius IH = 19.5 mm
External radius Re = 20.5 mm
Not F
R = 20.0 mm
Thickness
H = 1.0 mm
Height
L = 10.0 mm
R
Z
NR
K
Z
Q
J
P
D
C
y
H
M
+
L
X
With
B
R
F
Thermal dilation is worth: T () - Tref () = 05
. (T + T
S
I) + 2.(T + T
S
I) (R - R)/H

with:
·
T = 0 5
. °C, T = - 0 5
. ° C, T
= 0. °C
S
I
ref.
·
T = 01
. ° C, T = 01
. ° C, T
= 0. ° C
S
I
ref.
One tests the constraints, the efforts and bending moments in L and Mr. Les results of reference are
analytical. One obtains very good results whatever the type of element considered.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
47/54
5
Establishment of the elements of plate in Code_Aster
5.1 Description
:
These elements (of names MEDKTR3, MEDSTR3, MEDKQU4, MEDSQU4 and MEQ4QU4) rest on
meshs TRIA3 and plane QUAD4. These elements are not exact with the nodes and it is necessary to net with
several elements to obtain correct results.
5.2
Introduced use and developments:
These elements are used in the following way:
· AFFE_MODELE (MODELING: “DKT”.) for the triangle and the quadrangle of the type DKT
· AFFE_MODELE (MODELING: “DST”.) for the triangle and the quadrangle of the DST type
· AFFE_MODELE (MODELING: “Q4G”.)
One calls upon routine INI079 for the position of the points of Hammer and Gauss on the surface of
plate and the corresponding weights.
· AFFE_CARA_ELEM (HULL:(EPAISSEUR:“EP”
ANGL_REP: ('' '')
COEF_RIGI_DRZ: “CTOR”)
To make postprocessings (forced, generalized efforts,…) in a reference mark chosen by the user
who is not the local reference mark of the element, one gives a direction of reference D defined by two angles
nautical in the total reference mark. The projection of this direction of reference in the field of the fixed plate
a X1 direction of reference. The normal in the plan into fixed one second and the vector product of both
vectors previously definite makes it possible to define the local trihedron in which the efforts will be expressed
generalized and constraints. The user will have to take care that the selected reference axis is not found
not parallel with the normal of certain elements of plate of the model. By defect this direction of
reference is axis X of the total reference mark of definition of the grid.
Value CTOR corresponds to coefficent that the user can introduce for the processing of the terms of
rigidity and of mass according to normal rotation in the plan of the plate. This coefficient must be sufficiently
small not to disturb the energy balance of the element and not too small so that matrices of
rigidity and of mass are invertible. A value of 10­5 is put by defect.
· ELAS: (E:NAKED Young: ALPHA:. RHO:. )
For an elastic thermo behavior isotropic homogeneous in the thickness one uses key word ELAS
in DEFI_MATERIAU where the coefficients E are defined, Young modulus, Poisson's ratio,
thermal dilation coefficient and RHO density.
· ELAS_ORTH (_FO): (
E_L:ygl. E_T:ygt. .G_LT:glt. G_TZ:gtz. .NU_LT:nult.
ALPHA_L:L. ALPHA_T:T.)
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
48/54
For a thermo behavior elastic orthotropic whose axes of orthotropism are L, T and Z with
isotropy of axis L (fibers in the direction L coated by a matrix, for example) the seven should be given
independent coefficients ygl, longitudinal modulus Young, ygt, transverse modulus Young, glt,
modulus of rigidity in the plan LT, gtz, modulus of rigidity in plan TZ nult, coefficient of
Poisson in plan LT and the dilation coefficients thermal L and T for thermal dilation
longitudinal and transverse, respectively. The orthotropic behavior is not available
that associated the key word DEFI_COQU_MULT which makes it possible to define a multi-layer composite hull.
For only one orthotropic material, one will thus use DEFI_COQU_MULT with only one layer. If one
wish to use ELAS_ORTH with transverse shearing, it is necessarily necessary to employ
DST modeling. If modeling DKT is used, transverse shearing is not taken into account.
· ELAS_COQUE (_FO): (
MEMB_L:C1111. MEMB_LT:C1122. MEMB_T:C2222. MEMB_G_LT:C1212.
FLEX_L:D1111. FLEX_LT:D1122. FLEX_T:D2222. FLEX_G_LT:D1212.
CISA_L:G11…. CISA_T:G22…. ALPHA:. RHO:. )
This behavior was added in DEFI_MATERIAU to take into account matrices of rigidity
nonproportional in rigidity and inflection, obtained by homogenization of a multi-layer material.
The coefficients of the matrices of rigidity are then introduced with the hand by the user into the reference mark
user defined by key word ANGL_REP. The thickness given in AFFE_CARA_ELEM is only
used with the density defined by RHO. ALPHA is thermal dilation. If one wishes
to use ELAS_COQUE with transverse shearing, it is necessarily necessary to employ modeling
DST. If modeling DKT is used, transverse shearing is not taken into account.
· DEFI_COQU_MULT: (COUCHE: EPAISSEUR:“EP”
MATER: “material”
ORIENTATION: (''))
This key word makes it possible to define a multi-layer composite hull on the basis of the sub-base towards
roadbase starting from its characteristics sleeps by layer, thickness, type of material
constitutive and orientation of fibers compared to a reference axis. The type of constitutive material is
product by operator DEFI_MATERIAU under key word ELAS_ORTH. is the angle of the first
direction of orthotropism (longitudinal direction or direction of fibers) in the tangent plan with the element by report/ratio
with the first direction of the reference mark of reference defined by ANGL_REP. By defect is null, if not it must be
provided in degrees and must lie between ­ 90º and + 90º.
· AFFE_CHAR_MECA (DDL_IMPO: (
DX:. DY:. DZ:. DRX:. DRY:. DRZ:. DDL of plate in the reference mark
total.
FORCE_COQUE: (FX:. FY:. FZ:. MX:. MY:. MZ:. ) They are the efforts
surface on elements of plate. These efforts can be given in the total reference mark or
in the reference mark user defined by ANGL_REP.
· FORCE_NODALE: (FX:. FY:. FZ:. MX:. MY:. MZ:. ) It is about
efforts of hull in the total reference mark.
5.3
Calculation in linear elasticity:
The matrix of rigidity and the matrix of mass (respectively options RIGI_MECA and MASS_MECA)
are integrated numerically. It is not checked if the mesh is plane or not. Calculation holds account
owing to the fact that the terms corresponding to the DDL of plate are expressed in the local reference mark of
the element. A matrix of passage makes it possible to pass from the local DDL to the total DDL.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
49/54
Elementary calculations (CALC_ELEM) currently available correspond to the options:
· EPSI_ELNO_DEPL and SIGM_ELNO_DEPL which provide the strains and the stresses
with the nodes in the reference mark user of the element in lower skin, with semi thickness and in
higher skin of plate, the position being specified by the user. These values are stored
in the following way: 6 components of strain or stresses:
· EPXX EPYY EPZZ EPXY EPXZ EPYZ or SIXX SIYY SIZZ SIXY SIXZ SIYZ
· DEGE_ELNO_DEPL: who gives the deformations generalized by `element to the nodes to
to leave displacements in the reference mark user: EXX, EYY, EXY, KXX, KYY, KXY,
GAX, GAY.
· EFGE_ELNO_DEPL: who gives the efforts generalize by element with the nodes from
displacements: NXX, NYY, NXY, MXX, MYY, MXY, QX, QY.
· SIEF_ELGA_DEPL: who gives the efforts generalize by element at the points of Gauss with
to leave displacements: NXX, NYY, NXY, MXX, MYY, MXY, QX, QY.
· EPOT_ELEM_DEPL: who gives the elastic energy of deformation per element from
displacements.
· ECIN_ELEM_DEPL: who gives the kinetic energy by element.
Finally one calculates also option FORC_NODA of calculation of the nodal forces for operator CALC_NO.
5.4
Plastic design
The matrix of rigidity is there too integrated numerically. One calls upon the option of calculation
STAT_NON_LINE in which one defines in the level of the nonlinear behavior the number of
layers to be used for numerical integration. All laws of plane constraints available in
Code_Aster can be used.
STAT_NON_LINE (….
COMP_INCR: (RELATION:''
COQUE_NCOU:“A NUMBER OF LAYERS”)
….)
Elementary calculations (CALC_ELEM) currently available correspond to the options:
· EPSI_ELNO_DEPL which provides the deformations by element to the nodes in the reference mark
user starting from displacements, in lower skin, with semi thickness and in higher skin
of plate.
· SIGM_ELNO_COQU which makes it possible to obtain the stress field in the thickness by element
with the nodes for a given layer and a position requested (in lower skin, in the medium
or in higher skin of layer). These values are given in the reference mark user.
· SIEF_ELNO_ELGA which makes it possible to obtain the efforts generalized by element with the nodes in
the reference mark user.
· VARI_ELNO_ELGA which calculates the field of internal variables and the constraints by element
with the nodes for all the layers, in the local reference mark of the element.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
50/54
6 Conclusion
The finite elements of plate plans which we describe here are used in the structural analyzes
thin, in small displacements and deformations, of which the thickness report/ratio over characteristic length
is lower than 1/10. As these elements are plane, they do not take into account the curvature of
structures, and it is necessary to refine the grids if this one would be important.
It is elements for which strains and stresses in the plan of the element
vary linearly with the thickness of the plate. Moreover, the distortion associated with shearing
transverse is constant in the thickness of the element. Two types of elements of plate exist:
elements DKT, DKQ for which the transverse distortion is null and the elements DST, DSQ and Q4
for which it remains constant and nonnull in the thickness. One advises to use the second type
elements when the structure to be netted has a thickness report/ratio over characteristic length included/understood
between 1/20 and 1/10 and first in the remainder of the cases. When the transverse distortion is not
null, the elements of DST plate, DSQ and Q4 do not satisfy the equilibrium conditions 3D and them
boundary conditions on the nullity of stresses shear transverse on the faces higher
and lower of plate, compatible with a constant transverse distortion in the thickness of
plate. It results from it thus that on the level from the behavior a coefficient from 5/6 for a plate
homogeneous the usual relation between the constraints and the distortion transverses corrects in order to
to ensure the equality enters energies of shearing of the model 3D and the model of plate to distortion
constant. In this case, the arrow W has as an interpretation average transverse displacement in
the thickness of the plate.
The non-linear behaviors in plane constraints are available for the elements of plate
DKT and DKQ only. Indeed the rigorous taking into account of a constant transverse shearing
not no one on the thickness and the determination of the correction associated on rigidity with shearing by
report/ratio with a model satisfying the equilibrium conditions and the boundary conditions are not
possible and thus the use of the DST elements, DSQ and rigorously impossible Q4 return in
plasticity.
Elements corresponding to the machine elements exist in thermics; chainings
thermo mechanics is thus available except, for the moment, laminated materials.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
51/54
7 Bibliography
[1]
J.L. Batoz, G.Dhatt, “Modélisation of the structures by finite elements: beams and plates”,
Hermès, Paris, 1992.
[2]
D. Bui, “shearing in the plates and hulls: modeling and calculation”, Note HI-
71/7784, 1992.
[3]
J.G. Ren, “A new theory off laminated punt”, Composite Science and Technology, Vol.26,
p.225-239,1986.
[4]
T.A. Rock, E. Hinton, “A finite element method for the free vibration off punts Al for
transverse shear deformation “, Computers and Structures, Vol.6, p.37-44,1976.
[5]
T.J.R. Hughes, “The finite element method”, Prentice Hall, 1987.
[6]
E. Hinton, T. Rock and O.C. Zienkiewicz, “A notes one Mass Lumping and Related Processes in
the Finite Element Method “, Earthquake Engineering and Structural Dynamics, Vol4, p. 245-
249, 1976.
[7]
F. Voldoire “Modélisation by thermal and thermo homogenization elastic of components
mean mechanics ", CR MMN/97/091.
[8]
R3.11.01 P. MASSIN, F. VOLDOIRE, S. ANDRIEUX “Modèle of thermics for the hulls
thin ", Manuel of reference of Code_Aster.
[9]
V7.01.100 F. VOLDOIRE “thermoelastic hollow Cylindre”, Manuel of validation of
Code_Aster.
[10]
A.K. Noor, W.S. Burton, “Assessment off shear deformation theories for multilayered
composite punts “, ASME, Applied Mechanics Review, Vol.42, N°1, p.1-13,1989.
[11]
A.K. Noor, W.S. Burton, J.M. Peters “Assessment off computational models for multilayered
composite cylinders “in Analytical and Computational Models off Shells, Noor and Al Eds,
ASME, CED - Vol.3, P.419-442,1989.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
52/54
Appendix 1 Plaques orthotropic
For an orthotropic material like that represented below, composed for example of fibers of direction L
coated with a matrix, whose axes of orthotropism are L, T and Z with isotropy of axis L, the expression for
matrices H and H in the reference mark of orthotropism previously definite becomes:
H
H
0
L
LT

GLZ
0
HL = H
H
LT
TT
0 and HL =

0
GTZ

0
0
GLT
E
E
E
H
L
=
; H
T
L
L
TT =
G
=
1 -
LZ
LTTL
1 - LTTL
2 1
(+ LZ)
with
and
.
E
E
E
H
T LT
L TL
T
LT =
=
G
=
1 -
TZ
LTTL
1 - LTTL
2 1
(+ TZ)
The knowledge of the five independent coefficients E, E, G
, G
L
T
LT
TZ, LT is sufficient to determine
coefficients of the matrices H and H since:
AND

LT
TL =
and G
= G.
E
LZ
LT
L
If one indicates by the angle between the reference mark of orthotropism and the principal axis of the reference mark defined by the user in
means of ANGL_REP one establishes that:
H = TTH T
T
1
L 1 and H
= T H T

2
L
2
C2
S 2
CS



C
S
with: T
2
2
1 = S
C
- CS and T2 =

= cos, = sin
X L
-
where C
S
and = (,) like
S C

- CS CS C2 - S2
2
2


indicated on the figure below.
Z
L

X
T
In the case of forced initial of thermal origin, we have moreover:
L T




T
HT = - T H
1
L T T


0
where L and T are the dilation coefficients thermal in the directions L and T and T the variation of
temperature.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
53/54
Appendix 2 Facteurs of transverse correction of shearing for
orthotropic or laminated plates
The Hct matrix is defined so that the surface density of transverse energy of shearing obtained in
the case of the three-dimensional distribution of the constraints resulting from the resolution of balance is equal to that
model of plate based on the assumptions of Reissner, for a behavior in pure bending. One must
thus to find Hct such as:
1 +h/2
1
1

+h/2
xz
H-1
- 1
= TH T


ct
= Hct with =
and T =
dz = H

.
2
2
2

ct
yz
- H/2
- H/2
To obtain Hct one uses the distribution of following Z obtained starting from the resolution of the equations
of balance 3D without external couples:
Z
Z
= - (
with xz =
= 0 for z=±h/2.
, +
) D;
,
= - (, +
) D


xz
xx X xy y yz
xy X yy, y
yz
- H/2
- H/2
If there is no coupling membrane inflection (symmetry compared to z=0), constraints in
plan of element xx, yy, xy in the case of have as an expression a behavior of pure inflection:
= zA (Z) M with A Z = H Z H 1
()
() - F.
If (
H Z) and Hf do not depend on X and y one can determine Hct. Indeed:
Mxx, X - Mxy, y


Tx Mxx, X + Mxy, y
M xy, X - M


yy, y
(Z) = D (Z) T


1
+ D (Z)
2
where T =
=
and =
T

M

y
Mxy, X + M yy, y

yy, X


M xx, y

like:
Z
With
With
With
With
11 +
33
13 +
32
D = -

D,
1
2 A
With
With
With
H/2
31 +
23
22 +
33
-
Z
With
With
With
With
With
With
11 -
33
13 -
2
2
32
12
31
D = -
D.
2
2 A
With
With
With
With
With
H/2
31 -
23
33 -
2
2
22
32
21
-
+h/2
C =
DTH 1
- D
11
1 1dz;
- H/2
1 +h/2
1
+h/2
1

T C
C
11
12
T
It results from it that
-
H
T
- 1
=
T
with: C =
D H D

dz;
2
2 C
C
12
1

2
12
22
- H/2
- H/2
+h/2
C
=
DTH 1
- D
22
2 2dz
- H/2
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Code_Aster ®
Version
5.0
Titrate:
Elements of plate DKT, DST, DKQ, DSQ, Q4g
Date:
12/01/01
Author (S):
P. MASSIN
Key:
R3.07.03-A
Page:
54/54
1 +h/2
1
As in addition
- 1
-
H
1
1
-
=
TH T one proposes to take H = C to satisfy them as well as possible
2
2
ct
ct
11
- H/2
two equations whatever T and.
+h/2
By comparing Hct thus calculated with H
H
ct = dz one reveals the coefficients of correction of
- H/2
following transverse shearing: K = H11/H 11;K
= H12/H 12;K = H22/H 22
1
ct
ct
12
ct
ct
2
ct
ct.
For a homogeneous, isotropic or anisotropic plate, one finds as follows: Hct =kh H with k=5/6.
Note:
This method is valid only when the composite plate is symmetrical compared to z=0.
· For a multi-layer material, one establishes that:
NR
i-1
I
H
1
1
1
C
= I (H AT
2
2
p p
- Z AT) H-1
I
(

H With
p p
- Z A
11
p
I
p
I
I) +
4
2
2
i=1
p=1
p=1
i-1
I
1
1
3
3
-
1
1
2
2
-
(zi+1 - Z) [ATH 1
I
(

H With
p p
- Z A
I
) + (H AT
p p
- Z AT) H 1A
I
p
I
p
I
I

I]
24
2
2
p=1
p=1
1
+
(z5
5
T
- 1
i+1 - zi) A H
With
80
I

I
1
WITH + A
WITH + A
11
33
13
32
where: H = Z +1 - Z, = (Z +1 + Z
I
I
I
I
I
I) and A

for
2
I represents the matrix A + A
WITH + A
31
23
22
33
layer I.
· Validity of the choice H
= C-1
ct
11 can be examined a posteriori when one has an estimate of
solution (fields of displacements and plane constraints, in particular). One can then estimate
the variation enters the two estimates on energy. A step of calculation in two stages for
multi-layer plates and hulls (with Hct diagonal and two coefficients k1 and k2) was besides
developed by Noor and Burton [bib11] [bib12].
· In the case of an isotropic or anisotropic homogeneous plate the equality between two energies is
satisfied in a strict sense since D2 = 0. The choice makes above is then valid and no examination has
posteriori is not necessary.
Handbook of Référence
R3.07 booklet: Machine elements on average surface
HI-75/01/001/A

Outline document