Dans certaines conditions de vol, les aéronefs à voilure tournante souffrent parfois de l’émergence d’oscillations indésirables, phénomènes potentiellement instables connus sous le nom de Couplages Pilote-Aéronef aéroélastiques (CPA). Ces phénomènes affectent de manière critique la sécurité et la performance des aéronefs. Par conséquent, il est important d’être capable de prédire l’émergence de tels phénomènes dynamiques, le plus tôt possible dans le processus de conception des hélicoptères. Une revue de la littérature révèle que ces phénomènes sont le résultat d’interactions entre les comportements biodynamique du pilote et aéroélastique des hélicoptères. Afin d’avoir une plus grande modularité et granularité dans le processus de modélisation de systèmes complexes, une approche par bond graphs est adoptée. Un modèle aéromécanique d’hélicoptère et un modèle neuro-musculo-squelettique d’un des membres supérieurs du pilote sont développés en bond graphs. Parmi les représentations proposées, trois sont originales, notamment afin de modéliser : des efforts aérodynamiques quasi-statiques, la liaison traînée-battement-pas entre pale et moyeu rotor, et les efforts musculaires à partir d’un modèle de Hill qui tient compte d’une boucle de rétroaction neuromusculaire. Des résultats encourageants sont obtenus lorsque l’on compare la transmissibilité, entre l’angle de manche de pas cyclique imposé par le pilote et des accélérations latérales de la cabine, calculée à partir du modèle biodynamique, et à partir des résultats expérimentaux tirés de la littérature. Un modèle du système bioaéroélastique homme-machine est linéarisé, au voisinage d’un vol stationnaire, et analysé en termes de stabilité. L’étude révèle, comme conjecturé dans la littérature, que le mode régressif de traînée peut être déstabilisé. De plus, il apparaît que le mode progressif de traînée peut également être déstabilisé lors d’un CPA sur l’axe latéral-roulis. Un critère d’analyse de la stabilité d’un équilibre d’un système dynamique à partir d’un modèle linéaire limite la possibilité de prendre en compte certains comportements non-linéaires et donc réduit l’espace de conception. Les premières pierres vers une méthode basée sur des fonctions de Chetaev sont posées, afin de déterminer si l’équilibre d’un système dynamique est instable, directement à partir d’un modèle mathématique non-linéaire de grande dimension, à un coût de calcul potentiellement intéressant. Afin d’illustrer la pertinence de la proposition, le cas de la résonance sol d’un hélicoptère est présentée. / Under certain flight conditions, rotorcrafts might suffer from the emergence of undesirable oscillations, potentially unstable phenomena, known as aeroelastic Rotorcraft-Pilot Couplings (RPCs). These phenomena critically affect the safety and performance of rotorcraft designs. Therefore, there is an important interest in being able to predict the emergence of such dynamic phenomena, as soon as possible during the design process of helicopters. A review of the state-of-the-art reveals that these phenomena are the result of interactions between pilots’ biodynamics and helicopters’ aeroelastic behaviors. In order to provide more modularity and granularity in the modeling of complex systems, a bond graph based approach is used. A helicopter aeromechanical model and a pilot upper limb neuromusculoskeletal model are developed using bond graphs. Three original bond graph representations are proposed, to model: quasi-steady aerodynamic forces, lag-flap-pitch joint at blades’ roots, and a Hill-type muscle force model that accounts for muscle reflexive feedback. Encouraging results are found when comparing the pilot biodynamic model transmissibility cyclic lever angle to lateral cockpit accelerations computations to literature experimental results. A linear model of the coupled human-machine bioaeroelastic system around hover is analyzed in terms of stability. It reveals not only the regressing lag mode, as conjectured in literature, but also the advancing lag mode can be destabilized during a lateral-roll aeroelastic RPC. Furthermore, a criterion to assess the stability of the equilibrium of a dynamic system from a linear model limits the possibility to take into account nonlinear physical behaviors, reducing the design space. The first blocks towards a method based on Chetaev functions is proposed, to determine if an equilibrium is unstable, directly from its large nonlinear mathematical model, at a potentially interesting computational cost. The helicopter ‘ground resonance’ case illustrates the soundness of the proposal.

Für die elastische Mehrkörpersimulation bzw. die FEM-MKS-Kopplung sind reduzierte FE-Modelle von großer Bedeutung. Die Erstellung reduzierter Modelle mit hoher Abbildungsgüte im Rahmen einer Modellordnungsreduktion erfordert einerseits ein geeignetes Reduktions-verfahren und andererseits zuverlässige Korrelationsmethoden. Beides wird durch die Soft-ware MORPACK bereitgestellt. Die Korrelation reduzierter FE-Modelle basiert in MORPACK derzeit ausschließlich auf modalen Eigenschaften. Ausgehend von der Annahme, dass sich die Abbildungsgüte eines reduzierten FE-Modells erst im Rahmen einer Zeitbereichssimula-tion vollständig beurteilen lässt, ist eine dahingehende Erweiterung von MORPACK geplant. Für einfache Topologien muss die Möglichkeit bestehen, das dynamische Verhalten, redu-zierter Modelle, direkt in MORPACK zu simulieren. Mit Hilfe der resultierenden Zeitsignale werden die reduzierten Modelle bewertet. Für die Umsetzung dieser Idee muss in MORPACK zunächst ein eigenständiges EMKS-Programm implementiert werden. Die Implementierung des EMKS-Programms in MORPACK (bzw. MATLAB) stellt den Schwerpunkt dieser Arbeit dar. Es werden zunächst die Anforderungen an das EMKS-Programm formuliert. Nach der Behandlung aller erforderlichen theoretischen Grundlagen werden die Systemgleichungen hergeleitet. Anschließend wird ein Formalismus bereitgestellt, der den Aufbau der Systemgleichungen, auf Basis der Nutzereingaben ermöglicht. Nach der Implementierung des Formalismus wird das EMKS-Programm verifiziert und erprobt.:1 Einleitung 1 1.1 Motivation 1 1.2 Zielsetzung 2 1.3 Lösungsweg 3 2 Verifikation und Optimierung durch Zeitbereichssimulationen 5 2.1 Erweiterung von MORPACK 5 2.2 Anforderungen an das EMKS-Programm 10 2.3 Korrelation von Zeitsignalen 12 3 Grundlagen der elastischen Mehrkörpersimulation 16 3.1 Berücksichtigung elastischer Deformationen in Mehrkörpersystemen 16 3.2 Kinematik freier Einzelkörper 19 3.2.1 Räumliche Drehungen von Bezugssystemen 19 3.2.2 Methode des bewegten Bezugssystems 23 3.2.3 Diskretisierung und Variablen für die Zustandsbeschreibung 25 3.2.4 Kinematik der Schnittstellenknoten 28 3.3 Kinetik freier Einzelkörper 31 3.4 Wahl des Körperbezugssystems 40 3.4.1 Kinematische Zwangsbedingungen 40 3.4.2 Kinetische Zwangsbedingungen 42 3.5 Gebundene Mehrkörpersysteme 44 3.6 Daten von elastischen Körpern 48 4 Bewegungsgleichungen und EMKS Formalismus für zwei beliebig gekoppelte Körper 52 4.1 Modellbildung 52 4.2 Bewegungsgleichungen in einem Satz natürlicher Koordinaten 54 4.3 Transformation auf Minimalkoordinaten 62 4.3.1 Formalismus 63 4.3.2 Herleitung der notwendigen Vektoren und Matrizen 65 5 Erweiterung des EMKS-Algorithmus für die festgelegte Topologie 76 6 Implementierung in MORPACK 84 6.1 Struktur der Eingabe- und Definitionsdaten 84 6.2 Grafische Benutzeroberfläche und Einbindung in MORPACK 90 6.3 Implementierung des EMKS-Formalismus 92 7 Verifikation und Erprobung 98 7.1 Verifikation mit SIMPACK 98 7.2 Erprobung der Prozesskette 101 7.2.1 Erprobungsmodell 101 7.2.2 Ergebnisse der Zeitbereichssimulation im Vergleich zu modalen Korrelationskriterien 103 7.2.3 Optimierung durch Zeitbereichssimulation 108 8 Zusammenfassung und Ausblick 112 / Reduced FE-Models are very important for elastic multibody simulation and FEM-MKS-coupling. The generation of reduced FE-models with high approximation quality in a model order reduction requires on the one hand a suitable reduction method and on the other hand reliable correlation methods. Both are provided by the MORPACK software. In MORPACK the correlation of reduced FE models based currently only on modal properties. An extension of the MORPACK software is planned on the assumption, that the approximation quality of a reduced FE-model can be completely assessed only in a time domain simulation. For simple topologies, it must be possible to simulate the dynamic behavior of reduced models directly into MORPACK. With the correlation of resulting time signals, the reduced models are as-sessed. To realize this idea, an independent EMKS program must be implemented in MORPACK. The implementation of the EMKS program in MORPACK (respectively MATLAB) represents the focus of this thesis. The first part is to formulate the necessary requirements for the EMKS program. After handling of all the necessary theoretical foundations, the system equa-tions are derived. Subsequently, formalism is provided that allows a construction of the sys-tem equations based on the user input. After the implementation of the formalism, the EMKS program will verify and tested.:1 Einleitung 1 1.1 Motivation 1 1.2 Zielsetzung 2 1.3 Lösungsweg 3 2 Verifikation und Optimierung durch Zeitbereichssimulationen 5 2.1 Erweiterung von MORPACK 5 2.2 Anforderungen an das EMKS-Programm 10 2.3 Korrelation von Zeitsignalen 12 3 Grundlagen der elastischen Mehrkörpersimulation 16 3.1 Berücksichtigung elastischer Deformationen in Mehrkörpersystemen 16 3.2 Kinematik freier Einzelkörper 19 3.2.1 Räumliche Drehungen von Bezugssystemen 19 3.2.2 Methode des bewegten Bezugssystems 23 3.2.3 Diskretisierung und Variablen für die Zustandsbeschreibung 25 3.2.4 Kinematik der Schnittstellenknoten 28 3.3 Kinetik freier Einzelkörper 31 3.4 Wahl des Körperbezugssystems 40 3.4.1 Kinematische Zwangsbedingungen 40 3.4.2 Kinetische Zwangsbedingungen 42 3.5 Gebundene Mehrkörpersysteme 44 3.6 Daten von elastischen Körpern 48 4 Bewegungsgleichungen und EMKS Formalismus für zwei beliebig gekoppelte Körper 52 4.1 Modellbildung 52 4.2 Bewegungsgleichungen in einem Satz natürlicher Koordinaten 54 4.3 Transformation auf Minimalkoordinaten 62 4.3.1 Formalismus 63 4.3.2 Herleitung der notwendigen Vektoren und Matrizen 65 5 Erweiterung des EMKS-Algorithmus für die festgelegte Topologie 76 6 Implementierung in MORPACK 84 6.1 Struktur der Eingabe- und Definitionsdaten 84 6.2 Grafische Benutzeroberfläche und Einbindung in MORPACK 90 6.3 Implementierung des EMKS-Formalismus 92 7 Verifikation und Erprobung 98 7.1 Verifikation mit SIMPACK 98 7.2 Erprobung der Prozesskette 101 7.2.1 Erprobungsmodell 101 7.2.2 Ergebnisse der Zeitbereichssimulation im Vergleich zu modalen Korrelationskriterien 103 7.2.3 Optimierung durch Zeitbereichssimulation 108 8 Zusammenfassung und Ausblick 112

Indiana University-Purdue University Indianapolis (IUPUI) / In this thesis, a high- fidelity multibody dynamics model of a backhoe for simulating the digging operation is developed using the DIS (Dynamic Interactions Simulator)multibody dynamics software. Sand is used as a sample digging material to illustrate the model. The backhoe components (such as frame, manipulators links,track segments, wheels and sprockets) are modeled as rigid bodies. The geometry of the major moving components of the backhoe is created using the Pro/E solid modeling software. The components of the backhoe are imported to DIS and connected using joints (revolute, cylindrical and prismatic joints). Rotary and linear actuators along with PD (Proportional-Derivative) controllers are used to move and steer the backhoe and to move the backhoes manipulator in the desired trajectory. Sand is modeled using cubic shaped particles that can come into contact with each other, the backhoes bucket and ground. A cubical sand particle contact surface is modeled using eight spheres that are rigidly glued to each other to form a cubical shaped particle, The backhoe and ground surfaces are modeled as polygonal surfaces. A penalty technique is used to impose both joint and normal contact constraints (including track-wheels, track-terrain, bucket-particles and particles-particles contact). An asperity-based friction model is used to model joint and contact friction. A Cartesian Eulerian grid contact search algorithm is used to allow fast contact detection between particles. A recursive bounding box contact search algorithm is used to allow fast contact detection for polygonal contact surfaces and is used to detect contact between: track and ground; track and wheels; bucket and particles; and ground and particles. The governing equations of motion are solved along with joint/constraint equations using a time-accurate explicit solution procedure. The sand model is validated using a conical hopper sand flow experiment in which the sand flow rate during discharge and the angle of repose of the resulting sand pile are experimentally measured. The results of the conical hopper simulation are compared with previously published experimental results. Parameter studies are performed using the sand model to study the e ffects of the particle size and the orifi ces diameter of the hopper on the sand pile angle of repose and sand flow rate. The sand model is integrated with the backhoe model to simulate a typical digging operation. The model is used to predict the manipulators actuator forces needed to dig through a pile of sand. Integrating the sand model and backhoe model can help improving the performance of construction equipment by predicting, for various vehicle design alternatives: the actuator and joint forces, and the vehicle stability during digging.

The main objective of this work is to demonstrate some new computational methods for estimation, optimization and modeling of dynamical systems that use automatic differentiation. Particular focus will be upon dynamical systems arising in Aerospace Engineering. Automatic differentiation is a recursive computational algorithm, which enables computation of analytically rigorous partial derivatives of any user-specified function. All associated computations occur, in the background without user intervention, as the name implies. The computational methods of this dissertation are enabled by a new automatic differentiation tool, OCEA (Object oriented Coordinate Embedding Method). OCEA has been recently developed and makes possible efficient computation and evaluation of partial derivatives with minimal user coding. The key results in this dissertation details the use of OCEA through a number of computational studies in estimation and dynamical modeling. Several prototype problems are studied in order to evaluate judicious ways to use OCEA. Additionally, new solution methods are introduced in order to ascertain the extended capability of this new computational tool. Computational tradeoffs are studied in detail by looking at a number of different applications in the areas of estimation, dynamical system modeling, and validation of solution accuracy for complex dynamical systems. The results of these computational studies provide new insights and indicate the future potential of OCEA in its further development.

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