Simulator für die Clarke-Park-Transformation: abc nach dq0

Transformieren Sie abc-Ströme in αβ0- und dq0-Komponenten – mit animierten Zeigern, rotierenden Bezugssystemen, amplituden- oder leistungsinvarianter Normierung und Rücktransformation.

#abc to dq0 #Clarke transform #FOC #Park transform #power invariance #three-phase systems

Der Simulator für die Clarke–Park-Transformation überführt drei Phasenströme in ortsfeste αβ0-Komponenten und rotierende dq0-Komponenten. Er richtet sich an Ingenieure, die eine Konvention der Antriebsregelung überprüfen, Studierende, die Bezugssysteme kennenlernen, und Lehrende, die erklären möchten, warum sinusförmige Motorströme zu zeitlich konstanten Reglersignalen werden können. Verwenden Sie ihn als Rechner für die Clarke-Park-Transformation oder für die Umrechnung von abc nach dq0; die Zwischenschritte bleiben dabei sichtbar.

Wählen Sie Dreiphasiger Signalverlauf, um ein cosinusförmiges Mitsystem zu erzeugen. Geben Sie den Phasenstromscheitelwert, die elektrische Frequenz und die Anfangsphase des Signals ein. Voreingestellt ist ein symmetrisches System mit 10 A Scheitelwert und 50 Hz. Mit dem Regler für den elektrischen Phasenfortschritt wählen Sie einen Zeitpunkt innerhalb einer Periode; drücken Sie ▶, um die Animation zu starten. Die Wiedergabe ist bewusst verlangsamt; die physikalische Periodendauer und der gewählte Zeitpunkt sind weiterhin in den Rechenschritten angegeben.

Der Winkel der d-Achse wird gegen den Uhrzeigersinn von der ortsfesten α-Achse aus gemessen, die mit der Achse von Phase a zusammenfällt. Die q-Achse eilt der d-Achse um 90 elektrische Grad voraus. In der Voreinstellung dreht sich das Bezugssystem mit der Signalfrequenz und einem Winkelversatz von 30° zum Stromraumzeiger; daraus ergeben sich d ≈ 8,6603 A und q = 5 A. Ändern Sie den Anfangswinkel der d-Achse, um den Raumzeiger vollständig auf der d- oder q-Achse auszurichten. Bei einem Frequenzverhältnis von null bleibt das Bezugssystem ortsfest; bei −1 folgt es dem Gegensystem.

Wählen Sie abc-Momentanwerte eingeben, um drei gleichzeitig erfasste Stromwerte einzugeben. In diesem Modus legt der Regler für den Phasenfortschritt direkt den Systemwinkel θ fest; ▶ dreht das Bezugssystem um den festgehaltenen Stromraumzeiger. Die Kurven zeigen einen Winkeldurchlauf bei konstant gehaltenen Momentanwerten und keine Vorhersage ihres zeitlichen Verlaufs. Alle Stromeingaben erfolgen in Ampere; die elektrischen Einheiten sind im metrischen und im US-Anzeigemodus identisch.

Vergleichen Sie amplitudeninvariante und leistungsinvariante Normierung. Bei Amplitudeninvarianz entspricht der αβ-Betrag für ein symmetrisches System dem Phasenstromscheitelwert. Bei Leistungsinvarianz sind die Koordinaten in diesem Fall um den Faktor √(3/2) größer; das vollständige Skalarprodukt in dq0 stimmt dann ohne zusätzlichen Faktor mit dem in abc überein. Das Werkzeug rekonstruiert abc außerdem durch inverse Park- und Clarke-Transformation.

Erhöhen Sie das Amplitudenverhältnis von Gegen- zu Mitsystem, um Unsymmetrien zu untersuchen. Der αβ-Kreis wird zu einer Ellipse, und im synchron rotierenden dq-Bezugssystem tritt eine Welligkeit mit doppelter Frequenz auf. Fügen Sie allen Phasen denselben Gleichstromanteil hinzu, um das Nullsystem zu untersuchen; wählen Sie dazu zuerst Neutralleiter / unabhängige Signale. In einer Dreileiterschaltung muss die Summe der Ströme null sein. Das Nullsystem wird separat angegeben, da es in einer zweidimensionalen Raumzeigerdarstellung nicht sichtbar ist.

Dieser Simulator führt ideale algebraische Koordinatentransformationen durch. Er modelliert weder Motordrehmoment, magnetische Sättigung, Sensorverzögerungen, Stromregelkreise, PWM-Schaltvorgänge, Wechselrichtergrenzen noch einen Lagebeobachter. Der elektrische Winkel muss bereits bekannt sein; bei einem Motor entspricht er dem mechanischen Winkel multipliziert mit der Polpaarzahl, zuzüglich des gewählten Referenzversatzes. Gleichen Sie Achsenorientierung, Phasenfolge und Normierung ab, bevor Sie Ergebnisse mit Firmware oder einem anderen Lehrbuch vergleichen.

So wird gerechnet

Die amplitudeninvariante Clarke-Transformation lautet

iα=23(ia−ib+ic2),iβ=33(ib−ic),i0=ia+ib+ic3.i_alpha=frac23left(i_a-frac{i_b+i_c}{2}right),quad i_beta=frac{sqrt3}{3}(i_b-i_c),quad i_0=frac{i_a+i_b+i_c}{3}.

Die leistungsinvariante Form verwendet k=2/3k=sqrt{2/3} für die αβ-Zeilen und k0=1/3k_0=1/sqrt3 für die Summenzeile. Die beiden vollständigen Matrizen entsprechen den Gleichungen zur Clarke-Transformation von MathWorks (https://www.mathworks.com/help/simscape-electrical/ref/clarketransform.html). TI SPRA524, pp. 4–5, verwendet leistungsinvariante ortsfeste Koordinaten mit den Bezeichnungen d und q; diese entsprechen hier α und β.

Die Park-Transformation ist eine Drehung: α liegt entlang der Achse von Phase a, d liegt um θ gegen den Uhrzeigersinn gegenüber α verdreht, und q eilt d um 90° voraus:

id=iαcos⁡θ+iβsin⁡θ,iq=−iαsin⁡θ+iβcos⁡θ.i_d=i_alphacostheta+i_betasintheta,qquad i_q=-i_alphasintheta+i_betacostheta.

Diese Vorzeichenkonvention und die zugehörige Rücktransformation stimmen mit Microchip AN1078 Rev. B, pp. 8–10, überein. Die inverse Park-Transformation dreht um +θ; die inverse Clarke-Transformation stellt den gemeinsamen Anteil i0/(3k0)i_0/(3k_0) wieder her und dividiert die αβ-Projektionen durch 3k/23k/2.

Für amplitudeninvariante Größen gilt

∑iabc2=32(id2+iq2)+3i02.sum i_{abc}^2=frac32(i_d^2+i_q^2)+3i_0^2.

Für leistungsinvariante Größen sind beide Gewichte 1. Werden Spannung und Strom nach derselben Konvention transformiert, gelten dieselben Gewichte für die Momentanleistung: p=32(vdid+vqiq)+3v0i0p=frac32(v_di_d+v_qi_q)+3v_0i_0 bei amplitudeninvarianter Normierung oder p=vdid+vqiq+v0i0p=v_di_d+v_qi_q+v_0i_0 bei leistungsinvarianter Normierung. Die angezeigte Norm hat die Einheit A², nicht Watt.

Der Signalgenerator überlagert ein Mitsystem Icos⁡(u+φ)Icos(u+varphi), ein Gegensystem mit umgekehrter b–c-Folge und dem Scheitelwert rIrI sowie einen gleichen Gleichstromanteil IcI_c in allen Phasen. Er berechnet die Summe der Phasenströme analytisch als 3Ic3I_c, um Auslöschung bei Gleitkommarechnungen zu vermeiden. θ = Frequenzverhältnis des Bezugssystems × u + Anfangswinkel der d-Achse. Die Kurven über die gesamte Periode werden unabhängig vom gerade animierten Eingabewert berechnet.

Die Referenzfälle bilden den in JET veröffentlichten numerischen Vektor, die Stromwerte aus Abbildung 4 der UPC-Arbeit und die Koordinaten der TI-Grundraumzeiger nach. In Gleichung (7) der UPC-Arbeit steht Sinus, während Gleichung (1) und Abbildung 4 Cosinus verwenden; dieses Werkzeug folgt Gleichung (1) und Abbildung 4. In Gleichung (1.13) des JET-Beitrags fehlt in einer Zeile ein Minuszeichen bei Phase c; die Gleichungen (1.14)–(1.15) verwenden den korrekten negativen Wert. Die Transformationskoeffizienten wurden nicht an diese Druckfehler angepasst.

Annahmen

Ideale, gleichzeitig erfasste, reelle Strommomentanwerte; keine Filterung, Abtastverzögerung oder Quantisierung. Das Mitsystem des Signalverlaufs hat die Phasenfolge a → b → c mit um 120° versetzten Cosinusbezügen. Das Gegensystem hat dieselbe Anfangsphase von Phase a und die umgekehrte b–c-Folge; sein Amplitudenverhältnis beträgt 0–1. Die vom Generator hinzugefügte Nullkomponente ist ein Gleichstromanteil und keine Oberschwingung mit durch drei teilbarer Ordnungszahl.

Alle Winkel sind elektrische Winkel in Grad. Das Bezugssystem wird vom Benutzer festgelegt und muss nicht mit dem Rotor eines Motors übereinstimmen. Die 2D-Darstellung zeigt nur den αβ- oder dq-Anteil; 0 ist eine separate skalare Größe. Für die Längen im Diagramm gilt während der Animation ein fester Maßstab für die gesamte Periode; der helle Kreis dient der Orientierung, die graue Ortskurve zeigt den tatsächlichen Verlauf.

Rücktransformierte Ströme und Transformationsergebnisse behalten die volle numerische Genauigkeit. Nur die Rechenschritte zur Fehleranzeige werden auf 10⁻¹⁰ A bzw. A² gerundet, damit unbedeutende Auslöschungseffekte nicht als Abweichung in der Größenordnung einer Einheit erscheinen. Für aus Quelldiagrammen abgelesene Werte gelten ausdrücklich größere Toleranzen als für gedruckte numerische Referenzvektoren.

Anwendungsbereich

Die vollständige Transformation mit drei Komponenten gilt für symmetrische und unsymmetrische reelle abc-Ströme, auch bei einer von null abweichenden Summe. Eine isolierte Dreileiterschaltung erfordert die Summe null; die Eingabeprüfung lässt eine numerische Toleranz von 10⁻⁹ A zu. Verwenden Sie bei einer Summe ungleich null einen Neutralleiter als Stromrückweg oder betrachten Sie die Eingaben als unabhängige Signale.

Die wählbaren ganzzahligen Frequenzverhältnisse des Bezugssystems −2…2 ermöglichen einen stetigen Übergang in der Animationsschleife über eine Periode. Die erzeugte Signalfrequenz beträgt 0,01–100000 Hz; die Transformation ist algebraisch, und dieser Bereich lässt keinen Rückschluss auf die Bandbreite von Reglern oder Sensoren zu. Die eingegebenen Phasenströme und der gemeinsame Stromanteil sind auf ±100000 A begrenzt; der Scheitelwert des Mitsystems liegt im Bereich 0–100000 A. Dies sind Eingabegrenzen, keine Bemessungswerte von Geräten. Im Momentanwertmodus durchlaufen die Kurven θ bei konstant gehaltenen abc-Werten. Die Rechenschritte für Periodendauer und physikalische Zeit gelten nur im Signalverlaufsmodus.

Zeitlich konstante dq-Werte setzen sinusförmige Ströme eines symmetrischen Systems und ein mit diesem System synchron rotierendes Bezugssystem voraus. Ein Gegensystem erzeugt in einem mit dem Mitsystem synchron rotierenden Bezugssystem eine Welligkeit mit doppelter Frequenz. Motordynamik, PWM, FOC-Regler und Drehmomentberechnung sind nicht enthalten.

Quellen

  • Microchip, AN1078 Rev. B (2010), Sensorless Field Oriented Control (FOC) of a Permanent Magnet Synchronous Motor (PMSM), pp. 8–10. https://ww1.microchip.com/downloads/aemDocuments/documents/OTH/ApplicationNotes/ApplicationNotes/01078B.pdf ; landing page: https://www.microchip.com/en-us/application-notes/an1078
  • Z. Yu, Texas Instruments SPRA524 (March 1999), Space-Vector PWM With TMS320C24x/F24x Using Hardware and Software Determined Switching Patterns, pp. 4–5, Eq. 3 and Figure 3. https://www.ti.com/lit/pdf/spra524
  • MathWorks, Clarke Transform, complete amplitude- and power-invariant matrices. https://www.mathworks.com/help/simscape-electrical/ref/clarketransform.html
  • F. Filipović et al., Systematic Derivation of Clarke and Park Transformations through Vector Representation in Three-Phase Systems, Journal of Energy Technology 18(3), 145–160 (2025), pp. 156–157. https://doi.org/10.18690/jet.18.3.145-160.2025
  • I. Deiros Gras, Geometrically-Based Dynamical Modeling of Permanent Magnet Synchronous Motors, UPC bachelor's thesis (25 June 2024), printed pp. 14–18, Figure 4. https://upcommons.upc.edu/server/api/core/bitstreams/1124c6e7-d4a4-4070-83ae-c829c3764d64/content

Eingaben

EingabeSymbolEinheitVorgabe
Signalquelle
Signalverlauf: Phasenfortschritt animieren. Momentanwerte: abc konstant halten und das dq-Bezugssystem drehen.
signal_modeDreiphasiger Signalverlauf
Clarke-Normierung
Amplitude: Bei Symmetrie entspricht der αβ-Betrag dem Phasenstromscheitelwert. Leistung: Die vollständige αβ0-Transformation erhält das Skalarprodukt in abc.
normalizationAmplitudeninvariant (2/3)
Stromrückweg
Wählen Sie Neutralleiter / unabhängige Signale, um das Nullsystem zu untersuchen. In einem isolierten Dreileitersystem kann die Summe der Phasenströme nicht von null abweichen.
return_pathDreileiter (a + b + c = 0)
Elektr. Phasenfortschritt / Systemwinkel
Drücken Sie ▶. Im Signalverlaufsmodus werden Signale und Bezugssystem gemeinsam weitergeführt; im Momentanwertmodus wird dieser Winkel als θ verwendet. Die Wiedergabe ist für die Lehre verlangsamt.
angle_edeg40
Phasenstromscheitelwert des Mitsystems
Scheitelwert, kein Effektivwert. Geben Sie für einen symmetrischen sinusförmigen Verlauf den √2-fachen Phasenstrom-Effektivwert ein.
I_peakA10
Elektrische Signalfrequenz
Legt die physikalische Periodendauer und die verstrichene Zeit fest; die Wiedergabegeschwindigkeit mit ▶ dient der Veranschaulichung.
f_eHz50
Anfangsphase des Signals
Cosinusbezug: Phase a = Ipeak cos(progression + phase), vor Überlagerung von Gegen- und Nullsystem.
phase_ideg20
Amplitudenverhältnis Gegen- / Mitsystem
Bei 0 liegt ein symmetrisches Mitsystem vor. Fügt ein Gegensystem mit umgekehrter b–c-Folge und derselben Anfangsphase von Phase a hinzu; bei 1 ist die αβ-Ortskurve geradlinig.
negative_ratio0
Gleicher Gleichstromanteil je Phase
Wird a, b und c gleichermaßen überlagert. Werte ungleich null erfordern einen Neutralleiter oder einen anderen Stromrückweg.
I_commonA0
Frequenzverhältnis Bezugssystem / Signal
1 folgt dem Mitsystem; 0 hält das Bezugssystem ortsfest; −1 folgt dem Gegensystem. Ganzzahlige Verhältnisse ermöglichen eine nahtlose Animationsschleife.
frame_ratio1
Anfangswinkel der d-Achse
θ = Frequenzverhältnis × Phasenfortschritt + dieser Anfangswinkel. Alle Winkel sind elektrisch.
theta_offsetdeg-10
Momentaner Phasenstrom aa_entryA10
Momentaner Phasenstrom bb_entryA-5
Momentaner Phasenstrom cc_entryA-5

Ergebnisse

Mit den obigen Vorgabewerten liefert der Rechner:

ErgebnisWert
d-Achsenstrom8,66025 A
q-Achsenstrom5 A
Nullsystem 00 A
Ortsfeste Stromkomponente α5 A
Ortsfeste Stromkomponente β8,66025 A
Raumzeigerbetrag10 A
Elektrischer d-Achsenwinkel30 deg
abc-Rekonstruktionsfehler0 A

Mit Quellen abgeglichen

Jeder Fall unten wurde mit diesem Rechner berechnet und stimmt mit dem Wert aus seiner Quelle überein (innerhalb der Toleranz, wo eine angegeben ist).

JET 2025, gedrucktes Zahlenbeispiel zur Clarke/Park-Transformation

Quelle: Filipović et al. (2025), DOI 10.18690/jet.18.3.145-160.2025, printed pp. 156–157, Eqs. (1.13)–(1.15). Published voltage vector [141.42, −70.71, −70.71] V gives [α,β,0]=[141.42,0,0] V and [d,q,0]=[141.42,0,0] V at θ=0. Applied here with a 1 A/V scaling: these linear coordinate transforms do not depend on the physical unit. Phase-c sign is omitted in one line of Eq. (1.13), but cos(+120°) is negative and the later matrix input correctly uses −70.71. The correct sign is retained.

Eingaben: Signalquelle = abc-Momentanwerte eingeben; Elektr. Phasenfortschritt / Systemwinkel = 0 deg; Momentaner Phasenstrom a = 141,4 A; Momentaner Phasenstrom b = -70,71 A; Momentaner Phasenstrom c = -70,71 A; die übrigen wie oben.

ErgebnisQuelleDieser Rechner
Ortsfeste Stromkomponente α141,4 A141,4 A
Ortsfeste Stromkomponente β0 A0 A
Nullsystem 00 A0 A
d-Achsenstrom141,4 A141,4 A
q-Achsenstrom0 A0 A

UPC-Abschlussarbeit, Abbildung 4: 20 A Effektivwert, 50 Hz, 30°

Quelle: I. Deiros Gras, UPC thesis (2024), printed p. 14 Eq. (1), p. 18 Figure 4 (attached PDF pp. 16 and 20). The lower graph publishes constant levels d≈24.5 A and q≈14.1 A for a 20 A RMS balanced cosine set at 50 Hz and initial phase30°, with θ=100πt. Values are read from the published plot, not generated by this calculator; 0.3 A tolerance allows plot-reading resolution. Eq. (7), printed p.17, changes cos to sin inconsistently with Eq. (1) and Figure 4; follow Eq. (1) and the plotted positive d,q levels.

Eingaben: Phasenstromscheitelwert des Mitsystems = 28,284 A; Anfangsphase des Signals = 30 deg; Anfangswinkel der d-Achse = 0 deg; Elektr. Phasenfortschritt / Systemwinkel = 90 deg; die übrigen wie oben.

ErgebnisQuelleDieser Rechner
d-Achsenstrom24,5 ± 0,3 A24,495 A
q-Achsenstrom14,1 ± 0,3 A14,142 A
Nullsystem 00 A0 A

TI SPRA524 Grundraumzeiger U0 (100)

Quelle: TI SPRA524 (1999), Table 1 p.4 and Figure 3 p.5, https://www.ti.com/lit/pdf/spra524 . Switching state100 has phase vector [2/3, −1/3, −1/3] V/Vdc and published stationary coordinates (2/√6,0). The note calls its fixed axes d,q; these are α,β here. Test uses an equivalent 1 A per normalized unit and power-invariant normalization, with θ=0. Expected coordinates are the source's closed-form Figure 3 coordinates, evaluated as constants independently of the calculator.

Eingaben: Signalquelle = abc-Momentanwerte eingeben; Clarke-Normierung = Leistungsinvariant (√(2/3)); Elektr. Phasenfortschritt / Systemwinkel = 0 deg; Momentaner Phasenstrom a = 0,6667 A; Momentaner Phasenstrom b = -0,3333 A; Momentaner Phasenstrom c = -0,3333 A; die übrigen wie oben.

ErgebnisQuelleDieser Rechner
Ortsfeste Stromkomponente α0,8165 A0,8165 A
Ortsfeste Stromkomponente β0 A0 A
d-Achsenstrom0,8165 A0,8165 A
q-Achsenstrom0 A0 A

TI SPRA524 Grundraumzeiger U60 (110)

Quelle: TI SPRA524 (1999), Table 1 p.4 and Figure 3 p.5, https://www.ti.com/lit/pdf/spra524 . Switching state110 gives phase vector [1/3,1/3,−2/3] and the published fixed-axis space vector (1/√6,1/√2). Test scales normalized voltage values to equal numerical amperes and uses θ=0, so αβ=dq. Expected numbers are evaluations of the published closed-form Figure 3 coordinates, not results copied from this calculator.

Eingaben: Signalquelle = abc-Momentanwerte eingeben; Clarke-Normierung = Leistungsinvariant (√(2/3)); Elektr. Phasenfortschritt / Systemwinkel = 0 deg; Momentaner Phasenstrom a = 0,3333 A; Momentaner Phasenstrom b = 0,3333 A; Momentaner Phasenstrom c = -0,6667 A; die übrigen wie oben.

ErgebnisQuelleDieser Rechner
Ortsfeste Stromkomponente α0,4082 A0,4082 A
Ortsfeste Stromkomponente β0,7071 A0,7071 A
d-Achsenstrom0,4082 A0,4082 A
q-Achsenstrom0,7071 A0,7071 A

Related Resources

Responses

Deine E-Mail-Adresse wird nicht veröffentlicht. Erforderliche Felder sind mit * markiert

L ading...