Dimensionnement d'un SREC Mécanique
Dimensionnement d'un SREC Mécanique
Bureaux d’Études
Dimensionnement
Question 8.3 • Modifier le modèle de façon à prendre en compte une anticipation possible du freinage par lecture
des panneaux 100 m environ avant que le véhicule ne passe à leurs niveaux.
Question 1 • Déterminer l’énergie à dissiper pour freiner jusqu’à l’arrêt un véhicule de masse Mv = 1t alors qu’il est à
une vitesse de Vv = 110km/h. On supposera que 90% de l’énergie est dissipée dans le système de freinage, le reste
l’étant par pertes de roulement et de frottement.
1
Ec,v = MvVv2 = 4.67105 J = 1, 63kW h
2
Energie à dissiper par le système de freinage :
Il existe plusieurs Systèmes de Récupération d’Energie Cinétique (SREC ou KERS en anglais). On peut citer deux
grandes familles associées au type d’énergie stockée :
– le SREC mécanique : l’énergie cinétique du véhicule récupérée permet la mise en rotation d’un volant d’inertie ;
– le SREC électrique : l’énergie cinétique du véhicule récupérée est transformée en énergie électrique au moyen
d’un alternateur, puis stockée dans une batterie ;
Le SREC mécanique possède un rendement assez élevé, de l’ordre de 80%, alors que le SREC électrique, du
fait des multiples transformations d’énergie (mécanique - électrique - chimique, puis chimique - électrique - mécanique
lors de la restitution) est de l’ordre de 35%. On propose donc ici de retenir la solution du SREC mécanique.
Question 3 • Les contraintes mécaniques permettant le dimensionnement du SREC sont multiples. Mais un élément
limitant essentiel est la vitesse de rotation. Précisez quelles peuvent être les causes de dysfonctionnement liés à la
vitesse.
On ne peut donc pas emmagasiner une énergie supérieure à une valeur maximale. Cela veut dire qu’il faut conser-
ver un système de freinage classique en parallèle.
Pour des raisons de place disponible dans le véhicule on se limite à un système de 30cm de diamètre (±10%) et
de 20cm de long (±10%). Pour des raison de consommation de carburant, on souhaite réaliser un système de moins
de 100kg afin de ne pas trop augmenter la masse totale du véhicule.
Pour simplifier l’étude, on assimile le volant à un anneau circulaire de section carrée de rayon intérieur r1 , rayon
extérieur r2 et d’épaisseur e. Cet anneau est supposé libre d’effort sur ces bords et soumis à un mouvement de rotation
autour de son axe à une vitesse ω. On commencera le dimensionnement avec les valeurs suivantes : r1 = 10mm,
r2 = 150mm, e = 100mm et en considérant un disque en acier : E = 210MPa, ν = 0.3 et ρ = 7800kg/m3 .
Question 4 • Déterminer les caractéristiques inertielles du volant.
Question 5 • Déterminer la vitesse de rotation d’un tel volant afin qu’il puisse stocker l’énergie voulue.
2 Dimensionnement Mécanique
On cherche dans un premier temps à quantifier les risques d’éclatement du volant. Pour cela on cherche à
déterminer l’état de contrainte dans le disque en rotation. On ne s’intéresse pas à la liaison entre le volant et son
axe.
Question 6 • Poser le problème d’élasticité correspondant au disque seul en rotation autour de son axe à une vi-
tesse ω. Rappeler les différentes manière de résoudre un tel problème.
Question 7 • On fait une hypothèse de contraintes planes dans le plan perpendiculaire à l’axe du cylindre. Représenter
la forme du tenseur des contraintes dans un système d’axe cylindrique. Quelles est la validité de cette hypothèse ?
σrr σrθ 0
σ = σrθ σθθ 0
0 0 0 (e ,e
r θ ,ez )
Question 8 • On ne fait donc qu’une résolution du problème dans le plan (er , eθ ). On suppose que le champ de
déplacement plan est de la forme :
U(r, θ) = u(r)er
Justifier cette hypothèse.
u
εrr = u, r ; εθθ = ; εθθ = 0
r
Question 10 • Exprimer la relation de comportement (contraintes fonction des déformations) pour cette situation. E et
ν sont les caractéristiques élastiques du matériau.
Question 11 • Traduire l’équilibre en volume du volant sous forme de condition sur la composante de déplacement.
Question 13 • En appliquant les conditions aux limites sur les faces interne et externe de l’anneau, déterminer
complétement la solution en déplacement.
ω2 r12
µ
(1 + ν)λ − (1 − ν) − (3 + ν)ρ =0
σrr (r = r1 ) = 0 r12 8
⇒
σrr (r = r2 ) = 0 µ ω2 r22
(1 + ν)λ − (1 − ν) 2 − (3 + ν)ρ =0
r2 8
Donne :
(3 + ν) ω2 2
λ= ρ (r + r22 )
(1 + ν) 8 1
(3 + ν) ω2 2 2
µ=
ρ r r
(1 − ν) 8 1 2
Donc :
E(3 + ν) ω2 1 2 2 1 r12 r22
− r2
u(r) = 2
ρ (r1 + r2 ) − 2
1−ν 8 1+ν 1−ν r
r
1
σV M = (σI − σII )2 + (σI − σIII )2 + (σII − σIII )2
2
En contraintes planes : r
1 2
σV M = σrr + σ2θθ + (σrr − σθθ )2
2
Question 17 • Quelle est la zone la plus chargée. On propose différents types matériaux utilisables dans le tableau 5.
Un disque en acier est-il envisageable ?
3 Dimensionnement en Dynamique
On cherche maintenant une approximation de la fréquence fondamentale du système à l’aide de différents modèles
dynamique.
Question 19 • Rappeler pourquoi, en première approximation, l’obtention de la pulsation fondamentale est suffisante
pour un dimensionnement en dynamique.
Question 20 • Proposer une approximation de la pulsation fondamentale par équivalence avec un système à un degré
de liberté.
Question 21 • Déterminer une nouvelle approximation de la pulsation à l’aide du second modèle, dans lequel la partie
de l’axe située dans le volant est supposée indéformable.
avec :
πd 4 πr14
I= =
64 4
donc : s
12πEr14
ω2 =
m(a − e)3
On retrouve le premier modèle quand e = 0.
Question 24 • Comment améliorer la qualité des approximation ci-dessus. Proposer une amélioration pour le premier
modèle.
ω2 ≤ R(U(x))
pour toute forme U(x) cinématiquement admissible. Avec R(U) le quotient de Rayleigh :
Z a
EI(U”)2 dx
0
R(U) = Z a
ρSU 2 dx + mU(a/2)2
0
U(x) = x(x − a)
– une déformée statique (par exemple la déformée sous l’action d’une force au centre) :
F
y(x) = (3a2 x − 4x3 ) → U(x) = (3a2 x − 4x3 )
48EI
(valable sur la moitié gauche la poutre avec origine en A).
– un mode de vibration analytique d’un structure proche. Par exemple celui de la poutre sur deux appuis sans
masse :
πx
U(x) = sin( )
a
Par exemple, avec la déformée statique :
Z a/2
2 EI(24x)2 dx
0 48EI
R(U) = Z a/2
= 17 4 3
2 ρS(3a2 x − 4x3 )2 dx + ma6 35 ρSa + ma
0
Résultats
Pulsation (rd/s) r1 = 10mm r2 = 20mm
1DDL modèle 1 654 2632
1DDL modèle 2 2402 9672
Rayleigh modèle 1 654 2632
Cela valide le calcul réalisé avec l’équivalence avec un système 1ddl. Car l’inertie de l’arbre est négligeable devant
celle de la masse.
Question 26 • La liaison pivot étant composée de deux paliers de chaque coté du volant, en déduire les efforts dus à
l’inertie dans ces paliers.
mRθ̇2 mV 2
Finer = =
2 aR
Question 27 • Déterminer le moment cinétique au centre de gravité du volant en faisant l’approximation gyroscopique.
Question 28 • En appliquant le théorème du moment dynamique, déterminer le moment dans la liaison pivot induit
par la rotation du volant sur son axe (couple gyroscopique).
Question 29 • La liaison pivot étant réalisée par deux paliers distants d’un entraxe a. Déterminer les efforts corres-
pondants dans les paliers.
2Jωθ̇ 2JωV
Fgyro = =
a aR
Dans la suite, on ne prend en compte que les efforts dus au couple gyroscopique.
Question 33 • Calculer la pression de contact due aux efforts gyroscopiques. Pour cela, on utilise la surface projetée
égale au produit de la longueur du palier par le diamètre. Le critère en pression est-il respecté ?
Question 34 • Calculer la vitesse de glissement au contact. Le critère en vitesse maximale est-il respecté ?
5 Autres solution
Question 36 • Proposer d’autres solutions de guidage pour satisfaire au cahier des charges.
Question 37 • Proposer des solutions pour limiter les pertes dans le système.
∂f 1 ∂f ∂f
grad f = er + eθ + ez
∂r r ∂θ ∂z
1 ∂ur
∂ur ∂ur
− uθ
∂r r ∂θ ∂z
1 ∂uθ
∂uθ ∂uθ
grad U =
+ ur
∂r r ∂θ ∂z
1 ∂uz
∂uz ∂uz
∂r r ∂θ ∂z
∂2 f 1 ∂f 1 ∂2 f ∂2 f
∆f = + + + +
∂r2 r ∂r r2 ∂θ2 ∂z2
1 ∂uz ∂uθ ∂ur ∂uz ∂uθ uθ 1 ∂ur
rot U = − er + − eθ + + − ez
r ∂θ ∂z ∂z ∂r ∂r r r ∂θ
2 ∂uθ ur 2 ∂ur uθ
∆U = ∆ur − 2 − e + ∆uθ + 2 − e + ∆uz ez
r ∂θ r2 r r ∂θ r2 θ
∂σrr 1 ∂σrθ ∂σrz σrr − σθθ
div σ = + + + er +
∂r r ∂θ ∂z r
∂σθr 1 ∂σθθ ∂σθz σrθ
+ + +2 eθ +
∂r r ∂θ ∂z r
∂σzr 1 ∂σzθ ∂σzz σzr
+ + + ez
∂r r ∂θ ∂z r
$$1.--
$$5.--
$$0.--
$$4.--
$$/.--
$$3.--
$$-.--
73 7/
$!3.--
)#*
Figure 4 – Résultats éléments finis axi-symétrique : σzz et variations le long du rayon
!"#"$%&'"(
r1 = 10mm, r2 = 150mm, e = 100mm, acier, ω = 19000tr/mn
Modèles dynamiques
nertie fonctionne entre une vitesse maxi- ω
dϕ
σtdr
les capacités mécaniques du disque en ri dr dFc (σr + dσr) (r + dr) dϕ
minimale ωmin au-dessous de laquelle la r
erd son efficacité.
r0 σrrdϕ
’écrit donc : σtdr
1
= J (ω max
2 − ω min
2 )
h
2
Figure 6 – Premier
Figure 2 – Représentation mod
d’un disque enèle dynamique
rotation en fonction
des contraintes et de l’effort appliqué [11]
’inertie
Cette caractéristique se retrouve dans les matériaux composites
avancés à base de fibres de verre ou de carbone, mis au point
ons sur les matériaux pour les besoins des industries aéronautiques et spatiales.
eau tournant de rayon R, la masse est On note toutefois que, contrairement aux matériaux métalliques,
et le moment d’inertie s’écrit : qui présentent des caractéristiques homogènes indépendantes de
la direction de la contrainte subie, les matériaux composites
J = mR 2 présententFigure 7 – Second
leurs meilleures modèle dynamique
caractéristiques mécaniques dans la
direction des fibres, les caractéristiques mécaniques dans les
pression dans l’équation de l’énergie, on directions transverses étant beaucoup plus faibles.
1
Ec = mR 2ω 2
2 2.2 Équation différentielle de base
épend de la masse linéairement et de la
d’un disque en rotation
u carré. La contrainte dans l’anneau en
Le paragraphe 2.1 a abordé quelques considérations théoriques
Figure
sur l’énergie 8 – Troisidans
emmagasinée èmeunmod èle dynamique
anneau en rotation. Mais les
σ = ρω 2 R 2 volants d’inertie se présentent sous la forme de disques épais et
leur dimensionnement mécanique est plus complexe. Les
e du matériau. méthodes d’analyse par éléments finis s’appliquent bien entendu.
Modélisation du véhicule en courbe
inte maximale admissible du matériau,
: Dans ce paragraphe vont être donnés des éléments de dimen-
sionnement analytique permettant une première étude, tant en
1 σ max matériaux homogènes (métaux) que composites (anisotropes).
Ec = m
2 ρ La relation d’équilibre d’un disque en rotation en fonction des
Ec σ max contraintes et de l’effort appliqué est représentée figure 2 [11]. Elle
= peut s’écrire selon l’équation suivante (r : direction radiale ; t direc-
m 2ρ
tion circonférentielle) :
que l’énergie spécifique (rapportée à la
est d’autant plus importante que la
missible est élevée et la densité du dϕ
(σ r + dσ r )(r + dr )dϕ − σ r r dϕ − 2σ t dr sin + ρω 2 hr 2 dϕ = 0
montre le tableau 3. 2
Ce document a été délivré (30/01/2012 08:03) pour le compte de 7200039662 - ens cachan
(1) Les types NJ et NUP admettent une charge axiale faible (§ 4.4).
(2) N vitesse de rotation (en tr/min).
Dm diamètre moyen de roulement d +D
-------------- (en mm). ! 2 "
Le facteur de vitesse limite est donné pour une lubrification à la graisse. En lubrification à l’huile, la vitesse limite est augmentée d’environ 35 %.
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
© Techniques de l’Ingénieur, traité Génie mécanique B 5 370 − 5
tiwekacontentpdf_b5370 Ce document a été délivré pour le compte de 7200039662 - ens cachan // [Link]
et 5.1.4]
5.1.2 et ■ Capacités (dégrossissage des possibilités d’emploi) :
● p max = 10 à 20 MPa pour les matériaux métalliques,
classes CRITÈRES DE 50
CHOIX
MPaD’UN PALIER.
pour lesVALIDATION
plastiques.______________________________________________________________________________________________
(0) ● v max = généralement 2 à 3 m/s, voire 5 m/s en lubrification à
la graisse, 10 à 20 m/s en lubrification à l’huile.
● pv (MPa · m/s) en marche continue :
Fréquence gazeuse,
— ●ambiance du graissage :
conditions moyennes de ventilation de
M3
1 à 2— produits
m/s : métalliques : de 1 h à 1 jour ;
— produits plastiques : d’une semaine à quelques mois, voire une
pv max = 0,1 k D · k T1 /f
année ou à vie.
avec k D de
■ Calculs coefficient de (après
confirmation dimension (figure
choix du 1a ),: [B 5 330] § 4.3.
matériau)
k T1 coefficient de température (figure 1b ),
f Ambiance
3.1.1 coefficientgazeuse,
de frottement dynamique
graissage admis. ;
Sous liquide
— ambiance carter :mais hors projections d’huile
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
B 5 310 − 4 © Techniques de l’Ingénieur, traité Génie mécanique
tiwekacontentpdf_b5310 Ce document a été délivré pour le compte de 7200039662 - ens cachan // [Link]