Avec le développement des ordinateurs, la résolution numérique des équations de la physique est devenue un outil de calcul puissant pour établir des prédictions physiques. Mais le recours au calcul numérique entraîne des changements plus profonds pour les théories physiques. Le but de cette thèse est de montrer que le calcul numérique sur machine conduit à une véritable reformulation des théories physiques. Les lois et les principes fondamentaux formulés à l'aide d'équations différentielles sont reformulés de manière discrète. Pour cela, je me concentre sur le cas d'une théorie physique: la mécanique classique. Je montre que depuis les années 1980 une mécanique discrète a été développée. J'analyse cette approche et j'examine en particulier ce qu'elle nous apprend sur la représentation du temps comme continu dans les théories physiques. Dans une première partie, j'examine la résolution numérique sur machine en tant qu'outil pour la prédiction quantitative en physique. Je montre la nécessité pour les scientifiques d'y recourir et je propose une analyse des concepts fondamentaux de ce type de résolution. Dans une deuxième partie, j'examine dans quelle mesure le calcul numérique est un élément constitutif des théories physiques. Je défends la thèse selon laquelle la mécanique discrète est une nouvelle théorie du mouvement classique. Dans une troisième partie, je soutiens une thèse sur la représentation du temps comme continu dans les théories physiques. C'est une représentation dont les scientifiques peuvent se passer. J'examine ensuite en quel sens la représentation traditionnelle du temps comme continu est plus simple que la représentation discrète. / The numerical computation of the solutions of equations in physical theories enables scientists to make powerful predictions. But numerical computation also challenges physical theories in a more fundamental way. The aim of this dissertation is to show how numerical computation leads to a reformulation of physical theories. Fundamental laws and first principles usually formulated with differential equations are reformulated with discrete equations. To fulfill this goal, I focus on the case of classical mechanics. I study a discrete approach called discrete mechanics developed since the 1980's and I discuss its consequences on the usual continuous representation of time in physics. First, I study numerical computation as a means to make predictions in physics. The fundamental concepts of exact and numerical computations of differential equations are discussed. In the second part, I examine how numerical computation changes the fundamental principles of physical theories. I claim that discrete mechanics has to be considered as a new theory of classical motion. In a third part, I investigate the consequences of discrete mechanics on the continuous representation of time in physics. I claim that physicists do not have to necessarily represent time as continuous in their theories. The discrete representation is another possible choice. Finally, I compare the continuous representation of time and the discrete one according to criteria of simplicity.

The first and most imperative step when designing a biologically inspired robot is to identify the underlying mechanics of the system or animal of interest. It is most common, perhaps, that this process generates a set of coupled nonlinear ordinary or partial differential equations. For this class of systems, the models derived from morphology of the skeleton are usually very high dimensional, nonlinear, and complex. This is particularly true if joint and link flexibility are included in the model. In addition to complexities that arise from morphology of the animal, some of the external forces that influence the dynamics of animal motion are very hard to model. A very well-established example of these forces is the unsteady aerodynamic forces applied to the wings and the body of insects, birds, and bats. These forces result from the interaction of the flapping motion of the wing and the surround- ing air. These forces generate lift and drag during flapping flight regime. As a result, they play a significant role in the description of the physics that underlies such systems. In this research we focus on dynamic and kinematic models that govern the motion of ground based robots that emulate flapping flight. The restriction to ground based biologically inspired robotic systems is predicated on two observations. First, it has become increasingly popular to design and fabricate bio-inspired robots for wind tunnel studies. Second, by restricting the robotic systems to be anchored in an inertial frame, the robotic equations of motion are well understood, and we can focus attention on flapping wing aerodynamics for such nonlinear systems. We study nonlinear modeling, identification, and control problems that feature the above complexities. This document summarizes research progress and plans that focuses on two key aspects of modeling, identification, and control of nonlinear dynamics associated with flapping flight. / Ph. D.

<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.