Selbstfahrendes Motorad mit Arduino

Autoren: Junjie Lyu
Betreuer: Prof. Schneider
Art: PA
Projektlaufzeit: 04/2021-04/2022

Abb. 1: Arduino Engineering Kit - Selbstfahrendes Motorad

Thema

Aufbau eines selbsfahrenden Motorades, welches sich selbst während der Fahrt ausbalanciert.

Ziel

Das Arduino Engineering Kit ermöglicht den Aufbau dreier regelungstechnischer Herausforderungen. In diesem Projekt soll ein selbstfahrendes Motorad gebaut und programmiert werden.

Aufgabenstellung

1. Einarbeitung in das Thema, auch aus regelungstechnischer Sicht
2. Identifikation des Regelstrecke
3. Sichtung und Test des bestehenden Bausatzes
4. Aufbau des Systems (ggf. Platinenfertigung, etc.)
5. Vergleichen und bewerten Sie verschiedene Regleransätze (P, PI, PID und andere).
6. Modellbasierte Programmierung der Hardware via Matlab und Simulink
7. Test des Segway
8. Dokumentation nach wissenschaftlichem Stand
9. Erstellung von Gefährdungsbeurteilung und Betriebsanweisung

Anforderung

• Wissenschaftliche Vorgehensweise (Projektplan, etc.)
• Wöchentliche Fortschrittsberichte (informativ)
• Projektvorstellung im Wiki
• Machen Sie ein tolles Videos, welches die Funktion visualisiert.

Gantt - Chart

Bild1: Gantt Chart für Projektarbeit Selbstfahrendes Motorrad. Junjie Lyu (Diskussion) 12:54, 17. Okt. 2021

Datei:Gantt Chart für Projektarbeit.xlsx

Die Projektarbeit einschließlich Projektseminar wird im Sommersemester 2021 angemeldet. Das heißt, die Benotung der Projektarbeit muss bis zum Ende des Wintersemesters 21/22 erfolgt sein.

Bewertung des Bausatzes

Das Motorrad ist ein zweirädrig Roboter, der mit Hilfe einer rotierenden Scheibe (Trägheitsrad) das Gleichgewicht halten und sich bewegen kann, um zu kompensieren, wenn das Motorrad das Gleichgewicht verliert. Das Motorrad wird von einem Arduino MKR1000, dem Arduino MKR Motor Carrier, einem Gleichstrommotor zum Bewegen des Hinterrads, einem Encoder, einem Gleichstrommotor zum Steuern des Trägheitsrads, einer 6-Achsen-IMU, einem Standardservomotor zum Lenken des Motorradgriffs, einem Abstandssensor (Ultraschallsensor) und einem Drehzahlmesser (Hallsensor) gesteuert. Der Hardwareaufbau konnte mithilfe der Anweisung des Anleitungsvideos zusammengebaut werden. Das Kabel vom Gleichstrommotor zum Steuern des Trägheitsrads ist nicht lang genug. Versuchen Wir, durch das Innere des Motorrads zu führen, um dem Arduino MKR Motor Carrier anzuschließen. Der Akku rutscht leicht vom Motorkörper. Versuchen wir, mit einem Gummi den Akku mit Motorkörper befestigen.

Aufbau des Systems

Bild2: Darstellung des Motorrades. Abbildung zeigt die Komponenten des Mortorrades Junjie Lyu (Diskussion) 21:34, 11. Okt. 2021

Die Komponenten sind in der Bild2 dargestellt. Die Basis dieses Projekts bildet das Arduino-Board "MKR1000". MKR Motor Carrier ist eine MKR-Zusatzplatine für "MKR1000" zur Steuerung von Servo-, Gleichstrommotoren. Außerdem erweitert MKR Motor Carrier die Fähigkeiten von MKR1000 und vereinfacht die Anschluss zur anderen Aktoren und Sensoren über ein Reihe von 3-poligen Stiftleisten. IMU Sensor enthält verschiedene Sensoren, Beschleunigungsmesser, Gyrokope und Magnetometer, in einem einzigen Gehäuse. Mit IMU Sensor, der auf dem MKR1000 sitzt, wird die vertikale Position des selbstbalancierenden Motorrad gemisst und erkannt, wenn es fällt und der Kommuniziert mit MKR1000 via 122C. Hallsensor misst die Geschwindigkeit vom Trägheitsrad. Encoder misst die Geschwindigkeit vom Motorrad. Ultraschllsensor erkannt die Hindernisse vor dem Mortorrad. Servo-Motor ändert die Fahrrichtung des Motorrades. Das Motorrad nutzt Simulink, um die Sensoren, Aktuatoren und die Bewegung zu überwachen und zu steuern.

Identifikation der Regelstrecke

Bild 3: Beschreibung des Motorrades durch physikalische Größe. Junjie Lyu (Diskussion) 10:49, 12. Okt. 2021

Bild 4.1: Beobachtung der Regelstrecke.Junjie Lyu (Diskussion) 23:03, 22. Okt. 2021

Bild 4.2: Beobachtung der Regelstrecke.Junjie Lyu (Diskussion) 11:36, 23. Okt. 2021

Tabelle 1: Bedeutung dieser physikalischen Größe in Bild 3

physikalische Größe	Bedeutung
$\theta$	der Neigungswinkel $\theta$ ist 0 Grad , wenn das Motorrad vollkommen aufrecht steht, ist positiv, wenn das Motorrad von hinten gesehen gegen den Uhrzeigersinn geneigt ist, und negativ, wenn das Motorrad im Uhrzeigersinn geneigt ist.
$\phi$	die Rotationsverschiebung des Trägheitsrads relativ zum Rest des Motorrads ist $\phi$, wobei eine positive Verschiebung als gegen den Uhrzeigersinn definiert ist.
$h_{cm}$	Die Höhe des Massenschwerpunkts über dem Boden bei aufrechtem Motorrad ($\theta$ = 0) ist definiert als $h_{cm}$

In diesem Projekt selbstfahrendes Motorrad ist die Regelstrecke der Neigungswinkel $\theta$. Beobachten Wir unter den Bedingungen ohne externe Störung und Drehmoment vom DC Motor, wie die Regelstrecke reagiert. Im Bild4.1 und Bild4.2 werden die Ergebnisse beschreibt. Figur 1 und Figur 2 sind die Rotationsgeschwindigkeiten von Motorrad und Trägheitrad unter den ursprünglichen Bedingungen $\theta_0 = 0$ und $\dot\theta_0 = 0$. Figur 3 und Figur 4 sind die Rotationsgeschwindigkeiten von Motorrad und Trägheitrad unter den ursprünglichen Bedingungen $\theta_0 = 1$ und $\dot\theta_0 = 0$. Bild4.2 ist die Antwort von $\theta$ unter den ursprünglichen Bedingungen $\theta_0 = 0$ und $\dot\theta_0 = 0$. Wir können die folgende Ergebinisse bemerken: 1) Die Signals sind wegen des Gravitationsmomentes oszillierend. 2) $\dot\theta$ und $\dot\phi$ haben gleiche Frequenz und ähnliche Form aber invertiert. 3) Beim zweiten Durchlauf dauert es einige Zeit, bis ein stabiler Zustand erreicht wird, bei $\theta$ in einen rein oszillierenden Zustand übergeht, während sich das Signal beim ersten Durchlauf direkt in diesem Zustand befindet. 4) Bei stabilem Zustand ist das Signal vom zweiten Durchlauf oszillierend im Intervall ingefähr [360 720], aber das Signal vom ersten Durchlauf im Intervall ingefähr [0 360]. In der realen Welt wird dies aufgrund des physischen Bodens nicht passieren. wir konzentrieren uns auf das Verhalten des Systems, wenn der Neigungswinkel des Motorrads nahe 0 ist.

Mathematisches Modell des Systems

Das Drehmoment beim Motorrad $\tau_{net,M}$ hat 3 Hauptkomponenten, das auf das Motorrad um die Bodenachse des Rades wirkt:

$\rightarrow$ Gravitationsmoment $\tau_{g,M}$

$\rightarrow$ Trägheitsmoment des Trägheitsrads $\tau_{IW,M}$

$\rightarrow$ Externes Drehmoment $\tau_{ext,M}$

In dieser theoretischen Diskussion werden wir uns auf das ideale Szenario konzentrieren, bei dem $\tau_{ext,M} = 0$. Dissipative Kräfte wie Reibung und Widerstandseffekte zwischen beweglichen Teilen $\tau_{fric,IM} = 0$. Die Neigungswinkel $\theta$ ist sehr klein (nur wenige Grad), gilt für $\sin\theta \approx \theta$.

Nettodrehmoment beim Motorrad für Bodenradachsen: $\tau_{net,M} = I_M \cdot \ddot{\theta} = \tau_{g,M} + \tau_{IW,M} + \tau_{ext,M} \approx \tau_{g,M} + \tau_{IW,M} \qquad(1)$

Drehmoment beim Trägheitrad für Drehachse durch die Motorwelle des Trägheitsrades: $\tau_{net,IW} = I_{IW} \cdot \ddot{\phi} = \tau_{motor,IW} + \tau_{fric,IW} \approx \tau_{motor,IW} \qquad(2)$

Gravitationsmoment: $\tau_{g,M} = M_M \cdot g \cdot h_{cm} \cdot \sin\theta \qquad(3)$

Das von der Motorwelle auf das Trägheitrad ausgeübte Drehmoment ist gleich groß und entgegengesetzt zu dem vom Trägheitrad auf die Motorwelle ausgeübten Drehmoment: $\tau_{IW,M} = - \tau_{motor,IW} \qquad(4)$

Ersetze Gleichung (4) in Gleichung (1): $\tau_{net,M} = I_M \cdot \ddot{\theta} \approx \tau_{g,M} - \tau_{motor,IW} \qquad(5)$

Ersetze Gleichung (3) in Gleichung (5): $I_M \cdot \ddot{\theta} \approx M_M \cdot g \cdot h_{cm} \cdot \sin\theta - \tau_{motor,IW} \qquad(6)$

Ersetze Gleichung (2) in Gleichung (6): $I_M \cdot \ddot{\theta} \approx M_M \cdot g \cdot h_{cm} \cdot \sin\theta - I_{IW} \cdot \ddot{\phi} \qquad(7)$

Ersetze Gleichung $\sin\theta \approx \theta$ in Gleichung (7):$I_M \cdot \ddot{\theta} \approx M_M \cdot g \cdot h_{cm} \cdot \theta - I_{IW} \cdot \ddot{\phi} \qquad(8)$

wir programmieren den Motor so, dass er ein Drehmoment aufbringt, das proportional zum Neigungswinkel selbst ist. Das heißt: $\tau_{motor,IW} = K_p \cdot \theta$. Ersetze Gleichung $\tau_{motor,IW} = K_p \cdot \theta$ in Gleichung (8): $I_M \cdot \ddot{\theta} \approx M_M \cdot g \cdot h_{cm} \cdot \theta - K_p \cdot \theta \qquad(9)$
Von (9) erhalten wir: $I_M \cdot \ddot{\theta} \approx - (K_p-M_M \cdot g \cdot h_{cm})\theta \qquad(10)$

Wenn $K_p$ größer als $M_M \cdot g \cdot h_{cm}$ ist, dann geht die Winkelbeschleunigung des Motorrads um die Rad-Boden-Achse $\ddot{\theta}$ in die entgegengesetzte Richtung des Neigungswinkels $\theta$. Das heißt, wenn sich das Motorrad in eine Richtung neigt, wird es in die entgegengesetzte Richtung beschleunigt. Wenn sich das Motorrad auf die andere Seite des Gleichgewichtspunktes bewegt, wechselt die Winkelbeschleunigung die Richtung, um das Motorrad wieder in Richtung Gleichgewicht zu bewegen. Es handelt sich um ein stabiles Gleichgewicht, d. h., wenn das System aus der Gleichgewichtslage heraus gestört wird, stellt es sich selbst wieder ins Gleichgewicht zurück.

Bild5: Lösung der Differentialgleichung von theta. Abbildung zeigt die Prozesse der Lösungsfindung der Differentialgleichung von theta Junjie Lyu (Diskussion) 20:43, 22. Okt. 2021

die Differentialgleichung (10) löst: $\theta(t) = A\sin(\sqrt{\frac{K_p-M_M \cdot g \cdot h_{cm}}{I_M}}t+B)$

In der obigen Gleichung sind $A$ und $B$ Konstanten, die von den Anfangswerten von $\theta$ und $\dot{\theta}$ abhängen. $I_M$ ist das Trägheitsmoment des Motorradsystems um die Rad-Boden-Achse. Wenn die Konstante $K_p$ erhöht wird, wird die Schwingungsfrequenz schneller. Wir wollen dieses Verhalten so verbessern, dass die Amplitude der Schwingungen abnimmt und sich schließlich auf einen konstanten Neigungswinkel von 0 einpendelt und das Motorradsystem stabil auf zufällige Störeinflüsse anspricht.

Bewertung verschiedener Regleransätze

Bild 6: Regelkreis des Motorrades.Junjie Lyu (Diskussion) 12:19, 22. Okt. 2021

Tabelle 2: Bedeutung dieser physikalischen Größe in Bild 6

physikalische Größe	Bedeutung
$W$	Sollwert. Wenn das Motorrad ganz aufrecht steht, ist der Neigungswinkel "0". Das heißt: $W = 0$
$e$	Regelabweichung. Wir hoffen $e = 0$. Das heißt: $e = W - r = 0$.
$y$	Stellgröße des Reglers. Unter Berücksichtigung der dynamischen Eigenschaften der Regelstrecke wird Stellgröße des Reglers nach Regelabweichung bestimmt.
Fehler beim Parsen (Konvertierungsfehler. Der Server („cli“) hat berichtet: „[INVALID]“): {\displaystyle y^'}	Stellgröße. DC motor wirkt sich auf das Trägheitsrad zur Ausgleichung der Abweichung.
$z$	Störung. USB-Kabel gegen Motorräder, Luftwiderstand, Wind und verschiedene andere Ablenkungen
$y_{ges}$	Gesamtregelgröße. Die Gesamtregelgröße auf dem System.
$x$	Istwert. Der Neigungswinkel auf dem System
$r$	Rückführgröße. Die Größe des Neigungswinkels auf dem System wird von Sensoren gemisst.

Eine Bewertung erfolgt qualitativ auf Grundlage der Eigenschaften der drei Regler in Tabelle 3.

Tabelle 3: Bewertung verschiedener Regleransätze

Regler	Bewertung
P	Ein P-Regler erreicht in unserem Fall keine Regeldifferenz von 0. Die bleibende Regelabweichung kann durch Verringerung des Parameters $K_P$ verringert werden. Wenn $K_P$ zu niedrig eingestellt wird, reicht das Drehmoment des Motors nicht aus, um das Motorrad gegen das Gravitationsmoment im Gleichgewicht zu halten.
PD	Das Problem der propotionalen Regler, die bleibende Regelabweichung, ist beim PD-Regler allerdings weiterhin vorhanden.Entweder wird $K_P$ erhöht oder einen Integralanteil $K_i$ in den Regler eingefügt, um das Problem zu lösen.
PID	Der PID-Regler erweitert den PD-Regler durch einen Integralanteil, der die Regelabweichung über der Zeit aufsummiert und die Summe mit dem Faktor $K_i$ multipliziert. Dabei wird die Abweichung vollständig eliminiert.

In dem Projekt wurde ein PD-Regler verwendet.Die Anhäufung von Fehlern im Laufe der Zeit kann sich nachteilig auswirken,um das Motorrad später auszubalancieren. Deshalb wird auf ein I-Glied verzichtet.

Simulation des Reglers

Bild 7.1: Simulation des Reglers. Junjie Lyu (Diskussion) 12:57, 24. Okt. 2021

Bild 7.2: Kontroller des Systems.Junjie Lyu (Diskussion) 13:00, 24. Okt. 2021

In der folgenden Tabelle 4 sind die wichtigsten Parameter für die Ausführung des Modells erklärt. Einige Parameter gehen aus den Aufgaben der AEK-Seite hervor und werden mithilfe dieser berechnet. Die Bennung von theta0 und thetadot0 wurde in dem Fall beibehalten.

Tabelle 3: Parameterliste des Modells aus der start.m Datei

Parameter	Beschreibung	Einheit	Wert
theta0	Aus Aufgabe. Als t = 0, ist der Neigungswinkel theta0.	deg	1
thetadot0	Aus Aufgabe. Als t = 0, ist die Rotationsgeschwindigkeit des Motorrades thetadot0.	deg/s	0
rs_EM_s32	Externes Drehmoment z.B Wind	N*m	in der Simulation wird es durch Block "Random Number" erzeugt
rs_Theta_s32	der Neigungswinkel im Durchlauf	deg	-
rs_ThetaDot_s32	die Rotationsgeschwindigkeit im Durchlauf	deg/s	-
rs_MM_s32	das Moment vom DC-Motor	N*m	-
rs_TRGesch_s32	Rotationsgewschwindigkeit des Trägheitsrades	deg/s	-
par_Kp	Parameter für das P-Glied	-	0.5
par_Kd	Parameter für das D-Glied	-	0.08

Video

Verlinken Sie hier ein YouTube-Video zu Ihrem fertigen Projekt. Tipps zum Video finden Sie hier.

Weblinks

1. Arduino Engineering Kit
2. YouTube: Unboxing the Arduino Engineering Kit
3. YouTube: Motorcycle Self Balancing Using Reinforcement Learning
4. YouTube: Motorcycle Maneuvers
5. Arduino Store

Software

1. Arduino Engineering Kit Hardware Suppo​rt
2. Arduino Engineering Kit Hardware Suppo​rt für R2018b
3. Arduino Engineering​ Kit Project Files
4. Reinforcement learning with Self-balancing motorcycle

Siehe auch

1. Studentische Arbeiten bei Prof. Schneider
2. Anforderungen an eine wissenschaftlich Arbeit

---

→ zurück zum Hauptartikel: Studentische Arbeiten