CORINA OANCEA

ANALYSE DES DENTS D'ENGRENAGES DROITS

PAR LA MÉTHODE DES POTENTIELS COMPLEXES

Mémoire
présenté
à la Faculté des études supérieures
de l'Université Laval

ponr l'obtention
du grade de maitre ès sciences (M.Sc.)

Département de génie mécanique

FACULTÉ DES SCIENCES ET DE GÉNIE

Janvier 1997

O Corina Oancea, 1997

1*1 National Library
of Canada
Bibliothèquenationale
du Canada
Acquisitions and Acquisitions et
Bibliographie Services services bibliographiques
395 Wellington Street 395, nie Wellington
Ottawa ON K1A O N 4 OttawaON KtAONQ
Canada Canada

The author has granted a non- L'auteur a accordé une licence non
exclusive licence allowing the exclusive permettant a la
National Lïbrary of Canada to Bibliothèque nationale du Canada de
reproduce, loan, distribute or seii reproduire, prêter, distribuer ou
copies of this thesis in microform, vendre des copies de cette thèse sous
paper or electronic formats. la forme de microfiche/film, de
reproduction sur papier ou sur format
électronique.

The author retains ownership of the L'auteur consenre la propriété du

copyright in this thesis. Neither the droit d'auteur qui protège cette thèse.
thesis nor substantial extracts kom it Ni la thèse ni des extraits substantiels
may be printed or otherwise de celle-ci ne doivent être imprimés
reproduced without the author's ou autrement reproduits sans son
permission. autorisation-
 AVANT-PROPOS

Cette thèse est le résuitat d'un travail qui a pu naître et se dérouler grâce à la confiance
que mon directeur de thèse, le Professeur Alain Cardou a eue en moi. Fairnerais le remercier
pour son appui et ses conseils donnés tout le long de ce travaii, pour sa compréhension et sa
patience, et pour son support financier.

Je remercie également le Professeur Qaude Gosselin pour Faimable coilaboration que

nous avons eue dans le domaine des engrenages internes.

Je remercie toutes les personnes que j'ai rencontrées et avec lesquelles j'ai travaillé.
Je remercie particulièrement mes collègues de "trois-mille sept" pour leur accueil chaleureux,
pour toute l'amitié avec laquelle fis m'ont entourée, pour tous ies moments formidables que
nous avons partagés ensemble.

Je remercie Huguette et Michel Roy pour leur générosité et leur amitié.

Je remercie mon fils Tudor et mon époux Florin pour leur amour,pour le courage de
tenter l'aventure d'une toute nouvelle vie.

Je dédie cette thèse à mes chers parents Nusa et Ioan qui, même de loin, ont toujours
été A mes côtés avec leur soutien et leurs encouragements.
 TABLE DE MATIÈRES

m
UMÉ.............................................................................................. IZ
AVANT-PROPOS ....................................................................................
TABLE DE MATIÈRES .............................................................................IV
LETE DES SYMBOLES ...........................................................................Vm.
LISTE DES FIGURES ..,.......................................................................... .XI

CHAPITRE 1
LE C.~LCULDE LA RÉSISTANCEDES ENGRENAGESDROITS...................-1
1.1. Introduction ..................................................................................-1
1.2. Historique du calcul des contraintes dans Ies engrenages droits ........................1
1.3 . ~ ~ ~ l i c a tde
i olan héo one de l~asticité...................................................-6
1.4. Application de la théorie des fonctions d'une variable complexe ......................-7
1.S. Objectifs .......................................................................................11

CHAPITRE II
PE~ÉSENTATIONDE LA MÉTHoDE DES POTENTIELS COMPLEXES ............. -12
2.1. Variables complexes et fonctions de variables complexes .............................. -12
2.1 .1. Variables compkxes ................................................................. 12
2.1 .2. Fonction de variables complexes ..................................................-14
2.2. Transformation conforme................................................................ 15
2.2.1. Transformations ou Représentations ..............................................-15
2.2.2. Représentation conforme ............................................................16
2 .3 . La Méthode des Potentiels Complexes..................................................... 16
2.4. Équations de base dans la 'ïhéorie de lfiasticité plane .................................. 18
2.4.1. Équations de ~olossov.............................................................-18
. .
2.4.2. Condition limite ...................................................................1 9
2.5. Domaine semi-infii soumis à une force concentrée......................................20
2 .6 . Application de la Méthode des Potentiels Complexes ....................................21
2.6.1 . Expression du problème &lasticité après transformation conforme..........-23
2.6.2. TransfonnationdelaconditionIimite..............................................26
 2.6.3. Calcul des potentiek @((FIet Y(c) ............................................... 28
2.6.3.1. Calcul du potentieI @(c) .................................................-28
2.6.3.2. CalddupotentieI v(5) .................................................. 32
2.6.4. Transformation des équations de KoIossov ..................................... 35
2.6.4.1. Exphcitafion de @(O ...................................................... 36
2.6.4.2. Calcul des contraintes superficieiles.....................................-38
2.6.5. Modélisation ..........................................................................4
2.6.5.1. Type de dent à étudier ......................................................40
2.6 S.2. Ajustement des axes de coordomées.....................................41
2.6.6. Recherche de la transformation conforme par la programmation rion-
linéaire ...................................................................................42

CHAPITRE III
TRANSFORMATION CONFORME D'UN DOMADE SEMI-INFINI DANS
LE DEMIEMIPLAN. DANS LE CAS DES 3 DENTS .......................................... -44
3.1 . Translation de la dent centrale ............................................................... 44
3.2. Rotation d'un profil ............................. .......... . ........................*.....-47
3 .2. 1 . Position de l'axe de syrnéme du profil ........................................... 48
3.2.2. Position du profil .................................................................... -49
3.3. Transformation conforme .................................................................... 50
3.3.1. Transformation de la dent centrale ..................................................51
3.3.2. Transformation de la dem (1). dent translatée 2 droite ...........................51
3.3.3. Transformation de la dent (2). dent translatée à gauche..........................55
3.4. Recherche des coefficients de la n'ansformation conforme ..............................60

CHAPITRE TV
C O N T l U N ï E S SUPERFICIELLES DANS LES DENTS D'ENGRENAGES

4.2. Explicitation du potentiel @(CI pour la transformation conforrne..................... -68

.
..
4.3 . Cdcui des connaintes superfÏcieUes ....................................................... -71
4.4. L'infhence des dents adjacentes sur la contrainte superficielle .........................74
4.5. Calcul du facteur géométrique.......... . .....................
 CHAPITRE V
APPLICATION DE LA A&THODE DES POTENTIELS COMPLEXES AUX
ENGRENAGES INTERNES................................................................... 78
5.1. Définition d'un engrenage interne .........................................................-78
5.2. Géneration d'une dent à développante de cercle ......................................... -79
5.2.1 . Dimensions de gén6ration ..........................................................-80
5.2.2. Gdnhtion du prof2 actif............................................................82
5.2.3. Génération du profil de raccord............................................. 8 3
5 .2.4 . Assemblage des profils .............................................................-86
5.2.5. Influence du type de couteau sur le profil de raccord ............................ 87
5.2.6. Comparaison du profil complet avec le profil d é t e d é par une autre
m6thode ................................................................................. 88
5 .3 . Transformation d'un engrenage externe dans un engrenage interne ................... -89
5.3.1. Travaux publits ..................................................................... 89
5.3.2. L'application de la Méthode des Potentiels Complexes.........................-90
5 -3.2.1 . Coefficients de la transformation conforme..............................93
5.3.2.2. Calcul des contraintes pour les engrenages internes .................... 97
5.3.2.3. Section critique .............................................................. 98
5.3 -2.4. Formules pratiques pour calculer les contraintes superficielles
dans les engrenages internes.......... . . . ...............................100
5.3.2.5. Facteur géométrique ....................................................... -103
CONCLUSIONS ..................................................................................... 105
BIBLIOGRAPHIE .................................................................................. -106
ANNEXE 1 ......................................................................................... *.lIO
CALCUL DES CONSTANTES E&ELLES ................................................ 111
ANNEXE 2 ........................................................................................... 116
~ W D DEE CAS
ENGRENAGE EXTERNEDROK .PAWMÈIRES DE
GÉNÉRATIoN .............................................................117
&I'UDE DE CAS
. SNGRENAGE EXTERNE DROU .PARAMÈTRJ~s DE
GÉNÉRATION ............................................................ - 117
INDEX ................................................................................................ -168
 LISTE DES SYMBOLES

v2 =iapiacien
a = addendum
A(x,y) = afnxe de z
ajk, Aj, Bj,Xj, Yj = constantes réelles
ak = facteur de concentration de contraintes (Niemiuin et GIaubitz)
$ = angle de la force W avec la direction x

c, ak, bk = coefficients de la transformation conforme

t I

C,, CL = courbes dans le plan (u,v)

CI, = courbes dans le plan z
d = translation sur l'axe y
6 = translation sur l'axe x
E l 1, E12, Dl 1,D12 = fonctions auxiliaires
Q, ~2 = angles des rotations
F = largeur de la dent
$ = angle de pression
a, , Y = potentiels complexes
F(x,y) = fonction d'Airy
f(z) = fonction de variable complexe
FI1, F12,Tl 1, T l 2 = fonctions auxiliaires
+c = angle de pression
$i = angle de pression correspondant en n'importe quel point du profil
Qinvc = angle de pression correspondant à Rmvc
y, = constantes complexes
h = hauteur de la dent

Ind (z) = index d'une courbe par rapport ii z

J = facteur géométrique
kg. kv,km = facteurs AGMA
kb, k,ks= facteurs de concentration de contraintes (AUison et Hearn)
km = facteur de concentration de moment (Weilauer et Seireg)
kt = facteur de concentration de contraintes (Tiioshenko et Baud)
5 p = param&tresde Lamé
m = module de la dent
mc = rapport de contact
MPC = Mkthode des Potentiels Complexes
ps = ffactur de forme (Niermm et Glaubitz)
N = nombre de dents
v = coefficient de Poisson
n = direction n o d e à 5 = cst
P = pas diamétral
= angle entre les directions t et x dans le système de coordonnées curvilignes liée à la dent
d'engrenage
R = rayon primitif
r = rayon standard de I'outil
Rg = rayon de tête
Rb = rayon de base
Rdc = rayon de l'arrondi au sommet du couteau
Re = rayon externe
& = rayon externe du couteau
Reca = rayon externe du couteau
Rf= rayon correspondant au point d'intersection du profil actif et du pronl de raccord de la
dent
Ri = rayon interne
ri = rayon correspondant en n'importe quel point du profil
Rinvc= rayon de courbure du sommet du couteau pour le cas arrondi non-nul
Rr= rayon de pied
Rra = rayon pour un nouveau profil de raccord
rt = rayon d'arrondi de l'outil
s = épaisseur de la dent au niveau de la section critique
GC = contrainte superficielle en compression
G C ~ o, c ~ = c contraintes dans la section critique pour le côté compression
at = contrainte superficielle en tension
a t ~ ba, t ~ =ccontraintes dans la section critique pour le côté tension
a ~ , =qcontraintes
T = epaisseur circulaire sur le cercle de primitif
t = coordonnée curviIigne li6e ik la surface du domaine
m. = position de la force sur le domaine
Ti = épaisseur circulaire sur n'importe quel autre cercIe
Tinvc = 6p"seur circulaire à Rivc
Tx,Ty= composaates de la force T selon les directioos x et y variant avec la position
curviligne s
TV = cisaillement
u(x.y) = partie réelle d'une fonction de variable complexe
U(z) = fonction holomorphe complexe
v(x,y) = partie imaginaire d'une fonction de variable complexe
W = force appliquée à un point de la fkontiè~
o(Q= fonction de nan~fomiation
wj = coefficients de poids
Wt = force tangentielle entre les engrenages
x = déport
X,Y = coordomées cart6siemes de la dent
6, q = cordonnées liées à 6
x, y = coordonnées liées à z
xo. y0 = position de force sur la dent
Xi= I'abscisse d'un point du proffi réel
-
xi = I'abscisse d'un point du prafil tracé par la transformation conforme
Xf,Yf= coordonnées cartésiennes correspondantes au Rf
Xk,Yk = constantes réelles
Xt, Yt = coordomkes cartésiennes du profil de raccord de la dent
ma,YTa = coordonnées cartésiennes du pronl de raccord
y = ordonnée du point d'application de la force
Yt = angIe entre le rayon et la tangente à la trochoïde
Z = nombre de dents
z = variable complexe
-
z = variable complexe conjuguée
Zc = nombre de dents du couteau
 LISTE DES FIGURES

Figure 2.1. Diagramme d'Argand ou plan z .................................-13

Figure 2.2. Représentation conforme..........................................16
Figure 2.3. Force agissant sur un segment AB d'un solide................. 19
Figure 2.4. Force appliquée à un point quelconque z,, de la
frontiere du domaine...............................................21
+ 4

Figure 2.5. Deux dents en contact transformées en demi-plan chargé

par une force W .................................................... 22
Figure 2.6. Force appliquée en un point quelconque r, de la
frontière d'un espace semi4nhi.................................23
Figure 2.7. Espace semi-infini ................................................. 25
Figure 2.8. Système de coordonnées curvilignes Lié à la dent
d'engrenage ............................... . . ..................-36
Figure 2.9. Dent d'engrenage à modéliser ...................................-40
Figure 2.10. TransIation de la dent ..............................................41
CHAPITRE In
Figure 3.1. TmIation de la dent .............................................-44
Figure 3.2. L'influence de Ia translation sur "x" sur Ies dents
d'engrenage droit ............................... . . . . .....A 6
Figure 3.3. Rotation d'un profil .................................. . 47
Figure 3.4. Position de Paxe de symétrie du profil ..........................49
Figure 3.5. Position d'un profil MPC après rotation .......................-50
Figure 3.6. Transfomtion conforme proposée pour trois dents .........-50
Figure 3.7. Dent d'engrenage droit tracée par MPC ........................-51
Figure 3.8. Translation sur l'axe "x" de la dent centrale ....................52
Figure 3.9. Translation sur i'axe "y" de la dent centrale ...................-53
Figure 3.10. Dent d'engrenage droit MPC. qui a subi des translations
et une rotation ......................................................54
Figure 3.11. Deux dents d'engrenage droit représentées par MPC
(translation B droite) ...............................................54
Figure 3.12. Dent d'engrenage droit MPC. qui a subi des translations
et une rotation .......................................
 Figure 3.13. Deux dents d'engrenage droit représentées par MPC
(trans1ation à gauche) ............................................. -56
Figure 3.14. Représentation de trois dents par MPC..........................57
Figure 3.15. Position de chaque dent dans la transformation pour
trois dents .........,.................................................58
Figure 3.16. Dinikentes transformations appliquées aux dents
d'engrenage droit ................................................. -59
Figure 3.17. Profil approché d'une dent standard (pour le cas des
......
trois dents) ............. ..................................61
CHAPITRE IV
Figure 4.1. Contraintes superficielles dans un engrenage externe ........-74
CWXTRE: v
Figure 5.1. Parties d'me dent d'engrenage interne............... . .......78
Figure 5.2. Comparaison d'un engrenage externe et d'un engrenage
interne ..............................................................-80
Figure 5.3. Géométrie d'une dent d'engrenage interne ......................80
Figure 5.4. Gén6ration du profil actif ........................................-82
Figure 5.5. Génération du profil de raccord .................................-84
Figure 5.6. Assemblage des profils ............. ......... ................87
Figure 5.7. L'influence du couteau sur le profd de raccord .................88
Figure 5.8. Profil complet d'une dent d'engrenage interne ................-89
Figure 5.9. La transition de l'engrenage interne vers l'engrenage
externe .............................................................. -89
Figure 5.10. L'influence du nombre de termes de la transformation
conforme sur le profii de la dent ................................ -90
Figure 5.1 1. Transformaton des domaines proposée par Miyachika et
Oda .................................................................. 92
Figure 5.12. Engrenage interne. profil réel et profil tracé par MPC .........96
Figure 5.13. Contrainte superficielle dans un engrenage interne ...........-98
Figure 5.14. Section critique dans tui engrenage externe....................-99
Figure 5.15. Section critique dans un engrenage interne ....................-99
. Figure 5.16. Effort nominal dans une dent d'engrenage interne............-100
 Aujourd'hui les engrenages occupent une place spéciale dans les systèmes
mécaniques. C'est la façon la plus économique pour transmettre de la puissance et un
mouvement de rotation dans des conditions uniformes. Comme les exigences sont vastes et
avec des diff~cultésvariées, les engrenages sont très complexes et d'une grande diversité.
L'importance de l'engrenage, comme élément mécanique nécessaire et idéal, est démontrée
par la vaste gamme qu'on trouve dans toutes les industries. Le développement des nouvelles
technologies, comme l'électronique, a remplacé quelques applications de l'engrenage, mais il
reste toujours un élément mécanique dont l'utilisation croît continuellement

L'histoire des engrenages commence dans les civilisations antiques, avec les roues de
friction. Cependant, c'est seulement dans la période de la Renaissance que les mathématiciens
ont commencé à appliquer les principes géométriques pour déterminer le meilleur profil de la
dent d'engrenage. La conception d'un engrenage comprend des calculs mathématiques,
I'aspect géométrique, la détérioration, les materiaux, la fabrication et la vérification. Parmi
tous ces paramétres, il est essentiel de connaître avec précision les contraintes se trouvant
dans la dent d'engrenage pour prévenir certains risques de rupture. Par conséquent, plusieurs
méthodes théoriques et expérimentales ont été développées, ii partir de la fin du XM-ème
siècle.

1.2. Historique du calcul des contraintes dans Ies en-=nazes droits

Dans sa thèse, Paré 1281 présente un historique complet des methodes théoriques et
exp6rimentaIes développées pour le calcul des dents d'engrenages en flexion. Cette
présentation est rt5surn6e ci-dessous. Le lecteur trouvera le détail des références mentionnées
dans [28].
 C'est Lewis, en 1893, qui a analysé les dentures droites sym6triques par la théorie des
poutres, en inscrivant dans celles-ci une parabole d'isocontraintes représentant une poutre
d'egde résistance. L'extr6mit6 de celle-ci se trouve à i'intersection de la ligne d'action avec
Saxe de sym6trie de la dent. Lewis propose que la section de contrainte superficielle
maximale, ou la section critique, se situe au point de tangence de la parabole avec le profil de
la dent. Suite à cette hypothèse, la contrainte en tension se calcule comme pour une poutre en
flexion, avec un encastrement ii la section critique, sollicitée par la composante dans la
direction normale à la ligne de symetrie de la dent de Ia force appliquée sur la dent.

Un autre moyen d'6tude des contraintes dans les engrenages est la photoélasticité. En
1926, Timoshenko et Baud mesurent les contraintes et ils o b t i e ~ e ndest valeurs deux fois
plus grandes que celles calculées par la formule de Lewis. En pensant que cette différence est
due B la variation rapide de la section à la base de la dent, ils proposent un facteur de
concentration de contraintes (kt), facteur qui varie avec la largeur de la section et le rayon de
courbure au pied de la dent.

Comme la composante radiale de la force sur la dent crée une contrainte de

compression, en 1938 Memtt ajoute un terme supplémentaire à la formule de Lewis.

L'Amencan Gear Manufacturers Association (AGMA) utilise pour le calcul des

contraintes en tension dans les engrenages la formule proposée par Dolan et Broghamer en
1942. Après des études photoélastiques similaires à celles de Timoshenko et Baud, ils
trouvent que le facteur de concentration de contraintes O,qui multiplie la formule de Lewis
modifiée (flexionlcompression), est aussi lié à la position de la force sur la dent. Suite à cela
ils introduisent une formule empirique pour calculer le facteur de concentration de connaintes
en fonction de la section critique, de la hauteur de la charge, du rayon de courbure au pied de
la dent et de coefficients mesurés par photoélasticit6 (ceux-ci varient avec l'angle de pression
de la dent). Les résultats démontrent que la section critique est un peu moins élevée que celle
estimée par la parabole de Lewis. Leur étude a été réalisée pour des dents ayant des angles de
pression égaux B 14 1/20 et 200.

Toujours après une recherche photoélastique, en 1948, Heywood représente une dent
comme un trapèze équivalent et B partir de celui-ci développe une formule pour le calcul de la
contrainte maximale. Sa formule comprend en plus du facteur de concentration de contraintes
et du terme de contraintes en flexion, un facteur de correction de l'effet de la charge.
 En 1950, Niemann et Glaubitz font eux aussi des tests photoélastiques et ils proposent
que la contrainte maximale soit égale aux contraintes de. flexion, de compression et de
cisaillement, multiplikes par un facteur de concentration de contraintes ak.La contrainte de
cisaillement est reliée aux deux autres par un facteur ps dépendant de la forme de la dent.
Pour cette formule la section critique est située sous le cercle de tête, à deux fois le module de
Ia dent.

En 1957, KeUey et Pedersen réalisent eux aussi un test photodIastique, similaire à

celui de Dolan et Broghamer, avec l'intention d'étendre cette formule à I'angle de pression de
250. Avec les données recueillies, ils ont développé une formule analogue à celle de
Heywood et ajoutent un terne supplémentaire tenant compte de l'effet de I'angle de la
sollicitation, relativement aux direction principales.

Jacobson compare, en 1958, les résultats recueillis avec les formules proposees par:
Memtt, Dolan et Broghamer, Heywood, Niemann et Glaubitz, Kelley et Pederson appliquées
aux engrenages à un angle de pression de 200, toujours grâce à la photoelasticité. Suite à son
travail il conseille, pour les engrenages à 200 d'angle de pression, ayant entre 10 et 40 dents,
l'utilisation de la formule de Lewis modifiée, multipliée par un facteur de concentration de
contraintes dépendant du rayon de courbure du sommet de l'outil de taillage. Il propose aussi
une construction très simple pour localiser la section critique; au point de tangence du profii
de la dent avec une tangente à 300 par rapport à l'axe de symétrie de celle-ci.

En 1960, Wellauer et Seireg utilisent la théorie des plaques pour calculer le moment
de flexion dans une dent d'engrenage. Pour calculer les contraintes, ils suggèrent une formule
qui est en fonction d u . facteur de concentration de moment (km). Ils valident les résultats de
cette théorie avec un test expérimental en utilisant des jauges de contraintes. Ils proposent
l'utilisation de leur methode lorsque le rapport hauteurnargeur de la dent est faible.

En 1973, Wilcox et Coleman développent, à l'aide des rbsultats obtenus par les
éléments finis, une nouvelle formule pour le calcul de la contrainte au pied de la dent. Elle
peut s'appliquer aux dents symétriques ainsi qu'asymétriques, mais elle n'est pas fiable
lorsque la sollicitation se situe dans la partie inférieure de la dent.

Une ttude similaire est réaliske en 1974 par Chabert, Dang Tran et Mathis pour les
dents symétriques. Après comparaison avec les normes ISO {International Standardization
Organisation) et AGMA, ils suggèrent l'emploi de la norme ISO, celle-ci étant plus facile à
utiliser pour une pdcision comparable.

En 1974, Winter et Hirt mesurent expérimentalement les contraintes au pied de la dent

avec les jauges de contraintes. Suite à sa thèse de doctorat, Hirt conseille l'utilisation de la
norme ISO pour le calcul des contraintes car celle-ci tient compte d'un facteur de sensibilité
d'entaille en fatigue.

Comme le montre ce court historique, une des méthodes couramment suivie par les
chercheurs a &té d'obtenir des données par différentes techniques et de formuler des équations
empiriques permettant une détermination facile des contraintes. On doit souligner que la
plupan des formules ont été développées pour le côté en tension de la dent, car c'est de ce
côté qu'apparaissent les fissures.

Dans la plupart des rravaw présentés. les résultats n'étaient pas satisfaisants pour les
engrenages à grand nombre de dents. En 1980, Allison et H m suggèrent que cela est dû à
l'augmentation de l'influence de la contrainte en cisaillement sur ces engrenages. Plus un
engrenage a un nombre élevé de dents, plus celles-ci sont petites et ressentent le cisaillement
relativement aux autres contraintes. Suite à ces remarques, ils réalisent une étude
photoélastique, tenant compte aussi de la force de fnction et des différentes positions que la
force occupe lors de I'engrènement. Ils remarquent que la force de fnction affecte surtout la
contrainte du côté en tension, en la faisant varier avec la valeur du coefficient de fiction. De
même, lorsque la force est située au bas de la dent, la contrainte en tension devient plus
tlevée que celle en compression. Suite à ces r6sultats, ils proposent une fomule où chaque
terne de contrainte (flexion, compression, cisaillement) est multiplié par un facteur qui lui est
propre (kb , kc ,ks ). Comme cette force a un angle constant de 230 avec la surface de la
dent, cette formule est difficilement applicable pour les dents normalisées.

En 1980, ComelI modifie la formule de Heywood pour l'adapter à un logiciel qui

calcule les engrenages à grand rapport de contact (High Contact Ratio Gear, mc >2). Cette
modification est nécessaire pour que celle-ci s'applique à cette forme de dent.

On utilise des normes pour le calcul usuel des engrenages. Pour le calcul des
contraintes en flexion, la nome AGMA 218.01 est fkequernment utilisée au Canada. La
méthode proposée est celle développée par Dolan et Broghamer avec des facteurs
supplémentaires tenant compte du couple transmis h).de la variation de la charge avec la
vitesse (kv), et de la distribution de la force sur la dent (km). Présentement, les travaux
relatifs à cette norme consistent à l'informatisation des calculs, et à l'application d'algorithmes
pour identifier la section critique.

Depuis 1980, l'orientation des travaux de recherche a surtout 6té dirigée vers
I'application des méthodes num6riques: &men& Finis, Équations Intégrales, Méthode des
Potentiels Complexes, Bandes Finies.

La Méthode des É~émentsFis est une methode facile à appliquer et elle permet
d'knidier et de modéliser I'engrenage de façon plus globale. En effet avec cette mkthode, il est
possible de tenir compte de plusieurs paramètres comme le nombre de dents sur l'engrenage,
l'épaisseur de la jante, etc. Parmi les chercheurs qui ont utilisé la Méthode des Éléments F i s
on cite Oda, Nagamura et Aoki (198 1). Ils ont ddtemiiné l'influence de l'épaisseur de la jante,
sur laquelle repose la dent, ainsi que des dents adjacentes, sur la dent sollicitée. Ami, Harada
et Aida, en 1981 analysent l'influence de i'6paisseur de la jante pour une couronne dentée
supportant 9 dents. En 1985, Guillot et Tordion réalisent une étude similaire dans les
couronnes et jantes minces d'engrenages aux axes parallèles. En plus, Chang, Huston et Coy,
en 1983, étudient l'influence de la variation du profil de raccord. Ils ont trouvé que ce rayon
n'affecte que Ia contrainte en surface de la dent En 1984, Chong, Katayama, Kubo et Yabe
étudient par cette méthode les engrenages à denture interne.

En 1983 Rubenchick propose l'utilisation de la Méthodes des Équations Intdgrales

pour les engrenages droits. Ses résultats, comparés avec ceux de Cardou et Tordion 171
obtenus avec la Méthode des Potentiels Complexes, présentent une variation de 5% sur la
contrainte maximale.

En 1984, Gakwaya, Cardou et Dhatt ont comparé cette méthode avec ceile des
ÉI6rnent.sFinis. Selon leurs résultats, pour la même précision et pour des engrenages droits, la
Méthode des Équations Intégrales est moins dispendieuse à utiliser.

D'autre travaux ont été réalisés avec cette methode par les russes, comme Yakovlev en
1984, et par les japonais: Tanaka en 1984 et Oda,Miyachika, Koide et Minine en 1985.
1.3. &plication de la Théorie de ll%sticité

Le calcul des contraintes des dents d'engrenages sous charge peut se faire par diverses
mtthodes, analytiques ou num&iques, plus ou moins lourdes et plus ou moins exactes.
Puisqu'il s'agit la d'un problème d'élasticité plane, la Théorie de l'hsticité permet en principe
d'obtenir la solution exacte au problème.

La Théorie de ~hasticitéporte sur l'étude de la réponse des corps élastiques sous

l'action des forces extérieures. Un corps est élastique s'il possède la propriété de reprendre sa
forme originale après que les forces qui produisent les déformations cessent.

La propriété de l'élasticité est caractkrisée rnathkmatiquement par des relations liant

les forces et les déformations. Parmi les relations, la loi de Hooke est d'une importance
fondamentale. Elle exprime que les tensions produites sont proportîomelles aux forces qui
les produisent.

Une loi identique a été découverte indépendamment par Mariotte [3 11 et elle a été
utilisée pour étudier les efforts dans les poutres. Il déduit qu'une poutre résiste à la fiexion car
une partie des frbres sont en tension d'autres sont en compression. Mais Mariotte ne s'occupe
pas de la forme prise par l'axe de la poutre.

Cette h d e est réalisée par Jacob Bernoulli, qui combine les équations d'équilibre
avec la loi de Hooke et obtient l'équation différentielle de l'elastica (c'est à dire, la courbe de
déformation de l'axe de la poum).

Euler déduit I'tquation de Jacob Bernoulli en utilisant le principe de minimisation de

Daniel Bernoulli. Le travail d'Euler dans ce domaine a étd axé principalement sur l'étude de
l'elczstica et non sur la distribution des e f f dans
~ la section de la poutre.

L'étude de la distribution des efforts dans une poutre en tension et en flexion est due à
Coulomb, qui a aussi été le premier & étudier la résistance du fil mince à la torsion.

Le premier essai pour obtenir les équations d'équilibre dans les solides élastiques est
dû à Navier, et la date du 14 mai 1821 (date de son mémoire) marque la date de naissance de
la Théorie de l'~lasticit6.II deduit trois équations différentielles pour les composantes de
deplacement à l'intérieur d'un solide élastique isotrope. La forme de ses équations est
correcte, mais elles contiennent une seule constante d'6Iasticité. Navier obtient l'équation
d'équilibre à la surface d'un solide avec l'aide du principe de Lagrange.

La formulation mathématique de la Théorie de 1~lasticitt5qui n'a pas changé jusqu'à

présent, appartient 2 Cauchy. il démontre que I'état de tension dans un point d'un corps
déformable est complètement détermine par un système de 9 fonctions. Quand le corps est en
équili'bre les 9 fonctions doivent satisfaire 3 équations différentieIIes et leur nombre est réduit
à 6 à cause de certaines relations de symktne. L'état de déformation est aussi d6tefmin6 par 6
fonctions qui sont simplement liées aux composantes du vecteur de déplacement, quand le
déplacement est petit. Pour un corps élastique et des petites déformations, entre le système
des fonctions qui caractérise l'état de contrainte et celui qui caractérise l'état de déformation,
il y a une relation linéaire. Cette hypothèse représente une généralisation de la loi de Hooke.
Les équations de Cauchy sont semblables il celles de Navier mais elles ont deux constantes
d'elasticité. Cauchy a présenté ce résultat très important en 1822.

Depuis ce temps, la Théorie de l'Élasticité à été appliquée dans beaucoup de

problèmes d'ingénierie. La catégorie la plus ttudiée a 6té celle des probIèmes d'élasticité
plane, grâce à la grande variété des méthodes mathématiques disponibles.

1.4. A~rdicationde Ia théorie des fonctions d'une variabIe com~lexe

Les applications de la théorie des fonctions d'une variable complexe au problème plan
de la Theorie de l'Élasticité constituent un groupe de problèmes à part.

Les méthodes de KoIossov et Muskhelishvili [24] sont parmi les plus puissantes. Les
fonctions d'une variable complexe ont été introduites dans les problèmes d'élasticité plane en
1909 par Kolossov, et ont &té développées et utilisées par Muskhelishvili (1915). Ensemble
ils ont fondé une forte ecole russe. On peut mentionner les noms de Sokolnikoff [3 11, Savin
1291, Timoshenko [33], Ugodohikov et Kumetsov [35].

Dans ses études de 1961, Savin [29]analyse la distribution des contraintes autour de
trous de ciifferentes formes. Mentionnons en particulier le cas d'une plaque semi-infinie
perforée d'un trou circulaire, avec une tension appliquee parallèle à Ia frontière du semi-
infini.
 En 1960, Milne-Thomson [22] aborde les problèmes d'élasticité plane en utilisant la
théorie des fonctions d'une variable complexe. Il appIique la théorie de MuskhelishviIi, c'est-
M i r e la repr6sentation de la fonction biharmonique de deux variables la plus générale, au
moyen de deux fonctions arbitraires de la valeur complexe. Les contraintes et les
d6formations (en contrainte plane comme en déformation plane) se déduisent aisément de la
connaissance de ces deux fonctions arbitraires que I'auteur propose d'appeler "contraintes
complexes". Dans le cas d'un matériau isotrope, la solution des problèmes d'élasticité plane
consiste essentiellement à chercher les contraintes complexes. Dans le cas des corps qui ne
sont pas simplement connexes (anneaux par exemple), ü est nécessaire de supposer que les
contraintes complexes sont des fonctions uniformes. Dans ses travaux, Milne-Thomson
donne de très nombreux exemples d'application de la méthode. Ils sont relatifs aux corps
élastiques plans Iimités par des droites, des cercles ou même des courbes quelconques, grâce
à I'étude de la transformation des contraintes compIexes dms une représentation conforme.

Pour une étude plus approfondie des problèmes d'élasticité plane nous recommandons
les publications 114, 191.

En 1969, Sneddon et Lowengrub 1301, calculent la distribution des efforts au

voisinage d'une fissure dans le solide élastique (fissure connue sous le nom de "fissure de
Griffith") en utilisant les équations de Kolossov et Muskheiishvili [24].

En 1987, Hromadka et Lai [18] traitent la méthode des éléments finis par la MPC.
Leur méthodologie est la suivante: (1) définir un domaine polygonal simple et choisir une
fonction de variable complexe f (z) =u(x, y) + N(x,y) (avec z = x +iy), holomorphe dans le
domaine; (2) discrétiser la frontière du domaine dans rn déments de frontière; (3) définir en
chaque noeud une valeur nodale spécifique; (4) définir sur chaque élément nodal une fonction
linéaire d'essai, qui dépend des point nodaux extrêmes de 1161émentet des valeurs nodales
spécifiques en ces points; (5) définir une fonction linéaire globale d'essai; (6) définir une
fonction d'approximation comme l'intégrale de Cauchy pour la fonction linéaire globale
d'essai. Cette fonction est holomorphe daas le domaine et continue sur le contour et satisfait
les conditions iimites.

Grâce au thtorème de Riemann sur la transformation conforme qui etablit qu'il existe
une correspondance biunivoque entre un domaine quelconque et le disque unite, mais sans
expliciter la transformation conforme, plusieurs auteurs ont étudié les transformations
conformes des domaines compLigués dans des domaines plus simples à étudier. Parmi eux,
Wang [37Jmentionne la transformation d'une section avec cannelures dans le cercle unité.

En 1951, Hirano 1171 applique la Méthode des Potentiels Complexes (MPC) à la

résolution d'un problème d'élasticité plane, dans un domaine consistant en un plan semi-infini
sur lequel repose une courbe symétrique par rapport à l'axe des imaginaires. La méthode
repose sur la possibilité de passer par une transformation conforme adéquate du profil donné
(la "cloche" symétrique) au plan semi-infini.

En 1962, Ugodohikov et Ku~ietsovD5] appliquent la méthode de Muskhelishvili et

Kolossov pour calculer les efforts dans une dent d'engrenage, mais ils utilisent une
transformation conforme différente de celle de Hirano. Ils supposent que toute fonction qui
exprime le cercle unité dans la zone occupee du corps éiastique est transformée par une
fonction proposée par Ugodohikov.

Suite aux travaux de Hirano sur l'utilisation des potentiels complexes, en 1962, Aida
et Terauchi [l] reprennent l'idée et appliquent ses résultats aux profils d'engrenages droits,
mais la transformation conforme qu'ils utilisent est inadéquate pour représenter d'une manière
satisfaisante des pronls d'engrenages réels.

Avec la publication en 1973 de Baronet et Tordion 121, on a pu apprécier l'intérêt de la

methode pour le calcul des concentrations de contraintes au pied de dents normalisées, ou
encore, pour le calcul de leur flexibilité. En comparaison avec les méthodes proposées par
Dolan/Broghamer, Heywood, NiemandGlaubitz, KelleyRedersen, Coleman/Wilcox et
Chabert, la méthode de BaronetTordion donne de bons résultats aussi bien pour les
engrenages à grand nombre de dents (48 dents) qu'à faible nombre (14 dents). Cependant
leurs méthode n'est testée que pour de grands rayons de courbure au profd de raccord
(rt=0.38-m). Pourtant, le problème de trouver la transformation conforme d'un profil donné
n'est pas complètement résolu. Elle dépend de coefficients qui varient selon le type de dent à
étudier.

Ces travaux sont poursuivis par Cardou et Tordion [5,6,7],qui ajoutent en 1981 des
coefficients suppl6mentaûes à Ia transformation. Le problème d'application de cette méthode
est d'identifier les 2(0+1) coefficients (où n ddpend du nombre de termes dans la
transformation conforme). Une fois la transformation conforme réalisée, on emploie les
formules de la Théorie de l%lasticit6 appliquées au demi-plan, formules de Kolossov pi],
pour trouver les contraintes.

En 1985, dans un rapport substquent, Cardou et Tordion 181 compilent une série de
coefficients et de facteurs géométriques correspondant à plusieurs types d'engrenages
normalisés. Ils ont aussi étudie5 avec cette m4thode I'hfIuence du coefficient du frottement sur
les facteurs géométriques.

L'exactitude des rksultats fournis par la MPC a été mise en évidence par divers
travaux expérimentaux, notamment ceux de Gosselin, Gagnon, Cardou et Cloutier [I5]. Des
mesures de la déflexion d'une dent d'engrenage droit ont été effectuées au moyen d'un
interféromètre éIectronique à patrons de grandarités optiques. Les résultats expérimentaux
ont été comparés B ceux obtenus au moyen d'éléments n n i s de six types courants et six
densités de maillage ciiffikentes, et à ceux obtenus par la MPC. La comparaison des résultats
numériques et expérimentaux démontre que la Méthode des Potentiels Complexes calcule
précisément la déflexion d'une dent d'engrenage.

La MPC a été utilisée en 1993 par Nicoletto [25, 261 dans l'étude de la fissuration
dans b s engrenages, en se basant sur les travaux de Cardou et Tordion [6, 71. Dans ses
études, l'auteur analyse les contraintes dans la section critique et il fait une comparaison avec
la propagation de la fissure qu'il a déterminee expérimentalement. Le mode de fissuration est
caractérisé par un point d'origine situé sur la surface du profil de raccord, au point de
concentration maximale. Il trouve ainsi que la propagation de la fissure ressemble, au moins
au début, 2 la direction de cisaillement maximum. Nicoletto propose une méthode pour le
calcul du facteur d'intensité de contrainte basée sur l'utilisation d'une fonction de poids. La
fissure modifie la frontière du domaine qui modifie le contour d'intégration, en fonction de la
longueur de la fissure.

Suite aux travaux de Cardou et Tordion 16,71,Daniewicz 1121 présente en 1991 un

modé3le de prédiction analytique pour la courbe de fatigue pou. une dent d'engrenage. Le
modèle proposé comporte deux étapes: (a) l'origine de la fissure sur la surface du profil de
raccord au point de concentration maximale de contrainte, similaire B la methode adoptée par
I'AGMA; (b) la propagation de la fissure avec le calcul du facteur de concentration de
contrainte. Dans le calcul du facteur de concentration de contrainte, Daniewicz utifise les
composantes des contraintes déterminées par la MPC. Il modélise aussi analytiquement la
rigidité d'une dent fissurée; il utilise la MPC pour déterminer le déplacement du point O& la
force est appliquée, c'est 2 dire la flexibilité de la dent.

En 1991. Miyachika et O& [23] ont 6tudié l'effet de l'épaisseur de la jante d'un
engrenage interne sur les contraintes et la direction de propagation de la fissure sur la dent. Ils
traitent les engrenages internes par la même Méthode des Potentiels Complexes. La formule
de la transformation conforme est toujours celie proposée par Hirano, mais le nombre des
termes nécessaires est plus grand. Avec sept termes dans la transformation conforme, ils
analysent un engrenage avec module m=4, angle de pression $=20*, nombre de dents z=45,
déport x=û. Une comparaison avec la Méthode des Élkments Finis montre une différence de
seulement 0.8% pour le côté tension et 1.8% pour le côté compression.

Objectifs

Le but de ce travail est de poursuivre les travaux de Cardou et de Tordion concernant

l'application de la Mtthode des Potentiels Complexes aux engrenages droits. Comme
jusqu'ici, la MPC n'a été appliquée qu'au cas d'une dent isolée fixée à un massif s e m i - i n f ~ ,
on s'est proposé d'améliorer cette solution en modifiant la transformation conforme de Hirano
pour obtenir Ia représentation de 3 dents: ceci permettra d'établir de façon plus rigoureuse
l'influence des dents adjacentes sur la dent chargée. On s'est proposé aussi d'étudier les dents
d'engrenages internes.

Le deuxième chapitre présente la Théorie des Potentiels Complexes ainsi que son
application aux engrenages droits. Le troisième chapitre étudie la question de la
représentation de trois dents adjacentes. On en déduit la transformation conforme
correspondante. Le quatrième chapitre utilise cette transformation pour le calcul des
contraintes. Dans le cinquième chapitre nous étudions l'application de la MPC au cas des
engrenages internes. Ceci nous amène à la conclusion et à quelques recommandations.

Pour ne pas alourdir le texte, nous avons mis dans l'annexe 3 les programmes
informatiques utilisés. Ils sont écrits dans le langage MATLAB [The Math Works Inc.,
19891.
 PRÉSENTATION DE LA MI?TRODE DES POTENTIELS COMPLEXES

2.1. Variables com~lexeset fonctions de variables comr~lexes

Dans la Théorie de luasticité Plane, beaucoup de problèmes peuvent être résolus

avec des solutions math6rnatiques assez simpIes, en utilisant les variables complexes.

2.1.1. Variables compIexes

La notion de variable complexe est née avec le besoin de trouver la solution d'une
équation de la forme x2 + 1 = O. Comme il n'existait pas un nombre positif ou négatif dont le
carré soit égal à -1, Ies mathématiciens ont appelé le symboIe fl "nombre imaginaire". Les
nombres "imaginaires", en tant que racines carrées de nombres négatifs, sont mentionnés
pour la première fois en 1545 siècle par Cardan 1211. Jusqu'au milieu du XWI-e siècle, les
nombres complexes n'apparaissent qu'épisodiquement dans les travaux de certains
mathématiciens: Newton, Bernoulli [2 11.

Le mérite tout particulier en revient au grand mathématicien Leonhard Euler qui doit
être considéré comme l'un des créateurs de la théorie des fonctions d'une variable complexe.
Dans ses travaux, Euler étudia en détail les fonctions élhentaires d'une variable complexe, y
compris les fonctions logarithmique, exponentielle, trigonométriques, ainsi que les fonctions
trigonométriques inverses, donna les conditions de différentiabilité et les dl6ments du calcul
intégral des fonctions d'une variable complexe. Il d o ~ aussi a de nombreuses applications
des fonctions d'une variable complexe se rapportant a différents problèmes mathématiques. ïi
ouvrit l'ère de leurs applications en hydrodynamique et en cartographie. Le symbole "i"a été
introduit par Euler [21] et cette notation est utilisée encore aujourd'hui.
 On appelle variable complexe toute expression de la forme:

où: i représente l'unité imaginaire et i2 = -1 ;

x représente la partie réeue de la variable complexe;
y représente la partie imaginaire de la variable complexe.

L'interpr6btion géom6trique des nombres complexes remonte en 1799 et est connue

sous le nom de diagramme d'Argand [9].
La variable complexe peut être représentée par le point " A(x,y) ",dans Ie plan
" x - y ". A(x, y ) est I'affixe de 2.

Figure 2.1. Diagramme d'Argand ou plan z

La variable en coordonnées polaires a la forme:

On peut interpréter comme Ctant le vecteur OA. La variable représentée par le

vecteur OA' est donnee par

et eiie s'appelle la conjuguée de la variable 2.

2.1 -2. Fonction de variables comdexes

Une fonction de variable compIexe 2 est appelée fonction complexe. De façon

similaire au cas d'une variable complexe, une fonction complexe peut-être séparée en une
partie rée1.eet une partie imaginaire:

La fonction conjuguée de la fonction complexe f ( z ) est définie comme:

Une fonction complexe f(z) est appelée holomorphe ou analytique dans un domaine
a,si elle possède une dérivée unique en tout point du domaine D. Si l'on veut trouver les
dérivées partielles de f ( z ) on obtient, après la séparation des parties réelles et imaginaires les
équations:

Ces équations sont connues sous le nom de "conditions de Cauchy-Riemann" et elles

sont nécessaires et sufisantes pour que la fonction f ( z ) soit holomorphe.

En termes de dérivées partielles, la déride d'une fonction complexe peut s'exprimer

ainsi:

Si l'on elimine v dans les conditions Cauchy-Riemann, on a:

 De même, si l'on élimine u dans les équations on trouve:

Les équations ainsi déterminées démontrent que les parties réeiie et imaginaire d'une
fonction holomorphe de variable complexe sont des solutions de I'équation de Laplace.
Comme l'équation de Laplace est aussi appelée équation harmonique, les fonctions u et v
qui satisfont Féquation harmonique, sont appelées fonctions harmoniques.

2.2. Transformation conforme

La notion de transformation conforme est particulièrement importante en

mathématiques. Née de reprtsentations physiques, elle a de nombreuses applications
fondamentales dans les différents domaines: la méthode des représentations conformes
permet de résoudre les problèmes pratiques d'hydrodynamique et d'aérodynamique, des
théories des champs électrostatique, magnétique, thermique, de la Théorie de l'Élasticité.

Certains problèmes concernant les représentations conformes ont été résolus par
d'Alembert, Euler et Gauss [21]. En se référant à leurs travaux, Bernhard Riemann, dans sa
thèse intitulke " Fondements de la théorie générale des fonctions d'une variable complexe",
(1851) a posé la première pierre de la théorie géométrique des fonctions et a démontré
notamment (bien qu'avec des erreurs) le théorème fondamental de la possibilité de la
représentation conforme de tout domaine simplement connexe sur un autre. Les
représentations conformes sont largement appliquées comme instrument mathématique de
l'ktude de la mécanique des milieux continus. Parmi les promoteurs de ces applications,
citons pour l'hydrodynamique et l'aérodynamique Joukovsky et Tchaplyguine et pour la
Théorie de l'hsticité7 Kolossov et Muskhelishvili.

Transformations ou Représentations

L'ensemble des équations:

d é f i t en g6nkal une transformation ou une représentation, qui établit une correspondance
entre les points du plan des u - v et les points du plan des x - y. Si à chaque point du plan
des u - v correspond un point et un seul du plan des x -y on parlera d'une transformation
biunivoque.

2-2-2. Représentation conforme

Supposons que par les équations (2.10) le point (x,,,y,) du plan des x- y soit
transformé en le point (u,, v, ) du pIan des u - v tandis que les courbes Cl et C2 [se coupant
en (x,, y,) ] sont respectivement transformées en C 1 et C'2 [se coupant en (u,, v, ) J. Une
transformation telle que I'angle entre Cl et C2 en (x,, y,) est égal en grandeur et sens, à
l'angle entre C'1 et C'2 en (u,,~,), est dite conforme en (x,, y,). Donc f est la
transfonnation conforne de D sur D'si f est holomorphe et bijective.

Figure 2.2. Représentation conforme

2.3. La Méthode des Potentiels Com~lexes

La methode qui utilise les variables complexes s'appelle la méthode des potentiels
complexes. C'est une méthode mixte, qui a une partie analytique - associée 2i la Théorie de
d élasticité, et une partie numérique - associée à la détermination de la transformation
conforme.
 Dans I'élasticite plane, les contraintes peuvent être calcu16es partir de trois
équations:
- deux équations d'6quilr%re(avec des forces de volume nulles)

- une équation de compatibilité (avec des forces de volume nuIles)

où v est Ie coefficient de Poisson du matériau,

Si i'on combine Ies trois équations différentielles en une seule équation du quatrième
ordre, celle-ci dépend d'une fonction unique appelée fonction de potentiel F(x,y) ou
fonction d'Airy. Cette fonction a étk introduite en 1862 par l'astronome anglais Airy [3 11.

Avec la fonction d ' A j r les contraintes s'expriment ainsi:

et les équations d'équilibre sont satisfaites pour toutes les valeurs de x et y. L'équation de
compatibilit6 devient une équation biharmonique:

Le problème qui se pose, pour obtenir les contraintes en tout point d'un solide, est de
trouver la fonction d'Airy respectant les conditions aux limites et I'équation bihamionique.

Toute fonction biharmonique peut être représentée à l'aide de fonctions holomorphes

d'une variable complexe. Puisque les composantes de la contrainte et du déplacement
s'expriment par une fonction biharmonique, nous anivons ainsi à la représentation complexe
du problème plan de la Théorie de 16lasticité.

Le nom de MÉTHODE DES POTENTlELS COMPLEXES (MPC) vient du fait

qu'on associe la fonction de potentiel F ( x , y ) à une fonction holomorphe complexe U ( z ) (où
z=x+Q).
Ainsi, l'équation biharmonique devient

Le premier auteur à avoir démontré que l'on pouvait exprimer U à l'aide de deux
fonctions holomorphes de variable complexe @(z) et ~ ( z )a ,été Goursat en 1898 [3 11. La
solution de l'équation biharmonique a ainsi la forme:

où Re signifie la partie réelle de l'expression qui se trouve entre les crochets. Les
fonctions @ ( z )et ~ ( zdépendent
) des conditions limites du problème: géométrie et charges
appliquées.

On voit facilement l'avantage de la méthode des potentiels complexes: si dans

l'équation biharmonique en fonction de Fnous avons des dérivées partielles par rapport à x
et y, en utilisant les variables complexes on utilise une seule variable z , et donc des dérivées
ordinaires.

Pour un étude plus approfondie de la théorie des variables complexes nous

recommandons les publications [4,9,20].

2.4. Éauations de base dans la Théorie de l'Élasticité Plane

2.4.1. Eauations de Kolossov

Si les fonctions @(z) et ~ ( zsont

) connues, les composantes des contraintes a,,
a,,r, peuvent être détermin6es directement. Cest Kolossov (241 qui a appliqué la théorie
des variables complexes à un domaine quelconque, sirnpiement connexe. Il a ainsi obtenu
deux équations donnant les contraintes et une équation pour les déplacements.

OX + o,= 2[Qf(z) +W(z)J

q - uX+2ir, = 2[Z - W(z)+ Y f ( z )J
---
2p(u + b)= K - @(z)- z - @'(z)Y(z)

où K depend du type de problème:

A +3p
K= -
=3 4 v pour l'état de déformation plane, avec A , p = paramètres de Lamé
A+P

pour l'état de contrainte plane

Il faut donc trouver deux fonctions @(z) et Y(z) =x'(z) nomm6es potentiels
complexes, fonctions de variable complexe z = x + iy qui sont holomorphes dans le domaine
et satisfont les conditions aux limites.

2.4.2. Condition limite

On prend une abscisse curviligne s sur le contour d'un solide.

Figure 2.3. Force agissant sur un segment AB d'un solide

 Les fractions sont données par les composantes T,(s) et T,(s)

T, (s)= ox- cos(a) + ,z - sin(@

Pour un problème du premier type,la limite que l'on impose à none solution est une
force T appliquée sur le contour du domaine, le potentiel U(z) et les fonctions @(z) et Y(z)
doivent satisfaire à:

pour en être la solution.

2.5. Pomaine semi-infini soumis B une force concentrée

On démontre [24, 301 que les fonctions Q(z) et Y(z) pour un domaine semi-infini
quelconque, sollicité par une force concentrée, s'expriment sous la forme:

où y et y' sont des constantes complexes.

Il s'agit de trouver les fonctions W(z) et Yœ(z)holomorphes dans le domaine D, y

compris à I'înfïï (elles tendent vers zero Zi I'innni), ainsi que les constantes complexes y et
y' qui verifient les conditions aux limites.
 21

Pour une force W appliquée en un point z, de la frontière du domaine, les fonctions

potentieis deviennent [q:

Figure 2.4. Force appliquée à un point quelconque z, de la frontière du domaine

où @'(z) et YD(z)sont des fonctions holomorphes.

Pour certains domaines, la fonction d'Airy peut être dE

i de à trouver et i'on utilise
une transformation conforme qui modifie le domaine onginai et ses conditions limites en un
domaine plus simple à étudier. On peut se permettre d'utiliser la transformation conforme car
elle a la propriété de ne pas altérer la fonction biharmonique en passant d'un domaine à
l'autre, si la transformation est une fonction holomorphe.

2 -6.&plication de 13 Méthode des Potentiels Com@exesaux enmenages droits

Pour appliquer cette méthode aux engrenages droits, on considère une dent isolée sur
un demi-plan. Une transformation particuiière due à Hirano [17] la transforme dans un demi-
plan. Les équations de I'ékisticité plane classique peuvent ensuite être appliquées puisque le
profil de la dent est devenu une ligne droite.

Les principes de base des travaux effectues jusque là ont été:

 - de considérer que les dents adjacentes ont un effet négligeable;
- de considérer que I'engrenage est plein;
- de déterminer une transformaton conforme qui suive d'assez près le profd
de la dent,

Les hypothèses faites permettent la transformation du probIème réel, donné par deux
dents en contact, en un problème équivalent comportant une seule dent sous une charge W,
normaie à la surface. Sous l'effet de la transformation conforme la dent se transforme en un
demi-plan.

Figure 2.5. Deux dents en contact transformées en demi-plan chargé par une force W

La transformation de Hirano est la suivante:

où a,, bk et c sont des coefficients réels qui varient avec la forme de la dent. Cette
transformation dépend de 2n+l paramètres selon le nombre n de termes que I'on utilise.
Si I'on exprime cette relation pour chacune des coordom&s x et y on obtient:

La frontière du demi-plan inférieur ( 7 c 0) correspond ii 7 = O. Donc l'équation du

profil correspndant dans le plan ( z )est dom& par:

2.6.1. Ex~ressiondu problème d'élasticité a~rèstransfomiation conforme

On considère un espace semi-infiniborné par un contour ayant la forme d'une courbe

symétrique par rapport à I'axe Oy et asymptotique à Ox.

Figure 2.6. Force appliquée en un point quelconque z, de la frontière d'un espace semi-
inhi

On y applique une force concentrée W au point q, , faisant un angle avec I'axe Ox,
soit weiP.On a vu que les fonctions @(z) et Y(z) peuvent s'tcrire dans ce cas par les
équations (2.25) et (2.26).

On effectue la transformation z = w ( c ) ramenant ce profil à Paxe réel dans le plan

(C).La fonction o(&Jest holomorphe dans tout le demi-plan ( q < 0).
Soit Co I'image de q (& est sur I'axe réel).
 On trouve ainsi:

deuxième terme est une fonction holomorphe dans le domaine et tend vers zéro
lorsque (4 -. Par ailleurs, on peut négliger la constante Logc puisque, pour un état de
, fonction @(z) et Y(z) sont définies à une constante près (à une
contrainte d o ~ éles
fonction Iinéaire près pour @(z) ) .

On peut donc écrire:

Dans Y(c) (Cquation 2-26),on a le terne

La fonction

a un pôle à = goet tend vers zéro lorsque +a.

Eue peut s'écrire

1
ffU= +f 'KI
@'f C,)K - Co)

où f '(CI est une fonction holomorphe dans tout Ie plan, y compris à l'infini.
On peut donc écrire:
où la fonction Y'(c) est holomorphe, y compris à l'infini.
On pose:

Condition limite

Soit un point à l'intérieur du domaine.

Figure 2.7. Espace semi-infini

La condition limite (2.22) devient, lorsqu'on I'exprime dans le plan c:

Dans le cas présent, la frontière du domaine est transfom6e sur I'axe réel donc, si
l'on appelle t i'abscisse le long de cet axe (on pose t = 5):
 B
-
@(t)+ =W ( t )+F(t)= - i l (T,( t ) + iTy(t))dt+CS
@(O
+

(0 A

(le signe - appadt car le domaine est à droite quand on parcourt l'axe de gauche ii droite).

. .
2.6.2. Transformation de la condition h t e

On remplace @(t) et Y(t) daas l'équation (2.4 1), par leur valeurs (2.34) et (2.38):

t
ao(t)
W)-
@b(t)+ K(t)
(t)
+ @ * ( t )+--
M t )
?(t) +fi)
= -i 1( T ~+iTy (t) (t))dt +CSI
-m

@o(Q+yo(5-)=--*
weiB
LO& - t o )+-
we@
L O ~- t,,)
( ~ +- W O ) *
we-4
21t 2n: ' ) ( - t0) 2~

-
we" T-?O Mt0 we-$
qJC)+'u,(C>=-- 2x
hg-
@O
+-@'(&)(f -t,)
.-
~IL.

On analyse le premier terme:

Donc, si l'on prend un point de référence A à gauche de t, (', < t,), et un point
courant B:
 On voit donc que l'expression

donne la résuitante des tractions imposées sur l'arc AB (multipliée par le facteur -i), ce qui est
exactement la valeur de l'intégrale du membre de droite de I'équation (2.42).

Ainsi, i'équation (2.42)devient:

- J--..(T,( t )+iT,(t))dt

Mais
 Si I'on combine les équations (2.55) et (2.52) on obtient:

2.6.3. Calcul des potentiels Of (r) et y(r)

2.6.3.1. ÇalcuI du wtentiel

L'équation (2.52) s'&rit sous la forme:

On multiplie chaque terme de l'équation par.

p.
1
-

21n' t - c

où est un point à l'intérieur du domaine, et I'on intègre l'équation sur l'axe réel:

9 ' ( t ) est la valeur limite de <Om(r),

holomorphe dans le demi-plan inférieur (oii se trouve
C), et tend vers zéro à I'infini.
Les propriétés des intégrales de Cauchy sont connues:

où r est un contour fermé qui entoure une seule fois le point z du domaine D
L'index de Ia c o u .par rapport & est:

( z=1 si ZED

hd&) = O si zED

Donc

(le signe "-" apparafAtcar le domaine est à droite quant on parcourt I'axe de gauche à droite)
De même, \ r ( t ) est la valeur limite de F(0,holomorphe dans le demi-plan
supérieur puisque:

Cette fonction tend vers zéro à I'innni, donc:

(car se trouve dans le demi-plan inférieur et la fonction est holomorphe dans celui
supérieur).

A partir de l'équation (2.34) on obtient:

A partir de ï'équation (2.34) on peut aussi CalcuIer

weip
@'(t) = --- 1 + co" (t)
21r t-t*

-
a'(?)
,---we-'B 1 -( t )
+w
2r t-to

Si i'on remplace (2.67) dans la première intégrale de l'équation (2.64)on obtient:

Si l'on combine (2.68) et (2.64) on obtient:

Les deux integrales dans l'équation (2.70)deviennent avec la transformation de

 - -
La fonction --
W(t)
- c,
@'(O t
t
est la valleur limite de la fonction -- t
@'*(O
t-51
~'(0 '
holomorphe dans le demi-plan supérieur puisque @'*(Oet m'(O le sont dans le demi-plan
inférieur, et m'(O ne s'annule pas dans ce demi-plan. 5; est égaiement un point du demi-
plan inférieur. Cette fonction tend vers zéro lorsque C ++ =. Donc, la première integrale,
dans I'equation (2.72) est nulle. La seconde intégrale est du type Cauchy, - -
f(0
1 -dt,
2m'-t-ib,
aka"( t )
où f(t)=C- dt est une fonction holomorphe dans le demi-plan sup&ieur, y
M t - 5; )
compris à l'infini. D'après les propriétés des intégrales de Cauchy:

et on obtient

Pour la deuxième intégrde dans (2.70) on a:

L'intégra.de est la valeur limite de la fonction:

Le premier terme est holomorphe dans le demi-plan supérieur et son intégrale est
nulle. Si l'on remplace a>(T) par la transformation (2.27). le deuxième terme donne des
expressions de la forme:
 1
Comme la fonction f (0= - est holomorphe, y compris à
@ ' m C - 4, )K - CI 1
I'infini, dans le demi-plan supérieur, on obtient

we-""
r2 =--
2a &(ibk)(ibk -cl)(* - t,)
=k -l

ou dans une forme simplifi&x

Donc

où l'on a remplacé 5; par C, car <,est un point du demi-plan inférieur.

2.6.3.2. Calcul du potentiel

On part de l'équation conjuguée de (2.57):

On la multiplie par -1 1
-
21Cit-Cl
et on l'intègre sur I'axe réel ( cl est un point du demi-
plan inférieur):
 Y*(?) est la valeur limite de ~'(0 ,
holomorphe dans le demi-plan infdrieur (où se
trouve C,) et tend vers zéro à l'infini. D'après les propriétés des intégrales de Cauchy on a:

(2.84)
- -
Comme @"(t) est la valeur limite de ~ " ( c )holomorphe
, dans le demi-plan
supérieur et tendant vers zéro à l'infini, on a-

On remplace C, par 5 et l'équation (2.83) devient ainsi:

Si l'on combine (2.66) et (2.86) on obtient:

Pour la tramformation codonne de Hirano, la conjuguée prend la forme:

Comme cette fois les pôles C = -ib, sont dans le demi-plan inférieur, on a dans la
première int6grde des temies de la forme:
ce qui donne

Donc, 11t5quation(2.89) devient:

En même temps les deux autres intégrales dans (2.87)donnent:

Comme les pôles sont dans le demi-plan inférieur, en utilisant (2.88),on a des termes
de la forme:

1 1
La fonction f(0= -- est holomorphe dans le demi-plan infkrieur et
@'(CI c - 5o
I'équation (2.93) devient:

Le terme donné par c c étant une constante, on peut le négliger pour un état de
contrainte donné. Ainsi, si l'on combine (2.54), (2.88) et (2.94), le deuxième terme de
(2.92) est égal à:
 Le premier terme dans l'équation (2.92) est égal à:

Si Itonremplace les équations (2.91), (2.95) et (2.96) dans (2.86)on obtient:

-
we* O(co)--- 1 "(O[ab(n+w*cn]+
Y*(() = --
T-Co
21~.a'<To> a'(D

Comme on connaît la forme de y({) donnée par Iféquation(2.38), on remplace

Y9(r)
par sa valeur (298). et l'on obtient la valeur du potentiel Y (0:

Transformation des éauations de Kolossov

Les équations de Kolossov peuvent être refomulées dans le plan ( 0:

où t et n sont les directions tangentielles et n o d e s à la coordom6e curviligne 5, dans le
système de coordom6es curvilignes qu'on a lié à la dent d'engrenage.

Figure 2.8. Système de coordonnées curvilignes lié à la dent d'engrenage

Comme la transformation conforme conserve les angles, les courbes { s s t et q s s t

restent perpendiculaires entre eues dans les deux systèmes de coordonnées. C'est pour cela
qu'on les considère comme un système de coordonnées curvilignes lié à Ia dent.

2.6.4.1. Explicitation de

Si l'on combine les equations (2.8 1) et (2.39a) on obtient:

d'où
 et I'on a n equations complexes:

avec:

On pose:

ajk= . 1
(b, + 4 )' d ( - i b k )

Les potentiels @(ib,) et a'(-&), que l'on trouve dans I'équation (2.105) sont
remplacés dans les calculs subséquents par les constantes réelles X, et Y, (k=l, ...n). Ces
constantes sont fonction des coefficients de la transformation conforme, de la position et de
I'angle de la force sur la dent. Ainsi les potentiels @(il+) et a'(-ibk) s'expriment de cette
façon:

Pour calculer (-ibj ) on utilise la formule (2.54):

 avec

On est arrivé ainsi à un système de 2n équations linéaires:

Ce système dkquations étant résolu, on obtient:

2.6-4.2. Calcul des contraintes superf'icielles

Les contraintes superflcieiles sont situées à la surface de la dent, donc pour les points
à la frontière du domaine, c+= t

Pour simplifier l'écriture, nous dons noter

w cosp n [xI(f2 bt)+ 2ykbkt]

F,,( t ) = ---
+
(t2 b ~ ) ~
et i'équation (2.1 18) peut s'écrire ainsi

L'équation (2.106) on I'écrit

Pour simpliner I'écrinire,nous allons noter:

et I'équation(2.122)peut s'écrire:

~ ' ( t=)Tl,( t )+ iT12(t)

Mais, dans la contrainte supeficielle on a besoin de:

Comme la limite du domaine est une courbe k=cst, ia contrainte normale à la surface
a,(tangentielle à la courbe c=cst) est nulle, ainsi que le cisaillement r,. Les équations de
Kolossov (2.100)et (2.102 ) deviennent:
 or( t )= rm( t )= O

2 -6.S. Mode%sation de la dent d'engrena~eavec la transformation conforme

Pour c a l d e r les contraintes-dans une dent d'engrenage chargée au point (x, ,y,), par
la Méthode des Potentiels Complexes, il faut d'abord modéliser cette dent avec la
transformation conforme. La précision des formules donnant contraintes et déplacements,
analytiquement exactes, dépend essentieuement de la précision des transformatons obtenues.
La modélisation sera réussie lorsque le pronl de la transformation conforme, généré par les
équations (2.29), suivra parfaitement le profil de la dent. Le problème est donc de savoir
combien de termes n il faut prendre, et dans quel sens faire varier les paramètres (ak,bk),
pour s'approcher autant que possible du prof3 donné.

Figure 2.9. Dent d'engrenage à modéliser

 Dans la Figure 2.9, & est le rayon de tête, R est le rayon primitif. Rr est le rayon de
pied, 1est la partie développante, et 2 la partie trochoïde.

Le type de dent que nous d6sirons modéliser est une dent droite symétrique, générée
en deux parties: une développante de cercle et une trochoïde (engendrée par I'mondi de
I'outil-ccémaillère). Le système de coordomées X j , utilisé pour définir le profd a son
origine à I'intersection de son axe de symetrie et le cercle de pied de la dent. On remarque
qu'une partie de l'extrémité de la dent se trouve dans les j nkgatifs.

2.6.5 -2.Aiustement des axes de coordonnées

Comme la transformation conforme se situe toujours au dessus de I'axe Ox, il faut

ajuster le système d'axe de la dent (I, y), pour que celle-ci soit modélisable à partir du
système d'axes de la transformation conforme ( x , y). Cet ajustement se fait avec un
coefficient de translation d'axe, appelé d, un paramètre positif. qui ramène le système de
coordonnées de la dent au niveau de celui de la transformaton conforme.

Figure 2.10. Translationdeladent

Ainsi:
 Avec ce coefficient. I'extrémité de la dent devra tendre de façon asymptotique vers
I'axe Ox quand x tendra vers I'infi (km). Le param8tre de translation d doit être calcult
en même temps que les paramètres de la -ormation conforme.

2.6.6. Recherche de Ia transformation conforme ~ a Iar roma am mat ion non-linéaire

La recherche des paramètres consiste à minimiser une fonction U ,indicatrice de la

précision avec laquelle la fonction de Hirano représente le profil théorique de l'engrenage
étudié. Si i'on appelle xi I'abscisse d'un point du profd théorique, et I'abscisse donnée
par la transformation pour la m ê m e ordonnée, il s'agit de minimiser la distance entre les deux
profils pour obtenir les coefficients exacts. La distance entre les profils est représentée par
une fonction de type "moindres-carrés":

où N est le nombre de points de la transformation, compris entre le sommet et le pied de la

dent, et w, est un terme de poids permettant d'douer plus ou moins d'importance à certains
termes de la fonction des moindres-carrés.

La minimisation d'une fonction quadratique telle que U est un problème classique et

plusieurs méthodes sont disponibles: Gauss-Newton. Marquardt, Powell etc. Deux
méthodes de recherche directe ont été utilisées dans les travaux de Cardou et Tordion 151:
a) l'algorithme de Rosenbrock (avec un programme FORTRAN pour la recherche des
coefficients de la représentation conforme)
b) l'algorithme de Powell (avec un programme APL pour la recherche des
coefficients de la représentation conforme)

Dans l'algorithme de Rosenbrock qui agit sur les coefficients de la transformation

conforme, on explore xi dans p directions orthogonales, p étant un nombre de variables
indépendant. L'exploration dans chaque direction se fait avec un pas choisi. Après chaque
étape on vérifie la fonction à minimiser: si eiie est trop petite, alors on réalise une rotation
d'axes selon les directions de la plus grande décroissance trouvée dans cette étape.
 Dans la d e u x i h e méthode qui utilise i'aigorithme de Powell, on effectue une
recherche unidimensionnelle, dans chacune des directions, en utilisant une approximation
quadratique. La dent est gknérée par la transformation donnée par les équations (2.29) dont
on a soustrait le coefficient de translation d'axe d (pour la coordonnée y). L'algorithme de
Powell fait varier les coefficients, ce qui a pour effet de changer les coordonnées du pronl de
la transformation conforme, et par le fait même la valeur de la fonction II. Comme les
contraintes les plus élevées de la dent sont situées dans la partie inf'rieure (concave) de celle-
ci, on doit apporter un effort supplémentaire & la modélisation de cet endroit. C'est pour cette
raison que le terme de poids wj devient égal à 2. Pour la partie supérieure (convexe) de la
dent, le terme de poids a une valeur unitaire.

L'algorithme de Powell semble nettement pius efficace que celui de Rosenbrock en ce

qui concerne le résultat final. Les tract% des courbes obtenues sont indiscernables, ii l'échelle
1: 1, sauf dans la partie haute, qui présente peu d'intérêt pour Ifanalysedes contraintes. La
région de la section critique, au pied de la dent, est très bien rendue [a.

On sait que la transformation de Hirano dépend de 2ntl paramètres selon le nombre

n de termes que l'on utilise. Pour le profîl d'un engrenage standard, n=4 donne des très
bons résultats [SI. La transformation conforme des profils étudiés fait intervenir 10
paramètres (a,,4, c , d ). Les nombres de termes dans la transformation conforme dépend
de Ia forme du profd: pour un petit nombre de dents on a un profil courbe qui demande peu
de termes dans la transformation conforme; pour un grand nombre de dents, on s'approche
du proN d'une crémaillère avec des faces rectilignes, qui sont plus difficiles a représenter par
la transformation de Hirano.
 TRANSFORMATION CONFORME D'UN DOMAINE SEMI-INFINI DANS LE
DEMI-PLAN. DANS LE CAS DES 3 DENTS

Jusqu'ici la Méthode des Potentiels Complexes n'a été appliquée qu'au cas d'une
dent isolée fixée à un massif semi-infini. Pour un calcul plus précis des contraintes au pied de
la dent on s'est proposé d'établir Ifinnuencedes dents adjacentes sur la dent chargée et on a
essayé d'améliorer la solution utilisée jusqu'à présent. Pour cela on a modifié la
transformation conforme proposée par Hirano [17], pour une dent, afin d'obtenir la
représentation de 3 dents.

3.1. Translation de la dent centrale

Pour commencer, on supposera que les 3 dents sont obtenues par des simples
translations sur l'axe x de la dent centrale, à droite avec -8 et à gauche avec 6.

Y L " 4
Figure 3.1. Translation de la dent
 Pour simplifier. on considérera la transformation de Hirano (2.27). seulement à 1
texme:

La dent (O) seule correspond

La dent (1) d e correspond à:

La dent (2) seuie correspond à:

Les trois dents peuvent être représentées par:

z=q+ztf 22

On sait que (= Pour les points du contour 9 = 0. Si l'on note 5 = t , on

+itl.
reste avec = t . Comme z = x +j , après la séparation des parties réelle et imaginaire on
obtient:
NOTE
Les dents (1) et (2) seront Iégèrement différentes de la dent centrale (O). Cette
diffdrence dépendra de la vaieur 6 avec laquelle on a effectué la translation.

L INFLUENCE, SUR LES 3 DENTS, DE LA TRANSLATION SUR x

Figure 3.2. L' influence de la translation sur "x" sur les dents d'engrenage droit

Comme on peut Ie constater dans la figure 3.2, la valeur 8 influence aussi sur Ie
prohl de la dent centrale, surtout au niveau du creux de la dent:
Courbel: correspond à une translation &=O -5
Courbe2: correspond à une trardation 6=1-5
Courbe3: correspond à une translation 6=3.O

Pour la représentation des trois dents simplement translatées nous avons utilisé le
programme 1NFLU.M (annexe 3).

Les graphiques présentés dans les figures 3.5, et 3.7 à 3.16 sont obtenus avec le
programme ETAPE.M (annexe 3). Leur échelle est en pouces et les coefficients de la
transforation conforme 151 sont:
3 -2. Rotation d'un ~rofil

Suite aux observations faites, on va étudier la rotation d'un profil de dent

Figure 3.3. Rotation d'un profil

La ûansfonnation de Hirano (3.1 ), est déjà connue:

La rotation du profil se traduit, du point de vue mathématique, par une multiplication

par eie ,et ceci nous amène a
 Si l'on note:

et on connaît que:

on obtient:

Après la séparation des parties réelle et imaginaire on obtient

3.2.1. Position de l'axe de svmétrie du profil

O n doit maintenant vérifier la position de l'axe de symétrie du profd qui a subi la

rotation. Pour cela, on considère 5 = O dans les formules (3.17) et (3.18) et I'on obtient:
 X
-=- -Y
sine C O S ~

La formule (3.22) est la représentation d'une droite inclinée de 8 + a /2, ou

-
rr / 2 8 , par rapport à Ox.

Figure 3.4. Position de I'axe de symétrie du profil

3 -2-2. Position du mofil

On vérifie aussi la transformation donnée pour le profil. Pour cela, q = O et 6 = tet

les f o d e s (3.17) et (3.18) deviennent:

x= (ct+- r~~b')CoSe-tZ+bZ
ab
sin 8

Pour t +f=. x = ct cos 6 et y = ct sin 8. Cela donne la pente de I'asymptote:

 POSlTiON D UN PROFIL MPC N R E S ROTATiON
l
I '+ I
I 1

Figure 3.5. Position d'un profil MPC après rotation

3 -3. Transformation conforme oour les 3 den^

Pour avoir une continuité dans la transformation et aussi une meilleure approche avec
un engrenage, on propose, pour les dents adjacentes, des translations avec -6 tg8 1 2 sur
i'axe y des dents déjà translatées.
Pour les 3 dents, on propose donc, comme nouvelle transformation celle obtenue
ainsi:
- des trans1ations avec 8 et -6 & la dent centrale sur I'axe x ;
- des translations avec -6 - tg0 / 2 sur l'axe y des dents dejil translatées;
- des rotations de eie et e-i8 par rapport à i'axe y des dents translatées.

Figure 3.6. Transformation conforme proposée pour trois dents

 On analyse les points du contour, donc q = O et I'on garde la notation = t .

3 -3.1. Transformation de la dent centraie

kitransformation de &an0 nous donne:

a-t
Xo= c - t +
t2 +b2
a-b
Y. =-
t2 +b2
DENT CE-, MPC
1 I I 1 I
I

-- t------------

Figure 3.7. Dent d'engrenage droit tracée par MPC

3.3.2. Transformation de la dent (Il. dent translatée à droite

a) La dent subit une translation -6 sur l'axe x , et on note

 Dans ce cas. la transformation de Hirano, pour z, = x, + j,,
prend la forme:

a-b
YI=
t; +b2
En regardant la figure 3.8, on observe que la représentation au pied de la dent, après
la translation subie, n'est plus symétrique, fait justifié exactement par la translation S.

D E N T W S ~ A T É ESUR L AXE x
3 I I I I 1 t

Figure 3.8. Translation sur l'axe "x" de la dent centde

b) La dent translatée subit une translation -6 - tg8 1 2 sur i'axe y

Ainsi, Ia transformation devient:

 a-b 8
-6-tg-
'?l=t:+bz 2

Figure 3.9. Translation sur I'axe "y"de la dent centrale

c) La dent ainsi translatée subit une rotation de e-iû -

- cos 6 - isin 8 et la
transformation finale de la dent (1) est de la forme:

Si l'on représente ensemble les deux dents analysées, celle centrale et celle qui a subi
les transformationsprésentées, la transformation codonne correspondante donne le résultat
obtenu dans la figure 3.1 1.
 DENT MPC APRES TRANSLATlON SUR x. SUR y, ET ROTATlON

Figure 3.10. Dent d'engrenage droit MPC,qui a subi des translations et une rotation

DEUX DENTS MPC translation à droite

I 1 I
t
I

.
, .
8 8

-6 -4 -2 O 2 4

Figure 3.1 1. Dew dents d'engrenage droit représentées par MPC (translation à droite)

Le profil marqud par "*" represente le profü de la dent centrale; le proffi marqué par "O"
représente le profil de la dent translatée.
3.3.3. Transformation de la dent 12). dent translatée à gauche

En suivant Ies même étapes de calcul on obtient:

a) La dent subit une translation S sur I'axe x , et on note

t+S

Dans ce cas, la transformation de Hirano, pour 2, = x, + ly, prend la forme:

8
b) La dent translatée subit une translation -6 - tg- sur I'axe y
2
Ainsi, la transformation devient:

a-b 8
2 -&tg-
Y2 = t; + b 2

c) La dent ainsi translatée subit une rotation de eie=cos O+isin 8 e t la

transformation f d e de la dent (2) est de la forme:

x2= ( c - t 2 + tz2a +- t2b2 ) - ~ ~ s et:-a+-6b 6

-sinû+&g-sin0
2

Si I'on represente ensemble les deux dents analysées, celle centrale et celle qui a subi
les aanfonnations présentées on aura la representation de la figure 3.13.
 DENT MPC APRES TRANSLAllON SUR x, SUR y. ET ROTAllON

Figure 3.12. Dent d'engrenage droit MPC, qui a subi des translations et une rotation

DEUX DENTS MPC translation a gauche

3 1
I
1 1 I

8
,' #.

Figure 3.13. Deux dents d'engrenage droit représentées par MPC (translation à gauche)
Le profil marqué par "*" represente le profil de la dent centrale; le profil marqué par "O"
représente le profil de la dent translatée.
 Fimalement, pour les trois dents on a la tramformation conforme suivante:

On remplace (3.26), (3.33) et (3.40) dans (3.42) et l'on obtient:

On remplace (3.27),(3.35) et (3.41) dans (3.43) et I'on obtient:

ab a*~ 8
Y= ?2 +bl -(ctl + t; + b l -Ct2 - t:- + b2 cos 6 - 26 tg -cos 0
2
(3 -45)
TROIS DENTS MPC,TRANSLATION-SUR x ET y- ET ROTAION

S -6 -4 -2 O 2 4 6 8

Figure 3.14. ReprhentationdeeoisdentsparMPC

Pour valider les transformationsproposées pour chacune des trois dents, nous allons
étudier la position de chacune des dents dans la transformation finale:
 POSiiiON DE CHAQUE DENT
3.5
1 1 1 I 1 1 i a

Figure 3.15. Position de chaque dent dans la transformation pour trois dents

La courbe représente la tramifonnation conforme pour deux dents (avec

"-O"

translation à droite)
La courbe "-" représente la transformation conforme pour deux dents (avec
translation à gauche)
La courbe "-" représente la transformation conforme compIète pour trois dents
Le point "*"marque la postion de la dent centmie
Le point "O" marque la position de la dent adjacente à droite

Nous avons prouvé ainsi que les étapes que nous avons suivies pou.déterminer la
transformation conforme tiennent compte de la position de chaque dent par rapport aux dents
adjacentes, ainsi que du sens de rotation pour chacune des dents adjacentes.

Comme z = x + iy, on peut présenter la transformation conforme pour trois dents

dans une forme simplifiée:
 Nous avons choisi cette transformation (3.46) considérant que c'est ia meilleure
approche avec un engrenage. Si on la compare avec une autre tramformation, la Mi5rence
entre la transformation conforme compl2te que nous avons proposée et une autre, obtenue
par des simples translations de la dent centde est visible dans la figure 3.16.

DIFFERENTES TRANSFORMATIONS

Figure 3.16. Différentes transformations appliquées aux dents d'engrenage droit

La courbe "-" trace les pronls des dents simplement translatées; la courbe "-" trace les profiIs
des dents après la transformation conforme (3.46).

Donc, pour conclure toutes ces démarches, nous proposons une transformation
conforme (équation 3.46) qui peut représenter trois dents adjacentes d'un engrenage droit.
Maintenant, le problème est de trouver le nombre nécessaire de paramètres et leurs valeurs
pour que la courbe tracte par la transformation conforme proposée soit aussi proche que
possible du profil donné de la dent.
3 -4. Recherche des coefficients de Ia transformation conforme

Dans leurs travaux. Cardou et Tordion [5,6,7,8]trouvent que, pour modtliser une
dent isolée, le nombre suffisant de termes dans la transformation de Hirano est de quatre.
Pour commencer, nous avons aussi considéré une transformation à quatre termes:

Pour déterminer les coefficients de la transformation conforme il faut:

-définir le profil dome par un ensemble de points ( x i ,yi)
-appliquer l'algorithme de Powell avec des valeurs initiales pour les
paramètres de la transformation

Les détails des calcuis ainsi que le programme que nous avons utilisé se trouvent en
annexe 3.

Nous avons traité le cas d'un engrenage correspondant à un nombre de dents N=20,
angle de pression @=20°,pas diamétral P=l, rayon standard de l'outil r=0.38/P.

Pour trouver la représentation conforme, nous avons considéré les paramètres

constants suivants:

- angle de rotation de la dent: &3600/N 0=180

- translation de la dent sur I'axe x 6= 1.75

Les coefficients utilisés au départ et ceux obtenus se retrouvent dans le tableau

suivant:
 - - - - - - .. - - - -
Coefficients de départ Coefficients finaux
al 0. 1 0.09999
a2 0.3 0.27524

Figure 3.17. Profil approché d'une dent standard (pour le cas des trois dents)

La figure 3.17 donne les 3 courbes suivantes:

-courbe 1: le profil donné (obtenu par points avec GENERI)
-courbe 2: le profil de dépari (obtenu pour les coefficients de départ avec C0EFFI.M)
-courbe3: le profil approché (obtenu pour les coefficients optimisés avec C0EFFI.M)
Note pour la courbe 1: Pour définir le prof3 d'une dent d'engrenages droits, Cardou et
Tordion [SI présentent le programme GENERI. Ce programme définit le profil et,
pour chaque point il d o ~ les
e coordonnées (x, y), l'angle de la normale avec l'axe de
symétrie et le rayon du cercle sur lequel se trouve le point Pour le cas que nous avons
étudié, nous avons utilist les points tels que déterminés en [5] (tableau 1, page 6.4).
Pour tracer le profil de la dent adjacente, nous avons effectué la rotation du profil
central avec R0TAIION.M (annexe 3).

Comme on peut voir dans la figure 3.17, rapproche du profil tracé pour trois dents
se fait assez bien pour la dent centrale, spécialement pour le pied de la dent, zone d'effort
maximal, et moins bien pour celle adjacente. Tenant compte du fait que nous vouions voir
l'influence des dents adjacentes sur la distribution des contraintes daas la dent centrale, on
peut considérer la solution obtenue comme étant satisfaisante.

Si l'on compare les coefficients de la transfomiabon conforme pour une dent avec
ceux pour trois dents on observe la diminution du coefficient c qui passe de la valeur
c=0.94676 (pour une dent) à c=0.48775 (pour trois dents).
 CONTRAINTES SUPERFICIELLES DANS LES DENTS D'ENGRENAGES DROITS

Jusqu'ici, la Méthode des Potentiels Complexes n'a été appliquée qu'au cas d'une dent
isolee fïxfe à un massif semi-infi. Comme nous l'avons déjà mentionné, dans les travaux
effectués jusque Ià, on a considéré que les dents adjacentes ont un effet négligeable. Nous
nous sommes proposé de tenir compte des dents adjacentes et d'établir ainsi d'une façon plus
rigoureuse leur influence sur la dent chargée. Avec la transformation conforme de Hirano
pour obtenir la représentation de trois dents, présentée dans le chapitre précédent, nous allons
appliquer la Méthode des Potentiels Complexes dans le calcul des contraintes superficielles
dans les engrenages droits.

4.1. Calcul du potentiel @KI

Nous avons démontré dans le chapitre II que la fonction de potentiel @(apeut-être

exprimée par l'équation (2.34):

avec a*(5)donné par l'équation (2.64):

Nous ailons maintenant évaluer l'intégrale (4.2) pour la transformation conforme pour
trois dents, transformation qui est représentée sous forme simpW16e par l'équation (3.46):
 Pour simpliner les calculs on considérera la transformation seulement à un terme. Si
l'on exprime la transformation conforme sous une forme d6taill6e on obtient:

La première intégrale dans l'équation (4.2)devient avec la transformation conforme:

- -
.,
1
5
=- 1 " -~ d( ttW) =( t -)
21E o'(t) t -
i
2m'
' ( t ) 1
J --{ct+-+[c(t+6)+
o f ( t ) t - Cl
a
t-ib
a ]P kt +
--O -0 t+S-ib

- - -
a
+-23ri
-Q)
--
0 8 * ( t )eie dt
mr(t) t - c l t+&ib
++ce
1
21
m

-0
-ie '*( t )dt -t 6+ - J - a. - - t
a-
~'(t) ?-cl 2m' -0
ëie dt
a'(?)t-cl t - 6 - i b

La fonction --
' ( t ) t
- g,
m'(O f
est la valeur limite de Ia fonction
'
--
@"(O 5
@'K)c-
holomorphe dans le demi-plan supérieur puisque Q'*(r) et n'(O le sont dans le demi-plan
inférieur, et w ' ( 0 ne s'annule pas dans ce demi-plan. 5; est également un point du demi-plan
inf6rieur. Cette fonction tend vers zéro lorsque +=. La première intégrale dans l'équation
(4.6) est donc nulle:

En poursuivant le même raisonnement ont obtient pour les intégrales suivantes:

 La seconde intégrale dans l'équation (4.6) est du type Cauchy, -'1 -f@) dt, OC f(t)
2m-t-rb,
est une fonction holomorphe dans le demi-plan supérieur, y compris à l'infini. D'après les
propriétés des intégrales de Cauchy (2.68)on a:

Comme la fonction
-
f ( t )= --
' t )
g,
af(t) t
1
-

l'intégrale prend la valeur

Parce que les translations t + 6 et t - 6 ne changent pas le type de Sintégrale on peut

poursuivre avec le même raisonnement qui nous conduit aux résultats:
 Si l'on combine les équations (4.6). (4.7), (4.8), (4.9), (4.12), (4. I3), (4.14) on
obtient:

Pour la deuxième intégraie dans i'équation (4.2)on a conformément à (2.75)

L'intégrande est la valeur limite de la fonction (2.71)

Le premier terme est holomorphe dans le demi-plan supérieur et son intégrale est
nue:

(puisque w(to) et (3i(to)sont holomorphes dans le demi-plan supérieur)

Le deuxième terme dans (4.17) est:

La fonction f (c) = ---

w
1
3 c-c*
1
c-to
1
est holomorphe dans le demi-plan supérieur, y

compris à I'inf'i. Avec a>(t),la transformation conforme pour les trois dents (4.4), donne la
valeur de 4:
 Comme nous montre (2.54)

et donc:

Si l'on combine (4.20),(4.21), (4.22)et (4.23) on obtient:

On additionne (4.15), (4.24) et l'on obtient la valeur de @*(c)

c,
On remplace Cl par car C est un point du demi-plan inférieur.
Ainsi, la valeur de a(6) pour les trois dents est:

-
iû a W(ib- 6 ) - e a w ( i b+ 6 )
-e
c - (ib- 6) Z(ib- 6) r - (ib+ 6)g ( i b + 6 )
 La dérivée de la transformation conforme (4.4)est:

(2 .e-if'
a - a - eiû
m'(O = c(1+ 2cos 8 )-
(c-ib12 [(~+s)-ib]~-[(~-6)-ib]~

4.2. Explicitation du potentiel @(clDOW la transformation conforme compIète

D'après l'équation (4.26):

d'où

et I'on a 3n équations complexes:

 -
n
ak @'(ibk) -
- 6 )= Ob(-ibj - 6 )+ Ç
(-ibj - 8 ibk12 x ( i b k )
n - IL
- (4.3 lb)
-.i@ C, ak . -6 ) + e-ie ak Q'(ibk + 6 )
k=l (bj + bk l2 x(&k 6) k=i( -ib i -ibk -2q2 ' Z(ibk+6)

n -
a'(-ibj + 6)= 06 (-$ + 6)+ ak Wbk) +

1 -ib j +6-ibk)? ' =(ibk)

- n - (4.3 lc)
W(ibk+ 6 )
+eie ak
+akbk)2.W(ibk+6 )

D'après l'équation (4.28) on obtient:

m'(O = c(1+2 cos 8 )- c(C - -E[(c+ 6)-

n

l=i
"1

ibl)2 1=1
aleie

1
ib~
2 -z[(c
1 ~ 1
ap-"
- 6)- ibl l2
(4.32)

Comme o'(-ibk ) = Z ( i b k ) (4.33)

on calcule o'(-ibk) en tenant compte des propriétés des nombres complexes:

wr(-ibk) = c(1+ 2 cos 8 )+

n
a1
n

2+z1=1( bk + q - i b p
l=*(bk+bl)
o,ëi8
+zl=1(&+bl+i6)2

On exprime les potentiels W(ibk) et W(-ibk) par les équations (2.1 1O), (2.1 11)
et de la même façon:

a'(-ibk - 6)= XIk +i n k

-
- 8)= XIk
Q)'(ib& in

@'(-ib,+ 6)= X2, + iY2,

-
W(ibk+ 6 )= X2, - iYZk

En suivant les étapes de calcul pr6sentées en annexe 1(CALCUL DES CONSTANTES

&EUS) on obtient le système de 6n équations réelles:

Ce système d'équations étant résolu, on obtient bs valeurs des constantes réelles Xk,
Y,, X I k , YIk, X2& et YZk.Cela nous permet de déterminer @'(O
4.3. Calcul des contraintes su~erficielIq

Les contraintes superficielks sont situées à la sUTface des dents, donc pour les points à
Ia frontière du domaine

LZk(t-~+ibk)2X2k-iY2n
+eie
k=i [(t f 6)2+ bz] x ( i b k + 6)

Pour simplifier I'écriture, nous allons noter:

W ( t )= El ( t )+ iEi ( t )
avec

n
w 1
El ( t )= --- 2lrt-t* k=l
*2Ykbk?+~k(t2-bz)+
( t 2 + b;)2

'
k=~
z ( i b k-6 ) [(t + 612+ @]2

[(t + 612+ b:]2

 +'k=l
IL
=k
- 6)
w(ibk
kn,&(t+8)+xik[(t+ s>Z-b:]}sin
+ m2 +$12
[(t
e
+

D'après l'équation (4.28) on obtient:

Nous d o n s expliciter chaque terme de l'équation (4.44):

=k " ak[(t+ 612- bz + 2ibk(t+ 6)l

&=IN' + 6)- i b k 12 =
t=*z [(t+6)2+b:]2

"
k,J('-6)-&kl
ak
2 =z
,=l
-b: d i b k ( t - 6)]
tzk[(t-qZ
[(t - o ) +~ b;l2
On sait que:

eie = cos f
3 +isin 0
= cos 0 - isin 8
ëie
 Si l'on combine (4.49),(4.50),(4.51),(4.52), (4.53) et (4.54)on obtient:

-icKt
n

L=I
2ukW + al
+ a2+ $l2
0 - i~ a ~ [ ( t + 6 ) 2 - h ~ s h e"+ z2akbk(t+6> sin 8 -
n

bt
k=i Kt + 6)2 + l2
2

k,I Kt + 6)2 + 4 I

Pour simplifier I'écriture, nous allons écrire:

2akbk(t + 6 )
1 cos 8 -
(4.58)

Comme nous l'avons déjà mentionné en (2.126)dans la contrainte superficielle on a

besoin de:

La valeur de la contrainte superfcielle est donnée par:

 Si l'on combine (4.59) et (4.60) on obtient:

4.4. L'influence des dents adiacentes sur Ia contrainte su~eficielledans Ia dent chareée

Pour connaître l'influence des dents adjacentes sur la dent centrale. dent chargée. nous
avons programm6 la contrainte superficielle. équation (4.61). On a étudié le profù de la dent
d'engrenage droit standard, correspondant à:
- nombre de dents N=20 dents
- angle de pression 9=200
- pas diamétral P=l
- rayon standard de l'outil r=û.38/P
Contrainte superficielle SIGMAtO= 3.933,SIGlvlAt= 1-248

Figure 4.1. Contraintes superficieIIes dans un engrenage externe

 Courbe 1:représente les trois dents adjacentes
Courbe 2: représente la distribution des contraintes pour trois dents
Courbe 3: représente la distn'buton des contraintes pour une dent
L'échelle est en pouces.

La charge est appliquée sur le cercle de tête de la dent centrale, c'est-à-dire to=.Uû35.
Les coefficients de la transformation conforme que nous avons utilisés sont les
coefficients déteminés dans le chapitre III, c'est à dire:
. .

1 dent [8] 3 dents [chapitre m]

al 0.10074 0.09999 I

a2 0.26740 0.27524
a3 0.0121 1 0.01181

Le programme que nous avons utilisé s'appelle SIGMAN3.M et se trouve dans

l'annexe 3.
Si l'on analyse les résultats, on constate une grande différence entre les valeurs des
contraintes. Pour le cas d'une dent isolée, l'effort maximal a la valeur de 3.933 psi, en
comparaison avec 3.867 psi obtenue par Cardou et Tordion [8]. AGMA indique une valeur de
3.933 psi [8]. Si l'on considère les dents adjacentes, la valeur des contraintes est de 1.248 psi
(voir explications à la page suivante).

4.5. Calcul du facteur néornétrique

La formule de Lewis donne la contrainte superficieiie maximale pour un facteur de

forme dom& J.
où W,est la force tangentielle entre les engrenages
F est la largeur de la dent.

Donc, 3
F
= W - cos @ ,puisque W est la force appliqu6e par unit6 d'épaisseur.
La dation avec la contrainte caicuiée est:

où h, est la hauteur de dent étudide; dans l'exemple étudié &=2.2SW.

Dans le système américain, si W = l I b h , on a:

La charge étant appliquée sur le cercle de tête, nous avons trouvé que la contrainte
maximale de tension est de 1.248 psi pour W=l Iblm. On trouve donc comme facteur de
forme pour le côté en tension:

cos #
= -= 0.753
-

Le facteur geométrique obtenu quand on considère, dans la transformation conforme,

seulement la dent chargée est 6gd à:

- = 0.243 s'il est calculé par MPC [8]

- Jtmion= O. 239 s'il est calculé par AGMA 183

Si l'on &tudiela distribution des contraintes pour le cas des bois dents adjacentes, en
comparaison avec ceile pour une dent isolée (figure 4.I), on remarque une différence entre les
deux courbes. Elles gardent le même profil mais sont comme translatées une par rapport à
l'autre. On sait qu'au point d'application de la charge on a un point singulier dans la fonction
de potentiel @(t), qui intervient dans le calcul des contraintes superficielles. C'est justement
autour de ce point qu'on note le comportement différent des contraintes. On obtient pour trois
dents une zone de compression, sur le flanc de la dent, à la place d'une tension. Au pied de la
dent on a une tension mais sa valeur est trois fois plus petite que celle pour une dent, et cela
est du probablement a un terme singulier qui semble manquer dans le calcul des contraintes.
C'est cette même valeur qu'on considère dans le calcul du facteur géométrique, et le résultat
ainsi obtenu est aussi trois fois plus grand.
Les grandeurs géométriques employées pour le cas étudie se trouvent aussi dans le
premier tableau de I'annexe 2 .
 APPLICATION DE LA MI?IMODE DES POTENTELS COMPLEXES AUX
ENGRENAGES m m s

Le calcul des contraintes dans les engrenages internes est en générai détermini
comme pour les engrenages externes. Mais, comme la géométrie d'une dent d'engrenage
interne est différente de celle d'une dent d'engrenage externe, Mat de contrainte nécessite un
calcul adéquat. Nous présentons dans ce chapitre une méthode analytique pour la génération
dime dent d'engrenage interne droit et I'application de la Méthode des Potentiels Complexes
pour ce type d'engrenage.

5.1. Définition d'un ensenage interne;

Dans les chapitres précédents on a présenté l'engrenage externe, pour lequel les dents
se trouvent 2 l'extérieur de la surface d'un cylindre (pour les engrenages droits) ou d'un cône
(pour les engrenages coniques). L'engrenage interne est le complément, ou le négatif de
Ifengrenageexterne, car ses dents se trouvent à I'intérieur de la surface d'un cylindre ou d'un
cône, donc elles sont orientées vers le centre.

Figure 5.1. Parties d'une dent d'engrenage interne

 Parmi les caractéristiques d'un engrenage inteme on cite:

- pour le même rapport de transmission, la ligne d'action de l'engrenage

interne est pIus grande que pour l'engrenage externe;
- les dents sont plus solides en comparaison avec celles d'un engrenage
externe sindaire;
- I'engrenage interne forne sa propre protection autour des dents;
- c'est l'engrenagele plus recommandé pour certains types de machines du
point de vue sécurité;
- l'entraxe est égal à la différence entre les rayons primitifs du pignon et de Ia
roue;
- comme l'entraxe d'un engrenage interne avec le pignon conjugué est plus
petit que pour un engrenage externe similaire, I'aspect de I'ensemble est
plus compact;
- I'engrenage interne élimine l'utilisation d'un engrenage intermédiaire quand
on a besoin d'avoir deux arbres parallèles avec le même sens de rotation;
- I'engrenage interne ne peut pas être utilisé quand le nombre de dents du
pignon est proche du nombre de dents de la roue;
- un engrenage interne ne peut engrener qu'avec un pignon dont la denture est
externe. La possibilité d'interférence est donc plus grande que dans le cas
des engrenages externes;
- l'engrenage interne a comme inconvénient qu'il y a peu de machines-outils
qui peuvent le fabriquer.

Génération d'une dent à dévelo~~ante

de cercle

Comme dans Ie cas des engrenages externes, le prof2 le plus répandu est celui en
développante de cercle. Le profil est le même, mais la forme de la dent de I'engrenage
interne, est celle de l'espace entre les dents pour les engrenages externes. Ainsi, le contact se
fait à l'int&ieur, ou dans la partie concave, de la développante, au lieu de se faire dans la
partie convexe de la développante, comme dans i'engrenage externe. De plus le rayon de pied
est le rayon le plus grand, et les sommets des dents se trouvent sur le petit rayon.
 Figure 5.2. Comparaison d'un engrenage externe et d'un engrenage interne

5.2-1. Dimensions de génération

Figure 5.3. Géométrie d'une dent d'engrenage interne

Pour calculer les dimensions de genération, nous devons connaître certains

paramètres de génhtion:
 Z nombre de dents
m module (en millimètres)
qc anangle de pression
a addendm
b dedendum
r, rayon d'arrondi de l'outil
x facteur de déport de l'outil
Le problème traite dans ce chapitre considère le facteur de déport de Poutil x=û.

Pour les engrenages cylindriques aux axes parallèles ou gauches, AGMA a établi des
normes pour un profil normalisé.
Addendum (sailIie ):

Dedendum (creux):
b = -1.25
P
La hauteur de la dent:

Ainsi un engrenage normalisé aura les dimensions suivantes:

Rayon primitif:
rl

Rayon interne:

Rayon externe:
Re = R+b
Rayon de base:
l$, = R - C O S # ~
Pas diamétral:

où rn est toujours en rnillimiitres.

5.2.2. Génération du orofil acif

Le profil actif de la dent est la partie de celle-ci qui &srnet le mouvement de rotation
à i'engrenage. ou q o i t ce mouvement d'un autre engrenage.

Pour générer le prof3 actif d'un engrenage interne, il faut suivre les étapes suivantes
Pl:
1) Choisis le nombre de points N. pour tracer le profil;
2) Calculer le rayon ri correspondant en n'importe quel point du profil:

Figure 5.4. Génération du profil actif (fi3we tirée de [3])

3) Calculer l'angle de pression #,,correspondant à n'importe quel point du profil:

 & = arccos-Rb
ri

4) Wcuier l'angle polaire, à n'importe quel point du profil:

5) Calculer l'épaisseur circulaire T sur le cercle primitif:

6)Calculer 1'6paisseurcirculaire Ti sur n'importe quel autre cercle:

7) Calculer l'angle 19" formé avec l'axe de symétrie de la dent:

8) cdcder les coordonnées cadsiennes du proffi:

- I'abscisse du profd X :

- I'ordomee du profil Y:

5.2.3. Génération du profil de raccord

L'engrenage interne est d'habitude généré avec un couteau-pignon. La forme du

pignon suit, dans ce cas, le proN de raccord de la dent d'engrenage.
 r, = le rayon de raccord
8, = I'angle vectoriel de raccord
E, = l'angle de rotation de I'engrenage
interne
C = l'entraxe
Rm = le rayon exteme du couteau
E, = I'angie de rotation du couteau

Figure 5.5. Génération du profil de raccord

On analyse comme profd de raccord, une trochoide, et pour la générer il faut suivre
les étapes suivantes (pour le cas r, = 0):

1) Choisir le couteau avec lequel on va générer la trochoïde. Cela signifie qu'on

connaît ses caractéristiques:
Zc:nombres des dents du couteau
rn :module
$c :angle de pression
et qu'on peut calculer les dimensions du couteau avec les formules (5.I à 5.8)

2) Calculer I'entraxe C :

3) Calculer I'angle de rotation du couteau E, pour un rayon r, calculé:

 4) Calculer I'angIe de rotation de I'engrenage interne E, :

5) Calculer l'angle de la trochoïde 8, :

6) Calculer l'angle de la tmchoïde, par rapport à une ligne choisie comme référence
If

9, :

Si Iton utilise comme ligne de reférence l'axe de symttrie de la dent, 6 se calcule

avec la formule suivante:

avec

où Tcest l'épaisseur circulaire du couteau sur le cercle primitif, caicalée avec (5.12)
Pour obtenir @==, l'angle de pression au sommet du couteau, on utilise (5.10)

L'angle polaire se calcule avec (5.11)

 7) Calculer les coordonnées cartésiennes du prohl de raccord:
- l'abscisse de la trochoïde X,:

- I'ordomée de la trochoïde Y, :

Pour tracer les deux profils d'une dent d'engrenage interne, Ie profil actif et celui de
raccord, on peut utiliser le programme lhVI3WE.M (annexe 3).

5 .2 -4. Assemblage des profds

On peut assembler le profil actif et le prom de raccord, pour former une dent
complète, en trouvant graphiquement ou analytiquement leur point d'intersection (X,,Y,).

Méthode graphique: La méthode graphique est très simpie à employer, il s'agit de tracer les
deux profds dans le même plan et de localiser " à l'oeil " le point
d'intersection (X,. Yf) .

Méthode analytique: Cette méthode consiste à calculer le rayon du point où les deux profils se
rencontrent. Pour cela on utilise la formule:

Une fois le rayon d'intersection des deux profils calculk, on détermine les
coordonnées cartésiennes de celui-ci:

Pour calculer l'angle O;', formé avec l'axe de symétrie de la dent, on utilise les
formules (5.10) a (5.14)- pour le rayon Rf .
 Figure 5.6. Assemblage des profils

L'exemple présenté est pour le cas concret d'un engrenage interne avec 2=45 dents,
module m=4, angle de pression &=200, et un couteau avec Zc=17 dents. Le programme
utilisé est RACC0RD.M (annexe 3). Les figures 5.6 à 5.8,s. 10,512 et 5.13 sont en mm.

5 -2.5. Influence du t p e de couteau sur le profl de raccord

Pour faire cette analyse on utilise le programme C0UTEAU.M (annexe 3).

Avec le même engrenage, présenté dans l'exemple précèdent, et un deuxième couteau
avec Zcl=30dents, on obtient la figure 5.7 où:

Courbe la: représente la partie trochoïde d'une dent générée par un couteau à 30
dents
Courbe lb: représente la partie trochoïde d'une dent générée par un couteau à 17
dents
Courbe 2a: représente la partie développante d'une dent gtinérée par un couteau à 30
dents
Courbe 2b: représente la partie développante d'une dent générée par un couteau à 17
dents
 On observe que les couteaux déterminent des trajectoires différentes pour les parties
trochoïdes, mais cette différence n'affecte pas le profil de la dent car elle est en dehors du
profil de la dent. La dent garde la même représentation entre le point d'intersection des deux
parties (la partie trochoïde et Ia partie développante) et le rayon externe de la dent pour les
deux couteaux.

INFLUENCE DU COüTEAU SUR LE RACCORD

-85 rn a I I a 1 I I

Figure 5.7. L'infiuence du couteau sur le profil de raccord

5 -2.6. Com~araisondu orom complet avec le profil déterminé ~ aune

r autre méthode

On poursuit avec le même exemple d'engrenage interne de 2545 dents, et couteau de

ZC=17 dents. Le profil complet, généré en suivant le cheminement présenté dans les
paragraphes 5.2.2.. 5.2.3., et 5.2.4., a 6té comparé avec les coordonnées obtenues par une
autre mtthode [16]. L'exactitude des n5sultats est visible sur la figure 5.8, où le profil "-"
est le profil MPC et celui trac6 avec la mtthode[la est représenté par "+".
La figure 5.8. est tracée avec le programme INTERNE-M (annexe 3).
 ENGRENAGE INTERNE-PROFIL REEL

"r

Figure 5.8. Profil complet d'une dent d'engrenage interne

Transformation d'un enmen= externe dans un ewenaee interne

Travaux publiés

Suite à Ieun études, Von Eiff et al. [36], en 1990 ont montré que ia géométrie de la
dect change de façon continue quand on passe d'un type d'engrenage à I'autre.

pigmn type coumu

*<O

Figure 5.9. La transition de I'engrenage interne vers I'engrenage externe (figure adaptée
de 1361)
 De plus, ils ont montré que les engrenages internes et externes peuvent être décrits
par les xnêmes équations si I'on applique la règIe suivante: on considère un nombre de dents
positif dans le cas des engrenages externes et un nombre de dents négatif dans le cas des
engrenages internes; la c r é 4 2 r e a une infinité de dents.

5.3 -2. L'application de la Méthode des Potentiels Complexes au cas des enmena~esinternes

La Méthode des Potentiels Complexes peut aussi être utilisée dans le cas des
engrenages internes. Pou. cela on utilise encore la transformation propos6e par Hirano,
(équation 2.27):

Cette trausfonnatioo ramène la profil symétrique par rapport à l'axe des imaginaires
(dans le Plan (2)) à l'axe réel dans le Plan ( 6 ) .La condition qui s'impose, pour avoir une
parfaite modélisation de la dent interne, est un nombre plus grand des termes de la
transformation conforme. La transition de l'engrenage externe vers l'engrenage interne est
analysée avec le programme MSIT'M (annexe 3).

TRANSmON D UNE DENT EXTERNE VERS UNE INTERNE

Figure 5.10. L'influence du nombre de termes de la transformation conforme sur le

profil de la dent
 Courbe 1: correspond à une transformation conforme avec quatre termes; le
profil tracé est le profil d'une denture droite exteme.
Courbe 2: correspond à une transformation conforme avec cinq termes; on
observe que la développante de cercle commence changer de
concavité.
Courbe 3: correspond à une transformation conforme avec six termes.
Courbe 4: correspond à une transformation conforme avec sept termes; le profil
îracé est le profil d'une denture droite interne.

Comme on I'a déjà mentionné, la Méthode des Potentiels Complexes a aussi été
utilistk dans le cas des engrenages internes par Miyachika et Oda [23]. Ils ont aussi étudié
l'effet de l'épaisseur de la jante d'un engrenage interne sur les contraintes et la dùection de
propagation de la fissure sur la dent La transformation conforme qu'ils utilisent ûansforme
un domaine consistant en un plan semi-infini sur lequel repose une courbe symétrique par
rapport à l'axe des réels (Plan ( 2 ) ) dans l'axe imaginaire du plan semi-Mini (Plan (c)). La
forme de la transformation conforme proposée par Miyachika et Oda [23] est la suivante:

Comme on a voulu faire une comparaison entre l'équation (2.27), et celle proposée
par Miyachika et Oda, equation (5.25), on a ramené celle-ci dans le même plan. Si l'on
exprime la relation (5.25) pour chacune des coordonnées x et y on obtient:
 Maintenant il faut montrer que la transionnation proposée par Miyachika et Oda
(5.25) transforme l'espace semi-innni born6 d'une courbe symétrique par rapport à i'axe des
réels dans le demi-plan correspondant (54).

Figure 5.11. Transformation des domaines proposée par Miyachika et Oda

Si l'on prend la position de l'origine dans le plan z , c'est-à-dire le point de

coordonnées x = O , y = 0 , les équations (5.26) deviennent:

b
Dans Ia 2ème équation du systérne (Ç +a)* + q2 = --. En remplaçant ce résultat
C
a
dans la première équation du systéme, on obtient 5 = --. Comme la valeur du coefficient a
2
est toujours positive, le point correspondant à l'origine dans le plan z se trouve dans le demi-
plan (54) dans le plan 6. Nous avons utilisé un seul terme dans la transformation, étant
suffisant pour la démonstration.
 Puisque,les axes sont inversés, Ia frontière du demi-plan correspond à 5 = O. Donc,
l'équation du prof2 correspondant dans le plan(z) est donnée par:

On constate que la transformation (5.25) est la même que l'équation (2.27), avec les
axes inversés, c'est-à-dire x devient y et y devient x. Pour garder la même notation de
l'équation (2.27) nous avons chan@ la notation des coefficients ak et bk entre eux.

5.3.2.1. Coefficients de la transformation conforme;

Comme nous l'avons déjà mentionné, la transformation de Hirano (2.27) transforme

I'axe des réels du plan en une courbe symétrique par rapport à I'axe des imaginaires et
asymptotique à l'axe des réels. La forme de la courbe, ainsi obtenue, dépend des coefficients
c , ak, 4 , de la transformation conforme. La figure 5.10 nous montre que pour un
engrenage interne le nombre des termes dans la transformation conforme devrait être d'au
moins 7, cela implique la recherche de 15coefficients.
On powuit avec le même exemple d'engrenage interne de 2=45 dents, module m=4,
angle de pression 4, = 20°, mais cette foisci avec un rayon de l'arrondi de I'outil calculé
r,=1.4.
La transformation conforme correspondant à cet engrenage a les paramètres suivants

Pour le cas du profil de raccord, où I'arrondi de l'outil est non-nul, en se basant sur
la géometrie de la dent représentée dans les figures 5.4 et 5.5, il faut modifier les étapes
présentées dans 5.2.2. et 5.2.3, comme suit:
1) Choisir une valeur exptnmentale pour le rayon au sommet du profi en
déveIoppante pour le couteau IZ,
2) Calculer l'angle de pression, @ correspondant & (équation 5.10)

3) Calculer l'épaisseur T- du couteau à Ip,avec une formule derivée de (5.13)

4) Calculer Ie rayon de I'arrondi de PoutiI r,

r, =
2 cos, #

5) Calculer l'angle de pression #& à R,, rayon de I'arrondi au sommet du couteau

6)Calculer R, avec une formule dérivée de (5.7)

7) Calculer le rayon externe du couteau R

,

8) On suit les ttapes présent6es dans 5.2.3.

Calculer I'angle de rotation du couteau E, pour un rayon r, calculé ( r , prend des
valeurs entre ( R, + C) et Re, équation (5.18):
Calculer I'angIe de la trochoïde 8,, 6quation (5.20)

0, = arcsin
R, sin E,
--.R, €Ca
'=ta R

9) Calculer la variation de Tangle de la trochoïde par rapport à l'angle de rotation du

couteau:

10) Calculer la variation du rayon de I'arrondi par rapport à l'angle de rotation du

couteau:

\J*dJ)

rm

1 1) Calculer l'angle entre le rayon et la tangente à la trochoïde:

arctan

12) Calculer le rayon pour le nouveau profil de raccord Rm:

13) Calculer l'angle du nouveau profil de raccord:

Of = 0, +.rccos(
ru - rrsin yt
I* )
14) Calculer I'angle du nouveau profil par rapport à une axe de symétrie choisie
comme référence; si Pon utilise comme ligne de référence I'axe de sym6trie de la dent,
on obtient:
 15) CaIcuIer les coordonnées cartésiennes du profil de raccord:
- l'abscisse de la nouvelle trochoïde XT,:

- i'ordomée de la nouvelle trochoïde K :

Pour le profil actif de la dent, partie dbveloppante de cercle on suit les étapes
présentées en 5.2.2.

ENGRENAGE INTERNE, PROFIL REEL,PROFIL MPC='

T w 1 .

Figure 5.12. Engrenage interne. profil réel et profü tracé par MPC "-"

On constate facilement que la Méthode des Potentiels Complexes modélise

parfaitement la dent d'engrenage interne droit, comme elle le fait aussi pour l'engrenage
externe (figure 5.12).
5.3-2.2. Calcul des contraintes mur les enpnapes intemq;

Pour analyser les contraintes dans les engrenages internes, nous aiions utiliser la
Méthode des Potentiels Complexes. Dans 2.6.4.. nous avons reformulé les équations de
Kolossov dans le plan 5. Les équations

nous permettent de calculer les contraintes. Il faut exprimer les potentiels complexes dans le
plan 6 tenant compte de la transformation conforme utilisée. La force est appliquée à la
position &.

Le calcul des contraintes superficielles suit le cheminement indiqué dans 2.6.4.2. Le

cas traité est encore celui de l'engrenage interne étudié jusqu'ici: -5 dents, m=4,
#== 20°, et une charge PnlF=9.8 N/mm (F=largeur de la dent, 9.8 étant pour la
transformation du kgf en unités SI) appliquée sur le cercle de tête. Le programme
S1GMAN.M (annexe 3) permet de calculer la contrainte superficielle, d'abord le long de la
face de la dent (partie involute), puis dans la région critique du pied (partie de raccord). La
distribution des contraintes sur la face droite (contrainte en tension) et sur la face gauche
(contrainte en compression) de la dent est présentée sur la figure 5.13. Le signe "+" indique
la contrainte en tension et le signe "-" indique la contrainte en compression.
La distribution des contraintes démontre que la contrainte maximale se trouve dans la
partie trochoïde du pronl. On voit que la contrainte en compression est supérieure à celie en
tension. Si l'on fait une comparaison avec l'engrenage externe, l'effort au long du flanc de la
dent, partie involute, est plus élevé dans Ies engrenages internes que ceux externes.
Cependant la contrainte dans la section du raccord reste toujours plus grande qu'au niveau de
l'involute. On a obtenu pour la contrainte en tension a~a= 7.05 MPa et pour la contrainte
en compression une valeur de cc = -8.3 18 MPa.
I
I(

Si l'on regarde la connainte en tension autour du point d'application de la charge,

e0 = 0.158, on observe qu'elle prend une valeur négative. Cela veut dire que la dent subit
une compression proximité de la force appliquée (plus précistment jusqu'à 6 = 0.23).
 Contrainte supeficielle SIGMAc= 8.318. S I G W 7.05
12 1 B i
1 B 1 T

Figure 5.13. Contrainte supeficielle dans un engrenage interne

5.3.2.3. Sectioncritiaue;

La position de la section critique, dans une dent d'engrenage, est définie comme étant
la distance radiale entre le cercle primitif et le cercle passant par le point de contrainte
maximale. On peut évidemment d é f une ~ section critique à partir du côté tension et une
autre à partir du cet6 compression. On trouve que les deux côtés sont toujours tres proches.
sinon identiques (à la précision des calculs). La section critique, pour le cas étudié, où 1'0x1 a
trouvé les contraintes maximales. se trouve à 6 = 22.9.

Il existe deux méthodes approximatives pour la calculer l'une basée sur la parabole
de Lewis. et i'autre basée sur la théorie de la tangente. Dans les normes AGMA, la section
critique e s situte au niveau du point d'intersection de la parabole et du profil de raccord.
Hofer remplace la parabole par une tangente au profil de raccord. Les normes DIN (nonnes
allemandes) recommandent pour I'angle entre cette tangente et l'axe de symétrie de la dent
une vdeur de 300.
 a) Lewis b) tangente de Hoftr

Figure 5.14. Section critique dans un engrenage externe

Dans le cas des engrenages internes, pour une force appliquke sur le cercle de tête, la
parabole de Lewis rencontre Ifinvolute mais pas le profil de raccord. En appliquant la
méthode des tangentes à un engrenage interne, l'angle entre I'axe de s y d t n e de la dent et la
tangente au profil de raccord est plus grand que 300. Il y a peu de publications dans
lesquelles les engrenages internes sont étudiés [36].Dans ces travaux Ia section critique se
trouve pour des valeurs de I'angle entre 420 et 590. b
s angles recommandés sont de 4S0 et
500.

a) Lewis b) tangente de Hofer

Figure 5.15. Section critique dans un engrenage interne

 Dans le cas que nous avons étudié, on a considére les recommandations faites par
[23] et l'on a pris a=450. Ainsi les valeurs des contraintes dans la section critique
recommandée, calculées avec S1GUAN.M correspondent à' = 22-4 et sont: crr(a= 45")
=7.0432 MPa, q ( a= 45") = -8.3076 MPa

5.3.2.4. Formules ~ratiauesaour calculer les contraintes su~eficiellesdans les ewenaees

intemts

Figure 5.16. Effort nominai dans une dent d'engrenage interne (figure tirée de [23])

Les dimensions que l'on voit sont calculées à partir de la théorie développée dans (5.2.2):

s=10.28 mm y=1.827 mm

r,= 1.4 mm 0=700

L'angle O est choisi d'après les recommandations trouvées dans 1101, [23] et [36].
 Une force normale à la surface de la dent se décompose en composantes, normale
PR-sinO et paraU2Ie P, cos@ à l'axe de symétrie de la dent Les contraintes dans la section
critique sont données par:
- pour le côté tension

- pour le côté compression

Pn sin O
TN =
S-b

Pou.une dent d'engrenage interne droit, Miyachika et Oda 1231, ont déterminé des
équations pour faire le calcul des contraintes superficielles.
La contrainte en tension est calculée:

~,=(1+0.046s/r,)(O.67~,,+0.48JaZ mb +36<+1.140,,) (5.43)

Pour le côté en compression, la contrainte se calcule en fonction du rapport Ys

 LRS calculs faits avec les f o d e s présentées ont donné les valeurs suivantes:

- pour le côté en tension, avec (5.43) a,=7.29 MPa

- pour le c6té en compression avec (5.44) a, =-8.4855 MPa

Si Iton fait une comparaison entre les valeurs obtenues avec la Méthode des Potentiels
Complexes et celies obtenues à l'aide des calcds que nous vous avons présentés on observe,
pour un engrenage de W 5 dents et une force normale à la surface appliquée sur le cercle de
tete:
CONTRAINTES SUPERFICIELLES (MPa)
L ,
tension compression
W C ( a-) 7-05 -8.318
AGMA (o(a= 4 5 O ) ) 7.0432 -8.3076
Y 4 5 9 /,
~ ( C= ,
a 0.999 0.999
1 1

Formules pratiques ( a) 7-29 -8.4855 ,

. Cr/ 0- 1.O34 1-02

Analysons les résultats obtenus. La Méthode des Potentiels Complexes donne des
valeurs très proches pour la valeur maximale des contraintes et la valeur dans la section
critique indiquée par AGMA. Les rapports des valeurs de 0.999 pour le côté tension ainsi
que pour celui en compression en sont la preuve. Les résultats obtenus avec les formules
pratiques mentionnées donnent un écart de 3,4% pour le côté tension et de 2% pour le côté
compression. On peut donc remarquer pour le cas étudié une bonne corrélation des résultats
obtenus par la Méthode des Potentiels Complexes et les valeurs prédites par 1'AGM.A. De
plus les f o d e s pratiques proposées par [23] semblent passabIement exactes.

Il faut mentionner que la valeur maximale réelle de la contrainte dépend beaucoup de

I'épaisseur de la jante. La contrainte superficielle augmente du côté en tension autant que du
côté en compression, avec la diminution de I'épaisseur de la jante. La charge appliquée sur
une dent se fait sentir sur les dents non chargées situées d'une part et d'autre. Une jante plus
mince se déforme plus rapidement car Mat de contrainte est plus élevé [23]. Si I'bpaisseur
est plus petite que cinq modules (Sm),on doit tenir compte d'eile dans le calcul des
contraintes, en utilisant un facteur de concentration de connainte [36].Cependant, pour une
6paisseu.rsupérieure à Sm Mat de contrainte devient indkpendant de celle-ci et i'on peut donc
calculer l'état de contrainte à l'ai& de la Méthode des Potentiels Complexes. On notera que la
nonne ISO permet cette application pour une épaisseur supérieure 2 3.375rn.
On peut évaluer la contrainte -4edans les engrenages droits a jante mince avec
des formules simples, qui tiennent compte du module de la denture et de l'épaisseur de la
jante [2q.

On sait que la f o d e de Lewis donne la contrainte superficieue maximale pour un

facteur de forme donné, J:

où Wt est la force tangentielle entre les engrenages

F est la largeur de la dent.
P est le pas diamétral de la denture
J est le facteur géométrique.

Donc, wt = W
-
F
.
cos @ puisque W est la force appliquée par unité d'épaisseur. La
relation avec la contrainte calculée est:

Pour l'exemple que nous avons étudié, W = -=9.8 N/-, -5 dents, P =O -25.
F
# =200,on a obtenu le facteur géométrique:

pour le côt6 tension JICIUiO,,

=O.326
pour le côté compression Jcoq-ba =0276.

Pour faire la comparaison du facteur géométrique calcul6 avec les valeurs indiquées
par les normes AGMA, nous n'avons pas trouvé des formules nous permettant de faire le
calcul.
 Cependant, dans ses notes de cours de Tordion [34, figure 1A, page 7 Ibis], propose
une formule empirique pour représenter la courbe AGMA donnant le facteur géométrique
pour un engrenage droit externe:

Nous avons utilisé la suggestion de Von Eiff et al. [36] et nous avons utilisé un
nombre de dents négatif dans (5.48). Le résultat ainsi obtenu est J=0.341. La différence
entre cette valeur et celle obtenue par la Méthode des Potentiels Complexes est de 4% et nous
permet de déduire que le changement que nous avons apport6 à la formule (5.48) est valable.

Buckingham 13, tableau 22.4, page 4761, indique quelques valeurs pour le facteur
géométrique, pour le cas d'un engrenage chargé par une force appliquée sur le cercle de tête.
Pour un engrenage de 2 4 3 dents, le plus proche de celui que nous avons étudié, J=0.200.
Ceci est assez loin de J=0,326. Cependant, on peut dire que la valeur donnée par
Buckingham est du côté sécuritaire car elle donne une valeur de contrainte supérieure.
Les grandeurs géométriques employées pour le cas étudié se trouvent aussi dans le
deuxième tableau de l'annexe 2 .
 Le travail a porté sur l'analyse des dents d'engrenage droit par la Méthode des
Potentiels Complexes. Une =vue des travaux déjà accomplis dans ce domaine a d'abord été
effectuée. De cette revue, nous avons détermin6 les objectifs du présent travail, soit I'étude
de I'innuence des dents adjacentes sur la dent centrale chargée et l'étude des engrenages
internes par la même méthode.

On a d'abord présenté de façon détaillée la Méthode des Potentiels Complexes ainsi

que son application aux engrenages droits. Comme cette méthode repose sur la possibilité
de passer par transformation conforme du profïi donné au demi-plan nous avons déterminé
la transformation conforme dequate aux trois dents adjacentes. Une fois la transformation
trouvée, nous avons fait la recherche des coefficients nécessaires pour la modélisation des
dents. En utilisant par la suite la transformation conforme dans le calcul des contraintes,
nous avons noté dans la distribution des contraintes, pour le cas des trois dents adjacentes,
un comportement différent autour du point d'application de la charge, par rapport à une
dent isolée. On sait qu'au point d'application de la charge on a un point singulier dans la
fonction de potentiel @(t), qui intervient dans le calcul des contraintes superficielles.
Comme c'est justement dans cette zone que nous trouvons les changements, ceci est
probablement la cause pour laquelle la valeur maximale est différente de celle pour le cas
d'une dent. Comme les résultats qui concernent le calcul des contraintes superficielles ne
sont pas satisfaisants, il faudrait éclaircir la cause de ce problème.

Nous avons ensuite étudié les engrenages internes par la même mbthode. Tout
d'abord nous avons présenté une m6thode analytique pour la ghération du profil interne.
Le p r o g r m e que nous avons conçu peut être utilisé pour déterminer la géométrie de
n'importe quel engrenage interne droit. Nous avons appliqué la Methode des Potentiels
Complexes dans le calcul des contraintes superficielles. La valeur des contraintes
maximales ainsi obtenue est presque identique avec celle obtenu avec la norme AGMA, et,
pour le cas éîudié, a une différence de 4% de la valeur obtenue avec des formules pratiques.
Aida, T., Terauchi, Y., "On the Bending Stress of a Spur Gear", Bulletin of the
Japanese Society of Mechanical Engineen, Vol 5, No. 17, 1962, pp. 161-170.

Baronet, C.,Tordion, G. V., "Exact Stress Distribution in Standard Gear Teeth and
Geometry Factor", Trans. of AShE, Journal of Engineering for Industry, Vol. 95,
Nov. 1973, pp. 1159-1163.

Buckingham, E., "Analyticd Mechanics of Gears", New York, Dover Publications,

1988.

Carathéodory, C., "Theory of Funcrions of a Complex Variable", Chelsea Publishing

Company, New York, 1954.

Cardou, A., et Tordion, G. V., "Calcul des Contraintes dans les Engrenages Droits par
la Méthode des Potentiels Complexes", Rapport technique~o.EM-25. Département de
Génie Mécanique, Université Laval, Québec, Canada, sept. 1980.

Cardou, A.. et Tordion, G. V., "Numericd Implementation of Complex Potentials for

Gear Tooth Stress Andysis", ASME Journal of Mechanicd Design, Vol. 103, No. 2,
apr. 1981, pp. 516-528.

Cardou, A.. et Tordion, G. V., "Calculation of Spur Gear Tooth Flexibility by the
Complex Potential Method", ASME Journal of Mechanisms, Transmission, and
Automation in Design, Vol. 107, mars 1985.

Cardou, A., et Tordion, G. V., "Étude du facteur géométrique dans les engrenages
droits", Rapport technique No. SM-85-02, Dtpartement de Génie Mécanique,
Université Laval, Quebec. Canada, juin 1985.

Camer, G.F., Krook, M., Pearson, C. E., "Fonctions of a Complex Variable", New
York, Mc Graw Hi11 Book Co, 1966.
10. Chong, T. H., Kubo, A., "Simple Stress Formulae for a Thin-Rimmed Spur Gear",
ASME Journal of Mechanisms, Transmission, and Automation in Design, Papers No
84-DET-62,84-DET63,84-DET-64,1984.

11. Churchill,R V., "Introduction to Complex Variables and Applications", Mc Graw Hill
Book Co, New York, 1960.

12. Daniewicz, S., "Conception and Development of Improved Analyticai Prediction

Models for Fatigue Induced Tooth Breakage Due to Cyclic Bending in Spur Gear
Teeth", Ohi0 State University, 199 1.

13. Dudley, D. W., "Gear Handbook", Mc Graw Hill Book Co, 1962.

14. England, A., "Complex Variable Methods in Elasticity",Wiley Interscience, New

York, 1971.

15. Gosselin, C., Gagnon, P., Cardou, A., Cloutier, L., "Usage de l'Interférométrie par
Niveau de Granulosités dans Ia Mesure de la Rigidité de Dents d'Engrenage Droits",
3eme Congres Mondial des Engrenages et des Transmissions, dec.14, Paris, 2, 1992,
pp.765-774.

16. Gosselin, C., "CoILLrnunication Personnelle1'Université Laval, Québec, Canada, 1995

17. Hirano, F., " Research on Stress Concentration Factors for Two-Dimensional Elastic
Bodies", Transactions of the Japan Society of Mechanical Engineers, Vol. 17, No. 61,
1951, pp. 12-16.

18. Hromadka IX, TT.V., Lai, C., "The Complex Variable Boundary Element Method in
Engineering Analysis", Springeflerlag, New York, 1987.

19. Kdandiya, A. I., "Mathematical Methods of Two-Dimensionai Elasticity" Mir,

Moscou, 1973.

20. Kauffmann, A., Douriaux, R., "LesFonctions de la Variable Complexe", Ed. Eyrolles,
Paris, 1962.
21. Lavrentiev, M., Chabat, B., "Méthodes de la TMorie des Fonctions d'une Variable
Complexen,Mir, Moscou, 1977.

22. Milne-Thomson, L. M., "Systèmes Élastiques Plans", Dunod, Paris, 1968.

23. Miyachika, K., Oda, S., "Bending Strength of interna1 Spur Gears", ISME
International Conference on Motion And Power Transmissions, Hiroshima, Japon,
nov. 23-26, 1991.

24. MuskheIishviIi, N. I., Torne Basic Problems of the Mathematical Theory of

Elasticity", Noordhoff Ltd, Groningen, HoIland, 1953.

25. NicoIetto, G., "Gear Tooth Stress Analysis by the Complex Potentials Method",
Meccanica, No. 27, 1992, pp. 105-1 10.

26. Nicoletto, G., "Approximate Stress Intensity Factors for Cracked Gear Teeth",
Engineering Fracture Mechanics VOL44, No. 2, 1993, pp. 23 1-242.

27. OueIIet, De,Guillot, M., Cloutier, L., Cardou, A., "Formulation pour le calcul de la
contrainte maximale dans les engrenages droits à jante mince", Trans. CSME, Vol. 17,
NO.4A, 1993,pp. 687-698.

28. Paré,D., "Application de I'Infographie a la Méthode des Potentiels CompIexes pour le

Calcul des Engrenages Droits", These de maîtrise,Département de Génie Mécanique,
Université Laval, Québec ,Canada, juin 1987.

29. Savin, G. N., "Stress concentration Around Holes", Pergamon Press, London, 1961.

30. Sneddon, 1.. Lowengmb, M., "Crack Problems in the Classicai Theory of Elasticity",
John Wiley & Sons, New York. 1969.

3 1. Sokolnikoff, 1. S., "Mathematical Theory of Elasticity", Mc Graw Hill Book Co, New
York, 1956.

32. Terauchi, Y., Nagarnura, K., "Study on Deflection of Gear Teeth", Bulletin of the
JSME,Vol. 23, No. 1û4, oct. 1980.
33. Timoshenko, S. P., "Theoryof Elastic Stability", M c Graw Hill Book Co, New York,
1961.

34. Tordion, G., V., "Él6ments de Machines 11 (Engrenages, Engrenages Planétaires,

Cames)",Notes de cours, D6parternent de Génie Mécanique, Université Laval, Québec,
Canada, 1978

35. Ugodohikov, A. G., Kwetsov, A. M., "On Static Stress Cdculation of Gear Wheel
Teeth", 2nd Scientific and Technicd htemationd Conference on Electric Modeling of
Problems of Structural Mechanics, Rostov-on-Don, jan. 28- feb. 2, 1962, pp. 348-
354.

36. Von ES, H., Hirschmann, K. H.,Lechner, G.,"Influence of Gear Tooth Geometry
on Tooth Stress of External and Intemal Gears", Journal of Mechanical Design,
Vol. 112, N0.4, dec. 1990, pp.575-583.

37. Wang, C.T.,"AppïiedElasticity",Mc Graw HU Book Co, New York, 1953.

 CALCUL DES CONSTANTES &EUS

Explicitation de @'(-ibj)

w
B~= - 1 (tOsin - bj COQ)
2r tt +b,
Donc pour @'(-$) on a 2n équations réelles:
Donc pour @'(-3
- 6) on a Zn équations réelles:

Ojk = G,,sin 9 + Gu COS 9

Fjk= G3 sin û + G32COS 9

Exnlicitation de @'(-i% + 6)

W 1
A2, =-
2~ - S12+ b, [(to- 6)COSB+bj sinp]
(tO
Donc pour W(-ibj + 6) on a 2n équations réelles:

En regroupant les équations réelles, on obtient le système suivant:

ANNEXE 2
 1 Nombre de dents 1 20 1
1 Angle de pression (O) 1 20 i
Addendum (po) 1 1

Dedendum (po) 1.25

Diamètre primitif (PO) 1 20
Diamètre de base (po) 18.79
1 Diamètre interne (po) 1 18 1
1 Dim&treexterne (DO) 1 22.5 1
1 Ravon standard de l'outil (DO) 1 0.38 1
Force appliquée sur le cercie de tête (lb/po) 1
Inciinaison de la force par r a ~ ~ oàr l'axe
t "x" (0') - 150.5

[ ÉTUDE DE CAS: ENGRENAGE EXTERNE DROIT - PARA~M-S 1

DE G & ~ ? R A ~ ~ O N

Engrenage
interne Pignon .
Nombre de dents 45 17
Module (mm) 4 4

Pas diamétral (mm+ 0.25 0.25

Addendurn (mm) 4 4
Dedendum (mm) 5 5 ,
Diarnhe arimitif (mm) 180 68
1 Diamètre de base (mm) 1 169.14 1 63.8 1
Diamètre interne (mm) 172 60
Diamètre externe (mrn) 190 78
Ravon de l'outil (mm) 1.4 1.4
1 Force amliauée sur le cercle de tête ml-1 1 9.8 1 9.8 1
1 Inchaison de la force nar ramort à l'axe "x" (O1 1 -160 1 -160 1
Entraxe (mm) 1 56 .
46 L'influence de la translation sur I'axe x sur le profil des dents
% transformation conforme il quatre termes

ciear
95 Coefficients de la transformation conforme
A=[. 10128318670.24924230500 .O1 181172998 .O01711074321;
B=[.09959437338 .4352 1517630 .O2001055915.004242493981;
al=A( 1); a2=A(2); a3=A(3); a4=A(4);
b l=B(l); b2=B(2); b3=B(3); b4=B(4);
c=.94676929 150;
% Différentes translations
dl=1.5; d2=3; d3=.5;
% Angle de rotation de la dent, en radians
f=.3 14159;
% Coordonnée curviligne liée à la surface de la dent
t=-4: .O 1:4;
% Coordonnées de la dent réeiie "XR","YR
XR=.9467*t+. 1012*tJ(t.A2+.0995A2)+.2492*tJ(t.A2+.4352A2)+.0 118*tJ(~~2+.02Y)+
-0017*tJ(t."2+.0042Y);
YR=.1012*.0995J(t.AZ+.0995A2)+.2492*.435ZJ O 118*.02J(tA2+.0P2)+
-0017*.0042J(t.A2+.0042A2);
95 Coordom6es de la dent centrale "XO", "YO"
X k *t+al *t/(t?2+b 1A2)+a2*t/(t.A2+b2A2)+a3*tl(t.Y+b3A2)+a4*t/(t."2+b4"2);
Y k a l *b 1J(tA2+b1A2)+a2*b2J(t.A2+b2A2)+a3*b3J(t.Y+b3A2)+a4*b4J(t.A2+b4A2);
% L'influence de la première translation sur la coordonnée curviligne
t l=t-d1;
t2=tid 1;
% Explicitation de la transformation conforme
D l 1=t1.A2+b1A2; D12=t 1!2+bZA2; D13=t1.A2+b3A2; D 14=tl.A2+b4Y;
D2 1=tZA2+b lA2;D22=t2.A2+b2A2;D23=t2.A2+b3A2; D24=tZA2+MA2;
56 Coordoinées de la transformation conforme pour la premième translation
Xt=Xk*tl+al*tlJDl l + a 2 " t l . / D l 2 + a 3 * t l J D 1 3 2 1 +
a2*t2JD22+a3*t2.n>23+a4*t2JD24;
Yt=Yû+al*blJDl l-+a2*b2JD12+a3*b3JD13+a4*b4JDf4+al *blJD21+a2*bSJD22+
a3 *b3JD23+a4*b4JD24;
96 L'iauence de la deuxième translation sur la coordonnée curviligne
t l%t-d2;
t22=t+d2;
46 Explicitation de la transformation conforme
D 121=t12.A2+b1A2;
D 122=t12.A2+b2A2;
D 123=ti2.A2+b3A2;
D 124=t12.A2+b4A2;
D221=t22?2+b lA2;
D222=t22.A2+b2A2;
D223=t22.A2+b3A2;
D224=t22?2+MA2;
% Coordonnées de la transformation conforme pour la deuxième transIation
X~=Xk*t12+al*t12.D121+a2*t12~122+a3*ti2JD123+a4*t12JD124+c*t22+
al *t22D22 1+a2* t22JD222+a3*t22JD223+a4*t22JD224;
Yt2=YO+al*b lfD12l+a2*b2JD 122+a3*b3fD123+a4*b4JD124+al*b l./D221+
aZ*b2JD222+a3*b3 JD223+a4*b41D224;
46 L'influence de la troisième translation sur Ia coordomde curviligne
t 13=t-d3;
t23=t+d3;
Dl3 1=t13-A2+b1A2;
D 132=t13?2+bZA2;
D l33=tl3?2+b3^2;
D 134=t13.A2+b4A2;
D23 l=t23?2+b P2;
D232=t23.A2+b2A2;
D233=t23.A2+b3A2;
D234=t23 .A2+b4AZ;
% Coordonnées de la transformation conforme pour la troisième translation
XO=XO+c*t13+al*t131D131+a2*t13JD132+a3*t13JD133+a4*t13JD134+c*t23+
al *t23./D23 l+aZ*tU./D232+a3*tU./D233+a4*t23JD234;
YO=Yû+al*bl./D13 l+a2*bZID132+a3*b3JD133+a4*b4JD134+al*blJD23 1+
aî*bZJDU2+a3 *b3JD233+a4*b4./D234;
% Axes des coordom6es
Ya=[-1.5 31;Xa=[O O]; Xal=[-15 151; Yal=[O O];
 Xall ,'-.')
pi0t~t,Yt,Xt3,YU.'-', Xt2,Yt2,Y,Xa,Ya,'-.',a
tit.le(LINFLUENCE, SUR LES 3 DENTS, DE LA TRANSLATION SUR xt)
end
% Étapes de la transformation conforme pour trois dents

clear
% Coefficients de la transformation conforme pour une dent
A=[. 1012831867O.2492423O5ûû .O1 181172998 .O017 11074321;
B=[.09959437338 1517630 .O200 IO559 15 .Qû424249398];
al=A( 1); a2=A(2); a3=A(3); a4=A(4);
bl=B(l); b2=B(2); b3=B(3); b4=B(4);
c=.94676929 150;
% Valeur de la translation sur I'axe "x",-ci-
d= 1.5;
8 L'angle de rotation de la denken radians,"fl
f=.3 14159;
% Coordonnée curviligne liée à la surface du domaine
t=-2.5:.005:2.5;
% Coordonnées de la dent centrale,'XO","YO"
Xk* t+al *tJ(tmA2+b 1A2)+a2*t/(tA2+b2A2)+a3*tJ(t.A2+b3A2)+a4*t./(t.A2+MA2);
Y b 1*b 1./(P2+b 1A2)+a2*b2./(tA2+b2A2)+a3 *b3J(toA2+b3 A2)+a4*b4./(t."2+MA2);
% Rotation d'un profü
XOrot=XO*cos(f)-YO*sin(O;
YOrot=YO*cos(f)+XO*sin(f);
% Translation à droite de Ia dent centrale
t 1=Gd;
% Explicitation de la transformation corifonne pour la dent translatée à droite
D 1l=t 1."2+b lA2;D 12=t1.A2+b2A2;D 13=tl."Z+bY2; D M=t 1."2+MA2;
% Coordonnées de la dent translatée 5 droite
Xtd=c*tl+al*tlJDl l+a2*tlJDl2+a3*tlJD13+a4*tlJD14;
Ytd=al*blJDl lta2*bZJDl2+a3*b3./Df 3+a4*b4JDl4;
% Coordonnée de la dent translatée à droite sur I'axe x, sur I'axey et avec une rotation
Xd=(c*tl+al*tl./Dl l+a2*tl JD12+a3*tlJD13+a4*tlJD14)*cos(f)+(a~*b1./D1 l+a2*b2./
D 12+d *b3JD13+a4*b4./D14)*sin(f)-d*tan(V2)*sin(f);
Yd=-(c*tl+al*tlIDl I+a2*tlJD12+a3*tlJD13+a4*tlJD14)*sin(f)+
+(al*bl JD1 l+a2*b2JD12+a3*b3JD13+a4*b4JD14)*cos(f)-d*~(~2)*cos(f);
% Translation à gauche de la dent centde
t2=t+d;
% Explicitation de la transformation conforme pour la dent translatée à gauche
D21=t2?2+b lY; D22=tZA2+b2r; D23=t2.A2+b3A2;D24dLA2+b4Y;
% Coordom6es de la dent translatée à gauche
Xtg==*d+al *YLJD21+a2*rZJD22+a3*t2JD23+a4*t2./DZ+
Y t g a l *b 1./D21+a2*b2JD22+a3*b3JD23+a4*64JD;
945 Coordonnées de la dent transIatée à gauche sur i'axe x, sur I'axe y et avec une rotation
Xg=(cc*t2+al*t2JD21+a2*t2JD22+a3*t2./D23+a4*t2JD24)*cos(f)-
(al *b 1JD2 1+a2* b2JD22+a3 *b3./DU+a4*b4JD24)*sin(f)+d*tan(f/2)*sin(f);
Yg=(c*t2+al *t2JD2l+a2*t2JD22+a3*t2.fD23+a4*t2./D24)*sin(f)+(al *b 1JD2 l+a2*b2./
D22+a3*b3JD23+a4*bUDî4) *cos@)-d*tan(f72)*cos(f);
% Deux dents translatées: celle centmie et celle de droite
Xt 1=XO+Xtd;
Ytl=YO+Ytd;
% Trois dents obtenues par translation
Xt=XO+Xtd+Xtg ;
Yt=YO+Ytd+Ytg;
46 Transformation conforme complète des trois dents
XTC=Xo+(Xd+Xg) ;
YTC=YO+(Yd+Yg) ;
% Transformation conforme pour deux dents: dent centrale et celle de gauche
XTcg=xo+xg;
YTCg=YO+Yg;
46 Transformation conforme pour deux dents: dent centrale et celle de droite
XTCd=Xû+Xd;
YTCd=YO+Yd;
% Axes des coordom6es
Ya=[-3 31; Xa=[O O]; Xal=[-3 31; Yal=[O O]; Xar=-Xal*atan(f);
Y M a l *tan(f);Xard=Xal*atan(f>;Yard=-Xa( 0 ;Xarl=Xal *atan(f);
plot(Xa,Ya,'-.rt,Xa1 ,Yal ,'-.rl,Xal ,Yar,'-.b',Xar,Ya,'-.b',XO,YO,t-b',XUrot,YOrot,'r')
title('POSITI0ND UN PROFIL MPC APRES ROTATION')
pause
plot(XO,YO,X;iYa,'-.',Xal ,YaI ,Y)
title('DENT CENTRALE, WC')
pause
plot(Xtd,Ytd,Xa,Ya,'-.'Jral ,Yal,'-.')
 title('DENT TRANSLATÉE SUR L AXE x')
pause
piot@t&Ytd-d*tan(U2)~a,Ya,'-.'~al ,Yal,'-. ')
title(DENT TRANSLATÉES U R L AXE y')
pause
pI~t(Xd,Y~'r'~~Ya,'-~'$Cal ,Ya 1,'-.'Sard,Ya,'-.',XaI,Yard,'-.')
tide('DENT MPC APRES TRANSLATION SUR x, SUR y, ET ROTATION)
pause
Ya=[-1 31; Xa=[O O]; Xal=[-7 41; Yal=[O O];
pIot(Xa,Ya,'-.',Xal ,Yal ,'-.',Xd,Yd,'-
'Jrd(950),Yd(950)~o'$CTCd(95O)1YTCd(950)1'o',~C~~Cd~O,Y~,~~(~~),
XTCd(250),YTCd(zSO),'* ')
title('DEUXDENTS MPC translation à droite')
pause
Ya=[-3 331; Xa=[O O]; ?Cal=[-3 31; Yal=[O O]; Xar=-Xal*atan(f);
Y-Xa 1*tan(f) ;Xard=Xal *atan(f);Yard=-Xa 1( 0 ; Xar l=Xal *atm(&
plot(Xg,Yg~a,Y~'-.',Xa1 ,Yal,'.'Xar,Yq'-.',Xal ,Yar,'-.')
title('DENT MPC APRES TRANSLATION SUR x, SUR y, ET ROTATION')
pause
Ya=[-1 33; Xa=[O O]; Xal=[4 71;Yal=[O O];
plot(Xa,Ya,'-.',Xal,Yal ,'-.',Xg,YgT8--
'X~(~~),Y~(S~);~'XTC~(SO),YTC~(SO),'O'~C~,YTC~~O,YO*X~(
g(95O),YTCg(950),'*')
titIe('DEUX DENTS MPC translation à gauche')
pause
Ya=[-1 31; Xa=[O O]; Xal=[-7 71;YaI=[O O];
plot(XTC,YTC,Xa,Ya,'-.',Xal,Yal ,'-.',Xar-4.4,Ya,'-.',Xard+4.4,Ya,'-.')
title(' TROIS DENTS MPCJRANSLATION -SUR x ET y- ET ROTATION')
pause
pIot(X&Y~O,YO,X~Ya,'-*'ZaI ,Yal ,'-.'~a4.6,Y~'-.'$Ca+4.6,Ya,'-.')
title('TR0IS DENTS TRANSLATÉES MPC')
pause
~Io~(XI'C~,YTC~,'-.'JCTC~,YTC~,'--
',XTC,YTC,XTC~(~~O),YTC~(~~O),'O'$~~C(~~O),YTC(~~O),
'oZTCd(5ûû),YTCd(5Oû),
'*'Xrcg(500).YTCg(500),'*',xrc(500),YTC(),'*')
title('P0SïïION DE CHAQUE D E W )
pause
al,'-.')
pi0t~&Yk'~JCrC,YTC-.56978$ra,Ya,'-.'~al
Akinput (('amin4,amax4. inca6 ')
B l=input ('bmin 1, bmaxl, incb l= ')
BLinput ('bmin2,bmax2. incb2= ')
B3=input (%min3. bmax3, incb3= 3
B4input ('bmui4, bmax4?incb4= ')
C-input ('cmin,cmax, incc= ')
Einput ('emin, emax, ince= ')
nal=(Al(2)-Ai(l))/Ai(3);
na2=(A2(2)42(1))/A2(3);
na3=(A3(2)-A3(l))/A3(3);
na4=(A4(2)-A4(i))/A4(3);
nb l=(B 10)-B l(I))/B l(3);
nb2=@2(2)-B2(i))/B2(3);
nb3=@3(2)-B3(1))/s3(3);
n%-(B4(2)-B4(1))/B4(3);
nc=(C(2)-C(l))/C(3);
ne=(E(2)-E(l))/E(3);
il=l;i2=l;i3=l;ikl;j l=l;j2=lj3=1;j4=l;w=l;
om=1;
for v= l :ne
routine
end
[RGOPI?(l) RGOPTI(2) RGOETi(3) RGOPTI(4) RGOPTI(5) RGOP'ïI(6) ...
RGOPTI(7) RGOPTI(8) RGOPTI(9) RGOP'ïI( 10) OPTI];
v=RGOPTI( 1O);
om=1;
for w=l:nc
rouîine
end
~GOPTI(1)RGOPTI(2) RGOPTI(3) RGOPTI(4) RGOlTI(5) RGOPTI(6) ...
RGOPTI(7) RGQPTI(8) RGOPTI(9) RGUPTI( 10) O m ;
w=RGOPTI(9);
OPTI=l;
for j k l :nb4
routine
end
 IRGOm(1) RGOPTI(2) RmPTI(3) RGOPTI(4) RGOPI?(S) RGOPTI(6) ...
RGOPTI(7) RGOPTI(8) RGOPTI(9) RGOPTI(I0) O m ;
jkRGOFTI(8);
OPfl=I;
for j3=l :nb3
routine
end
[RGOPTI(l) R o P T I ( 2 ) RGOPTI(3) RCOPTI(4) RGOPII(5) RGOPTI(6) ...
RGOPTI(7) RGûPTI(8) RGOPTI(9) RGOPTI(l0) OPTI];
j3=RGOPTI(7);
om=1;
for j2=l :nb2
routine
end
[RGOPTI( 1) RGûPTI(2) RGûPTI(3) RGOPTI(4) RGOPTI(5) RGOPTI(6) ...
RGOFTI(7) RGUPTI(8) RGOPTI(9) RGOPTI(I0) OPTQ ;
j2=RGOFTI(6);
om=1;
for j l=l:nb l
routine
end
[RGOPTI(l) RGûP'ïI(2) RGOFTI(3) RGOPTI(4) RGOPTI(5) RGOPTI(6) ...
RGOPTI(7) RGOPTT(8) RGOPTI(9) RGOETI(1O) OPTI];
j I=RGOPTI(S);
Orn=l;
for i4=1 :na4
routine
end
mGOFTI(i) RGOPTI(2) RGOPTI(3) RGOPTI(4) RGOPII(5) RGOFTI(6) ...
RGOPTI(7) RGOPTI(8) RGOPTI(9) RGûPTI(10) OPTI];
ikRGOPTI(4) ;
OPTI=l;
for i3= 1 :na3
routine
end
[RGOFTi(l) RGOPTI(2) RGOPTI(3) RGOFTI(4) RGOPTI(5) RGOPTI(6) ...
RGoPT?:(7)RGOPTI(8) RGOPTI(9) RGOPTI(10) OPTI];
i3=RGûPTI(3);
Orn=l;
for ï2= 1 :na2
routine
end
[RGOm(l)RGOm(2) RGOPTI(3)RGOP'ïI(4) RGOPTI(5) RGOm(6) ...
RGOPTI(7) RGOPII(8) RGOFTI(9) RGOPTI(10) OPTI];
i2=RGOPI1(2);
OPTI=l;
for il=l:nal
routine
end
LRGOPTI(1)RGûP'i'I(2)RGOFïI(3) RGOFTi(4) RGOPTI(5) RGûPTI(6) ...
RGOPTI(7) RGOPTI(8) RGOPTI(9) RGOPTi(l0) OPTI];
% CALCUL DE LA DENT M T I L E
D 1 10=(ttA2)+@1(1)"2);
D 120=(~~2)+(B2(1)~2);
D l3û=(~~2)+@3(1)~2);
D 140=(t.A2)+(B4(1)A2);
D2Iû=(t 1.A2)+(B1(1)A2);
D220=(t 1 .A2)+(B2(1)A2);
Dî3O=(t l.A2)+(B3(l)A2);
D24û=(t 1.A2)+(B4(l)A2);
D3 10=(t2.A2)+(B1( 1)A2);
D32û=(t2.Y)+(B2(l)A2);
D33û=(t2.A2)+(B3( 1)A2);
D34û=(t2.A2)+(B4( 1)^2);
xlû=C(l)*t+Al(l)*t/D1lO+A2(1)*tJD120+A3(1)*tJD 130+A4(1)*tJD140;
x20=(C(1)*tl+Al(1)*tlJD2lO+A2(l)*tlJD220+A3(1)*tl/DUO+A4(1)*tlID240)*...
cos(teta)-(Al(l)*B l(l)JD210tA2(1)*B2(1)JD220+A3(1)*B3(l)JD230+A4(1)*...
B4(1)JD240)*sin(teta);
x3û=(C(l)*t2+Al(l)*t2./D310+A2(1)*t2JD320tA3(1)*t2JD33û+A4(1)*t2JD340)*...
cos(teta)+(Al(l)*B 1(1)JD3 lChA2(1)*B2(1)D320+A3(1)*B3(1)JD33û+A4(1)*B4(1)j...
D34O)*sin(teta);
XTco=x10+x20+~30;
ylO=AI(I)*B 1(1)JD110+A2(1)*B2(1)ID120+A3(l)*B3(1)JDI30+A4(1)*B4(l)JD1~;
y~0=(~(l)*tl+Al(l)*tlJD2lQeA2(l)*tlJD220+A3(l)*t~lD230tA4(l)*tlJD240)*...
sin(teta)+(Al(l)*B 1(1)~210+A2(1)*B2(1)ID220cA3(1)*B3(1)/D230+A4(1)*B4(1)./...
D24o)*cos(teta);
~3Q-(C(l)*t2+Al(l)*t2JD3lO+A2(1)*i2JD32U+A3(1)*t2JD330çA4(l)*t21D340)*...
sin(teta)+(Al(1)*B l(l)JD3 10+A2(1)*B2(1)A320+A3(1)*B3(1)JD33û+A4(1)*B4(l)J...
D340)*cos(teta);
YTCû=ylû+y20+y30-2*d*sin(teta)-E(1);
% COEFFICIENTS OPTIhriALS
a1=Ai(i)+(RGOfTI(l)-l)*A1(3);
a2=A2(1)+(RGOETI(Z)- 1)*A2(3);
a3=A3(1)+(RGOPI1(3)- L)*A3(3);
a4=A4(1)+(RGOFïl(4)- 1)*A4(3);
b l=B 1(l)+(RGOPTI(S)-l)*B 1(3);
b2=B2(l)+(RGOPTI(6)- l)*B2(3);
b3=B3(1)+(FtGOFTI(7)- l)*B3(3);
bkB4(1)+(RGOPII(8)- l)*B4(3);
c=C(I)+(RGOPTI(9)-l)*C(3);
e=E(I)+(RGOPTI(lO)- l)*E(3);
% CALCUL DE LA DENT OPTIlMALE
D 1l=(teA2)+(b1A2);
D 1Z=(t."2)+(bZY);
D 13=(t."Z)+(bY2);
D 14=(~~2)+@4"2);
t l=t+d;
t2= t-d;
D2 1=(tl ?2)+(b lA2);
D22=(t 1.Y)+(b2A2);
D23=(t l .A2)+(b3A2);
D24=(t 1.Y)+(MA2);
D3 1=(t2.Y)+@ lA2);
D32=(t2.Y)+(b2A2);
D33=(t2.A2)i(ô3A2);
D34=(tP2)+@tA2);
xl=c*t+al*tJDl l+a2*t/D12+a3*tJD13+a4*t./Dl4;
x2=(c*t l+al *t1JD2l+a2*tl ./D22+a3*t 1JD23+a4*t1JD24)*cos(teta)-...
 (al *b1JD2I+a2*b21D22+a3*b3JD23ia4*b4JD24)*sin(teta);
x3=(c*t2+al *t2JD3l+a2*t2JD32+a3*t2JD33+a4*t2JD34)*cos(~(a *b 1JD3 i+a2*...
62JD32+a3*b3JD33+a4*b4JD34)*sin(teta);
XTCo=xl+X2+x3;
yl=al*blJDll+a2*b2JDl2+a3fb31D13+a4*b4JD14;
y2=(c*tl+al*tl/D2l+a2*tlJD22+a3*tlJD23+a4*tlID24)*sin(teta)~al *bl ID2l+a2*...
b21D22+a3*b3JD23+a4*bQJD24)*cos(teta);
y3=-(c*t2+al *t21D3 l+a2*t2JD32+a3*t2JD33+a4*WD34)*sin(teta)+...
+(al *blJD3 I+aZ*b2JD32+a3*b3JD33+a4*b4JD34)*cos(teta);
YTCo-y l+y2+y3-2*d*sin(teta)e;
XR=1.O375171*t+. 10850075*t/(t.A2+.098935225A2)+. 19659848*tJ(t.A2+.40113083A2)+...
.010490963*tl(t.A2+.019692ûû3A2)+.ûû 19275009*tl(t.A2+.004u1 l24W2);
YR=.lO8SOO75*.O98935225J(t.A2+.09893S225A2)+.19659848*.40113083J(t.A2+...
.4O 1 13083A2)+.010490963*.019692003J(t.Y+.O 19692003A2)+...
.O019275009*.004231 1243J(~"2+.0042311243"2)-.16550389;
pIot(xR1 ,YR1,xR2?YR2XTcO,YTcO,'-.',xTco,YTCo)
titie(r al= ',num2str(al),'a2= ',num2str(a2),'
a3= ',niim2str(a3),
a4=
' ',num2str(a4),'
bl=

pause
format long
[al a2 a3 a4Jt, ml b2 b3 b4]'
[C e]',[d teta]

pause
 % Sous-routine pour le calcul des coeffkients de la transformation conforme

96 CALCUL DE LA DENT
global t d f i l i2i3 i4jl j Z j 3 j 4 w v A l A2A3 ~ 4 BBl 2 ~ B3 ~ C E
al=Al(l)+(il-l)*Ai(3);
a2=A2( 1)+(i2-i)*A2(3) ;
a3=A3(1)+(3- l)*A3(3);
a4=A4( l)+(i4- l)*A4(3);
bl=B 1 ( l ) N 1-l)*B 1(3);
bZ=BZ(l)+(j2-1)*B2(3);
b3=B3(l)+(j3-l)*B3(3);
b4=B4( l)+(j4- l)*B4(3);
c=C(l)+(w- l)*C(3);
e=E( l)+(v- I)*E(3);
D l l=(t."2)+@ lA2);
D l2=(t."22)+@2y);
D 13=(t.Y)+(b3Y);
D 14=(t.A2)+@4A2);
t 1=t+d;
t2=t-d;
D2 l=(t P2)+(blA2);
D22=(t 1."2)+(bZA2);
D23=(t 1.A2)+(b3A2);
DX=(t 1.^2)+@tA2);
D3 1=(t2.A2)+@1A2);
D32=(t2.A2)+(b2A2);
D33=(t2.A2)+(b3A2);
D34=(t2!2)+(MA2);
xl=c*t+al*t./Dl l+a2*t/D12+a3*tlD13+a4*t/D14;
x2=(cf t lia1 WJD2 1+a2*t1JD22+a3*t 1JD23+a4*tIJD24)*cos(teta);
x3=(c*t2+al *t2JD3l+a2*t2JD32+a3*t2JD33+a4Irt2JD34)*sin(teta);
x4=(-(al *b 1lDZl+a2*b2./D22+a3*b3JD23+a4*b4./D24)+..-
al *bl JD3 1+a2*b2./D32+a3*b3JD33+a4*b4JD34)*sin(teta);
XTC=x l+xî+fi+x4;
yl=al*bl./D1l+a2*b2./D12+a3*b3JD13+a4*b41D14;
y2=(c*tl+al*tlJD21+a2*tlJD22+a3*tlJDu+a4*tlJD24-...
(c*tZ+al*t2JD3l +a2*QJ'D32+a3*t21D33+a4*4/D34))*sin(teta);
y3=(al *b lJM l + a 2 * b 2 1 D 2 2 + a 3 * b 3 J D 2 3 + a 4 * b 4 J D . .
D32+a3*b3JD33+a4*b41D34)*cos(teta);
YTC= l+y2+y3-2*d*sin(teta)s;
dif=diff(YTC);
for k=I:60
if (dif(k)>O).break,end
end
LIM=XTC(k);
% calcul de la fonction d'optimisation
1=0;
for k=1:60
if YTC(k)<=max(YR1)
if YTC(k)>=min(YR1)
if XTC(k)<=LiM
1=1+ 1;
XC(l)=XTC(k);
YC(I)=YrC(k);
end
end
end
end
pas=YRl(Z)-YRl(1);
for k=l:l
n 1(k)=((YC(k)-min(YR1))/pas) ;
end
n=max(n 1);
xi=XR 1(n);
xs=XR 1(n+ 1);
yi=min(YRl)+pas*(n-1);
for k= 1:1
xint(k)=xi+(xs-xi)*(YC(k)-yi)/pas;
end
FOP=~~m(((xint-xC) ."2)A);
if (FOP<OPTI)
OPTI=FOP;
RGOKi(l)=i 1;
RGOPTI(Z)=i2;
RGOPTI(3)=i3;
RGOETI(4)=i4;
RGOPTI(S)=j 1;
RGUI?TI(6)=j2;
RGOPTI(7)=j3;
RGOITI(8)=j4;
RGOPTI(9)=w;
RGOPTI( 1O)=v;
else
end
% Rotation d'un profil donné

clear
% LE PROFIL DE LA DENT cALCULÉ AVEC GENERl[5]:
% XO,YO: coordonnées du chaque point
% R: le rayon sur lequel se trouve le point
% N: nombre de points pour tracer le profil
N=30;
XO=[1.3 13 1.117 1.029 -972 -933 .907 390 380.877 -878,876 A69 -859.846 330.8 12
-791 -768 .743 .717.688 .657.625.590.554.516 -476 -435.392 -3471;
Yb[-.O99 -.O18 .O63 -143 ,224 -305-386 -467 -547 .628 .709 -790 -871 .951 1.032 1.1 13
1.194 1.275 1.356 1.436 1.517 1.598 1.679 1.760 1.840 1.921 2.002 2.083 2.164 2.2451;
R=[8.75 8.803 8.872 8.946 9.023 9.100 9.179 9.259 9.339 9.419 9.500 9.579 9.659 9.738
9.8 17 9.896 9.975 10,054 10.133 10.212 10.290 10,369 10,448 10.526 10.605 10.684
10.763 10.842 10.921 11.ûûûj;
W CALCUL DES COORDONNÉESPOUR LE PROFIL QUI A SUBI UNE ROTATION
46 L'angle de rotation du prohl (en radians)
teta=18*piJ180;
% Angle d'inclinaison de chaque point du profil par rapport à l'axe de symétrie de la dent
% centrale
for i= 1:N
alfa(i)=asin(XO(i)JR(i));
end
W Angle d'inclinaison de chaque point du profilaprès rotation, par rapport à l'axe de
8 symktrie de la dent centde
for i=l :N
alfa1(i)=teta-alfa(i);
% Coordonnées du profil après rotation
Xrot(i)=R(i) *sin(alfa1(i));
Yrot(i)=R(i)*cos(alfa1(i))-R( 1);
end
pIot(XO,YO,Xrot,Yrot)
title('R0TATION D UNE DENT DoNNÉE')
end
46 Calcul des contraintes superficieï.Ies pour trois dents

clear
% Paramètres connus: force appliquée "W"
%6 point de l'application de la force "to"
% angle d'inclinaison de la force "beta"en degrées
W=l;
to=.0035;
km=-I5O.S;
% Coefficients de la transformation conforme pour une dent: CARDOU et TORDION [8]
AC=[.10074 -26740 -01211 .ûû175];
BC=[.09864.4663 1 .O205 1.00413];
cc=.% 118;
ac l=AC(l); ac2=AC(2); ac3=AC(3); ac4=AC(4);
bcl=BC(l); bc2=BC(2); bc3=BC(3); bc4=BC(4);
% Coefficients de la transformation conforme pour trois dents dents
A=[.09999.27524 .O t 181 .ml911;
B=[.09689 S052 I .O190 1 -003941;
c=.48775;
al=A( 1); a2=A(2); a3=A(3); a4=A(4);
bl=B(l); b2=%(2);b3=B(3); b4=B(4);
% Valeur de la translation sur I'axe "xW,-d-
d= 1.75;
% L'angle de rotation de la dent, en radians. "f'
f=18*piA80;
% Coordonnée curviligne liée à la surface du domaine
t=-2.5:.0 1:2.5;
% Coordonnées de la dent CARDOU,"XC","YC"
XC=cc*t+ac I *t./(~~2+bc 1A2)+ac2*t./(tn2+bc2A2)+ac3*tJ(t.A2+bc3 A2)+ac4*t./.-.
(tA2+bc4A2);
YC=ac 1*bc 1J(t.A2+bc1A2)+ac2*bc2J(t.A2+bc2A2)+ac3 *bc3J(t,A2+bc3A2)+ac4*b~4./...
(tAZ+bc4AZ)+.203#;
% CALCUL DE LA CONTRAINTE SUPERFICIELLE DANS UNE DENT ISOLÉE
% Calcul de la dérivée de la msfommtion conforme pour la dent centrale
for k=1:4
for j=lA
ECAR(k)=ac I ./(BC(k)+bc 1).A2+acU(BC(k)+bc2).A2+~3
Y..
(BC(k)+bc3)Î2+a~4J(BC(k)+bc4).~2;
OMPCAR(k)=cc+sum(ECAR(k));
% Calcul de ALFACAR(i,k)
OMCAR(k)=AC(k)./OMPCAR(k);
ALFACAR(j ,k)=OMCAR@)I((BC(j)+BC(k)).A2);
% Calcul des coefficients réels Ru), Sa)
RC(j)=((to*cos(beta*pi/l 80)+BC(i).*sin(beta*pi/l80))./...
(toA2+BC(i)!2))* W/(2*pi);
SC(j)=((to*sin(beta*pi/l80)-BC(j).*cos(bel go))./ ...
(toA2+BC(j).^2))*W/(2*pi);
MATALFACAR(j,k)=[ALFACAR(j ,k)1;
end
end
for k= 1:4
for j= 1:4
if k=j
U(j,k)= 1;
else
U(j,k)=O;
end
end
end
% Calcul des variables réelles XCAR(j)
for k=1:4
for j=l:4
MATXCAR(j,k)=MATALFACAR(i ,k)+U(j,k);
end
end
XCAR=(MATXCAR)\[RCP;
% Calcul des variables réelles YCARU)
for k=1:4
for j= l:4
MAnCAR(j,k)=MATALFACAR(i&)-UUa);
 end
end
YCAR=(MATYCAR)\[(-Sq';
46 Calcul de la contrainte superfkielle
ElCAR=(OMC~(~)*gCCAR(l)*(tA2-bcl.A2)+2*YCAR(1)*bcI *t)./(tA2+bc1.A2).A2+...
0MCAR(2)*(XCAR(2)*(t.A2-bc2.A2)+2*YCAR(2) *bc2*t)./(t/"'+b~2.~2)."2+.. .
OMCAR(~)*(XCAR(~)*(P~-~C~ .A2)+2*YCAR(3)*bc3*t)./(t.A2+bc3.A2).A2+...
OMCAR(4)*~CAR(4)*(t.A2-bc44A2)+2*YCAR(4)*bc4*t)./(~A2+bc4.A2) ?2);
E 1lCAR=-(cos(beta*pi/l8O)J(t-to))*W/(2*pi)+E ICAR;
E2CAR=(OMCAR(l) *(2*XCAR( I)*bc 1*t-YCAR( 1)*(t. A 2 - bl~.A2))./(t.A2ib~ 1."2). "2+.. .
O M C A R ( ~ ) * ( ~ * X C A R ( ~ ) * ~ C ~ * ~ - Y C A R ( 2 ) * ( t .. A2-b2+.
OMCAR(3)*(2*XCAR(3)*bc3*t-YCAR(3)*(t Y-bc3 ."2))J(t A2+b~3.A2).A2+...
OMCAR(~)*(~*XCAR(~)*~C~*~-YCAR(~)*(~~~-~C~.~~))./( .A2);
E 12CAR=-(sin@eta*pill80)./(t-to)) *W/(2*pi)+E2CAR;
D lCAR=acl *(tA2-bc1A2)./(t.A2+bclA2).A2+ac2*(tY-bc2A2)./(t.A2+bc2A2).A2+...
ac3*(t?2-bc3 A2).l(tA 2 + b c 3 A 2 ) . A 2 + a c 4 * ( t z 2 ;
Dl lCAR=cc-DICAR;
D 12CAR=-2*t.*(ac l *bc 1./(t!Z+bc P2).Y+a~2*bc2./(t~2+bc2~2) ?2+...
ac3*bc3J ( ~ ~ 2 i b c 3.A2+ac4*bc4./(t.A2+bc4A2)
~2) !2);
SIGMA-*(E 1lCAR.*D 1lCAR+E12CAR.*D 12CAR)J(D 11CARA2+D12CARA2);
SIGMAMIGMAO ./4;
% Calcul des coordonnées cartésiennes pour la contrainte superficielle
% SIGMAcO est la contrainte en compression pour une dent isolée
% SIGMAtO est la contrainte en tension pour une dent isolée
DIM=length(XC) ;
ü=i;ij=i;
for i=2:@IM)
if XC(i)c=û
H lcCAR=sqrt((XC(i- 1)-XC(i))A2+(YC(i)-YC(i-
1)IA2);
XSIGMAcO(ii)=XC(i)-SIGMAO(i)*(YC(i)-YC(i- 1))./H1cCAR;
YSIGMAcO(ii)=YC(i)-SEMAO(i)*OLC(i- 1)-XC(i))./H 1cCAR;
ü=ü+l;
else
H 1tCAR=sqrt((XC(i)-XC(i- 1))A2+(YC(i-1)-YC(i))A2);
XSIGMAtO(ü)=XC(i)-SIGMAO(i) *(YC(i- 1)-YC(i))./H1tCAR;
YSIGMAto(jj)=YC(i)-SIGMAO(i)*(XC(i)-XC(i1))JHl tCAR;
 j+j+l;
end
end
% CALCUL. DE LA TRANSFORMATION CONFORME POUR TROIS DENTS
% Coordomées de Ia dent centrde,"XO","YOt'
Xû=c*t+aI *tJ(tA2+b1A2)+a2*t/(~A2+b2A2)+a3*tJ(t.A2+b3A2)+a4*U(t-A2+MA2);
Yû=al *b 1J(tmA2+b 1A2)+a2*b2J(~A2+b2A2)+a3*b3J(tA2+b3A2)+a4*b4./(~A2+W2);
% Translation à droite de la dent centrale
t 1=t-d;
Dl l=tloA2+blA2;
D 12=t1!2+b2"2;
D 13=tl.AZ+b3AZ;
D M=t 1?2+bP2;
% Dent translatée a droite sur I'axe x, sur l'axe y et avec une rotation
Xd=(c*tl+al*tlJDl l+a2*tl JD12+a3*tlJD13+a4*tlJD14)*cos(f)+(al*b1./Dl1+...
a2*b2JD 12+a3*b3JD13+a4*b4JD14)*sin(f)-d*tan(fl2)*sin(f);
Yd=-(c*tl+al*tlJDl l+a2*tlJD12+a3*tlJD 13+a4*tlJD14)*sin(f)+...
(al*b lJDl l+a2*b2JD12+a3*b3./D13+a4*b4JD14)*cos(f)-d*tan(Y2)*cos(f);
95 Translation à gauche de la dent centrale
t2=t+d;
D2 1=QOA2+b lA2;
D22=t2.A2+b2A2;
D23=tP2+b3"2;
D24=t2.A2+b4A2;
Q Dent translatée à gauche sur i'axe x, sur i'axe y et avec une rotation
Xg=(c*Q+al *t2JD2I+a2*t2JD22+a3*t2JD23+a4*t2JD24)*cos(f)-(al *b 1JD2 l+...
a2*bZ./D22+a3*b31DW+a4*b4JD24)*sin(f)+d*tan(f72)*sin(f);
Yg=(c*Q+al*t2JD2l+a2*t2JD22+a3*iL.lD23+a4*t2JD24)*sin(f)+(al*b 1JD21+...
a2*b2.~22+a3*b3oD23+a4*b41D24)*cos(f)-d*~(fl2)*cos(Q;
% Transformation conforme complete des trois dents
XTC=XO+(Xd+Xg);
YTC=Yû+(Yd+Yg)--36634;
% CLCUL DES CONTRAINTES SUPERFICIELLES POUR 3 DENTS ADJACENTES
% Calcul de la dérivée de la transformation conforme
for k=1:4
for j=l:4
E(k)=a i J(B (k)+b1).^2+a2J(B(k)+b2).^2+a3J0.
a4./(B(k)+b4).A2+(al*(2*((J3(k)+b l).^2dA2)*c0s(f)+...
4*d*(B(k)+b1).*sin(f)))J((B(k)+b 1).A2+d.A2).A2+...
*--.
(â2*(2*((B(k)+b2).A2-daA2)*~~~(f)+4*d*(B0.
sin(f)))J((s (k)+b2)-n2+d,A2)hA2+(a3 *(2*((B(k)+b3).A2-.--
d.A2)*cos(f)+4*d*(B(k)+b3).*sin(f))./(@(k)+b3).A2td.A2) ?2+...
(a4*(2*((B(k)+b4).A2d.A2)*cos(f)+4*d*@@)+b4).**..
sin@))) J((B @)+!A).A2+d.A2).Y;
OMP(k)s*(1+2*cos(f))+sum(E(k));
9% Calcul de ALFA(j,k)
OM&)=A(k)./OMP(k);
ALFA(j,k)=OM(k)J((B(i)+B(k))?2);
96 Calcul de L(j,k),M(j,k)
H1 1(j,k)=-Am)*((B ( j ) + B ( k ) ) . A 2 d - A 2 ) . / ( ( B ( 2 ;
H 12ÿ,k)=2*A(k)*d*(B(j)+B(k))./(@ (j)+B(k)) ."2+;
T1(k)=al*((B (k)+b1).A2-d.A2)./((B(k)+b 1).A2+d.A2).A2+a2*((B (k)+b2)!2-. ..
d.^2)J((B(l~)+b2).~2+d-~2).~2+a3 *((BOc)+b3).A2d.A2)./-..
((B( k ) + b 3 ) . n 2 i d . A 2 ) . A 2 + a 4 * ( ( B ( k ) + b 4 ) . Y - d . 2 ;
=@)=al *(((B(k)+b 1).A224*d.A2)*cos(f)+4*d*(B(k)+b l)*sin(f))./...
((B(k)+b 1).A24*d.Y).A2+a2*(((B(k)+b2).A2-4*d.A2)*cos(f)+...
4*d*(B(k)+b2)*~in(f))./((B(k)+b2).~294*d.~2).~2+a3 *(((B(i~)+b3).~2-...
4*d-A2)*cos(f)+4*d*(B(k)+b3)*sin(f)) J((B(k)+b3)."2+4*d."2)?2+...
a4*(((B(k)+b4) . A 2 - 4 * d . A 2 ) * c o s ( f ) + 4 * d * ( B o . / . . .
((B(k)+t14).~2+4*d.~2)."2;
T3(k)=al *cos(f)./(B(k)+bI) .A2+a2*cos(f)./(B(k)+b2).A2+...
a3*~os(f)J(B(k)+b3).~2+a4*cos(f) J ( B ( k ) + b 4 ) ;
T1 1(k)=c*(1+2*cos(f))+sum(T 1(k))+sum(T2(k))+sum(T3);
T4(k)=Z*al*d*(B(k)+bl)J((B(k)+b 1) .A2+d.A2).A2+2*a2*d*(B(k)+b2)./...
((B(k)+b2).A2+daA2).A2+2*a3*d*(B(k)+b3) j(@(k)+b3) .A2d.A2).A2+...
2*a4*d*(B(k)+b4)./((s(k)+M).A2+d-A2).A2;
TS(k)=al*(4*d*(B(k)+b l)*cos(f)-((B(k)+b 1).A24*d.A2)*...
sin(f))J((B(k)+b l).A2+4*d.Y).^2+a2*(4*d*(BOc)+b2)*cos(f)- ...
((B (k)+b2).Y-4*d.Y)*sin(f))./((B ( k ) + b 2 ) . E w - .
a3*(4*d*(s(k)+b3)*co~(f)-((B(k)+b3).^24*d.~2)*sh(~)./(@&)+b3).~2+...
4*d.A2).A2+a4*(4*d*(B(k)+b4)*cos(f)-((B0..
((B(k)+b4).A2+4*d.A2).A2;
T6@)=al*sin(f)J(B(k)+b 1).A2+a2*sin(f)l(B(k)+b2).A2+...
a3*~in(f)J(B(k)+b3).~2+a4*sin(f)J(B(k)+b4)?2;
120r)=s=CrW))-s=V5(k))-sum(T6(k)) ;
(31 l(j&)=@l l(j,k).*TI l(k)+H12(j,k).*T12(k))J(TI 1(k).A2+T12(k).A2);
Gl2(jJr)=(H12(j,k).*T1l(k)-Hl l(j,k)-*T12(k))J(Tll(I~).~2+T12(k).Y);
L(jJc)=Gll (j,k)*cos(f)G12(j&)*sin(£);
M ( j , l ~ ) ~(j,k)*sin(f)+G12(jJc)*cos(£);
ll
OMl(k)=A(k)J(T 11(k)+i*T12(k));
OMZ(k)=A(k)l(T1l(k)-i*T12(k));
96 CaIcd de Z(j,k)).V(j,k)
Z(j,k)=H 11(j,k)JOMP@);
V(j,k)=H lZ(j,k)JOMP(k):
% chhlde N(jJc),O(jJr), P(jJc)?F(j.k)
G2 i (i,k)=(A(k)J(B(j)+B(l~).~2)). *T11(k)./(T11(k).^2+T12(k)?2);
G22(j,k)=(A(k)J(B(i)+B(k).Y)).*Tl20() J(T1l(k).A2+T12(k).A2);
EU 1(i,k)=-A@)*((Bu)+BO()).A2-4*ddA2) J((B(j)+B(k)).^2+4*dA2).A2;
H22(j,k)=PA(k)*d*(B (j)+B(k))J(@(i)+B(k)).A24*d.A2).^2;
G3 l(j,k)=(H2l(j,k).*Tll(k)+H22(j,k).*Tl2(k))J(Tl I(k).A2+T12(k).A2);
G32(j,k)=(H22(j9k).*Til(k)-H2i(j,k).*T12(k))./(Tl l(k).A2+T12(k).A2);
N(j,k)=G2l(j,k)*cos(f)+G22(j,k)*sin(f);
O(j,k)=G2 1(j,k)*sin(f)-G22(j,k)*cos(f);
% CaIcui de P(j,k)).Fÿ,k)
P(j,k)=G3 l(j,k)*~os(f)-G32(j,k)*~in(f);
F(jJc)=G31(j,k)*sin(f)+G32(j,k)*cos(f);
% Calcul des coefficients réels Rÿ), S(i),Rl(j), S l(i)JU(j), S2(j)
R(j)=((to*cos(beta*pi/I8O)+B(i). *sin(beta*pi/l80))./(toA2+B(j)!2))* ...
W/(2 *pi);
S(j)=((to*sin@eta*pi/i 80)-B(j).*cos(beta*pi/180))J(t0~2+B(j).~2)) *...
W/(2 *pi);
R 1(j)=(((to+d)*cos(beta*pi/l80)+Bÿ). *sin(beta*pi/lSO))J...
((to+d)^2+B(j). ^2))*W/(Z*pi);
S 1(j)=(((to+d)*sin(beta*pi/l80)-B (i).*cos@eta*pi/180))J...
((t~id)~Z+B(j).~Z))* W/(2*pi);
R2(j)=(((to-d) *cos(beta*pi/l80)+Bÿ).*sin(beta*pillSO))J...
((to-d)^2+B(j).^2))*W/(2*pi);
S2(j)=(((to-d)*sin(beta*pi/l80)-B(j).*cos(beta*pi/I 80))J...
 ((to-d)^2+B (j) .^2))*W/(2pi);
end
end
96 CALCUL DES COEFFICIENTS RÉELS
MATRICE=[l+ALFA(I,1) ALFA(1,2j ALFA(1,3) ALFA(1,4) 0 0 0 0 -L(1,1) -L( 1,2)
-L(1,3) -L(1,4) -M(l,l) -M(1,2) -M(1,3) -M(1,4) -L(l,l) -L(1,2) -L(1,3) -L(1,4) +M(l,I)
+M(1,2) +M( 1,3) +M(I ,4);
ALFA(2-1) 1+ALFA(2,2) ALFA(2-3) ALFA(2,4) O O O 0 -L(2,1) -L(2,2) -L(2,3)
-L(2,4) -M(2,1) -M(2,2) -M(2,3) -M(2,4) 442-1) -L(2,2) -L(2,3) -L(2,4) +M(2,1) +M(22)
+M(2,3) +M(2.4);
ALFA(3.1) ALFA(3,2) I+ALFA(3,3) ALFA(3.4) O 0 0 0 -L(3,1) -L(3,2) -L(3,3)
-L(3,4) -M(3,1) -M(3,2) -M(3,3) -M(3,4) -L(3,1) -L(3,2) -L(3,3) -L(3,4) +M(3,1) +M(3,2)
+M(3,3) +M(3,4);
ALFA(4,l) ALFA(4,2) ALFA(4,3) 1+ALFA(4,4) O O O O -L(4,1) -L(4,2) -L(4,3)
-L(4,4) -M(4,1) -M(4,2) -M(4,3) -M(4,4) -L(4,1) -L(4,2) -L(4,3) -L(4,4) +M(4,1) +M(4,2)
+M(4,3) +M(4,4);
O O O O 1-ALFA(1,l) ALFA(1,2) ALFA( 1,3) ALFA( 1,4) -M(1,1) -M(1,2) -M(1,3)
-M(1,4) +L(1,1) +L( 1,2) +L(1,3) +L(1,4) +M(1, 1) +M(1,2) +M(1,3) +M(1,4) +L(l, 1)
+L(1,2) +L( 13 ) +L( 1.4);
O O O O ALFA(2,l) 1-ALFA(2.2) ALFA(2-3) ALFA(2,4) -M(2,1) -M(2,2) -M(2,3)
-M(2,4) +L(2,1) +L(2,2) +L(2,3) +L(2,4) +M(2,1) +M(2,2) +M(2,3) +M(2,4) +L(2,1)
+L(23) +L(2,3) +L(2,4);
O O O O ALFA(3,I) ALFA(3.2) 1-ALFA(3,3) ALFA(3,4) -M(3,1) -M(3 2) -M(3,3)
-M(3,4) +L(3,1) +L(3,2) +L(3,3) +L(3.4) +M(3,1) +M(3,2) +M(3-3) +M(3,4) +L(3,1)
+L(3,2) +L(3,3) +L(3,4);
O O O O ALFA(4.1) ALFA(4-2) ALFA(4,3) 1-ALFA(4,4) -M(4,1) -M(4,2) -M(4,3)
-M(4,4) +L(4,1) +L(4,2) +L(4,3) +L(4,4) +M(4,1) +M(4,2) +M(4,3) +M(4,4) +L(4,1)
+L(4,2) +L(4,3) +L(4,4);
-Z(l, 1) -Z(1,2) -2(1,3) -2(1,4) +V(1,1) +V( 12) +V( 193) +V(1,4) 1+N(1,1) +N( 9 2 )

+N( 1,3) +N(1,4) +O( 1,I) a ( 1 . 2 ) +O( 1.3) +O( 1.4) -p( 1,l) -p( 1.2) -p( 1.3) -p( 1,4) +F( 1-1
+F(1,2) +F(1,3) +F( 1-4);
-Z(2,1) -2(2,2) -2(2,3) -2(2,4) +V@,1) +V(23) +V(2,3) +V(2,4) +N(Z 1) 1+N(2,2)
+N(2,3) +N(2,4) +O(Z, 1) +0(2,2) +0(2,3) +O(2,4) -p(2,1) -p(22) -p(2,3) -p(2,4) +F@, 1)
+F(2,2) +F(2,3) +F(2,4);
 -V(3,1)-V(33)-V(3,3)-V(3,4) +Z(3,l) +2(32)+Z(3,3) +2(3,4) -F(3,1)-F(3,2)
-F(3,3) -F(3,4)+P(3,1) +p(32) +P(3,3) +P(3,4) -0(3,1)-0(3,2)-0(3,3) -0(3,4)-N(3,1)
-N(3,2) 1-N(3,3)-N(3,4);
-V(4,1)-V(42)--V4,3)-V(4,4)+Z(4,1) +Z(43) +2(4,3)+2(4,4) -F(4,1)OF(4.2)
-F(4,3)-F(4,4)+P(4,1) +P(42) +P(4,3) +P(4,4) -0(4,1)-0(4,2)-0(4,3) -0(4,4)-N(4,1)
-N(4,2)-N(4,3) 1-N(4,4)];
RESULTAT=[R(I); R(2); R(3); R(4); S(1); S(2); S(3); S(4);
Ri(1); Rl(2); Rl(3); Rl(4); Sl(1); Sl(2); Sl(3); Sl(4);
R2(1); R2(2); R2(3); R2(4); S2(1); S2(2); S2(3); S2(4)];
% CALCUL DES VARIABLES RÉELLES
VR=(MATRIcE)\[RESULTAT];
46 IDENTIFICATION DES VARL4BLES R&LLES
X(l)=VR(l); X(2)=VR(2); X(3)=VR(3); X(4)=vR(4);
Y(1)=VR(S); Y(2)=VR(6); Y(3)=VR(7); Y(4)=VR(8);
X1(1)=VR(9); X1(2)=VR(lO): XI(3)=VR(ll); X1(4)=VR(12);
Y1(1)=VR(13); Y1(2)=VR(14); Y1(3)=VR(15); Y1(4)=VR(16);
X2(1)=VR(17); X2(2)=VR(18); X2(3)=VR(19); X2(4)=VR(20);
Y2(1)=VR(21); Y2(2)=VR(22); Y2(3)=VR(23); Y2(4)=VR(24);
% Calcul de la contrainte superficielle
el=OM(l)*(X(1)*(t.A2-b 1 .A2)+2*Y(l)*bl*t).f(t.A2+b1.A2).A2+OM(2)*(X(2)*(t.Y-
b2.A2)+2*Y(2)*b2*t)J(t.A2+b22A2).A2+...
OM(3)*(X(3)*(t.A2-b3.A2)+2*Y(3)*b3*t)J(t.Y+b3.Y).Y+OM(4)*0[(4)*(t.A2-
bP2)+2*Y(4)*b4*t)J(t.A2+b4.A2)."2;
e2=OM1(1)*(Xl(l)*((t+d).Y-
blDA2)+2*Y l(l)*bl*(t+d))*co~(f)J((t+d).~2+bl.Y).~2+OM1(2)*(X1(2)*((t+d).~2-
b2?2)+2*Y 1(2)*b2*(t+d))*co~(f)J((t+d)~~2+b2.~2). A2+..D
OM1(3)*(X1(3)*((t+~l).~2-
b3."2)+î*Y 1(3)*b3*(t+d))*co~(f)J((t+d).~2+b3.Y).~2+OM1(4)*(X1(4)*((t+d).~2-
bV2)+2*Y l(4)*b4*(t+d))*~os(f)J((t+d).~2+b4.~2). A2;
e3=OMl(l)*(-Y l(l)*((t+~i).~2-
bl.Y)+2*X1(1)*bl*(t+d))*~in(f)l((t+d)."2+b1.~2).~2+OM1(2)*(-Y 1(2)*((t+d).Y-
b2."2)+2*X 1(2)*b2*(t+d))*~in(f)J((t+d).~2+b2.~2).~2+..~
OM1(3)*(-Y1(3)*((t+d)."2-
b3 DA2)+2*X1(3)*b3*(t+d))*sin(f)./((t+d).A2ib3 .Y).A2+0M1(4)*(-Y1(4)*((t+d).Y -
b4.^2)+2*X l(4)*b4*(t+d))*sin(f)J((t+d).A2ib4.A2).A2;
 e4=OM2(1)*(X2(1)*((td).A2-b1.A2)+2*Y2(1)*b1*(t-d))*cos(f)J((t-
d)."2+b1."2)."2+OM2(2)*(X2(2)*((ta) ."2-b2."2)+2*Y 2(2)*b2*(t-d))*~0~(f)
./((tg
d).A2+b2.A2).A2+...
OM2(3)*(X2(3)*((t-d) ."Z-b3 . A 2 ) + 2 * Y Z ( 3 ) * 7
d)."2+b3 ."2)."2+OM2(4) *(X2(4)*((t4) -"2-b4?2)+2*Y2(4)*M*(t-d))*~0~(f) ./((t-
d).A2+b4.A2).A2;
e5=OM2(1)*(-Y2(1)*((td).A2-b 1.A2)+2*X2(1)*bI*(t-d))*sin(f)J((t-
d)."2+b l .A2).A2+OM2(2)*(-Y2(2)*((t-d).A2-b2-A2)+2*~(2)*b2*(t-d))*~h(QJ((t-
d)."2+b2.Y)?2+...
OM2(3)*(-Y2(3)*((t-d).A2-b3.A2)+2*X2(3)*b3*(t4))*~h(QJ((t-
d)X + b 3 .AZ).A2+OM2(4)
*(-Y2(4) *((t-d) .A2-b4.A2)+2*X2(4)*b4* (t-d))*sin(f)J((t-
d).A2+b4.AZ).A2;
E 11=-(cos(beta*pill80)J(t-to))*W/(2*pi)+e l+eZe3+e4+e5;
e6=0M(1)*(2*X(l)*b 1*t-Y(1)*(t?2-
b 1.A2))J(t.A2+b1.A2). A2+OM(2)*(2*X(2)*b2*t-Y(2)*(tY-b2.". ..
0M(3)*(2*X(3)*b3*t-Y(3)*(t.A2-b3.Y))J(t.A2+b3.A2).A2+OM(4)*(2*X(4)*b4*t-
Y(4)*(t."2-b4."2)) J(t,A2+b4.A2).Y;
e7=OM1(1)*(-Y l(l)*((t+~i).~2-
b 1.A2)+2*X1(l)*bl*(t+d))*cos(f)J((t+d).A2+b P2)."2+OM1(2)*(-Y 1(2)*((t+~l).~2-
b2!2)+2*X l(2)*b2*(t+d))*cos(f)J((t+d) ?2+b2."2) !2+..
OMl(3)*(-Y 1(3)*((t+d)."2-
b3 !2)+2*X 1(3)*b3*(t+d))*~os(f)J((t+d).~2+b3 ."2) ?2+OM 1(4)*(-Y 1(4)*((t+d) ."2-
WA2)+2*X1(4)*b4*(t+d))*cos(f)J((t+d). AZ+b4.A2) ?2;
e8=oMl(l)*(Xi(i)*((t+d).A2-
b 1."2)+2*Y l(l)+bl*(t+d))*~in(f)J((t+d).~2+b 1.A2).Y+OM1(2)*(X1(2)*((t+d).Y-
b2. A2)+2*Y1(2)*b2*(t+d))*sh(f) J((t+d) .A2+b2.A2).A2+...
OMl(3)*(X l(3)*((t+d)."2-
b3 ?2)+2*Y 1(3)*b3*(t+d))*sin(f)
-2 ' . )2
A.3b+2
A. )di t ( (J
b4.A2)+2*Y 1(4)*b4*(t+d))*sin(f) J((t+~i).~2+b4.~2).Y;
e9=OM2(1)*(-Y2(1)*((td).A2-bl.Y)+2*X2(1)*b 1*(t-d))*cos(f)J((t-
d). "2+b l."2) ."2+0M2(2)*(-YZ(2) *((t-d) . " 2 - b 2 A
d).A2+b2.A2).A2+...
Om(3)*(-Y2(3)*((t-d).Y-b3.A2)+2*X2(3)*b3*(t-d))*~~~(QJ((t-
."2) ."Z+On(12(4)*(-Y
d).~2+b3 2(4) *((
d).A2+b4.Y).A2;
 e1WM2(1)*~(l)*((t-d).A2-bl.A2)+2*Y2(1)*b 1*(td))*sin(f)J((t-
d)."2+b 1-9)."2+OM2(2)*(X2(2)*((t-d).A2-b2.r)+2*Y2(2)* b2*(t-d)) *sh(f)./((t-
d).A2+b2.A2).A2+...
oM2(3)*(X2(3)*((t-d).A.2-b3-"2)+2*Y2(3)*
d).^2+b3 .A2).AZ+OM2(4)*(X2(4) *((t-d).A2-b4.A2)+2*Y 2(4)*M*(t-d))*sb(f) J((t-
d).A2+b4.A2).A2;
E12=-(sin(beta*pi/l80)J(t-to)) *W/(Z*pi)+e6+e7+e8+e9-e10;
Dl=al*(tA2-b 1A2)l(t.A2+b1A2).A2+a2*(t.A2-b2A2)J(t.A2+b2A2).A2+a3*(~A2-
b3A2)J(tA2+b3Y).A2+a4*(t.A2-b4'L2) J ( ~ ~ 2 + b 4?2; ~2)
D2=al *((t+d) !2-b 1A2)J((t+d)?2+b 1A2).A2+a2*((t+d) !2-
b2A2)J((t+d).Y +b2A2)?2+a3 *((t+~i)?2-b3~2)J((t+à). Yl+b3 AZ).A2+a4*((t+d). A2-
MA2)J((t+d) .A2+b4A2).A2;
D3=al *((t-d)?Z-b 1A2)J((t-d).A2+b 1A2).A2+a2*((t-d).A2-b2A2)J((t-
d).A2+b2A2).A2+a3*((t-d) ?2-b3 A2)J(J
d).A2+b4"2).A2;
D4=2*(t+d).*(al *b 1J((t+~l).~Z+b lA2).A2+a2*b2J((t+d).A2+b2A2).A2+a3*b3J...
((t+d).A2+b3"2). A2+a4*b4J((t+d) .A2+b4A2).A2);
D5=2*(t-d).*(al *b 1J((t-d)!2+b 1Y).A2+a2*b2J((t-d).A2+b2A2).A2+a3*b3 ./((t-
d) .l2+bY2) .A2+a4*b4./((t-d) !2+MA2) !2);
D 1 l=c*(1+2*cos(f))-D 1-@2+D3)*cos(f)+(DeD5)*sin(f);
D6=2*t.*(al *b 1J ( ~ ~ 2 + lA2)
b .A2+a2*b2J(t.A2+b2A2).A2+a3*b3J(t.A2+b3A2).A2+...
a4* b4J(t. A2+b4A2)!2);
D 12=-D6-@4+D5)*cos(f)-@2-D3)*sin(f);
SIGMA=4*(E11.*D 11+E12.*D12)J(D 1leA2+D12!2);
SIGMA=SIGMAJ4;
% Calcul des coordonnées cartésiennes pour la contrainte superficielle pour trois dents
% SIGMAc est la contrainte en compression
% SIGMAt est Ia contrainte en tension
DlM=length(XTC);
ii=l;J=l;
for i=2:@IM)
if XTC(i)<=û
H 1c=sqrt((XïC(i- 1)-XTC(i))"2+(YTC(i)-YTC(i- 1))A2);
XSIGMAc(ii)=XTC(i)-SIGMA(i)*(YTC(i)-YTC(i-
1))JH1c ;
YSIGMAc(Ü)=YTC(i)-SIGMA(i)*C(i-1)-XTC(i))JH 1c;
ü=ü+I ;
 else
H1t=sqrt((XTC(i)-XTC(i-
l))A2+CYTC(i-1)-YTC(i))A2);
XSIGMAt(jj)=XTC(i)-SIGlAA(i)*(YTC(i-1)-YTC1~
YSIGMAt(ij)=YTC(i)-SIGM.A(i)*mC(i)-XTC(i- 1)) .AI1t;
jj=jj+l;
end
end
% Axes des coordom&s
xa=[O O]; ya=[-.5 3 1; xal=[-6 61; yal=[O
O];
~IO~~SIGMA~,YSIGMAC~J~SIGMA~O,YSIG~~O~C,YTC~SIG
X S 1 G M A G Y S I G U ~ x ~ y ~ ' -ya
. ' 1d,Y)
,
title(E'Contrainte superficielle','SIGMAtk ',num2str(rnax(4*SIGMAû)),', SIGMAt=
',num2str(max(4*SIGMA))])
pause
96 Génération d'une dent d'engrenage interne

cIear
96 CARACTÉRISTIQUES DE L'ENGRENAGE
% Parametres connus de L'engrenage interne
96 Lnombre de dents
% m=rnoduie
96 PA=angle de pression
% N=nombre de points que i'on désire pour tracer I'involute
"P. CARAC~ÉRISTIQUES DU COUTEAU-PIGNON
% Parametres connus du couteau pignon:
% Zc=nombre de dents
% mc= modul
9% PAc= angle de pression
% M-liombre de points que l'on désire pour tracer la tmchoïde
=input ('Z=')

mc=input ('mc=')
PAc=input (TAC=')
M=input('M=')
% CALCUL DU PROFIL DE RACCORD (PARTIE EN TROCHOÏDE),
% PARTIE G~NERÉEPAR LE COUTEAU
% Pas diarnetral
Pc=l/mc;
% Addendum
ac=l/Pc;
% Dedendum
bc= 1.25Pc;
Sb Diametre primitif
Dc=Zc/Pc;
% Diametre interne
% Diametre externe
Dec=Dc+2*bc;
% Rayon de base
Rbc=Dd2*cos(PAc*pi/180);
% CALCUL DU PROFTL ACTIF (PARTIE EN DÉVELOPPANTE DE CERCLE)
95 Pas diarnetral
P=l/m;
% Addendum
a=1/P;
% Dedendum
b= l.25/P;
% Diametre primitif
D=m;
% Diametre exterieur
De=D+Z*b;
% Diametre interne
Di=D-2*a;
4% Rayon de base
Rb=D/2*cos(PA*pill80);
% Calcul de I'entraxe
C=(Z-Zc)/(2*P);
% Calcul de la fonction "ulvolute"pour l'angle de pression au sommet du couteau
Beta=2*Rbc/Dec;
PAec=acos(Beta);
inv(PAec)=tan(PAec)-PAec;
% L'épaisseur circulaire de la dent sur fe cercle primitif
T=.S*pi/P;
Q Calcul du rayon d'intersectionde la partie en développante de cercle et de la partie en
trochoïde
W=(C*sin(PA*pi/l80)+@e~/2)*sin(PAec))~2+Rb~2;
Rf=sqrt(W);
% Calcul de la fonction "involute"pour l'anglede pression PA
inv(PA)=tan(PA*pifl80)-PA*pi/l80;
% Calcul de la développante de cercle
for i=l :N
 R2(i)=Rb+(Rf-Rb) *i/N;
alfa(i)=Rb/R2(i);
PA1(i)=acos(aKa(i));
INV(i)=(tan(PAi (i)))-PAl (i);
Teta(i)=(T'/D-inv(PA)+INV(i));
X(i)=R2(i)*sin(Teta(i));
Y(i)=FU(i)*cos(Teta(i));
end
KI';Mt;
45 Calcul des coordonnées cartdsiennes pour le rayon interne
Ri=Di/2;
fi=(acos(Rb/Ri));
INVi=(tan(fi))-fi;
teta=(T/D-inv(PA)+Wi);
x=Rï*sin(teta);
y=Ri*cos(teta);
8 Calcul des coordonnées cartésiennes pour le rayon de base
teta l=(T/D-hv(PA));
xb=Rb*sin(teta1);
yb=Rb*cos(tetal );
% ' o c d des coordonnées cartésiennes pour le rayon Rf
alfaf=Rb/Rf;
fif=(acos(alfaf));
INVf=tan(f@-ffi,
teta2=Cr/D-inv(PA)+INVf);
xf=Wsin(teta2);
yf=Rf*cos(teta2);
% Calcul de la trochoïde
for j=l :M
rt(j)=Rf+@e/2-Rf)*j/M;
E l=rt(j)A2-CA2-DecA2/4;
E2=C*Dec;
E3=E 1E2;
% L'angle de rotation du couteau
epsilonc=acos(E3)*18O/pi;
46 L'angle de rotation de I'engrenage
 epsilonZ=@dD)*epsilonc;
% Calcd de l'angle vectoriel de la trochoïde "anglet"
E4=-@ec/2)*sin(epsilonc+pi/180);
E5=E4Jrt(j);
ES=(asin(ES))* i Solpi;
angIet=E6epsilon2;
% L'épaisseur circulaire de la dent-pignon au sommet du couteau, en radians
Tec=T/Dc+inv(PA)-hv(PAec);
% Calcul de I'angle vectoriel de la trochoïde "tetasecond" par rapport a l'axe de symetrie de
% la dent
deIta=(pin-(ZclZ) *Tec)* 18O/pi;
tetaseconde=delta-andet;
Xt(j)=rt(j)*sin(tetaseconde*pill80);
Yt(j)=rt(j)*cos(te*iseconde*pi/lûO);
end
lm]';rytl';
% CALCUL DU FLANC GAUCHE DE LA DENT-PARTIE DÉVELOPPANTE
% L'épaisseur de la dent sur le cercle de base
Fi=acos(Rb/Rb);
INVb=tan(Fi)-Fi;
Tb=2*Rb*(T/D-invPA)+INVb);
X(l)=Ri*sin(teta);
X(N)=Rf*sin(teta2)
for i=l:N
dx(i)=X(i)-xb;
Xr(i)=xb-Tb-dx(i) ;
end
[Wl)I; kwN)I;[Xrl';
% CALCUL DU FLANC GAUCHE DE LA DENT-PARTIE TROCHOIDE
% L'épaisseur de la dent au niveau du Rf
inv(PAc)=tan(PAc*pi/l80)-PAc*pi/l80;
T~=~*W((T/D)-~V(PAC)+INV~);
o=Tb+2*dX(N);
for j=l:M
dXt(i)=Xt(i)-xf;
Xtr(j)=xf-o-&tu);
end
Fm
Xd=[x Xr(l)];Yd=@ y];
plot(Xt,-Yt&Y a,-Y,Xtr,-Yî,Xd,-Yd)
titie('ENGRENAGE INTERNE DROIT')
pause
disp('N0MBRE DE DENTS (a')
W?(CrnD
disp('N0MBRE DE DENTS POUR LE COUTEAU (Zc)')
disp(r[Zcll)
disp('ANGLE DE PRESSION (PA)')
dispuPAII)
disp(FACTEUR ADDENDW)
&SP([[111)
disp('FACTEURDEDENDUM')
disp(K 1.2513)
disp(TRANSITI0N DEVELOPPANTE-TROCHOIDE (xf,yf)')
disp([ rxfl MI)
disp(TROCH0IDE-FLANC DROIT DE LA DENT (rt(i),Xp,Yp)')
disp([ f d ' WI' C-Ytl' 1)
disp('1NVOLUTE-FLANCDROïï DE LA DENT (R2(i),X,Y)')
disp(f WI' IX]' [-YIt 1)
disp('INV0LUTE-FLANC GAUCHE DE LA DENT (R2(i),Xr,Y)')
~isp(CW1' ml'C-YI'I)
disp('TEt0CHOIDE-FLANCGAUCHE DE LA DENT(rt(j),Xtr,Yt)')
~ s ~ ( [ rml'
w [-ym
% LA DENT D'ENGRENAGE INTERNE TRACEE PAR LA METHODE CLAUDE
% GOSSELIN, pour le cas: 2545 dents, moduled, angle de pression=20
XCl=[-5.6278 -5.55349 -5.48364 -5.42246 -5.37352 4.36969 -3.56359 -2.93555 -2.46488
-2.13003 - 1.90865 1.90852 2.12992 2,46480 2.93547 3 S6348 4.36952 5.37326 5.422 19
5.48338 5.55322 5.627531;
YC 1=[-94.6768 -94.65221 -94.6 1 1 16 -94.55537 -94.48757 -92.63499 -90.94164 -89.42254
-88.09 1 12 -86.95907 -86.03636 -86.03687 -86.9596 1 -88.09 164 -89.42300 -90.94199
-92.63520 -94.48759 -94.55537 -94.6 1 1 18 -94.65224 -94.6768 11;
plot(Xt,-Yu,-Y&-Y,Xbt,-Yt,Xd,-Yd,XC 1,YC1,'o')
title('ENGRENAGE INTERNE-PROFIL &EL')
% L'assemblage des deux parties de la dent

clear
% CARA~~ÉRISTIQUES DE L'ENGRENAGE
% Parametres connus de l'engrenage interne
% Z=nombre de dents
% m=module
% PA=angie de pression
% N=nombre de points que l'on désire pour tracer l'involute
% CARACTÉRI:STIQUESDU COUTEAU-PIGNON
% Parametres connus du couteau pignon:
% Zcaombre de dents
% mc=moduI
% PAc=angle de pression
% Maombre de points que I'on désire pour tracer la trochoïde
-put ('Z-)
m=input (lm=')
PA=input ('PA=')
N=input ('N=')
zc=input ('('Zc=')
mc=input ('mc=')
PAc=input ('PA&)
M=input('M=')
% CALCUL DU PROFIL DE RACCORD (PARTIEEN TROCHOIDE),
96 PARTIE GÉNERÉE PAR LE COUTEAU
96 Pas diametral
Pc=l/mc;
% Addendum
ac=l/l?c;
% Dedendum
bc= 1.îS/Pc;
% Diametre primitif
Dc=Zc/Pc ;
% Diametre interne
 Dic=Dc-2*ac;
% Diametre externe
Dec=Dc+Z*bc;
46 Rayon de base
R b oW2*cos(PAc*pill80);
95 CALCUL DU PROFIL ACTIF (PARTIE EN DI%ELOPPANTE DE CERCLE)
% Pas diametrd
P=I/m;
% Addendum
a=l/P;
% Dedendm
b= 1Z/P;
8 Diametre primitif
D=Z/P;
% Diametre exterieur
De=D+2*b;
% Diametre interne
Di=D-2*a;
% Rayon de base
Rb=D/2*cos(PA*pi/l80);
% Calcul de l'entraxe
C=(Z-Zc)/(2*P);
% Calcul de la fonction "involute"pour i'angle de pression au sommet du couteau
Beta=2*RbclDec;
PAec=acos(Beta);
inv(PAec)=tan(PAec)-PAec;
4% L'epaisseurcirculaire de la dent sur le cercle primitif
T=S*piP;
% Calcule du rayon d'intersection de la partie en développante de cercle et de la partie en
% trochoïde
W=(C*sin(PA*pi/l80)iDec/2) *sin(PAec))"2+RbA2;
Rf=sqrt(W);
9% Calcul de la fonction "involute"pour l'angle de pression PA
inv(PA)=tan(PA*pi/l80)-PA*pi1180;
% Calcul de la développante de cercle
for i=l :N
 R2(i)=Rbt(Rf-Rb)*iM,
a.Ka(i)=RbM(i);
PA 1(i)=acos(aKa(i));
INV(i)=(tan(PA 1(i)))-PAI(i);
Teta(i)=(T/D-inv(PA)+INV(i));
X(i)=RZ(i)*sin(Teta(i));
Y(i)=R2(i)*cos(Teta(i));
end
BI';M';
% Calcul des coordonnées cartésiennes pour le rayon interne
Ri=Di/2;
fi=(acos(Rb/Ri));
Wi=(tan(fi))-fi;
teta=(Tn>-Uiv(PA)+Wi);
x=Ri*sin(teta);
y=fi*cos(teta);
% Calcul des coordonnées cartésiennes pour le rayon de base
teta1=(TD-inv(PA));
xb=Rb*sin(tetal);
yb=Rb*cos(tetal);
% Calcd des coordonnées cartésiennes pour le rayon Rf
a.ifaf=Rb/Rf;
fif=(acos(aifaf));
INVf=tan(fif)-fif;
teta2=(T/D-inv(PA)+wvf);
xf=Wsin(teta2);
yf=Rf*cos(teta2);
% Calcul de Ia trochoïde
for j=l:M
rt(j)=D/2+(Del2-D/2)*jM
El =1-t(j)~2-C~2-Dec~W4;
EZ=C*Dec;
E3=E 1E2;
% L'angle de rotation du couteau
epsilonc=acos(E3)*180/pi;
% L'angle de rotation de l'engrenage
 epsilon2=~)*epsiIonc;
% Calcul de l'angle vectoriel de la trochoïde "anglet"
E4=@ed2)*sin(epsiIonc*pi/L 80);
ES=EQJrt(j);
E6=(asin@S))* i Solpi;
angIet=E6epsilon2;
% L'épaisseur circdaire de la dent-pignon au sommet
% du couteau, en radians
Tec=T/Dc+inv(PA)-invCpAec);
46 CaIcul de l'angle vectoriel de la trochorde "tetasecond"par rapport a l'axe de symeme de
% Ia dent
deIta=(pin-(zcn)*Tec)* 18O/pi;

end
Cxtl'; IYtl';
% CALCUL DU FLANC GAUCHE DE LA DENT-PARTIE DÉVELOPPANTE
9% L'épaisseur de la dent sur le cercle de base

for i= 1:N

end

titIe(rASSEMBLAGE DES PROFILS','Xf=',numîstr(-xf),' Yf=',num2str(-yf)])

pause
disp(TRANSITI0N DEVELOPPANTE-TROCHOIDE (Xf,YQt)
disp([[-xfl [-Yfll)
% L'innuence du type de couteau sur le profil de raccord

cIear
% CARACTÉRISTIQUES DE L'ENGRENAGE
% Paramètres COMUS de l'engrenage interne
% Z=nombre de dents
% m=rnodule
% PA=angle de pression
% N=nombre de points que l'on désire pour tracer I'involute
% CARA~RISTIQUESDU COUTEAU-PIGNON
% Paramètres connus du couteau pignon:
% Zc=nombre de dents
% rnc=modul
% PAc-angle de pression
% M=nombre de points que l'on désire pour tracer la trochoïde
z==iIlput('Z)
m=input ('In=')
PA=input ('PA=')
N=input (Tu=')
Zc=input ('Zc=')
mc=input ('mc=')
PAc=inpu t ('PAc=')
M=input ('M=')
zc l=ïnput ('Zc 1=')
9% CALCUL DU PROFIL DE RACCORD (PARTLE EN TROCHOIDE),
% PARTIE GÉNERÉE PAR LE COUTEAU
% Pas diametrd
Pc=1/mc;
% Addendum
ac=l/Pc;
% Dedendum
bc= 1.25Pc;
% Diametre primitif
Dc=Zc/Pc ;
% Diametre interne
Dic=Dc-2*ac;
% Diametre externe
Dec=Dc+2*bc;
96 Rayon de base
Rbc=Dc/2*cos(PAc~il180);
96 CALCUL DU PROFILACIIF (PARTIEEN DÉVELOPPANTE DE CERCLE)
% Pas diametrd
P=l/m;
% Addendum
a=l/P;
% Dedendum
b= 1.25P;
% Diametre primitif
D=Z/P;
% Diametre exterieur
De=D+2*b;
% Diametre interne
Di=D-2*a;
% Rayon de base
Rb=D/2*cos(PA*pill80);
% Calcul de l'entraxe
C=(Z-Zc)/(2*P);
% Calcul de la fonction "involute" pour l'angle de pression au sommet du couteau
Beta=2*Rbc/Dec;
PAec=acos(Beta);
inv(PAec)=tan(PAec)-PAec;
% L'6paisseur circulaire de la dent sur le cercle primitif
T=.S*pilP;
% Calcule du rayon d'intersection de la partie en développante de cercle
% et de la partie en trochoïde
W=(C*~in(PA*pi/l80)+(D&)*sin(PAec))~2+RbY;
Eu=sqrt(W);
% Calcalcul de la fonction "involute" pour l'angle de pression PA
inv(PA)=tan(PA*pil180)-PA*pi/ 180;
% Calcul de la développante de cercle
for i=I :N
R2(i)=Rb+(Rf-Rb)*iM;
alfa(i)=Rb/R2(i);
PA1(i)=acos(alfa(i));
INV(i)=(tan(PAi(i)))-PA l(i);
Teta(i)=fl/D-inv(PA)+W(i)) ;
X(i)=R2(i)*sin(Teta(i));
Y(i)=RZ(i)*cos(Teta(i));
end
[XI'; M';
% Cdcd des coordomées cart6siemes pour le rayon interne
Ri=Di/2;
fi=(acos(Rb/Ri)) ;
INVi=(tan(fi))-fi;
teta=(T/D-inv(PA)+INVi);
x=Ri*sin(teta);
y=Ri*cos(teta);
% Calcul des coordonnées cartésiennes pour le rayon de base
tetal=(T/D-inv(PA));
xb=Rb*sin(teta1);
yb=Rb*cos(tetal);
96 Calcul des coordonnées cartésiennes pour le rayon Rf
alfaf=Rb/Rf;
fif=(acos(dfaf));
INVf=tan(fif)-f~
teta2=Cr/D-inv(PA)iINVf);
xf=Wsh(teta2) ;
yf=Wcos(teta2);
% Calcul de Ia trochoïde
for j=l:M
rt(j)=D/2i(De/2-D/2) *jM;
E 1=rt(j)A2-CA2-DecA2/4;
E2=C*k;
E3=E IE2;
% L'angle de rotation du couteau
epsilonc=acos(E3)* 18O/pi;
% L'angle de rotation de I'engrenage
epsilon2=@c/D)*epsilonc;
% Calcul de l'angle vectoriel de la trochoïde "anglet"
E4=-@ed2)*sin(epsiIonc+Pi/l80);
E5=E4Jrt(j);
ESC(asin(E5))* 18O/pi;
anglet=Ekpsilon2;
% L'épaisseur circulaire de la dent-pignon au sommet du couteau, en radians
Tec=T/Dc+inv(PA)-inv(PAec);
% Calcul de I'angle vectoriel de la trochoïde "tetasecond"par rapport a l'axe de symetne de
la dent
deita=(pi/Z-(zcn)*Tec) * 18O/pi;
tetaseconde=anglet-delta;
XtQ)=rt(j)*sin(tetaseconde*pi/ 180);
Yt(j)=rt(j)*cos(tetaseconde*pill80);
end
ml'; ml';
96 CALCUL DU FLANC GAUCHE DE LA DENT-PARTE DÉVELOPPANTE
% L'épaisseur de la dent sur le cercle de base
Fi=acos (RbfRb);
rNvb=-tan(Fi)-Fi;
Tb=2*Rb*(T/D-inv(PA)+nUVb);
X(l)=Ri*sin(teta);
X(N)=Wsin(teta2);
for i= l :N
dx(i)=X(i)-xb;
Xr(i)=xb-Tb-dx(i) ;
end
l)l;kw')I;
Dl';
pWX&-YtXr,-Y)
% PROFIL DE RACCORD GÉNERÉ PAR UN AUTRE COUTEA
% Diametre primitif du couteau "Dc1"
Dcl=Zcl/Pc;
% Diametre interne du couteau "Dic1"
Dic l=Dc 1-2*ac;
% Diametre externe du couteau "Decl"
Dec l=Dc l+2*bc;
46 Rayon de base du couteau "Rbc1"
Rbc l=Dc lR*cos(PAc*pÿ180);
% Calcul de I'entraxel
C1=(Z-Zci)/(Z*P);
% Calcul de la fonction "involute" pour I'angle de pression au sommet du couteau1
Betal=2*Rbcl/Dec 1;
PAec l=acos(Betal);
inv(PAec l)=tan(PAec1)-PAec 1;
% Calcule du rayon d'intersection de la partie en développante de cercle 1
% et de la partie en trochoïde 1
W i=(Cl*sh(PA*pi/l8O)+(Dec lR)*sin(PAec1))A2+RbA2;
Rf l=sqrt(Wl);
5% Calcul de la développante de cercle 1
for i= 1:N
R21(i)=Rb+(Rf 1-Rb)*i/N;
alfa1(i)=Rb/RZ1(i);
PA1 l (i)=acos(alfal(i));
INV1(i)=(tan(PA11(i)))-PA11(i);
Tet a 1(i)=(T/D-inv(PA)+INV 1(i));
X 1(i)=R21(i)*sin(Tetal(i));
Y 1(i)=R2 1(i)*cos (Teta1(i));
end
Cxll'; Dl]';
% Calcul des coordonnées cartésiennes pour le rayon Rf
alfaf1= R b M ;
£31 =(acos(alfaf1));
INVl f=tan(fif1)-NI;
teta2 1=(Tm-inv(PA)+INV1f);
xfl=Rfl *sin(tetaZ1);
yfl=Rf l *cos(teta21);
% Calcul de Ia trochoïdel
for j= l :M
rt 1(i)=D/2+@e/2-D/2)*j/M;
E 1l=rt 1(j)--C I "2-Dec l W4;
% L'angle de rotation du couteau
epsiloncI =ZEOS(E~ 1)* 180/pi;
% L'angle de rotation de l'engrenage
epsilon2l=(Dc l/D)*epsilonc1;
% Calcul de l'angle vectoriel de la trochoïde "anglet"
E4I=@ec 1/2)*sin(epsilonc1+pi11 80);
ESl=E4lJrt10);
E6 l=(asin(ES l)* MO/pi;
angletl=E61epsilon2 1;
Q L'épaisseurcirculaire de la dent-pignon au sommet du couteau, en radians
Tecl=T/Dc 1+inv(PA)-inv(PAec 1);
% Calcul de I'angle vectoriel de la trochoide "tetasecond"par rapport a l'axe de symetrie de
% la dent
delta l=(pÿz-(Zc ln)% 1)* 180/pi;
tetaseconde1=anglet 1-delta1;

end
[ x t il'; [Yt il';
for i=l :N
dx 1(i)=X l (i)-xb;
X 1r(i)=xb-Tb-dx l (i);
end

title('INm,UENCE DU COUTEAU S U R LE RACCORD')

% Transition d'une dent externe droit vers une dent interne droit (transition d'une
% transformation conforme à quatre termes vers une transformation conforme à sept
% fermes)

clear
% COEFFICIENTS DE LA TRANSFORMATION CONFORME
A=p 1-8225.02 8.242 2.8 12 -891 -2078 -09811;
B=[14.8 12.8 4.6 1.7.58 -145 .073];

% Coordonnée curviligne liée à la surface du domaine

t=-27: .O2:27;
% Transformation conforme,quatre termes
X4=c*t+al*tJ(t?2+b 1A2)+a2*tJ(t.A2+b2A2)+a3*t./(~~2+b3~2)+a4*tJ(t.~2+b4~2);
Y4=al*b lJ(t.A2+b1A2)+a2*b2J(t-A2+b2A2)+a3*b3J(t-A2+b3A2)+a4*I .7J(t.A2+b4A2);
46 Transformation conforme, cinq termes
X5=X4+a5 *tJ(t.A2+b5A2);
Y5=Y4+~iS*b5J(t.~2+b5~2);
% Transformation conforme, six termes
X6=X5+a6*tJ(~~2+b6~2);
Y6=YS+a6*b6J(t .A2+bb6A2);
% Transformaton conforme, sept termes
X7=X6+a7*~/(t.~2+b7Y);
Y7=Y6ta7*b7J(t.A2+b7A2);
plot(X4,Y4,'r',X5,Y5,*rt$C6,Y6~'b',X7,Y7,%')
title(TRANSm0N D UNE DENT EXTERNE VERS UNE INTERNE')
end
% Calcul des contraintes superficielles engrenage interne

clear
% Parametres connus:
% force appiiquke 'W"
9% point de l'application de la force "to"
9% angle d'inclinaison de Ia force "beta"en degrées
W=input(W=')
to=input('to=')
beta=input('beta=')
% Coordonnée curviligne liée à la surface du domaine
tl=45:.1:-.16;t2=0.185:.1:45;
t=CM It213;
% Coefficients de la transformation conforme
A=[3 1-8225.02 8.242 2.8 12.891 .2078 -09811;
B=[14.8 12.8 4.6 1.7.58 -145 .073];
C=[.1241;
al=A( 1); a2=A(2); a3=A(3); a4=A(4); &=A@); a6=A(6); a7=A(7);
bkB(1); b2=%(2); b3=B(3); b4=B(4); b5=B(5); b6=B(6); b7=B(7);
1);
CS(
% Calcul de la dent par la MPC
XTCs*t+a 1 *tJ(t A2+b1.A2)+a2*t/(tA2+b2.A2)+a3 *t./(t.A2+b3.A2)+a4*t./(t-A2+b4.A2)+
a5*tJ(t. A2+b5.A2)+a6*~(t. A2+b6.A2)+a7*t;
YTClal*b lJ(teA2+b1 .Y)+a2*b2J(tA2+b2.A2)+a3*b3J(t.A2+b3.A2)+a4*b4./(t.A2+b4."2)+
a5*b5J(t.A2+b5.A2)+a6*b6 J(~~2+b6.~2)+a7*b7J(t.~2+b7.~2);
% CALCUL DE LA CONTMINTE SUPERFICIELLE
% Calcul de la dérivée de la transfonnation conforme
for k=1:7
for j=l:ï
E(k)=[al J(B(k)+b 1).A2+a2.1(B(k)+b2).A2+a3 ./(B(k)+b3) ?2+
a4.f(B(k)+b4).A2+aS J(B(k)+b5)."2+a6J(B(k)+b6).A2+a7.(k)+b7).A2];
OMp(k)=t+s~('Oc));
% Cdcui de ALFA(j,k)
OM(k)=A(k) ./OMP(k);
 ALFA(j,k)=OM(k)/(((sQ)+B (k))X)
;
% CalcuI des coefficients nkis Ru), Sÿ)
R(j)=((to*cos(beta*pi/I 80)+Bü).*sin(beta*pi/I 80))./(toA2+Bü)."2)) *
W/(Z*pi);
S(i)=((to*sin@eta"pi/180)-B (j)-*cos(beta*pi/80)) J(toA2+B (i)."2))*
W/(î*pi);
MATALFAQ&)=[ALFA(jJc) J;
end
end
for k=l:7
for j=l:7
if k==j
U(j,k)=l ;
else
Uÿ&)=o;
end
end
end
46 Calcul des variables réelles X(j)
for k=1:7
for j=l:7
MATX&.k)=hrIATALFA(j,k)+U(i
,k);
end
end
X=(MATX)\[RIf;
% Calcul des variables delles Yu)
for k=1:7
for j=l:7
MATY(j,k)=MATALFA(jA)-Uÿ,k);
end
end
Y=WATY)\[(-S)I f ;
% Calcul de la contrainte superficielle
FL=OM(l)*(X(l)*(t."2-b IA2)+2*Y(1)*b1 *t)J(LA2+b 1.A2).A2+0M(2)*(X(2)*(t.A2-b2A2)+
2*Y(2)*b2*t)./(t-A2+b2.A2).A2tOM(3)*(X(3)*(t.A2-b3A2)+2*Y(3)*b3*t)./(t.A2+b3.A2).A2+
OM(4)*~(4)*(t.n2-b4n2)+2*Y(4)*b4*t)J(t-A2+M.A2).A2+OM(5)*(X(5)*(t.A2-b5A2)+
2*Y(s)*b5*t)J(tA2+bS.A2).A24M(@*(x(6)*(ttA2-b6A2)+2*Y(s)*b6*t).(ttA2+b6~A2).A2+
OM(7)*m(n*(tA2-b7A2)+2*Y(7)*b7*t)J(~A2+b7.A2).A2;
F1l=(cos(betaf pi/l80)J(t-to))* W/(2*pi)+Fl;
FZdM(l)*(2*X(l)*b I*t-Y(1)*(~~2-b
1"2))./(~~2+b
1A2).A2-i-C)M(2)*(2*X(2)*b2*t-Y(2)*
./(~~2+b3~2).~2+
(tA2-b2A2))J(~A2+b2A2)oA2+OM(3)*(2*X(3)*b3*t-Y(3)*(~A2-b3A2))
OM(4)*(2*X(4)*b4*t-Y(4)*(tA2-MA2))J(~A2+MA2)oA2~M(~*(2*X(5)*b5*t-Y(5)*
* ( 2 * X ( 6 ) 2 +
(tA2-b5"2))J(t0A2+b5A2).A2+OM(Q
0M(7)*(2*X(7)*b7*t-Y(7)*(t"2-b7A2))./(t.A2+b7A2)
."2;
Fi 2=-(sin(beta*pill80)l(t-to) ) *W/(Z*pi)+F2;
T l 1=c-(al *(tA2-blA2)J(tA2+b1A2).A2+a2*(tA2-b2A2)./(~A2+b2A2).A2+a3*(t.A2-b3A2)./
(t.A2+b3A2).A2+a4*(~A2-b4A2)J(t.A2+b4A2)*A2+a5*(t.A2-b5A2)./(t.A2+b5A2).A2+a6*
.f(t.AZ+b7A2).A2);
(tA2-b6A2)./(t.A2+b6A2).A2+a7*(t.A2-b7A2)
T12=-2*t *(al *b 1./(tA2+b1A2).A2+a2*b2./(t*A2+b2A2).A2+a3*b3
J(t.A2+b3A2).A2+
a4*b4J(tOA2+MA2).
A2+a5*b5./(t.A2+b5A2)
.A2+a6*06./(~A2+b6A2).A2+a7*b7./
(~~2+bF2).~2);
SIGMA=4*(F1l.*T1 l+FlZ.*T12)J(Tl 1.A2+T12.A2);
SIGMA=SIGMA/2;
96 Calcul des coordonnées catésiennes pour la contrainte superficielle
% SIGMAc est la contrainte en compression
% SIGMAt est la contrainte en tension
DIM=length(XTC)
fi=2;jj= 1;
for i=2:@LM)
if XTC(î)<=û

else

end
end
x=[-1.8957 1.8957]; y=[9.8154 9.8l541;xf=[O O];yl=[12 -31;
~ ~ O ~ ~ C , Y T C X S T G ~ C , Y S I G ~ C , ' O ' ~ S I G1~,y~ ,-')
1,'- YSIG~~~[A~
title(rlontrainte superficielle','SIG- ',num3_ntr(&(2*SIGMA)), ', SIGMAt-
',num2str(max(2+SIGMA))1)
pause
~~SP(FORCEAPPLIQUÉE~)
&sp(C[9-811)
disp('P0INT D APPLICATION DE LA FORCE')
disp(K0-158ll)
disp('ANGLE D INCLINAISON DE LA FORCE PAR RAPPORT À L AXE x')
disp([r-16011)
XSIGMA=[BS IGMAc]' [XSIGMAt]']
YSIGMA=[CySIGMAc]' FSIGMAt]']
 INDEX
Condition limite 19,25,26
conditions Cauchy-Riemann 14
conjugué 13, 14,32,33,79
contrainte superficieue 2, 39,73,74,75,97, 102, 103
engrenage droit 10,46,59,74
équation biharmonique 17, 18
equations de Kolossov 8, 18,35,39,97
facteur g6ométrique 75,76, 103, 104
fonction complexe 14
fonction d ' A j r 17,2 1
fonction de potentiel 17,18,63
holomorphe 8, 14, 15, 16, 17, 18, 19,20,21,23,24,25,28,29,31,33, 34,64,65,66
Méthode des Potentiels Complexes 5.9, 10, 11, 16,18,21,40,44,63,78,90,91,96,97,
102, 103, 104
Modélisation 40,43,90
section critique 2,3,5, 10,43,98,99, 100, 101,102
transformation conforme 8,9, 1 1, 15, 16,2 1,22,23,33, 36,37,40,4 1,42,43,44,50,53,
57,58,59,60,63,64,66,68,75,90,91,93,97
variable complexe 7, 8, 12, 13, 14, 15, 18, 19
ILSI IAKGET (QA-3)

APPLIED I M G E .lnc
-
=-
--
.-.
--
t 653 East Main Streer
-
-
--
--
- --
Rochester. NY 14609 USA
Phone: 716/482-0300
F a 7161288-5989