Guide complet — Simulation du Moteur Asynchrone Triphasé en LTspice
Ce guide explique comment compléter et utiliser le schéma LTspice du moteur asynchrone triphasé
(3ph_ACMotor). Chaque partie du circuit et chaque ligne du sous-circuit est expliquée de façon
simple.
1 Lecture du schéma — Ce qui est déjà présent
Élément Nom Rôle À compléter ?
Alimentation triphasée 400 V / 50 Hz du OUI — tension et
Sources Va, Vb, Vc Va, Vb, Vc
réseau fréquence à définir
3ph_ACMoto Symbole du moteur asynchrone triphasé
Moteur U1 Déjà présent
r avec ses 6 bornes
Résistance de mesure sur la borne W
Résistance R1 (1kΩ) R1 Déjà présent
(vitesse)
Interrupteur sur la borne W — peut être
Switch S S À adapter
remplacé par un fil
Flèche (GND) 0 Référence de masse (neutre du réseau) Déjà présent
2 Les 6 bornes du symbole 3ph_ACMotor
Le sous-circuit .subckt 3ph_ACMotor a b c N W J possède 6 bornes. Voici leur rôle exact :
Borne Type Rôle Connexion dans LTspice
Phase A du réseau triphasé
a Entrée → Va (source triphasée phase 0°)
(stator)
Phase B du réseau triphasé
b Entrée → Vb (source triphasée −120°)
(stator)
Phase C du réseau triphasé
c Entrée → Vc (source triphasée −240°)
(stator)
Neutre du réseau (référence 0
N Entrée → GND (masse 0)
V)
Vitesse angulaire Ω en rad/s (1 → R1 (1kΩ) vers GND V(W)×60/(2π) = N
W Sortie
V = 1 rad/s) tr/min
Couple moteur C en Nm (1 V = 1 → Résistance de charge (couple
J Sortie
Nm) résistant) ou masse via résistance 1 MΩ
Comment lire V(W) en tr/min dans LTspice ?
Après simulation, dans le Waveform Viewer, clic droit sur l'axe Y → "Add Plot Plane"
Puis saisir l'expression : V(W)*60/(2*3.14159) => donne la vitesse en tr/min
À vide (Cr=0), le moteur 1.1kW/2 paires de pôles atteint environ 1470-1490 tr/min
La vitesse de synchronisme est Ns = 60×50/2 = 1500 tr/min (jamais atteinte)
3 Complétion 1 — Sources triphasées Va, Vb, Vc
Les trois sources de tension doivent délivrer un système triphasé équilibré 400 V / 50 Hz. Dans
LTspice, on utilise des sources SINE avec un décalage de phase de 120° entre chaque phase.
Formule de la tension de phase :
Réseau 400 V = tension composée (entre deux phases). La tension de phase (entre phase et neutre)
vaut :
V_phase = 400 / √3 ≈ 230.9 V (valeur efficace)
V_crête = 230.9 × √2 ≈ 326.6 V (valeur de crête pour la source SINE)
Attributs des sources dans LTspice :
Sources triphasées — attribut Value à copier
* Paramètres à définir (en haut du schéma via .param)
.param fn=50 ; fréquence réseau 50 Hz
.param Vn=400 ; tension composée réseau 400 V
.param Vph={Vn/sqrt(3)} ; tension de phase = 230.9 V
* Source Va — Phase A (référence 0°)
Va a 0 SINE(0 {Vph*sqrt(2)} {fn} 0 0 0)
* Source Vb — Phase B (retard de 120°)
Vb b 0 SINE(0 {Vph*sqrt(2)} {fn} 0 0 -120)
* Source Vc — Phase C (retard de 240°)
Vc c 0 SINE(0 {Vph*sqrt(2)} {fn} 0 0 -240)
Syntaxe SINE : SINE(offset amplitude fréquence délai amortissement phase)
• Offset = 0 : pas de composante continue
• Amplitude = Vph×√2 : valeur de crête (LTspice attend la crête, pas l'efficace)
• Fréquence = fn = 50 : fréquence réseau
• Phase = 0 / -120 / -240 : décalages temporels entre phases (en degrés)
Retenir :
SINE(0 326.6 50) pour Va. SINE(0 326.6 50 0 0 -120) pour Vb. SINE(0 326.6 50 0 0 -240) pour Vc. Ou utiliser .param
pour les valeurs.
4 Complétion 2 — Borne J : couple résistant
La borne J du moteur représente le couple électromagnétique développé (en Nm, sous forme de
tension en V). Pour simuler une charge, on doit appliquer un couple résistant sur cette borne. Trois
méthodes selon le cas :
Méthode A — Démarrage à vide (aucune charge)
Connecter J à la masse via une grande résistance (1 MΩ). Le moteur démarre et accélère sans
charge → atteint la vitesse de synchronisme moins le glissement à vide (très faible).
Méthode A — À vide
* Borne J connectée à GND via 1 MΩ = charge nulle (à vide)
R_J J 0 1Meg
Méthode B — Couple résistant constant (levage, convoyeur)
On injecte un courant constant dans le noeud J. Ce courant représente le couple résistant en Nm.
Une source de courant DC suffit.
Méthode B — Couple constant
* Couple résistant constant = 7.3 Nm (moteur 1.1kW à 1440 tr/min)
* Calcul : Cr = Pn / Omega_n = 1100 / (1440*2*pi/60) = 7.3 Nm
I_Cr J 0 7.3
* Ou en paramètre :
.param Cr=7.3
I_Cr J 0 {Cr}
Méthode C — Couple résistant proportionnel à la vitesse² (pompe/ventilateur)
Le couple d'une pompe centrifuge est proportionnel au carré de la vitesse : Cr = k × Ω². On utilise une
source B non-linéaire.
Méthode C — Couple quadratique (pompe/ventilateur)
* Couple résistant quadratique : Cr = k * Omega^2
* k = Cr_nominal / Omega_nominal^2
* Omega_nominal = 1440 * 2*pi/60 = 150.8 rad/s
* k = 7.3 / 150.8^2 = 3.21e-4
.param k_pompe=3.21e-4
B_Cr J 0 I={k_pompe*V(W)*V(W)}
* Note : V(W) = Omega en rad/s
Retenir :
À vide : R_J 1Meg. Charge constante : I_Cr = Cr en Nm. Pompe : B_Cr I={k*V(W)^2}. Le noeud J ne doit JAMAIS rester
flottant.
5 Commandes de simulation et mesures
.tran — Simulation transitoire
On utilise une simulation transitoire pour observer le démarrage du moteur. Il faut simuler au moins
300 à 500 ms pour que le moteur atteigne son régime.
Directive .tran à coller dans le schéma
* Simulation transitoire de 0 à 500 ms, pas maxi 100 µs
.tran 0 500m 0 100u
* Options importantes pour la convergence (moteur = circuit difficile)
.options method=gear reltol=0.01 abstol=1u vntol=1m
* Si problème de convergence : augmenter le pas maxi
.tran 0 500m 0 500u
Ce que signifient les options :
• method=gear : méthode d'intégration numérique plus stable pour les circuits avec inductances
(comme un moteur). Évite les oscillations parasites.
• reltol=0.01 : tolérance relative de 1% (au lieu de 0.1% par défaut). Accélère la simulation sans
perdre en précision visible.
• abstol=1u : tolérance absolue sur les courants. Évite les faux problèmes de convergence autour
de zéro.
Mesures à observer dans le Waveform Viewer
Expression LTspice Ce que ça mesure Valeur attendue (régime)
V(W) Vitesse angulaire Ω en rad/s ≈ 150-153 rad/s
V(W)*60/(2*3.14159) Vitesse N en tr/min ≈ 1430–1470 tr/min
V(J) Couple électromagnétique en Nm ≈ 7.3 Nm (charge nominale)
I(Va) Courant absorbé phase A en A ≈ 2.6 A nominal
Puissance instantanée phase A en
V(a)*I(Va) ≈ 400 W (1/3 de Pa)
W
(1-V(W)/(2*3.14159*50/2)
Glissement g (sans unité) ≈ 0.03 à 0.05 (3-5%)
)
6 Netlist complète prête à utiliser
Voici la netlist complète à copier dans un fichier .net ou à recréer dans LTspice. Elle intègre toutes
les compléments expliqués ci-dessus.
Netlist principale — moteur_3ph_simulation.net
* Simulation moteur asynchrone 1.1kW / 400V / 50Hz
.param fn=50 Vn=400 Vph={Vn/sqrt(3)}
* ■ Sources triphasées ■
Va a 0 SINE(0 {Vph*sqrt(2)} {fn} 0 0 0)
Vb b 0 SINE(0 {Vph*sqrt(2)} {fn} 0 0 -120)
Vc c 0 SINE(0 {Vph*sqrt(2)} {fn} 0 0 -240)
* ■ Moteur ■
X1 a b c 0 W J 3ph_ACMotor Zp=2 J=2.1m
* ■ Mesure vitesse (borne W) ■
R1 W 0 1k
* ■ Charge sur borne J (choisir UNE méthode) ■
R_J J 0 1Meg ; méthode A : à vide
* I_Cr J 0 7.3 ; méthode B : couple constant 7.3 Nm
* B_Cr J 0 I={3.21e-4*V(W)*V(W)} ; méthode C : pompe
* ■ Simulation ■
.tran 0 500m 0 100u
.options method=gear reltol=0.01 abstol=1u vntol=1m
Paramètres internes du sous-circuit (déjà dans le .subckt)
Paramètre Valeur Signification physique
Résistance
Rs = 7.5 Ω Résistance des bobinages en cuivre du stator
statorique
Résistance
Rr = 5.5 Ω Résistance équivalente de la cage d'aluminium
rotorique
Inductance de fuite
Ls = 26.43 mH Flux qui ne traverse pas l'entrefer (fuite)
stator
Inductance de fuite
Lr = 18.8 mH Flux de fuite côté rotor
rotor
Inductance
Lm = 351.05 mH Flux principal traversant l'entrefer = couplage
magnétisante
Rf = 2 kΩ Résistance fer Pertes fer (courants de Foucault + hystérésis)
Nb de paires de
Zp = 2 2 paires de pôles → Ns = 1500 tr/min à 50 Hz
pôles
J = 2.1 mkg·m² Moment d'inertie Inertie mécanique du rotor (résistance à l'accélération)
7 Fonctionnement interne du sous-circuit — Modèle dq
Le sous-circuit implémente le modèle d-q du moteur asynchrone. C'est la représentation
mathématique la plus utilisée pour simuler les moteurs dans un repère fixe (stator). Voici les grandes
étapes :
Bloc Éléments LTspice Ce que ça fait
Convertit les 3 tensions abc en courants vers les axes
1 — Filtrage Fa, Fbd, Fbq, Fcd, Fcq
d et q. Les résistances 1MΩ évitent les noeuds
entrée Ra, Rb, Rc (1MΩ)
flottants.
Circuit équivalent de la composante "d" : stator
2 — Axe d Gab, Gac, Rds, Lds Ldm
(Ls_+Rs_), magnétisation (Lm_■Rf_), rotor (Lr_+Rr_).
(direct) (Rf_), Ldr, Ad
Ad = OTA qui couple avec la vitesse Ω.
3 — Axe q Gq, Rqs, Lqs Lqm (Rf_), Identique à l'axe d mais décalé de 90°. As = OTA qui
(quadrature) Lqr, As couple avec Ω dans l'autre sens.
4 — Couple Hdm, Hqm, Hdr, Hqr Calcule le couple C = (3/2)×Zp×Lm×(Id_m×Iq_r −
électromag. A_dm_qr, A_qm_dr Iq_m×Id_r). Les OTA font le produit croisé des flux.
5 — Vitesse Intègre le couple net (C − Cr) divisé par J pour obtenir
Gw, Cj (=J)
mécanique Ω. Cj = condensateur d'inertie : C×dΩ/dt = couple net.
Pourquoi le modèle d-q ?
Le courant triphasé abc est difficile à analyser (3 variables sinusoïdales).
La transformation de Clarke (abc → αβ) puis de Park (αβ → dq) simplifie tout.
Dans le repère dq fixe au stator, les équations deviennent du courant continu.
LTspice simule alors un circuit électrique classique (R, L, C, sources contrôlées).
Les OTA (Operational Transconductance Amplifier) réalisent les multiplications.
Équation mécanique du moteur
La dynamique mécanique est modélisée par le condensateur Cj :
J × dΩ/dt = C_em − C_résistant
C_em = (3/2) × Zp × Lm × (I_dm × I_qr − I_qm × I_dr)
• J = moment d'inertie (2.1 mkg·m²) — comme une masse en translation
• C_em = couple électromagnétique produit par l'interaction des flux d et q
• C_résistant = couple de la charge mécanique (ce qu'on connecte à J)
Retenir :
V(W) = Omega rad/s = vitesse. V(J) = couple em. La borne J collecte le courant de charge. Cj intègre pour donner Om
ega.
Synthèse
Récapitulatif — Tout en un coup d'œil
Élément à ajouter Valeur / Expression LTspice Pourquoi
Va (source phase A) SINE(0 {Vph*sqrt(2)} {fn} 0 0 0) Phase A du réseau 400V/50Hz
SINE(0 {Vph*sqrt(2)} {fn} 0 0
Vb (source phase B) Phase B décalée de −120°
-120)
SINE(0 {Vph*sqrt(2)} {fn} 0 0
Vc (source phase C) Phase C décalée de −240°
-240)
.param fn=50 Vn=400
Paramètres globaux Facilite les calculs automatiques
Vph={Vn/sqrt(3)}
R_J J 0 1Meg (vide) ou I_Cr J 0 Couple résistant de la charge
R_J sur borne J
{Cr} mécanique
Simulation 500 ms pour voir le
.tran .tran 0 500m 0 100u
démarrage
.options method=gear reltol=0.01 Stabilité numérique du moteur
Conversion rad/s → tr/min dans le
Mesure vitesse V(W)*60/(2*3.14159)
viewer
Mesure couple V(J) Couple électromagnétique en Nm
Les deux fichiers joints contiennent la netlist complète (moteur_3ph_simulation.net) prête à être importée dans
LTspice.