Co-Simulationsverfahren können nicht nur im Bereich der Simulation multiphysikalischer Systeme, sondern auch zur Parallelisierung der Berechnung großer monodisziplinärer Modelle verwendet werden. Diese Arbeit beschäftigt sich mit der Reduktion der Rechenzeit von Simulationen großer mehrkörperdynamischer Systeme, die durch das Aufteilen eines monolithischen Modells in gekoppelte Subsysteme erzielt werden kann. Die Kopplung zwischen den Subsystemen ist problemabhängig und kann auf verschiedene Arten beschrieben werden. Bei den im Rahmen dieser Arbeit verwendeten Verfahren, sind die Subsysteme durch nichtlineare konstitutive Gleichungen miteinander verbunden. Informationen in Form von Kopplungsvariablen werden nur an bestimmten Kommunikations- oder Makrozeitpunkten zwischen den Subsystemen ausgetauscht. Zwischen den Makrozeitpunkten werden die Kopplungsvariablen in den Subsystemen durch Extra-/Interpolationspolynome approximiert. Diese Vorgehensweise erlaubt es, die Subsysteme, bzw. deren Bewegungsgleichungen, in jedem Makrozeitschritt unabhängig voneinander zu integrieren. Wenn ein Co-Simulationsverfahren vom Jacobi-Typ verwendet wird, werden alle Subsysteme parallel gelöst. Ein wesentlicher Nachteil vieler verfügbarer Co-Simulationsansätze ist die Verwendung eines äquidistanten Makrozeitgitters. In praktischen Anwendungsfällen, insbesondere wenn hoch nichtlineare Systeme oder Modelle mit stark veränderlichen Eigenschaften simuliert werden, ist die Anwendung einer konstanten Makrozeitschrittweite oftmals nicht rechenzeiteffizient. Im Rahmen dieser Arbeit wird ein explizites und ein implizites Co-Simulationsverfahren vorgestellt, die beide über einen Algorithmus zur Steuerung der Makroschrittweite verfügen. Neben der Makroschrittweite wird durch den Algorithmus auch der Polynomgrad der Approximationspolynome für die Kopplungsvariablen geregelt. Die Schrittweitenkontrolle basiert auf verschiedenen Fehlerschätzern für den lokalen Fehler eines Makrozeitschritts, die Ordnungssteuerung basiert auf dem Wert einer Norm der skalierten Ableitungen der Approximationspolynome. Diese Vorgehensweise orientiert sich an den Methoden, die in numerischen Differentialgleichungslösern, wie z. B. BDF-Implementierungen, verwendet werden. Die parallele Implementierung der Co-Simulationsverfahren erfolgt durch einen hybriden MPI (Message Passing Interface) und OpenMP (Open Multi-Processing) Ansatz. Jedes Subsystem wird von einem MPI-Prozess ausgeführt; falls mehrere Integrationen desselben Subsystems notwendig sind, beispielsweise zur Fehlerschätzung oder im Rahmen der Gradientenberechnung für das implizite Co-Simulationsverfahren, werden diese mittels OpenMP parallelisiert. Die Rechenzeiteffizienz des expliziten und des impliziten Co-Simulationsverfahrens mit und ohne Schrittweitensteuerung wird ausführlich anhand zahlreicher numerischer Beispiele untersucht. Die Simulationen wurden alle vollständig parallelisiert auf dem Lichtenberg Hochleistungsrechner der Technischen Universität Darmstadt ausgeführt. Der Einfluss verschiedener Verfahrensparameter auf die Rechenzeit der maßgeblichen Prozesse, die innerhalb der Co-Simulationsverfahren ausgeführt werden, wird detailliert analysiert. Die Testsimulationen verdeutlichen die Steigerung der Effizienz der Co-Simulationsverfahren, die durch die Implementierung des Makroschrittweiten- und Ordnungssteuerungsalgorithmus erzielt wird. Außerdem wird die drastische Rechenzeitersparnis aufgezeigt, die durch die Verwendung von Co-Simulationsverfahren im Vergleich zu einer monolithischen Simulation großer Mehrkörpersysteme ermöglicht wird. Die erzielte Rechenzeitersparnis lässt sich durch eine Kombination aus Parallelisierung und vorteilhafter Ausnutzung des Multirate Attributs der Co-Simulationsverfahren erklären.