REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE L’ENSEIGNEMENT SUPERIEUR ET DE LA RECHERCHE SCIENTIFIQUE

UNIVERSITÉ DE SAAD DAHLAB, BLIDA
FACULTE DES SCIENCES DE L’INGENIEUR
DEPARTEMENT D’AERONAUTIQUE

MEMOIRE DE FIN D’ETUDES

EN VUE DE L’OBTENTION DU DIPLOME D’INGENIEUR D’ETAT EN AERONAUTIQUE

OPTION : STRUCTURE

THEME

VIBRATION DE TORSION DES

ARBRES DE MACHINES

Présenté par: Encadré par:

Doumi Noureddine MR Kirade Abedelkader

Année universitaire 2010-2011

 ‫ﻣﻠﺧص‬
‫ﯾﮭدف ھدا اﻟﻌﻣل إﻟﻰ دراﺳﺔ اﻻھﺗزازات اﻟﻧﺎﺟﻣﺔ ﻋن اﻵﻻت اﻟدوارة‬
‫وطرﯾﻘﺔ ﺗﻔﺎدي اﻟﺣﺎﻻت اﻟﺧطﯾرة)اﻟﺗﺣطم( و ذﻟك ﺑﺈﯾﺟﺎد ﺗواﺗر اھﺗزاز اﻟﻔﺗل‬
‫ﻟﻠدوار وﻣﻘﺎرﻧﺗﮫ ﺑﺎﻟﺗواﺗر اﻟﺣدي اﻟﻣﺳﻣوح ﺑﮫ‬

Résumé
Le but de notre travail est d’étudier les vibrations créent
par les machines tournantes et comment éviter les cas
dangereux (destruction de la machine) à partir de
détermination des fréquences propres de vibration de
torsion et prévenir le risque dangereux (fréquence
critique)

Summary
The traget of this work is for study of the vibrations
resulted from the rotative machines and avoiding the
dangerous cases (disaster) through determining of
frequency of torsion vibration and preventing the
dangerous cases (critical frequency)
Je dédie ce travail de fin d’étude avant tout a ma mère et
mon père qui m’aident pendant toute ma vie

Puis à mon seul frère ABDERRAZAK, à mes chères sœurs, a

toute ma famille.

A mes chers amis : SALEM, SAMIR, YOUCEF, DJAMEL,

HOUCINE, ABDELNASSER, AZZEDINE.

A ma promotion d’étude : HASSENE BELLALA, SABRI,

TAYEB BOUNADRI, WALID BOUSSAHA, SAID, FAROUK,
NASSIM……..etc.

Mon promoteur : KIRAD ABEDELKADER.

DOUMI NOUREDDINE
 SOMMAIR
-LISTE DES FIGURES

-RESUME

INTRODUCTION GENERALE ………………………………………………………………………………1

CHAPITRE I : Introduction sur les vibrations des arbres

1. GENERALITES SUR LA VIBRATION……...............................................................2

1.1 La recherche actuelle sur les vibrations ………………………………………………..4

1.2 Bref historique ……………………………………………………………………………………..4
1.3 Importance de l’étude des vibrations………………………………………..………….5

1.4. L’origine de la vibration……………………………………………………………………………..6

a) Les Mauvais alignements…………………………………………….………………6

b) Les balourds……………………………………………………….……………………….6
c) Les fixations……………………………………………………….……………………….7

1.5. Les Différents types de vibration……………………………………………………………....7

a)Vibration relative d’arbre…………………………………………………………………..7

b) Vibration absolues d’arbres…………………………………………..………………….8

c) vibration dues au déséquilibrage……………………………………………………….8

d) vibration dues aux excentricités………………………………………………………..9

1.6. L’effet des vibrations sur une structure………………………………..……………………9

2. ARBRES ET APPUIS D’ARBRES…………………………………………..……………………...12

2.1. Arbres…………………………………………………………………………………………………….12

2.1.1. Généralités………………………………………………………………….……………………….12

2.1.2. Constructions des arbres……………………………………………………..……………….15

2.2. Appuis d’arbres ……………………………………………………………………………………….16

 a)-paliers lisses…………………………………………………………………………………………....16

b)-palier a roulements………………………………………………………………..…………..…16

Chapitre II calcul statique des arbres (torsion)

Généralités…………………………………………………………………………..…………….18

1- Contraintes et déformations en torsion……………………………………………………19

2-Valeurs de la contrainte……………………………………………………………..…………….19

3-Moment quadratique polaire…………………………………………………..……………….22

a) section circulaire pleine………………………………………………..………………….22

b) section circulaire creuse…………………………………………………….…………….22

4- Calculs de résistance en torsion………………………………………………..……………..22

5-condition de déformation……………………………………………………………..………….24

Chapitre III modélisation dynamique des arbres

Vibration de torsion……………………………………………………………………………26

1-Modélisations des lignes d’arbres………………………………………………..……………..26

2- Modélisation des arbres coudes……………………………………………..………………….28

2-1- Rigidités torsionnelle du coude………………………………………………..……………..28

2-2 -Longueurs équivalentes d’un coude…………………………………………..……………29

2-3- Moment d’inertie d’un coude………………………………………………………….……..30

3- Moment d’inertie Equivalent de l’équipage mobile……………………………………31

a) masse équivalente du piston et de son axe……………….……………..32

b) Masse équivalente de la bielle………………………………….………………33
c) Moment d’inertie équivalente……………………………………….…………33
d) Modèle d’un coude et des pièces mobiles………………...…..…………34

4-Modélisation des arbres non coudes (mise en évidence des vibrations de

torsion)…………………………………………………………………………………………………………..34
4-1-Arbre a une section constante portant un disque…………………………………….34

4-2-Arbre a une section variable portant un disque………………………………………..37

4-3-Arbre a une section constante portant deux disques…………………..…………..39

5-Modélisation des accouplements élastiques…………………………………..…………..42

6-Modélisation des réducteurs……………………………………………………………………….42

CHAPITRE V RECHERCHE DES PULSATIONS PROPRES ET DEFORMES MODALES

INTRODUCTION…………………………………………………………………………………………….44

1-HYPOTHESES……………………………………………………………………………….………………45

2- Equations du mouvement du disque………..………………………………………………..45

3- Méthodes de résolution du système…………………………………….…………………….48

3-1-ALGORITHME DE RUTSHAUSER………………………………………………………………..48

a)-Propriétés des matrices tri diagonales………………………………………………..……..48

b)-Algorithme…………………………………………………………………………………………………48

3-2-Application de l’algorithme de RUTISHUSER…………………………………………….50

a)-Détermination des matrices [R], [S]…………………………………………....……………50

b)-Calcul des éléments de la matrice [C]=[S] ∗ [R]…………………………….…………..50

c)-Relation entre les éléments de [C]et[C° ]………………………………..……….…………51

3-3-TEST DE CONVERGENCE QUADRATIQUE…………………………………………………..51

a)Rapidité de la convergence………………………………………………………………………….51

b)-Algorithme de la méthode de RUTISHAUSER accélérée…………………..………...52

3-4-DEFORMES MODALES……………………………………………………………………..……….52

4-VIBRATION DE TORSION DE L’ARBRE (PRISE EN COMPTE LA MASSE DE

L’ARBRE)…………………………………………………………………………………………………………54

4-1-MOUVEMENT DE TORSION…………………………………………………..………………….54
a)-Equation du mouvement………………………………………………………………..………….54

b) -Résolution mathématique de l’équation de Mouvement…………………….…...56

4-2-MATRICES DE TRANSFERT DUES A LA VIBRATIONDE TORSION….................58

a)-Etablir la matrice de transfert…………………………………………………………………….59

b) Matrice de transfert d’un disque…………………………………………………...............60

4-3-DETERMINATION DES FRQUENCES PROPRES DE TORSION DU SYSTEME

(ARBRE + DISQUE)……………………………………………………………................................61

a)Relation de transfert de premier élément……………………………………................62

b) Relation de transfert de second élément…………………………………………………..62

c)Relation de transfert globale du système…………………………………………………….62

d) Equation de fréquence……………………………………………………………………………….63

5-PARTIE CALCUL…………………………………………………………………………................66

5-1-ORGANISATION DE CALCUL……………………………………………………………………66

5-2-DISCRIPTION DU PROGRAMME RUTHISHAUSER……………………………………..66

a)Lecture des donnés………………………………………………………………………………….....65

b) Calcul des éléments…………………………………………………………………..……………….65

c)Calcul des valeurs propres……………………………………………………………………..…...65

d) Calcul des pulsations propres……………………………………………………………………..66

e)Calcules des déformes modales…………………………………………………………………..67

5-3-ORGANIGRAME………………………………………………………………………………………..68

5-4-APPLICATION DE CALCUL………………………………………………………………………….73

5-4-1Description de l’installation…………………………………………………..……………….73

5-4-2 Modélisation de l’installation………………………………………………………..……..74

5-4-3Caractéristique inertielles et torsionnelles du modèle …………………….…...76

5-4-4-Détermination des pulsations propres et déformes modales…………….....76

5-5-PROGRAMME : calcul des pulsations propres et déformées modales des

Lignes d’arbres non ramifiées en vibration de torsion……...............................77

Affichage des résultats………………………………………………………………………….……….81

Conclusion………………………………………………………………………………………………………83

-Liste des symboles

-BIBLIOGRAPHIE
LISTE DES SYMBOLES :

A l’air de section

b largeur de bras

d diamètre

E Module de Young

G Module d’élasticité transversale

h hauteur du bras

V vitesse du piston

I Moment d’inertie du disque

K Rigidité torsionnelle du coude

N Nombre de disques

T temps

X position
Mb Masse rotative de la bielle

Mp Masse du piston

M Masse du maneton

IC Moment d’inertie du coude

It Moment d’inertie de tourillon

I Moment d’inertie du maneton

Ieq Moment d’inertie équivalent

I Moment d’inertie du bras

Moment quadratique de la section du

maneton

Moment quadratique de la section du

tourillon

J Moment quadratique polaire

L Longueur équivalente
L Longueur du tourillon

L Longueur du maneton

Km Rigidité torsionnelle du maneton

Kt Rigidité torsionnelle du tourillon

K Rigidité torsionnelle du bras

w Moment résistant polaire

[M] Matrice carrée d’inertie

[K] Matrice carrée des rigidités

[C ° ] Matrice de transfert de la poutre

négligent sa masse

D Matrice de transfert du système

globale

[D ] Matrice de transfert du disque

 SYMBOLE GREQUE :

Angle de rotation de la manivelle

comte à partir de
l’axe longitudinale
Angle de déviation de la fibre
Μ moment d’inertie de masse d’arbre
par rapport à l’axe zz’ par unité de longueur
angle de torsion produit par
l’application de deux couples
ξ masse volumique
σ contrainte normale
Contrainte normale admissible
τ contrainte tangentielle
Contrainte tangentielle admissible
Pulsation propre
π 3.141592654
LISTE DES FIGURES :
FIGURE : PAGE :
1-1 Vibration absolus d’arbres 11
1-2 Différentes types d’arbres 14
1-3 Quelques formes d’arbres 14
2-1 Diagrammes de contrainte tangentielle 20
2-2 Schémas géométrique de la déformation de l’arbre 21
2-3 Section d’arbres 22

3-1 Exemples de modèle de la ligne d’arbre 27

3-2 Coudes de vilebrequins soumis à deux couples 29

3-3 Dimensions d’un coude de vilebrequins 30

3-4 Modèle à un disque 31

3-5 systèmes de bielle manivelle 32

3-6 Arbre perpendiculaire porte un disque 35

3-7 Arbre perpendiculaire à section variable porte un disque 37

3-8 Arbre perpendiculaire à plusieurs sections variables portent un

disque 38

3-9 Arbre portant deux disques 39

3-10 Arbre portant deux disques (schéma de dimensions) 41

3-11 Modélisation d’un accouplement 42

4-1 Notation du modèle 45

 Introduction générale

Le but de notre présent mémoire intitulé "vibration de torsion des arbres de

machines" est d'analyse le problème crée par la vibration de machine.
Le premier chapitre de ce projet se limite à deux parties, le premier, comportant une
introduction consacrée aux phénomènes vibratoires et leur importance; Nous avons exposé
une étude détaillée et essentielle sur les causes ou l'origine de ces vibrations.
La deuxième donne les définitions principales des arbres et ses appuis et aussi les différents
types d'arbres que l'on trouve dans les machines tournantes.
D’une manière générale avant de passer aux calculs statiques et dynamiques, ce là facilité la
compréhension des considérations théoriques.
Les arbres doivent être calculés pour remplir un certain nombre de conditions:
- résistance mécanique: en aucun point la contrainte doit dépasser la résistance pratique du
métal.
- Rigidité à la torsion: elle doit être telle que les déformations élastiques, dues aux charges
appliquées, restent inférieures aux déformations maximales admissibles sous peine
d'engendrer des vibrations et un guidage défectueux.
Ces calculs statiques sont décrits dans le deuxième chapitre
Dans le troisième chapitre modélisation dynamique des arbres des machines, On se
propose pour cela d’établir les modèles dynamiques et les équations mathématiques.
Dans le quatrième chapitre : recherche des pulsations propres et déformées modales
Est le vif du sujet proposé, étant un sujet très vaste, qui demande un temps assez important et
des moyens considérables, pour cerner et traiter dans son ensemble. Pour notre part, nous
avons limité le sujet sur: une étude des vibrations de torsion des arbres de machines.
Dans ce chapitre, nous avons été guidés par deux cas: D'une part, on néglige la masse
de l'arbre, et d'autre part, tenant compte la masse de l'arbre, on aborde le problème avec une
autre méthode approchée et ce, en utilisant la méthode des matrices de transfert.
A la fin de ce chapitre on fait une application sur une installation ; moteur, arbre,
hélice d’un avion léger.
Un programme informatique pour le calcule des fréquences propres et déformes modales.

Page 1
1. GENERALITES SUR LA VIBRATION :

Depuis le début du vingtième siècle, les études et les recherches sur les machines, ont
participés à l'évolution de la science et des techniques modernes, ces machines sont
destinées à rendre le travail de l'homme plus productif et plus facile.
Sur le plan mondial la conception, la réalisation, et la commercialisation de ces
machines tournantes telles que: les moteurs, les compresseur, les turbines, réducteurs, rotor
d'hélicoptère, et hélices, ..., etc., deviennent de plus en plus importantes, leur utilisation
dans les différents. Secteurs à dépassé les prévisions de la statique, grâce à leur avantage,
leur développement et leur précision, bien que leur prix atteigne une valeur très élevée par
rapport à d'autres moyens.
Parmi les composantes essentielles des machines tournantes, les appuis d'arbres, les
arbres représentent les organes les plus vulnérables de ces machines, et exigent une
attention plus particulière, afin de maîtriser leur choix, leur paramètre d'exploitation, leur
comportement dynamique ... etc.
Donc pour un bon fonctionnement de ces machines, il faut que les arbres soit sous un
contrôle continu, mais malheureusement peuvent être touchés par des problèmes pires, qui
ne sont pas faciles à résoudre, parce qu'ils influent directement sur leur fonctionnement et
peuvent les détériorer.
Généralement, toute machine animée d'un mouvement de rotation est plus exposée
aux problèmes d'exploitation, alors il est plus haut intérêt pour le constructeur et l'utilisateur
de connaître les causes pour supprimer ou tout au moins limiter ces problèmes, afin de
prolonger la durée de vie de l'arbre et d'éviter tout incident grave (la résonance).
Pour augmenter le rendement et la durée de vie de ces machines, il est nécessaire d'éliminer
tous les problèmes d'exploitation, qui peuvent surgir au cours de fonctionnement.
Parmi les nombreux problèmes, qui peuvent influer sur la durée de vie et le mode de
fonctionnement des machines tournantes, on peut citer quelques problèmes qui sont plus au
moins importants:
- Problèmes de vibrations
- Problèmes de température

Page 2
- Mode de graissage et les incidents
Mais ce qui nous intéresse dans notre projet sont les problèmes de vibrations.

Page 3
1.1. LA RECHERCHE ACTUELLE SUR LES VIBRATIONS:

1.2. Bref historique:

L'homme s'est intéressé au phénomènes vibratoires lorsqu'il a construit les

premiers instruments de musique ou la recherche de lois de la production et transmission
du son a permis la réalisation de ces instruments.
C'est ainsi que PYTAHAGORE (582-507 AV JC) à prouvé expérimentalement que
Si deux cordes sont identiquement tendues, les tons qu'elles produisent différent d'un
Octave quand la longueur de l'une est le double de l'autre.
Malgré la connaissance acquise par les anciens, il faut attendre le début du 17éme
siècle pour que GALILEE (1564-1642) démontre que le ton d'un son est déterminé par la
fréquence de sa vibration.
Ces résultats furent confirmés pour la première fois par BUICK TAYLOR
(1685-1731)
Plusieurs mathématiciens renommés ont également étudies le problème de la corde
vibrante parmi les quels nous citons, Bernoulli, d'ALEMBERT, EULER, LAGRANGE,
ROURIER, ... etc.
Leurs études ont montrés qu'une corde peut vibrer latéralement de différentes façon
appelées" mode de vibration "le premier mode correspondant a la fréquence la plus
basse.
L'équation différentielle régissant les vibrations d'une corde fut établie par
ALEMBERT vers 1750, elle porte aujourd'hui le nom de " équation d'onde"
De son coté FOURIER établi en 1822 son mémoire célèbre sur la théorie des séries
harmoniques.
A partir de la loi de HOOKE, EULER, et BERNOUILLI ont étudié les vibrations des
poutres. Leurs méthodes de calcul, basée sur la conservation d'énergie fut reprise par LORD
RAYLEIGH (1842-1919), elle porte depuis son nom. Parmi les chercheurs contemporains
mentionnons Stodola (1859-1942) qui a établi une méthode d'analyse des vibrations des
poutres lors de ses travaux sur les vibrations des aubes de turbines.
Mais l'histoire des phénomènes vibratoires qui concernait au début la création et la
Page 4
transmission du son a pris un nouvel essor lorsque le courant alternatif fut adopté comme
agent de transport d'énergie; Des lors, l'étude des phénomènes oscillatoires acquit un intérêt
exceptionnel qui n'est donc pas étonnant que les électriciens qui prenant alors la relève des
mécaniciens et physiciens, aient mis au point les méthodes de calcul les plus commodes et
les plus fructueuses.
Ces méthodes ont été transposées depuis à l'étude des vibrations mécaniques et
acoustiques. Il est ainsi dégagé progressivement une théorie générale des vibrations
indépendantes de leur support physique.
1.3. IMPORTANCE DE L'ETUDE DE VIBRATION :

Les vibrations sont des phénomènes gênants apparaissent sur les machines, elles
s'influent directement sur leur durée de vie et sur leur fiabilité et diminuent aussi très souvent
leur précision, leur performance et également leur prix de vente.
Les vibrations mécaniques ont ordinairement un caractère parasitaire, qu'on traitait
auparavant comme un inconvénient mineur, ainsi pour éviter les résonances on se contentait
de se tenir assez loin des conditions dites " critiques ". Il est d'ailleurs reconnu aujourd'hui
que même tenues loin de la résonance, les vibrations altèrent la structure des matériaux, et
par suite leur endurance. Il est donc très important de savoir étudier les vibrations et de
connaître les méthodes de réduction de leur effet nuisible.
C'est à l'occasion de la rupture causées par les vibrations sur les vilebrequins des
moteurs, des ailettes de turbines, des avions en vol que les chercheurs ont d'abords élaboré
des méthodes particulières à chaque situation et c'est l'accumulation de l'expérience dans ce
domaine qui a rendu possible la création d'une- théorie générale sur les vibrations.
Le moteur d'avion comprend un grand nombre de masses liées par des éléments de
différentes rigidités, représentant ainsi un système compliqué, soumis aux vibrations
torsionnelles.
Les vibrations à la torsion du moteur diminuent la longévité des éléments et provoquent par
fois leur détérioration. De même, le confort de l'appareil se voit altérer par ces vibrations et
les bruits qui les accompagnent.
De nombreuses machines ou élément de machines sont sujet durant leur
fonctionnement à des vibrations. Les causes de ces vibrations sont assez variées: Mode De
Page 5
fonctionnement de la machine, imprécisions d'exécution ou de montage, usures, défaut de
fonctionnement...
L'étude des problèmes de vibrations a acquis une importance de plus en plus grande dans
les calculs des pièces de machines.
C'est un fait bien connu, que seule la théorie des vibrations permet de comprendre
entièrement nombres de problèmes, dont l'importance pratique est considérable, tels que:
équilibrage des machines, vibrations de torsion des arbres et des trains d'engrenages,
vibrations des disques et des Aubes des turbines, mouvement tourbillonnaire de vibration
des arbres.
L'étude des vibrations dans les moteurs des avions ct la recherche des moyens,
pour leur diminution, contribuent à l'amélioration de la fiabilité, il la longévité ainsi que du
confort, exigences contemporaines auxquelles on accorde beaucoup d'importance à
l'heure actuelle, lorsque les vibrations dépassent les limites admises il se pose la question
de l'origine de ces vibrations ce n'est que lorsque les excitateurs et les causes sont connues
que l'on peut éliminer efficacement les vibrations, c'est pour cette raison que le point de
départ de toute décision est la recherche des causes ou l'origine de ces vibrations.

1.4. L'ORIGINE DE LA VIBRATION:

Toute machine en fonctionnement est source de vibrations et de bruit, les

vibrations observées sur une machine peuvent être provoques par de nombreux
phénomènes, les plus importants sont les suivants:
a)- les mauvais alignements: en pratique plus de 50 % des machines ne sont pas
alignées, les raisons de ce mauvais alignement sont multiples:
- Echauffement différentiel des éléments de la machine
- Mauvais serrage
- Dilatation et usure des paliers, des arbres, ... etc.
Le mauvais alignement se traduit par des vibrations axiales de fréquence double de la
fréquence de rotation.
b)-les balourds : parmi les principales causes des balourds, on cite: - Les dissymétries de
conception des pièces

Page 6
- Les hétérogénéités de matières
- Les défauts d'usinage, les assemblages dissymétriques
- Les conditions de fonctionnement
Les balourds ne peuvent jamais être éliminés entièrement, compte tenue de
vitesses de rotation actuelles, le problème de l'équilibrage devient plus important encore.
c)- les fixations : la fixation des machines joue un rôle très important dans l'élimination ou
la production des vibrations.
Toute mauvaise fixation de la machine peut être due à une mauvaise assise, à un
châssis inadéquat, à une absence de support antivibratoires qui permettre à la transmission
des vibrations des machines voisines, on peut avoir d'autre causes de vibrations qui sont plus
au moins importantes.

1.5. LES DIFFERENTS TYPES DE VIBRATIONS DES MACHINES:

On définit ainsi trois types de vibrations des machines:

- Vibrations relatives d'arbres: ce sont le mouvement relatif de l'arbre par rapport à son palier
- Vibrations absolues de paliers: ce sont les mouvements extérieurs du palier mesurés par
rapport à un point fixe
-Vibrations absolues d'arbres: ce sont les mouvements de l'arbre mesurés par rapport à un
point fixe

a)- Vibrations relatives d'arbres:

On entend par vibrations relatives d'arbres, les mouvements rapides de l'arbre par
rapport à son palier, elles représentent les réactions du rotor aux forces alternatives agissant
sur lui par exemple un rotor balourdé est le siège de vibrations d'arbres, son centre de gravité
se déplace sur une orbite appelé " orbite cinétique " Les vibrations relatives sont jugées
d'après les deux critères:
- critère 1: Elles ne doivent pas dépasser un seuil de vibrations prédéterminé.
- critère 2_: la vibration de leur évolution dans le temps ne doit pas dépasser un
Seuil de gradient prédéterminé.

Page 7
b)- Vibrations absolues d'arbres:

JI est évident que des rotors légers tournant dans les paliers très rigides et très lourds ne
produiront que des vibrations de paliers faibles (FIG-1-1, a) c'est le cas de la plupart des
compresseurs rotatifs.
Lorsque la masse du rotor, dans une même structure d'accueil, augmente, le rapport,
entre les vibrations relatives d'arbres ct les vibrations absolues de paliers (FIG-1-1, b)
change. Dans le cas de très grosses machines, turbines par exemple.
Le rapport entre la masse du rotor et la masse du palier augmente encore (FIG-1-1, c).
Si les vibrations absolues de paliers sont plus grandes que vibrations absolues
d'arbres, les vibrations relatives d'arbres deviennent alors négligeables, on mesure alors de
préférence les vibrations absolues d'arbres.

c). Les vibrations dues aux déséquilibrages:

Equilibrer veut dire améliorer la répartition des masses du rotor de telle manière,
que les forces centrifuges libres f= m.r. Autour de l'axe de rotation imposées par la
construction ne dépassent pas les tolérances admises, de façon générale un rotor est
parfaitement équilibré, si son axe d'inertie principal se confond avec son axe de rotation.
Dans tous les autres cas, l'équilibrage dynamique du rotor est perturbé et le rotor est
affecté d'un balourd suivant la répartition de balourd, on distingue sur le rotor deux
balourds:
- défaut de centre de gravité : il provoque un déplacement parallèle à l'axe principal
d'inertie par rapport à l'axe de rotation
- l'axe principal d'inertie fait un angle par rapport à l'axe de rotation
NOTA:
1- Les disques ou tambours compresseur sont généralement équilibrés avant le montage
des ailettes
2- L'arbre compresseur est également équilibré dynamiquement avant le montage avec le
rotor (généralement par enlèvement de métal à des endroits bien définis)

Page 8
d). Vibrations dues aux excentricités:
Les rotors n'étant jamais parfaitement équilibrés, les arbres peuvent prendre en
rotation une certaine courbure appelée "excentricité"
L'augmentation de la valeur d'excentricité provient en général:
- Soit d'un mauvais alignement des paliers
- Soit de l'apparition d'un balourd sur un rotor
Les forces excitatrices agissantes sur le rotor contiennent les forces
aérodynamiques (hélicoptère), magnétiques (moteur) et mécaniques (machine tournante),
la plus importante de ces forces et la force centrifuge du balourd reproduisant le
déséquilibré Du rotor. La qualité de la balance a été déterminée par des expériences pour
différentes types de machines en fonction des tours et donnée par des tableaux ou des
graphes.
L'augmentation du niveau de vibration entraîne successivement une augmentation
de la valeur d'excentricité.
Les faibles jeux radicaux qui existent entre les parties fixes et les parties mobiles des
machines rendent indispensables le contrôle de l'excentricité de l'arbre.
En effet, lorsque la valeur de l'excentricité est élevée des vibrations importantes
prenant naissance des frottements entre les parties fixes et les parties mobiles peuvent se
produire.
* Effets de frottements:
Les frottements peuvent être dus à plusieurs causes dont les deux principales sont:
1- Frottements par manque de lubrification, il y a frottement de l'arbre avec le Palier, ce
phénomène produit des vibrations à hautes fréquence et à large bande qui ne sont pas liées
au régime de rotation de l'arbre.
2 - Les pièces tournantes présentent des aspérités ducs au mauvais usinage, ces aspérités
provoquent des frottements qui donnent naissance à des vibrations de hautes fréquences.

1.6. L'EFFET DES VIBRATIONS SUR UNE STRUCTURE:

Comme on la vu précédemment, le phénomène vibratoire est très souvent
préjudiciable aux structures lorsqu'elles sont soumises à des sollicitations dynamiques crées

Page 9
par des forces extérieures, ou intérieures, ou par des déplacements.
L'amplitude du déplacement engendrée est souvent beaucoup plus importante que
celle des sollicitations statiques et les contraintes induites peuvent entraîner la rupture.
Sur les aérodynes (avions, hélicoptères), les problèmes de vibrations se présentent de façon
très particulière; En effet (le moteur d'avion, le rotor d'hélicoptère) constitue un puissant
générateur de vibrations et sa présence créé des problèmes spécifiques à ce type d'engin.
Généralement pour résoudre le problème que créent les vibrations d'une structure
quelconque, procède comme suit:
- Déceler les sources de vibrations
- Rechercher les possibilités de résonance, et ce par la détermination des Caractéristiques
dynamiques de la structure c'est à dire les fréquences et modes propres.
L'estimation des fréquences propres permet d'éviter la coïncidence entre ces
dernières et les fréquences du système ce qui pourrait conduire au phénomène de résonance
qui entraîne le plus souvent la rupture de la structure.
Ainsi, tout signal de vibrations que peut avoir un système peut être inter
prête comme combinaison linéaire des modes propres de la structure.
En fonction des résultats trouvés on pourra imaginer les moyens d'atténuer
le phénomène résonance.
Nul n'imaginera donc, les conséquences désastreuses qu'engendraient les
vibrations sur une structure volante ct ces répercussions dans le domaine
aéronautique.
Aussi, il n'est plus possible d'ignorer les phénomènes vibratoires, il faut
donc essayer de les prévoir pour dimensionner les structures au stade des bureaux
d'études et les modifier judicieusement, si une fois construites, elles ont des
caractéristiques dynamiques gênantes.

Page
10
FIG -1-1 –VIBRATION ABSOLUS D’ARBRE

Page
11
2. ARBRES ET APPUIS D'ARBRES:
Les montages mécaniques le plus courant et les plus anciens comportent deux
parties, un élément interne, l'arbre et un élément externe, le palier, on nomme tourillon la
partie de l'arbre qui tourne dans le palier.
L'utilisation des paliers dans les machines tournantes sert guider les arbres en
rotation, et pour reporter leur charge sur un appui fixe, tant que l'arbre joue le rôle d'une
colonne vertébrale dans la machine, ils sont utilisés dans des nombreux moteurs ou autre
dispositifs mécaniques.
2.1. ARBRES:

2.1.1. GENERALITES:

Les arbres son des pièces cylindriques, destinés a [Link] ces pièces
viennent se fixer des organes de transmissions (roues dentées, disques, ..., ex
L'arbre constitue le plus souvent la partie mobile de la machine
(Arbre d'un moteur par exemple).
La position importante des arbres dans les machines tournantes permet de transmettre
des puissances mécaniques reçues par l'intermédiaire des éléments moteurs, ou éléments
récepteurs.
Il existe plusieurs types d'arbres, les plus utilisés sont:
-arbres moteurs
- arbres de transmission
- arbres récepteurs (voir figure. 1.2)
a)-arbres moteurs: Il constitue l'organe essentiel de machines motrices ou
moteur dont le mouvement principal est un mouvement de rotation (arbre de turbine par
exemple)
b)-arbres de transmission : Ce sont des organes intermédiaires, recevant l'énergie
mécanique d'un moteur, ils la transmettent à une machine réceptrice (par l'intermédiaire
d'organes de transmission: accouplements, réducteurs) et ils peuvent porter engrenages,
roues dentées, etc.
c)-arbres récepteurs : Constituant l'organe principal des machines réceptrices,

Page
12
(compresseurs par exemple), ils reçoivent l'énergie mécanique, soit directement du
moteur, soit par l'intermédiaire des arbres de transmission et transmettent cette énergie
aux organes récepteurs.
Parmi les différentes formes d'arbres, on peut citer quelques exemples (voir figure 1.3)

Page
13
FIG-1-2 DIFFERENTES TYPES D’ARBRES

FIG - 1-3 QUELQUES FORMES D’ARBRES

Page
14
2.1.2. CONSTRUCTION DES ARBRES :

* Diamètre et vitesse de rotation :

Les diamètres des arbres de transmission sont normalisés (norme NF E 22.100) et égaux à
(en mm) 20, 25, 30, 35, 40, 45, 50, 60, 70, 80, 90, 100, 120 et 150
Leurs fréquences de rotation le sont aussi (norme NF E 22.003) et ont les valeurs suivantes
(en tours par minutes)

25 63 160 400 1000 2500

31.5 80 200 500 1250 3150

40 100 250 630 1600

50 125 315 800 2000

*Martiaux utilises :

a)-Aciers non allié :

Les aciers doux étirés (R=40 à 50 dan/mm2) sont employés sans usinage, grâce à leur faibles
tolérances, pour les arbres de transmission, les aciers mi-durs et durs, souvent trempés
(R=55 à 80 dan/mm2) sont utilisés, usinés, ou forgés, pour les arbres incorporés aux
machines.

b)- fontes ductiles a graphite sphéroïdal :

Alliées ou nom (R=65 à 100 dan/mm2), elles permettent l'obtention économique d'arbres de
formes complexes (vilebrequins creux à contre poids incorporés de moteur d'automobile
obtenus par fonderie)

c)-Aciers allies:

Au chrome ou au nickel chrome trempés et revenus (R ≥120daN/mm2) .il sont réservés

pour les fortes charges.

Page
15
2.2. APPUIS D'ARBRES:

Dans un ensemble mécanique, certaines pièces opèrent un mouvement relatif par

rapport à d'autre pièces, elles mêmes fixes ou mobiles, il importe que le guidage assure la
précision et la fiabilité de la trajectoire assignée à chaque pièce pour qu'elle remplisse sa
fonction. La qualité du guidage choisi est élément déterminant du fonctionnement correct et
du bon rendement de l'ensemble d'une générale, on distingue les guidages par glissement
(paliers lisses) ct les guidages par roulement (paliers à roulements).
Le choix de la solution la mieux adaptée n’est pas toujours simple car de
nombreux paramètres doivent être pris en compte, notamment la charge de base
dynamique de fonctionnement, l’éventuelle perturbation prévisible (vibrations,
chocs…etc.)

a). PALIERS LISSES:

Suivant l'orientation de l'arbre et le sens de la charge à la quelle, ils sont soumis, on

Distingue:
- Les paliers porteurs (axe horizontal, charge radiale)
- Les paliers de butée (axe horizontal, charge axiale)
- Les boitards(axe vertical, charge radiale)
- Les crapaudines (axe vertical, charge axiale et radiale)

b). paliers à roulement:

Le roulement est un organe mécanique destiné à permettre la transmission des efforts

entre deux pièces en rotation l'une par rapport à l'autre avec un frottement minimal.
JI se distingue par un faible encombrement axial, la précision des centrages, la
faiblesse des couples de résistance au démarrage et des coefficients de frottement en service.
Le roulement est normalement constitue de quatre éléments de base :
-une bague intérieur (solidaire a l’arbre)
- une bague extérieure
-des éléments roulants (permettant la rotation relative d’une bague par rapport à l’autre)
-une cage maintenant l’écartement des éléments roulants et les empêchant de frotter d’un

Page
16
contre l’autre.
On classe généralement le roulement en trois familles principales :
- Les roulements à billes
- Les roulements à rouleaux
- Les roulements à aiguilles
*Les différents montages de roulements:
Est utilisé deux types de roulements:
-Roulement à billes
- Roulement à rouleaux
Le roulement à bille: est destiné à encaisser les efforts axiaux de la turbine et de
Compresseur, c'est le roulement de la butée.
Le roulement a rouleaux : est destiné à encaisser les efforts radiaux de la turbine et de
compresseur.

*Exemples de montages:

a)- Réacteurs simples attelages:

b)- Réacteurs doubles attelages:

Page
17
Généralité :

Tous les arbres qu'on rencontre changent plus ou moins de forme (se
déforment) sous l'influence de l’action extérieure. L'intensité de ces déformations
dépend de la matière des arbres, de leurs formes géométriques, de leurs dimensions
et des efforts appliqués.
Lorsque l'on calcule la solidité des parties d'une construction ou d'une
machine, les déformations des arbres acquièrent une importance considérable. Ces
questions sont étudiées dans les cours de résistance des matériaux.
La résistance des matériaux est la science la plus générale qui étudie la solidité des
machines, la création de différents types de machines et mécanismes, avions,
hélicoptère, turbines et moteurs, ... etc.
Du point de vue physique, il' est clair que le matériau ne puisse pas supporter
des contraintes aussi grandes qu'on veut Aussi, les valeurs des contraintes
maximales doivent elles, en portant de la condition de fiabilité de la pièce au cours
de l'exploitation, être limitées par certaines quantités permises, ces quantités sont
dites contraintes admissibles.
En pratique, on rencontre 3 cas de calcul de résistance:
1°_ Les sollicitations étant connues, l'on demande de trouver, pour le matériau
choisi, les dimensions correspondantes de la section transversale de l'arbre
2°_ Les dimensions de la pièce et le matériau sont connus, l'on demande de vérifier
si l'arbre pourra supporter une sollicitation donnée.
3°_ Sont connus le matériau, la dimension de la pièce, ainsi que la façon dont elle est
chargée. L'on demande de déterminer la valeur admissible de cette charge.
A la base de tous ces calculs est opposée la condition de résistance:

Page
18
* Torsion: La torsion a lieu lorsque sur un arbre agissent des forées extérieures qui
engendrent un moment par rapport à son axe. La déformation de torsion
s'accompagne d'une rotation des sections transversales de l'arbre l'une par rapport à
l'autre autour de l'axe de cette dernière.
1. CONTRAINTES ET DEFORMATIONS EN TORSION:
On considère la fig. ci dessous, deux sections droites très voisines s et s' situées à la
distance dz. Une portion de fibre telle mn située à la distance φ que après déformation

élastique, en m'n', soit φ l'angle dont la section s’à tourné par rapport à S

Comme la déformation élastique est très petite, la figure m n' n" peut être assimilée à un
triangle rectangulaire par ailleurs, le glissement de la section s' par rapport s peut être
caractérisé par l'angle de déviation y de la fibre m'n'.
La déviation a donc pour mesure:

≃ = = . …………………………..( − )

2-VALEUR DE CONTRAINTE:

Le glissement dû à la torsion est analogue à celui provoqué par le cisaillement et la

contrainte tangentielle correspondante est donnée par l'expression:

=G.
D’après (3-1), on a:
=G. … … … … … … . (3 − 2)
Le rapport ⁄ représente l'angle de torsion rapporté à l'unité de longueur, il est appelé

Page
19
angle unitaire de torsion et il est désigne par θ. Exprimant le cisaillement' dans y les fibres
extérieures de l'arbre par les contraintes, on trouve en accord avec la loi de HOOKE, pour
le cisaillement, le rapport entre les contraintes tangentielles dans les fibres extérieures L et
l'angle de torsion θ

=G. .θ ……………………. (3-3)

D'après la formule (3-2), on voit qu'en torsion la déformation de cisaillement et les

contraintes tangentielles sont directement proportionnelles à la distance jusqu'au
centre de gravité de la section.

FIG -2-1

Le diagramme des contraintes tangentielles suivant la section droite de l'arbre est

représenté sur la (fig.:2-1)
Au centre de gravité d'une section circulaire les contraintes tangentielles sont nulles.
Elles sont maximales aux points de la section voisine de la surface de l'arbre.
Le moment de torsion en une section quelconque de l'arbre qui est la résultante de la
section simultanée des contraintes tangentielles agissant dans les plans
Élémentaires dA se trouvant à la distance ~ du centre de la section peut être exprimé
par l'équation:

…………………. (3-4)
La distribution des contraintes tangentielles L dans la section se détermine d'après
les schémas géométrique des déformations de l'arbre en torsion reproduit sur la
(FIG- 2-2)
Page
20
 FIG - 2-2

En mettant (3-3) dans (3-4), il vient que:

De là, on obtient une formule permettant de déterminer l'angle de torsion relatif de l'arbre:

= = = …………………..( − )
.

Où GJP - est la rigidité de la section transversale de l'arbre en torsion qui se

Mesure en Kgf.cm2
L'angle de rotation total d'un arbre de longueur L est égal à

=∫ . = . = .L………………… (3-6)
. .

Où ⁄ est la rigidité de l'arbre à la torsion qui se mesure en [Link]

En mettant la valeur de tirée de (3-5) dans (3-3), trouvons maintenant la contrainte
tangentielle 't en un point quelconque de la section de l’arbre :

.
= ……………………………………………….( − )

La contrainte tangentielle maximale est de toute évidence égal à : max= .

Ou bien à:
max= ……………………………………………..( − )

WP= ⁄ - est le moment résistant polaire.

Page
21
 3. MOMENT QUADRATIQUE POLAIRE :

a)- Section circulaire pleine:(fig.: 2.3,a)

On a: =
= = ≃ . . ……………( − )
.

b)section circulaire creuse : (fig2.3, b)

.
On a: =

( ) ( )
WP= =

WP=0.2.D3 (1-C4)…………………… (3-10)

Ou : =
FIG - 2-3

Donc, pour un arbre cylindrique plein, de diamètre d, le moment résistant polaire est
donnée par la formule (3-9) et la contrainte tangentielle maximale est

…… ……………. (3-11)
Pour un arbre cylindrique creux, wP détermine d'après la formule (3-10), la contraint
tangentielle maximale est :
.
=
( )
……………….( − )

4. CALCUL DE RESISTANCE EN TORSION:

Pour que la pièce résiste en toute sécurité à la torsion, il faut que la contrainte
maximale soit inférieure à la contrainte admissible du métal.
D'où la condition de la résistance statique de l'arbre à la torsion est sous forme

Page
22
 = ≤ … .. ..…………..( )

max - est la contrainte tangentielle admissible, sous l'action d'une charge statique;
On adopte (sans tenir compte de la concentration des contraintes et d'autres facteurs
qui Altérant la résistance).

=( . à . )

D'après la formule (3-13), le moment résistant de l'arbre en torsion doit être:

≥ …………….( − )

En plus vérification à la résistance, les autre quantités étant connues, on peut

également d'après cette formule (3-14) choisir le diamètre de l'arbre ou calculer le
moment de torsion.
Si on tient compte que pour une circulaire pleine w p =0,2d2

Il vient : ≥ = . . …………….( − )
. .

Pour les arbres creux, on l'avantage de réduire beaucoup le Poids de l'arbre moyennement
une légère augmentation du diamètre (arbre de turbine par exemple), ≌0.2.
Il vient:

………….………. (3-16)

Si l'on exprime le moment de torsion par la puissance N en ch. et le nombre de tours par
minute par n, on aura:

= . ( . . ) …………………..( − )

Et la formule (3-15) se transforme en:

≥ . . …………………..( − )
.

Page
23
Alors que la formule (3-16) s'écrit:

≥ . …………..………………….( − )
. ( )

La puissance N est donnée en Kilowatts (1 ch = 0,736 KW) le moment de torsion peut

exprimer par:

Mt= . . ( . . )……………………..( − )
.

5-Condition de Déformation:

Le diamètre d'un arbre de grande longueur, pour éviter une brusque détente en cas de
variation brutale du couple, sera calculé pour que la déformation angulaire entre les sections
extrêmes soit inférieure à une valeur maximale; si non l'arbre se comporterait comme un
ressort de torsion, telle que l'angle de torsion unitaire θ.
Reste inférieure à une valeur maximale; par exemple: 0,25 degré par mètre de
longueur, valeur qui est couramment adoptée en pareil cas :

θmax= ≤ …………………………………..( − )
.

D'ou le moment d'inertie polaire garantissant une rigidité admissible:

JP≥ ………………………………………………...( − )
.

Par conséquent le diamètre d'un arbre cylindrique pleine d, doit être:

.
d≥ …………………………………….( − )
. .

Tan disque le diamètre extérieur d'un arbre cylindrique creux D, pour un C donné:

.
d≥ ( ). .
………………. ………….. (3-24)

Page
24
Remarque:

Le calcul à la condition de déformation impose une valeur de diamètre nettement plus grande
que le calcul à la condition de sécurité. C'est évidement cette deuxième valeur du diamètre
qu'il faudra adopter. Surtout si l'arbre a une grande longueur.

Page
25
VIBRATIONS DE TORSION :

Lors delà conception, il est nécessaire de prévoir le comportement vibratoire des

structures en cours de fonctionnement afin d'augmenter la fiabilité et de prévenir les
ruptures.
L'ensemble mécanique dans l'avion (moteur - arbre - hélice par exemple) représente
un système déformable sous l'action des couples de torsion un tel système admet plusieurs
modes de vibrations libres.
Chacun de ceux ci étant caractérisé par une fréquence naturelle et des amplitudes
relatives pour chaque élément du système.
Le danger principal provient cependant d'une résonance qui risque de se produire
entre les vibrations provoquées par le couple moteur et la fréquence de vibration propre de
torsion de l'arbre; les vitesses de rotation aux quelles ces phénomènes de résonances se
produisent sont dites vitesses critiques de torsion.
Dans tel cas, le fonctionnement peut être très dangereux ct entraîne la rupture de
l'arbre ou l'usure rapide des autres éléments de machine.

1. MODELISATION DES LIGNES D'ARBRES:

Les recherches ayant trait aux vibrations des machines et des mécanismes, doivent
comprendre une étude théorique suivie d'une vérification expérimentale de la conclusion
théorique.
L'expérience à montré que dans la recherche théorique concernant un certain aspect
du système physique réel, les calculs peuvent être considérablement simplifies, ne
concernant qu'une partie des propretés (celles faisant l'objet de l'étude) et l'on néglige les
autre. Le résultat de cette simplification du système physique réel s'appelle modèle L'étude
des vibrations de torsion dans les transmissions des machines procède en général des
modèles linéaires par différentes méthodes à masses concentrées, reliées par les élasticités
non pesantes.
L'un des problèmes fondamentaux est la détermination des pulsations propres du
système, afin d'éviter les résonances aux régimes de fonctionnement.

Page
26
Caractéristique du modèle:
La complexité de la chaîne de transmission (moteur, réducteur, arbre de
transmission, hélice, par exemple) nécessiterait un nombre très élevé de degrés de liberté, ce
qui entraînerait un coût du calcul excessif. Ce qui par contre nous parait faisable est de
calculer.
En utilisant le modèle classique "masse ressort". Le modèle le plus couramment
utilisé pour les transmissions de mouvement de rotation, consiste en un ensemble de
solide de révolution, supposés infiniment régi des, réunis entre eux par des ressorts
de torsion sans masses, de rigidité torsionnelle Ki, comme indiqué sur la figure (3-1).
Ainsi, il est nécessaire de créer, pour une ligne d'arbre, différents modèles pour:
-Réaliser l'équilibrage d'un vilebrequin.
-Connaître les fréquences propres des vibrations de torsion, de flexion ou Axiales.
-Déterminer les amplitudes en quelques points particulières des oscillations forcés
de torsion; de flexion ou axiale.
-Donner une image de l'état de contrainte de la matière ... etc.
Nous nous limiterons à la caractérisation des modèles traduisant les N
premières fréquences naturelles de torsion et les amplitudes des oscillations forcées.
Les éléments de ces modèles doivent satisfaire à des règles bien connues:
Proportionnalité entre efforts et déformations; conservation des masses, pour que les
lois de la mécanique leur soient applicables et que la résolution numérique des
équations conduise à l'utilisation de schémas mathématiques classiques.
Ces modèles représentent une distribution de paramètres et de degrés de liberté qui
dépendent de l'expérimentation.

FIG-3-1-Exemple de modèle de ligne d’arbre

Page
27
2. MODELISATION DES ARBRES COUDES :

Le vilebrequin constitué par des coudes qui ont un mouvement de rotation

reliés à des pièces (bielles - pistons) animées d'un mouvement alternatif.
Si l'0n Veut remplacer le vilebrequin par un arbre cylindrique équivalent, il faut
commencer par étudier la rigidité ainsi que le moment d'inertie à la torsion.

2.1. RIGIDITE TORSIONNELLE DU COUDE :

Si l'on considère un coude de vilebrequin soumis à deux couples opposés "C" portés
par l'axe du tourillon et appliqués à chaque extrémité, la déformation globale en torsion est la
somme des déformations.
Cela se traduit par la rotation rotative des deux sections extrêmes du coude d'un angle 0
(rad).
La rigidité torsionnelle K du coude est: K = CI θ

Ou encore: = + +

Kt : rigidité torsionnelle d'un tourillon.

Km : rigidité torsionnelle d'un magnéton.
Kb : rigidité équivalente d'un bras.
Si dans un premier temps, nous faisons l'hypothèse que les formules simples de RDM
sont applicables aux coudes de vilebrequin, nous pouvons calculer la rigidité K en fonction
des dimensions:

. . .
Kt= , Km= , k m=

Page
28
 FIG - 3-2 -Coude de vilebrequin soumis à 2 couples C.

.
Jt= ( − ) , Jm= ( − ) , Jb=

dim: diamètres des alésages dans les tourillons et les manetons s'ils existent.
dt ,dm : diamètres respectifs des manetons et tourillons
Lb; b : dimensions de la section du bras.
D'ou:
. .
= + +
.

2.2. Longueur équivalente d'un coude:

Selon la formule classique de la R.D.M, réduisant le coude de vilebrequin à une suite

de poutre, la longueur équivalente Leq de l'arbre du diamètre équivalent en
Torsion Deq est égale:

Leq= .

.
Soit : Leq= [ + + . . ]
.

Un coude est un système complexe ct seule la prise en compte des formes réelles permet de
déterminer avec précision la longueur équivalente. Pour cela, plusieurs semi-
expérimentales ont été proposées pour la détermination de cette dernière, en prenant comme
diamètre de l'arbre équivalent au coude, le diamètre du tourillon de vilebrequin.

Page
29
- Formule de CARTER :

. . . .
Leq= .[ + + ]
. . .

- Formule de KER-WILSON:

Formule de TIMOSHENKO:

. . . .
Leq= .[ + + .

FIG -3-3- Dimension d’un coude de vilebrequin

2.3. Moment D'inertie D'un Coude:

Le moment d'inertie IC d'un coude de vilebrequin est donnée par:

IC=Im+2It+2Ib

Im : Moment d'inertie du maneton par rapport à l'axe du vilebrequin.

=( ζ )[ + . ]

ζ en (Kg/m3) : masse volumique du matériau.

It : Moment d'inertie d'un demi - tourillon par rapport à l'axe du vilebrequin
Ib: Moment d'inertie d'un bras par rapport à l'axe du vilebrequin.
On effectue un découpage du bras et de son contre poids en éléments simples; puis pour
Page
30
chaque élément, parés avoir calculé son inertie, par rapport à un axe Passant par son centre
de gravité, et parallèle à l'axe du vilebrequin, en appliquant lé théorème de
HUYGHENS.
La somme de ces inerties élémentaires est égale à l'inertie du bras par rapport à l'axe
du vilebrequin.

a)-Modèle a un disque:

Le modèle à un disque est constitué par un arbre sans masse de diamètre Deq et de
longueur Leq portant un disque d'épaisseur nul et d'inertie le

IC = It + Im+ Ib1+ Ib2

Avec : Ib1 ct Ib2 : moment d'inertie des bras 1 et 2
Leq1 et Leq2 : sont les longueurs équivalentes des deux
demi-coudes.

FIG-3-4 modèle a un disque

3. Moment d'inertie Equivalent de L'équipage Mobile:

Pour tenir compte de l'effet de l'attelage mobile sur le système élastique

(vilebrequin), nous allons calculer une masse équivalente constante liée au maneton
du coude.

Page
31
 FIG-3-5- SYSTEME DE BIELLE-MANIVELLE

a)- Masse équivalente du piston et de son axe:

L'énergie cinétique du piston à tout instant au cours de son de placement, doit être égale à
l'énergie cinétique de la masse équivalente attachée au mouvement du maneton. La fonction
de force due à la pesanteur est identique dans les deux cas :

= [ ( )⃗ ]

Avec: Mp(Kg) - masse du piston et de son axe

( )⃗ - vecteur vitesse du point A par rapport au repère fixe
L'énergie cinétique de la masse équivalente attachée au maneton est:

= ( ). .

Avec: Meq(α) (Kg) - masse équivalente du piston, fonction de l'angle de

Rotation (rd) de la manivelle.
R - rayon de la manivelle.
w = dα/dt (rd/s) - vitesse de rotation constante de la manivelle
D'ou:

La masse équivalente moyenne cherchée sera obtenue en calculant la moyenne arithmétique

de la fonction Mel' (a) pour a variant de 0 à 1t.

Page
32
 = ∫ ( ).

b)- Masse équivalente de la bielle:

On posant:

Mb(Kg) - masse de la bielle

α (rd)- angle de rotation de la manivelle.

w= (rd/s) - vitesse de rotation constante de la manivelle.

°
- vecteur rotation de la bielle par rapport au Repère fixe

( ⁄ ) -Tenseur d'inertie de la bielle par rapport à Repère lie à la bielle et d'origine

( )⃗ -vecteur vitesse de G par rapport au repère fixe

L'énergie cinétique de la bielle est :

= . ( ⃗) + ⃗ . (̅ ⁄ ). ⃗°
L'énergie cinétique de la masse équivalente ponctuelle Meb attachée au maneton est:

= ( ). .

La masse équivalente de la bielle, en fonction de est donnée par l'expression.

. ( )⃗ + °
. ( ⁄ ). °⃗
( )= .
.

La masse équivalente cherchée sera obtenue en calculant la moyenne arithmétique de la

fonction M ( ) pour variant de 0 à π.

. ( )⃗ ⃗ . ( ⁄ ). ⃗
M = ∫ d
.

Page
33
c)- Moment d'inertie équivalent:

Avec :
- Ieq : moment d'inertie équivalent (N.m. s2)
- Meq : masse équivalente du piston (kg)
-Meb : masse équivalente de la bielle (kg)
-R : rayon de la manivelle (m)

d). Modèle D'un Coude et des Pièces Mobiles Associées:

Modèle à un disque: Le modèle est identique à celui que nous avons défini au paragraphe
([Link], a). Dans ce cas, le disque a pour moment d'inertie
I =IC+Ieq

IC - Moment d'inertie du coude (N.m.s2)

Ieq- Moment d'inertie équivalent de l'équipage mobile.

4. MODELISATION DES ARBRES NON COUDES :

(Mise en évidence des vibrations de torsion)

4.1. Arbre a une section constante portant un disque:

Un arbre vertical AB dont l'extrémité B est fixé perpendiculairement au centre d'un disque
(FIG:3-6) peut également servir d'exemple d'un système à un degré de liberté,
Les vibrations de torsion se produit dans l'arbre AB quand le disque tourne autour de l'axe
AB

FIG : 3-6

Page
34
Le moment de torsion engendre dans l'arbre sera:
= − . …………………………..( − )
Où, K - la rigidité de torsion de l'arbre
Le moment d'inertie du disque par rapport à l'axe de rotation
= . ……………………….( − )
De l'équation (4-1) et (4-2), on peut établir l'équation du mouvement vibratoire en torsion de
ce système:

………………... (4-3)
ω - La fréquence propre tic vibration libre en torsion de l'arbre et disque

………………….(4-4)
I - moment d'inertie du disque
On peut trouver la solution générale de l'équation (4-3) sous la forme:
= . ( − )………( − )
ϕ - Amplitude angulaire tic la vibration de torsion déterminée par les
Conditions initiales.
Pour éviter les phénomènes résonance, il faut calculer les fréquences propres ω du système.
Si l'arbre examiné comme un système continu en comprenant sa masse, dans
Ce cas le système possède infiniment ω (ω1,ω2,ω3,…...ωn)
Si nous considérons l'arbre comme un système sans tenir compte sa masse, le
système n'existe qu'une fréquence propre de torsion œ calculée d'après la formule (4-4)
Pour calculer ω, il faut connaître les paramètres suivants:
………….……………. (4-6)

Jp - Moment d'inertie polaire tic la section tic l'arbre.

.
=

I= (Moment d'inertie du disque)

Page
35
Enfin, on peut déterminer la fréquence ω sous la forme :
.
ω= . ………………………………….( − )
.

Avec: G- Module de cisaillement

L- Longueur de l'arbre
d- Diamètre de l'arbre.

Page
36
4.2. Arbre a une section variable portant un disque:

En réalité, dans la pratique toujours on rencontre l'arbre constitué sous la forme

variable (fig.:3-7)

FIG -3-7

Dans ce cas, l'angle de torsion de ce système peut être déterminé

. .
θ= + ………………….( − )
. .

J , J - Sont les moments d'inertie des sections de l'arbre de diamètres

Respectivement.

= , = …….………( − )

De (4-8) et (4-9), on peut former:

. .
θ= ……………….……( − )

Où, Leq - est la longueur équivalente de l'arbre de deux tronçons de diamètre d1 et d2:

.
= = …………….( − )

Nous constatons qu'on peut transformer un arbre composé de deux tronçons

différents du diamètre à un arbre de diamètre d, et de longueur Leq

Page
37
La rigidité de torsion de l'arbre équivalente est déterminée:

= ……………………..( − )
.

Dans ce cas la fréquence propre du système (arbre + disque) calculée selon la

formule:
= =
. .

ω=d . . .
… … … … … … (4 − 13)

En général, si l'arbre est composé de n tronçons de longueurs 11etl2,….. ,ln et de

diamètres correspondants d1,d2,…… ,dn dans ce cas on peut établir la longueur
équivalente de ce système sous la forme:
. . .
= = + + ⋯+ ……( − )

De ce la, on peut facilement calculer la fréquence propre du système compliqué.

(FIG-3-8)

FIG- 3-8

Page
38
4.3. Arbre a une section constante portant deux disques:

Première Méthode :

Considérons un arbre portant deux disques (FIG-3-9), si l'on applique sur ces
disques deux couples de torsion opposés et qu'on les supprimés brusquement.
L'arbre sera soumis à des vibrations de torsion pendant les quelles les disques
tournent toujours en sens contraire (système à deux degrés de liberté)

FIG-3-9

En désignant: K - la rigidité torsionnelle de l'arbre.

I1 ,I2- moments d'inerties de deux disques.
, − La position angulaire de chaque disque par Rapport à leurs
positions de repos.
Supposons que le matériau composant le dernier, suit la loi de " HOOKE".
Nous avons: =− ( − )……………..( − )
Avec: Mt - couple de torsion de l'arbre
θ - angle de torsion de l'arbre
Négligeons l'inertie propre de l'arbre, appliquons la loi fondamentale de la
dynamique, nous permet d'écrire l'équation du mouvement de la vibration libre de
torsion:
"
1) . + ( − )……………….( − )
2) . + ( − )……………….( − )
- étant l'accélération angulaire du disque.
= − -L'accélération angulaire de déformation de l'arbre.

Page
39
Par différence des équations (4-16) et (4-17) nous obtenons l'équation différentielle
homogène du second ordre en (θ)
. − . =−

.( + )− . = .

. +( + ). =− ⇒ + = .

− −
⇒ + . =

⇒ +K.θ + = ⇒
.
+K.θ =
.

( + )
+ . =
.

En posant: k =
.

( )
D’ou, ω= . pulsation propre du système
.

De l'équation (4-16) et (4-17), on peut transformer sous la forme matricielle:

−
. + . =
−
[M]. [ θ"] + [K]. [θ] = [0]
Où, [ ]- matrice des inerties du système.
[ ]- matrice des rigidités
Supposons que la vibration en torsion du système est libre et harmonique, c'est pourquoi la
solution peut être formée sous la forme matricieJ1e suivante
[ ] = [ ]. sin … … … … … . . (4 − 20)
[ ]- La matrice d'amplitude de l'angle de deux disques.

La période des oscillations du système est : T=

Page
40
Et la fréquence des oscillations par seconde sera: f = =

Deuxième Méthode:

Si le système subit une vibration en torsion, il existe sur le long de l'arbre une
section intermédiaire (nn) qui reste fixe pendant la vibration en torsion. Cette section
est appelé section neutre ou section nodale.
IL s'agit de déterminer sa position, c'est à dire calcul a et b.
Maintenant considérons le système de deux disques, un disque avec la longueur a et
l'autre disque avec la longueur b la période de deux systèmes est identique
(FIG:3-10)

FIG-3-10

On a :T = 2 , =

Puisque : T = T , on peut tirer:

= ⇒ =

.
Nous calculons enfin: = , k =
. .

La fréquence propre du système déterminée soit avec a, soit avec b:

ω = =d .
. .

Page
41
 .
Soit : = = .
. .

5. MODELISATION DES ACCOUPLEMENTS ELASTIQUES:

Quand les nœuds de vibrations se situent à leur niveau, les accouplements ont une grande
souplesse en torsion de là, il faut donner une grande importance à leur Modélisation
(FIG:3-11).

FIG-3-11-Modélisation d'un accouplement

K - Rigidité torsionnelle
I1 et I2 - Moments d'inertie des deux parties de l'accouplement
Remarque: on peut négliger I1 et I2 car ils sont très faibles

6. MODELISATION DES REDUCTEURS:

C'est un ensemble de roues dentées s'engrenant mutuellement. (Pignons, trains D’engrenages)

On a :

T = I .ω + I .ω = ω = ω I + .I

Soit : T = ω .I avec :I =I + .I
T - énergie cinétique du système équivalent
I - moment d'inertie du système équivalent
= ⁄ - rapport de réduction du couple
- vitesse angulaire de chaque roue.

Page
42
FIG-3-12-MODELISATION D’UN REDUCTEUR

Page
43
Introduction :

Un mode de vibration est un état privilège de vibration d'une structure où les forces de
déformation élastique sont compensées exactement par les forces d'inertie crées le
mouvement des masses.
Les efforts excitateurs ne sont alors équilibrés que par les forces d'amortissement
engendrées par la vitesse des mouvements vibratoires de la structure: on dit que cette structure
vibre à l'état de résonance.
Un mode vibratoire d'une machine est caractérisé par:
- Une fréquence de résonance (ou fréquence propre)
- Un déformée vibratoire modale (ou amplitude relative)
- Un Coefficient d'amortissement.
Les deux premières caractéristiques (fréquence propre et déformée) sont
indispensables à connaître pour la résolution qualitative d'un problème de vibration
de torsion.
La complexité de forme de la déformée modale (nombre de nœuds) augmente
avec le rang du monde vibratoire.

Page
44
1. HYPOTHESES :

- Le système obéit à la loi de HOOKE (les déformations sont proportionnelles aux

Efforts)
- Seules les déformations de torsion sont prises en compte
- Le système est conservatif
- L'axe de rotation est un axe de différentiel fixe.
- Les centres de gravité des disques appartiennent à l'axe
De rotation.
2. Équation du mouvement du disque:

Considérons (FIG-4-1) le modèle constitué par une ligne d'arbre comportant n

volant de moment d'inertie Il,12, ….In où sont supposées concentrées toutes les masses,
et relies entre eux par (0-1) portions d'arbre, OU sont concentrées les élasticités, et
présentant des rigidités torsionnelle K1 ,K2 ….kn-1 soient θ1,θ2, ….θn les angles de
rotation des disques au cours de la vibration.

FIG -4-1 –Notation des modèles

La mise en équations est faite à partir de la FIG (4-1)

Energie cinétique : 2E = ∑ I .

Energie potentielle :
= ( − ) + ( − ) + ( − )

Page
45
D’où l’équation de LAGRANGE: − =
On peut donc déduire le système différentiel du second
ordre a coefficients constants

. + ( − )=

. − . +( + ) − . = ……( − )

∀ ∊[ , − ]

[M][ ] + [K][ ] … … … … … … . ( − ) =0

. − . + . =

Avec : =

Le système différentiel ci dessus peut se mettre sous la forme matricielle suivante:

Avec :[ ] = -matrice carrée des inerties

[ ] – matrice carrée des rigidités

Et [ ]= ; [ ]=

Page
46
Le déplacement étant proportionnel aux efforts, la solution de cette équation
matricielle est de la forme:

….... (4-23)
Avec: [ ]= (i=1,…, n)
Donc, en reportant la solution (4-23) dans l'équation matricielle, nous obtenons:
-[ ][ ] + [ ][ ] =
−[ ] +[ ] [ ]=
Ou encore: −[ ] . [ ]. + [ ] .[ ] [ ] =

C'est un système à un n équations linéaires et homogènes en A1,A2, .... ,An qui n'est
résoluble (c'est à dire que le système ait une solution en Ai non nulle) ; d'après la
théorie des équations linéaires que si son déterminant est nul.
Si l'on désigne par [1] la matrice unité, les valeurs de sont données par:

Où encore, les valeurs de cherchées sont les valeurs propres de la matrice:

(En posant ( 2 = C)
[ ] = [ ] .[ ]

Page
47
3. METHODE DE RESOLUTION DU SYSTEME:

Il existe plusieurs méthodes de résolution parmi ces méthodes on a la méthode

Directe de résolution du système:
- Méthode de ULRICH – GIRAUDEAU
- Méthode de HOLZER
- Méthode de RUTISHAUSER
Nous avons choisi parmi toutes ces méthodes, la méthode de RUTISHAUSER car
elle présente l'avantage d'être précise et rapidement convergence, lorsque les va1curs
propres sont éloignées.
De plus la convergence peut être accélérée par une méthode de déflation.
3.1. ALGORITHME DE RUTISHAUSER:

a)- Propriétés des matrices tri diagonales: soit [ ] une matrice tri diagonale

Et (λ) le déterminant caractéristique de cette matrice

Si di j est le facteur de l'élément situé sur la ligne j la colonne; nous obtenons la

formule de récurrence:

Page
48
Les termes symétriques par rapport à la diagonale n'interviennent que par leur
produit. On peut donc les modifier à condition que leur produit reste constant. La
matrice [C·] définie ci après a donc les mêmes valeurs propres que [Cl mais par les
mêmes vecteurs propres:

b)- Algorithme: Il s'applique plus particulièrement aux matrices

GAUSS a démontré que l'on pouvait décomposer une matrice quasi-triangulaire en

un produit [R] x [ ] d'une façon unique.
Si l'une des matrices [ ] ou [S] possède des (1) sur sa diagonale principale sous
cette condition, on peut donc construire la suite matricielle suivante:

Les propriétés des éléments de cette suite sont:

 [ ][C.] ont les mêmes valeurs propres
 [ ][ ]··,[ ]sont des matrices quasi-triangulaires
 Lorsque n tend vers l'infini, [ ] tend vers une matrice triangulaire inférieure.
Les valeurs propres apparaissent donc sur la diagonale principale dans l'ordre
décroissant à partir du haut.

Page
49
3.2. APPLICATION DE L'ALGORITHME DE RUTISHAUSER:

a)- Détermination des matrices [R], [S]:

Nous obtenons les relations:

Les éléments blctf31 sont obtenu à partir des éléments C1eto.1 de la matrice [C·] par
les relations suivantes:

b)- Calcul des éléments de la matrice [ ]= [S] x [R]

Les éléments , De la matrice sont donnés par les relations:

, = + ∀ i ∊[ , − ]

, = ∀ i ∊[ , − ]
Page
50
 , , = ∀ i ∊[ , − ]

, = .

> +2
, =
> +2

La matrice [ ] obtenue est un tri diagonale possédant des (1) sous la diagonale
principale comme la matrice [ ° ]

c)- Relation entre les éléments de [ ] ET [ ° ]:

Elles sont obtenues à partir des relations précédentes que nous résumons:

C = C° C = +

b = C° - et = .

3.3. TEST DE CONVERGENCE QUADRATIQUE:

a)- Rapidité de la convergence:

- La convergence est d'autant plus rapide que la différence entre les valeurs propres
est grande. Pour les racines rapprochées, la convergence est médiocre.
- La méthode de RUTISHAUSER accélérée à convergence quadratique permet
d'obtenir, dans cc cas ct pour les matrices tri diagonales, des temps de calcul très
faibles avec une précision remarquable.
Page
51
- Lorsqu'on connaît une valeur propre; il est possible de se ramener à une matrice
d'ordre (n-l) qui contient les autres racines, le calcul est fait en même temps que la
déflation (diminution de la taille de la matrice) tout en conservant la forme tri
diagonale.
b)- Algorithme de la méthode de RUTHISHAUSER accélères:

La convergence quadratique est obtenue comme suit:

- Après chaque itération, on retranche la valeur (l'itération précédente) à tous les
termes diagonaux Ci j
- Lorsque Cn est nul, ou encore lorsque an est nul, la valeur propre VPn est la somme
des Cn retranchés, le test d'arrêt portera sur la valeur de .
- On recommence le processus en supprimant la dernière ligne et la dernière colonne
de la matrice [Cl, et en ajoutant VPn à tous les termes diagonaux de la matrice
d'ordre (n-l) on obtient la valeur propre VPn-1
- Le calcul s'arrête lorsque la valeur propre VP2 est calculée, VPI est alors égale à

Les valeurs propres ne sont pas obtenues par ordre croissant, elles apparaissent dans
un ordre quelconque

3.4. DEFORMEES MODALES:

Les valeurs propres de la matrice aux pulsations étant connues. IL est possible de
déterminer les vecteurs propres correspondants, ces derniers caractérisent les
déformées modales (ou les amplitudes relatives des disques) des divers modes de
vibration.
La matrice [Cl à utiliser, pour le calcul des composantes des vecteurs propres, est
celle défini par le produit : [M] [K]

Page
52
Les composantes des vecteurs propres étant définies à un coefficient prés, il est donc
nécessaire de fixer une composante, par exemple l'amplitude du disque n° l sera
prise égale 1 rad.
Les composantes des vecteurs propres sont obtenues à partir de l'équation:

[ ]− [ ] ( )=(0)

Avec: - Pulsation propre N ° i, [1] matrice unité

(Xi) - Matrice colonne des composantes du vecteur propre associé à la
Valeur propre
(0) - Matrice colonne composée d'éléments tous nuls Nous obtenons ainsi les n
relations suivantes:

() ∀ ∊[ , − ] - . , + − , − . , =

(n) - . , + − , =

En posant: , = , nous obtenons:

, = − . .

, = , − , +( − . ) ,

Page
53
∀ j ∊[ , − ]

La relation (n) inutilisée est automatiquement vérifiée. Elle permet de contrôler la

précision des calculs numériques.

4. VIBRATIONS DE TORSION SE EN COMPTE LA MASSE DE L'ARBRE) :

Pour illustrer ce cas, nous considérons le cas général de la torsion d'un élément
de poutre.
En appliquant les théorèmes généraux de la dynamique, on obtient l'équation aux
dérivées partielles qui régie ce mouvement. Cette étude va nous permettre de
déterminer les matrices de transfert de notre système afin d'aboutir- à l'équation de
fréquence recherchée.
Nous nous proposons donc, de donner les grandes lignes de cette méthode de calcul.

4.1. MOUVEMENT DE TORSION :

a)- Equation du mouvement:

Le mouvement de torsion d'une poutre est défini par les grandeurs suivantes:
M - Couple de torsion
- Angle de torsion
μ - Moment d'inertie de masse de l'inertie de l'arbre par rapport à l'axe Z par
Unité de longueur.
J - Moment d'inertie polaire
G - Module de COULOMB ou de cisaillement
Schématiquement, on représente un élément de poutre longueur dz soumis à la torsion
comme suit:

Page
54
Si l'on suppose qu'aucun couple extérieur ne sollicite cet élément. L’application du
théorème du moment cinétique:

∂M
.d . =M+ −M
∂Z

Soit : . = …………….( − )

Par ailleurs, le couple M et l'angle sont lies pat ta relation:

= . . ………………( − )

Cette dernière résulte de la formule de résistance des matériaux où le déplacement est

donné par:

. .
=∫ = …………….( − )
.

Avec: MZI - moment dû à un couple unitaire.

En reportant (4-28) dans (4-27) on obtient l'équation aux dérivées partielles qui régit le
mouvement de torsion qui s'écrit

. = . . ……….( − )

Dans le cas d'une section constante on a:

. = . . ……….( − )

D'autre part, pour une poutre de section droite circulaire I et J sont liés par la relation:

= . …………….( − )
Page
55
Cette dernière est obtenue en supposant que dans le mouvement de torsion, les sections
droites restent des sections droites.
L’équation du mouvement devient alors, pour une poutre de section circulaire:

Ou encore : = .

………… (4-33)
= .

Avec a= - célérité de l'onde de torsion

Remarque: l'utilisation de la relation = . Peut entraîner des erreurs

Importantes dans le cas d'Une section quelconque.
b)- Résolution mathématique de l'équation de torsion:

Pour résoudre J'équation (4-31) nous recherchons des solutions qui sont des produits de
fonctions de la forme:
( , ) = ( ). ( )
L'équation de torsion s'écrit alors:

μ. . = . . .

La séparation de la fonction de la variable d'espace de celle de ta variable du temps

Conduit a:

L'égalité entre ces deux fonctions est valable quelque soient Z et t ce qui entraîne que la
valeur commune des deux termes soit une constante que l'on prendra négative afin que les
solutions restent "bornées" dans te temps.
La constante est prise égale à (- ) ; il s'ensuit que:

Page
56
Et :

On forme alors le système suivant:

+ . =
.
+ . =
.

Dont les solutions sont données par:

f(t)= [Link] ( t) +[Link] ( t)

( )= . . + . . .
. .

À une valeur de ( ) est associée la solution:

(Z,t) = ([Link] t + [Link] t).([Link] Z + [Link] Z)

.
Avec : =
.

Remarque : Dans le cas d'une poutre de section circulaire devient égal à:

.
= =

Page
57
4.2. MATRICES DE TRANSFERT DUES A LA VIBRATION DE TORSION:

Nous déterminons la matrice de transfert d'un élément de poutre soumis à la torsion, en

établissant la relation de transfert liant les points j et j+ 1 de l'élément suivant:
Au point j c'est à dire en Z=O s'exerce un couple de torsion M j
Au point j+ 1 c'est à dire en Z=L s'exerce un couple de torsion M j+l

En considérons que la poutre présente une section droite circulaire l'équation qui
régit le mouvement de torsion s'écrit alors:

. ( )
+ . . ( )= ……….( − )

Avec :
G - module de cisaillement
–masse volumique

.
En posant: =

L’équation (4-34) devient:

. ( )
+ . ( ) ………..( − )

La solution de l'équation (4-35) est formée sous la forme:

(z)=[Link] Z + [Link] Z …………. (4-36)

Les constantes A, B déterminées par les conditions aux limites de l'arbre.

Page
58
a)- Etablir la matrice de transfert:

En dérivant (4-36) par rapport à Z, il vient que:

ϕ (z) = − . + . ……..( − )

En utilisant la rc1ation (4-37). L'expression de (4-28) devient:

= . . ( )=− . . . + . . . … … … … … … … ( − 8)

L'expression (4-36) et (4-38) s'écrivent sous la forme matricielle:

ϕ A
= .
M − . . . . B

L'application de cette écriture aux extrémités du tronçon nous donne: - au point j

correspondant à Z=O, la relation suivante:
ϕ 1 0 A
= . … … … … . . . (4 − 39)
M 0 G. Jβ

- au point j+ 1 correspondant à Z=L, la relation suivante:

ϕ A
= . … . . (4 − 40)
M − . . . . B

De (4-40), on peut calculer:

…………………….…. (4-41)

En substituant (4-41) dans (4-40), on trouve:

Page
59
 1 0
ϕ 1 . ϕ
= . 1
M − . . . .
G. Jβ M

Après l'opération de produit de deux matrices, on a:

= . . . ……( − )
− .
.. .
La formule (4-42) est la relation de transfert de l'arbre de point j+ 1 à j, on peut
désigner la matrice:
⁄ .
[ ]=
− . .

[ ] - est la matrice de transfert d'un élément de poutre soumis à la

Torsion.
Remarque: Dans le cas où l'on néglige la masse du tronçon c'est à dire lorsque
( =0) la matrice de transfert notée:
° ⁄ .
=

b)- Matrice de transfert d'un disque:

Pour un disque de moment d'inertie 1 et de masse M, les conditions de part ct d'autre

du disque sont:

- égalité des déplacements: =

- à discontinuité des couples: = −

Page
60
La relation de transfert se traduit par la relation matricielle suivante:

Où encore:
[D] - matrice de transfert du disque.

4.3. APPLICATION DE LA METHODE DES MATRICES DE TRANSFERT

A LA DETERMINATION DES FREQUENCES PROPRES DE TORSION
DU SYSTEME (arbre + disque):

Dans l'étude des vibrations de torsion, on suppose un disque et un tronçon

d'arbre fixé, comme indiqué sur les schémas ci-dessous:
On subdivise notre système en deux éléments, l'un comportant un disque de masse
M ct de moment d'inertie I, l'autre représentant un tronçon, de tronçon constant d et
de longueur L (figure ci-contre)

On établit alors la relation de transfert de chaque élément puis la matrice de

transfert globale résultant du produit des matrices.
En explicitant les conditions aux limites, on aboutit à l'équation de fréquence

Page
61
a)- Relation de transfert du premier élément:

Cet élément est constitué de la poutre la longueur L et de masse M. La relation de

transfert entre les points:
[ ] =[ ] … … … … … … … .. … … … … ( − )

⁄ .
= .
− . .

b)- Relation de transfert du second élément:

Cet élément est représente par le disque, de moment d'inertie 1, La relation de

transfert reliant les points (3) et (2) s'écrit, connaissant la matrice [Dt]
Ainsi :

Où encore : = .
− .

c)- Relation de transfert globale du système:

La relation de transfert du système, lie les points (3) et (1) par l'intermédiaire de la
matrice globale de la relation (4-43)

On a: [ ] = [ ]. [ ]

Que l'on reporte dans (4-44) :

[ ] =[ ]. [ ] = [ ]. [ ]. [ ] ……….( − )

On pose : D =[ ]. [ ]

En explicitant cette dernière relation, on obtient:

1 sin βL⁄G. Jβ
[D ] =
− I cosβL − βG. J. sinβL cosβL − . I. sin βL⁄G. Jβ

La relation de transfert (4-45) devient:

Page
62
 ⁄ .
= − − . . −
.. . ……..( − )
.

d)- Equation de fréquence:

L'équation de fréquences est obtenue en introduisant les conditions aux

Limites qu'impose notre système, ainsi;
En Z=O correspondant à l'encastrement (point 1)

≠ =
En Z=L correspondant à l'extrémité libre (point 3)

M3=0 ≠

En explicitant ces conditions dans la relation (4-46), il vient que:

=[ ].

En revenant à l'écriture algébrique, on forme le système suivant:

+ ⁄ .

..
= ( − )
.
Comme est différent de zéro, cette dernière équation n'est valable que si:

. .
− = …………….……..( − )
.

Cette relation est l'équation de fréquences à la torsion du système soumis à la torsion (lorsque
l'on prend en considération la masse de l'arbre)
Donc; d'après la relation (4-47), on a:

. = (L’équation que l'on résous graphiquement)

En posant L=X
, on obtient l'équation:

tg x = OU =

Page
63
 = ⇒ = ⇒ = .

= ⇒ = ⇒ = .

Donc : = .

Remarque:

Lorsqu’on néglige la masse de l'arbre, la relation (4-46) précédemment établie, reste

valable, à la seule différence, que la matrice de transfert de la poutre (premier élément) est
notée [D° ]

Il vient que alors: [r] = [D ]. [D° ]. [r]

Ou encore en explicitant les conditions aux limites:

⁄ .
= . .
− −
.

= .
On forme le système : .
= − .

.
La dernière équation n'est vérifiée que si: =0

Page
64
La fréquence propre est alors immédiate et donnée par:

= …….. (4-48)
.

Avec: G- module de cisaillement

I- moment d'inertie du disque

J- moment d'inertie polaire de

L- longueur de l’arbre.

Page
65
5- partie de calcul :

5.1. ORGANISATION DES CALCULS:

Pour notre étude de vibration de torsion, en se basant sur les données

théoriques précédentes, et avec la méthode de RUTHIS, nous avons réussi à réaliser
un programme qui permet de calculer les pulsations propres, ainsi que les déformes
modales de n'importe quelle ligne d'arbre non ramifiée.

5.2. Description du programme RUTISHAUSER :

a)-lecture des données :

Les noms des variables utilisées sont:

N : nombre de disque des schémas équivalents
EPSI : valeur permettant de contrôler la convergence de la suite matricielle de
RUTISHAUSER; valeur courante:
INERT (I): moment d'inertie par rapport à son axe, du disque "i" (N valeurs à lire;
en N.m.s)
RAID (I): rigidité torsionnelle du tronçon séparant les disques i et i+ 1 (N-l) valeur à
lire en (Nm/rd)

b)- calcul des éléments (C) et ( ) de la matrice [ ]:

Les éléments (C) de la diagonale principale de la matrice [ ] sont rangés dans le

tableau C(I) puis dans D(I), afin d'être conservés jusqu'au test de convergence.
Les éléments ( ) situés au dessus de la diagonale principale de la matrice [ ] sont
rangés dans le tableau ALPHA(I)
c)-calcul des valeurs propres, test de convergence, test d’arrêt :

Les différentes séquences de calcul apparaissent dans l'ordre suivant:

- Calcul des éléments (C) et ( ) de la nouvelle matrice [C] de la suite matricielle de
RUTISHAUSER
- Calcul de la variation relative des éléments [ ] entre cette itération et la

Page
66
précédente
- Rangement dans le tableau [ ] des éléments (C) de cette itération, en vue du
calcul de la variation relative des éléments (C) au cours de la prochaine itération.
- Test de convergence sur la variation relative des éléments (C)
- Test de sécurité sur le nombre d'itération (nous
Considérons que 100 itérations suffisent généralement
Noms de variables utilisées:

NB: nombre d'itérations

A: élément (b) de la matrice [R]
B: élément (B) de la matrice [S]

=∑ − : Éléments CI à l'itération 1

Q=s1/s2 : Pour test d'arrêt ou de convergence

d)- calcul des pulsations propres :

Les valeurs propres cherchées sont les éléments du tableau C(I) lorsque le test de
convergence est satisfait.
Les valeurs propres apparaissent par ordre décroissant; C(l) est la plus grande,) est la
plus petite.
Les systèmes à (N) disques étudies sont libres à leurs deux extrémités, il n'existe
donc que (N-l) pulsations différentes de zéro.
Noms de variables utilisées:

OM(J) [rd/s] - Pulsation propre qui correspond au mode (N-J)

TM(J) [tr/min] - Pulsation propre qui correspond au mode (N-J)
e)- Calcul des déformées modales :

Le nom de variable utilisé est:

X (I, J) : amplitude relative du disque au mode J

Page
67
Page
68
Page
69
Page
70
Page
71
Page
72
5.4. APPLICATION DE CALCUL :

5.4.1. Description de l'installation : L’étude qui va suivre portera sur une

installation motrice arbre hélice pour un avion léger.
*Le moteur:
Utilisé est de puissance moyenne, il comporte 6 cylindres en ligne attachés à un
vilebrequin et qui porte à son extrémité un volant de grande inertie utilisé pour
homogénéiser le mouvement des pistons.
Le passage du moteur à l'arbre porte hélice se fait par l'intermédiaire qui permet le
changement de la vitesse de rotation du vilebrequin à l'hélice.
*L’arbre port hélice : est un tronçon plein de l, 5 m de longueur construit
Généralement en acier A-32, de rigidité G =8.10-10 N/m2 , attaché à la 2eme roue
Du réducteur et porte hélice à son extrémité.
* un accouplement : grand vitesse et disposé entre l'extrémité de l'arbre et le
moyen de hélice pour prévenir un fonctionnement à vide, et pour concentrer les
nœuds de vibrations.
*hélice : est un disque de grande inertie; elle peut être par réglable au sol ou en vol
suivant la mission de l'avion.
Le schéma de principe de l'installation est le suivant:

Page
73
Donc, il s'agit de l'installation d'un moteur 6 cylindres en V à cycles à 4 temps à
allumage commandé.
5.4.2. Modélisation de L'installation:
Réellement on a une ligne d'arbre essentiellement composés de coude, d'arbre,
manivelles, arbres épaules, accouplement, rotor, réducteur, hélice, ..., etc. seront
modélisés par des disques ayant différents inerties (1) relies entre eux par des
raideurs (K).
L'ensemble mécanique est modélisé comme suit:
a)- Modélisation du moteur:
Le vilebrequin avec ses manivelles, bielles, pistons ct volants, est une pièce trop
compliquée pour que l'on puisse déterminer exactement les fréquences naturelles en
torsion.
Ainsi, devrons-nous pour les déterminer, simplifier le système, en remplaçant les
éléments du mécanisme bielles-manivelles par des disques équivalents, et le
vilebrequin par un arbre rectiligne ayant la même flexibilité à la torsion.
Le moteur tourne à la vitesse N1, à la sortie du moteur on rajoute un volant dans le
but d'emmagasiner ct d'absorber l'énergie cinétique.

Page
74
b)- Modélisation de liaison:

On entend par liaison; les accouplements qu'il faut leur donner une grande
importance puisque c'est à leur niveau que se situent les nœuds de vibrations. Le
choix d'un bon accouplement donne une meilleure tenue de la ligne d'arbre.
Ils sont généralement modélisés par deux disques séparés par une rigidité, les AGV
et APV se sont des accouplements a grande et petite vitesse respectivement, leur
intérêt est de transmettre l'énergie admissible à l'état couplé.
.
c)- Modélisation de l'hélice:
Par l'intermédiaire de l'APV, en vue du régime de l'hélice nominale, le moteur par
l'intermédiaire de l'arbre transmettre une partie de son énergie cinétique à travers le
vilebrequin permettant à l'hélice d'avoir une bonne tenue aérodynamique, vibratoire,
et un bon rendement; elle est modélisée par une inertie In
NOTA: Le réducteur et le volant seront rassemblés et modélisés par un seul disque
"14 "

Page
75
5.4.3. Caractéristiques inertielles et torsionnelle du modèle:

- Les moments d'inerties (1), les rigidités à la torsion (K) sont indiquées dans le tableau suivant:

Tableau des moments d'inerties et des rigidités au modèle du moteur V6

.
N° de disque Désignation Moment N° de Tronçon Rigidité
d'inerties 1 en torsionnelle K
K/m2 en N.m rd
1 Moteur 0,0184 1 4,88 .105
2 0,0184
3 0,0184
4 Volant d'inertie 0,2000 2 4,88 .105
+ réducteur
5 Hélice 0,5000 3 8,10 .105
4 0,20 105

5.4.4. Détermination des pulsations propres et déformées modales :

A laide du programme de RUTHISHAUSER, on obtient directement les pulsations

propres et déformées modales des différents modes.
Les résultats obtenus sont consignés dans les tableaux.

Page
76
Page
77
I/PARTIE DONNEES
.................
r------------------T
I MATRICE INERTIES I
L------------------J
T--------------T------------------T
I DISQUE N°: I INERTIES(kgm2) I
I--------------L------------------J
I 1 I .0184 I
L---------------------------------L
I 2 I .0184 I
L---------------------------------L
I 3 I .0184 I
L---------------------------------L
I 4 I .2000 I
L---------------------------------L
I 5 I .5000 I
L---------------------------------L
r------------------T
I MATRICE RAIDEURS I
L------------------J
T-------------T-----------------T
I TROCONS N°:I RAIDEURS (Nm/rd) I
L-------------L-----------------L
I 1 I 4.88E+5 I
L-------------------------------L
I 2 I 4.88E+5 I
L-------------------------------L
I 3 I 8.10E+5 I
L-------------------------------L
I 4 I .20E+5 I
L-------------------------------L
II/PARTIE CALCUL ET DETERMINATION
---------------------------------
T----------------------------------------T
DETERMINATION DES ELEMENTS(C) ET(ALPHA) I
L------L----------------L----------------I
I (I) I C(I) I alpha(I) I
L------L----------------L----------------L
I 1 I .2652174E+8 I .000000E+14 I
L------L----------------L----------------L
I 2 I .0000000E+8 I 7.034026E+14 I
L------L----------------L----------------L
I 3 I .0000000E+8 I 7.034026E+14 I
L------L----------------L----------------L
I 4 I .0000000E+8 I 1.782881E+14 I
L------L----------------L----------------L
I 5 I .0097600E+8 I .000040E+14 I
L------L----------------L----------------L

T------------------------------------------I
DETERMINATION DES PULSATIONS PROPRES EN I
L----------L-------------L-----------------L
I MODE N°: I OMEGA(rad/s)I VITESSE(tr/mIN I
L----------L-------------L-----------------L
I 1 I 343.6625 I 3283.399 I
L----------L-------------L-----------------L
I 2 I 2856.2206 L 27288.731 I
L----------L-------------L-----------------L
I 3 I 7154.2304 I 68352.514 I
L----------L-------------L-----------------L
I 4 I 9738.5615 I 93043.577 I
L----------L-------------L-----------------L

Page 80
T-----------------------------------------------------------T
I CALCUL DES DEFORMES MODALES DES QUATRES(4) PREMIER MODES I
L-----------------------------------------------------------L
T--------------------------------------------------T
I MODE N°: 1 SA PULSATION PROPRES w=343.6625rad/s I
T--------------T-----------------------------------T
I DISQUE N° I AMPLITUDE RELATIVE I
L--------------L-----------------------------------J
I 1 I 1.0000000 I
L--------------------------------------------------J
I 2 I 0.9955469 I
L--------------------------------------------------J
I 3 I 0.9786595 I
L--------------------------------------------------J
I 4 I 0.9786595 I
L--------------------------------------------------J
I 5 I -0.5012140 I
L--------------------------------------------------J
T--------------------------------------------------T
I MODE N°:2 SAPULSATION PROPRES w=2856.2206 rad/s I
T--------------T-----------------------------------T
I DISQUE N° I AMPLITUDE RELATIVE I
L--------------L-----------------------------------J
I 1 I 1.0000000 I
L--------------------------------------------------J
I 2 I 0.6924034 I
L--------------------------------------------------J
I 3 I 0.1718258 I
L--------------------------------------------------J
I 4 I -0.1726484 I
L--------------------------------------------------J
I 5 I 0.0008773 I
L--------------------------------------------------J
T--------------------------------------------------T
I MODE N°:3 SA PULSATION PROPRES w=7154.230rad/s I
T--------------T-----------------------------------T
I DISQUE N°: I AMPLITUDE RELATIVE I
L--------------L-----------------------------------J
I 1 I 1.0000000 I
L--------------------------------------------------J
I 2 I -0.9298515 I
L--------------------------------------------------J
I 3 I -0.652277 I
L--------------------------------------------------J
I 4 I 0.0917269 I
L--------------------------------------------------J
I 5 I -0.0002019 I
L--------------------------------------------------J
T--------------------------------------------------T
I MODE N°: 4 SA PULSATION PROPRES w=9728.561rad/s I
T--------------T-----------------------------------T
I DISQUE N° I AMPLITUDE RELATIVE I
L--------------L-----------------------------------J
I 1 I 1.0000000 I
L--------------------------------------------------J
I 2 I -2.5759190 I
L--------------------------------------------------J
I 3 I 3.0594390 I
L--------------------------------------------------J
I 4 I -0.1366283 I
L--------------------------------------------------J
I 5 I 0.0003650 I
L--------------------------------------------------J

Page 81
 CONCLUSION

Le présent mémoire traite les vibrations de machines et vise comme objectifs:

- L'analyse des éléments constituant la machine
- La modélisation dynamique
- La détermination des équations mathématiques
- Utilisation des méthodes de calculs des coefficients, des fréquences propres et
amplitudes
- Un programme d'analyse numérique.
La détermination des fréquences propres de vibration de torsion de la ligne
d'arbre permet de prévenir les risques dangereux (fréquences critiques) que
peuvent causer ces vibrations du point de vue de la résistance des matériaux et les
préjudices qu'elles peuvent causer au bon fonctionnement de l'appareil.
Nous entendons par fréquences critiques (vitesses dangereuses) la fréquence
pour la quelle, la vitesse de rotation est égale à la fréquence propre du système
(résonance).
Ainsi, grâce à la détermination des fréquences propres, on peut modifier ces
fréquences de telle sorte que la réponse des efforts que transmet la machine
(moteur, rotor) au fuselage soient minimaux.
Toute les techniques de calcul sont facilement introduites dans les
ordinateurs actuels qui permettent non seulement de résoudre rapidement les
systèmes d'équations les plus complexes mais surtout d'obtenir les solutions
exactes en limitant le nombre des hypothèses simplificatrices au minimum.
A l'aide de programme présenté dans ce mémoire, on obtient directement les
pulsations propres et les déformées modales des différents modes.
Donc, nous avons constaté que le calcul des vibrations était indispensable à
la machine, afin d'amortir les efforts transmis à la cellule de l'appareil.

Page 82
 BIBLIOGRAPHIE :
-Hélicoptère théorie pratique [Link] FORTS ET [Link]

7éme édition

-Mécanique des vibrations MICHAEL DELPEDRO édition

1985

-Dynamique des machines alternatives R .OUZIAUX DUNOD

Paris

-Techniques de l’ingénieur

-mini projet 4émeannée ingénieur Etude des vibrations de

torsion d’une installation d’un avion léger composées moteur-
arbre-hélice

-Aide mémoire de résistance des matériaux [Link]

HUEE édition 1975

-mécanique générale .1erpartie [Link] CIACO

éditeur Louvain-lanéure 1987.

-cours 5éme année ingénieur module vibrations

 CHAPITRE I
INTRODUCTION SUR LES
VIBRATION DES ARBRES
 CHAPITRE II
Calcul statique des arbres
(torsion)
Généralités
 CHAPITRE III
MODELISATION DYNAMIQUE
DES ARBRES
 CHAPITRE V
RECHERCHE DES PULSATION
PROPRE ET DEFORMES
MODALES