REPUBLIQUE DU SENEGAL

UNIVERSITE CHEIKH ANTA DIOP

ECOLE SUPERIEURE POLYTECHNIQUE

Centre de Thiès
Département Génie Civil

Projet de fin d'études

En vue de l'obtention du Diplôme d'Ingénieur de Conception

Titre: MODELISATION DES COQUES DE REVOLUTION A L'AIDE

DES ELEMENTS FINIS .APPLiCATION : CAS DES RESERVOIRS
D'ASSAINISSEMENT

Auteur: Mamadou DIOP

Directeur interne: Dr Moustapha NDIAYE

Année: 2002 - 2003

\~J
-',
Projet de Fin d'Etudes Ecole Supérieure Polytechnique

DEDICACES

Je dédie ce projet à mon défunt père, à ma mère, à mes

frères et sœurs, à mes amis, à toute la promotion et à tout
ce qui me sont chérs.

Mamadou DIOP l Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

REMERCIEMENTS

Mes remerciements vont à l'endroit de tout ceux qui ont participé de prés ou de loin à la
réussite de ce projet . Mais une mention spéciale est réservée à Monsieur Moustapha
NDlAYE qui n'a ménagé aucun effort pour la réussite de ce travail. Car tout au long de
l'étude il a montré une disponibilité sans faille et un enthousiasme sans commune mesure
pour mener à bien ce projet.

Mamadou DIOP II Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

SOMMAIRE

Dans le cadre de notre projet de Fin d'Etude intitulé modélisation des coques de révolution
par éléments finis, nous avons eu à définir la méthode des éléments finis qui a été
développée pour la première fois dans les années 60 .Nous avons aussi établi les relations
géométriques d'une coque d'épaisseur variable avant de le spécifier sur les coques
axisymétriques, objet de notre étude. Les relations cinématiques et mécaniques basées sur
les hypothèses de Reissner/Mindlin et de Kirchhoff pour les comportements
axisymétriques ont été aussi présentés. Mais comme toute étude par la méthode des
éléments finis doit comporter les conditions aux limites , c'est pourquoi elles ont été
formulées ainsi que la formulation variationnelle.
Dans la deuxième partie nous avons étudié trois éléments dont les deux à savoir l'élément
CAXI_K et CAXCQ sont basés sur la théorie de Reissner/Mindlin et l'élément CAXCL
sur la théorie de Kirchhoff. Ainsi pour chaque élément on a isolé sa matrice de rigidité, de
masse et les charges équivalentes aux nœuds.
Dans la dernière partie on a dimensionné un réservoir à l'aide du logiciel Microfe afin de
mesurer la puissance de nos éléments par rapport à ce que nous offre le marché.

Mamadou DIOP III Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Table des matières

Introduction 1

1. Théorie 2
1.1 Définition de la méthode des éléments finis 2
1.1.1 Qu'est que c'est la méthode des éléments finis 2
1.1.2 Comment cette méthode fonctionne 2
[Link] Discrétisation géométrique 2
[Link] Construction des approximations nodales 2
[Link] Calcul des matrices élémentaires 3
1.1.3 Assemblage 2
1.1.4 Prises en compte des conditions aux limites 3
1.1.5 Evaluations des grandeurs élémentaires 3
1.2 Description géométrique 3
1.2.1 Aspects généraux - Définition 3
1.2.2 Géométrie de la surface 3
1.2.3 Vecteurs de base tç, t 11 et n relatif au trièdre fondamental .4

1.2.4 Vecteurs de la base duale auxiliaire PI; , P'Il et n 6

1.2.5. Courbure et torsion 8

[Link] Evaluation de la courbure K 9
[Link] Evaluation de la torsion T Il
1.2.6 . Courbure et torsion dans les directions tl; et t'Il 12

1.2.7 Spécialisation dans le cas d'un repère cartésien et cylindrique 13

1.3 Relations cinématiques 15
1.4 Conditions d'équilibre 17
1.5 Loi constitutive 19
1.5.1 Cas anisotrope 20
1.5.2 Cas orthotrope 21
1.5.3 Cas isotrope 21
1.6 Conditions de bord 22
1.7 Formulation variationnelle 23
1.7.1 Equivalence des travaux virtuels 23

Mamadou DIOP N Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

1.8 Fonctions d'interpolation 25

1.8.1 Elément tronconique sans contrainte transversale (CAXI_ K ) 25
[Link] Matrice de Rigidité 25
[Link] Matrice de masse 29
[Link] Forces équivalentes aux nœuds 29
1.8.2 Elément tronconique linéaire (CAXC L ) 30
[Link] Matrice de Rigidité 34
[Link] Matrice de masse ,.. , 34
[Link] Forces équivalentes aux nœuds 35
1.8.3 Elément isoparamétrique quadratique (CAXCQ) 36
[Link] Matrice de Rigidité 36.
[Link] Matrice de masse 37
[Link] Forces équivalentes aux nœuds 37
2 Application pratique 38
Conclusion , 41

Mamadou DIOP v Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Listes des figures et tableaux

Figure 1.1: Représentation paramétrique de la surface moyenne d'une coque: P 4

Figure 1.2: Angle œentre les directions tl; et t11 du repère local: P5
Figure 1.3: vecteurs de la base duale: P7
Figure 1.4: Changement de direction du vecteur normal; P8
Figure 1.5: Repère cylindrique et cartésien: PlI
Figure 1.6: Tronc de cône: P13
Figure 1.7: Rayon de courbure R, et Re: P 14
Figure 1.8: Rayon de courbure Rs et Re: P14
Figure 1.9 : : Rayon de courbure R s et Re: P14
Figure 1.10: cas d'une coque à épaisseur variable: P15
Figure 1.1 1 : Portion de coque: P23
Figure 1.1 2 : Moment Me et effort résultant de membrane Ne : P24
Figure 1.1 3 : Moment Ms et effort résultant de membrane Ns et effort tranchant Ts : P25
Figure 1.1 4 : Elément de référence : P26
Figure 1.1 5 : Elément tronconique CAXI_L : P26
Figure 1.1 6 : Géométrie d'un élément génératrice d'un tronc de cône: P30
Figure 1.1 7: Elément tronconique linéaire CAXI_L : P30
Figure 1.18: Elément quadratique CAXI_Q : P 35
Figure 1.19: Photo d'un bassin de décantation: P 38

Tableau 1: Conditions aux limites P 22

Mamadou DIOP VI Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Liste des symboles et abréviations

f intégrale

{} matrice colonne

[] matrice
T
[M] = transposée de la matrice M

W, W e formes intégrales globales et élémentaires

8 symbole de calcul des variations
TI fonctionnelle d'énergie
[B] matrice reliant les déformations aux variables nodales

[H] matrice de comportement élastique

[Hm], [Hmf] et [Hf] matrices de comportement homogénéisées de membrane,

flexion et couplage membrane-flexion.

He = rigidité de cisaillement

E, El ,E2 = module de Young

CT = cisaillement transversal
[J] = jacobienne de transformation

[k lae ,[k l,ob =matrices de rigidité globale élémentaire en repère local et global

Mamadou DIOP VII Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

z, z' = coordonnée suivant n, suivant nç

e = vecteur rotation
pm, p r, prof = inerties homogénéisées de membrane, flexion et couplage membrane-
flexion

Ms Me = moments de flexion

NsNe = efforts de membrane

R , (Rs , Re) =rayon de courbure (d'un méridien, d'un parallèle)

Mamadou DIOP VIII Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

N oms des principaux éléments utilisés.

CAXI L élément de coque axisymétrique linéaire de type Mindlin
CAXI K élément de coque axisymétrique de type Kirchhoff
CAXCQ élément de coque axisymétrique quadratique de type Mindlin

Mamadou DIOP IX Année académique 2002/2(~H

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Introduction
Lorsque correctement appuyée, une coque devient un élément structural de grande
efficacité (considérant le rapport poids/portée) qui transmet les charges par action
membranaire.
En effet, à l'exception des zones adjacentes aux discontinuités de géométrie ou de
chargement (bordures, ouvertures, changement brusque d'épaisseur, charges concentrées
ou sectorielles) , les effets flexionnels sont négligeables. Les propriétés structurales
intéressantes des coques ont été mises à contribution depuis l'antiquité dans la construction
navale.
De nos jours, les coques sont souvent utilisées en Génie civil dans la construction des
voûtes cylindriques ou sphériques en maçonnerie. Avec l'avènement du béton armé, les
formes et les applications se sont diversifiées vers la construction des barrages voûtes, des
châteaux d'eau et réservoirs et des toitures de grandes portées. Les structures en forme de
coque constituent les éléments structuraux de base dans la construction aéronautique et
aérospatiale. Le comportement structural des plaques et coques est décrit par des équations
aux dérivées partielles complexes. Des solutions de ces équations ne sont obtenues , de
façon précise, que pour quelques rares cas, souvent d'ordre académique. Les difficultés
rencontrées dans la résolution de tels cas de structures sont contournées avec l'usage de la
méthode des éléments finis. Cette méthode donne une solution certes approximative, mais
suffisamment précise pour des applications pratiques.
Bien qu'étant élaborée depuis le début des années soixante, la recherche est en cours
encore aujourd'hui, les objectifs s'inscrivant dans l'amélioration des solutions obtenues.
On distingue deux catégories de coques: les coques de révolution et les coques à
géométrie quelconque. Leur analyse repose sur la description de leur surface avec un
système de coordonnées le plus approprié, la détermination des paramètres géométriques
déterminants et l'écriture des équations aux dérivées partielles fondamentales qui décrivent
leur comportement mécanique.

Mamadou DIOP 1 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

1 Théorie générale

1.1 Définition de la méthode des éléments finis

Dans cette partie où nous essayons de voir de façon générale la méthode des éléments finis,
il s'avère important afin d'éviter toute mauvaise compréhension de définir les principes de
base de la méthode des éléments finis

1.1.1 Qu'est que c'est la méthode des éléments finis

La méthode des éléments finis est une technique d'analyse numérique qui permet d'obtenir
des solutions approchées dans une large variété de problèmes d'ingénieries .Tout au début
cette méthode était développée pour étudier les structures complexes, mais par la suite
elle a était étendue et appliquée dans le champs des mécaniques continues. Aujourd'hui on
l'utilise de plus en plus pour la résolution de beaucoup de problèmes.

1.1.2 Comment cette méthode fonctionne

La procédure de discrétisation en éléments finis réduit le problème à un nombre finis
d'inconnus en divisant la région de solution en éléments et en exprimant le champs
variable d'inconnu par des fonctions d'approximation Les fonctions d'interpolation sont
définies en termes de valeurs du champs de variable à des points spécifiques appelés
nœuds. Les valeurs nodales du champs variable et les fonctions d'interpolation des divers
éléments définissent de manière complète le comportement de ce champs.
La nature de la solution et le degré de d'approximation ne dépend pas de la taille et du
nombre d'éléments utilisés mais plutôt des fonctions d'interpolation choisies

[Link] Discrétisation géométrique

La première étape consiste à diviser le continuum ou la région de solution en éléments.
Une variété d'élément doit être utilisé et les différent d'éléments doivent être employés
dans la même région de solution. Néanmoins lorsqu'on analyse une structure élastique
composée de différents type d'éléments telles les poutres et plaques, il est non seulement
souhaité, mais plutôt nécessaire d'utilisé différents éléments dans la même solution.

[Link] Construction des approximations nodales

L'étape suivante est l'affectation de nœuds à chaque éléments et de choisir les fonctions
d'interpolation pour représenter la variation du champs variable à travers l'élément. Ce

Mamadou DIOP 2 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

champs peut être un scalaire, un vecteur, ou un tenseur. Souvent les polynômes sont
choisis comme fonction d'interpolation pour le champs variable à cause de sa facilité
d'intégration et de différenciation. Le degré du polynôme choisi dépend du nombre de
nœuds affecté à l'élément, la nature et le nombre d'inconnus à chaque nœuds, ct à une
certaine exigence de continuité au niveau des nœuds et le long des limites de l'élément.

[Link] Calcul des matrices élémentaires

Une fois le modèle d'élément fini établi, nous sommes prêt pour déterminer l'équation
matricielle définissant les propriétés des divers éléments composant l'ensemble. Pour cela
on utilise l'approche variationnelle.

1.1.3 Assemblage
Pour trouver les propriétés du système modélisé par un réseau d'éléments nous devons
assembler tous les propriétés des éléments. En d'autres termes nous combinons les
équations matricielles celles définissant les comportements de l'ensemble du système

1.1.4 Prises en compte des conditions aux limites

Avant que les systèmes d'équation soient prêts pour la solution ils doivent être modifiés en
tenant compte des conditions aux frontières du problème. A ce stade on impose les valeurs
nodales connues.

1.1.5 Evaluations des grandeurs élémentaires

La procédure d'assemblage donne un ensemble d'équation simultané que nous résolvons
pour obtenir les valeurs nodales inconnues du problèmes.

1.2 Description géométrique

1.2.1 Aspects généraux - Définitions

Une coque de révolution est définie géométriquement par une surface moyenne A et une
épaisseur h =2t .La surface moyenne est obtenue par rotation d'une courbe plane ou
méridien autour d'un axe de révolution Z.L'épaisseur h, supposée petite par rapport aux
autres dimensions de la coque ( longueur du méridien, circonférence) est définie suivant la
direction z normale à la surface moyenne. Suivant l'ordre de grandeur de h par rapport aux
autres aux autres dimensions on introduit parfois l'adjectifmince ou épais' aux coques. Ce

Mamadou DIOP 3 Année académique 2002/2:>OJ

Projet de Fin d'Etudes Ecole Supérieure-Polytechnique

qualificatif n'implique pas seulement une caractéristique géométrique mais sous -entend
également un rôle particulier des déformations dites de cisaillement transversal (CT ).

1.2.2 Géométrie de surface

On définit la géométrie de la coque d'épaisseur t en coordonnées paramétriques en

exprimant les coordonnées cartésiennes du vecteur position üM =X i + Y j + Z k de la

surface moyenne comme:

X = f(S~ u)
y = g(~, YI) (1.1)

Z = h(é. 11)

1.2.3 Vecteurs de base tl;, tT) et n relatif au trièdre fondamental

Les valeurs de ç = constante et " = constante sont les courbes de surface CI; et CT). Un point
M de la surface moyenne peut toujours être défini à l'intersection d'une courbe CI; et d'une

courbe CT)' On appelle trièdre fondamental relatif à cette surface, le trièdre formé par les
trois vecteurs unitaires
tl; tangent à CI; dirigé vers les ç croissants
tT) tangent à CT) dirigé vers les" croissants et

n normal à la surface moyenne et formant avec tl; et un tT) trièdre droit.

OM
t
"
=--"
B
et (1.2 a)

Mamadou DIOP 4 Année académique 200212003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

OM ç = (X,ç Y'ç Z'ç)' OMT] = (X'T] Y'T] Z'T]) (1.2 b )

OM ç ® OMT] = (Y,ç Z'T] - Z'ç Y'T]: Z'ç X'T] - X'ç Z'T] X'J: Y, -
':> T]
Y,~; X,T] )
avec

B= IOM'lll et (1.2 c)

Comme montré à la figure suivante, en général les vecteurs t et tT] ne sont pas
ç

perpendiculaires. Ils font entre eux un angle 0) tel que

Figure 1.2: Angle 0) entre les directions

t ç et t ll du repère local

H/AB

(1.3 a)

cos œ = tç • tT] = (X'çXT]' + Y'çY'T] + Z'çZ'T])/AB

Suite à une variation de dç et dn, le vecteur position varie comme

(1.3 b)

et le carré de la longueur de l'arc correspondant sur la surface moyenne est

2
ds =[Link] = A 2 dç 2 + 2 AB cos 0) dçdT] + B 2 dn 2 = 1 (1.4 a)

Mamadou DIOP 5 Année académique 2U02/2003

Projet de Fin d'Etudes Ecole Supérieure flolytedmique

Si tl; et t sont orthogonaux H = 0 et œ = rr/2 comme dans le cas des coques à symétrie
11
de révolution. 1 est la première forme fondamentale de la surface. Développer autrement,

A2
1 = [dl; d11 ] AB cos Cûl {dl;} (lA b)
[ B
2 d11
BAcos œ

Le tenseur métrique de la surface est

AB cos Cûl (lAc)

2
B

La matrice de transformation géométrique  ( repère local ~ repère global) est construite

à partir de la relation:

les composantes du vecteur fi dans le repère globaL

X'I;IA X, lB a
11
Â= Y'I;IA Y, lB b (1.5 b)
11
Z'I;IA Z, lB c
11

Les vecteurs étant unitaires 1Â 1= sin 00, l'aire du parallélogramme que forment tl; et

t .
11

Mamadou DIOP 6 AImée académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

1.2.4 Vecteurs de la base duale auxiliaire PS ' P'll et n

Les vecteurs tangents t et t'Il sont unitaires mais souvent non orthogonaux. Par
s

commodité d'autres vecteurs tangents Ps ' P'Il sont construites telle que P soit
s

perpendiculaire à t et p à t): (voir figure suivante)

'Il 'Il ~

Pll

Figure 1.3: vecteurs de la base duale

Suivant la base (té. t ) les vecteurs P): ,p s'expriment comme

_ 'Il ~ 'Il

Ps == PSSt s + PS'll t'Il (1.6 a)

P'Il -P t +p t
- 'ilS S 'Il'll 'Il
(1.6 b)

t • P): == (t t n) == 0
'Il ~ 'Il 'Il

(1.6 c)

== p ): + cos œP
'Il~ 'Il'Il

t • P == (t t):n) == sin œ == p ): cos œ + p

11 11 11 ~ 'Il~ 1111
De (1.6 c) on tire

Mamadcu DIOP 7 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure flolytechnique

p _p _ 1
SS - l'Pl - sin CD
(1.6 d)
p = P =-cot CD
ST] T]S
Soit
_1_ -cot a
P
s sm CD
_l_
CD
\
PT] = -cot CD
sm CD
a t
T]
(1.6 e)

n a a 1 n

1.2.5 Courbure et torsion

La courbure K et la torsion T par définition, sont des variations de la direction du vecteur

normal n. Pour une variation dS et dn le point M décrit une trajectoire MM' (fig. 1.4).

D
////'t

.1'.1'.1'.1',,,,

t
Figure 1.4: Changement de direction du vecteur normal

On suppose que t est le vecteur unitaire tangent à la surface dans la direction MM' et b est
un vecteur unitaire b =n x t et tel que ( t, b, n ) forme un trièdre droit. On peut alors écrire
que

Mamadou DIOP 8 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

t- dOM _dOM b =0 x t =_1 (0 x dOM) (1.7 a, b)

-/dOMI- ds ds

De la figure 1.4 on déduit que

do = K ds t - T ds b do = K t - Tb (1.8 a)
ds

Les composantes de do / ds sont

K = t _ do = dOM - do et T =_ b _ do =_ (o,dOM,do) (1.8 c, d)

ds ds2 ds ds2

( 0, dOM, do) représente le produit mixe de trois vecteurs.

[Link] Evaluation de la courbure K

dOM e do = (OM'ç dç + OM,y]dn) e (o'ç dç + 0'Y] dn) (1.9 a)

2 2
dOM e do=(o,): eOM,): )dç +(0,): eOM, +0, eOM,): )dçdY]+(o, eOM, )dY] (1.9 b)
':> ':> ':> Tl Y] ':> Y] Y]

Se rappelant que, 0 e dOM = 0 e (A t dç + B ty] dn) = 0 alors,

ç

d( 0 edOM)=(o edOM)'ç dç +( 0 edOM),y] dY]=O

(llOa,b)
d(o edOM) =(o'ç edOM+oedOM'ç )dç +( n,y] edOM+ 0 edOM,y] )dY] =0

o'ç - dOM = - 0 - dOM'ç

(1 Il a, b)
0, -dOM =-0 -dOM,
y] y]

Le développement donne

Mamadou DIOP 9 Purrlée acadénique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

n, • (OM,): dç + OM, dn) == -n. (OM,): dç + OM, d11)

11 ~ 11 ~11 1111

n, .OM, =-[Link],
11 11 1111

n'Ç .A\ =-[Link]'çç =L

n'Ç • B t = -n. OM'Ç11 = M

11

n, • A t): = -n • OM,): =M
11 ~ ~ll

n, .Bt =-[Link], =N
11 11 1111

Par substitution de (1.11) dans (1.9 b) on obtient la seconde forme fondamentale

2 2
II =[Link]=Ldç + 2Mdçd11 + Nd11 (1.12 a)

avec

L == -( n e OM'çç ) = -aX'ÇÇ - bY'çç - cZ'çç

M = -( n • OM'Ç11 ) = -aX'Ç11 - bY'Ç11 - cZ,Ç11 (1.12 b)

N = -( n • 0 M ) = -aX - bY - cZ,
'1111 '1111 '1111 1111

(1.12 c)

La courbure est

Mamadou DIOP 10 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

(1.13)

On a aussi

]=[~M
~J
D,,; • \ D,,; • t Tl
t
r Tl
D, • t
,;
D,
Tl
•
Tl A

Les courbures sont maximums ou minimums sur les courbes e = constantes et ,; =

constante. Le rayon de courbure Rs relatif à la courbe méridienne (de = 0) et celui relatif à

la courbe circonférentielle Re (ds = 0 ), lorsque l'épaisseur est constante sont représentés à
la figure 1.5 avec t l1 = i a et t E, = t

Courbe méridienne

P (r, Z)

o
ir r
Figure 1.5 : Rayon de courbure R s et Re

[Link] Evaluation de la torsion T

L'expansion du produit mixte donne

Mamadou DIOP Il Année académique 2002/2003

 Projet de Fin d'Etudes Ecole Supérieure Polytechnique

(n, dOM, dn) = ( n x dOM). dn = [n x (OM'ç dç + OM'l1 d11]. do (1.14 a)

Avec les considérations géométriques de la figure 1.3

(n,dOM, dn)= (A Plldç- B P dn (1.14 b)

ç).

Observant (8.6 e)

(n, dOM, dn) = [A (-._I_ t - cot œ t): ) dl; - B (-._I_t ): - cot œ t ) d11]. dn
sm œ 11 ':> sm œ ':> 11

(1.14c)

La torsion est alors

2 2)]
T = _1_ [-L( fdl;d11+COS (j) dç 2)+M(t dl; -Xd112)+N("*dl;d11+COS (j) d11
sin œ A 2 dç2 + 2 AB dl;d11 + B 2 dr?

(1.15)

1.2.6. Courbure et torsion dans les directions t et t

ç 11

Dans la direction \ ( d11 = 0 )

(1.16)
T =
ç
N+cot w Kç
Dans la direction t (dç = 0)
11

(1.17)

Marnadou DIOP i2 Année académique 2002/2003

---------------_._---~-------~._-~.~--~_._ , .. ..
, ~ ~ ..
Projet de Fin d'Etudes Ecole Supérieure Polytechnique

1.2.7 Spécialisation dans le cas d'une coque de révolution

Dans le cas d'une coque de révolution à épaisseur constante les vecteurs \ et

t correspondent respectivement à i e et t
11

Un repère cartésien est défini tel que l'axe Z coïncide avec l'axe de symétrie de révolution
du solide ou de la coque. Le vecteur position d'un point quelconque est défini par ses
coordonnées cartésiens X, Y, Z

ie

i
y

x
Figure 1.6 : Repère cylindrique et cartésien

Les composantes du vecteur x dans les bases cartésienne et cylindrique sont reliées par la

matrice orthogonale [Q x r] :

x cosê -sinê 0 r

y sinû cosê 0 o (1.18)

z o o 1 z

Mamadou DIOP 13 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

7.

Figure 1.7 Tronc de cône

z ..

a
A
• p t

Figure 1.8 couronne circulaire R, = Re = 00

B
t

D -------------------------~------
p
fi

a A r
Ir

Figure1.9 cylindre Cylindre Rs = 00 Re= RI

Mamadou DIOP 14 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Lorsque l'épaisseur est variable , on introduit le paramétre S (Figure 1.10 ) suivant

l'épaisseur tel que

Xq( ç ,8, Ç)=xp( ç, 8) + z( Ç) il ( ç ,8 )

z( ç )= ~ ç -1 < ~< 1
set e const.

z' ç= const

z
) Courbe méridienne
Moyenne ç = 0
fi

p
k

Figure 1.10 cas d'une coque à

épaisseur variable

1.3 Relations cinématiques

1.3.1 Hypothèse de Mindlin (Hypothèse des sections droites)

L'hypothèse des sections droites permet d'exprimer les déplacements virtuels d'un point
quelconque q en fonction des déplacements virtuels du point p situé sur la surface moyenne
( ou de référence) et en fonction d'un accroissement de déplacements virtuels dû à la
rotation de la section. Cette hypothèse nous a permis d'établir les déformations.
dans le repère curviligne.

Mamadou DIOP 15 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

a
e
es
ee cos~
as
__
1
Rs
_ sin~
0

0
l 1

li
r r 1

X = X, = 0 0 a w
as 1

1
(1.3.1 )

y Xe 0 0 CO{$ J ~

y 1 ~
Rs as

T
u = [u w ~] (1.3.2)

x y] (1.3.3 )

~ cos~ 1
0 0
as r Rs
1 _ sin~ ~
L= 0 0 (1.3.4)
Rs r as
0 0 a cos~
1
as r

1.3.2 Hypothèse de Kirchhoff ( Hypothèse de la conservation des normales )

La théorie des coques dites de Kirchhoff est basée sur l'hypothèse de la conservation des
normales : «les points matériels situés sur une normale n à la surface moyenne A avant
déformation restent sur une normale à la surface moyenne déformée .
Les déformations de membranes ainsi que les courbures ont été données dans le repère
cylindrique. On admet que les déformations de CT sont négligeables. Ce qui implique que
'FO

Mamadou DIOP . 16 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

es a __1 0
as Rs
li
cos~ _ sin~
ee 0
r r
{:} -

XS
_1Q_ aRs_1
Rsas as 2 R~
0 a
as
w (1.3.5)

coscp _ coscp Q
xe --- 0
r Rs r as

T
U = [u w ~] ( 1.3.6)

(1.3.7)

»
as
cos~
r
1 a aRs 1 _coscp
Rsas as 2Rl r Rs l
_1 sinè _coscp
L= Rs
0
Rs
(13.8)
r

0 0 eas cos~
r

1.4 Conditions d'équilibre

Les équations d'équilibre ont été établies à partir de l'expression des travaux virtuels
(P.T.V. )
W =Wint - W ext = 0 (1.4.1

a(r N ) rT
s
U a2~ a2
·--a-s----"s'---+ COS<p Ne + -R--r (f - Pm -a-t-2 - Pmf-a-t-2 )
s
=0
s
rT a(r T ) aw
2
(1.42.)
- R s - sin <p Ne as s - r (f - Pmf
s
at 2) =0
a2U a2~
s
a(r M )
a s
S e s
+ coscp N +rT - r (m - P
s mf
-a ~ -pf -a
t-' t
2) = 0

Cette relation traduite sous forme matricielle donne

Mamadou DIOP 17 Aimée académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

- 1-
--- a ar coscP 0 0 _1 N
OS ras r R e

1 siné _ ~ _lèr M
--- 0 0
R r as ros
"" s

0 0 a -
coscP
- 1 M
as r e

T
s

(3
r(f-p - -Z p - )
aZu aZ
s m t mf tZ a a
r ( f -p
s mf
a2'i)
at (lA.3)

(
r m -p
azu az~
--p--
s mf at Zat Z

avec

M TJ (104.4)
e s

aZu Pf--
r (ms- pmf - -
ê2~]
at è t Z
Z

(1.4.5)
cos~ 1
r
o o Rs
_sin~
LT =
Rs
1
r
o o (1.4.6)

o o a cos~
1
as r

Mamadou DIOP 18 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polyteclmique

Habituellement on aL' = LT mais dans le cas présent tel n'est pas le cas dans:
l 'hypothèse de Mindlin à cause du terme - aa _1 aar
s r s
l 'hypothèse de Kirchhoff à cause de négligence de la déformation de CT

1.5 Loi constitutive

Nous allons considérer une coque constituée de matériaux élastiques subissant de petites
déformations et de petits déplacements. Les relations contraintes - déformations sont
établies avec l'hypothèse de contraintes planes et d'anisotropie plane. On suppose

également que les propriétés élastiques sont axisymétriques. D'où les relations suivantes:

(J= [H] {E}+ l(J

o
J (l.5.1)

avec
o
H = Hl2 Hn o <cr >= <(Js (JS (Jsz' > <E >= <Es ES Esz' >

o o Gsz·
(1.5.2 a, b, c)

z' est une direction d'orthotropie

En considérant une coque constituée d'un empilement de couches orthotropes, les axes
d' orthotropie étant s et e .En considérant les définitions des efforts résultants et en utilisant
les relations contraintes-déformations ( équation ci-dessous) et les relations cinématiques
basées sur l'hypothèse de Mindlin ( équation 1.3.1) on obtient:
t t
N s= J 0,2 (Js dz Ns = f al (JS dz ( 1.5.3 a, b )
-t -t

t t
Ms = J 0,2 (Js z dz Ms = J al (JS Z dz (1.5.4 a, b)
-t -t
ona

Mamadou DIOP 19 Armée académique 200212003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

{N} = lH j{e} + lH j{X} +

m mf
r0
J

{M}=lH j{e}+lH j{X}+lM J ( 1.5.5)

mf f 0

T =H y+T
5 c 0

Ces relations traduites sous forme matricielle donne :

N H H 0 e N
m mf 0

M H H 0 X+ M ( 1.5.6)
mf f 0

TJ 0 0 H yJ T
l
1
5 L C 0

t t t
où Hm = f [HA] dz , Hmf
-t
= f [HA] z dz,
-t
Hf= f [HA] Z2 dz
-t
(1.5.7)

sont les rigidités homogénéisées avec

[ HA ] = l~:HII H12 ]
t a
2G
2 HI2 a HII
et H =kf
c -t a sz
dz (1.5.8a,b)
al 1

1.5.1 cas anisotrope

Dans le cas anisotrope la matrice H est donnée par

HII H12 0

[H] = HI2 H22 0 (1.5.9)

o 0 Gsz'

1.5.2 Cas orthotrope

Mamadou DIOl' 20 Année académique 1002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Dans ce cas précis la matrice H est donné ci-dessous

El V2lEI
0
1-v12Y21 1- VI2V2l

Y21El E2 (1.5.10)
[H ] 0
1- Y12Y2l 1 - YI2Y21

0 0 Gsz

< cr 0 ) = <~ T [ Hl ( afh + a1h Y2l ) Hl ( afh Y12 + a1h ) 0 ]

(1.5.11)
est un vecteur de contraintes initiales d'origines thermiques où ~T représente la variation
de la température et al th et a2 th sont des coefficients de dilatation thermique dans les
directions 1et 2

1.5.3 Cas isotrope

Dans le cas d'un matériau isotrope la matrice H est la suivante

E Ey
0
1- y2 1- y 2
Ey E (1.5.12)
[H] 1- y2 1- y 2
0

0 0 Gsz

Eath~T
< (JOt> =
1-u
< 1 1 a> (1.5.13)

est un vecteur de contraintes initiales d'origines thermiques où ~T représente la variation

th th
de la température et al et a2 sont des coefficients de dilatation thermique dans les

directions 1et 2 .

Mamadcu DIOP 21 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

1.6 Conditions au bords

Les différents types de conditions aux limites sont représentés dans le tableau suivant

Type de conditions Variables connues Variables inconnues Représentations

Encastrement Su

Appui fixe
Up=0 fi =0

U p M=O
N s r, Ms

Ns Ts fi
nt s

1
s
( rotule) : Su et Sr

1 1 \ Il

Appui roulant suivant U=O [z=O M=O Nsw Ms

n Su et Sr t

/ / .
;
4
,
f3 =0
Idem sans rotule : Su U =0 [z=O
et Sr
Nsw Ms

n
Ve>
.

~
,i
Appui roulant suivant F, w =0 M =0 U r, Ms
t: Su et Sr :
n
t
,
-,
Idem sans rotule r, w=O fi =0 UTs Ms

'# t ;

Appui roulant U=O Fz=M=O Fr W fi

" 1
1

Fr=W=M=O U r, i
fi ,1
1

Liaisons internes Ms=O Ns=O 1

- rotules
- joint au
dilatation
Symétrie W=u=O Nsw Ms
Ts=O fi =0

Tableau 1

Mamadou DIOP 22 AImée académique 2002/2003

 Projet de Fin d'Etudes Ecole Supérieure Polytechnique

1.7 Formulation variationnelle

Cette partie est relative aux aspects théoriques des coque;s de révolution dont le
comportement est axisymétrique .les champs de déplacement virtuels ont été défini avec
l'hypothèse des sections droites ( modèle de Reissner/Mindlin ), puis avec l'hypothèse
d'un champs de déplacement indépendant de e, on établit les expressions de déformations
virtuelles en utilisant la description cylindrique ou curviligne. Nous donnons l'expression
du principe des travaux virtuels dans le cas d'une épaisseur constante pour des
sollicitations volumiques, surfaciques, statiques et dynamique

. Principe des travaux virtuels en axisymétrique

Le principe des travaux virtuels dans le cas d'un problème axisymétrique s'écrit
W =Wint - Wext = 0 (1.7.1)
Lorsque l'épaisseur est constant on a

W = 2n
int
f C(e*){N}+(X*){M}+y*T )rds
S S
(1.7.2)

avec (e*) = (e; eê) (N* ) = (N* N*) ( M* ) =( M* M* '> 1

(1.7.3)
ses e
Les efforts résultants de membrane suivant s et e de flexion (moments )autour de io et
de t et l'effort tranchant sont représentés à la Figure 1.11
z

Ne
Me
-. k

1r

le
Figure 1.11 Portion de coque

Mamadou DIOP 23 AImée académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

b ) Plan i e. n

Figure 1.12 Moment Me et effort résultant de membrane Ne

Me

a )Plan t, n
Figure 1.13 Moment Ms et effort résultant de membrane Ns et effort tranchant Ts

Le travail virtuel externe est la somme du travail virtuel externe des forces de volume
et d'inertie.

W (inertie) = 211:
ext
J (u* f
S S
+ w* f + ~* f i ) r ds
z s

* ( 8zu 8z~ A* ( 8 zu p f8z~) * è 2w

J
-211: su. pm àt Z - pmf àt z t + 1--' • pmf àt z + p àt z t + w . pm ct2

+ 211: ( r u*F + r w*F + r ~* M ) s (1.7.4)

S Z f

Mamadou DIOP 24 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

t t
Avec Pm = fP
-t
al a2 dz Pmf = fP
-t
al a2 z dz (1.7.5a, b)

P -
f
J p aa
-t 1 2
Z2 dz (1.7.6)

,
OU U1=1~ <X2=1~ (1. 7.7)
Rs Re

1.8 Fonctions d'interpolation

1.8.1 Elément tronconique sans contrainte transversale (CAXI_ K )
L'élément CAXI_ K est un élément tronconique mais il est basé sur la théorie de Kirchhoff
avec approximations linéaires de la géométrie des déplacements tangentiels de la surface
moyenne et une approximation cubique du déplacement transversal .Il a trois degrés de
liberté par nœud

[Link] Matrice de Rigidité

La matrice de rigidité a été établi à l'aide du principe des travaux virtuels (P.T.V ), des
expressions des déformations de membranes et des courbures.

-1 ~ ; ~ 1

o
-~-------I---------F----+. ~
1

Figure 1.14 Elément de référence

Marnadou DIOP 25 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

z
S,u

2
L
z,
w

Figure 1.15 Elément tronconique CAXCL

Approximation de u (~) -1::;; ~::;; 1

L = 1- ç L 1+ ç
u (~) = Nr Ul + N~ U2 NI N
1 2
(1.8.1a, b, c)
2 2
AU =1(U2 -ud a U = a U o~ = 2.- ô u (l.8.2a, b)
o~ 2 as o~ as L o~

approximation de w ( ~ ):

(1.8.3a, b)

NI =l (1- ç)2 ( 2 + ç) N 2 =- L ( 1- ~ 2) ( 1- ~ )
4 8

Mamadou DIOP 26 Année académique 2002/2003

 Projet de Fin d'Etudes Ecole Supérieure Polytechnique

"

(1.8.4)

_ 6 ):
Nl,ss- 2'"
L

_ -6 ~
N 3 ,ss-- L2

or les déformations de membrane sont donnés pal'

a 0
as
{e}=
:'f=
sJ
c
r
_8
r
1:) (1.8.5)

li
1

L L L L
u=N li +N li li N 0 N 0 w
1 1 2 2 1 2 1

=> = (1.8.6a, b)

w N li
w= NI w 1 + N 3 w 2 0 NI 0 2
3

,
t
MamaJou DIOP 27 Année académique 200212003
J
Projet de Fin d'Etudes Ecole Supérieure Polytechnique

U
1

a 0 N
L
0 N2
L
0 w
as 1 1

d'où {e}=
c S
0 N 0 N u
2
r r 1 3

w
2

aNr 0
aNr 0 1 il 0 l 0
as as L

=lc~f
(1.8.7)
c N }" _NIS CNr _N3S _NIS CNr _N3S
r r r r r r r

donc dans le repère cylindrique la matrice B m est donnée par

0 0 l 0 0
L
[Bm J=[ -{
INLC _lN,S -lNzS INLC -lN3S -lN.8
Lr 1 r r r z r r
(1.8.8)

0 - NI,ss - Nz,ss 0 -[Link] -N-+.ss

[Bf] = (1.8.9)
0 - C N' s - C N zs 0 -C N3.s - C N4s
r ' r ' r r '

La matrice de rigidité est obtenue à partir de la relation suivante:

L
Wi~t = 2n J( <e*) ([Hm ]{e}+ [Hmr] {X} ) + <X*) QHmr] {e}+ [Hf ]{X})) r ds
o
(1.8.10)

Mamadou DIOP 28 AImee académique 2Ù02/21)(J3

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

=(u*)
n lac
[k] lac ~ n Jlac (1.8.11)

[k]loc =2n fo [kç]dS=2n -1f [kç] ~ dç (1.8.12) :

(1.8.13)

[Link] Matrice de masse

La matrice de masse est obtenue, elle à partir de la relation suivante:

L 8 2u 8 2w
828 828 8 2u
Wi~t = 2n Je U*pm àt2 + W*pm àt2 + 8*pm 81: 2 + U*pmf àt2 + 8*pmf àt 2 ) r ds
o
(1.8.14)

2
8U }
=(U*) [m] {
àt 2n (1.8.15)
n lac lac
lac

2 2
ae aw
avec e*= - w*,s -=-~(-) (1.8.16)
2 as 2
at at
d'où

[m] =2n
lac 0
f [m ]dS=2n J[m ] ~ dç
ç -1 ç ( 1.8.17)

Mamadou DIOP 29 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

+ {NS} P f <NS ) - {N LJP mf <Nt, ) - tNt, J p m <NL ) ) r (1.8.18)

<NL)=<NL 0 0 NL 0 0>
1 2

<NC) = <0 N NON N) (1.8.19)

2 3 4

<NC ) = <0 N NON N)

,S 1,s 2,s 3,5 4,5

[Link] Forces équivalentes aux nœuds

Le vecteur des forces équivalentes aux nœuds est:

{fJ=2n J {f} ~ ds (1.8.20)

{fJ={Nc)f, +{Nc)f, +{N~}m, (1.8.21)

Les matrices {N L} {NC} {NS} sont définies en (1.8.19)

1.8.2 Elément tronconique linéaire (caxi_L )

2
s

22 ----------------------------
2
L

o r

Mamadou DIOP 30 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polyteclmique

Figure 1,16 Géométrie d'un élément génératrice d'un tronc de cône

L'élément CAXC L est un élément tronconique mixte à deux nœuds avec approximations
linéaires de la géométrie ( coordonnées r et Z ), des déplacements et rotations et des
approximations constantes des efforts tranchants .L' élément CAXI L a trois degrés de
liberté par nœuds.

s, u
Z,W

'---'r-------------.
o B r,U
Figure 1.17 Elément tronconique linéaire CAXCL

On pose C = cos <p

S = sine
La relation liant le repère global au repère local est la suivante:

Mamadou DiOP 31 Année académique 200212003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

( 1.8.22)

C
1

U=N U +N U U o o \\T
1
1 1 2 2

=> = (1.8.23)

W=N W +N W W o o
1 1 2 2

Vi
2

or les déformations de membrane sont donnés par

r
C~ s~
aS as fU) (1.8.24)
1
Lr
o lw
U
1

C2- S2- N 0 N2 0 W
as as 1 1

d'où {e}= (1.8.25)

1 0 0 N 0 N U
r 1 2 2

Mamadou DIOP 32 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

aN aN
r-c -C

;,
cas-I 0 cas-1 0 0
L
0

N
1
0
N
_1
0
~l 0
N
r
1
0
(1.8.26)

r r

donc dans le repère cylindrique la matrice Bm est donnée par:

C s 0 C s 0
L L L L
N N (1.8.27)

r
o 0
r
2
o 0

(u (U W ~ U W ~ ) (1.8.28 )
n) 1 1 2 2 2

les courbures sont données par:

{X} =
Xs

Xe
ras

l?~
or ~ = ~IN 1+ ~2N 2 ::::::> {X}=
C
-L~l -L~2

l-~lN1 +-~2~2
r r
C
(1.8.29)

donc dans le repère cylindrique la matrice Br est donnée par

o o o o i-
L
(1.8.30)
NIC
o o o o
r

les déformations de contraintes transversal (CT) sont données par:

y= ~ + \V = ~ - sau
as + C aw
as or ~ = ~ N + ~ N
1122
,U = N U + NU,
1122
(1.8.31)

W=N I WI +N 2 W

donc dans le repère cylindrique la matrice Be est donnée par

Mamadou DIOP 33 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

(B c )=(~ -i: (1.8.32)

[Link] Matrice de Rigidité

L
Wi~t = Tn Je (e*)( [Hm]{e} + [Hmf]{X} ) + (X*) QHmf ]{e} + [Hr]{z}))
o
(1.8..33)

(1.8.34)

(k mft)

L
avec [k mf l- 2n J e[BmY e[Hm] [B m]+[Hmf] [Br ])
o

(1.8.35)

k mft = 2 1t fo ~Jr ds (1.8.36)

[Link] Matrice de masse

la matrice de masse est obtenu à partir de l'équation suivante

a2u a2fJ a2u
W~t = 2n
L
Je
o
U; (Pm----f
Dt
+ Pmf-2 t) + fJ* t epmf----f + Pr t) ) r ds
Dt Dt
(1.8.37)

Mamadou DIOP 34 Année académique 20~2:2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Dans le cas d'un matériau homogène, elle se présente sous la forme suivante

al 0 -2bC rm 0 -be

al -2bS 0 rm -bS
h2 al -bC -bS h2 rm
U U
[m]=27tph ~ (1.8.38)
a2 0 -2bC

a2 -2bS

SYM h2 a2
12

(1.8.39)

[Link] Forces équivalentes aux nœuds

On distingue 4 types de sollicitations:

Les sollicitations surfaciques, volumiques, linéaires et celles qUI sont dues à la
température.
• Sollicitations volumiques
Il s'agit principalement des forces de gravité et des forces découlant des pressions
moyennes exercées sur la surface moyenne dont leurs expressions sont les suivantes:
- forces de gravité

f = - 2n p g h L <0 (2r + r ) - ( hS)2 0

nv 6 1 2 2
(2r + r )
2 2
-(hi )2 ) (1.8.40)

- Si la surface moyenne est soumise à une pression uniforme on a

f =- 2n p L <- Sa Ca 0 -Sa Ca 0) (1.8.41)

nv 6 1 1 2 2

• Sollicitations linéiques au nœud 2 :

f = - 2n <0 0 0 F FM) (1.842)

ns r2 z2 2

• Sollicitations dues à la température:

Mamadou DIOP 35 Année académique 2002/2003

 Projet de Fin d'Etudes Ecole Supérieure Polytechnique

E athh r L T
m 0
cr 1- V

(l.8.43)

1.8.3 Elément isoparamétrique quadratique (CAXI_Q)

Z,W

_1 1

-1

E1ément de référence

E1ément réel

Figure 1.18 E1ément quadratique CAXI_Q

[Link] Matrice de rigidité

La forme intégrale élémentaire Wint est

(1.8.44 )

1 1
k = 2 1t JJk
-1 -1 ~
dÇ dç (1.8.45)

Mamadou DIOP 36 Année académique 200212003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

avec (1.8.46)

B est la matrice qui lie le vecteur des déformation réelles au vecteur Un

et J est le Jacobien. H est la matrice définie en (1.5 .2a )

[Link] Matrice de masse

Pour définir la matrice de masse cohérente, on considère

(1.8.4ï )

a2U- a2U
avec ab" =( .... at2; .... ) (1.8.48 )

1 1

[m)=27t f f [N)rp [N)JdçdÇ

-1 -1
(1.8.49)

[Link] Forces équivalentes aux nœuds

Pour les sollicitations volumiques les charges équivalentes sont données par :

(1.8.50)

Mamadou DIOP 37 Année académique 2002/2003

Projet de Fin d'Etudes Ecole Sup érieure Polytechnique

2 Application Pratique

Figure 1.19 : Photo d'un bassin de décantation

A l'aide du logiciel Microfe on a dimensionné un réservoir de type cylindrique comme

montré ci-dessus. Ainsi les efforts internes, le mémoire disque et le temps de calcul sont
évalués. Le nombre d'éléments, de nœuds et d'équation utilisé est considérable

MamadouDJOP 38 Année académique 200212003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

ANALYSE STATIQUE DU RESERVOIR

(proFE-80999/20699)

--Contrôler le traitement des sous-structures: 0 sec

--Valeurs système globales

Éléments Noeuds Équations Contraintes Espace disque

1580 1525 9150 4477881 34 Mbyte

Cas de charge :3

--Valeurs système Sous-/Superstructure

Sous Structure/Superstructure Principale- Interne
Elément Espace mémoire Nœud Equation Nœud Equation

1 42876 334 Kbyte 24 144 48 288

2 1350942 10 Mbyte 108 648 396 2376
3 18666 145 Kbyte 24 144 12 72
4 65901 514 Kbyte 36 216 41 246
5 3178812 24 Mbyte 120 720 872 5232
DECANT 334836 2615 Kbyte 0 0 156 936

--Espace disque nécessaire

Mémoire nécessaire disponible

Calculatrice d'équation 26 Mbyte OUI

Espace disque nécessaire disponible Répertoire:\

Résultats 4333 Kbyte 8510 Mbyte 'C:\MB\PROJEKTE\P...'

Marnadou DIOP 39 AImée académique 2002/2003

Projet de Fin d'Etudes Ecole Supérieure Polytechnique

Bloc 42 Mbyte 8510 Mbyte 'C:\MB\MicroFe\BL...'

--Résolution de la fonction statique

Sous-struc./ Temps calcul
Superstructure sec

1 1
2 25
3 1
4 1
5 39
DECANT 8

--Charge totale / Réaction d'appui totale

Cas Px i Ax Py ! Ay Pz / Az
[Kn] [Kn] [Kn]

1 -0.00 / 0.00 -7.07/7.07 -14420.82/ 14420.82

2 -24.88 / 24.88 -45.20 / 45.20 -25116.46/25116.46
3 0.00 / -0.00 0.00; -0.00 -1674.43/1674.43

--Résultats
pour la Sous structure: 3 sec

--Fin de l'analyse statique

Temps de calcul: 103 sec

*** Calcul correct ***

Mamadou DIOP 40 Année académique 2002/2003

 Projet de Fin d'Etudes Ecole Supérieure Polytechnique

CONCLUSION ET RECOMMANDATIONS

Afin de modéliser la coque de révolution plusieurs éléments sont utilisés: l'élément

CAXI_L, CAXI_K pour les coques à épaisseur constante et l'élément CAXI_Q pour les
coques à épaisseur variable. Pour chaque élément on a déterminé la matrice de rigidité, de
masse et les forces équivalentes aux nœuds .Ces trois facteurs une fois connu permettent
d'évaluer les déplacements aux niveaux des nœuds, les efforts internes qui comparés aux
valeurs de référence nous donne une idée sur la fiabilité de chaque l'élément. Des patchs-
tests auraient pu être faits pour évaluer les performances et la précision des modèles. Après
de telles vérifications, une comparaison entre le logiciel Microfe et les résultats découlant
des éléments testés nous édifiera plus amplement sur la validité du modèle.
Enfin nous pensons que le travail qui a été fait sera poursuivi dans le cadre des projets de
fin d'étude pour déboucher sur une plate forme qui exploitée donnerait de bons résultats.

Mamadou DIOP 41 Année académique 2OV2/20Û3

... _..,-
~. - ..
 BIBLIOGRAPHIE

[1] Jean -Louis BA TOZ and Gouri DATT. Modélisation des structures par éléments finis.
Vol 3. Mars 1992.564 p.

[2] Kenneth H. HUEBNER , Earl A. THORHTON and Ted G. BYROM. The finite
element method .Third edition . 1995.627 P

[3] O.C ZIENKIEWICZ, FRS and R. L TAYLOR .The finite element method. Fourth
edit ion . Vol 1. 1988.648 P
[4] Thomas J.R. HUGHES. The finite element method.1987.803 p

[5] Jean -Louis BATOZ and Gouri DATT. Modélisation des structures par éléments finis.
Vol2. septembrel995. 483 p.

[6] Jean -Louis BATOZ and Gouri DATT. Modélisation des structures par éléments finis.
Vol 1. Mars 1992.564 p.
 ANNEXES

Résultais du dhnenslonnement du réservoir à l'aide du logiciel Microfe

Plan de repérage

Elé D-2 Plancher

Structure x = -2.04 -3.54 -4.09 -3.54 -2.04 0.00 2.04 m

y= -3.54 -2.04 0.00 2.04 3.54 4.09 3.54m
z= 1.74 1.74 1.74 1.74 1.74 1.74 1.74 m
x= 3.54 4.09 3.54 2.04 0.00 -2.04 m
y= 2.04 0.00 -2.04 -3.54 -4.09 -3.54 m
z= 1.74 1.74 1.74 1.74 1.74 1.74 m

Matériau Coque isotrope

Epais 15.0cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kJ~/m2
Nu 0.20
Surface = 50.18 m"2
 Volume = 7.53 m A3

Etage [Link]

Elé D-l Plancher

Structure x = -14.53 -8.39 0.00 8.39 14.53 16.77 14.53 m

y = 8.39 14.53 16.77 14.53 8.39 0.00 -8.39 m
z= 0.00 0.00 0.00 0.00 0.00 0.00 0.00 m
x= 8.39 0.00 -8.39 -14.53 -16.77 -14.53 m
y= -14.53 -16.77 -14.53 -8.39 0.00 8.39 m
z= 0.00 0.00 0.00 0.00 0.00 0.00 m

Matériau Coque isotrope

Epais 45.0 cm
Poids = 25.00 kN/m3
Mad E = 3.00e+007 kN/m2
Nu 0.20
Surface = 837.22 m A2
Volume = 376.75 m A3

Etage Fondation

Elé W-25 Mur

Structure x = -3.54 -4.09 m

y = -2.04 0.00 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids = 25.00 kN/m3
Mad E = 3.00e+007 kN/m2
Nu = 0.20
 Surface = 3.79 m"2
Volume = 0.57 m /\3

Etage Déversoir

Elé W-12 Mur

Structure x = -16.48 -14.27 m

y = 0.00 8.24 m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 ml\2

Volume = 10.52 m"'3

Etage Murs

Elé W-ll Mur

Structure x = -14.27 -16.48 m

y = -8.24 0.00 m
Niveau inf. = 0.00 rn, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 ml\2
Volume = 10.52 ml\3
Etage Murs

Elé W-10 Mur

Structure x = -8.24 -14.27 m

y= -14.27 -8.24m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 m"2
Volume = 10.52 ml\3

Etage Murs

Elé W-9 Mur

Structure x= 0.00 -8.24 m

y= -16.48-14.27m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 m"2
Volume = 10.52 ml\3

Etage Murs
Elé W-S Mur

Structure x= 8.24 0.00 m

y= -14.27-16.48m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
MNu 0.20
Surface = 30.06 m A2
Volume = 10.52 m A3

Etage Murs

Elé W-7 Mur

Structure x = 14.27 8.24 m

y = -8.24 -14.27 m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu = 0.20
Surface = 30.06 m A2
Volume = 10.52 m A3

Etage Murs

EléW-6 Mur
Structure x= 16.48 14.27 m
y= 0.00 -8.24 m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 ml\2
Volume = 10.52 ml\3

Etage Murs

Elé W-S Mur

Structure x = 14.27 16.48 m

y= 8.24 0.00 m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 ml\2
Volume = 10.52 ml\3

Etage Murs

Elé W-4 Mur

Structure x= 8.24 14.27 m

 y = 14.27 8.24 m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 35.0 cm
Poids 25.00 kN/m3
Mad E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 m A2
Volume = 10.52 m A3

Etage Murs

Elé W-2 Mur

Structure x = -4.09 -3.54 m

y= 0.00 2.04 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids = 25.00 kN/m3
Mad E = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 m A2
Volume = 0.57 m A3

Etage Déversoir

Elé W-l Mur

Structure x= 0.00 8.24 m

y= 16.48 14.27m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m
Matériau Coque isotrope
Epais 35.0 cm
Poids 25.00 kN/m3
Mad E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 m A2
Volume = 10.52 m A3

Etage Murs

Elé W-13 Mur

Structure x= -0.75 -1.30m

y =-= 1.30 0.75 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 f i

Matériau Coque isotrope

Epais = 20.0 cm
Poids = 25.00 kN/m3
Mad E = 3.00e+Oü7 kN/m2
Nu 0.20
Surface = 1.35 m A2
Volume = 0.27 m A3

Etage cylindre_inf

Elé W-14 Mur

Structure x = -1.30 -1.50 m

y= 0.75 0.00 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

 Epais 20.0 cm
Poids 25.00 kN/m3
Mad E = 3.00e+007 kN/m2
Nu 0.20
Surface =
Volume =

Etage cylindre_inf

Elé W-15 Mur

Structure x= -1.50 -1.30m

y= 0.00 -0.75 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais 20.0 cm
Poids 25.00 k..~/m3

Mad E = 3.00e+007 kN/m2

Nu 0.20
Surface =
Volume =

Etage cylindre_inf

Elé W-16 Mur

Structure x= -1.30 -0.75m

y= -0.75 -1.30m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais 20.0 cm
Poids 25.00 kN/m3
 Mad E = 3.00e+007 kN/m2
Nu = 0.20
Surface =
Volume =

Etage cylindre_inf

Elé W-17 Mur

Structure x = -0.75 0.00 m

y = -1.30 -1.50 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais 20.0 cm
Poids 25.00 lu~/m3

Mad E = 3.00e+007 kN/m2

Nu = 0.20
Surface =
Volume =

Etage cylindre_inf

Elé W-18 Mur

Structure x= 0.00 0.75 m

y = -1.50 -1.30 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais = 20.0 cm
Poids = 25.00 kN/m3
Mad E = 3.00e+007 kN/m2
Nu 0.20
1
Surface = 1.35 m A2
Volume = 0.27 m A3

Etage cylindre: inf

Elé W-19 Mur

Structure x= 0.75 1.30 m

y= -1.30 -0.75 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais 20.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 1.35 m A2
Volume = 0.27 m A3

Etage cylindre_inf

Elé W-20 Mur

Structure x= 1.30 1.50 m

y = -0.75 0.00 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais = 20.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu = 0.20
Surface = 1.35 m A2
Volume = 0.27 m A3
Etage cylindre_inf

'Elé W-21 Mur

Structure x= 1.50 1.30 m

y= 0.00 0.75m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais = 20.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 1.35 ml\2
Volume = 0.27 ml\3

Etage cylindrejnf

Eié \"-22 Mur

Structure x= 1.30 0.75 m

y= 0.75 1.30 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais = 20.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 1.35 ml\2
Volume = 0.27 ml\3

Etage cylindre_inf
 Elé W-23 Mur

Structure x= 0.75 0.00 m

y= 1.30 1.50 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais 20.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 1.35 ml\2
Volume = 0.27 ml\3

Etage cylindrejnf

Elé W-24 Mur

Structure x= 0.00 -0.75 m

y= 1.50 1.30 m
Niveau inf. = 0.00 m, Niveau sup. = 1.74 m

Matériau Coque isotrope

Epais 20.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 1.35 ml\2
Volume = 0.27 ml\3

Etage cylindre_inf

Eléw-l Mur

_ _ _ _,_~" _ _.... ~ ~ , _ ~ - . - -••• •_ _ J.. •.••••

Structure x = -3.54 -2.04 m
y= 2.04 3.54 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 m A2
Volume = 0.57 m A3

Etage Déversoir

Elé w-2 Mur

Structure x = -2.04 0.00 m

y= 3.54 4.09 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu Surface = 3.79 m A2
Volume = 0.57 m A3

Etage Déversoir

Elé w-3 Mur

Structure x= 0.00 2.04 m

y= 4.09 3.54 m
 Niveau inf. 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 m A2
Volume = 0.57 m A3

Etage Déversoir

Elé w-4 Mur

Structure x= 2.04 3.54 m

y= 3.54 2.04 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 m A2
Volume = 0.57 m A3

Etage Déversoir

Elé w-S Mur

Structure x= 3.54 4.09 ID

Y= 2.04 0.00 ID

Niveau inf. = 1.73 ID, Niveau sup. = 3.52 ID

Matériau Coque isotrope
Epais 15.0 cm
Poids 25.00 kN/m3
ModE = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 m"2
Volume = 0.57 m"3

Ètage Déversoir

Elé w-6 Mur

Structure x= 4.09 3.54 m

y= 0.00 -2.04 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 rn" 2
Volume = 0.57 ml\3

Etage Déversoir

Elé w-7 Mur

Structure x= 3.54 2.04 m

y = -2.04 -3.54 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais = 15.0 cm
 Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 m A2
Volume = 0.57 m A3

Etage Déversoir

Elé w-S Mur

Structure x= 2.04 0.00 m

y = -3.54 -4.09 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

lviatériau Coque isotrope

Epais 15.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface =

Volume =

Etage Déversoir

Elé w-9 Mur

Structure x= 0.00 -2.04 m

y = -4.09 -3.54 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids = 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
 Nu 0.20
Surface =

Volume =

Etage Déversoir

Elé w-IO Mur

Structure x = -2.04 -3.54 m

y= -3.54 -2.04 m
Niveau inf. = 1.73 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 15.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 3.79 m A2
Volume = 0.57 m A3

Etage Déversoir

Elé w-ll Mur

Structure x = -14.27 -8.24 m

y= 8.24 14.27 m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais 45.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 m A2
 Volume = 13.53 m A3

Etage Murs

Elé w-12 Mur

Structure x = -8.24 0.00 m

y= 14.27 16.48 m
Niveau inf. = 0.00 m, Niveau sup. = 3.52 m

Matériau Coque isotrope

Epais = 45.0 cm
Poids 25.00 kN/m3
Mod E = 3.00e+007 kN/m2
Nu 0.20
Surface = 30.06 m A2
Volume = 13.53 m A3

Etage Murs

Elé F-l App surfacique

Structure x = 14.53 16.77 14.53 8.39 0.00 -8.39 -14.53 m

y = -8.39 0.00 8.39 14.53 16.77 14.53 8.39 m
x = -16.77 -14.53 -8.39 0.00 8.39 14.53 m
y= 0.00 -8.39 -14.53 -16.77 -14.53 -8.39 m
Niveau = 0.0 m

Appui Appui Compr/Trac Trans suivant axe x = 2.00e+005 kN/m3

Appui Compr/Trac Trans suivant axe y = 2.00e+005 kN/m3
Appui Compr Trans suivant axe z = 2.00e+005 kN/m3

Ètage Fondation
Moments ms, toutes les surfaces
Cas 1 Min = -16.64kNm/m
Max = 25.82kNm/m
Pas = 2.00kNm/m

Moments Ms
 il! :O:::JO.S::
200.C::
L ::'OO.c::

-~00.::::

l" -:0 :::J 0 . :: ::

b- -,;.:::JO - :::

ï-=:::Jo.::c
--::=oo.c:':

Cont Normale sr ,toutes les surfaces

Cas 1 Min = -590.1 OkN/ml\2
Max = 293.71kN/m1\2
Pas = 100.00kN/m1\2

Contrainte nonnale sr