Multibody Dynamics on Differentiable Manifolds

2021 ◽  
Vol 16 (4) ◽  
Author(s):  
Edward J. Haug

Abstract Topological and vector space attributes of Euclidean space are consolidated from the mathematical literature and employed to create a differentiable manifold structure for holonomic multibody kinematics and dynamics. Using vector space properties of Euclidean space and multivariable calculus, a local kinematic parameterization is presented that establishes the regular configuration space of a multibody system as a differentiable manifold. Topological properties of Euclidean space show that this manifold is naturally partitioned into disjoint, maximal, path connected, singularity free domains of kinematic and dynamic functionality. Using the manifold parameterization, the d'Alembert variational equations of multibody dynamics yield well-posed ordinary differential equations of motion on these domains, without introducing Lagrange multipliers. Solutions of the differential equations satisfy configuration, velocity, and acceleration constraint equations and the variational equations of dynamics, i.e., multibody kinematics and dynamics are embedded in these ordinary differential equations. Two examples, one planar and one spatial, are treated using the formulation presented. Solutions obtained are shown to satisfy all three forms of kinematic constraint to within specified error tolerances, using fourth-order Runge–Kutta numerical integration methods.

Author(s):  
Edward J. Haug

Abstract The manipulator differentiable manifold formulation presented in Part I of this paper is used to create algorithms for forward and inverse kinematics on maximal, singularity free, path connected manifold components. Existence of forward and inverse configuration mappings in manifold components is extended to obtain forward and inverse velocity mappings. Computational algorithms for forward and inverse configuration and velocity analysis on a time grid are derived for each of the four categories of manipulator treated. Manifold parameterizations derived in Part I are used to transform variational equations of motion in Cartesian generalized coordinates to second order ordinary differential equations of manipulator dynamics, in terms of both input and output coordinates, avoiding ad-hoc derivation of equations of motion. This process is illustrated in evaluating terms required for equations of motion of the four model manipulators analyzed in Part I. It is shown that computation involved in forward and inverse kinematics and in evaluation of equations of manipulator dynamics can be carried out in real-time on modern microprocessors, supporting on-line implementation of modern methods of manipulator control.


2005 ◽  
Vol 128 (2) ◽  
pp. 359-364
Author(s):  
C. Q. Liu ◽  
Fang Li ◽  
R. L. Huston

Governing dynamical equations of motion for a basketball rolling on the rim of a basket are developed and presented. These equations form a system of five first-order, ordinary differential equations. Given suitable initial conditions, these equations are readily integrated numerically. The results of these integrations predict the success (into the basket) or failure (off the outside of the rim) of the basketball shot. A series of examples are presented. The examples show that minor changes in the initial conditions can produce major changes in the subsequent ball motion. Shooting and coaching strategies are recommended.


2014 ◽  
Vol 2014 ◽  
pp. 1-8 ◽  
Author(s):  
Zhaolin Jiang ◽  
Tingting Xu ◽  
Fuliang Lu

The skew-circulant matrix has been used in solving ordinary differential equations. We prove that the set of skew-circulants with complex entries has an idempotent basis. On that basis, a skew-cyclic group of automorphisms and functional equations on the skew-circulant algebra is introduced. And different operators on linear vector space that are isomorphic to the algebra ofn×ncomplex skew-circulant matrices are displayed in this paper.


2011 ◽  
Vol 141 ◽  
pp. 98-102
Author(s):  
Zong Yu Chang ◽  
Lei Wang ◽  
Da Xiao Gao ◽  
Zhong Qiang Zheng

Deployment of buoy systems is one of the most important procedures for buoy system’s operation. In this paper, buoy system with surface buoy, cable segments with components, anchor and so on is modeled by applying multibody dynamics method. Then numerical method is used to solve the ordinary differential equations and dynamics simulations are achieved while anchor is casting from board. The trajectories, velocity and acceleration of different nodes in buoy system are obtained. The transient tension force of each part of cable is analyzed in the process of deployment. This work is helpful for design and deployment planning of buoy system.


1998 ◽  
Vol 5 (5-6) ◽  
pp. 277-288 ◽  
Author(s):  
Ali H. Nayfeh ◽  
Haider N. Arafat

The nonlinear planar response of cantilever metallic beams to combination parametric and external subcombination resonances is investigated, taking into account the effects of cubic geometric and inertia nonlinearities. The beams considered here are assumed to have large length-to-width aspect ratios and thin rectangular cross sections. Hence, the effects of shear deformations and rotatory inertia are neglected. For the case of combination parametric resonance, a two-mode Galerkin discretization along with Hamilton’s extended principle is used to obtain two second-order nonlinear ordinary-differential equations of motion and associated boundary conditions. Then, the method of multiple scales is applied to obtain a set of four first-order nonlinear ordinary-differential equations governing the modulation of the amplitudes and phases of the two excited modes. For the case of subcombination resonance, the method of multiple scales is applied directly to the Lagrangian and virtual-work term. Then using Hamilton’s extended principle, we obtain a set of four first-order nonlinear ordinary-differential equations governing the amplitudes and phases of the two excited modes. In both cases, the modulation equations are used to generate frequency- and force-response curves. We found that the trivial solution exhibits a jump as it undergoes a subcritical pitchfork bifurcation. Similarly, the nontrivial solutions also exhibit jumps as they undergo saddle-node bifurcations.


Author(s):  
Abdolreza Pasharavesh ◽  
Y. Alizadeh Vaghasloo ◽  
A. Fallah ◽  
M. T. Ahmadian

In this study vibration amplitude, frequency and damping of a microbeam is controlled using a RLC block containing a capacitor, resistor and inductor in series with the microbeam. Applying this method all of the considerable characteristics of the oscillatory system can be determined and controlled with no change in the geometrical and physical characteristics of the microbeam. Euler-Bernoulli assumptions are made for the microbeam and the electrical current through the microbeam is computed by considering the microbeam deflection and its voltage. Considering the RLC block, the voltage difference between the microbeam and the substrate is calculated. Two coupled nonlinear partial differential equations are obtained for the deflection and the voltage. The one parameter Galerkin method is employed to transform the equations of motion to a set of nonlinear coupled ordinary differential equations. Differential quadrature method (DQM) is implemented to solve the governing nonlinear ordinary differential equations. The effect of the controller parameters such as capacitance, resistance and inductance on the amplitude, frequency and damping is studied. Also the internal resonance between the electrical and mechanical parts of the system is studied. Results indicate using these elements, amplitude, frequency and damping can be controlled as desired by the user.


Author(s):  
Edward J. Haug

A method is presented for formulating and numerically integrating ordinary differential equations of motion for nonholonomically constrained multibody systems. Tangent space coordinates are defined in configuration and velocity spaces as independent generalized coordinates that serve as state variables in the formulation, yielding ordinary differential equations of motion. Orthogonal-dependent coordinates and velocities are used to enforce constraints at position, velocity, and acceleration levels. Criteria that assure accuracy of constraint satisfaction and well conditioning of the reduced mass matrix in the equations of motion are used as the basis for updating local coordinates on configuration and velocity constraint manifolds, transparent to the user and at minimal computational cost. The formulation is developed for multibody systems with nonlinear holonomic constraints and nonholonomic constraints that are linear in velocity coordinates and nonlinear in configuration coordinates. A computational algorithm for implementing the approach is presented and used in the solution of three examples: one planar and two spatial. Numerical results using a fifth-order Runge–Kutta–Fehlberg explicit integrator verify that accurate results are obtained, satisfying all the three forms of kinematic constraint, to within error tolerances that are embedded in the formulation.


2021 ◽  
Vol 83 (2) ◽  
pp. 151-159
Author(s):  
E.A. Korovaytseva

Results of hyperelastic soft shells nonlinear axisymmetric dynamic deforming problems solution algorithm testing are represented in the work. Equations of motion are given in vector-matrix form. For the nonlinear initial-boundary value problem solution an algorithm which lies in reduction of the system of partial differential equations of motion to the system of ordinary differential equations with the help of lines method is developed. At this finite-difference approximation of partial time derivatives is used. The system of ordinary differential equations obtained as a result of this approximation is solved using parameter differentiation method at each time step. The algorithm testing results are represented for the case of pressure uniformly distributed along the meridian of the shell and linearly increasing in time. Three types of elastic potential characterizing shell material are considered: Neo-hookean, Mooney – Rivlin and Yeoh. Features of numerical realization of the algorithm used are pointed out. These features are connected both with the properties of soft shells deforming equations system and with the features of the algorithm itself. The results are compared with analytical solution of the problem considered. Solution behavior at critical pressure value is investigated. Formulations and conclusions given in analytical studies of the problem are clarified. Applicability of the used algorithm to solution of the problems of soft shells dynamic deforming in the range of displacements several times greater than initial dimensions of the shell and deformations much greater than unity is shown. The numerical solution of the initial boundary value problem of nonstationary dynamic deformation of the soft shell is obtained without assumptions about the limitation of displacements and deformations. The results of the calculations are in good agreement with the results of analytical studies of the test problem.


1999 ◽  
Vol 121 (1) ◽  
pp. 110-113 ◽  
Author(s):  
E. Wong ◽  
J. W. Zu

The dynamic behavior of a simply-supported spinning Timoshenko shaft with coupled bending and torsion is analyzed. This is accomplished by transforming the set of nonlinear partial differential equations of motion into a set of linear ordinary differential equations. This set of ordinary differential equations is a time-varying system and the solution is obtained analytically in terms of Chebyshev series. A beating phenomoenon is observed from the numerical simulations, which is not observed for shaft systems where only bending vibration is considered.


Sign in / Sign up

Export Citation Format

Share Document