Energy-shaping techniques are used to expand the range of autonomous motion of unmanned aerial systems without prohibitively {color{black}increasing the computational cost of the resultant controller}. Passivity-based control presents a method to implement a static, nonlinear state feedback control law that stabilizes the motion of an aircraft with a large region of attraction. {color{black} The energy-based control scheme is applied to both multirotor and fixed-wing aircraft}. Multirotor aircraft dynamics are cast into a port-Hamiltonian System and the concept of trajectory tracking using canonical feedback transformation is implemented to construct a cross-track controller. Fixed-wing aircraft dynamics are cast in port-Hamiltonian form and a passivity-based nonlinear control law for steady, wings-level flight of a fixed-wing aircraft to a specified inertial velocity (speed, course, and climb angle) is constructed. Results in simulations and experiments suggest robustness, and a large region of attraction of the controller. The control law extended to support time-varying inertial velocity tracking that incorporates banking to turn. The results are extended by including a line-of-sight guidance law and varying the direction as a function of position relative to a desired path, rather than as a function of time. The control law and the associated proof of stability follow similarly to that of the time-varying directional stabilization problem. The results are supported with simulations as well as experimental flight tests. / Doctor of Philosophy / This dissertation presents an alternative but intuitive approach to regulate unmanned aerial vehicles' flight that would allow for more maneuverability {color{black} than conventional methods}. This scheme relies on modifying the energy of the system to achieve the desired motion and leverages the properties of the aircraft rather than eliminating them and imposing different properties. This approach is applied to both fixed-wing and aircraft and quadcopters. Simulations and experimental flights have show the efficacy of this approach compared to other more established methods.

Cette thèse est motivée par un problème aéronautique: le ballottement du carburantdans des réservoirs d’ailes d’avion très flexibles. Les vibrations induites par le couplagedu fluide avec la structure peuvent conduire à des problèmes tels que l’inconfort des passagers,une manoeuvrabilité réduite, voire même provoquer un comportement instable. Cette thèse apour objectif de développer de nouveaux modèles d’interaction fluide-structure, en mettant enoeuvre la théorie des systèmes Hamiltoniens à ports d’interaction (pHs). Le formalisme pHsfournit d’une part un cadre unifié pour la description des systèmes multi-physiques complexeset d’autre part une approche modulaire pour l’interconnexion des sous-systèmes grâce auxports d’interaction. Cette thèse s’intéresse aussi à la conception de contrôleurs à partir desmodèles pHs. Des modèles pHs sont proposés pour les équations de ballottement du liquide en partantdes équations de Saint Venant en 1D et 2D. L’originalité du travail est de donner des modèlespHs pour le ballottement dans des réservoirs en mouvement. Les ports d’interaction sont utiliséspour coupler la dynamique du ballottement à la dynamique d’une poutre contrôlée par desactionneurs piézo-électriques, celle-ci étant préalablement modélisée sous forme pHs. Aprèsl’écriture des équations aux dérivées partielles dans le formalisme pHs, une approximation endimension finie est obtenue en utilisant une méthode pseudo-spectrale géométrique qui conservela structure pHs du modèle continu au niveau discret. La thèse propose plusieurs extensionsde la méthode pseudo-spectrale géométrique, permettant la discrétisation des systèmesavec des opérateurs différentiels du second ordre d’une part et avec un opérateur d’entrée nonborné d’autre part. Des essais expérimentaux ont été effectués sur une structure constituéed’une poutre liée à un réservoir afin d’assurer la validité du modèle pHs du ballottementdu liquide couplé à la poutre flexible, et de valider la méthode pseudo-spectrale de semi-discrétisation.Le modèle pHs a finalement été utilisé pour concevoir un contrôleur basé surla passivité pour réduire les vibrations du système couplé. / This thesis is motivated by an aeronautical issue: the fuel sloshing in tanksof very flexible wings. The vibrations due to these coupled phenomena can lead to problemslike reduced passenger comfort and maneuverability, and even unstable behavior. Thisthesis aims at developing new models of fluid-structure interaction based on the theory ofport-Hamiltonian systems (pHs). The pHs formalism provides a unified framework for thedescription of complex multi-physics systems and a modular approach for the coupling ofsubsystems thanks to interconnection ports. Furthermore, the design of controllers using pHsmodels is also addressed. PHs models are proposed for the equations of liquid sloshing based on 1D and 2D SaintVenant equations and for the equations of structural dynamics. The originality of the workis to give pHs models of sloshing in moving containers. The interconnection ports are used tocouple the sloshing dynamics to the structural dynamics of a beam controlled by piezoelectricactuators. After writing the partial differential equations of the coupled system using thepHs formalism, a finite-dimensional approximation is obtained by using a geometric pseudospectralmethod that preserves the pHs structure of the infinite-dimensional model at thediscrete level. The thesis proposes several extensions of the geometric pseudo-spectral method,allowing the discretization of systems with second-order differential operators and with anunbounded input operator. Experimental tests on a structure made of a beam connected to atank were carried out to validate both the pHs model of liquid sloshing in moving containersand the pseudo-spectral semi-discretization method. The pHs model was finally used to designa passivity-based controller for reducing the vibrations of the coupled system.

Dans cette thèse, nous abordons le problème de la stabilisation des systèmes non linéaires. En particulier, nous nous concentrons sur les modèles où l'énergie joue un rôle fondamental. Ce cadre énergétique est adapté pour capturer les phénomènes de plusieurs domaines physiques tels que les systèmes mécaniques, les systèmes électriques, les systèmes hydrauliques, etc. Le point de départ des contrôleurs proposés sont les concepts de système passif, des sorties passives et des fonctions d'énergie (ou stockage). Dans ce travail, nous étudions deux classes de systèmes dynamiques, à savoir les Hamiltoniens à ports (PH) et les Euler-Lagrange (EL), qui conviennent pour représenter de nombreux processus physiques. Une première étape vers la construction des contrôleurs est de montrer la passivité des systèmes PH et la caractérisation de leurs sorties passives. Par la suite, nous explorons l'utilisation des différentes sorties passives dans deux techniques bien connues de contrôle par passivité (PBC), c'est-à-dire le contrôle par interconnexion (CbI) et l'équilibrage énergétique (EB), et nous comparons les résultats obtenus dans les deux approches. De plus, nous proposons une nouvelle méthodologie dans laquelle la loi de commande est composée d'un terme proportionnel (P), un terme intégral (I) et, éventuellement, un terme dérivatif (D) de la sortie passive. Dans cette stratégie, l'énergie du système en boucle fermée est façonnée sans qu'il soit nécessaire de résoudre des équations différentielles partielles (PDE). Nous analysons le scénario du régulateur PID à l'aide des différentes sorties passives précédemment caractérisées. Enfin, nous appliquons un schéma PID-PBC récemment proposé dans la littérature à un système mécanique complexe, à savoir un pendule inversé ultra flexible, représenté sous la forme d'un modèle contraint EL. La conception du contrôleur, la preuve de stabilité, ainsi que les simulations et les résultats expérimentaux sont présentés pour montrer l'applicabilité de cette technique aux systèmes physiques. / In this thesis we address the problem of stabilization of nonlinear systems. In particular, we focus on models where the energy plays a fundamental role. This energy-based framework is suitable to capture the phenomena of several physical domains, such as mechanical systems, electrical systems, hydraulic systems, etc. The starting point in the proposed controllers are the concepts of passive system, passive outputs and energy (storage) functions. In this work we study two classes of dynamical systems, namely port-Hamiltonian (PH) and Euler-Lagrange (EL), which are suitable to represent many physical processes. A first step towards the controller design is to show the passivity of the PH systems and the characterization of their passive outputs. Thereafter, we explore the use of the different passive outputs in two well-known passivity-based control (PBC) techniques, that is control by interconnection (CbI) and energy balancing (EB), and we compare the obtained results in both approaches. In addition, we propose a novel methodology in which the controller consists in a proportional (P), an integral (I) and, possibly, a derivative (D) term of the passive output. In this approach the energy of the closed-loop system is shaped without the necessity of solving partial differential equations (PDEs). We analyze the scenario of the PID controller using the different passive outputs previously characterized. Finally, we apply a PID-PBC scheme recently proposed in the literature to a complex mechanical system, namely an ultra flexible inverted pendulum, which is represented as a constrained EL model. The controller design, the stability proof, as well as simulations and experimental results are presented to show the applicability of this technique to physical systems.

Cette thèse porte sur la conception de commandes robustes pour les systèmes non linéaires, mettant l'accent sur les systèmes mécaniques. Des résultats concluants sont présentés pour résoudre deux situations très abordées dans la théorie du contrôle : 1) La stabilité des systèmes non linéaires perturbés ; 2) Le suivi global de trajectoire dans les systèmes mécaniques en ayant seulement connaissance de la position. Nous avons commencé par donner une méthode de conception des commandes robustes pour assurer une régulation de sortie non passive. En outre, si le système est perturbé (pas appariés), des preuves rigoureuses pour les rejeter sont fournies. Ce résultat est principalement inspiré d'un changement de coordonnées et de l'action intégrale dynamique. Si le scénario à traiter concerne des systèmes mécaniques avec des perturbations variant dans le temps, nous dotons le système de propriétés comme IISS (Integral Input- State Stable) et ISS (Input-State Stable). Ce résultat est obtenu en modifiant la procédure de conception de manière à rejeter les perturbations constantes (pas appariés). Cependant, en raison de la non-linéarité du système, les commandes qui en résultent ont une grande complexité. Pour le même problème, un deuxième et élégant résultat est donné au cas où un changement préalable de variable (impulsions) est réalisé. Finalement, une réponse convaincante au problème de suivi de trajectoire pour les systèmes mécaniques est donnée en tenant compte uniquement des informations de position. Nous résolvons ce problème en deux étapes. Premièrement, quelques modifications sont apportées à la preuve de stabilité d'un observateur de vitesse basée sur la théorie de l'invariance et l’Immersion récemment publié. Notez que ceci est un observateur satisfaisant la convergence exponentielle de vitesse dans les systèmes mécaniques. Deuxièmement et sur la base du changement de coordonnées (impulsions), un contrôleur de suivi avec stabilité exponentielle, tenant compte de la position et de la vitesse, est proposé. De telle sorte qu'avec la combinaison des deux résultats, le suivi de trajectoire exponentielle avec retour de position est donné. / This thesis focuses on the design of robust control for nonlinear systems, mainly on mechanical systems. The results presented are to two situations widely discussed in control theory: 1) The stability of nonlinear systems disturbed; 2) The global tracking trajectory in mechanical systems having only knowledge of the position. We started giving a design method of robust controls to ensure regulation on non-passive output. In addition, if the system is perturbed (constant unmatched), rigorous proof to its rejection is provided. This result is based mainly on change of coordinates and integral dynamic control. When the scenario to deal are mechanical systems with time-varying matched and unmatched, disturbance, the system is endowed with strong properties as IISS (Integral Input-State Stable) and ISS (Input-State Stable). This is achieved based on the design method to rejection of constant disturbances (unmatched). However, due to the nonlinearity of the system, the controllers have a high complexity. For the same problem, a second and elegant result is given making a initial change of coordinate on the momenta variable, such that the controller significantly simplifies, preserving the aforementioned robustness properties. Finally, a convincing answer to the problem of global exponential tracking of mechanical systems is given taking into account only the position information. We solve this problem in two steps. First, some slight variation is presented to the proof of stability of a speed observer based on Immersion and Invariance theory recently published. Note that this is a speed observer satisfying the exponential convergence speed in mechanical systems. Secondly, and based on the change of coordinates (momenta), a globally exponentially stable tracking controller with position and velocity known is proposed. The combination of both results give the first global exponential tracking controller of mechanical systems without velocity measurements.

Cette thèse porte sur l’étude des systèmes de lois de conservation couplés par une interface mobile. Un modèle dynamique d’un procédé d’extrusion obtenu à partir des bilans de masse, de taux d’humidité et d’énergie est proposé. Ce modèle exprime le transport de la matière et de la chaleur dans une extrudeuse par des systèmes d’équations hyperboliques définis sur deux domaines complémentaires variant dans le temps. L’évolution des domaines est dictée par une Equation aux Dérivées Ordinaires (EDO) issue du bilan de masse total dans une extrudeuse. Par le principe des applications contractantes l’existence et l’unicité de la solution pour cette classe de système sont prouvées. Le problème de stabilisation de l’interface mobile est aussi abordé en utilisation le formalisme des systèmes à retard. La méthode des caractéristiques permet de représenter le système composé des équations issues du bilan de masse par un système à retard sur l’entrée. Au moyen d’un contrôleur prédictif la position de l’interface est stabilisée autour d’un point équilibre. La dernière partie de ce travail est dédiée à l’étude des systèmes Hamiltoniens à ports frontière couplés par une interface mobile. Ces systèmes augmentés de variables couleur qui sont des fonctions caractéristiques du domaine peuvent s’exprimer comme des systèmes Hamiltoniens à ports frontière / This thesis is devoted to the analysis of Partial Differential Equations (PDEs) which are coupled through a moving interface. The motion of the interface obeys to an Ordinary Differential Equation (ODE) which arises from a conservation law. The first part of this thesis concerns the modelling of an extrusion process based on mass, moisture content and energy balances. These balances laws express heat and homogeneous material transport in an extruder by hyperbolic PDEs which are defined in complementary time-varying domains. The evolution of the coupled domains is given by an ODE which is derived from the conservation of mass in an extruder. In the second part of the manuscript, a mathematical analysis has been performed in order to prove the existence and the uniqueness of solution for such class of systems by mean of contraction mapping principle. The third part of the thesis concerns the transformation of an extrusion process mass balance equations into a particular input delay system framework using characteristics method. Then, the stabilization of the moving interface by a predictor-based controller has been proposed. Finally, an extension of the analysis of moving interface problems to a particular class of systems of conservations laws has been developed. Port-Hamiltonian formulation of systems of two conservation laws defined on two complementary time-varying intervals has been studied. It has been shown that the coupled system is a port-Hamiltonian system augmented with two variables being the characteristic functions of the two spatial domains

Cette thèse répond à une problématique de commande frontière d’une conduite acoustique dont l’actionnement est assuré par un réseau d’actionneurs/capteurs co-localisés constituant une peau active. Pour faire face au caractère intrinsèquement multiphysique de ce problème vibro-acoustique, nous avons choisi dans cette thèse d’employer une approche hamiltonienne à ports, approche structurée basée sur la représentation des échanges entre différents domaines énergétiques au sein d’un système et entre différents systèmes. Nous avons proposé une modélisation hamiltonienne à ports de l’équation d’onde interconnectée à la frontière au système d’actionnement distribué, correspondant à une formulation 2D du problème physique. Nous avons développé une méthode de discrétisation spatiale basée sur l’utilisation de différences finies sur plusieurs grilles en quinconce qui préserve la structure hamiltonienne à ports de l’équation d’onde. Cette méthode permet en outre d’interconnecter facilement le système discrétisé avec d’autres sous-systèmes, dans le but de mettre en place un actionnement par exemple. Son principal avantage sur d’autres méthodes préservatives de structure réside dans sa simplicité de mise en œuvre qui découle de l’utilisation de différences finies. Concernant la commande du système vibro-acoustique, nous avons proposé une méthode de synthèse de loi de commande distribuée pour les systèmes régis par deux lois de conservation en 1D. L’originalité de cette méthode réside en le fait qu’elle repose sur le calcul d’invariants structuraux (fonctions de Casimir) exploités afin de modifier la structure du système en boucle fermée. Les conditions d’application sur un système 2D sont étudiées et des résultats numériques valident les lois de commande synthétisées. / This thesis deals with the boundary control of an acoustic by a network of co-localised sensors/actuators which constitutes a smart skin. In order to cope with this multiphysical problem, we chose to place our study in the framework of port-Hamiltonian systems, a structured approach based on the representation of energy exchanges between different energy domains between different systems of subsystems. We proposed a port-Hamiltonian model of the wave equation interconnected through its boundary to the distributed actuation system, which corresponds to a 2D formulation of the physical problem. We developed a spatial discretization method based on the use of finite differences on several staggered grids that preserve the port-Hamiltonian structure of the wave equation. This method also permits to easily interconnect the discretized system with other subsystems, which is convenient for instance for control purposes. Its main advantage over other structure preserving methods is its simplicity of implementation which stems from the use of finite differences. In order to control the vibro-acoustic system, we proposed a control law synthesis method for systems governed by two conservation laws in 1D. The originality of this method lies in the fact that it relies on the computation of structural invariants (Casimir functions) exploited in order to modify the structure of the system in closed loop. The conditions of application of these laws on a 2D system are studied and numerical results validate the synthesized control laws.

L'objectif principal de la thèse est d'établir un modèle sous forme hamiltonienne à ports pour la dynamique du plasma dans les réacteurs de fusion de type tokamak, puis de démontrer le potentiel de cette approche pour aborder les problèmes d'intégration numérique et de commande non linéaire. Les bilans thermo-magnéto-hydrodynamiques, écrits sous forme hamiltonienne à ports à l'aide de structures Stokes-Dirac, conduisent à un modèle 3D “ multi-physique ” du plasma. Ensuite, un modèle 1D équivalent au modèle de diffusion résistive est obtenu en supposant les mêmes hypothèses d'équilibre quasi-statique et de symétries. Un schéma symplectique de réduction spatiale de ce modèle 1D qui préserve la structure du modèle et ses invariants est établi. Il ouvre la voie à des travaux ultérieurs de commande non linéaire fondés sur la structure géométrique d'interconnexion et les bilans du modèle. La commande IDA-PBC (Interconnection and Damping Assignment - Passivity Based Control) basée sur la passivité du modèle est d'abord synthétisée pour ce système en dimension finie. Finalement, une commande IDA-PBC associée avec la commande à la frontière est proposée pour le système en dimension infinie. Les controlleurs sont testés et validés avec les simulateurs des tokamak (METIS pour le Tore Supra de CEA/ Cadarache, et RAPTOR pour le TCV de l'EPFL Lausanne, Suisse). / The modelling and analysis of the plasma dynamics in tokamaks using the port-Hamiltonian approach is the main project purpose. Thermo-mMagnetohydrodynamics balances have been written in port-Hamiltonian form using Stokes-Dirac interconnection structures and 3D differential forms. A simplified 1D model for control has been derived using quasi-static and symmetry assumptions. It has been proved to be equivalent to a classical 1D control model: the resistive diffusion model for the poloidal magnetic flux. Then a geometric spatial integration scheme has been developped. It preserves both the symplecticity of the Dirac interconnection structure and physically conserved extensive quantities. This will allow coming works on energy-based approaches for the non linear control of the plasma dynamics.An Interconnection and Damping Assignment - Passivity Based Control (IDA-PBC) , the most general Port-Hamiltonian control, is chosen first to deal with the studied Tokamak system. It is based on a model made of the two coupled PDEs of resistive diffusion for the magnetic poloidal flux and of radial thermal diffusion. The used TMHD couplings are the Lorentz forces (with non-uniform resistivity) and the bootstrap current. The loop voltage at the plasma boundary, the total external current and the plasma heating power are considered as controller outputs. Due to the actuator constraints which imply to have a physically feasible current profile deposits, a feedforward control is used to ensure the compatibility with the actuator physical capability. Then, the IDA-PBC controllers, both finite-dimensional and infinite-dimensional, are designed to improve the system stabilization and convergence speed. The proposed works are validated against the simulation data obtained from the Tore-Supra WEST (CEA/Cadarache, France) test case and from RAPTOR code for the TCV real-time control system (CRPP/ EPFL, Lausanne, Switzerland).

Die Anwenung von Operator-Splitting-Methoden auf gewöhnliche Differentialgleichungen ist gut etabliert. Für Differential-algebraische Gleichungen und partielle Differential-algebraische Gleichungen unterliegt sie jedoch vielen Einschränkungen aufgrund des Vorhandenseins von Nebenbedingungen. Die räumliche Diskretisierung reduziert PDAEs und lenkt unseren Fokus auf das Konzept der DAEs. Um eine reibungslose Übertragung des Operator-Splittings von ODEs auf DAEs durchzuführen, ist es wichtig, eine geeignete entkoppelte Struktur für das gewünschte Differential-algebraische System zu haben. In dieser Arbeit betrachten wir ein Modell, das partielle Differentialgleichungen für elektromagnetische Bauelemente - modelliert durch die Maxwell-Gleichungen - mit Differential-algebraischen Gleichungen koppelt, die die elementaren Schaltungselemente beschreiben. Nach der räumlichen Diskretisierung der klassischen Formulierung der Maxwell-Gleichungen mit Hilfe der finiten Integrationstechnik formulieren wir das resultierende gekoppelte System als Differential-algebraische Gleichung. Um eine geeignete Entkopplung zu bekommen, verwenden wir den zweigorientierten Loop-Cutset-Ansatz für die Schaltungsmodellierung. Daraus folgt, dass wir in der Lage sind, eine geeignete Operatorzerlegung so zu konstruieren, dass wir eine natürliche topologisch entkoppelte Port-Hamiltonsche DAE-Struktur erhalten. Wir schlagen einen Operator-Splitting-Ansatz für die Schaltungs-DAEs und gekoppelten Feld-Schaltungs-DAEs in entkoppelter Form vor und analysieren seine numerischen Eigenschaften. Darüber hinaus nutzen wir das Hamiltonsche Verhalten der inhärenten gewöhnlichen Differentialgleichung durch die Verwendung expliziter und energieerhaltender Zeitintegrations-methoden. Schließlich führen wir numerische Tests, um das mathematische Modell zu illustrieren und die Konvergenzergebnisse für das vorgeschlagene DAE-Operator-Splitting zu demonstrieren. / Le equazioni algebriche differenziali e algebriche alle derivate parziali hanno avuto un enorme successo come modelli di sistemi dinamici vincolati. Nella modellazione matem- atica, spesso si desidera catturare diversi aspetti di una situazione come le leggi di conservazione della fisica, il trasporto convettivo o la diffusione. Queste aspetti si riflettono nel sistema di equazioni del modello come operatori diversi. La tecnica dell’Operator Splitting si è rivelata una strategia di successo per affrontare problemi così complicati. L’applicazione dei metodi di Operator Splitting alle equazioni differenziali ordinarie (ODE) è ormai una tecnologia ben consolidata. Tuttavia, per equazioni algebriche differenziali (DAE) e algebriche differenziali parziali (PDAE), l’approccio è soggetto a molte restrizioni dovute alla presenza di vincoli e alla proprietà di indice. La discretizzazione spaziale riduce le PDAE e indirizza la nostra attenzione al concetto di DAE. Le DAE emergono in problemi dinamici vincolati come circuiti elettrici o reti di trasporto di energia. Al fine di generalizzare agevolmente la tecnica dell’Operator Splitting dalle ODE alle DAE, è importante avere una struttura disaccoppiata adeguata per il sistema algebrico differenziale desiderato. In questa tesi, consideriamo un modello che accoppia equazioni differenziali alle derivate parziali per dispositivi elettromagnetici -modellati dalle equazioni di Maxwell- con equazioni algebriche differenziali che descrivono gli elementi base del circuito. Dopo aver discretizzato spazialmente la formulazione classica delle equazioni di Maxwell usando la tecnica di integrazione finita, formuliamo il sistema accoppiato risultante come una equazione algebrica differenziale. Interpretando il dispositivo elettromagnetico come un elemento capacitivo, l’indice dell’intero sistema di circuito e campo accoppiato può essere specificato utilizzando le proprietà topologiche del circuito e non supera il valore di due. Per eseguire un disaccoppiamento appropriato, utilizziamo l’approccio loop-cutset per la modellazione dei circuiti. In tal modo siamo in grado di costruire una opportuna decomposizione dell’operatore tale da ottenere una naturale struttura disaccoppiata port-Hamiltonian DAE. Proponiamo un approccio di suddivisione dell’operatore per i DAE a circuito disaccoppiato e a circuito di campo accoppiato utilizzando gli algoritmi di divisione Lie-Trotter e Strang e per analizzare le proprietà numeriche di questi sistemi. Inoltre, sfruttiamo il comportamento hamiltoniano del sistema di equazioni differenziali ordinarie mediante l’utilizzo di metodi di integrazione temporale con esatta conservazione dell’energia. Poggiando sull’analisi di convergenza del metodo di suddivisione dell’operatore ODE, deriviamo i risultati di convergenza per l’approccio proposto che dipendono dall’indice delsistema e quindi dalla sua struttura topologica. Infine, eseguiamo prove numeriche di sistemi circuitali, nonchè sistemi accoppiati a circuito di campo, per testare il modello matematico e dimostrare i risultati di convergenza per la proposta Operator Splitting DAE. / The application of operator splitting methods to ordinary differential equations (ODEs) is well established. However, for differential-algebraic equations (DAEs) and partial differential-algebraic equations (PDAEs), it is subjected to many restrictions due to the presence of constraints. In constrained dynamical problems as electrical circuits or energy transport networks, DAEs arise. In order to perform a smooth transfer of the operator splitting from ODEs to DAEs, it is important to have a suitable decoupled structure for the desired differential-algebraic system. In this thesis, we consider a model which couples partial differential equations for electro- magnetic devices -modeled by Maxwell’s equations- with differential-algebraic equations describing the basic circuit elements. After spatially discretizing the classical formulation of Maxwell’s equations using the finite integration technique, we formulate the resulting coupled system as a differential-algebraic equation. To perform an appropriate decoupling, we use the branch oriented loop-cutset approach for circuit modeling. It follows that we are able to construct a suitable operator decomposition such that we obtain a natural topologically decoupled port-Hamiltonian DAE structure. We propose an operator splitting approach for the decoupled circuit and coupled field- circuit DAEs using the Lie-Trotter and Strang splitting algorithms and analyze its numerical properties. Furthermore, we exploit the Hamiltonian behavior of the system’s inherent ordinary differential equation by the utilization of explicit and energy-preserving time integration methods. Based on the convergence analysis of the ODE operator splitting method, we derive convergence results for the proposed approach that depends on the index of the system and thus on its topological structure. Finally, we perform numerical tests, to underline the mathematical model and to demonstrate the convergence results for the proposed DAE operator splitting.