Even Order Explicit Symplectic Geometric Algorithms for Quaternion Kinematical Differential Equation in Guidance Navigation and Control via Diagonal Pad\`{e} Approximation and Cayley Transform

The Quaternion kinematical differential equation (QKDE) plays a key role in navigation, control and guidance systems. Although explicit symplectic geometric algorithms (ESGA) for this problem are available, there is a lack of a unified way for constructing high order symplectic difference schemes with configurable order parameter. We present even order explicit symplectic geometric algorithms to solve the QKDE with diagonal Pad\`{e} approximation and Cayley transform. The maximum absolute error for solving the QKDE is $\mathcal{O}(\tau^{2\ell})$ where $\tau$ is the time step and $\ell$ is the order parameter. The linear time complexity and constant space complexity of computation as well as the simple algorithmic structure show that our algorithms are appropriate for realtime applications in aeronautics, astronautics, robotics, visual-inertial odemetry and so on. The performance of the proposed algorithms are verified and validated by mathematical analysis and numerical simulation

Geometrijski integratori za numeričku simulaciju dinamike letjelica

U fokusu istraživanja je prilagodba i primjena novoprojektiranih geometrijskih integracijskih algoritama na numeričke zadaće dinamike letjelica te je ocijenjena njihova pogodnost i kvaliteta za primjenu u okviru numeričkih simulacijskih procedura za izračunavanje dinamičkih odziva zračnih i svemirskih letjelica. Za razliku od standardnih numeričkih metoda integracije običnih diferencijalnih jednadžbi koje operiraju u vektorskim prostorima, ciljani geometrijski algoritmi projektirani su na mnogostrukostima i Lievim grupama. Polazne geometrijske metode uobličene su na način da algoritamski zadovoljavaju određena kinematička i dinamička ograničenja koja proistječu iz diferencijalno-geometrijske strukture sustava. Nakon pojedinačnog opisa izvoda matematičkih modela odabranih algoritama i njihove prilagodbe za ciljane zadaće integracije dinamike letjelica, prezentirani su primjeri primjene geometrijskih integracijskih metoda te njihova usporedba s klasičnim procedurama vektorskih prostora

Biokinematic analysis of human body

Thesis (Doctoral)--Izmir Institute of Technology, Mechanical Engineering, Izmir, 2011Includes bibliographical references (leaves: 118-123)Text in English; Abstract: Turkish and Englishxiii, 123 leavesThis thesis concentrates on the development of rigid body geometries by using method of intersections, where simple geometric shapes representing revolute (R) and prismatic (P) joint motions are intersected by means of desired space or subspace requirements to create specific rigid body geometries in predefined octahedral fixed frame. Using the methodical approach, space and subspace motions are clearly visualized by the help of resulting geometrical entities that have physical constraints with respect to the fixed working volume. Also, this work focuses on one of the main areas of the fundamental mechanism and machine science, which is the structural synthesis of robot manipulators by inserting recurrent screws into the theory. After the transformation unit screw equations are presented, physical representations and kinematic representations of kinematic pairs with recurrent screws are given and the new universal mobility formulations for mechanisms and manipulators are introduced. Moreover the study deals with the synthesis of mechanisms by using quaternion and dual quaternion algebra to derive the objective function. Three different methods as interpolation approximation, least squares approximation and Chebyshev approximation is introduced in the function generation synthesis procedures of spherical four bar mechanism in six precision points. Separate examples are given for each section and the results are tabulated. Comparisons between the methods are also given. As an application part of the thesis, the most important elements of the human body and skeletal system is investigated by means of their kinematic structures and degrees of freedom. At the end of each section, an example is given as a mechanism or manipulator that can represent the behavior of the related element in the human body

Ineslatic analysis of geometrically exact rods

In this work, a formulation for rod structures, able to consider coupled geometric and constitutive sources of nonlinearity in both the static and the dynamic range, is developed. It is extended for allowing the inclusion of passive energy dissipating elements as a special rod element and geometric irregularities as a full three-dimensional body connected to the framed structure by means of a two-scale model. The proposed formulation is based on the Reissner-Simo geometrically exact formulation for rods considering an initially curved reference configuration and extended to include arbitrary distribution of composite materials in the cross sections. Each material point of the cross section is assumed to be composed of several simple materials with their own thermodynamically consistent constitutive laws. The simple mixing rule is used for treating the resulting composite. Cross sections are meshed into a grid of quadrilaterals, each of them corresponding to a fiber directed along the axis of the beam. A mesh independent response is obtained by means of the regularization of the energy dissipated at constitutive level considering the characteristic length of the mesh and the fracture energy of the materials. Local and global damage indices have been developed based on the ratio between the visco-elastic and nonlinear stresses. The consistent linearization of the weak form of the momentum balance equations is performed considering the effects of rate dependent inelasticity. Due to the fact that the deformation map belongs to a nonlinear manifold, an appropriated version of Newmark's scheme and of the iterative updating procedure of the involved variables is developed. The space discretization of the linearized problem is performed using the standard Galerkin finite element approach. A Newton-Raphson type of iterative scheme is used for the step-by-step solution of the discrete problem. A specific element for energy dissipating devices is developed, based on the rod model but releasing the rotational degrees of freedom. Appropriated constitutive relations are given for a wide variety of possible dissipative mechanisms. Several numerical examples have been included for the validation of the proposed formulation. The examples include elastic and inelastic finite deformation response of framed structures with initially straight and curved beams. Comparisons with existing literature is performed for the case of plasticity and new results are presented for degrading and composite materials. Those examples show how the present formulation is able to capture different complex mechanical phenomena such as the uncoupling of the dynamic response from resonance due to inelastic incursions and suppression of the high frequency content. The study of realistic flexible pre-cast and cast in place reinforced concrete framed structures subjected to static and dynamic actions is also carried out. Detailed studies regarding to the evolution of local damage indices, energy dissipation and ductility demands are presented. The studies include the seismic response of concrete structures with energy dissipating devices. Advantages of the use of passive control are verified

Constitutive and Geometric Nonlinear Models for the Seismic Analysis of RC Structures with Energy Dissipators

Nowadays, the use of energy dissipating devices to improve the seismic response of RC structures constitutes a mature branch of the innovative procedures in earthquake engineering. However, even though the benefits derived from this technique are well known and widely accepted, the numerical methods for the simulation of the nonlinear seismic response of RC structures with passive control devices is a field in which new developments are continuously preformed both in computational mechanics and earthquake engineering. In this work, a state of the art of the advanced models  for the numerical simulation of the nonlinear dynamic response of RC structures with passive energy dissipating devices subjected to seismic loading is made. The most commonly used passive energy dissipating devices are described, together with their dissipative mechanisms as well as with the numerical procedures used in modeling RC structures provided with such devices. The most important approaches for the formulation of beam models for RC structures are reviewed, with emphasis on the theory and numerics of formulations that consider both geometric and constitutive sources on nonlinearity. In the same manner, a more complete treatment is given to the constitutive nonlinearity in the context of fiber-like approaches including the corresponding cross sectional analysis. Special attention is paid to the use of damage indices able of estimating the remaining load carrying capacity of structures after a seismic action. Finally, nonlinear constitutive and geometric formulations for RC beam elements are examined, together with energy dissipating devices formulated as simpler beams with adequate constitutive laws. Numerical examples allow to illustrate the capacities of the presented formulations

Inelastic analysis of geometrically exact rods

Postprint (published version

New Directions for Contact Integrators

Contact integrators are a family of geometric numerical schemes which guarantee the conservation of the contact structure. In this work we review the construction of both the variational and Hamiltonian versions of these methods. We illustrate some of the advantages of geometric integration in the dissipative setting by focusing on models inspired by recent studies in celestial mechanics and cosmology.Comment: To appear as Chapter 24 in GSI 2021, Springer LNCS 1282

Proceedings of the ECCOMAS Thematic Conference on Multibody Dynamics 2015

This volume contains the full papers accepted for presentation at the ECCOMAS Thematic Conference on Multibody Dynamics 2015 held in the Barcelona School of Industrial Engineering, Universitat Politècnica de Catalunya, on June 29 - July 2, 2015. The ECCOMAS Thematic Conference on Multibody Dynamics is an international meeting held once every two years in a European country. Continuing the very successful series of past conferences that have been organized in Lisbon (2003), Madrid (2005), Milan (2007), Warsaw (2009), Brussels (2011) and Zagreb (2013); this edition will once again serve as a meeting point for the international researchers, scientists and experts from academia, research laboratories and industry working in the area of multibody dynamics. Applications are related to many fields of contemporary engineering, such as vehicle and railway systems, aeronautical and space vehicles, robotic manipulators, mechatronic and autonomous systems, smart structures, biomechanical systems and nanotechnologies. The topics of the conference include, but are not restricted to: ● Formulations and Numerical Methods ● Efficient Methods and Real-Time Applications ● Flexible Multibody Dynamics ● Contact Dynamics and Constraints ● Multiphysics and Coupled Problems ● Control and Optimization ● Software Development and Computer Technology ● Aerospace and Maritime Applications ● Biomechanics ● Railroad Vehicle Dynamics ● Road Vehicle Dynamics ● Robotics ● Benchmark ProblemsPostprint (published version