Seminaraufgabe SoSe 2021: Einspurmodell Gruppe B

Autoren: Björn Schlottke und Dennis Schleicher

Einleitung

Dieser Artikel dokumentiert die Ausarbeitung der Gruppe B im Modul Systems Design Engineering im Sommersemester 2021. Die Aufgabe beinhaltet die Entwicklung des linearen Einspurmodells nach dem V-Modell als Tool zur Komplexitätsbeherrschung. Hierfür wurden die Inhalte des Skriptes Datei:Skript Fahrwerkmanagent ProfGöbel.pdf in Submodule und zugehörige Komponenten aufgeteilt und in einem MATLAB 2020a Simulink Modell umgesetzt.

Vorgehensweise nach V-Modell

Wie bereits erläutert wurde die Entwicklung nach dem V-Modell durchgeführt. Dies beinhaltet die in die folgenden Schritte, welche an dieser Stelle noch einmal kurz aufgelistet sind.

1. Anforderungsdefinition
2. Funktionaler Systementwurf
3. Technischer Systementwurf
4. Komponentenspezifikation
5. Programmierung
6. Komponententest
7. Integrationstest
8. Systemtest
9. Abnahmetest

Die Schritte 6 - 9 werden in dem vorliegenden Artikel nicht behandelt, weil sie die Testung der Ausarbeitungen von Gruppe J sowie den Abnahmetest des entwickelten Modells durch Prof. Göbel, betreffen. Somit wird lediglich die Entwicklung des Simulinkmodells behandelt.

Anforderungen

Zunächst wurden alle Anforderungen zu denen sich Gruppe B im Rahmen dieser Seminararbeit verpflichtet hat, in einem Lastenheft festgehalten. Hierfür wurde das Lastenheft mit der folgenden Struktur versehen.

1. Formale Anforderungen
   1. Nachhaltigkeit des Programmcodes
   2. Abgaben des letzten Meilensteins im HSHL Wiki
2. Termine und Fristen
3. Aufbau
   1. Teilmodul Mannöversteuerung
   2. Teilmodul Modell
   3. Teilmodul Auswertung
4. Modellgrundlagen

Das gesamte Lastenheft kann im SVN eingesehen werden. In Tabelle 1 wird ein Auszug aus dem Lastenheft dargestellt. Die Abkürzungen in der Spalte Typ stehen für Inhalt (= I) und Anforderung (= A).

Tabelle 1: Formale Anforderungen aus dem Lastenheft

ID	Typ	Kapitel	Inhalt
001	I	1	Formale Anforderungen
002	A		Es muss die Software MATLAB 2020a benutzt werden.
003	A		Die Vorgehensweise muss dem V-Modell nach Vorlage von Prof. Göbel entsprechen.
004	A		Es muss ein lineares Einspurmodell in MATLAB Simulink erstellt werden.
005	A		Für die Erstellung des Modells müssen die Parameter extrahiert und in MATLAB angelegt werden.
006	A		Die Parameter müssen nachträglich anpassbar sein.

Technischer Systementwurf

Im folgenden Kapitel wird der gesamte Entwurf des Systems erläutert. An dieser Stelle wird der funktionale Systementwurf aus dem V-Modell übergangen, da er keinen Teil der Seminaraufgabe dargestellt hat. Grundlegend ist das Simulink Modell so aufgebaut, dass es drei Submodule besitzt, welche wiederum weitere Komponenten beinhalten. Innerhalb der Komponenten ist die eigentliche Funktionalität hinterlegt. Diese Struktur wurde für die Bearbeitung der Seminaraufgabe vorgegeben und wurde ebenfalls unter den Anforderungen im Lastenheft festgehalten. Die oberste Ebene dieser Struktur lässt sich in Abbildung 1 erkennen.

Abb. 1: Submodule des Simulinkmodells

Das Submodul Mannöversteuerung besitzt nur eine Komponente, nämlich die Mannöverauswahl. Unter dem Submodul Modell hingegen befinden sich die Komponenten Reifen, Karosserie, Gierdynamik und Schwimmwinkel. Als Letztes besitzt das Submodul Auswertung keine Komponenten. Weitergehend kann erkannt werden, dass die Submodule untereinander Verbindungen besitzen. Somit werden die Parameter aus Tabelle 2 zwischen den Submodulen übergeben.

Tabelle 2: Übergebene Parameter der Submodule

Parameter	Beschreibung	Einheit
$\delta_v$	Lenkwinkel vorne	$rad$
$\delta_h$	Lenkwinkel hinten	$rad$
$_Kv_{Cx}$	Fahrzeuglängsgeschwindigkeit im Fahrzeugschwerpunkt C	$\frac{m}{s}$
$_Ka_{Cy}$	Fahrzeugquerbeschleunigung im Fahrzeugschwerpunkt C	$\frac{m}{s^2}$
$_Ka_{Cx}$	Fahrzeufglängsbeschleunigung im Fahrzeugschwerpunkt C	$\frac{m}{s^2}$
$_KF_{Cy}$	Fahrzeugquerkraft im Fahrzeugschwerpunkt C	$N$
$_KF_{Cx}$	Fahrzeuglängskraft im Fahrzeugschwerpunkt C	$N$
$\beta$	Schwimmwinkel	$rad$
$\beta_v$	Schwimmwinkelanteil vorne	$rad$
$\beta_h$	Schwimmwinkelanteil hinten	$rad$
$_R\alpha_v$	Schräglaufwinkel vorne	$rad$
$_R\alpha_h$	Schräglaufwinkel hinten	$rad$
$_K\dot{\psi}$	Gierwinkelrate (Gierwinkelgeschwindigkeit)	$\frac{rad}{s}$

Weitergehend wurden die Inhalte der Submodule einmal dargestellt. Hierbei handelt es sich um die Ebene auf denen die Komponenten erkennbar sind.

Abb. 3: Submodul Modell

Es lässt sich direkt erkennen, dass das Submodul Modell die eigentliche Funktionalität des linearen Einspurmodells besitzt und diese entsprechend hinterlegt.

Komponentenspezifikation

Die Komponentenspezifikation definiert die Aufgabe, das Verhalten, den inneren Aufbau und die Schnittstellen der Komponenten. In den Modulen "Mannöversteuerung" und "Auswertung" werden keine Mathematischen Formeln verwendet (bei der "Mannöversteuerung" erfolgt eine Zuweisung des hinteren Lenkwinkels über eine Multiplikation mit 0, diese ist aber trivial). Das Modul "Modell" wird, wie im technischen Systemplan beschrieben, in vier Komponenten unterteilt. Deren Funktionsweise soll im Folgenden mithilfe mathematischer Zusammenhänge erläutert werden. Bei der Herleitung wurden die in beschriebenen Formeln verwendet.

Schwimmwinkel

Die Formel der Winkelgeschwindigkeit des Schwimmwinkels wird aus dem vereinfachten Schwerpunktsatz abgeleitet. Dabei ist $_KF_{Cy}$ die im Schwerpunkt der Hinterachse angreifende Kraft, m die Fahrzeugmasse und $_K\dot{\psi}$ die Gierwinkelgeschwindigkeit.

$_K\dot{\beta} = \frac{_KF_{Cy}}{m \cdot _Kv_{Cx}} - _K\dot{\psi}$

Auf den Hintergrund der Vereinfachung des Schwerpunktsatzes wird im Rahmen der entsprechenden Komponente näher eingegangen.

Durch Integration von $_K\dot{\beta}$ ergibt sich der Schwimmwinkel des Fahrzeugs $_K\beta$. Der Schwimmwinkel der Vorder- und Hinterachse $\beta_v$ und $\beta_h$ kann jeweils aus dem Schwimmwinkel des Fahrzeugs $_K\beta$, dem Schwerpunktabstand zur Vorder- bzw. Hinterachse $l_v$ bzw. $l_h$, der Gierwinkelgeschwindigkeit $_K\dot{\psi}$ und der Fahrzeuglängsgeschwindigkeit $_Kv_{Cx}$ berechnet werden.

$\beta_v = \beta + \frac{l_v \cdot _K\dot{\psi}}{_Kv_{Cx}}$

$\beta_h = \beta - \frac{l_h \cdot _K\dot{\psi}}{_Kv_{Cx}}$

Räder

Der Achs-Schräglaufwinkel ist der Winkel zwischen der Reifenlängsachse und dem Geschwindigkeitsvektor des Radmittelpunkts. Daher ergeben sich die Achs-Schräglaufwinkel der Vorder- und Hinterachse $\alpha_v$ und $\alpha_h$ aus der Differenz des jeweiligen Lenk- und Schwimmwinkels.

$\alpha_v = \delta_v - \beta_v$

$\alpha_h = \delta_h - \beta_h$

Basierend auf den Achs-Schräglaufwinkeln lassen sich zusammen mit der jeweiligen Achssteifigkeit $c_v$ und $c_h$ die Querkräfte an den Rädern bzw. Achsen berechnen.

$_RF_{yv} = c_v \cdot \alpha_v$

$_RF_{yh} = c_h \cdot \alpha_h$

Karosserie

Nach dem Skript wird beim linearen Einspurmodell von kleinen Winkeln ausgegangen und damit ergeben sich damit folgende Formeln für die Kräfte in x- und y-Richtung.

$_KF_{xv} \approx \delta_v \cdot _RF_{yv}$

$_KF_{xh} \approx \delta_h \cdot _RF_{yh}$

$_KF_{yv} \approx _RF_{yv}$

$_KF_{yh} \approx _RF_{yh}$

Aus den Kräften in x- und y-Richtung werden jeweils die Kräftesummen $_KF_{Cx}$ und $_KF_{Cy}$ gebildet. Durch einsetzen können die folgenden Formeln gebildet werden.

$_KF_{Cx} = \delta_v \cdot _RF_{yv} + \delta_h \cdot _RF_{yh}$

$_KF_{Cy} = _RF_{yv} + _RF_{yh}$

Durch das Dividieren der jeweiligen Kräfte (in x- und y-Richtung) mit der Masse des Fahrzeugs ergeben sich die Fahrzeuglängsbeschleunigungen in x- und y-Richtung.

$_Ka_{Cy} = \frac{_KF_{Cy}}{m}$

$_Ka_{Cx} = \frac{_KF_{Cx}}{m}$

Der Schwerpunkt Abstand der Hinterachse $l_h$ ergibt sich aus der Differenz von dem Radstand $l$ minus dem Schwerpunkt Abstand zur Vorderachse $l_v$.

$l_h = l - l_v$

Mithilfe des Drallsatzes lässt sich die Gierwinkelbeschleunigung $_K\ddot{\psi}$ berechnen. Dabei wird die Differenz der Produkte aus dem Schwerpunktabstand zur Vorder- bzw. Hinterachse $l_v$ bzw. $l_h$ und den jeweiligen Querkräften durch die Massenträgheit $J_{zz}$ geteilt.

$_K\ddot{\psi} = \frac{l_v \cdot _KF_{yv} - l_h \cdot _KF_{yh}}{J_{zz}}$

Gierwinkel

Durch Integration der Gierwinkelbeschleunigung $_K\ddot{\psi}$ ergibt sich die Gierwinkelgeschwindigkeit $_K\dot{\psi}$.

Programmierung

Die Programmierung beinhaltetet zwei unterschiedliche Komponenten. Die Simulation erfolgt über das Simulink Modell und die Parameter werden vom Parameterskript aus MATLAB übergeben.

Parameterskript

Das Parameterskript übergibt zunächst alle im Modell benötigten Parameter:

Abb. 4: Übergebene Parameter aus dem Skript

Es werden 4 Fahrmannöver definiert, welche die Nummer 1-4 zugewiesen wird. Diese sind:

1. Impulsartiges Einlenken (Lenkwinkelsprung)
2. Slalomfahrt
3. Keine Lenkbewegung
4. Gleichmäßiges Einlenken

Zudem wird zwischen Links-und Rechtskurve unterschieden.

Es werden Benutzerabfragen eingebaut, welche die Interaktion mit dem Skript bestimmen.

Die Eingaben vom Benutzer werden nur dann übernommen wenn es sich um Eingaben innherb des Definitionsbereiches handelt.Bei den Fahrmannövern ist dieser Definitionsbereich, wie oben zu sehen 1-4. Der Benutzer kann eine freie Eingabe tätigen, wenn nun eine 5 oder eine andere Zahl/Buchstabe außerhalb des Definitionsraumes eingeben wird bekommt der Nutzer eine Fehlermeldung und die Ausführung des Skriptes wird abgebrochen.

Die gleiche Vorgehensweise wurde für die Auswahl einer Links- oder Rechtskurve implementiert.

Außerdem wird eine dritte Eingabe vom Benutzer gefordert, welche Abfragt ob er eine Darstellung des Simulationsergebnisses wünscht oder nicht. Der gesamte MATLAB-Code des Parameter Modells kann im SVN eingesehen und ausprobiert werden.

Simulink Modell

Im Folgenden wird die Umsetzung der im Kapitel Komponentenspezifikation beschriebenen Gleichungen erläutert. Außerdem wird ebenfalls auf die Simulationsmöglichkeiten der unterschiedlichen Fahrmannöver eingegangen, die bereits im Parameterskript abgefragt werden. Um die Funktionsweise des Simulinkmodells verstehen zu können, werden dementsprechend die Komponenten abgebildet und anschließend erläutert.

Komponente Mannöverauswahl

In der Abbildung lässt sich dementsprechend erkennen, dass die Komponente eine Switch Abfrage für die Variable Mannöver ID beinhaltet. Aufgrund dieser ID werden dann die unterschiedlichen Lenkwinkel Verläufe simuliert. Außerdem wird der Lenkwinkel hinten $\delta_h$ gebildet und zwar mit einem relativen Anteil von dem Lenkwinkel vorne $\delta_v$. Im Nachgang wurde dieser realtive Anteil gleich 0 % gesetzt, da kein Lenkwinkel an der Hinterachse existiert.

Die nächste Komponente ist der Reifen, welche die Umsetzung der Gleichung für die Achsschräglaufwinkel $\alpha_v$ und $\alpha_h$ beinhaltet. Außerdem werden den Reifen und somit Achsen, wirkenden Querkräfte $_RF_{yv}$ und $_RF_{yh}$ berechnet. Für die Umsetzung wurden einfache Additions- und Multiplikationsblöcke verwendet.

Komponente Reifen

Weitergehend beinhaltet die Komponente Karosserie die Berechnungen für die Fahrzeugquer- und längskraft $_KF_{Cy}$ und $_KF_{Cx}$ sowie die Fahrzeugquer- und längsbeschleunigung $_Ka_{Cy}$ und $_Ka_{Cx}$ im Fahrzeugschwerpunkt C. Außerdem wird noch die Gierwinkelbeschleunigung $_K\ddot{\psi}$ berechnet.

Komponente Karosserie

Die Komponente Gierdynamik hingegen beinhaltet nur eine einzige Berechnung, nämlich die Integration mit der Gierwinkelbeschleunigung $_K\ddot{\psi}$ zur Gierwinkelrate, bzw. Gierwinkelgeschwindigkeit $_K\dot{\psi}$.

Komponente Gierdynamik

Als letzte Komponente wurde die Komponente Schwimmwinkel verwendet. Diese berechnet den Schwimmwinkel $\beta$ mit seinen Anteilen für die Vorder- und Hinterachse $\beta_v$ und $\beta_h$.

Komponente Schwimmwinkel

Dokumentation

Link zu der Ausarbeitung in SVN

---

→ zurück zum Hauptartikel: "Systems Design Engineering" - Seminaraufgabe SoSe 2021: Einspurmodell“