Diese Arbeit befasst sich mit Differential-algebraischen Gleichungen (DAEs). DAEs spielen eine wichtige Rolle in der Modellierung, der Simulation und der Optimierung von Netzwerken und gekoppelten Problemen in vielen Anwendungsgebieten. Es werden in Bezug auf die Modellierung und die numerische Simulation von DAEs bereits bestehende Ergebnisse diskutiert und erweitert. Des Weiteren wird die globale eindeutige Lösbarkeit und die Sensitivität der Lösungen mit Hinsicht auf Störungen der DAEs untersucht. Häufig wird die Modellierung von multiphysikalischen Anwendungen durch die Kopplung mehrerer einzelner DAE Systeme realisiert. Diese Herangehensweise kann hochdimensionale DAEs erzeugen, welche aufgrund von Instabilitäten nicht von klassischen numerischen Methoden, simuliert werden können. Angesichts dieser Herausforderungen werden drei Ziele formuliert: Erstens wird ein globales Lösungstheorem formuliert und bewiesen, welches auf gekoppelte Systeme angewandt werden kann, um deren Kopplungsansatz mathematisch zu rechtfertigen. Zweitens werden numerische Methoden vorgestellt, welche unter wesentlich schwächeren Strukturannahmen stabil sind und sich daher für die Simulation von gekoppelten Systemen eignen. Drittens wird eine Strategie präsentiert, die es ermöglicht, explizite Methoden auf gekoppelte Systeme anzuwenden. Um diese Ziele zu erreichen, braucht man ein Entkopplungsverfahren für DAEs, welches die folgenden drei Eigenschaften erfüllt: Die Komplexität des Entkopplungsverfahrens sollte nicht die Komplexität der DAE überschreiten. Das Entkopplungsverfahren sollte Eigenschaften wie Symmetrie, Monotonie und positive Definitheit erhalten. Das Entkopplungsverfahren sollte durch einen Schritt-für-Schritt Ansatz mit unabhängigen Schritten realisiert werden. Sowohl das Konzept des Tractability Index als auch das des Strangeness Index liefert kein solches Entkopplungsverfahren. Daher wird hier ein neues Index Konzept eingeführt, das diesen Anforderungen entspricht. / This thesis addresses differential-algebraic equations (DAEs). They play an important role in the modeling, simulation and optimization of networks and coupled problems in various applications. The main application in this thesis are coupled problems in electric circuit simulation. We discuss and extend existing results regarding the modeling and numerical simulation of DAEs. Furthermore, we investigate the global unique solvability and the sensitivity of solutions with respect to perturbations of DAEs. Nowadays the modeling of multi-physical applications is often realized by coupling systems of DAEs together with the help of additional coupling terms. This approach may yield high dimensional DAEs which cannot be simulated, due to instabilities, by standard numerical methods. Regarding these challenges we formulate three objectives: First we provide a global solvability theorem which can be applied to coupled systems to mathematically justify their coupling approach. Second we introduce numerical methods which are stable without needing any structural assumptions. Third we provide a way to apply explicit methods to coupled systems to be able to handle the size of the coupled systems by parallelizing the algorithms. To achieve these objectives, we need a decoupling procedure which fulfills the following three properties: The complexity of the decoupling procedure has to reflect the complexity of the DAE, i.e. the decoupling procedure should be state-independent if possible. The decoupling procedure should preserve properties like symmetry, monotonicity and positive definiteness. The decoupling procedure should be realized by a step-by-step approach with independent stages. Both the Tractability Index concept and the Strangeness Index concept do not provide such a decoupling procedure. For this reason we introduce a new index concept which provides such a decoupling procedure.

Industrial robot manipulators are general-purpose machines used for industrial automation in order to increase productivity, flexibility, and quality. Other reasons for using industrial robots are cost saving, and elimination of heavy and health-hazardous work. Robot motion control is a key competence for robot manufacturers, and the current development is focused on increasing the robot performance, reducing the robot cost, improving safety, and introducing new functionalities. Therefore, there is a need to continuously improve the models and control methods in order to fulfil all conflicting requirements, such as increased performance for a robot with lower weight, and thus lower mechanical stiffness and more complicated vibration modes. One reason for this development of the robot mechanical structure is of course cost-reduction, but other benefits are lower power consumption, improved dexterity, safety issues, and low environmental impact. This thesis deals with three different aspects of modeling and control of flexible, i.e., elastic, manipulators. For an accurate description of a modern industrial manipulator, the traditional flexible joint model, described in literature, is not sufficient. An improved model where the elasticity is described by a number of localized multidimensional spring-damper pairs is therefore proposed. This model is called the extended flexible joint model. This work describes identification, feedforward control, and feedback control, using this model. The proposed identification method is a frequency-domain non-linear gray-box method, which is evaluated by the identification of a modern six-axes robot manipulator. The identified model gives a good description of the global behavior of this robot. The inverse dynamics control problem is discussed, and a solution methodology is proposed. This methodology is based on a differential algebraic equation (DAE) formulation of the problem. Feedforward control of a two-axes manipulator is then studied using this DAE approach. Finally, a benchmark problem for robust feedback control of a single-axis extended flexible joint model is presented and some proposed solutions are analyzed.

<p>Industrial robot manipulators are general-purpose machines used for industrial automation in order to increase productivity, flexibility, and quality. Other reasons for using industrial robots are cost saving, and elimination of heavy and health-hazardous work. Robot motion control is a key competence for robot manufacturers, and the current development is focused on increasing the robot performance, reducing the robot cost, improving safety, and introducing new functionalities. Therefore, there is a need to continuously improve the models and control methods in order to fulfil all conflicting requirements, such as increased performance for a robot with lower weight, and thus lower mechanical stiffness and more complicated vibration modes. One reason for this development of the robot mechanical structure is of course cost-reduction, but other benefits are lower power consumption, improved dexterity, safety issues, and low environmental impact.</p><p>This thesis deals with three different aspects of modeling and control of flexible, i.e., elastic, manipulators. For an accurate description of a modern industrial manipulator, the traditional flexible joint model, described in literature, is not sufficient. An improved model where the elasticity is described by a number of localized multidimensional spring-damper pairs is therefore proposed. This model is called the extended flexible joint model. This work describes identification, feedforward control, and feedback control, using this model.</p><p>The proposed identification method is a frequency-domain non-linear gray-box method, which is evaluated by the identification of a modern six-axes robot manipulator. The identified model gives a good description of the global behavior of this robot.</p><p>The inverse dynamics control problem is discussed, and a solution methodology is proposed. This methodology is based on a differential algebraic equation (DAE) formulation of the problem. Feedforward control of a two-axes manipulator is then studied using this DAE approach.</p><p>Finally, a benchmark problem for robust feedback control of a single-axis extended flexible joint model is presented and some proposed solutions are analyzed.</p>

[pt] Esse trabalho apresenta uma abordagem baseada em programação matemática para a solução de problemas de valor inicial de contorno constitutivo elastoplástico. Considerando a plasticidade associativa, as equações constitutivas locais, em sua forma discreta, são formuladas como problemas de otimização cônica. Especificamente, é demonstrado que métodos implícitos de retorno a superfície para os critérios mais conhecidos da literatura, como o de Rankine, von Mises, Tresca, Drucker-Prager e Mohr Coulomb, podem ser expressos como problemas de otimização cônica de segunda ordem e semidefinida. Além disso, um novo método numérico para a determinação do operador elastoplástico consistente, baseado na derivada paramétrica de primeira ordem das soluções ótimas, é proposto. / [en] This work presents a mathematical programming approach for elastoplastic constitutive initial boundary value problems. Considering associative plasticity, the local discrete constitutive equations are formulated as conic programs. Specifically, it is demonstrated that implicit return-mapping schemes for well-known yield criteria, such as the Rankine, von Mises, Tresca, Drucker-Prager, and Mohr-Coulomb criteria, can be expressed as secondorder and semidefinite conic programs. Additionally, a novel scheme for the numerical evaluation of the consistent elastoplastic tangent operator, based on a first-order parameter derivative of the optimal solutions, is proposed.

This research effort develops a comprehensive computational framework to support the parametric optimal design of uncertain dynamical systems. Uncertainty comes from various sources, such as: system parameters, initial conditions, sensor and actuator noise, and external forcing. Treatment of uncertainty in design is of paramount practical importance because all real-life systems are affected by it; not accounting for uncertainty may result in poor robustness, sub-optimal performance and higher manufacturing costs. Contemporary methods for the quantification of uncertainty in dynamical systems are computationally intensive which, so far, have made a robust design optimization methodology prohibitive. Some existing algorithms address uncertainty in sensors and actuators during an optimal design; however, a comprehensive design framework that can treat all kinds of uncertainty with diverse distribution characteristics in a unified way is currently unavailable. The computational framework uses Generalized Polynomial Chaos methodology to quantify the effects of various sources of uncertainty found in dynamical systems; a Least-Squares Collocation Method is used to solve the corresponding uncertain differential equations. This technique is significantly faster computationally than traditional sampling methods and makes the construction of a parametric optimal design framework for uncertain systems feasible. The novel framework allows to directly treat uncertainty in the parametric optimal design process. Specifically, the following design problems are addressed: motion planning of fully-actuated and under-actuated systems; multi-objective robust design optimization; and optimal uncertainty apportionment concurrently with robust design optimization. The framework advances the state-of-the-art and enables engineers to produce more robust and optimally performing designs at an optimal manufacturing cost. / Ph. D.

La recherche et les questions abordées dans cette thèse sont de type inverse : la reconstitution d'une source ponctuelle ou la complétion d'une donnée à la limite inconnue à l'extrémité du domaine dans les modèles paraboliques de transport de contaminants. La modélisation mathématique des problèmes de pollution des eaux fait intervenir deux traceurs, l'oxygène dissous (OD) et la demande biochimique en oxygène (DBO) qui est la quantité d'oxygène nécessaire à la biodégradation de la matière organique. En effet, au cours des procédés d'autoépuration, certaines bactéries aérobies jouent un rôle principal. Ces micro-organismes décomposent les matières organiques polluantes en utilisant l'oxygène dissous dans le milieu. Afin de compenser ces données manquantes, les champs, solutions du problème, sont observés directement ou indirectement. Les problèmes inverses qui en résultent sont quasi certainement mal-posés voire même sévèrement mal-posés pour la plupart. Dans cette thèse, nous proposons justement une analyse aussi poussée que possible sur la question de l'identifiabilité pour les deux problèmes inverses décrits ci-dessus. Nous avons démontré un résultat d'unicité pour des sources fixes dans le cas d'observations décalées. La réalité pour l'observation est nuancée et l'idéal n'est pas acquis ; des mesures directes sur la DBO sont difficiles à obtenir. En revanche collecter des données sur l'OD est possible en temps réel et avec un faible coût. La DBO est donc observée de façon indirecte, grâce au couplage dans le système de Streeter et Phelps, l'information passe de l'OD à la DBO. Pour ce problème aussi, nous avons produit un résultat d'unicité pour la reconstruction de la source ou puits ponctuel qui serait présent dans l'équation de transport sur l'OD. Nous avons ensuite examiné des questions annexes à l'identifiabilité telles que le degré d'instabilité des équations à résoudre. De ce type d'informations dépendent le comportement des méthodes numériques et des algorithmes de calcul à utiliser. / The research and the questions approached on this thesis are inverse type : the reconstruction of point-wise source or the data completion problem in parabolic models of transport of contaminants. The mathematical modelling of the problems of water pollution includes two tracers, the dissolved oxygen (DO) and the biochemical demand in oxygen (BDO) which is the quantity of oxygen necessary for the biodegradation of organic matter. Indeed, during the biodegradation process, aerobic bacteria play a leading part. These micro-organisms decompose polluting organic matters by using the dissolved oxygen in the middle. To compensate these missing data, fields, solutions of the problem, are observed directly or indirectly. The resulting inverse problems are ill-posed. Their mathematical study rises big complications and their numerical treatment isn't easy. We demonstrated a uniqueness result for fixed sources in the case of moved observations. The reality for the observation is qualified and the ideal is not acquired; direct measures on the BOD are difficult to obtain. On the Other hand to collect data on the DO is possible in real time With a moderate cost. The BOD is thus observed in indirect way, thanks to the coupling in the system of Streeter and Phelps, the information passes from the DO to the BOD. For this problem, we produced a uniqueness result for the reconstruction of source. Then, we examined the degree of instability of the equation to be solved. The behaviour of numerical methods depend on this type of information.

<p>A time-discrete formulation of the variational principle of mechanics is used to provide a consistent theoretical framework for the construction and analysis of low order integration methods. These are applied to mechanical systems subject to mixed constraints and dry frictional contacts and impacts---machines. The framework includes physics motivated constraint regularization and stabilization schemes. This is done by adding potential energy and Rayleigh dissipation terms in the Lagrangian formulation used throughout. These terms explicitly depend on the value of the Lagrange multipliers enforcing constraints. Having finite energy, the multipliers are thus massless ghost particles. The main numerical stepping method produced with the framework is called SPOOK.</p><p>Variational integrators preserve physical invariants globally, exactly in some cases, approximately but within fixed global bounds for others. This allows to product realistic physical trajectories even with the low order methods. These are needed in the solution of nonsmooth problems such as dry frictional contacts and in addition, they are computationally inexpensive. The combination of strong stability, low order, and the global preservation of invariants allows for large integration time steps, but without loosing accuracy on the important and visible physical quantities. SPOOK is thus well-suited for interactive simulations, such as those commonly used in virtual environment applications, because it is fast, stable, and faithful to the physics.</p><p>New results include a stable discretization of highly oscillatory terms of constraint regularization; a linearly stable constraint stabilization scheme based on ghost potential and Rayleigh dissipation terms; a single-step, strictly dissipative, approximate impact model; a quasi-linear complementarity formulation of dry friction that is isotropic and solvable for any nonnegative value of friction coefficients; an analysis of a splitting scheme to solve frictional contact complementarity problems; a stable, quaternion-based rigid body stepping scheme and a stable linear approximation thereof. SPOOK includes all these elements. It is linearly implicit and linearly stable, it requires the solution of either one linear system of equations of one mixed linear complementarity problem per regular time step, and two of the same when an impact condition is detected. The changes in energy caused by constraints, impacts, and dry friction, are all shown to be strictly dissipative in comparison with the free system. Since all regularization and stabilization parameters are introduced in the physics, they map directly onto physical properties and thus allow modeling of a variety of phenomena, such as constraint compliance, for instance.</p><p>Tutorial material is included for continuous and discrete-time analytic mechanics, quaternion algebra, complementarity problems, rigid body dynamics, constraint kinematics, and special topics in numerical linear algebra needed in the solution of the stepping equations of SPOOK.</p><p>The qualitative and quantitative aspects of SPOOK are demonstrated by comparison with a variety of standard techniques on well known test cases which are analyzed in details. SPOOK compares favorably for all these examples. In particular, it handles ill-posed and degenerate problems seamlessly and systematically. An implementation suitable for large scale performance and accuracy testing is left for future work.</p>

A time-discrete formulation of the variational principle of mechanics is used to provide a consistent theoretical framework for the construction and analysis of low order integration methods. These are applied to mechanical systems subject to mixed constraints and dry frictional contacts and impacts---machines. The framework includes physics motivated constraint regularization and stabilization schemes. This is done by adding potential energy and Rayleigh dissipation terms in the Lagrangian formulation used throughout. These terms explicitly depend on the value of the Lagrange multipliers enforcing constraints. Having finite energy, the multipliers are thus massless ghost particles. The main numerical stepping method produced with the framework is called SPOOK. Variational integrators preserve physical invariants globally, exactly in some cases, approximately but within fixed global bounds for others. This allows to product realistic physical trajectories even with the low order methods. These are needed in the solution of nonsmooth problems such as dry frictional contacts and in addition, they are computationally inexpensive. The combination of strong stability, low order, and the global preservation of invariants allows for large integration time steps, but without loosing accuracy on the important and visible physical quantities. SPOOK is thus well-suited for interactive simulations, such as those commonly used in virtual environment applications, because it is fast, stable, and faithful to the physics. New results include a stable discretization of highly oscillatory terms of constraint regularization; a linearly stable constraint stabilization scheme based on ghost potential and Rayleigh dissipation terms; a single-step, strictly dissipative, approximate impact model; a quasi-linear complementarity formulation of dry friction that is isotropic and solvable for any nonnegative value of friction coefficients; an analysis of a splitting scheme to solve frictional contact complementarity problems; a stable, quaternion-based rigid body stepping scheme and a stable linear approximation thereof. SPOOK includes all these elements. It is linearly implicit and linearly stable, it requires the solution of either one linear system of equations of one mixed linear complementarity problem per regular time step, and two of the same when an impact condition is detected. The changes in energy caused by constraints, impacts, and dry friction, are all shown to be strictly dissipative in comparison with the free system. Since all regularization and stabilization parameters are introduced in the physics, they map directly onto physical properties and thus allow modeling of a variety of phenomena, such as constraint compliance, for instance. Tutorial material is included for continuous and discrete-time analytic mechanics, quaternion algebra, complementarity problems, rigid body dynamics, constraint kinematics, and special topics in numerical linear algebra needed in the solution of the stepping equations of SPOOK. The qualitative and quantitative aspects of SPOOK are demonstrated by comparison with a variety of standard techniques on well known test cases which are analyzed in details. SPOOK compares favorably for all these examples. In particular, it handles ill-posed and degenerate problems seamlessly and systematically. An implementation suitable for large scale performance and accuracy testing is left for future work.

Cette thèse aborde les formalismes pour la modélisation multi-physique en support au cycle en V deconception. Ce travail a été réalisé dans le cadre du projet ANR–MoCoSyMec, selon la méthodologie duprototypage virtuel fonctionnel (PVF) et illustré sur des systèmes électromagnétiques.Nous nous sommes principalement intéressés au langage VHDL-AMS, en tant que support aux différentsniveaux de modélisation apparaissant dans le cycle en V de conception. Cela nous a conduits à traiter laportabilité et l’interopérabilité en VHDL-AMS de diverses méthodes et outils de modélisation. Nous avonsproposé et validé, via le formalisme des composants logiciels ICAr, des solutions aux limites de l’utilisation deVHDL-AMS pour modéliser certains phénomènes physiques reposants sur des calculs numériques.Nous avons étendu la norme ICAr pour supporter des modèles dynamiques décrits par des équationsdifférentielles algébriques (DAE) ; et pour des besoins de co-simulation, nous pouvons également y associer unsolveur. Ces développements sont désormais capitalisés dans le framework CADES.Enfin, nous avons proposé une architecture pour le portage de modèles d’un formalisme à un autre. Elle a étédéfinie et mise en oeuvre plus particulièrement pour des modèles magnétiques réluctants (Reluctool) et desMEMS magnétiques (MacMMems) vers le VHDL-AMS.Ces formalismes et méthodologies sont mis en oeuvre autour du PVF d’un contacteur électromagnétique. / This PhD report deals with modeling formalisms for multi-physical systems in the design V- cycle. Thiswork was carried out within the French ANR-MoCoSyMec project, according to the methodology of functionalvirtual prototyping (PVF) and illustrated with electromagnetical systems.The work focuses on the VHDL-AMS modeling language, as a support for several modeling levels appearingin the design V-cycle. In this work, the portability and interoperability problems have been studied, usingVHDL-AMS, for various modeling methods and tools. Solutions have been proposed and validated for use limitsof VHDL-AMS language, specifically for the modeling of some physical phenomena using numericalcomputations, through the software component formalism called ICAr.The ICAr software component standard has been extended to support dynamic models described throughdifferential algebraic equations (DAE). It has also been extended for co-simulation purposes in which a solver isassociated to the dynamic model inside the ICAr component. These developed solutions are now available in theframework CADES.Finally, architecture has been proposed for the transforming of models from a professional formalism intoanother, specifically into VHDL-AMS. It has been designed and implemented for reluctant magnetic models(RelucTool) and magnetic MEMS (MacMMems).These formalisms and methodologies are implemented around the functional virtual prototyping (PVF) of anelectromagnetic contactor.

Chaque année, plus de 5 millions de personne à travers le monde deviennent hémiplégiques suite à un accident vasculaire cérébral. Ce soudain déficit neurologique conduit bien souvent à une perte partielle ou totale de la station debout et/ou à la perte de la capacité de déambulation. Dans l’optique de proposer de nouvelles solutions d’assistance situées entre le fauteuil roulant et le déambulateur, cette thèse s’inscrit dans le cadre du projet ANR TECSAN VHIPOD « véhicule individuel de transport en station debout auto-équilibrée pour personnes handicapées avec aide à la verticalisation ». Dans ce contexte, ces travaux de recherche apportent des éléments de réponse à deux problématiques fondamentales du projet : l’assistance au passage assis-debout (PAD) des personnes hémiplégiques et le déplacement à l’aide d’un véhicule auto-équilibré à deux roues. Ces problématiques sont abordées du point de vue de la robotique avec comme question centrale : peut-on utiliser l’approche Takagi-Sugeno (TS) pour la synthèse d’une commande ? Dans un premier temps, la problématique de mobilité des personnes handicapées a été traitée sur la base d’une solution de type gyropode. Des lois de commande basées sur les approches TS standard et descripteur ont été proposées afin d’étudier la stabilisation des gyropodes dans des situations particulières telles que le déplacement sur un terrain en pente ou le franchissement de petites marches. Les résultats obtenus ont non seulement permis d’aboutir à un concept potentiellement capable de franchir des obstacles, mais ils ont également permis de souligner la principale difficulté liée à l’applicabilité de l’approche TS en raison du conservatisme des conditions LMIs (inégalités matricielles linéaires). Dans un second temps, un banc d’assistance au PAD à architecture parallèle a été conçu. Ce type de manipulateur constitué de multiples boucles cinématiques présente un modèle dynamique très complexe (habituellement donné sous forme d’équations différentielles ordinaires). L’application de lois de commande basées sur l’approche TS est souvent vouée à l’échec compte tenu du grand nombre de non-linéarités dans le modèle. Afin de remédier à ce problème, une nouvelle approche de modélisation a été proposée. À partir d’un jeu de coordonnées bien particulier, le principe des puissances virtuelles est utilisé pour générer un modèle dynamique sous forme d’équations algébro-différentielles (DAEs). Cette approche permet d’aboutir à un modèle quasi-LPV où les seuls paramètres variants représentent les multiplicateurs de Lagrange issus de la modélisation DAE. Les résultats obtenus ont été validés en simulation sur un robot parallèle à 2 degrés de liberté (ddl) puis sur un robot parallèle à 3 ddl développé pour l’assistance au PAD. / Every year more than 5 million people worldwide become hemiplegic as a direct consequence of stroke. This neurological deficiency, often leads to a partial or a total loss of standing up abilities and /or ambulation skills. In order to propose new supporting solutions lying between the wheelchair and the walker, this thesis comes within the ANR TECSAN project named VHIPOD “self-balanced transporter for disabled persons with sit-to-stand function”. In this context, this research provides some answers for two key issues of the project : the sit-to-stand assistance (STS) of hemiplegic people and their mobility through a two wheeled self-balanced solution. These issues are addressed from a robotic point of view while focusing on a key question : are we able to extend the use of Takagi-Sugeno approach (TS) to the control of complex systems ? Firstly, the issue of mobility of disabled persons was treated on the basis of a self-balanced solution. Control laws based on the standard and descriptor TS approaches have been proposed for the stabilization of gyropod in particular situations such as moving along a slope or crossing small steps. The results have led to the design a two-wheeled transporter which is potentially able to deal with the steps. On the other hand, these results have also highlighted the main challenge related to the use of TS approach such as the conservatisms of the LMIs constraints (Linear Matrix Inequalities). In a second time, a test bench for the STS assistance based on parallel kinematic manipulator (PKM) was designed. This kind of manipulator characterized by several closed kinematic chains often presents a complex dynamical model (given as a set of ordinary differential equations, ODEs). The application of control laws based on the TS approach is often doomed to failure given the large number of non-linear terms in the model. To overcome this problem, a new modeling approach was proposed. From a particular set of coordinates, the principle of virtual power was used to generate a dynamical model based on the differential algebraic equations (DAEs). This approach leads to a quasi-LPV model where the only varying parameters are the Lagrange multipliers derived from the constraint equations of the DAE model. The results were validated on simulation through a 2-DOF (degrees of freedom) parallel robot (Biglide) and a 3-DOF manipulator (Triglide) designed for the STS assistance.