LA FUSÉE ARIANE 5

1. MISE EN SITUATION

Suite à l’échec du tout premier vol (501) du lanceur Ariane 5, deux vols de
qualification (les vols 502 et 503) ont été réalisés avec succès les 30 octobre 1997 et
21 octobre 1998. Ceux-ci ont permis de faire quelques modifications et le premier vol
commercial a ensuite eu lieu le 10 décembre 1999. Depuis cette date, de nombreux
lancements ont été réalisés avec un taux de succès qui place cette fusée parmi les
plus compétitives actuellement.

Avec une technologie à corps central à deux étages et deux propulseurs latéraux,
Ariane 5 a marqué une évolution notable dans la famille des lanceurs européens. Ce
choix permet d’obtenir une grande source de puissance au décollage et une
meilleure gestion de la position et de la répartition des masses. En revanche, il induit
des problèmes de couplage importants qui nécessitent une orientation des axes des
tuyères situées en extrémité du corps central et des deux propulseurs latéraux.

1.1. Structure générale de la fusée Ariane 5

La fusée Ariane 5 est globalement composée d’un corps central CC à deux étages et
deux propulseurs latéraux PG et PD (voir figure 1) :
er
• Le 1 étage, commun à tous les types de lancement, est appelé EPC comme
Étage Principal Cryogénique, car il transporte de l'oxygène et de l'hydrogène
liquides refroidis à des très basses températures. Ces fluides permettent
d’alimenter le moteur cryogénique Vulcain situé en extrémité basse de cet
étage. La tuyère de ce moteur est orientable par un montage de type Cardan
selon deux directions orthogonales et elle est refroidie en permanence par la
circulation de l’oxygène et de l’hydrogène liquide.
Le moteur Vulcain 2 est une évolution du moteur Vulcain initial introduite depuis
2002. D’une part, il peut fournir une puissance supérieure d’environ 20 % et,
 d’autre part, son coût de fabrication est inférieur de près de 30 %. Son
introduction a permis d’accroître sensiblement les capacités du lanceur Ariane.
Actuellement, les deux moteurs coexistent jusqu’à la disparition annoncée du
moteur Vulcain initial.
Ce moteur Vulcain est mis en route avant le décollage du lanceur et fonctionne
jusqu’à séparation des deux étages du corps central, soit un peu moins de 10
minutes.
• La configuration du 2ème étage, appelé EPS comme Étage Propulsif
Stockable, peut varier selon la mission. Il contient le satellite ou les éléments à
propulser. Cet étage est proposé en plusieurs configurations modulables de
façon à s’adapter à toutes les demandes des clients. Au début, ce deuxième
étage était appelé à être remplacé par l’avion spatial Hermès, permettant
d’envoyer des hommes dans l’espace : ce projet a été abandonné.
• De chaque côté du corps central CC à deux étages sont fixés les propulseurs
PG et PD ou EAP comme Étages Accélérateurs à Poudre. Les tuyères TG et
TD placées respectivement aux extrémités de ces deux EAP sont aussi
orientables selon deux directions orthogonales par un montage de type Cardan.
La propulsion des EAP, réalisée par la combustion de Propergol solide, est
particulièrement efficace car les deux EAP fournissent 90 % de la poussée
totale du décollage. Allumés peu après le décollage, ils fonctionnent à pleine
puissance pendant une durée de 130 secondes puis, lorsque leurs réservoirs
sont vides, ils se détachent du corps central CC et retombent dans la mer où
leurs carcasses sont ensuite récupérées.

1.2. Le moteur Vulcain

Le moteur Vulcain (voir figures 2 et 3) fonctionne avec un mélange oxygène -

hydrogène liquides (ergols) dont la proportion peut légèrement varier en fonction des
caractéristiques du vol : il fonctionne d’ailleurs toujours avec une légère présence
d’hydrogène dans les gaz brûlés. L’injecteur permet de mélanger et d’enflammer le
mélange réactif à l’entrée de la chambre de combustion. Celle-ci, de forme
convergente puis divergente, est en alliage de cuivre. L’évacuation des gaz brûlés se
poursuit par une tuyère à fort évasement permettant d’accélérer ceux-ci et de créer
la poussée du moteur. L’échauffement induit par cet écoulement nécessite un
refroidissement des parois qui seraient détruites sans cela. À cet effet, une
circulation d’hydrogène liquide super critique est établie autour de la chambre de
combustion et de la tuyère. Une alimentation commune est réalisée par un tore à la
limite chambre – tuyère (voir figure 3) ; le circuit entourant la chambre est constitué
de 360 canaux droits, contenus dans un plan passant par l’axe du moteur, à section
(largeur, profondeur) variable usinés dans le bloc d’alliage de cuivre. Cette évolution
de section est rendue nécessaire par la forte variation de profil de la pièce due à la
présence du col sonique située à l’entrée du divergent de la tuyère. Le circuit
remonte jusqu’à l’injecteur où l’hydrogène participe alors à la combustion. C’est le
circuit régénératif.
La tuyère, située après le col sonique, est refroidie par 456 canaux de section
constante et carrée en Inconel (voir figure 4) dont le soudage bords à bords avec un
mouvement spiralé constitue la tuyère elle-même. Cette technique permet d’avoir un
écoulement d’hydrogène dans des tubes lisses et permet un très fort évasement de
la tuyère. C’est le circuit dump cooling.
L’hydrogène parcourant ce circuit est perdu et ne participe, à sa sortie, qu’à une
poussée minime puisqu’il n’est pas mis en réaction avec de l’oxygène.
 Figure 1

Corps central CC à deux étages

er
Le 1 étage, ou EPC (Étage Principal Cryogénique), est
situé dans la partie basse de ce corps central ; la tuyère
de cet étage est allumée avant le décollage et le reste
jusqu’à séparation des deux étages.
ème
Le 2 étage, ou EPS (Ensemble Propulsif Stockable),
est situé dans la partie haute de ce corps central ; la
tuyère de cet étage est allumée au moment de la
séparation des deux étages.

Moteur Vulcain et sa tuyère centrale TC

La tuyère centrale TC du moteur Vulcain du 1er étage du
corps central CC est située en extrémité basse de ce
→ →
corps. Elle est orientable autour des directions y F et z F

de ± 7°.
→
xF

→ →
yF zF

Propulseurs latéraux PG et PD
Les deux EAP (Étage Accélérateurs à Poudre) latéraux sont utilisés au décollage
pour apporter de la puissance et sont ensuite largués.
Le propulseur EAP gauche ( y 〉 0 ) sera noté PG ; en extrémité de ce propulseur se
→ →
trouve la tuyère gauche TG, orientable autour des directions y F et z F de ± 6° .

Le propulseur EAP droit ( y 〈 0 ) sera noté PD ; en extrémité de ce propulseur se

→ →
trouve la tuyère droite TD, orientable autour des directions y F et z F de ± 6° .
 Figure 2 : Vue en perspective du moteur Vulcain

Turbopompe H2

Turbopompe O 2

Divergent tuyère
 Figure 3 : Schéma descriptif du moteur Vulcain

O2

INJECTEURS

CHAMBRE DE
COMBUSTION

Canaux droits à
section variable

H2 Canaux obliques à
section constante

DIVERGENT

→
xT
Figure 4 : Structure du divergent

a ( x 1)
R 0 = 293,65

→
dx dl x1 y
x2

a( x 2 )

dx

1796

R1 = 880,12

Longueur en millimètres →
xT
2. TRAVAIL DEMANDÉ

L’épreuve proposée s’articule autour de six parties parfaitement indépendantes qui

devront être rédigées sur trois copies différentes, en regroupant les parties 1 et 2
d’une part, les parties 3 et 4 d’autre part et enfin les parties 5 et 6. Ces parties sont :

1. Étude de la phase de décollage de la fusée et validation de la nécessité de

pouvoir orienter les axes des tuyères TG et TD du moteur Vulcain.
2. Étude des modes de vibrations de la tuyère centrale TC, des pulsations
propres, des facteurs d’amortissement et des conditions de résonance.
3. Étude des contraintes dans la tuyère centrale TC par une étude analytique
sur un modèle de type «enveloppe mince».
4. Étude de la déformée de la tuyère centrale TC par une étude analytique sur
un modèle de type « poutre ».
5. Étude thermodynamique des écoulements dans la tuyère.
6. Étude des transferts thermiques à la paroi dans le circuit régénératif et dans
le dump cooling.
 ― Partie 1 ―
Étude de la phase de décollage de la fusée et validation de la nécessité
de pouvoir orienter les axes des tuyères TG et TD du moteur Vulcain

Le but de cette partie est de vérifier qu’il est nécessaire d’avoir un système
d’orientation pour les axes des tuyères TG et TD afin que la trajectoire de la fusée
puisse être pilotée.

La phase de vol étudiée ici correspond à la phase de montée avant largage des
propulseurs latéraux. La fusée est donc complète, à savoir qu’elle est composée :
• de son corps central CC à deux étages ;
• de ses deux propulseurs latéraux à gauche et à droite PG et PD ;
• de la tuyère centrale TC orientable du moteur cryogénique du corps central
CC ;
• des deux tuyères orientables TG et TD des propulseurs latéraux PG et PD.

Pour simplifier l’étude, nous ferons les trois hypothèses suivantes :

Hypothèse 1 : la terre est supposée parfaitement fixe pendant la phase étudiée.
S’il est possible de négliger la rotation de la terre autour du soleil, l’influence de la
rotation de la terre autour de son axe est non négligeable ; cependant, une étude
plus poussée intégrant cette rotation permettrait de montrer que son influence ne
change rien au but de l’étude proposée ici et aboutirait aux mêmes conclusions.
Hypothèse 2 : bien que la masse varie sur la phase de vol étudiée dans cette partie,
nous considérerons que l’ensemble F est à masse conservative.
Cette hypothèse simplificatrice n’est pas restrictive sur la nécessité d’orienter les
axes des tuyères TG et TD.
Hypothèse 3 : les tuyères des propulseurs latéraux sont bloquées en position
centrale et ne sont donc pas orientables.
La fusée Ariane 5 est la première fusée européenne à intégrer les rotations des trois
tuyères. Si l’orientation d’une seule tuyère est théoriquement possible (comme on va
le montrer ici), il est possible que son angle d’inclinaison atteigne les valeurs
maximales admissibles dans les cas extrêmes ; pratiquement, lors des vols, les
inclinaisons n’ont jamais dépassé 2° sur chacune des trois tuyères grâce au travail
combiné de celles-ci.

Paramétrage
→ → →
On associe à la terre le repère R g (O, x g , y g , z g ) où O est le centre de la terre. Ce

→
repère est supposé galiléen. La direction x g correspond à la normale locale

terrestre au point de décollage à Kourou (voir figure 5).

On associe à l’ensemble F = {corps central CC + propulseur gauche PG + propulseur

droit PD + tuyère gauche TG + tuyère droite TD} de la fusée le repère
→ → →
R F (G, x F , y F , z F ) où G est le centre d’inertie de l’ensemble F considéré. La
→
direction x F correspond à l’axe longitudinal de la fusée (voir figure 5). La masse de
→ → →
cet ensemble F est mF et sa matrice d’inertie en G dans la base B F ( x F , y F , z F )

=
A F 0 0 
s’écrit : I (G,F) =  0 BF 0  .
 0 0 CF 

→
zg

Figure 5

→
O yg

Kourou

→
→ xF
xg
G
Le paramétrage de la position de l’ensemble F est réalisé par les six paramètres x, y,
z, ψ, θ et ϕ tels que :
→ → →
• la position relative des origines O et G des deux repères R g (O, x g , y g , z g ) et

→ → →
RF (G, x F , y F , z F ) est réalisée par les trois paramètres linéaires x, y et z tels
→ → → →
que OG = x( t ) x g + y( t ) y g + z( t ) z g .

→ → → → → →
• l’orientation relative des deux bases B g ( x g , y g , z g ) et B F ( x F , y F , z F ) est

repérée par les trois rotations correspondant aux angles d’Euler ψ, θ et ϕ :

→ →
→ →
yg → z1
→ y2
y1 → z2 → yF →
x1 y2 xF
Ψ θ ϕ
→ → →
→ → xg → → y1 → → x
z1 = zg x2 = x1 z2 = zF

On associe à la tuyère centrale TC (la seule qui reste orientable dans notre étude) du
→ → →
corps central CC le repère R T (O F , x T , y T , z T ) où OF est le point tel que
→ →
OF G = x G x F . La masse et les éléments d’inertie de la tuyère centrale TC sont
supposés négligeables face à ceux de l’ensemble F correspondant au reste de la
fusée.
L’orientation de l’axe de la tuyère centrale TC est réalisé par un montage de type
Cardan (voir figure 6) :
Ensemble F → →
y3 = yF
Figure 6

OF 3

Tuyère TC
→
z3

→
xT
Les deux angles δ et ε sont tels que :

→ →
→ xF → y3
x3 → yT →
z3 xT

δ → ε
→
zF → → x3
→ →
y F= y 3 z3 = zT

L’orientation de la tuyère centrale TC est réalisée à l’aide d’un système à vérins. Les
actions mécaniques qui participent à l’orientation de la tuyère centrale TC sont
→ →
modélisées dans cette partie par deux couples purs M(OF ,F → 3) = CMδ y 3 et
→ →
M(OF ,3 → TC) = CMε z 3 .

Le croisillon du Cardan, notée 3 (voir figure 6), et le système à vérins qui provoque la
rotation de la tuyère centrale TC sont supposés respectivement de masse
négligeable. Les deux liaisons pivot qui réalisent ce montage de type Cardan sont
considérées comme parfaites.

Actions mécaniques s’appliquant sur l’ensemble Σ = {F + 3 + TC}

L’ensemble Σ = {F + 3 + TC} est soumis (voir figure 7) :
• aux actions mécaniques de l’environnement, modélisées par le torseur suivant

→   → → → 
R (E → F )  X x + Y y + Z E zF
{T(E → F)}G = →  E F E F
 = → → →
 .
   
M(G,E → F) G L E x F + ME y F + NE z F  G

• à l’action de la pesanteur dirigée selon la verticale locale et assimilée à un

glisseur pur au point G ; la valeur de la norme de l’accélération de la pesanteur
g varie légèrement pendant la montée mais sera considérée comme constante
ici.
• à la poussée sur la tuyère gauche TG du propulseur gauche PG ; cet effort est
→ →
assimilé à un glisseur de résultante F G = FG x F au point A G tel que
→ → →
OF A G = −a x F + b y F .
 • à la poussée sur la tuyère droite TD du propulseur droit PD ; cet effort est
→ →
assimilé à un glisseur de résultante F D = FD x F au point A D tel que
→ → →
OF A D = −a x F − b y F .

• à la poussée sur la tuyère centrale TC du corps central CC ; cette poussée est

→ →
assimilée à un glisseur de résultante F C = FC x T au point O T tel que
→ →
OF O T = −c x T .

→
xF
Figure 7
Corps central CC

→
yF G
→
Propulseur gauche PG xT Propulseur droit PD

Tuyère droite TD
Tuyère gauche TG OF
AG AD
OT
Tuyère centrale TC

Question 1.1. Déterminer l’énergie cinétique de l’ensemble Σ en mouvement par

rapport au repère galiléen R g . Il est demandé de mettre l’expression obtenue sous la

. . . . . . . . .
forme 2T(Σ / R g ) = mF f ( x, y , z) + A F f1(θ, Ψ, θ, ϕ) + BF f2 (θ, Ψ, θ, ϕ) + CF f3 ( Ψ, ϕ, θ) .

Question 1.2.
a. Pour chaque action mécanique agissant sur l’ensemble étudié Σ en
mouvement par rapport au repère galiléen R g , déterminer la puissance

développée.
b. Pour chaque action mécanique agissant à l’intérieur de l’ensemble étudié Σ,
déterminer la puissance développée.
Nota : pour chaque action mécanique, préciser clairement les étapes du calcul et
l’expression obtenue pour le calcul de la puissance. Il n’est pas conseillé de faire la
somme de toutes ces puissances.
Question 1.3. En déduire les huit équations de Lagrange associées aux huit
paramètres x, y, z, Ψ, θ, ϕ, δ et ε qui caractérisent les mouvements de l’ensemble Σ
par rapport à la terre. Il est demandé de ne pas développer les dérivations par
rapport au temps.
Question 1.4. Pour chaque équation trouvée précédemment, préciser l’équivalence
avec les théorèmes généraux de la mécanique (système isolé, théorème utilisé,
éventuel point d’écriture et direction de projection).
Question 1.5. On se place dans le cas particulier d’un vol pour lequel le plan
→ → → →
(G, x F , y F ) reste au voisinage du plan (O, x g , y g ) .

a. Justifier rapidement que dans cette configuration de vol les paramètres

z et θ restent petits.
b. Les angles Ψ et ϕ jouant des rôles similaires, on considère que Ψ est un
paramètre angulaire qui donne l’orientation de la fusée et ϕ une petite
perturbation autour de cette position moyenne.
Vérifier, en considérant les paramètres z, θ et ϕ et leurs dérivées très petits,
que les équations de Lagrange obtenues à la question 1.3. et
correspondant aux paramètres Ψ, θ et ϕ deviennent, après linéarisation
autour des positions z ≈ 0 , θ ≈ 0° et ϕ ≈ 0° :
d
Lψ → [CF (ϕ& + ψ& )] = (FD − FG )b − FC x G (sin ε + θ cos ε sin δ) + ME θ + NE
dt
d
Lθ → [ A F θ& ] = FC x G ϕ cos ε sin δ + L E − ME ϕ
dt
d
Lϕ → [C F (ϕ& + ψ& )] = (FD − FG )b − FC x G sin ε + NE
dt
c. En déduire, à partir de ces équations, que les deux rotations δ et ε sont
nécessaires pour contrecarrer les composantes ME et NE du moment des
actions mécaniques de l’environnement sur la fusée en vol, dans la mesure où
les efforts FC , FD et FG ne peuvent être régulés. Justifier alors la nécessité
d’orienter les axes des tuyères TG et TD.
 ― Partie 2 ―

Étude des modes de vibrations de la tuyère centrale TC, des pulsations

propres, des facteurs d’amortissement et des conditions de résonance

L’objectif de cette partie est de vérifier que les fréquences de résonance de la tuyère
centrale TC du moteur Vulcain restent inférieures à 20 Hz lors de la phase de
décollage sur une hauteur de 2 000 mètres. Cette limitation est choisie pour éviter
les couplages vibratoires avec la structure de la fusée et la commande du système
d’orientation par les deux vérins hydrauliques.

Pour cette étude, le modèle de comportement utilisé pour la tuyère centrale TC est
représenté sur la figure 8. Il s’agit d’un ensemble construit autour d’une liaison pivot
→ →
d’axe (O, z 0 = z ) montée sur un ressort axial et un ressort latéral de raideurs
respectives K A et K L et sur deux amortisseurs de coefficients respectifs λ A et λ L .
La rigidité et l’amortissement de la timonerie liée au bâti 0 sont modélisés
respectivement par un ressort de raideur K B et un amortisseur de coefficient λ B .
Ces trois ressorts et ces trois amortisseurs sont liés au bâti 0.
Bâti 0
λB
Figure 8 B
O →
KB y
λL
→
Ff →
O0 y0
KL
λA
KA Tuyère TC →
n
G →
FC
Bâti 0
β α →
→ y
mg
C
θ
A

→
x0 →
x
Paramétrage
→ → → →
• Le repère (O 0, x 0 , y 0 , z 0 = z ) lié au bâti est supposé galiléen pendant l’étude.
→→→ →
• Le repère (O, x , y , z 0 = z ) est lié à la tuyère centrale TC.

→ → →
• O0O = x( t ) x 0 + y( t ) y 0 .

→ → → →
• ( x 0, x ) = ( y 0, y ) = θ .

→ →
• ( g , x ) = β avec β variable de valeur maximale βMaxi = 10° .
→ →
• ( y , n ) = α avec α variable de valeur moyenne α = 80° .

• AC = L1 = 1,35 m OA = L2 = 1,62 m .

• OB = L3 = 1,10 m OG = L 4 = 0,75 m .

Caractéristiques des ressorts et des amortisseurs

• K A = 58.10 6 N.m −1 λ A = 5.10 5 N.s .m −1 .

• K B = 24.10 4 N.m −1 λ B = 18.10 4 N.s .m −1.

• K L = 55.10 6 N.m −1 λ L = 3.10 5 N.s .m −1 .

Caractéristiques cinétiques de la tuyère centrale TC

La tuyère centrale TC est de masse m = 6 400 kg et la matrice associée à
→ → →
l’opérateur d’inertie de la tuyère exprimée en O dans la base ( x , y , z 0 ) est

=
 A 0 0
 
I (O, TC) =  0 B 0  . On donne : C = 11.10 3 kg.m 2 .
 0 0 C
 

Actions mécaniques s’appliquant sur la tuyère centrale TC

La tuyère centrale TC est soumise (voir figure 8) :
 • à la résultante des forces de pression appliquées sur le flanc de la tuyère ; cet
→ →
effort est assimilé à un glisseur de résultante F f = −Ff y au point O avec

Ff = 25.105 N .
• à l’action de deux vérins d’orientation de la tuyère ; cet effort est assimilé à un
→ →
glisseur de résultante F C = FC n au point C avec FC = 1.10 5 N (valeur

moyenne). Pour l’étude on prendra α = 80° .

→
• à l’action de la pesanteur qui est inclinée par rapport l’axe (O, x ) d’un angle
β ; en phase de décollage vertical β = 0° et, pour l’étude, on prendra
β = βMaxi = 10° , ce qui correspond à la valeur à la fin de la montée avant

séparation des deux étages du corps central. L’accélération de la pesanteur g

est considérée comme constante et vaut : g = 9,81 m.s −2 .

L’action de la timonerie sur la tuyère centrale TC sera modélisée par le glisseur :
→ → . . → 
   R ( timonerie → TC ) = −K B ( y − L 3 sin θ ) y 0 − λ B ( y − L 3 θ cos θ ) y 0
   
T( timonerie → TC) = 
   → 
 B  0 
 B

Question 2.1. Déterminer les trois équations différentielles, sous forme littérale,
faisant intervenir les paramètres x, y et θ , obtenues par projection dans la base
→ → →
( x 0 , y 0 , z ) du Principe Fondamental de la Dynamique appliqué en O au modèle de

la tuyère centrale TC proposé sur la figure 8.

Question 2.2. En considérant que les variations des paramètres x, y et θ restent très
petites en amplitude et en vitesse, linéariser alors ces trois équations différentielles.
Écrire ces trois équations linéarisées sous la forme : [M](&q&) + [C](q& ) + [K ](q) = [A ],

 &x&   x&  x

     
avec (&q&) =  &y&  , (q& ) =  y&  , (q) =  y  , [M] matrice de masse, [C] matrice
 &θ&   θ&   θ
     
d’amortissement et [K ] matrice de raideur. L’application numérique, non demandée,
donne :
  ..  .
64.102 0 0   x  5.105 0 0   x  58.106 0 − 249.10 4   x   − 366.102 
  ..   .   
 0 64.102 48.102   y  +  0 48.104 − 2.105   y  +  0 55.106 − 227.103   y  =  − 249.104 
 0      
 48.102 11.103   ..   0 − 2.105 22.10 4   .   0 26.104 29.10 4   θ   153.103 
  θ    θ    
   
Question 2.3. Calculer les pulsations propres non amorties et les modes propres
associés.
Question 2.4. Soit [P] la matrice des vecteurs propres exprimée dans la base
→ → →
( x 0 , y 0 , z ) . En déduire sous forme numérique les matrices [M * ]= [P]T [M] [P] ,

[C *]= [P]T [C] [P] et [K *]= [P]T [K ] [P] .

Question 2.5. Les éléments des matrices [M *] , [C *] et [K *] sont notés

respectivement mij* , c ij* et k ij* . On supposera en première approximation que les

c ii*
facteurs d’amortissement s’écrivent : ξ ii = . Calculer ces facteurs
2 k ii* mii*

d’amortissement.
Question 2.6. Conclure quant au respect du critère relatif aux fréquences de
résonance imposé par le cahier des charges.
 ― Partie 3 ―

Étude des contraintes dans la tuyère centrale TC par une étude analytique sur
un modèle de type « enveloppe mince »

L’objectif de cette partie est de déterminer les contraintes dans la partie divergente
de la tuyère centrale TC sous l’action de la pression intérieure de l’écoulement.

Pour l’étude proposée le divergent de la tuyère centrale TC est modélisé comme un

cône de révolution d’épaisseur e constante (voir figure 9 ci-dessous).
Le matériau constitutif de la tuyère est supposé élastique linéaire, homogène et
isotrope dans cette partie.

Figure 9
→
er

α
s

rmini
e →
en

→
p(z ) = p0 = cste es
Tuyère TC

rMaxi

→
ez
La forme réelle du divergent n’est pas conique mais la modélisation retenue donne
une très bonne approximation pour l’étude qui va suivre.
L’écoulement est stabilisé. Pour simplifier les calculs, on supposera que :
• la pression à l’intérieur de la tuyère est constante et égale à p 0 ,

• la pression à l’extérieur de la tuyère est nulle,

• les effets dus à la pesanteur sont négligeables devant ceux dus à la pression
à l’intérieur de la tuyère.
Le problème étant de révolution, on va utiliser un paramétrage cylindrique (r, θ, z) .
→ → → →
On rappelle que pour un vecteur u = ur (r, θ, z ) e r + u θ (r, θ, z) e θ + u z (r, θ, z) e z , le

 ∂ur 1 ∂ur u θ ∂ur 

 ∂r −
r ∂θ r ∂z 
= →  ∂u 1 ∂u θ ur ∂u θ 
gradient s’écrit : g rad u =  θ +  .
 ∂r r ∂θ r ∂z 
 ∂u z 1 ∂u z ∂u z 
 ∂r r ∂θ ∂z  ( → → →
 e r , e θ, e z )

Question 3.1. On décide de prendre comme vecteur déplacement infinitésimal d’un

→ → →
point M de la structure déformée le vecteur u (M) = ur (r, z) e r + u z (r, z) e z . Justifier
rapidement la validité de ce choix.
Question 3.2. Justifier que la matrice associée au tenseur des contraintes s’écrit
 σ rr
=
0 σ rz 
sous la forme : σ =  0 σ θθ 0  . Les expressions des différents
σ rz 0 σ zz  ( e , e , e )
→ → →
r θ z

éléments de cette matrice ne sont pas demandées.

Pour la suite de l’étude, comme l’épaisseur du divergent de la tuyère centrale TC est
très petite devant les autres dimensions, le modèle d’enveloppe mince peut être
retenu : dans ce cas, nous ferons l’hypothèse que la matrice des contraintes dans le
→ → →
milieu continu peut être exprimée, dans la base locale (en , e θ , e s ) , par :

=
σnn (s) 0 0 

σ= 0 σ θθ (s) 0 
 0 0 σ ss (s) ( → → →
e n, e θ, e s )
avec les contraintes σ nn (s) , σ θθ (s) et σ ss (s) ne dépendant que de l’abscisse

curviligne s (voir figure 9) et donc supposées uniformes dans une tranche et selon
l’épaisseur.
Question 3.3. La contrainte σ nn (s) est approchée par la moyenne des contraintes
σ nn (s) sur les parois intérieure et extérieure du divergent de la tuyère. En appliquant
cette hypothèse, déterminer la valeur de σ nn (s) .
Question 3.4. Pour déterminer l’expression de la contrainte σ ss (s) , on isole un

volume élémentaire de la tuyère d’épaisseur e, de hauteur ds et de largeur rdθ

autour d’un point M. Faire le bilan des actions mécaniques agissant sur cet élément
→
et écrire le théorème de la résultante statique en projection sur la direction e z . On
fera un schéma clair de l’élément isolé sur lequel on précisera les actions
σ ss (s) dσ ss (s) p 0
mécaniques extérieures à cet élément. En déduire que + = tan α .
s ds e
Sachant que l’extrémité basse du divergent est libre, déterminer l’expression de
σ ss (s) en fonction de p 0 , e, α , r(s) et rMaxi .

Question 3.5. Pour déterminer l’expression de la contrainte σ θθ (s) , on isole une

« demi-tranche » de la tuyère d’épaisseur e, de hauteur ds et de longueur r π . On

fera un schéma clair de l’élément isolé sur lequel on précisera les actions
mécaniques extérieures à cet élément. Appliquer le théorème de la résultante
p0
statique à cette «demi-tranche» et montrer que σ θθ (s) = r ( s) .
e cos α
Question 3.6. Vérifier que la contrainte σ nn (s) reste très inférieure, en valeur
absolue, aux deux autres contraintes σ θθ (s) et σ ss (s) . Pour la suite de l’étude, la

contrainte σ nn (s) sera alors négligée.

Question 3.7. Exprimer la contrainte équivalente de Von Mises en fonction de
σ θθ (s) et σ ss (s) puis en fonction de p 0 , e, α , r(s) et rMaxi .

En déduire la zone la plus contrainte et donner l’expression de la contrainte

équivalente de Von Mises dans cette zone en fonction de p 0 , e, α , rmin i et rMaxi .
Question 3.8. Les caractéristiques géométriques du divergent de la tuyère centrale
TC sont :
• rmin i = 0,29365 m ;

• rMaxi = 0,88012 m ;

• zMaxi − z min i = 1,796 m ;

• e = 30 mm .

La pression maximale à l’intérieur de la tuyère est : p 0 = 110.10 5 Pa .

Vérifier que la limite élastique σ 0,2 = 600.10 6 Pa n’est pas atteinte.
 ― Partie 4 ―

Étude de la déformée de la tuyère centrale TC par une étude analytique sur un

modèle de type « poutre »

L’objectif de cette partie est de vérifier que, lors de la phase de décollage sur une
hauteur de 2 000 mètres, la déformation de la tuyère centrale TC ne perturbe pas
l’asservissement d’orientation de celle-ci.
Pour cela, on va utiliser un modèle simplifié pour lequel les hypothèses de la
Résistance des Matériaux sont supposées vérifiées.

L’asservissement en orientation est conçu pour que la tuyère centrale TC ne s’incline

→
pas d’un angle supérieur à ± 3 ° autour de l’axe (O, z ) , défini sur la figure 10 ci-
dessous, lors de la phase de décollage sur 2 000 mètres. Pour ne pas perturber
l’asservissement, le cahier des charges impose que la déformée angulaire au

voisinage du point O reste inférieure à ± 0,5 ° .

Pour faire cette étude, la tuyère centrale TC est modélisée comme une poutre
cylindrique, de section constante, articulée en O et en appui au point B avec la fusée
considérée comme le bâti (voir figure 10).

B
Figure 10
→
y

O →
Ff

Tuyère TC
G →
n
→
FC
→→→
Le repère (O, x , y , z ) est lié
α →
à la tuyère centrale TC. β y
→
mg C
A
→
x
Sur cette figure, la tuyère centrale TC est représentée dans une position d’orientation
quelconque par rapport à la fusée.

Paramétrage
→→
• ( g , x ) = β avec β variable de valeur maximale βMaxi = 10° .
→ →
• ( y , n ) = α avec α variable de valeur moyenne α = 80° .
• AC = L1 = 1,35 m OA = L 2 = 1,62 m .

• OB = L3 = 1,10 m OG = L 4 = 0,75 m .

Caractéristiques de la tuyère centrale TC

• m = 6 400 kg .

• module d’Young : E = 20.1010 Pa .

→ →
• moment quadratique par rapport à l’axe (G, z ) : I (G, z ) = 1.10 − 2 m 4 .

Actions mécaniques s’appliquant sur la tuyère centrale TC

Elle est soumise (voir figure 10) :
• à la résultante des forces de pression appliquées sur le flanc de la tuyère ; cet
→ →
effort est assimilé à un glisseur de résultante F f = −Ff y au point O avec
Ff = 25.105 N .
• à l’action des vérins d’orientation de la tuyère ; cette action est assimilée à un
→ →
glisseur de résultante F C = FC n au point C avec FC = 1.10 5 N . Pour l’étude,

on prendra α = 80° .
→
• à l’action de la pesanteur qui est inclinée par rapport l’axe (O, x ) d’un angle
β ; en phase de décollage vertical β = 0° , et pour l’étude, on prendra
β = βMax = 10° , ce qui correspond à la valeur à la fin de la montée avant

séparation des deux étages du corps central. L’accélération de la pesanteur g

est considérée comme constante et vaut g = 9,81 m.s −2 .

Question 4.1. En utilisant la méthode de votre choix, déterminer littéralement la
déformée angulaire au voisinage du point O.
Question 4.2. Faire l’application numérique et conclure quant au respect du critère
relatif à la déformée angulaire au voisinage du point O imposé par le cahier des
charges.
 ― Partie 5 ―

Étude thermodynamique des écoulements dans la tuyère

L’objectif de cette partie est de caractériser l’écoulement isentropique des gaz brûlés
dans le moteur Vulcain et particulièrement la température de ceux-ci et de
déterminer la poussée de ce moteur.

Il est conseillé de se reporter à l’annexe pour les définitions générales de

thermodynamique classique.
Il est demandé de respecter scrupuleusement dans les développements les
notations fournies dans cette annexe et dans le texte ainsi que les valeurs
numériques fournies.

Question 5.1. On considère un écoulement permanent de gaz parfait dans un tube

de section variable d’aire A tel que représenté sur la figure 11. La vitesse considérée,
dans une section de pression P uniforme, est la vitesse moyenne w de l’écoulement.
Avec ces notations exprimer le travail Wp des forces de pression en fonction de P1 ,

v1, P2 , v 2 et de la masse dm contenue dans le volume A 1B1C1D1 lors d’un petit

déplacement A1C1 du fluide.

Figure11

D1
T1
B1
w1 B2
P1 T2
D2
C1
A1
P2
A2
C2 w2
z1

z2
Question 5.2. En faisant les hypothèses suivantes :
P
• le gaz est parfait ( = rT ), non visqueux,
ρ
• l’écoulement se fait sans échange de chaleur avec les parois,
• les forces de pesanteur sont négligeables, pour le gaz,
appliquer le principe général de la conservation de l’énergie lors de ce déplacement
et montrer que la relation différentielle liant la vitesse w de l’écoulement dans une
section à pression constante et l’enthalpie H du gaz dans cette même section s’écrit
w dw + dH = 0 .

Ces hypothèses seront conservées pour les questions 5.3. à 5.10..

Question 5.3. En s’aidant de la définition de la vitesse du son et du nombre de Mach

M de l’écoulement dans une section, montrer que l’on peut relier la vitesse w à la
dw 1 dP
pression P dans une section par : =− .
w γ M2 P
Question 5.4. En exprimant la conservation du débit masse, écrire la relation liant
dA dw dv dA dw
, et et en déduire la relation : = (M2 − 1) .
A w v A w
dA dP
Question 5.5. Déduire des relations précédentes la relation liant à et
A P
caractériser graphiquement, en fonction de l’abscisse x suivant l’axe de la tuyère, les
variations de A, P, M et w lorsque l’on veut obtenir une vitesse supersonique
d’écoulement des gaz brûlés en sortie de tuyère.
Question 5.6. On s’intéresse à la caractérisation des grandeurs de l’écoulement
dans une section de la tuyère en fonction des caractéristiques de l’écoulement dans
la section d’entrée notée 1 ( P1 , T1 , w1 , γ , r). Exprimer la vitesse locale w à l’aide de
ces paramètres et de la pression locale P en se servant de la forme intégrée de la
relation établie à la question 5.2. et de l’expression de la variation de l’enthalpie dans
une évolution adiabatique.
Question 5.7. Lorsque la pression de sortie tend vers zéro (ce qui est le cas lors du
vol d’Ariane 5), montrer que la vitesse maximale w s atteinte par les gaz brûlés dans
1/ 2
 2γ 
la section de sortie de la tuyère s’exprime par : w s =  w 12 + r T1  .
 γ −1 
Question 5.8. En reprenant la valeur numérique du débit masse des ergols introduits
dans le moteur, calculer la force de poussée du moteur Vulcain dans le vide.
Question 5.9. Exprimer la température moyenne Tc des gaz brûlés au col en
fonction des caractéristiques de l’écoulement en entrée, puis faire l’application
numérique. On repartira éventuellement de l’expression obtenue à la question 5.6.
en faisant intervenir la vitesse sonique de l’écoulement des gaz.
Question 5.10. Montrer que dans la section de sortie de la tuyère, la température
des gaz peut s’exprimer, lorsque l’on néglige le nombre de Mach en entrée, par la
2T1
relation : Ts = .
2 + M2s ( γ − 1)
 ― Partie 6 ―

Étude des transferts thermiques à la paroi dans le circuit régénératif et dans le

dump cooling

L’objectif de cette partie est de définir des canaux de refroidissement capables de

maintenir la température interne de la chambre et de la tuyère à des valeurs
compatibles avec la tenue mécanique de ces pièces.

Il est conseillé de se reporter à l’annexe pour les définitions générales de

thermodynamique classique.
Il est demandé de respecter scrupuleusement dans les développements les
notations fournies dans cette annexe et dans le texte ainsi que les valeurs
numériques fournies.

En réalité l’écoulement des gaz n’est pas tout à fait adiabatique car les parois de la
chambre et de la tuyère sont refroidies. Nous allons nous intéresser maintenant aux
transferts thermiques à la paroi. La partie 5 donne une condition limite en
température et vitesse du coté des gaz brûlés et permet alors d’obtenir à partir de
considérations de mécanique des fluides la valeur du coefficient de transfert
thermique hg entre la paroi et ces gaz. Le coefficient de transfert entre la paroi et

l’hydrogène est également issu de corrélations appropriées et est supposé connu. Il

faut maintenant définir des canaux de refroidissement capables de maintenir la
température interne de la paroi de la chambre et de la tuyère à des valeurs
compatibles avec la tenue mécanique de ces pièces ( T 〈 750 K ). On va pour cela
considérer une tranche axiale de tuyère proche du tore d’injection et contenant le
circuit régénératif. Les figures 12 et 13 montrent un détail du cheminement thermique
entre les gaz brûlés et l’hydrogène liquide parcourant les canalisations.
 Figure 12 : Cheminement du flux entre gaz brûlés et hydrogène

a c
L
Th

e
Tg Lignes de flux
Ka

Figure 13

1/2 Ailette

1/2 Canal

Flux thermiques

Distinction des flux d'ailette

et de fond de canal
Phénomène de pied d'ailette

Le transfert conductif est régi par la seule loi de Fourier, reliant le flux au gradient de
température, définie en annexe.
La partie centrale séparant 2 canaux est assimilable à une ailette que nous allons
étudier en premier.

Question 6.1. Considérons une ailette de température moyenne à la base T0 (voir

figure 14) et baignée par un fluide de température uniforme Tf . Localement elle

échange avec celui-ci un flux de chaleur dφ sur une surface dS caractérisée par la
température locale de la surface externe de l’ailette T et un coefficient d’échange

constant h d’unité W.m −2 .K −1 tel que dφ = h dS(T − Tf ) .

Figure 14 : Schéma et notations concernant une ailette

Température T
Section A

→
x
T0 x x+dx

Périmètre P
Surface dS

Le nombre de Biot B i , caractérisant le rapport entre convection externe et conduction

Ah
interne, est défini ici, avec les notations de la figure 14, par Bi = .
Pλ
On ne peut considérer la température dans une section d’ailette comme uniforme
que si le nombre de Biot B i est inférieur ou égal à 0,1.
Compte tenu des valeurs numériques communiquées en annexe, peut-on considérer
qu’une paroi inter canaux (jouant le rôle d’ailette dans le circuit régénératif) a une
température transversale uniforme ? (On se placera au niveau du tore d’introduction
de l’hydrogène liquide).
Question 6.2. Dans l’affirmative, écrire l’équation de la conservation du flux entre les
sections x et x + dx de l’ailette définie figure 14 en respectant ses notations.

Les caractéristiques thermophysiques du matériau seront considérées comme

constantes spatialement et indépendantes de la température pour cette
question et les suivantes.

Question 6.3. Intégrer l’équation différentielle du second ordre obtenue à la question

6.2. et montrer qu’elle peut se mettre sous la forme : θ = C1e mx + C 2 e −mx .

On précisera ce que représentent θ , C1 , C 2 et m.
Question 6.4. En supposant nul le flux sortant par l’extrémité de l’ailette, montrer
que le flux φ évacué latéralement par l’ailette de périmètre P, de section A, de
longueur L, de conductivité thermique λ, de coefficient d’échange avec le fluide h,
s’écrit φ = (T0 − Tf ) h p λ A tanh (mL ) avec tanh : tangente hyperbolique.
Question 6.5. On décompose le flux thermique passant des gaz brûlés vers
l’hydrogène en deux chemins (voir figures 12 et 13). Le premier passe à travers la
paroi et rejoint l’hydrogène par le fond du canal. Le second alimente l’ailette
constituée par la paroi séparant deux canaux et atteint l’hydrogène par le coté latéral
du canal (on montrera plus tard que le flux s’écoulant en extrémité de paroi inter
canaux est négligeable). La partition de ces deux flux est une ligne de flux
(assimilable à une ligne droite) s’accrochant au coin inférieur du canal et partageant
la base du canal dans un rapport K. Un examen plus approfondi des conditions
thermiques conclut à l’égalité des températures de surface de la paroi coté gaz
brûlés.
a. Exprimer, sous la forme φ c = f (K ) (Tg − Th ) , le flux (par canal et par unité de

longueur) transféré par la partie de paroi située sous le canal en fonction de

l’écart Tg − Th des températures de gaz brûlés et d’hydrogène, des coefficients

d’échange thermique h g et hh du coté des gaz et de l’hydrogène, de la

conductivité λ du cuivre, de la partition K du flux et de la géométrie du canal

définie figure 12. II est conseillé de s’aider de la formulation du flux à travers
une section trapézoïdale donnée en annexe.
b. Exprimer de la même façon, sous la forme φ a = g(K ) (Tg − Th ) , le flux

transféré des gaz brûlés vers l’hydrogène par la paroi de séparation entre
canaux ( par unité de longueur de canal).
Question 6.6. Exprimer l’équation en K avec le symbolisme f(K) et g(K) assurant
l’uniformité de température de paroi coté gaz brûlés sous le canal et sous la paroi
séparant deux canaux.
Question 6.7. La condition résultant de la question 6.6. est vérifiée pour K = 0,51.
Évaluer alors φ a et φ c et vérifier la validité de cette valeur de K. Évaluer ensuite la

température de paroi obtenue coté gaz brûlés et le flux transféré à la paroi de la

tuyère sur toute sa circonférence par mètre de longueur de canal lorsque l’on se
situe en entrée du circuit régénératif au niveau du tore d’introduction de l’hydrogène
liquide.
Question 6.8. À la même cote axiale du circuit régénératif, calculer la température
de la paroi inter canaux en son extrémité supérieure. Peut-on en conclure que le flux
éventuellement perdu par l’extrémité de cette paroi est négligeable ?
Question 6.9.
a. Exprimer littéralement la température d’hydrogène au bout d’un parcours de
longueur l dans les canaux ?
b. Déterminer la température d’hydrogène en entrée d’injecteur du moteur
Vulcain sachant que le cumul du flux échangé à la paroi sur toute la longueur
du circuit régénératif est de 20 MW .
 ― Annexe ―

1. RAPPELS DE THERMODYNAMIQUE ET NOTATIONS

• P : pression d’un fluide en Pa .

• v : volume massique d’un fluide en m 3 .kg −1 .
• T : température d’un fluide en K.
• Cp : chaleur spécifique à pression constante d’un fluide en [Link] −1.K −1 .
• C v : chaleur spécifique à volume constant d’un fluide en [Link] −1.K −1 .
• r = Cp − C v : constante massique d’un fluide en [Link] −1.K −1 .
Cp
• γ= : coefficient d’évolution adiabatique d’un fluide.
Cv
• a s = ( γ P v )1 / 2 : célérité du son dans un fluide en m.s −1 .
• w = vitesse d’écoulement d’un fluide en m.s −1 .
w
• M= : nombre de Mach de l’écoulement.
as
• Pv = rT : équation d’état des gaz parfaits.
• Pv γ = cste : équation caractérisant l’évolution adiabatique réversible d’un fluide.
• U : énergie interne massique d’un fluide en [Link] −1 .
• H = U + P v : enthalpie massique d’un fluide en [Link] −1 .
w2
• Ht = H + : enthalpie totale massique d’un fluide en écoulement en [Link] −1 .
2
• Q : chaleur échangée avec l’extérieur lors d’une évolution d’un fluide en [Link] −1 .
• W : travail échangé avec l’extérieur lors d’une évolution d’un fluide en [Link] −1 .
• dU = dW + dQ : relation caractérisant la conservation de l’énergie d’une évolution
sans écoulement échangeant dQ et dW avec l’extérieur.

2. RAPPELS DE THERMIQUE ET NOTATIONS

• λ : conductivité thermique d’un milieu solide ou fluide en W.m −1.K −1 .

• h : coefficient d’échange entre la paroi d’un solide et son milieu environnant en

W.m −2 .K −1 .
 dT
• φ = −λ A , loi de Fourier traduisant que le flux thermique s’écoulant dans un
dx
→
solide à travers une surface A perpendiculaire à l’axe des x est :
- proportionnel à la conductivité thermique du solide ;
- proportionnel à la surface considérée ;
- proportionnel au gradient de température local ;
- de même sens que celui des isothermes décroissantes.
• Flux thermique passant à travers une zone trapézoïdale délimitée par ses bases de
λ L (a − b )
températures Ta et Tb (avec les notations de la figure 15) : φ = (Ta − Tb ) .
a
e ln
b

Figure 15
Base de température Ta

e Base de température Tb

Rappel mathématique :
e x − e−x
• tanh ( x ) = avec tanh : tangente hyperbolique.
ex + e− x
3. DONNÉES NUMÉRIQUES CONCERNANT LE MOTEUR VULCAIN
HM-60

Circuit régénératif Circuit dump cooling

Extrémité Tore
Col Tore Sortie
Repère coté d’injection
sonique d’injection tuyère
chambre entrée H2

Développement axial (m) 0 0,426 0,700 0,72 2,53

Rayon intérieur R i (m) 0,207 0,131 0,293 0,300 0,88

Épaisseur paroi e (mm) 1,5 1 2,5 0,4 0,4

Largeur canal a (mm) 2 1,3 2,6 3,2 3,2

Hauteur canal L (mm) 11 10,4 12 3,2 3,2

Épaisseur paroi entre

1,61 1 2 0,4 0,4
canaux c (mm)
Température des gaz
3 600 3 317 1 941 1 930 1 300
brûlés (K)

Température du H2 (K) 35 35

Coefficient d’échange
paroi - gaz brûlés hg 17 300 31 600 7 550

( W.m −2 .K −1 )
Coefficient d’échange
paroi canal - Hydrogène 38 000
hh ( W.m −2 .K −1 )

Conductivité de la paroi
380 380 380 20 20
métallique ( W.m −2 .K −1 )

Pression d’alimentation
130.105
(Pa)
• Débit H2 dans circuit régénératif (alimentant l’injecteur) : 36 kg.s −1.

• Débit H2 dans circuit dump cooling (Hydrogène perdu) : 2,2 kg.s −1.

• Débit O 2 dans injecteur : 224 kg.s −1.

• rgazbrûlés = 565,5 [Link] −1.K −1 .

Cp
• γ= = 1,23 (gaz brûlés).
Cv

• Cp moyen hydrogène = 12 800 [Link] −1.K −1 .