PLAN DE TRAVAIL

1 5
INTRODUCTION CONCLUSION

2 PROBLEMATIQUE

4 METHODOLOGIE
& RESULTATS

3 OBJECTIF
 01

INTRODUCTION
 Processus
Le concasseur à mâchoires est un équipement
incontournable dans les industries extractive,
miniers et de production des granulats, il occupe
une position prépondérante en tant que premier
maillon de la chaine de réduction granulométrique.
Sa mission consiste à transformer des bloc de
roche pouvant atteindre des dimensions métrique
en fragments exploitables (inferieur à 200mm) pour
les étapes ultérieures des traitements (broyage,
etc.)

CONCASSEUR SANDVIK CJ815

SCHEMA CINEMATIQUE DU CONCASSEUR A
MACHOIRS SANDVIK CJ815
INTRODUCTION
 Fréquence de
service
𝑵
𝒇= = 𝟑, 𝟑𝟑 𝑯𝒛
𝟔𝟎
INTRODUCTION
CHOIX DE LA THEORIE DE
MODELISATION
(TIMOSHENKO VS EULER)

INTRODUCTION
 Théorie d’Euler-Bernoulli Théorie de Timoshenko
(Modèle classique) (Modèle Avancé)

Hypothèse : le déplacement et la rotation

Hypothèse : les sections droites restent
sont traités comme deux variables
perpendiculaires à la fibre neutre.
indépendantes.

Néglige totalement le cisaillement Intègre un facteur de correction de

transversal et l’inertie rotative. cisaillement (K=0,886).

Echec scientifique : sous-estime la

Choix validé : empêche le verrouillage en
déformation de 8 à 12% pour des
cisaillement (Shear Locking).
géométries épaisses.

CONCLUSION : Parce que L/D=5,2 (Poutre courte <10), la théorie de Timoshenko

INTRODUCTION
est la seule mathématiquement valide pour ce rotor.
 CARACTERISTIQUES
TECHNIQUES DE L’ARBRE
EXCENTRIQUE

INTRODUCTION
 Propriétés du matériau LE PARAMETRE CLE
Acier 42CrMo4 Rapport L/D=5,2
Ce faible élancement classe l’arbre dans la
Limite élastique : 750 MPa catégorie absolue des « poutres courtes »,
imposant le rejet des théories classiques.
Module de Young : 210 GPa
Masse total : 3000 Kg
INTRODUCTION
 02

PROBLEMATIQUE
Contexte :
o Charges cycliques : 200 tr/min
o Pics d’impact : x 2,3 la charge normale
o Risque : rupture par fatigue = arrêts de production

PROBLEMATIQUE
Question de
recherche « Utilisation de Timoshenko +
EF sous SolidoWorks »

?
« Dans quelle mesure un modèle EF basé
sur Timoshenko permet-il de prédire le
comportement dynamique et la tenue en
fatigue de l’arbre, et quelles sont ces
marges de sécurités ? »
PROBLEMATIQUE
 03

OBJECTIF
 ORIGINALITE
Objectifs Application de la théorie de Timoshenko (L/D=5,2)
Spécifiques • Poutres courtes ⇒ prise en compte du cisaillement

ANALYSE STATIQUE
• Déterminer les contraintes, déplacements, déformations
• Vérifier le domaine élastique
ANALYSE MODALE
• Déterminer les fréquences propres
• Vérifier l’absence de résonance
• Evaluer la stabilité dynamique

ANALYSE DE FATIGUE
• Estimer la durée de vie
• Identifier les zones critiques
OBJECTIF
 04

METHODOLOGIE
 METHODOLOGIE

1 5
ÉTABLISSEMENT DU VALIDATION
MODÈLE DYNAMIQUE NUMÉRIQUE ET
RÉSOLUTION

04

2 DISCRÉTISATION
DES ÉLÉMENTS
FINIS 4 APPLICATION DES
CHARGEMENTS
VARIABLES

3
ASSAMBLAGE
ET CONDITIONS METHODOLOGIE
AUX LIMITES
Equation générale du mouvement
(Formule de Lagrange)

𝑴 𝒒ሷ + 𝑪 𝒒ሶ + [𝐊]{𝐪ሽ = {𝐅(𝐭)ሽ

C’est l’équation qui régit les comportements

de tous les systèmes mécaniques

METHODOLOGIE
JUSTIFICATION DE LA THEORIE
DE TIMOSHENKO
• Théorie de Timoshenko (L/D=5,2)
• Facteur de Couper (k=0,886)
 1
EQUATION GOUVERNANTE POUR
UNE POUTRE DE TIMOSHENKO

 Déplacement transversale
𝛛𝟐 𝛚 𝛛𝟐 𝝎 𝛛𝝋
Contrairement à Bernoulli, cette théorie
𝝆𝑨 𝟐 − 𝑲𝑮𝑨
𝛛𝒕
−
𝛛𝒙𝟐 𝛛𝒙
= 𝒇(𝒙, 𝒕)découple le déplacement transversale et la
rotation, elle introduit deux effets
• Rotation de la section 𝝋 physiques majeurs : cisaillement
transversale (KGA) et l’inertie de rotation.
𝛛𝟐 𝚿 𝛛𝟐 𝝋
𝝆𝑰 𝟐 − 𝑬𝑰 𝟐 − 𝑲𝑮𝑨
𝛛𝝎
−𝝋 =𝟎
K=0,886
𝛛𝒕 𝛛𝒙 𝛛𝒙

2 MATRICE DE RIGIDITÉ ÉLÉMENTAIRE 𝑲𝒆

𝑲𝒆 = 𝑲𝒆𝒇 + 𝑲𝒆𝒄
METHODOLOGIE
 DISCRÉTISATION PAR LA MEF
• Interpolation linéaire (pour Timoshenko
Équations (sur un élément de longueur L)
o 𝝎 𝒙 = 𝑵𝟏 𝒙 𝝎𝟏 + 𝑵𝟐 (𝒙)𝝎𝟐
o ∅ 𝒙 = 𝑵𝟏 𝒙 ∅𝟏 + 𝑵𝟐 (𝒙)∅𝟐

Avec les fonctions de forme linéaire :

𝒙 𝑥
𝑵𝟏 𝒙 =𝟏 − ; 𝑁2 𝑥 =
𝒍 𝑙

Bernoulli (Cubique), on utilise une

interpolation linéaire où 𝝎 et ∅ sont
indépendants pour éviter le phénomène de
verrouillage en cisaillement
METHODOLOGIE
• Matrice de masse élémentaire 𝑴𝒆 𝟒 × 𝟒
𝑴𝒆 = 𝑴𝒆𝒕 + 𝑴𝒆𝒓

2 0 1 0
𝐼 𝐼 0 0 1 0
𝜌𝐴𝐿 0 𝐴 0 𝑝𝑙 0 1 0 −1
𝑀𝑡 = 𝐴 Translation 𝑀𝑟𝑒 = Rotation
6 1 0 2 0 𝐿 0 0 0 0
𝐼 𝐼 0 −1 0 1
0 0
𝐴 𝐴

• Matrice de rigidité élémentaire 𝑲𝒆 𝟒 × 𝟒

𝑲𝒆 = 𝑲𝒆𝒇 + 𝑲𝒆𝒄 𝐾𝐺𝐴 𝐾𝐺𝐴 𝐾𝐺𝐴 𝐾𝐺𝐴
− −
𝐿 2 𝐿 2
−𝐾𝐺𝐴 𝐸𝐼 𝐾𝐺𝐴𝐿 𝐾𝐺𝐴 𝐸𝐼 𝐾𝐺𝐴𝐿
En développant, on 𝐿 𝐿
+
3
−
2
− +
𝐿 6
𝐾=
trouve la matrice de 𝐾𝐺𝐴 𝐾𝐺𝐴 𝐾𝐺𝐴 𝐾𝐺𝐴
− − −
rigidité : 𝐿 2 𝐿 2
𝐾𝐺𝐴 −𝐸𝐼 𝐾𝐺𝐴𝐿 𝐾𝐺𝐴 𝐸𝐼 𝐾𝐺𝐴𝐿
+ + METHODOLOGIE
2 𝐿 6 2 𝐿 3
ASSEMBLAGE GLOBAL
ET CONDITIONS AUX LIMITES
• Assemblage
𝒏 𝒏
L’assemblage suit la correspondance entre les
𝑴 =෍ 𝑴𝒆 ; 𝑲 =෍ 𝑲𝒆 degrés de liberté locaux et globaux.
𝒆=𝟏 𝒆=𝟏 Les conditions aux limites essentielles pour
reproduire le comportement réel des paliers.
• Conditions aux limites

Les paliers sont modélisés en appui simples, de position

x=200mm et x=1800mm.
Blocage des déplacement transversaux :
𝑢 = 0, 𝜔 = 0,
𝑟𝑜𝑡𝑎𝑡𝑖𝑜𝑛 𝑙𝑖𝑏𝑟𝑒 ∶ 𝜃𝑥 ; 𝜃𝑦 𝑒𝑡 𝜃𝑧

METHODOLOGIE
Charges
variables appliquées

1 Effort de flexion
• Calcul de la vitesse
𝜋𝑋𝐷𝑋𝑁 𝜋 × 0,480 × 200
𝑉= = = 5,025 m/s
60 60

• Calcul de l’effort de flexion

200000 × 2,3
𝐹𝑓 = = 91500 𝑁
5,026 METHODOLOGIE
2 Effort de Torsion
• Calcul de la vitesse angulaire
2𝜋𝑁 2𝜋×200
𝑷𝒎𝒐𝒕𝒆𝒖𝒓 W= = = 20, 944 rad/s
C= 60 60

𝒘
• Calcul du couple de torsion
200000
C= = 9550 𝑁𝑚
20,944

• Effort de torsion équivalent F

𝐷
9549
R= = 𝑂, 240𝑚 𝐹𝑡𝑜𝑟𝑠𝑖𝑜𝑛 = × 2,3 = 91511 𝑁
2 0,240
3 Effort sur le paliers
• Calcul de l’effort sur les paliers

𝑭𝒇𝒍𝒆𝒙𝒊𝒐𝒏 × 𝑳𝒂𝒑𝒑𝒍𝒊𝒄𝒂𝒕𝒊𝒐𝒏
𝑭𝒑𝒂𝒍𝒊𝒆𝒓 =
𝒍𝒆𝒏𝒕𝒓𝒆 𝒑𝒂𝒍𝒊𝒆𝒓

𝟗𝟓𝟓𝟎𝟎×𝟎,𝟖𝟎𝟎
𝑭𝒑𝒂𝒍𝒊𝒆𝒓 = = 𝟒𝟕𝟕𝟓𝟎 𝐍 par palier
𝟏𝟔𝟎𝟎

METHODOLOGIE
 VALIDATION NUMÉRIQUE

• Etude de convergence
du maillage

Maillage Éléments Contrainte max (MPa)

Grossier 15000 4,3
Moyen 45237 4,30
Fin 120 000 4,33

L’écart entre le maillage moyen et fin est inferieur à 1%

METHODOLOGIE
 04

RESULTATS
DONNÉES DU CONCASSEUR ET DE
L’ARBRE EXCENTRIQUE
Paramètres valeurs Unités
Largueur entre paliers 1600 Mm
Diamètre nominale 480 Mm
Vitesse de rotation 2OO 200 tr/min
Module de Young 210 GPa
Coefficient de poisson 0.3
Limite élastique 750 MPa
Puissance du moteur 200 KW
Masse linéique µ 1667 Kg/m
Masse totale de l’arbre 3000 Kg
Capacité max 568 t/h
Matériaux Acier 42CrMo4
Poids volant 3752 Kg
Longueur 2500 Mm
Excentricité 20 Mm
 MODELE DE L’ARBRE
Modèle géométrique de
l'arbre excentrique
Sandvik CJ815

ARBRE EXCENTRIQUE
RESULTATS
Trois analyses complémentaires ont été menées
sur l’arbre excentrique :

• Statique : pour vérifier la résistance sous charge

maximale
• Modale : pour vérifier l’abscence de résonnance
• Fatigue : pour estimer la durée de vie
• Matériau : Acier 42 CrMo4 (Re=750MPa
• Vitesse de service : 200tr/min (3,33 Hz)

RESULTATS
 TABLEAU RESULTATS DE
L’ANALYSE STATIQUE
Grandeur Résultats Noued critique Seuil critique
11070 (congé de
Von Mises 4,30 MPa 730 MPa
raccordement)

0,060 mm 18026
Déplacement De 0,1 à 0,3 mm
(extrémité libre)

Déformation 1,67 × 10−5 Élément 2602 Limite élastique

750
• Coefficient de sécurité statique 𝑆𝑠𝑡𝑎𝑡𝑖𝑞𝑢𝑒 = ≈ 174
4,30

RESULTATS
 TABLEAU RESULTATS DE
L’ANALYSE MODALE

Indicateur valeur
Première fréquence propre 83,8 Hz
Vitesse critique 5028 tr/min
Marge de sécurité 96% (facteur 25)
Amplitude modale max 0,05mm
Risque Aucun

L’amplitude modale confirme la grande rigidité de l’ensemble. Il s’agit d’une

déformée modale normalisée, pas d’une prédiction quantitative de l’amplitude
vibratoire en service réel (force d’excitation et amortissement)

RESULTATS
 TABLEAU RESULTATS DE
L’ANALYSE EN FATIGUE
Analyse Indicateur Résultat Interprétation

Estimation
Fatigue Durée de vie 9000000 cycles
conservative

Facteur de charge
77 au nœud Aucun risque de
Fatigue en min
11070 rupture

Zones très
Facteur de charge 655700 au
Fatigue faiblement
en max nœud 65O3
sollicités
Cohérence des
résultats

«La cohérence entre les trois analyses valide

la fiabilité du modèle numérique.
La zone critique identifiée par l’analyse
statique est exactement celle où l’analyse de
fatigue prédit la durée de vie minimale.»

RESULTATS
 05

CONCLUSION
Contributions

Ce mémoire apporte une contribution originale sur plusieurs plans :

• L’application de la théorie de Timoshenko à un arbre de concasseur à
mâchoires, justifiée par le rapport L/D de 5,2
• Le développement d’une démarche d’analyse couplée associant les
analyses statiques, modale et de fatigue
• La quantification précise des marges de sécurité : statique (174),
dynamique (facteur 25) et en fatigue (facteur de charge 77)
• La mise en évaluation du surdimensionnement de l’arbre, qui pourrait
être allégé de 30 à 40% sans compromettre sa tenue mécanique

CONCLUSION
Limites

Les résultats présents doivent être interprétés

à la lumière des hypothèses suivantes :
• matériau homogène et isotrope,
• absence de balourd
• paliers modélisés comme des appuis
simples
• chargement cyclique simplifié (R=-1)
• absence d’effets thermiques et de
corrosion.

CONCLUSION
Perspectives

Ce travail ouvre plusieurs perspectives pour

des recherches futures :
• L’optimisation topologique pour alléger
l’arbre de 30 à 40%
• La validation expérimentale par jauges de
déformation et accéléromètres
• L’analyse tribologique approfondie des
surfaces de contact
• L’intégration dans un outil de naissance
prédictive dans le cadre de l’industrie 4.0
• L’extension de la méthodologie à d’autres
composants critiques du concasseur CONCLUSION
Ce mémoire a porté sur l’analyse numérique de l’arbre excentrique CJ815.
Sous les charges critiques nous avons appliqué la théorie de Timoshenko,
justifiée par le faible rapport L/D, et mener trois analyses couplées les résultats
montrent une contrainte maximal de 4,30 MPa, une première fréquence propre
à 83,3 Hz et une durée de vie de 𝟗 × 𝟏𝟎𝟔 de cycles.
Les marges de sécurité exceptionnelles :
174; 24 et 77 confirment un surdimensionnement important offrant un
potentiel d’allégement de 30 à 40% et ouvrant la voie à une maintenance
prédictive.

CONCLUSION
 MERCI POUR
VOTRE ATTENTION

CONCLUSION