Co-Simulation wird in vielen Bereichen verwendet, um zwei oder mehrere Subsystem-Solver im Zeitbereich zu koppeln. Bei mechanischen Systemen gibt es im Wesentlichen zwei Kopplungsarten, nämlich Kopplungsverfahren auf der Basis von Konstitutivgleichungen und Kopplungsansätze basierend auf algebraischen Bedingungsgleichungen. In dieser Dissertation werden die numerische Stabilität und das Konvergenzverhalten von verschiedenen Co-Simulationsverfahren analysiert. Für die Stabilitätsanalyse muss ein Testmodell definiert werden. In Anlehnung an die Stabilitätsdefinition für numerische Zeitintegrationsschemata (Dahlquistsche Stabilitätstheorie) werden lineare Co-Simulationstestmodelle definiert. Als Co-Simulationstestmodelle werden hier Zweimassenschwinger eingesetzt, wobei die Massen über ein Feder-Dämpfer-Element oder über ein starres Gelenk verbunden sind. Durch Diskretisierung des Testmodells mit einem Co-Simulationsverfahren können Rekurrenzen-Gleichungen abgeleitet werden, die die zeitdiskrete Co-Simulationslösung beschreiben. Das Stabilitätsverhalten des linearen Zweimassenschwingers ist bei Anwendung von Kopplungsverfahren auf Basis eingeprägter Kopplungskräfte/-momente durch 7 unabhängige Parameter charakterisiert. Um verschiedene Co-Simulationsansätze zu vergleichen, sind 2D-Stabilitätsplots praktikabel. Daher werden 5 der 7 Parameter fixiert, so dass der Spektralradius als Funktion der verbleibenden 2 Parameter dargestellt werden kann. Die Ergebnisse zeigen erwartungsgemäß, dass implizite Kopplungsschemata eine wesentlich bessere numerische Stabilität als explizite Methoden aufweisen. Darüber hinaus kann ein verbessertes Stabilitätsverhalten durch Erweiterungen der Kopplungsbedingungen erreicht werden, d.h. durch Berücksichtigung von Ableitungen und Integralen der konstitutiven Gleichungen. Ein sehr gutes Stabilitätsverhalten ist mit dem D-erweiterten Kraft/Kraft-Kopplungsansatz in Verbindung mit quadratischen Approximationsfunktionen erreichbar. Die Analyse der numerischen Stabilität von Co-Simulationsverfahren mit algebraischen Bedingungen ist das zweite zentrale Thema dieser Arbeit. Hier müssen 5 unabhängige Parameter für das lineare Testmodell definiert werden. Der dimensionslose Real- und Imaginärteil des Eigenwerts von Subsystem 1 werden als Achsen in 2D-Stabilitätsplots verwendet; die anderen 3 Parameter werden konstant gehalten. Drei klassische Stabilisierungsmethoden, nämlich Baumgarte-Stabilisierung, Methode der gewichteten Multiplikatoren und Koordinatenprojektion, werden für unterschiedliche Approximationsordnungen diskutiert. Alternativ kann der Lagrange-Multiplikator zwischen den Makrozeitpunkten diskretisiert werden (erweiterter Multiplikatorenansatz), um Co-Simulationsansätze auf Index-2 und auf Index-1-Ebene zu erhalten. Dadurch können gleichzeitig die Kopplungsbedingungen und deren Zeitableitungen an den Makrozeitpunkten erfüllt werden. Es werden mehrere Mehrkörpermodelle verwendet, um die praktische Anwendung der erwähnten Co-Simulationsverfahren aufzuzeigen.