REPUBLIQUE DU CAMEROUN

Paix - Travail - Patrie

UNIVERSITE DE DOUALA
ECOLE DOCTORALE DES SCIENCES FONDAMENTALES ET APPLIQUEES
Unité de Formation Doctorale des Sciences de l’Ingénieur

Laboratoire Materiaux et Mecanique

MEMOIRE
Présenté et soutenu en vue de l’obtention du diplôme de
Master 2 Recherche en Sciences de l’Ingénieur
PARCOURS TYPE : GENIE MECANIQUE ET MATERIAUX
Par :
Ing. TAM LIKENG Christian

PROBLEMES DE STABILITE DES STRUCTURES

COMPRIMEE A COQUE MINCE

Cas du FLAMBEMENT
 PLAN DE LA PRESENTATION

⸕ Contexte et problématique

⸕ Revue de littérature

⸕ Matériels et méthodes

⸕ Résultats et discussions
 Mots Clés

⸕ Flambement

⸕ Valeur critique

⸕ Théorie des coques minces

⸕ Conditions aux limites

⸕ Charges critiques

⸕ Méthode expérimentale
 Contexte et problématique
⸕ Contexte

Génie Civil

Une coque est un élément d'épaisseur faible par rapport aux

dimensions globales et possédant un rayon de courbure. C’est
un élément surfacique tridimensionnel prenant en compte
simultanément les effets de membrane et de flexion.

Construction industrielle

Aéronautique et naval
 Contexte et problématique
⸕ Problématique
Ces structures sont soumises à des sollicitations diverses, telle que:
• La compression
• La traction P1: Comment cerner le
• La flexion phénomène d’instabilité de
2nd ordre?
Effet engendré

• Des déplacements
• Des contraintes P2: Comment les résorber
afin d’optimiser les sections
utilisées?
2Nd Ordre

• Flambement
• Deversement / voilement
 Contexte et problématique
⸕ Problématique (suite)

la géométrie différentielle

La théorie des coques fait intervenir

La mécanique
plusieurs domaines de la science

Les EDP
 Contexte et problématique
⸕ Problématique (suite)

La théorie des coques est basée sur

deux principales hypothèses

• L’hypothèse de Kirchhoff Love pour les coques minces

est moins bien vérifiée pour les coques épaisses

• L’hypothèse Hencky – Mindlin – Reissner – Nagdhi

pour les coques épaisses bien que prenant en compte
le cisaillement transverse, une discrétisation directe de
ces modèles par des éléments finis conformes conduit
à des instabilités numériques dès que l’épaisseur
devient trop petite.
 Objectifs de l’étude
⸕ Objectif général

Etudier et analyser les modes de ruines des ouvrages de type

coques, et proposer des solutions pour limiter ces ruines
structurelles

⸕ Objectif Spécifique

Déterminer pour une coque soumise à une sollicitation de type compression axiale
l’expression de l’instabilité de 2nd ordre tel que le flambement,
Justifier en fonction des matériaux l’impact de cette instabilité et proposer des
solutions pour résorber ce phénomène,
 Revue de littérature

Les équations locales

les hypothèses cinématiques de KIRCHOFF-Love et Reissner-Mindlin[1]

la loi de comportement

Loi de Hooke, Effort membranaires

Effet de second ordre du fait

les déformations ne sont pas proportionnelles aux

contraintes

nos éléments de coques soumis au

flambement comme des éléments non
lineaires
 Revue de littérature

Formulation variationnelle

Charges critiques

Influence des conditions aux limites

Choix du matériaux et analyse

 Revue de littérature
Définitions
Coques minces

Rapport épaisseur/dimensions : L'épaisseur ( h ) est très faible par rapport au rayon de courbure ( R )
ou aux dimensions principales (généralement h/R<1/20h/R < 1/20h/R < 1/20)

Hypothèses : Basée sur la théorie des coques minces (par exemple, théorie de Love-Kirchhoff). Les
contraintes transversales (cisaillement transverse, compression normale à la surface) sont négligées. Les
déformations sont principalement dues à la flexion et à l'étirement de la surface moyenne

Comportement : La coque résiste principalement par des efforts membranaires

(tractions/compressions dans le plan de la coque) et, dans une moindre mesure, par des moments de
flexion)

Exemple: Réservoirs à paroi fine, coques d'avions, dômes architecturaux légers

 Revue de littérature
Définitions
Coques épaisses

Rapport épaisseur/dimensions : L'épaisseur est significative par rapport au rayon de courbure ou aux
dimensions (généralement h/R>1/20h/R > 1/20h/R > 1/20)

Hypothèses : Utilise la théorie des coques épaisses (par exemple, théorie de Reissner-Mindlin),
qui prend en compte les cisaillements transversaux et les contraintes normales à travers
l'épaisseur

Comportement : Les efforts membranaires, les moments de flexion et les cisaillements

transversaux jouent un rôle important. Les déformations à travers l'épaisseur ne sont pas
négligeables

Exemple: Coques de réacteurs nucléaires, tunnels en béton, coques de navires à paroi épaisse
 Revue de littérature
La formulation variationnelle pour
l’études ces coques minces

formulation variationnelle en déplacement P. Mandal [2]

2000
lorsque des structures de coque à paroi mince sont chargées en
compression, leur résistance est limitée par le flambement et comme ce
flambement est souvent catastrophique, une prédiction fiable de la
résistance au flambement des coques est une aspiration forte des
ingénieurs Il faudra donc dans le cadre de notre étude, spécifier le type de
matériaux ainsi que la nature du matériau utilisé.

Etudes des instabilités et charge critique

Chou [3] 1997, première incidence d'importance pratique d'importance pratique est associée à la
conception du pont Britannia du pont Britannia en 1845 par Robert Stephenson

Bryan [4] 1998 , a publié un article fournissant une solution pour le flambement initial d'une plaque
carrée d'une plaque carrée simplement supportée, qui a donné lieu à un grand nombre d'articles sur le calcul
des charges critiques des ponts .

Timoshenko [6] 1938, calcul des charges critiques des plaques présentant une grande variété de
conditions aux limites.
 Revue de littérature
Cependant, dans structures aéronautiques, les tôles minces sont souvent
utilisées au-delà de leurs charges critiques, et la charge qui peut supporter la
structure est déterminée par la résistance ultime en compression

Essais réalisés par Schuman et Back [13] 2003,.

ont indiqué que pour les plaques suffisamment larges et minces pour se
déformer sous la charge, la charge ultime pouvant être supportée
n'augmentait pas proportionnellement à la largeur, et était en fait
relativement insensible à l'augmentation de la largeur. Ils ont conclu
qu’après le flambage, une plaque se comporte comme si seule une partie
de de sa largeur est efficace pour supporter la charge

Hubert Debski [14] 2018,

Des colonnes à paroi mince en matériaux composite de carbone-epoxy
avec une section transversale en forme de chapeau suivant le system
HexPly M12, Hexcel. Cette étude se fait à paritr des essais sur une tole
pliée dans une section transversale en forme de chapeau avec les
dimensions particulières, l’analyse du modele de calcul sur la base de la
méthode des éléments finis a été utilisé pour l’analyse du flambage et du
comportement après flambage des éléments à paroi mince
 Revue de littérature
Pour rappel,

Jamal [18], 1999

démontre que les instabilités de second ordre tels que les flambements
peuvent être du à des imperfections localisées, le cas des coques
cylindriques soumises à la compression axiale et dont l’étude s’est basée
sur une analyse à double échelle incluant les modes d'interaction,
Ceci conduit à un système d'équations différentielles non linéaires
complexes associés à des dérivées discontinues, ces équations ont permis
la réduction de la charge critique de flambement proportionnellement à
l’imperfection localisée

Analyse des problèmes de flambement des cylindres sous compression axiale

Legrognec [19] en 2003
parmi l’ensemble des modèles de coque existants, la théorie de Donnell
[19] est la plus fréquemment utilisée pour l’analyse de problèmes de
flambement de cylindres sous compression axiale. Elle permet aisément, le
plus souvent par le biais de méthodes énergétiques, de calculer la valeur
critique du chargement en élasticité pour un mode axisymétrique ou non
axisymétrique faisant apparaître simultanément des ondes axiales et
circonférentielles ; en abordant les aspects dans le domaine plastique il
met en évidence le fait que pour une coque assez épaisse, le flambement a
toutes les chances de se produire dans le domaine plastique, de plus il
aborde les questions liées aux conditions aux limites avec un encastrement
des bases doit fournir la même valeur critique et le même type de mode
 Revue de littérature
Influence des conditions aux limites

Batterman [20] 1999

Etudie analytiquement le cas de bords libres, se plaçant ainsi dans la
situation la plus défavorable parmi toutes les conditions aux limites
envisageables. Il propose donc une solution analytique du flambement
élastoplastique d’un cylindre sous compression axiale il a par ailleurs
travaillé sur une formulation lagrangienne totale où la déformation est
décomposée de manière additive et les équations constitutives écrites sur
le modèle des matériaux standard généralisés dont le cadre
thermodynamique a fait ses preuves en plasticité de petites
transformations
Debski et Katarzyna Falkowicz [16],

la structure de coque mince dans sa structure détaillée peut être assimilée à

une poutre, poteau nervure constitué d’éléments mince entremêlé les uns
sur les autres, ces éléments (plaque) sont soumis à la compression ou au
cisaillement et peuvent conduire à des flambements et donc à une ruine
accélérée de la structure
 Matériels et Méthodes
Matériels

1. Analyse des différentes formes fondamentales

Formes de surfaces quadratiques

Deuxième forme fondamentale

Troisième forme de surface quadratique

Courbure de Gauss

2. Théorie Bidimensionnelle des coques

 Matériels et Méthodes
Matériels

2. Théorie Bidimensionnelle des coques

Tenseurs de déformations
𝜀= 𝜀 𝑖𝑗 . (𝑔
¿ 𝑖 𝑗
𝑔 )

Tenseurs de contrainte dans la membrane

Formulation variationnelle
tridimensionnelle lié à l’équilibre se présente comme suit :

(2.63)

Où est l’espace de Sobolev et représentent respectivement les produits scalaire des tenseurs et
vecteurs.
 Matériels et Méthodes
Matériels
En laboratoire

Des viroles en aciers dont les dimensions varient en hauteur, diamètre, épaisseur

Une presse de compression

 Matériels et Méthodes
Méthodes
Geometrie Differentielle
Système des coordonnées de surface

La représentation paramétrique de la surface peut être donnée comme suit :

(1- 25)

L'équation paramétrique (1.22) peut également être écrite sous la forme

surface moyenne et système de coordonnées

(1- 26)
Principales directions et lignes de courbures
Considérons une petite région d’une surface lisse S près d’un point typique M. « Lisse » signifie que la surface est continue
Ω et ne contient aucune discontinuité de pente, c’est-à-dire des plis ou des sommets. Si nous dessinons diverses courbes le
long de la surface S à travers point M, les tangentes à ces courbes sont placées sur un plan appelé plan tangent à la surface
en M (plan P sur la). Une ligne perpendiculaire au plan tangent et passant par le point M est appelé la normale à la surface S
au point M, et elle est noté ou . Puisque la surface est lisse (comme défini ci-dessus), le plan tangent et par conséquent, la
normale est déterminée. Une section normale de la surface S au point M peut être définie
Plan tangent au vecteur noraml
 Matériels et Méthodes
Méthodes
Geometrie Differentielle
Formes de surfaces quadratiques

La première forme fondamentale

(1- 30) Les quantités et sont appelées

surface moyenne et système de coordonnées
paramètres de Lamé
(1- 35)

La deuxième forme fondamentale

(1- 43)
La formule de Frenet

Vecteur unitaire tangent

(1- 46)

Vecteur normaux unitaires Deuxieme forme quadratique de surface Coefficient de forme

 Matériels et Méthodes
Méthodes
Geometrie Differentielle
Formes de surfaces quadratiques

La troisième forme fondamentale

(1- 50) Les vecteurs unités normales dirigé
vers le centre de courbures O, avec un
rayon de courbure noté R

En admettant l’égalité pour 𝜕 𝒓 ,1 𝜕𝒓,2 𝜕 𝒓

2
𝜕 𝒓
2
déterminées les dérivées des vecteurs = =
unités tangentes 𝒅 𝒙2 𝒅 𝒙1 1
𝜕𝒙 𝜕𝒙
2 2
𝜕𝒙 𝜕𝒙
1

𝑪 𝜕
2 (
Variation du vecteur unitaire suivant une coordonnée La prise en compte de la première 𝑨1 𝒂 1 )= 1 (
𝑨 2 𝒂2 )
donnée 𝜕𝒙
forme fondamentale nous donne: 𝜕𝒙

𝒒 =𝒂2 ⋅
L’introduction du symbole de surface
de Christoffel et la courbure 1
𝜕 𝒂1
𝜕𝒙
1
=
𝜕
𝜕𝒙
1(
𝒂 1 ⋅ 𝒂2) −
𝜕 𝒂2
𝜕𝒙
1
⋅
( )
𝒂1 =− 𝒂 1 ⋅
𝜕 𝒂2
𝜕𝒙
1
=− 𝒂 1 ⋅ 𝒂2 , 1

L’expression du tenseur de courbure 𝒕 1 =− 𝒂1 ⋅ − ( 𝑨1

𝑹1 1
𝒂 =
𝑨1
𝑹1 )
 Matériels et Méthodes
Méthodes
Geometrie Differentielle
Formes de surfaces quadratiques

La troisième forme fondamentale

𝜕 𝒂1
( 𝜕 𝑨1
)
𝑨1
Par substituons des équations présentées
1
plus haut on obtient l’expression de la =¿ − 𝒂 2 + 𝒂3
dérivée des vecteurs unitaires , ,
𝜕 𝒙1 𝑨2 𝜕𝒙 2
𝑹1

[ ][ ]
𝝏 𝒂𝟏
𝝏 𝒙𝟐
𝟎
(
𝟏 𝝏 𝑨𝟐
𝑨𝟏 𝝏 𝒙 𝟏 ) 𝟎

[]
𝒂𝟏
( )
𝝏 𝒂𝟐 𝝏 𝑨𝟐 𝑨𝟐
𝟐
= − 𝟏 𝟎 − 𝒂𝟐
On obtient les deux matrices et 𝝏𝒙 𝑨𝟏 𝝏 𝒙 𝟏 𝑹𝟐
𝒂𝟑
𝝏 𝒂𝟑 𝑨𝟐
𝝏 𝒙𝟐 𝟎 𝟎
𝑹𝟐

𝝀 𝝀
Avec leur forme contractuelle mathématique 𝒂3 , 𝜷 =− 𝒃𝜶𝝀 𝒂 =− 𝒃𝜶 𝒂 𝝀
1
La courbure de Gauss (équation caractéristique)
𝜕 𝒙1
obtenue en laissant déduire les identités du symbole
de Christofel
[ (
1
𝑨1
𝜕 𝑨2
𝜕 𝒙1 )] +
1
𝜕 𝒙2 [ (
1
𝑨2
𝜕 𝑨1
𝜕 𝒙2 )] =
𝑨1 𝑨2
𝑲

𝑨𝜶 𝑨𝜷 𝜶
La troisième forme quadratique fondamentale 𝒄 𝜶𝜷 = =𝒃 𝜷 𝒃𝝁𝜷
Selon R, NZENGWA, 𝑲
 Matériels et Méthodes
Méthodes
Théorie bidimensionnel des coques
Les théories de coques 3-D sont basées sur des hypothèses qui sont dans la majorité des cas en physique ou en mécanique faites pour
pouvoir modéliser ou simplifier les calculs. Dans le cas des théories sur les coques A. KHALED (31) and M. S. Qatu; R. W.
Sullivan; W. Wang (32)ont présenté une représentation mathématique de certaines théories afin d'inclure l'influence des hypothèses
sur les première et deuxième équations fondamentales. Loin de ces considérations, le modèle bidimensionnel de la N-T pour les
coques épaisses élastiques linéaires a été déduit du problème tridimensionnel sans hypothèses spécifiques de nature géométrique ou
mécanique R., NZENGWA; B.H., TAGNE (2). Les équations à deux dimensions sont déduites en appliquant une analyse
asymptotique à une famille d’équations variationnelles obtenues à partir d’une coque abstraite mise à l’échelle par échelonnage
multiple des équations tridimensionnelles initiales.
Précédemment nous donnions la définition géométrique de la surface moyenne de la coque. Nous définirons d’abord une coque
comme un solide tridimensionnel dans l’espace avant de présenter quelques théories célèbres qui sont utilisées.

Surface moyenne S Nous choisissons dans ce cadre, la coque à épaisseur constante

Dont une coque complète peut être définie comme:

ℎ ℎ
Ω= ൜𝑀∈ℝ3,𝑂𝑀= 𝜑ሺ 𝑥1,𝑥2ሻ + 𝑧𝑎3,𝜑ሺ 𝑥1,𝑥2ሻ ∈𝑆,− ≤ 𝑧≤ ൠ
M

2 2
𝒉
𝟐 ⃗3
𝑒
𝒉
⃗
𝜑
−
𝟐
O
⃗2
𝑒
⃗1
𝑒
Figure : Coque 3D (à gauche) et surface moyenne de la coque ( à droite)
 Matériels et Méthodes
Méthodes
Théorie bidimensionnel des coques

⃗2 Tenseur métrique tridimensionnel

3𝑎
M
⃗
𝑎
𝑔 𝛼𝛽 =𝑔 𝛼 . 𝑔 𝛽 =𝜇𝛼𝜆 𝜇 𝜆𝛽 𝑎 𝜆𝛾=( 1 − 𝑥 3 𝐻 + ( 𝑥3 ) 𝐾 ) 𝑎𝛼𝛽
2
S
m
⃗1
𝑎
Tenseur mixte du changement de surface
⃗3
𝑒
⃗2
𝑒 Ω
e Et sa réciproque vu
⃗1o
𝑒 qu’il est inversible
Domaine de la coque Ω

et
Rapport entre l’espace métrique et la surface métrique
𝝆 ( 𝒙 3 ) =1 − 𝒛𝑯 + 𝒛 2 𝑲
𝟏 𝟏 𝟏 𝟏 𝟏 𝟏 𝟏
𝑯= ൬ + ൰ = 𝒃𝜶𝜶 = 𝒕𝒓(𝒃) 𝑲= ×
𝟐 𝑹𝟏 𝑹𝟐 𝟐 𝟐 𝑹𝟏 𝑹𝟐
Avec Et

Courbure moyenne Courbure de Gauss

on peut négliger la courbure de Gauss dans la théorie des coques minces en raison de sa faible épaisseur (en supposant que la
Base locale et membranaire métrique de la surface moyenne est presque égale à celle de toute autre surface le long d'une épaisseur de coque)
 Matériels et Méthodes
Méthodes
Théorie bidimensionnel des coques
Tenseur de déformation

𝜀= 𝜀 𝑖𝑗 . (𝑔 𝑖 𝑗
𝑔 ) 𝜀 𝛼𝛼 =𝜇𝛼 ( ∇𝛼 ~
𝑢 𝜆 −𝑏 𝜆𝛼 ~
𝜆
𝑢3 ) 𝜀 𝛼𝛼 =𝜇𝛼 ~
𝜆
𝑢 𝜆 ,𝛼 +
(
𝐴𝜆 , 𝛽 ~
𝐴𝛽
𝑢 𝛽 − 𝑏 𝜆𝛼 ~
𝑢3
)
Les déformations de cisaillement
transverse
Sont calculées et obtenues comme suit:

Modèle de coque Kirchhoff-Love

Sur la base des hypothèses que Kirchhoff a introduites dans le but de dériver la
théorie d'une plaque mince, Love 0. C. Zienktewicz and Y K Cheung (35) a été 1.La coquille est mince.
le premier à dériver un ensemble d'équations de base décrivant le comportement
d'une coque mince et élastique.. 2. Les contraintes et les déplacements sont suffisamment faibles pour que les quantités de
Généralement appelée première approximation de Love, cette théorie classique
de la petite déformation d’une coque mince est basée sur les hypothèses magnitude de second ordre et d'ordre supérieur dans les relations contrainte-déplacements
suivantes, également connues sous le nom d’hypothèses de Kirchhoff-Love puissent être négligées par rapport aux termes du premier ordre.

3. La contrainte normale transversale est faible par rapport aux autres composantes de la
contrainte normale et peut être négligée.
4. Une normale à la surface de référence avant la déformation reste droite et normale à la surface
configuration initiale et configuration en déformation
de référence déformée et ne subit aucune extension.
 Matériels et Méthodes
Méthodes
Théorie bidimensionnel des coques
Modèle de coque Nzengwa-Tagne
Un déplacement de la surface moyenne
Cette théorie des coques remonte à 1999, où sans hypothèse ad hoc sur les déformations
des coques ni sur la texture des matériaux, un modèle bidimensionnel des coques
~
𝑢 ( 𝑥1 , 𝑥2 , 𝑧 )
élastiques linéaires est déduit du problème tridimensionnel d’une coque d’épaisseur [26].
À partir de différentes échelles sur les composantes tangentes et normales du composantes de la contrainte de cisaillement transverse
déplacement largement utilisées dans les travaux de Ciarlet et al. [95], le déplacement
limite semble être un déplacement de Kirchhoff-Love d'un type différent 𝜀33 ¿ 𝑢 3 ,3 ( 𝑥 1 , 𝑥2 , 𝑧 ) =~
𝑢3 , 3 ( 𝑥1 , 𝑥 2 , 𝑧 ) =0 ⇒ ~
𝑢 3 ( 𝑥1 , 𝑥 2 , 𝑧 ) =𝜉 3 ( 𝑥1 , 𝑥
La déformation de cisaillement
2 𝜀 𝛼 3 =𝜇
𝜆 ~
𝑢 +~
𝑢3 , 𝛼 +𝑏 𝛼 ~
𝜆
𝑢𝜆
𝛼 𝜆, 3

La solution générale qui satisfait l’équation se trouve dans Nzengwa et al [26] et consiste à Le déplacement limite trouvé dans le modèle en coque de N-T contient le
déplacement classique de Kirchhoff-Love
trouver dans la base qui satisfait

Le tenseur de déformation déduit de

la cinématique des coques N-T par
rapport à l'état de déformation plan se
lit
, sont les tenseurs de déformation de la membrane e
𝜺 𝜶𝜷 ( 𝒖 ) = 𝒆𝜶𝜷 ( ~
𝒖 ) − 𝒛 𝒌𝜶𝜷 ( ~ ~
𝟐
𝒖)+ 𝒛 𝑸 𝜶𝜷 ( 𝒖 )
changement de tenseur de courbure et le changement de la
troisième forme fondamentale
 Matériels et Méthodes
Expérimentation

Nous allons expérimenter les théories énoncées plus haut à travers les essais de compression à réaliser sur certains tubes ronds et autres
viroles. Les essais consistent à appliquer une charge sur une virole dont la base est soudée. De ce modèle nous pourrons observer les
déplacements (conforme au champ de déplacement selon Grognec, Philippe Le (17) ).

Etant donné que nous sommes en compression simple nous pouvons dire que les composantes du tenseur tangent matériel peuvent être
explicités. En effet on considère un cylindre circulaire de longueur l, de rayon R et d’épaisseur t. l’une des extrémités du cylindre L’une
des extrémités du cylindre reste en contact avec le plan x = 0 (les déplacements radiaux et les rotations sont libres) tandis qu’en x = l
s’applique une contrainte de compression axiale Σ répartie uniformément. On prendra λ = -Σ > 0. Enfin, l’écrouissage est pris purement
isotrope linéaire, de module h constant.
Charge critique et mode de bifurcation
Valeur critique sous la forme
mode de bifurcation sinusoidal
𝝀
𝒂( 𝒏)
𝒄 =𝜶
𝒕 2 𝒏2 𝝅 2
12 𝒍 ²
+ 𝜸−
𝜷2
𝜶 ( ) 1 𝒍²
𝑹 ² 𝒏2 𝝅 2
Le mode de bifurcation recherché s’obtient en minimisant avec pour un cylindre
suffisamment long, la minimisation sur donne une bonne approximation de la

√
valeur critique
𝑬𝒕 4𝒉
𝝀 𝒂𝒄 =2 𝝀𝒍𝒄 =
𝑹 3 [ ( 5 − 4 𝝂 ) 𝑬 +4 𝒉 ( 1−𝝂 2 ) ]
 Matériels et Méthodes
Expérimentation
Description du dispositif

Grandeurs à mesurer Résultats attendus

Effort à la compression : valeur donnée Courbes de déformation/déplacement

en kN
Courbe Contrainte-déformation
Contrainte à la compression : valeur
donnée en Mpa Normes
NF EN 12-390-3
Déformation/Déplacement à l’aide des
comparateurs Materiaux
Presse hydraulique Acier S235
charge maximale applicable est de 4500kN
Principe: insérer une éprouvette dans un bâti dans lequel une presse Capteurs
hydraulique qui effectue un mouvement de translation (du haut vers le bas) ; Deux Comparateurs étalonnés par le
cette presse applique sur l’éprouvette une charge et la lecture se fait dès qu’il
y a contact entre la presse et l’éprouvette LABOGENIE (précision 0.01mm)
 Matériels et Méthodes
Expérimentation
Description du dispositif

Installation du dispositif de mesure

Plat en L fixé sur le bâti de la presse

Fixation des comparateurs par

l’entreprise d’un ruban adhésif

Nous avons eu a effectuer les

chargements sur les viroles, et nous
avons fait les mesures des déformations,
grâce à un chronomètre nous avons
saisit le temps mis pour les déformations

Presse hydraulique
charge maximale applicable est de 4500kN
Principe: insérer une éprouvette dans un bâti dans lequel une presse
hydraulique qui effectue un mouvement de translation (du haut vers le bas) ;
cette presse applique sur l’éprouvette une charge et la lecture se fait dès qu’il
y a contact entre la presse et l’éprouvette
 Matériels et Méthodes
Expérimentation
Description des echantillons
N° Virole Image eprouvette Dimensions N° Virole Image eprouvette Dimensions
(eprouvette) (eprouvette)
Chargement sur éprouvette
metallique
V1 Hauteur : 330 mm Hauteur : 280 mm
V2
Diametre exterieur :
P= charge presse
120 mm
Diametre exterieur :
125 mm

h
d

N° Virole Image eprouvette Dimensions N° Virole Image eprouvette Dimensions

D (eprouvette) (eprouvette)

Hauteur : 190 mm Hauteur : 330 mm

V3 V4
Diametre exterieur : Diametre exterieur :
Lieu des essais 125 mm 50 mm
Laboratoire BHYGRAPH
GEOTECHNIQUE S.A
 Matériels et Méthodes
Expérimentation
Description des échantillons
N° Virole Image eprouvette Dimensions N° Virole Image eprouvette Dimensions
(eprouvette) (eprouvette)

Hauteur : 190 mm Hauteur : 280 mm

V5 V2
Diametre exterieur : Diametre exterieur :
50 mm 50 mm

Num hci(mm) dc(mm) hc/dc Type d'élancement Charge

Virole appliquée

V1 330 125 2,64 Silo elancé 455

V2 280 120 2,33 Silo elancé 455
V3 190 125 1,52 Silo d'élancement 455
intermédiaire

V4 330 50 6,60 Silo tres elancé 455

V5 190 50 3,80 Silo elancé 455
V6 280 50 5,60 Silo tres elancé 455
Représentation des élancements des silos Tableau récapitulatif des échantillons
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Déformations en fonction du temps
Virole 1 Virole 2
Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée
V1 330 125 2 2,64 Silo elancé V2 280 120 2 2,33 Silo elancé 455

Déformation en fonction du temps Déformation en fonction du temps

90 90

80 80

70 70

60 60

50 50

40 40

30 30

20
20
10
10
0
0 1 4 7
1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 67 70 73 76 79 82 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 67 70 73 76 79 82

Temps (s) Deformation (mm)

Temps (s) Deformation (mm)
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Déformations en fonction du temps
Virole 3 Virole 4
Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée
V3 190 125 1,5 1,52 Silo d'élancement intermédiaire 455 V4 330 50 2 6,60 Silo tres elancé 455

Déformation en fonction du temps Déformation en fonction du temps

60 90

80
50
70

40 60

50
30
40

20 30

20
10
10

0 0
1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 67 70 73 76 79 82

Temps (s) Deformation (mm) Temps (s) Deformation (mm)

 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Déformations en fonction du temps
Virole 5 Virole 6
Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée
V5 190 50 2 3,80 Silo elancé 455 V6 280 50 2 5,60 Silo tres elancé 455

Déformation en fonction du temps Déformation en fonction du temps

90 80

80 70
70
60
60
50
50
40
40
30
30
20
20

10 10

0 0
1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 67 70 73 76 79 82 1 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58 61 64 67 70 73

Temps (s) Deformation (mm) Temps (s) Deformation (mm)

 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Courbes contraintes déformations
Virole 1 Virole 2
Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée
V1 330 125 2 2,64 Silo elancé V2 280 120 2 2,33 Silo elancé 455

𝜀 (𝑚𝑚)
𝜀
Courbe Contrainte-déformation
Courbe Contrainte-déformation
(𝑚𝑚)
3.5
25
3
20
2.5

2 15

1.5 10

1
5

𝜎 𝜎
0.5

0 (𝑀𝑃𝑎) 0 (𝑀𝑃𝑎)

Series1
Series1
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Déformations en fonction du temps
Virole 3 Virole 4
Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée
V3 190 125 1,5 1,52 Silo d'élancement intermédiaire 455 V4 330 50 2 6,60 Silo tres elancé 455

𝜀 𝜀 Courbe Contrainte-déformation
(𝑚𝑚)
Courbe Contrainte-déformation
(𝑚𝑚)
14 30

12 25

10
20
8
15
6
10
4

𝜎
5
2

𝜎 (𝑀𝑃𝑎
(𝑀𝑃𝑎 )
0 0

Series1 Series1
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Déformations en fonction du temps
Virole 5 Virole 6
Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée Num Virole hci(mm) dc(mm) Ep (mm) hc/dc Type d'élancement Charge appliquée
V5 190 50 2 3,80 Silo elancé 455 V6 280 50 2 5,60 Silo tres elancé 455

𝜀 𝜀
Courbe Contrainte-déformation
(𝑚𝑚)
Courbe Contrainte-déformation
40 (𝑚𝑚)
30

35 25
30
20
25

20 15

15 10
10
5

𝜎 𝜎
5

(𝑀𝑃 (𝑀𝑃𝑎
0

𝑎) )
0

Series1
Series1
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Regroupement statistique
Modele Regroupement des valeurs statistiques
J.El bahaoui (5) nous allons nous rendre compte qu’il faut mettre en Virole tres élancée (silo tres Virole élancée (silo tres Virole avec elancement
élancé) élancé) intermédiaire (silo avec
avant les aspects liés aux défaut localisés sur les coques cylindriques.
elancement intermédiaire)
En effet le modele fiabiliste qu’il propose est défini par G(X) = R(X) –
V4, V6 V1, V2, V5 V3
S(X) où R(X) est la ressources disponible et S(X) la sollicitation.
Le modele elements finis va nous permettre d’étudier la variation de la
charge critique dans le problème du flambage des coques cylindriques
minces. Seulement dans le cadre de la présente étude, il nous sera
difficile d’utiliser ces méthodes, mais juste de les constater sur les
résultats obtenus après expérimentation.
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats
Observations sur echantillons
Images Observations Déformation en tete observée aussi sur virole 2 Déformation en tete observée aussi sur virole 3.
Cette déformation est relativement faible,
Déformation en tete de virole observée. La valeur
comparée au modele numérique.
est présentée plus haut dans les tableaux qui
précedent.

Virole V1 Virole V2 Virole V3

Déformation en tete observée aussi sur virole 4. Déformation en tete observée aussi sur virole 5. Déformation en tete observée aussi sur virole 6.
Ce modele étant élancé, l’on s’attendrait à une
présence du flambage au niveau du tronc. Les
résultats obtenus sur l’éprouvette sont importants
comparé au modele numérique.

Virole V4 Virole V5
 RESULTATS ET DISCUSSSIONS
Expérimentation
Description du modèle numérique

Type de maillage:
Logiciel utilisé: - La méthode de Coons Definition
Cette méthode bien qu’ancienne permet d’obtenir des
ROBOT STRUCTURAL résultats (maillage) des volumes. D’apres Cuillière, J C
(36) qui présente cette méthode comme la plus efficiente
ANALYSIS pour la modélisation des objet 3D dans le cadre de la CAO
(Conception Assistée par Ordinateur). Il démonter a cet
effet que la génération du maillage revêt trois aspects
différents :
Le remaillage adaptatif
Le contrôle de l’erreur géométrique lors de la discrétisation
La discrétisation à priori.

Parametrage du maillage:

Pour une charge appliquée de

3 699 187.00 daN/m²
Taille des éléments fixée à 1cm
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats du modèle numérique
Déformations en fonction du temps
Virole D (mm) Virole D (mm) Virole D (mm)

V1 V2 V3

Nous avons un max de 3.9mm Nous avons un max de 4.06mm Nous avons un max de 51.169mm

Virole D (mm) Virole D (mm) Virole D (mm)

V5 V6
V4

Nous avons un max de 6.515mm Nous avons un max de 1.848mm Nous n’avons pas de déformation obtenu
 RESULTATS ET DISCUSSSIONS
Expérimentation
Résultats du modèle numérique
Analyse des résultats obtenus

Cette déformation observée est

conforme en partie aux observations
faites par J.El bahaoui (11), qui lui
met en relation une intégration
numérique sur l’épaisseur de la
coque et utilise ainsi une
linéarisation cohérente de la
cinématique de la coque et du
modele elasto-plastique qui conduit
au modele de NEWTON-
RAPHSON

Image des déformations obtenues par PRANK (1997)

 CONCLUSIONS ET PERSPECTIVES

But Limites des lois usuelles de

Rendre au mieux le comportement des comportement
structures intégrant dans leurs géométries
 le modèle de Kirchhoff-Love est directement difficile à
des courbures et ainsi d’en optimiser la
conception et le dimensionnement, a mettre en œuvre sous forme numérique et néglige le
travers les modèles numériques et
cisaillement transverse ;
physiques utilisés dans le cadre de notre
travail,
 pour des éléments finis conformes le modèle de Reissner –
Mindlin conduit à des instabilités numériques dès que
Difficultés et perspectives l’épaisseur devient trop petite.
 les défauts d’ordre structurel,
 les points d’application des charges (les charges étaient uniquement appliquées en tete de poteaux),
 le choix des matériaux (en fonction des limites d’élasticité),
 l’ajustement du modele numérique dans un logiciel dédié.
 La loi de comportement associée aux résultats obtenus.
Ces aspects feront l’objet de recherches complémentaires à la suite de ce travail
 Merci pour votre attention!!!

VOS OBSERVATIONS SONT VIVEMENT

SOUHAITEES