Conception de réservoirs à hydrogène composites
Conception de réservoirs à hydrogène composites
MÉMOIRE DE MASTER
Domaine : Sciences et Techniques
Filière : Génie Mécanique
Spécialité : Construction Mécanique
Jury :
Merci à tous et que Dieu vous bénisse et vous garde, si Dieu le veut.
Sommaire
Sommaire
Introduction …………………………………………………………………………... 4
I.1 L’hydrogène gazeux ……………………………………………………………… 5
I.2 La Production d’hydrogène ………………………………………………………. 5
I.3 Intérêt énergétique …………………………………………………………………. 6
I.4 Stockage d’hydrogène ……………………………………………………………... 6
I.4.1 Stockage gazeux …………………………………………………………….. 6
I.4.2 Stockage liquide …………………………………………………………….. 7
I.4.3 Stockage Solid……………………………………………………………….. 7
I.5 Stockage et transport de l’hydrogène………………………………………………. 8
I.6 Stockage du l’hydrogène comprimé ……………………………………………….. 8
I.6.1 réservoirs sous pression…………………………………………………………... 8
I.6.1.1Les réservoirs de « type I »………………………………………………….. 9
I.6.1.2Les réservoirs de « type II »………………………………………………… 9
I.6.1.3 Les réservoirs de « type III »……………………………………………….. 9
I.6.1.4 Les réservoirs de « type IV »……………………………………………….. 10
Introduction ……………………………………………………………………………. 13
II.1Définition………………………………………………………………………….... 13
II.2Classification des matériaux composites…………………………………………… 14
II.3 Composants ……………………………………………………………………….. 14
II.3.1 Renforts……………………………………………………………………….. 15
II.3.1.1 Différents types de renfort………………………………………………… 16
I
II.[Link] Fibres de verre………………………………………………………… 16
II.[Link] Fibres de carbone……………………………………………………… 16
II.[Link] Fibres aramides………………………………………………………... 18
II.[Link] Fibres céramiques……………………………………………………... 18
II.[Link] Fibre Naturelle………………………………………………………… 18
II.3.1.2 Caractéristiques mécaniques des fibres……………………………………. 19
II.3.2 La matrice……………………………………………………………………... 19
II.3.2.1 Résines thermodurcissables……………………………………………….. 20
II.3.2.2Résines thermoplastiques………………………………………………….. 20
II.3.2.3 Les matrices naturelles……………………………………………………. 21
II.3.2.4 Résines époxydes…………………………………………………………. 21
II.3.2.5 Les résines thermostables…………………………………………………. 21
II.3.3 Les charge……………………………………………………………………… 21
II.3.4 Les additifs…………………………………………………………………….. 22
II.3.5 L’interface……………………………………………………………………... 22
II.4 Stratifiés……………………………………………………………………………. 22
II.5 Sandwichs………………………………………………………………………….. 23
II.6 Procédé de moulage par enroulement filamentaire………………………………... 24
II.7Comportement mécanique de la monocouche……………………………………… 26
II.7.1 Loi de comportement………………………………………………………….. 26
II.7.2 Loi de mélange……………………………………………………………....... 28
II.7.3 Expression des relations d’élasticité hors axes principaux…………………… 29
Introduction…………………………………………………………………………….. 32
III-1 structure des réservoirs de type IV ………………………………………………. 32
III.1.1 Architecture d’un réservoir type IV…………………………………………….. 32
III.1.1.1 L’enveloppe étanche (ou liner)……………………………………………. 32
III.1.1.2 La structure résistante……………………………………………………... 33
[Link] protection externe……………………………………………………… 33
III.2 Fabrication du Réservoir type IV par enroulement filamentaire…………………. 33
III.3 Contrôle de réservoir fabriqué……………………………………………………. 34
III-4 Modèle d’enroulement filamentaire (netting theory)…………………………….. 35
I
III.4.1 Conception de la dôme………………………………………………………. 38
Introduction…………………………………………………………………………… 41
IV .1 Détermination des paramètres du réservoir (MATLAB)……………………….. 41
IV.1.1 Paramètres d’entrée………………………………………………………… 41
IV.1.2 Paramètres de sortie………………………………………………………… 41
IV.1.3 Organigramme……………………………………………………………… 41
IV .2 Modélisation de la géométrie du réservoir (SOLIWORKS)……………………. 43
IV .3 Simulation par la méthode des éléments finis (ANSYS)……………………….. 43
IV.4 Application………………………………………………………………………. 43
IV.4.1Donnée………………………………………………………………………. 43
IV.4.2 Résultats……………………………………………………………………... 44
IV.4.2.1 Paramètres du réservoir (MATLAB)…………………………………… 44
IV.4.2.2 Modèle géométrique (SOLIDWORKS)………………………………... 46
IV.4.2.3 Modèle par la méthode des éléments finis (ANSYS) 46
IV.[Link] Maillage……………………………………………………………. 46
IV.[Link] Matériaux…………………………………………………………... 47
IV.[Link] condition aux limites ………………………………………………. 47
IV.4.2.3.4charge……………………………………………………………….. 47
IV.[Link] Nombre de couche…………………………………………………. 47
IV.4.2.3.6Les graphes de ANSYS…………………………………………….. 48
Références Bibliographiques…………………………………………………………... 52
I
Liste des Figures
CHAPITRE I : Généralités sur l’hydrogène gazeux
Figure(I.1) : la pollution………………………………………………………………. 4
comprimé
II
Figure (II.10) : les différant types de sandwiche……………………………………… 24
Figure (II. 12) : Réservoir sous pression constitué d'un tuyau d'acier horizontal……... 26
II
Liste des tableaux
Thermoplastiques
III
Nomenclature
IV
الملخص
، ألنه يمكن تخزينه في أربعة أنواع من خزانات الضغط، الغاز الغازي هو الخيار األفضل.يعتبر غاز الهيدروجين بديلا جيداا للوقود األحفوري
على المواد.وخاصة النوع الرابع من الخزانات (موضوع دراستنا) بفضل األداء الفني والقدرة التنافسية من حيث التكلفة وتطويره في هذا العمل
. حل لتصنيع الخزانات.المركبة كحل واعد
الهدف من هذا العمل هو دراسة الخزان من النوع الرابع ليؤدي إلى تصميم األجزاء المختلفة المكونة لخزان الهيدروجين في الحالة الغازية عند
سيتم استخدام نظرية لف الخيوط (نظرية المعاوضة) لتحديد زوايا االتجاه وسماكة.) قاعدة معدنية، غلف مركب،الضغط العالي (بطانة بوليمر
(ANSYS). سيتم استخدام هذه البيانات التي تم الحصول عليها لمحاكاة السلوك الميكانيكي للخزان باستخدام برنامج.ثنايا الطبقات
، وسمك وعدد طبقات غطاء الخزان، وهي زاوية اتجاه األلياف، لتحديد المسند المتعلق بأوعية الضغطMATLAB تم تطوير برنامج على
. الخزان.فضلا عن الشكل الهندسي لقبتين
.ANSYS النمذجة بطريقة العناصر المحدودة تتم باستخدام برنامج.)SOLIDWORKS) تم تصميم الشكل الهندسي باستخدام البرنامج
Résumé
Le gaz hydrogène est considéré comme une bonne alternative aux combustibles fossiles. Gaz gazeux est la
meilleure option, car elle peut être stockée dans les quatre types de réservoirs sous pression, notamment le quatrième
type de réservoir (objet de notre étude) grâce à des performances techniques et une compétitivité des coûts et de le
développer dans ce travail sur les matériaux composites comme une solution prometteuse. Solution pour la fabrication
de réservoirs.
L’objectif de ce travail est l’étude de réservoir de type IV pour aboutir à la conception des différentes
parties constituants le réservoir de stockage d’hydrogène à l’état gazeux à haute pression (liner polymère, coque
composite, embase métallique). La théorie d’enroulement filamentaire (netting theory) sera utilisée la détermination des
angles d’orientation et l’épaisseur des plis de la stratification. Ces données obtenues seront utilisées pour simuler le
comportement mécanique du réservoir à l'aide du programme (ANSYS) .
Un programme sur MATLAB a été développé pour déterminer la donnée relative au réservoir sous pression, qui est
l'angle de direction des fibres, l'épaisseur et le nombre de couches du couvercle du réservoir, ainsi que la géométrie
forme des deux dômes du réservoir.
La géométrie a été modélisée avec SOLIDWORKS. La modélisation par la méthode des éléments finis est opérée en
utilisant le logiciel ANSYS.
Abstract
Hydrogen gas is considered a good alternative to fossil fuels. Gaseous gas is the best option, because it can be stored in
the four types of pressure tanks, especially the fourth type of tank (object of our study) thanks to technical performance
and cost competitiveness and to develop it in this work. On composite materials as a promising solution. Solution for the
manufacture of tanks.
The objective of this work is the study of type IV tank to lead to the design of the different parts constituting the hydrogen
storage tank in the gaseous state at high pressure (polymer liner, composite shell, metal base). The filament winding
theory (netting theory) will be used to determine the orientation angles and the thickness of the layers of the stratification.
These data obtained will be used to simulate the mechanical behavior of the reservoir using the program (ANSYS).
A program on MATLAB has been developed to determine the datum relating to the pressure vessel, which is the angle of
direction of the fibers, the thickness and the number of layers of the tank cover, as well as the shape geometry of the two
domes of the tank. .
The geometry was modeled with SOLIDWORKS. Modeling by the finite element method is carried out using the ANSYS
software.
Introduction générale
Introduction général
Introduction générale
Pour que l’hydrogène puisse réellement devenir le vecteur énergétique de demain, il faut qu’il
soit disponible à tous moment. Dans ce contexte mettre au point des modes de transport, de stockage
et de distribution efficaces, représente un enjeu crucial. Aujourd’hui, le stockage de l’hydrogène reste
la solution la plus satisfaisante dans l’attente de spécifications techniques et économiques des
industriels qui soient stables. Il existe trois stratégies différentes de stockage selon l’état de
l’hydrogène (solide, liquide et gazeux). Son conditionnement sous forme gazeuse reste néanmoins
l’option la plus prometteuse. Le stockage sous forme de gaz comprimé à de très hautes pressions
apparaît à l’heure actuelle comme la solution présentant le meilleur compromis en termes de densité
massique et volumétrique. [1]
L'hydrogène peut être stocké dans les quatre types de récipients sous pression. Le choix du
stockage se fait en fonction de l'application finale qui nécessite un compromis entre performances
techniques et compétitivité coût. L'H2 en tant que gaz industriel est stocké dans des réservoirs de type
I, dont la pression est de 150 à 300 bars (généralement 200 bars). Ce sont les réservoirs à haute
pression les plus répandus aujourd'hui et les moins chers. Lorsque seules des pressions plus élevées
sont requises - principalement pour des applications stationnaires - les réservoirs de type II sont
préférés. Les cuves de type III et de type IV sont destinées à des applications portables, pour lesquelles
le gain de poids est essentiel. Cependant, ces réservoirs sont beaucoup plus chers. [2].
1
Introduction général
Type II : Récipient sous pression constitué d'un cerceau de revêtement métallique épais
enveloppé d'un composite fibre-résine.
Type III : Récipient sous pression constitué d'un revêtement métallique entièrement
enveloppé d'un composite fibre-résine.
type IV Récipient sous pression constitué d'un revêtement polymère entièrement
enveloppé d'un composite fibre-résine Le port est métallique et intégré dans la structure
(boss).
Le développement de nouveaux matériaux est soumis aux besoins des utilisateurs. Les matériaux
composites stratifiés se présentent actuellement comme une solution prometteuse pour la fabrication des
réservoirs et des tubes répondant aux exigences environnementales et à la contrainte du poids. Le type IV de
réservoir à une structure en matériaux composites. Ces structures de renforcement constitué de fibres (verre,
aramide, carbone) et d’une résine (thermodurcissable ou thermoplastique) [3].
L’objectif de ce travail est l’étude de réservoir de type IV pour aboutir à la conception des
différentes parties constituants le réservoir de stockage d’hydrogène à l’état gazeux à haute pression
(liner polymère, coque composite, embase métallique). La théorie d’enroulement filamentaire
(netting theory) sera utilisée la détermination des angles d’orientation et l’épaisseur des plis de la
stratification. Ces données obtenues seront utilisées pour simuler le comportement mécanique du
réservoir à l'aide du programme (ANSYS).
Le premier chapitre : Où nous abordons dans ce chapitre l'étude de l'hydrogène gazeux à travers
la définition, son intérêt énergétique, les solutions nécessaires de production et comment le
stocker et le transporter.
Le deuxième chapitre : Dans ce chapitre, nous abordons les généralités des matériaux
composites, à travers la définition des matériaux composites, leurs composants, et leur ingénierie
(structure).
Le troisième chapitre : Dans ce chapitre, nous présentons les détails de la modélisation du
comportement mécanique de la structure du réservoir. Les logiciels utilisés dans notre étude
seront aussi présentés dans ce chapitre, en expliquant le travail de ces logiciels et de la simulation.
Le quatrième chapitre : Les modélisations et simulations du réservoir type 4 en utilisant les
logiciels MATLAB, SOLIDWORKS et ANSYS et extraire les résultats.
2
CHAPITRE I
Généralité sur
l’hydrogène gazeux
CHABITRE I/ Généralité sur l’hydrogène gazeux
CHABITRE I
Généralité sur l’hydrogène gazeux
Introduction
Confronté à l’augmentation des besoins énergétiques mondiaux (doublement des besoins d’ici
à 2050), à la raréfaction des ressources fossiles et à la nécessité de réduire les émissions de gaz à effet
de serre issues de l’utilisation de ces ressources, il convient de développer et diversifier nos ressources
en énergie. Partant de ce constat, on assiste depuis quelques années à la montée en puissance
d’énergies et de ressources dites alternatives. Bien sûr, toutes ces ressources et les technologies qui y
sont associées, n’en sont pas toutes au même stade de maturité : certaines n’en sont qu’au stade de
recherche alors que d’autres en sont au stade de développement, de diffusion voire même au stade
d’industrialisation à grande échelle.
On peut citer à ce titre :
– les ressources actuelles et à court terme : éolien, solaire, biomasse, stockage électrique. . .
– les ressources à moyen terme (20 ans) : hydrogène notamment ;
– les ressources à long terme (50 ans) : fusion nucléaire
4
CHABITRE I/ Généralité sur l’hydrogène gazeux
Le dihydrogène est la forme moléculaire de l’élément hydrogène qui existe à l’état gazeux
aux conditions normales de température et de pression. La molécule comporte deux
atomes d’hydrogène ; sa formule chimique est H2 . Le dihydrogène est également
appelé « hydrogène moléculaire. » ou, à l’état gazeux, « gaz (d’) hydrogène. »
5
CHABITRE I/ Généralité sur l’hydrogène gazeux
La voie gazeuse est la plus simple et la plus répandue des technologies de stockage.
L’hydrogène est contenu dans des réservoirs sous pression. Plus la pression est élevée, plus la quantité
d’hydrogène stocké est grande. Les réservoirs actuellement développés fonctionnent avec des
pressions maximales aux alentours de 300-350 bar (30-35 MPa) voir 480 bar. Le travail actuel
consiste à créer des réservoirs gazeux sous 700 bar afin d’améliorer la capacité volumique, principal
défaut de cette technologie. [8/9]
6
CHABITRE I/ Généralité sur l’hydrogène gazeux
Le stockage solide semble être la voie la plus prometteuse pour l’avenir. Les capacités
volumiques sont très intéressantes puisque supérieures à celles du stockage liquide. Certaines formes.
7
CHABITRE I/ Généralité sur l’hydrogène gazeux
Le réseau de gaz naturel peut aussi contenir de l’hydrogène jusqu’à 20% en volume mais cela
nécessite de le purifier en aval du réseau. Notons qu’il faut dix fois plus d’espace pour stocker de
l’hydrogène gazeux par rapport à l’essence.
L'hydrogène à été stocké et distribué sous pression il y a plusieurs années, avec des
bouteilles ou assemblages de bouteilles cylindriques, en acier, gonflées à 20 ou 25 MPa (types
I et II). L’inconvénient le poids qui résulte de l’utilisation d’aciers à bas niveaux de contraintes
pour éviter les problèmes de fragilisation par l’hydrogène. La situation a radicalement changé
avec l'apparition de la technologie des réservoirs composites dits de type III ou IV. Leur
principe de base est de séparer les deux fonctions essentielles que sont
8
CHABITRE I/ Généralité sur l’hydrogène gazeux
L’étanchéité et la tenue mécanique pour les gérer l'une indépendamment de l'autre. Dans ce type
de réservoir on associe une vessie en résine (thermodurcissable ou thermoplastique) à une
structure de renforcement constituée de fibres (verre, aramide, carbone) permettent de travailler
à des pressions beaucoup plus élevées tout en réduisant la masse et en évitant les risques de
rupture explosive en cas d’agressions externes sévères. C’est ainsi que 70 MPa est pratiquement
devenu le standard actuel. Du fait des lois de compressibilité, cette valeur semble être une limite,
même si des réservoirs de pression nominale égale ou supérieure à 100 MPa ont été réalisés à
titre expérimental. [13]
Les réservoirs de stockage gazeux hyperbare peuvent être classés en quatre catégories selon
leurs matériaux constitutifs mineurs
9
CHABITRE I/ Généralité sur l’hydrogène gazeux
10
CHABITRE I/ Généralité sur l’hydrogène gazeux
Les réservoirs dans le cadre de cette étude se limitent aux réservoirs de type IV constitué, d’une part d’un liner
en polymère assurant la fonction d’étanchéité et d’autre part, d’un matériau composite assurant la fonction
structuration à la pression interne
11
CHAPITRE II
Généralités sur
Matériaux
Composites
CHAPITRE II / Généralités sur Matériaux Composites
CHAPITRE II
Généralités sur Matériaux Composites
Introduction
Pour l’utilisation intensive de matériaux composites, nous présentons dans ce chapitre ces
matériaux ainsi que leurs caractéristiques générales et leurs composants. Ensuite, seront exposés
les avantages et les inconvénients des matériaux composites et les matériaux composites
structuraux, suivi par les procédé de fabrication des composites.
II.1Définition
Un matériau composite est constitué de l'assemblage de deux ou plusieurs matériaux de
natures différentes. Leur association est complémentaire est permet d'aboutir à un matériau dont
les performances recherchées seront supérieures à celles des composants pris séparément. Un
matériau composite est constitué dans le cas le plus général d'une ou plusieurs phases
discontinues réparties dans une phase continue. La phase discontinue, appelée renfort ou
matériau renforçant, est habituellement plus dure avec des propriétés mécaniques supérieures à
celles de la phase continue appelée matrice. [15]
13
CHAPITRE II / Généralités sur Matériaux Composites
14
CHAPITRE II / Généralités sur Matériaux Composites
15
CHAPITRE II / Généralités sur Matériaux Composites
Les fibres de carbone sont pour la plupart obtenues à partir d’une fibre précurseur de
type poly acrylonitrile (PAN). Ces fibres subissent une étape d’oxydation pour faire réticuler le
PAN et le rendre ainsi moins « fusible » en vue des traitements thermiques ultérieurs.
Elles sont ensuite soumises à une seconde étape, qui peut être une carbonisation à
1500°C ou à une graphitisation à 3000°C. Pendant cette étape, les fibres sont maintenues sous
tension afin de conserver l’orientation de la structure moléculaire.
Dans le premier cas, on obtient des fibres haute résistances ou HR, et dans le deuxième
cas, on obtient des fibres haut module ou HM (Tableau 3). Il existe également des fibres
16
CHAPITRE II / Généralités sur Matériaux Composites
module intermédiaire, appelée inter médiate modules ou IM. Elles peuvent ensuite
passer dans un électrolyte ou subir une oxydation pour créer des groupements polaires, destinés
à favoriser l’ennoyage. Il existe différents traitements de surface pour les fibres de carbone [22].
Les fibres de carbone continues sont disponibles sous diverses formes de mèches de
3000 à 24000 filaments, de fibres coupées, de rubans ou bien encore de tissus, qu’ils soient uni
ou multidirectionnels. Elles présentent des caractéristiques spécifiques élevées comme leur haut
module en traction. Elles ont également un coefficient de dilatation très faible, une conductivité
thermique élevée. En revanche, certains inconvénients demeurent comme leur résistance limitée
aux chocs, une faible résistance à l’abrasion.
17
CHAPITRE II / Généralités sur Matériaux Composites
Fibres aramides. Des renforts hybrides de type verre kevlar ou carbone kevlar sont
largement utilisés dans le domaine des loisirs (ski, raquette de tennis) [21]
Fibres naturelles
Les fibres naturelles sont des biomatériaux d’origine végétale ou animale, tissées, tricotées ou,
tressées. Elles ont été initialement utilisées pour la production de vêtements, tissus, tapis et cordages.
Actuellement, les fibres naturelles sont également utilisées dans l’industrie, dans les applications
automobiles, le bâtiment …etc. Les fibres naturelles permettent de réduire le coût et le poids dépièces
et d’améliorer les perspectives de recyclage et l’impact environnemental.
Les fibres naturelles peuvent être classées en trois grands groupes selon leur origine :
Les fibres végétales qui comprennent : des fibres provenant des poils séminaux de graines
18
CHAPITRE II / Généralités sur Matériaux Composites
(coton, kapok), des fibres libériennes extraites de tiges de plantes (lin, chanvre, jute, ortie),
des fibres extraites de feuilles (sisal), de troncs (palmier), d’enveloppes de fruits (noix de coco).
Les fibres animales
Les fibres minérales t
II.3.1.2 Caractéristiques mécaniques des fibres
Il existe différents types de fibres. Elles peuvent être scindées en deux groupes, les fibres à
haut module et les fibres à haute résistance. Les fibres à haut module ont une résistance faible et celles
à forte résistance ont un module faible. [21]
II.3.2 La matrice
Le rôle de la matrice est principalement de lier les fibres, mais aussi, d’assurer une
répartition spatiale est homogène du renfort. D’un point de vue mécanique, la matrice transmet
et répartie les efforts extérieurs vers le renfort. D’un autre côté, elle apporte à la structure une
tenue chimique et donner la forme désirée au produit [24], (figure II.7).
19
CHAPITRE II / Généralités sur Matériaux Composites
II.3.2.2Résines thermoplastiques
Les résines thermoplastiques ont des propriétés mécaniques faibles. Ces résines sont
solides et nécessitent une transformation à très haute température. Les polychlorures de vinyle
(PVC), les polyéthylènes, polypropylène, polystyrène, polycarbonate polyamide sont quelques
exemples de ces résines thermoplastiques. Les résines thermoplastiques classiquement
rencontrées sont PEEK, K3B.
De même que pour les résines thermodurcissables, les matériaux les plus performants
20
CHAPITRE II / Généralités sur Matériaux Composites
ont des caractéristiques mécaniques élevées et une masse volumique faible : ces dernières sont
présentées dans le tableau (II.4). [21]
21
CHAPITRE II / Généralités sur Matériaux Composites
pigments et colorants ;
agents anti
retraits ;
agents anti
ultraviolets ;
accélérateur ;
catalyseur.
II.3.5 L’interface
En plus de ces fibres et la matrice, il faut rajouter : une interface qui assure la
compatibilité renfort/matrice, qui transmet les contraintes de l’un à l’autre sans déplacement
relatif. Bonne adhérence en couche fine.
Des produits chimiques entre aussi dans la composition du composite, l’interphase etc...
qui peut jouer sur le comportement mécanique, mais n’interviennent pratiquement jamais dans
le calcul de structure composite. [20]
II.4 Stratifiés
Un stratifié est constitué d'un empilement de monocouches ayant chacun une orientation
propre à un référentiel commun aux couches et désigné comme le référentiel du stratifié.
En jouant sur l'ordre et l'orientation de ces couches. Il est possible d'adapter finement les
22
CHAPITRE II / Généralités sur Matériaux Composites
II.5 Sandwichs
Matériaux composés de deux semelles (ou peaux) de grande rigidité et de faible épaisseur
enveloppant une âme (ou cœur) de forte épaisseur et faible résistance.
L'ensemble forme une structure d'une grande légèreté. Le matériau sandwich possède une
grande légèreté en flexion et c'est un excellent isolant thermique.
23
CHAPITRE II / Généralités sur Matériaux Composites
24
CHAPITRE II / Généralités sur Matériaux Composites
25
CHAPITRE II / Généralités sur Matériaux Composites
Figure (II. 12) : Réservoir sous pression constitué d'un tuyau d'acier horizontal. [30]
26
CHAPITRE II / Généralités sur Matériaux Composites
𝜎5 0 0 0 0 𝑐55 0 𝜀5
[𝜎6 ] [ 0 0 0 0 0 𝑐66 ] [𝜀6 ]
27
CHAPITRE II / Généralités sur Matériaux Composites
Avec : Vf + Vm = 1 (II.6)
a. Cas général :
Module d'élasticité longitudinal :
E11 = EmVm + EfVf (II.7)
Module d'élasticité transversal :
𝐸𝑚 𝐸𝑓
𝐸22 = 𝐸 (II.8)
𝑓 𝑉𝑚 +𝐸𝑚 𝑉𝑓
Coefficient de poisson :
V12 = 𝜗𝑚 Vm + 𝜗𝑓Vf (II.9)
il faut un schéma :
28
CHAPITRE II / Généralités sur Matériaux Composites
Figure(II.14) : Relations entre système les coordonnées cylindrique et système les coordonnées des
références du composite.
Dans cette partie, les relations entre les composantes de la matrice de rigidité ou de souplesse hors axes et
les composantes exprimées dans les axes principaux de la couche k sont développées dans le cas d’une
rotation autour de l’axe 𝑍̅ = ̅̅̅
𝑁3 .on note θ l’angle de rotation autour de l’axe 𝑍̅ = ̅̅̅
𝑁3 tel que 𝜃 =(𝑥̅ ; ̅̅̅
𝑁1)
.la matrice de rigidité s’exprime alors en fonction des composantes dans la base d’orthotrope sous la forme
suivante :
𝑄̅̂ 𝐾 = 𝑇̂ −1 𝑄̂ 𝐾 𝑇̂ (II.11)
Avec𝑇̂ = 𝑇̂ −𝑇 . En notation vectorielle, le changement de base s’écrit pour les contraintes et les
déformations sous la forme suivante :
𝐶2 𝑆2 −√2𝑆𝐶
−1
𝑇 = [ 𝑆2 𝐶2 √2𝑆𝐶 ]
√2𝑆𝐶 −√2𝑆𝐶 𝐶2 − 𝑆2
suivante :
29
CHAPITRE II / Généralités sur Matériaux Composites
𝑄̅̂ 𝐾 = 𝑇̂ −1 𝑄̂ 𝐾 𝑇̂ (II.13)
Pour simplifier l’écriture, l’indice k est omis dans la suite. Les composantes dans la base
globale s’écrivent en fonction des composantes dans la base d’orthotrope sous la forme
suivante :
̅̅̅̅̅
𝑄11=𝑄11 𝐶 4 +𝑄22 𝑆 4 +2(𝑄12 + 𝑄66 )𝑆 2 𝐶 2
̅̅̅̅̅
𝑄12=(𝐶 4 + 𝑆 4 )𝑄12 + (𝑄11 + 𝑄22 + 2𝑄66 )𝑆 2 𝐶 2
̅̅̅̅̅
𝑄16=(𝑄11 − 𝑄12 − 𝑄66 )√2𝑆𝐶 3 + (𝑄12 − 𝑄12 + 𝑄66 )√2𝐶𝑆 3 (II.14)
̅̅̅̅̅
𝑄22 =𝑄22 𝐶 4 + 𝑄11 𝑆 4 + 2(𝑄12 + 𝑄66 )𝑆 2 𝐶 2
̅̅̅̅̅
𝑄26 =(𝑄11 − 𝑄12 − 𝑄66 )√2𝑆 3 𝐶 + (𝑄11 − 𝑄12 + 𝑄66 )√2𝐶 3 𝑆
̅̅̅̅̅
𝑄66 =(𝐶 4 + 𝑆 4 )𝑄66 + (𝑄11 + 𝑄22 − 2𝑄12 − 𝑄66 )2𝑆 2 𝐶 2
La notation vectorielle retenue permet d’obtenir à partir des relations précédentes les
composantes de la matrice de souplesse 𝑆̂ 𝐾 dans la base globale (ẋ, ẏ, ż) en fonction des
composantes Ŝ k dans la base d’orthotrope (Ṅ1 , Ṅ2 , Ṅ3 ). Pour obtenir l’expression de la relation
entre les composantes de Ŝ k dans la base globale noté S̄ IJ et celles dans la base d’orthotrope noté
SIJ, il suffit de remplacer formellement la lettre Q par la lettre S dans l’expression (II.13).
30
CHAPITRE III
Analyse du réservoir de
Type IV
CHAPITRE III/ Analyse du réservoir de Type IV
CHAPITRE III
Analyse du réservoir de Type IV
Introduction
Ce chapitre s’intéresse au type de réservoir choisi pour cette étude. Il s’agit de présenter
l’architecture et le comportement mécanique des réservoirs de type IV pour la conception de ses
différentes parties.
Figure (III.1) : Schéma des composants d’un réservoir de type IV. [1]
On notera encore que cette enveloppe d’étanchéité sert de mandrin pour le bobinage de
la fibre et doit donc être suffisamment rigide ou pouvoir être rigidifiée, par exemple par une
légère pression interne ou par un autre procédé, pendant la durée de cette opération tout en
conservant une géométrie précise. On notera également que la polymérisation des polymères
utilisés comme résine nécessite en général des températures pouvant atteindre 180°C pour les
résines thermodurcissables, voire plus pour les résines thermoplastiques. Le choix du polymère
de liner ne peut donc s’affranchir de celui du composite.
Les résines les plus utilisées jusqu’à présent sont des thermodurcissables qui exigent
de les polymériser dans des fours.
III.1.1.2 structure résistante.
Les torons de fibres sont déposés par enroulement filamentaire : ce procédé totalement
automatisé permet, couche par couche, de choisir les angles d’enroulement. Ce sont ces
angles qui vont donner à la structure finale son aptitude à résister au chargement de pression
interne. Ainsi dans l’absolu, un angle longitudinal permettrait de bien résister aux efforts
axiaux dus aux effets de fond, tandis qu’un angle transversal assurerait la tenue
circonférentielle.
[Link] externe.
33
CHAPITRE III/ Analyse du réservoir de Type IV
Figure (III.2) : Procédé de fabrication des réservoirs composites par enroulement filamentaire
On distingue trois types d’enroulement pour les réservoirs de stockage Figure (III.3)
l’enroulement « circonférentiel », pour lequel le renfort fait un angle d’environ 90° avec l’axe
de rotation du mandrin, l’enroulement « polaire» utilisé pour le bobinage de sphère et de
réservoirs et enfin l’enroulement « hélicoïdal ou hélicoïdal croisé». Ce dernier permet d’ajuster
l’angle des renforts par rapport à l’axe de la structure en fonction des caractéristiques
mécaniques désirée.
Ce bobinage conduit à une structure composite stratifiée équilibrée du type 90°/±θ par rapport
à l’axe du réservoir. Les couches à 90° sont destinées à reprendre les efforts transversaux induits par
la pression interne sur les parois tubulaires alors que les couches croisées reprennent principalement
les efforts longitudinaux induits par la pression sur les dômes et les embases du réservoir
34
CHAPITRE III/ Analyse du réservoir de Type IV
effectuée en plusieurs points avant le bobinage. La mesure du diamètre du réservoir est
effectuée après la dépose de chaque couche circonférentielles au cours du bobinage. Une série
de mesures est également effectuée après la polymérisation du réservoir afin d’obtenir le
diamètre externe final du réservoir. La mesure de la masse de l’ensemble embase-liner et du
réservoir est également réalisée.
L’ensemble embases-liner est pesé avant le bobinage et le réservoir final est pesé après
polymérisation. La différence des deux mesures donne la masse de composite. La dernière
méthode de contrôle de réservoir « radiographie X » est effectuée après polymérisation. Ces
radiographies permettent de visualiser les différents éléments du réservoir ainsi que leur
épaisseur. Elles permettent également de valider la tenue des interfaces des matériaux traduisant
à la présence de jeu de chaque interface. [33]
III-4 Modèle d’enroulement filamentaire (netting theory) [34]
Il s'agit d'une technique analytique simple utilisée dans la conception de récipients sous
pression composites, en particulier de récipients sous pression composites à enroulement filaments.
Cette théorie d'analyse traite composite comme des stratifiés. Elle établit la relation entre les
contraintes résultant des plis récipient sous pression et la pression interne, les propriétés du matériau
et les paramètres de traitement.
Les hypothèses adoptées pour cette théorie d’ d’enroulement filamentaire sont :
1. Toutes les charges sont supportées uniquement par les fibres.
2. Toutes les fibres sont uniformément sollicitées en tension.
3. La contribution de la matrice à supporter la charge est négligée.
4. La matrice sert uniquement à maintenir la position géométrique des fibres.
5. Toute interaction entre les fibres est également négligée.
6. Les contraintes transversales et les contraintes de cisaillement dans les plis composites
sont faibles, comparées à la résistance ultime à la traction des fibres.
7. La condition de compatibilité des souches est négligée.
8. La contribution au partage de charge du liner est minime ou inexistante.
Considérons la section cylindrique d’un réservoir obtenu par enroulement filamentaire à
extrémité fermée de rayon R, soumis à une pression interne de p (figure III.4). x désigne la direction
longitudinale et y désigne la direction circonférentielle Il est supposé que le récipient sous pression
est par enroulement filamentaire en utilisant uniquement un enroulement hélicoïdal ou une
combinaison d’enroulement circonférentielle et hélicoïdal.
L'enroulement hélicoïdal est réalisé avec un angle α avec la direction axiale x.
35
CHAPITRE III/ Analyse du réservoir de Type IV
𝜎𝑐𝑖𝑟𝑐𝑜𝑛𝑓
𝜎𝑐𝑖𝑟𝑐𝑜𝑛𝑓
𝜎𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑𝑖𝑛𝑎𝑙𝑒
𝜎𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑𝑖𝑛𝑎𝑙𝑒 𝜎𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑𝑖𝑛𝑎𝑙𝑒
Enroulement
hélicoïdal à ± 𝛼 𝜎𝑐𝑖𝑟𝑐𝑜𝑛𝑓
Enroulement
Circonférentiell
e
Figure (III.4) : Les Contraintes dans le réservoir
Supposons que l'épaisseur totale des plis hélicoïdaux soit 𝑡𝛼 . La moitié des plis hélicoïdaux
ont un angle d’orientation +𝛼 L'autre moitié des plis ont un angle – 𝛼 Supposons aussi que l'épaisseur
totale des plis soit 𝑡ℎ .
L'aire des fibres hélicoïdales coupant la longueur unitaire de la surface le long de la direction
circonférentielle est 𝑡𝛼 𝑠𝑖𝑛𝛼 .
L'aire des fibres hélicoïdales coupant la longueur unitaire de la surface le long de la direction axiale
est 𝑡𝛼 𝑐𝑜𝑠𝛼.
Les objectifs de cette analyse seront donc :
1. Si le réservoir est enroulé uniquement avec des fibres hélicoïdales, ±𝛼 avec une contrainte
de fibre admissible, déterminer l'épaisseur du pli hélicoïdale et l'angle d'enroulement α
2. Si le réservoir est enroulé avec à la fois des fibres hélicoïdales ±𝛼 et des fibres
circonférentielles, avec une contrainte de fibre admissible, déterminantes l'épaisseur du pli hélicoïdale
et l'épaisseur du pli circonférentielle..
Tout d'abord, l'équation d'équilibre des forces le long de la direction axiale est considérée.
Circonférentiell 𝑡
e 𝜎𝛼 𝛼cos𝛼
2
1sinα
1cos α
𝑡
𝜎𝛼 𝛼cos𝛼
2
36
CHAPITRE III/ Analyse du réservoir de Type IV
La contribution des contraintes résultantes des fibres circonférentielle dans la direction axiale
peut être négligée.
La composante de la force résultante des fibres hélicoïdal es dans la direction axiale sera
donc𝜎𝛼 𝑡𝛼 𝑐𝑜𝑠 2 𝛼. Où 𝜎𝛼 est la résistance à la traction hélicoïdale
Cela doit être en équilibre avec la force résultante dans la coque composite due à la pression interne
dans la direction axiale.
𝑡𝑓
(𝜎⏟𝑓 1𝐶𝑂𝑆 𝛼 ⏟
𝑐𝑜𝑠𝛼 × ⏟
2 ) ⏟
2𝜋𝑅 =𝑝 ⏟2
𝜋𝑅
⏟ 2 𝑐𝑜𝑚𝑝𝑜𝑠𝑎𝑛𝑡𝑒 𝑝𝑜𝑢𝑟 𝑙𝑒𝑠 𝑝𝑜𝑢𝑟 𝑡𝑜𝑢𝑡𝑒 𝑧𝑜𝑛𝑒 𝑑𝑒 𝑙𝑎 𝑐𝑜𝑞𝑢𝑒
𝑠𝑡𝑟𝑒𝑠𝑠
⏟ 𝑟é𝑔𝑖𝑜𝑛 𝑎𝑥𝑖𝑎𝑙𝑒 𝑓𝑖𝑏𝑟𝑒𝑠±𝛼 𝑙𝑎 𝑐𝑖𝑟𝑐𝑜𝑛𝑓é𝑟𝑒𝑛𝑐𝑒 𝑑′𝑒𝑥𝑡𝑟é𝑚𝑖𝑡é
𝑝𝑜𝑢𝑟 𝑢𝑛𝑖𝑡é 𝑑𝑒 𝑙𝑜𝑛𝑔𝑢𝑒𝑢𝑟
𝑃𝑅 𝑡𝛼
2
= 2𝜎2 2
𝑐𝑜𝑠 2 𝛼 𝑡𝛼 = 2𝜎 𝑃𝑅
𝑐𝑜𝑠2𝛼
(III.1)
2
𝑡 Circonférentiel 𝑡
𝜎𝛼 𝛼cos 𝜎𝛼 𝛼cos
2 2
𝛼 𝛼
S'il n'y a pas de fibres le long de la direction circonférentielle, alors la force résultante dans
la direction circonférentielle du pli du composite due à la pression interne est 𝑝𝑅.
𝑡𝑎
𝑝2𝑅𝐿 = 2𝜎𝑓 𝑠𝑖𝑛𝛼. 2𝐿𝑠𝑖𝑛𝛼 ===> 𝑝𝑅 == (𝜎𝛼 𝑡𝛼 𝑠𝑖𝑛2𝛼 )
2
𝑝𝑅
𝑡𝛼 = (III.2)
𝜎𝛼 𝑠𝑖𝑛2 𝛼
37
CHAPITRE III/ Analyse du réservoir de Type IV
Circonférentiel
le
𝑡
𝑡
𝜎𝛼 𝛼cos𝛼 𝜎𝛼 𝛼cos𝛼
2 2
𝜎𝐻 𝑇𝐻
𝑡𝛼
𝑝𝑅 = (𝜎𝐻 𝑡𝐻 ) + (2𝜎2 𝑠𝑖𝑛2 𝛼) = (𝜎𝐻 𝑡𝐻 ) + (𝜎𝛼 𝑡𝛼 𝑠𝑖𝑛2 𝛼)
2
Où 𝜎𝐻 est la résistance à la traction circonférentielle.
On tenant compte de la valeur de 𝑡𝛼 donnée par la formule (III. 2) On trouve :
𝑝𝑅 𝑡𝑎𝑛2 𝛼
𝑡𝐻 = (1 − ) (III.3)
𝜎𝐻 2
Par conséquent, pour un réservoir sous pression interne cette théorie estime l'épaisseur totale
du pli du composite dans sa section cylindrique à (𝑡𝛼 + 𝑡𝐻 ) .
III.4.1 Conception de la dôme
Le dôme isotensoïde est la forme la plus stable pour les filaments enroulés récipient sous
pression en composite fabriqué par enroulement filamentaire. Il fournit à la structure de poids
minimum et la capacité de charge maximum L’angle d’enroulement entre la trajectoire de la fibre
et le méridien est donné par la condition de Clairaut (figure III.8) [35] :
𝑟0
𝑠𝑖𝑛𝛼 = (III.4)
𝑟
1 𝜌3 𝑑𝜌
𝜉= ∫𝜌 2 2 1/2 (III.5)
𝜌 −𝜌0 6]
[ 2 −𝜌
1−𝜌0
39
CHAPITRE IV
Programmation et
Application
CHAPITRE IV /Programmation et Application
CHAPITRE IV
Programmation et Application
Introduction :
Ce chapitre vise la réalisation de deux objectifs :
IV.1.3 Organigramme
41
CHAPITRE IV /Programmation et Application
Début
Tracé de graphe
𝜉 = 𝑓(𝜌)
Tracé de graphe
𝛼 = 𝑓(𝑟)
Fin
42
CHAPITRE IV /Programmation et Application
Les données de simulation sur ANSYS seront :le maillage de la structure (réservoir), les
caractéristiques mécaniques du matériau, les condition aux limites, la charge (pression interne),
nombre de couche,
Les résultats de la simulation seront les valeurs et les graphes d’iso valeurs des contraintes Von
Mises, les déplacements et les déformations.
IV.4 Application
Dans cette partie du travail, nous sommes appuyés sur les valeurs du type III de réservoirs
d'hydrogène et les avons appliqués au type IV, car les réservoirs sont géométriquement similaires et
pour faciliter notre étude.
43
CHAPITRE IV /Programmation et Application
-
- Figure (IV.2) : Conception et dimensions du réservoir type IV
-
- IV.4.2 Résultats
r
µ
Figure(IV.3) : variation de l’angle d’orientation en fonction du rayon
44
CHAPITRE IV /Programmation et Application
𝝃 = 𝒛/𝒓𝒄
𝝆 = 𝒓𝟎 /𝒓𝒄
- En tenant compte des nombres des couches déterminés les épaisseurs des couches
deviennent :
- L’épaisseur des couches hélicoïdales : thef = npe*ehe = 11 mm
- L’épaisseur des couches circonférentielles : thof = npo*eho = 9.5 mm
- L’épaisseur totale des couches : ttot = thef+thof = 20.5 mm
45
CHAPITRE IV /Programmation et Application
46
CHAPITRE IV /Programmation et Application
Composite (carbone/époxy)
En tenant compte de ces nombres de couches les épaisseurs des couches deviennent :
- L’épaisseur des couches hélicoïdales : thef = 12 mm
- L’épaisseur des couches circonférentielles : thof = 9.5 mm
- L’épaisseur totale des couches : thof = thef+thof= 21.5 mm
La répartition des couches du matériau composite stratifié est indiquée par la figure (VI.7) ci-après.
47
CHAPITRE IV /Programmation et Application
𝟗𝟎𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
±𝟗. 𝟒𝟒𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
∓𝟗. 𝟒𝟒𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
±𝟗. 𝟒𝟒𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
∓𝟗. 𝟒𝟒𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
±𝟗. 𝟒𝟒𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
∓𝟗. 𝟒𝟒𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
±𝟗. 𝟒𝟒𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
∓𝟗. 𝟒𝟒𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
±𝟗. 𝟒𝟒𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟏 𝒍𝒂𝒚𝒆𝒓)
∓𝟗. 𝟒𝟒𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
±𝟗. 𝟒𝟒𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓) 𝟗𝟎𝟎 (𝒄𝒊𝒓𝒄𝒐𝒏𝒇é𝒓𝒆𝒏𝒕𝒊𝒆𝒍𝒍𝒆; 𝟏𝟗 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓)
±𝟗. 𝟒𝟒𝟎 (𝒉𝒐𝒍𝒊𝒄𝒐𝒊𝒅𝒂𝒍; 𝟔 𝒍𝒂𝒚𝒆𝒓)
𝟎
∓𝟗. 𝟒𝟒 (𝟐 𝒍𝒂𝒚𝒆𝒓)
𝟗𝟎𝟎 (𝟐 𝒍𝒂𝒚𝒆𝒓) ±𝟗. 𝟒𝟒𝟎 (𝒉𝒐𝒍𝒊𝒄𝒐𝒊𝒅𝒂𝒍; 𝟔 𝒍𝒂𝒚𝒆𝒓)
IV.4.2.3.6Les graphes :
Le modèle éléments finis sur ANSYS devrait donner les résultats suivants :
Faute de temps alloué à ce travail de recherche cette partie n’a pas été achevée. Nous espérons
qu’elle sera abordée dans des travaux futurs du laboratoire.
48
Conclusion Générale
Conclusion Générale
Conclusion Générale
Cette étude s’est inscrite dans le cadre des travaux de recherche et de développement sur
l’analyse et la conception des réservoirs en composites type IV destinés au stockage de l’hydrogène.
Dans un premier temps, Il agissait de mettre en valeur les connaissances sur l’hydrogène
gazeux, ses propriétés, sa formation, son stockage et ses méthodes de transport. Puis des généralités
des matériaux composites ont été exposées.
Ce travail a permis de développer les méthodes et les théories nécessaires pour la conception
des réservoirs de type VI pour le stockage de l’hydrogène gazeux. Il s’agit de la théorie d’enroulement
filamentaire (netting theory) qui a permis de déterminer l’épaisseur des plis de la stratification et les
angles d’orientation en fonction du rayon de la dôme et dans la partie cylindrique du réservoir. Ceci
a permis aussi de déterminer la forme géodésique de la dôme.
Ces données ont été ensuite utilisées pour modéliser la forme géométrique à l’aide du logiciel
SOLIDWORKS. Puis ce modèle géométrique a été transféré au logiciel ANSYS pour simuler le
comportement mécanique du réservoir sous la pression interne.
Notre visite à une usine (Maghreb Pipe Industries-M’sila) au cours du mois de mai a été d’un
grand intérêt pour découvrir et comprendre les méthodes de fabrication des tubes en matériaux
composites par le procédé d’enroulement filamentaire.
La nouveauté de cet axe de recherche au laboratoire et le temps alloué à ce travail n’ont pas
permis d’exploiter les résultats de la modélisation par élément finis sur ANSYS.
Nous espérons que ce travail ouvre une fenêtre et sera un support pour la continuation de ce
type d’étude sur les réservoirs de stockage de l’hydrogène.
50
Références Bibliographiques
Références Bibliographiques
Références Bibliographiques
52
Références Bibliographiques
[15] Bahlouli, N., « Comportement mécanique en température du pli élémentaire d'un stratifié
carbone/époxy », 1994, Cachan, Ecole normale supérieure.
[16]Jean Marie [16]Berthelot, « matériaux composites, comportement mécanique et analyse des
structures », 4éme édition, TEC &DOC, Le Mans Novembre 2010.
[17] Damien Ducret, « Elasticité anisotrope et endommagement des matériaux composites,
caractérisation ultrasonore et modélisation micromécanique » Thèse de doctorat, institut national
des sciences appliquées, Lyon, 2000.
.[16/06/2022]
nantes.[16/06/2022]
2016 /2017.
[24]Laurent gormet, « Généralités sur les matériaux composites », ECN, Centre Nantes 2008.
[25]A. R. Blythe, « Electrical resistivity measurements of polymer materials », Poly Test, 4 195-
209, 1984.
53
Références Bibliographiques
[35] Gunyoung Park, Hyoseong Jang and Chul Kim Design of composite layer and liner for
Structure safetyofhydrogenpressure vessel (type 4) Journal of Mechanical Science and Technology 35
(8) 2021
[36]: Lei Zu *, Sotiris Koussios, Adriaan Beukers Design of filament–wound domes based on
continuum theory and non-geodesic roving trajectories Composites: Part A 41 (2010) 1312–
1320
[37]Pranjali Sharma , Tapan Bera , Kaladhar Semwal , Rajesh M. Badhe Alok Sharma ,
S.S.V. Ramakumar , Swati Neogi Theoretical analysis of design of filament woundtype 3 composite
cylinder for the storage of compressed hydrogen gas. i n t e r n a t i o n a l j o u r n a l o f h y d r o g
e n e n e r g y 4 5 ( 2 0 2 0 ) 2 5 3 8 6 e2 5 3 9 7
54