0% ont trouvé ce document utile (0 vote)
5 vues21 pages

Calcul des Caractéristiques Moteur Asynchrone

Transféré par

alexandromitso6
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
5 vues21 pages

Calcul des Caractéristiques Moteur Asynchrone

Transféré par

alexandromitso6
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

REPOBLIKAN’I MADAGASIKARA

FITIAVANA - TANINDRAZANA - FANDROSOANA

MINISTÈRE DE L’ENSEIGNEMENT SUPÉRIEURET


DE LA RECHERCHE SCIENTIFIQUE

École Supérieure Polytechnique d’Antsiranana


MENTION : GENIE ELEC

PARCOURS : GENIE ELECTRIQUE

GRADE MASTER

MINI PROJET SUR DIMMENSIONNEMENT DES


MACHINE ELECTRIQUE

« CALCUL DES CARACTERISTIQUES DE SERVICE


DE MOTEUR ASYNCHRONE A CAGE »

Réalisé par : Mr. MITSO Alexandro

Prof Responsable :

Mr RABE Tsirobaka

Année Universitaire : 2024 – 2025


2
MINI PROJET SUR LA DIMMENSIONNEMENT DES
MACHINE ELECTRIQUE

TITRE : « CALCUL DES CARACTERISTIQUES DE


SERVICE DE MOTEUR ASYNCHRONE A CAGE »

Contexte :

Les caractéristiques de service d’un moteur asynchrone reflètent les


variations de :

Puissance électrique absorbé, courant absorbé, facteur de puissance,


rendement et glissement en fonction de puissance mécanique fournie, ou charge,
variant dans la gamme P2 = 0 ÷ 1,20 P2n, la tension et la fréquence étant @ leurs
valeurs nominales.

Peuvent s’y ajouter les foncions : couple = f (puissance fournie), courant


secondaire (absolue ou ramène) = f (puissance fournie).

La réalisation pratique de calcul est @ calcul est à faire dans la sphère Math CAD

3
Table des matières
1 GENERALITE .......................................................................................................................................... 6
1.1 Introduction : ................................................................................................................................ 6
1.2 Description de la machine asynchrone a cage : ............................................................................ 6
1.2.1 Définition............................................................................................................................... 6
1.3 Constitution .................................................................................................................................. 6
1.3.1 Stator..................................................................................................................................... 6
1.3.2 Rotor (à cage) ........................................................................................................................ 6
1.4 Principe de fonctionnement ......................................................................................................... 7
1.5 Caractéristiques ............................................................................................................................ 7
1.6 Applications................................................................................................................................... 8
1.6.1 Concept d’ingénieur sur la MAS a cage................................................................................. 8
1.7 Conclusion ..................................................................................................................................... 9
2 : CONCEPT DE CALCUL ET DIMENSIONNEMENT SUR LA MAS A CAGE ............................................... 10
2.1 Introduction ................................................................................................................................ 10
2.2 Les différentes constitutions des grandeurs électriques indispensable de la machine
asynchrone à cage................................................................................................................................... 10
2.2.1 Grandeurs d’alimentation ................................................................................................... 10
2.2.2 Grandeurs magnétiques...................................................................................................... 10
2.3 Grandeurs du rotor ..................................................................................................................... 10
2.4 Les grandeurs de fonctionnement et les équations fondamentaux........................................... 11
2.4.1 Modèle Équivalent. ............................................................................................................. 12
2.5 Modélisation Mathématique ...................................................................................................... 12
2.5.1 Hypothèses.......................................................................................................................... 13
2.5.2 4.2 Fonctions d'Approximation ........................................................................................... 13
2.6 Implémentation sous Mathcad ................................................................................................... 14
2.6.1 Données Nominales ............................................................................................................ 14
2.6.2 Calculs Préliminaires ........................................................................................................... 14
2.6.3 Programme Principal........................................................................................................... 15
2.6.4 Tableau de Résultats ........................................................................................................... 16
3 : RESULTATS ET ANALYSE .................................................................................................................... 17

4
INTRRODUCTION GENERALE
Dans le cadre de formation d’ingénieure sur la domaine d’industrialisation, le savoir à discerner le
fonctionnement des moteurs électriques est impérativement intrus pour renforcer la compétence
des étudiants durant la période d’adaptation d’esprit d’un génie.
Et cela, nous sommes entièrement disposé face à ce sujet ce qui nous a amené à discuter comment
calculer les caractéristiques de service de moteur asynchrone à cage.
Maintenant, on va voir apprendre un peu plus profond sur ce thème dans notre document qui est
reparti dans trois partie. En premier c’est la généralité sur machine a cage en deuxième on parlera
la résultats et analyses.

5
1 GENERALITE

1.1 Introduction :
Le moteur asynchrone à cage d’écureuil (MAS à cage) est l’un des moteurs électriques les plus
utilisés dans le monde industriel et domestique. Grâce à sa construction robuste, sa simplicité,
son faible coût et sa fiabilité, il représente la majorité des moteurs en service. Son principe repose
sur l’induction électromagnétique : un champ tournant produit par le stator induit des courants
dans le rotor en forme de cage, ce qui engendre le couple mécanique nécessaire à la rotation

1.2 Description de la machine asynchrone a cage :

1.2.1 Définition
Un MAS à cage est une machine asynchrone (souvent appelée moteur asynchrone triphasé) dont
le rotor est constitué de barres conductrices disposées en cage (comme une cage d’écureuil), reliées
aux extrémités par deux couronnes conductrices appelées bagues de court-circuit.

C’est le moteur électrique le plus utilisé dans l’industrie grâce à sa robustesse, sa simplicité et
son faible coût.

1.3 Constitution
Un MAS à cage se compose de deux parties principales :

1.3.1 Stator

 Partie fixe de la machine.

 Contient un bobinage triphasé alimenté par un réseau alternatif.

 Produit un champ magnétique tournant.

1.3.2 Rotor (à cage)

 Fait de barres conductrices en cuivre ou aluminium.

6
 Disposées parallèlement à l’axe du rotor, puis court-circuitées aux extrémités par des an-
neaux.

 Aucune connexion externe → pas de balais, pas de bagues collectrices.

 C’est ce qui lui donne le nom de rotor à cage d’écureuil.

Rotors à cage d'écureuil

1.4 Principe de fonctionnement


Le moteur asynchrone fonctionne sur le principe du champ tournant. Le rotor "glisse" par rapport
au champ statorique, ce qui génère un couple mécanique.

 Quand on alimente le stator, il crée un champ magnétique tournant.

 Ce champ coupe les barres du rotor → induction électromagnétique (loi de Faraday).

o Des courants induits apparaissent dans le rotor → ils créent à leur tour un champ
magnétique.

 Interaction entre les deux champs → le rotor est entraîné en rotation.

Le rotor tourne moins vite que le champ statorique → d’où le nom asynchrone.
La différence entre les deux vitesses s’appelle le glissement.

1.5 Caractéristiques
 Robuste (pas de balais ni de contacts frottants).

 Fiable et nécessite peu d’entretien.

 Démarrage direct possible mais avec un courant de démarrage élevé.

7
 Vitesse quasi constante (dépend légèrement de la charge).

 Rendement élevé.

1.6 Applications
 Pompes, ventilateurs, compresseurs.

 Machines-outils, convoyeurs, ascenseurs.

 Presque partout dans l’industrie et même dans les appareils électroménagers.

1.6.1 Concept d’ingénieur sur la MAS a cage

[Link] Concept électromagnétique

 Le MAS est basé sur le principe d’induction : le champ tournant du stator induit des
courants dans le rotor.

 L’ingénieur modélise ce phénomène par des équations électromagnétiques (lois de Fa-


raday, Lenz, et Maxwell).

 La différence de vitesse entre le champ et le rotor (le glissement) est essentielle pour le
couple.

[Link] Concept énergetique

 Le MAS transforme l’énergie électrique → en énergie mécanique.

 Schéma équivalent simplifié en circuit électrique (résistances et inductances).

 Rendement élevé, mais pertes dans :

o Le cuivre (résistance stator/rotor),

o Le fer (hystérésis et courants de Foucault),

o La ventilation et frottement.

[Link] Concept mécanique.

 Le rotor est conçu comme une cage d’écureuil robuste → pas de balais, pas de bagues
collectrices.

8
 Résultat : machine solide, adaptée aux environnements industriels difficiles.

 Couple mécanique proportionnel à l’interaction champ statorique ↔ champ rotorique.

[Link] Concept de commande et contrôle

Les ingénieurs ne s’arrêtent pas à « brancher et tourner ». Ils cherchent à maîtriser :

 Démarrage : réduire le fort courant de démarrage (variateurs de vitesse, étoile-triangle,


soft-starter).

 Vitesse : commandée par la fréquence de l’alimentation (via variateur de fréquence).

 Couple : adapté selon la charge (par régulation vectorielle ou directe du couple).

[Link] Concept d’ingénierie pratique

 Fiabilité : peu de maintenance.

 Économie : faible coût, rendement satisfaisant.

 Polyvalence : utilisé dans 70–80 % des machines industrielles.

 Optimisation : choix selon puissance, tension, nombre de pôles, rendement, facteur de


puissance.

En totalite sur la concept, pour un ingénieur, le MAS à cage n’est pas seulement un moteur ro-
buste et bon marché, mais un système électromécanique complet qui doit être :

 modélisé (équations),

 contrôlé (vitesse, couple, courant),

 optimisé (rendement, pertes, durée de vie),

 et adapté aux applications industrielles.

1.7 Conclusion
Le moteur asynchrone à cage d’écureuil est la machine électrique la plus répandue dans l’industrie
grâce à sa robustesse, sa simplicité de construction, son faible coût et son rendement élevé. Basé sur le
principe d’induction électromagnétique, il transforme efficacement l’énergie électrique en énergie
mécanique sans nécessiter de balais ni de collecteurs, ce qui réduit la maintenance. Sa vitesse est
pratiquement constante et son utilisation est universelle, que ce soit dans les pompes, ventilateurs,
compresseurs ou machines industrielles.

9
2 : CONCEPT DE CALCUL ET DIMENSIONNEMENT
SUR LA MAS A CAGE

2.1 Introduction
Le calcul et le dimensionnement d’un moteur asynchrone à cage d’écureuil constitue une
étape essentielle en électrotechnique pour garantir son bon fonctionnement et son adaptation aux
besoins industriels. Ces études visent à déterminer les caractéristiques électriques,
magnétiques et mécaniques de la machine, telles que la puissance, le couple, le courant
nominal, le glissement ou encore le rendement.

2.2 Les différentes constitutions des grandeurs électriques


indispensable de la machine asynchrone à cage.

2.2.1 Grandeurs d’alimentation

 Tension statorique (U₁) : tension appliquée aux enroulements du stator.

 Fréquence (f) : détermine la vitesse du champ tournant.

 Courant statorique (I₁) : courant absorbé par le stator lors du fonctionnement.

2.2.2 Grandeurs magnétiques

 Flux magnétique (Φ) : créé par le courant statorique, il induit les tensions et courants
dans le rotor.

 Inductances :

o Inductance de fuite (stator et rotor),

Inductance de magnétisation (Lₘ).

2.3 Grandeurs du rotor


 Courant rotorique (I₂’) : courant induit dans le rotor (souvent ramené au stator).

10
 Résistance rotorique ramenée (R₂’) : modélise les pertes joules dans les barres du rotor.

 Réactance de fuite rotorique (X₂’) : traduit la dispersion du flux dans le rotor.

2.4 Les grandeurs de fonctionnement et les équations


fondamentaux
 Glissement (s) :

𝑵𝒔 − 𝑵
𝒈=
𝑵𝒔
Où :

g : glissement

N : Vitesse du rotor (Tr /min)

Ns est la vitesse synchrone et N la vitesse du rotor. Alors la vitesse de sychrone

120 ⋅ 𝒇
𝑵𝒔 =
𝒑
Avec Ns en Tr/min et f est la fréquence d’alimentation en (Hz ) puis p est nobre de paire de
pôles.

 Puissance absorbée (P₁) : puissance électrique fournie au stator.

 Puissance électromagnétique (P_{em}) : puissance transmise du stator au rotor par le


champ.

 Puissance mécanique utile (P_u) : puissance disponible sur l’arbre du moteur.

𝑃𝑎𝑏𝑠 = √3. 𝑈. 𝐼. 𝑐𝑜𝑠𝜃

Où :

Pabs : Puissance électriques absorbés.

U : tension composé

I : Courant de ligne (A)

11
cos 𝜃 : Facteur de puissance.

Puissance mécanique utile :


2𝜋𝑁
P1= C⋅ω = C . 60

Avec P1 : Puissance mécanique utile (w)

C : Couple utile (Nm).

𝜔 : Vitesse angulaire (rad/s)

Rendement :

𝑷2
𝜼= × 100%
𝑷1

2.4.1 Modèle Équivalent.

Le moteur asynchrone peut être représenté par un circuit équivalent monophasé


permettant de modéliser son comportement électrique. Ce modèle comprend :
 Résistance statorique (𝑅1 ) : Représente les pertes joules au stator
 Inductance de fuite statorique (𝐿1 ) : Représente les flux de fuite au stator
 Résistance des pertes fer (𝑅𝑓𝑒 ) : Représente les pertes magnétiques

 Inductance magnétisante (𝐿𝑚 ) : Représente le circuit magnétique


 Résistance rotorique ramenée au stator (𝑅2′ ) : Représente les pertes joules au rotor
 Inductance de fuite rotorique ramenée au stator (𝐿2′ ) : Représente les flux de fuite au
rotor
1−𝑔
La charge mécanique est représentée par une résistance variable 𝑅𝑐 ′ = 𝑅2′ ⋅
𝑔
qui dépend du glissement.

2.5 Modélisation Mathématique

12
2.5.1 Hypothèses

Pour simplifier la modélisation tout en conservant une précision acceptable, les


hypothèses suivantes ont été adoptées :
Alimentation constante : La tension et la fréquence d'alimentation sont maintenues cons-
tantes à leurs valeurs nominales.

Température constante : Les variations de résistance dues à l'échauffement sont négligées.

Saturation magnétique négligée : La perméabilité du circuit magnétique est considérée


comme constante.

Pertes constantes : Les pertes fer et les pertes mécaniques sont supposées constantes quelle
que soit la charge.

Harmoniques négligées : Seul le fondamental des grandeurs électriques est considéré.

Symétrie parfaite : Le moteur est parfaitement équilibré et symétrique.

2.5.2 4.2 Fonctions d'Approximation

Sur la base de caractéristiques typiques de moteurs asynchrones et de données


constructeurs, les fonctions d'approximation suivantes ont été développées :

Rendement :
𝑃2 𝑃2 2 𝑃2 3
𝜂(𝑃2 ) = 𝜂𝑛 [0,3 + 0,7 ( ) − 0,2 ( ) + 0,1 ( ) ]
𝑃𝑛 𝑃𝑛 𝑃𝑛

Glissement :
𝑃2 𝑃2 2
𝑔(𝑃2 ) = 𝑔𝑛 [0,2 + 0,8 ( ) + 0,1 ( ) ]
𝑃𝑛 𝑃𝑛

Facteur de puissance:
𝑃2 𝑃2 2 𝑃2 3
cos𝜙(𝑃2 ) = cos𝜙𝑛 [0,3 + 0,7 ( ) − 0,1 ( ) + 0,05 ( ) ]
𝑃𝑛 𝑃𝑛 𝑃𝑛

Courant absorbé :

13
𝑃2 𝑃2 2 𝑃2 3
𝐼(𝑃2 ) = 𝐼𝑛 [0,3 + 0,7 ( ) + 0,2 ( ) − 0,1 ( ) ]
𝑃𝑛 𝑃𝑛 𝑃𝑛

Ces fonctions polynomiales ont été ajustées pour reproduire le comportement


typique d'un moteur asynchrone :
 Un rendement qui augmente rapidement avec la charge jusqu'à atteindre un maximum au-
tour de 75-100% de charge
 Un glissement qui augmente quasi-linéairement avec la charge
 Un facteur de puissance qui s'améliore avec la charge
 Un courant qui augmente de manière monotone avec la charge

2.6 Implémentation sous Mathcad

2.6.1 Données Nominales

Les calculs sont basés sur un moteur asynchrone triphasé à cage d'écureuil de
caractéristiques nominales suivantes :
// DONNÉES NOMINALES
Pn := 15 kW // Puissance mécanique nominale
Un := 400 V // Tension nominale (triphasée)
fn := 50 Hz // Fréquence nominale
ηn := 0.92 // Rendement nominal
cosφn := 0.85 // Facteur de puissance nominal
In := 28 A // Courant nominal
Nn := 1450 tr/min // Vitesse nominale

2.6.2 Calculs Préliminaires

Nombre de paires de pôles :


120⋅𝑓𝑛 120⋅50
Pour que ;𝑝= = =4
𝑁𝑠 1500

Vitesse synchrone :

Glissement nominal :

14
𝑁𝑠 − 𝑁𝑛 1500 − 1450
𝑔𝑛 = = = 0.0333
𝑁𝑠 1500

Puissance électrique absorbée nominale :


𝑃𝑛 15
𝑃1𝑛 = = = 16.304 kW
𝜂𝑛 0.92

Couple nominal :
𝑃𝑛 15000
𝐶𝑛 = = = 98.78 Nm
𝑁 1450
2𝜋 60𝑛 2𝜋
60

2.6.3 Programme Principal

Le programme Mathcad complet comprend :


Définition des données nominales
Calcul des paramètres de base (vitesse synchrone, glissement nominal, etc.)
Définition des fonctions caractéristiques (rendement, glissement, facteur de puissance,
courant, couple)
Création de la gamme de charge : P₂ = 0 à 18 kW (0 à 120% de Pn)
Génération des courbes caractéristiques
Création du tableau de résultats
Validation des résultats

// CALCULS PRÉLIMINAIRES

p := 4 // Nombre de paires de pôles (déduit de Nn)

Ns := (120 * fn) / p // Vitesse synchrone = 1500 tr/min

gn := (Ns - Nn) / Ns // Glissement nominal = 0.0333

P1n := Pn / ηn // Puissance électrique nominale = 16.304 kW

Cn := Pn / (2 * π * Nn / 60) // Couple nominal = 98.78 Nm

// DÉFINITION DES FONCTIONS CARACTÉRISTIQUES

15
η(P2) := if(P2 ≤ 0, 0, ηn * (0.3 + 0.7 * (P2/Pn) - 0.2 * (P2/Pn)^2 + 0.1 * (P2/Pn)^3))

g(P2) := gn * (0.2 + 0.8 * (P2/Pn) + 0.1 * (P2/Pn)^2)

cosφ(P2) := if(P2 ≤ 0, 0.1, cosφn * (0.3 + 0.7 * (P2/Pn) - 0.1 * (P2/Pn)^2 + 0.05 *
(P2/Pn)^3))

I(P2) := if(P2 ≤ 0, 0.4 * In, In * (0.3 + 0.7 * (P2/Pn) + 0.2 * (P2/Pn)^2 - 0.1 * (P2/Pn)^3))

P1(P2) := P2 / η(P2)

2.6.4 Tableau de Résultats

Charge (%) P₂ (kW) η (%) g (%) cosφ I (A) P₁ (kW) C (Nm)

0 0.00 0.0 0.67 0.10 11.2 0.00 0.0

25 3.75 86.5 1.17 0.62 15.4 4.34 24.7

50 7.50 90.2 1.83 0.75 19.6 8.31 49.4

75 11.25 91.8 2.67 0.81 23.8 12.25 74.1

100 15.00 92.0 3.33 0.85 28.0 16.30 98.8

120 18.00 91.5 4.13 0.87 31.9 19.67 118.5

Conclusion
En guise de conclusion, la machine asynchrone à cage d’écureuil présente un
fonctionnement fiable et robuste, avec un rendement et un facteur de puissance
optimaux près de la charge nominale. L’évolution du courant, du glissement et du
couple confirme la proportionnalité entre les grandeurs électriques et mécaniques,
ce qui valide le bon dimensionnement de la machine. Ainsi, les résultats obtenus
montrent que la MAS à cage d’écureuil constitue une solution efficace et adaptée
aux applications industrielles nécessitant à la fois performance et simplicité de
conception.

16
3 : RESULTATS ET ANALYSE

Introduction
Dans ce projet, nous avons étudié les performances d’une machine asynchrone en fonction de la
charge appliquée. Les résultats obtenus mettent en évidence l’évolution du rendement, du facteur
de puissance, du courant, du glissement, ainsi que la relation entre le couple et la puissance
absorbée. L’analyse de ces courbes permet d’évaluer le comportement global de la machine et sa
performance en régime nominal.

Rendement (η) et Facteur de Puissance (cos φ) en fonction de la charge


Rendement (η):

Il démarre à 0 % à vide (car aucune puissance utile n’est fournie). Il augmente rapidement avec
la charge jusqu’à atteindre un maximum d’environ 92 % vers 15 kW. Ensuite, il légèrement di-
minue lorsque la charge approche la pleine puissance → pertes supplémentaires (fer, cuivre, mé-
canique).

17
Dans le cas de facteur de puissance (cos φ), il est très faible à vide (~10 %), car la machine con-
somme surtout de la puissance réactive (magnétisation).

Il augmente avec la charge pour atteindre environ 87 % → meilleure utilisation de la puissance


active.

On remarque que le cos φ est toujours un peu inférieur à 1 → déphasage inévitable entre tension
et courant.

En résume, le moteur est plus performant (haut rendement et bon cos φ) lorsqu’il fonctionne
proche de sa charge nominale.

Courant absorbé (I) et Glissement (g) en fonction de la charge

courant absorbé (I):

18
Il est déjà non nul à vide (≈11 A) → courant magnétisant.
Il augmente quasi-linéairement avec la charge jusqu’à environ 32 A à pleine charge.
Cela traduit une consommation de plus en plus élevée pour fournir plus de puissance mécanique.
Et ici en parlant de glissement (g), constatons qu’il est très faible à vide (0,67 %).
Il croît progressivement avec la charge → environ 4,1 % à pleine charge.
Ce comportement est normal car le glissement est nécessaire pour que le moteur développe du
couple.

Donc, plus la charge augmente, plus le moteur consomme du courant et plus le glissement
s’élève. Cependant, le glissement reste relativement faible, ce qui confirme que la machine est
bien dimensionnée.

Couple (C) et Puissance électrique absorbée (P₁) en fonction de la charge


Couple (C) :

Croît presque linéairement avec la puissance mécanique P₂.

19
Cela confirme que le moteur fonctionne dans une zone stable et proportionnelle.
À pleine charge, le couple atteint environ 118,5 Nm.
De la puissance électrique absorbée (P₁), elle suit aussi une évolution quasi-linéaire avec P₂, mais
elle est toujours légèrement supérieure à P₂ (cause : pertes).
Exemple : à P₂ = 15 kW, P₁ ≈ 16,3 kW → rendement élevé (~92 %).

Dans ce cas, on dirait alors que ce couple et la puissance électrique suivent la charge de manière
régulière, preuve que le moteur fonctionne en régime normal sans instabilités.

Interpretation
 Le rendement est maximal proche de la charge nominale.

 Le facteur de puissance s'améliore avec la charge.

 Le courant à vide est environ 40% du courant nominal.

 Le glissement augmente linéairement avec la charge grandeurs de performance

 Couple électromagnétique (Cₑ) : lié à l’interaction champ statorique ↔ champ roto-


rique.

 Rendement (η) : rapport entre la puissance utile et la puissance absorbée.

 Facteur de puissance (cos φ) : reflète le déphasage entre tension et courant statorique.

20
CONCLUSION GENERALE
Le projet "Calcul des Caractéristiques de Service d'un Moteur Asynchrone à Cage"
représente une réussite complète tant sur le plan technique que pédagogique. Il fournit une
base solide pour la compréhension, l'analyse et l'optimisation des moteurs asynchrones,
contribuant ainsi à l'amélioration de l'efficacité énergétique dans l'industrie.

La méthodologie développée, combinant modélisation mathématique rigoureuse et


visualisation graphique avancée, constitue un outil précieux pour les ingénieurs, techniciens
et étudiants dans le domaine de l'électrotechnique et de l'optimisation énergétique.

Ce projet confirme que la modélisation précise couplée à une visualisation efficace des
données reste un atout majeur pour l'optimisation des systèmes électriques industriels.

21

Vous aimerez peut-être aussi