Quaternion Regularization of the Equations of the Perturbed Spatial Restricted Three-Body Problem: I

We develop a quaternion method for regularizing the differential equations of the perturbed spatial restricted three-body problem by using the Kustaanheimo–Stiefel variables, which is methodologically closely related to the quaternion method for regularizing the differential equations of perturbed spatial two-body problem, which was proposed by the author of the present paper.A survey of papers related to the regularization of the differential equations of the two- and threebody problems is given. The original Newtonian equations of perturbed spatial restricted three-body problem are considered, and the problem of their regularization is posed; the energy relations and the differential equations describing the variations in the energies of the system in the perturbed spatial restricted three-body problem are given, as well as the first integrals of the differential equations of the unperturbed spatial restricted circular three-body problem (Jacobi integrals); the equations of perturbed spatial restricted three-body problem written in terms of rotating coordinate systems whose angular motion is described by the rotation quaternions (Euler (Rodrigues–Hamilton) parameters) are considered; and the differential equations for angular momenta in the restricted three-body problem are given.Local regular quaternion differential equations of perturbed spatial restricted three-body problem in the Kustaanheimo–Stiefel variables, i.e., equations regular in a neighborhood of the first and second body of finite mass, are obtained. The equations are systems of nonlinear nonstationary eleventhorder differential equations. These equations employ, as additional dependent variables, the energy characteristics of motion of the body under study (a body of a negligibly small mass) and the time whose derivative with respect to a new independent variable is equal to the distance from the body of negligibly small mass to the first or second body of finite mass.The equations obtained in the paper permit developing regular methods for determining solutions, in analytical or numerical form, of problems difficult for classicalmethods, such as the motion of a body of negligibly small mass in a neighborhood of the other two bodies of finite masses.     

Reibung, Selbsthemmung, Selbsterregung: Überlegungen zum Painlevéschen Paradoxon

Equations of motion of certain rigid body mechanisms which contain Coulomb friction may fail for particular parameter regions, Painlevé's paradoxical results emerge. When small deformations are permitted, a set of singularly perturbed algebro-differential equations governs the motions. The number characteristic for the paradox discriminates between decaying small high-frequency oscillations and oscillations which possibly grow. In order to decide about stability in the second case, it is necessary to smoothen the frictional characteristic and to study the coupled high-frequency oscillations of the complete system. The paper outlines the individual steps of such an investigation for a simple example.     

On the rotational motion of a gyrostat about a fixed point with mass distribution

The perturbed rotational motion of a gyrostat about a fixed point with mass distribution near to Lagrange’s case is investigated. The gyrostat is subjected under the influence of a variable restoring moment vector, a perturbing moment vector, and a third component of a gyrostatic moment vector. It is assumed that the angular velocity of the gyrostat is sufficiently large, its direction is close to the axis of dynamic symmetry, and the perturbing moments are small as compared to the restoring ones. These assumptions permit us to introduce a small parameter. Averaged systems of the equations of motion in the first and second approximations are obtained. Also, the evolution of the precession angle up to the second approximation is determined. The graphical representations of the nutation and precession angles are presented to describe the motion at any time.     

A comparison between a collocation and weak implementation of the rigid‐body motion constraint on a particle surface

In this work, we implemented and compared two different methods to impose the rigid‐body motion constraint on a solid particle moving inside a fluid. We consider a fictitious domain method to easily manage the particle motion. As the solid as well as the fluid inertia are neglected, the particle can be discretized through its boundary only. The rigid‐body motion is imposed via Lagrange multipliers on the boundary. In the first method, such constraints are imposed in discrete points on the boundary (collocation), whereas in the second the constraint is imposed in a weak way on elements dividing the particle surface. Two test problems, that is, a spherical and an ellipsoidal particle in a sheared Newtonian fluid, are chosen to compare the methods. In both cases, the analysis is carried out in 2D as well as in 3D. The results show that for the collocation method an optimal number of collocation points exist leading to the smallest error. However, small variations in the optimal value can generate large deviations. In the weak implementation, the error is only mildly affected by the number of elements used to discretize the particle boundary and by the Lagrange multiplier's interpolation space. A further analysis is carried out to study the effect of an approximated integration of weak constraints. A comparison between the two methods showed that the same accuracy can be achieved by using less constraints if the weak discretization is used. Finally, the rigid‐body motion imposed via weak constraints leads to better conditioned linear systems. Copyright © 2009 John Wiley & Sons, Ltd.     

Accelerating motion of a sphere in a maxwell fluid

The translatory accelerating motion of a sphere due to an arbitrarily applied force in an unlimited Maxwell fluid is considered. The exact solutions for the velocity of the sphere for three particular types of accelerating motion are presented. The first is for a falling sphere; the second is for the decelerating motion of a sphere after the force which maintains the sphere with a constant velocity is removed; the third is for the motion of the sphere subjected to an impulsive force. The exact solutions are expressed in terms of real, regular, definite integrals which can be evaluated by numerical technique. Also presented are the asymptotic solutions for the velocity of the sphere in all three cases which are valid for small values of time.     

Stability of high-frequency periodic motions of a heavy rigid body with a horizontally vibrating suspension point

The motion of a heavy rigid body one of whose points (the suspension point) executes horizontal harmonic high-frequency vibrations with small amplitude is considered. The problem of existence of high-frequency periodic motions with period equal to the period of the suspension point vibrations is considered. The stability conditions for the revealed motions are obtained in the linear approximation. The following three special cases of mass distribution in the body are considered; a body whose center of mass lies on the principal axis of inertia, a body whose center of mass lies in the principal plane of inertia, and a dynamically symmetric body.     

Melnikov's method for rigid bodies subject to small perturbation torques

Summary In this paper, the global motion of rigid bodies subjected to small perturbation torques, either conservative or dissipative, is investigated by means of Melnikov's method. Deprit's variables are introduced to transform the equations of motion into a standard form which is rendered suitable for the application of Melnikov's method. The Melnikov method is used to predict the transversal intersections of stable and unstable manifolds for the pertubed rigid-body motion. The chosen examples are a self-excited rigid body subject to a small periodic torque in a viscous medium, and the heavy rigid body. It is shown in both cases that there exist transversal intersections of heteroclinic orbits for certain ranges of parameter values.     

Impact of a spherical rigid body on the surface of a cavity in a compressible liquid: An axisymmetric problem

The shock interaction of a spherical rigid body with a spherical cavity is studied. This nonstationary mixed boundary-value problem with an unknown boundary is reduced to an infinite system of linear Volterra equations of the second kind and the differential equation of motion of the body. The hydrodynamic and kinematic characteristics of the process are obtained __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 1, pp. 11–19, January 2008.     

Motion of a small sphere in a nonhomogeneous flow of an incompressible liquid

We consider the motion of a small sphere in an arbitrary potential flow of an ideal liquid. For the general case we obtain an integral of the equations of motion and a particular solution. We find flows in which the force acting on the sphere is central. We also obtain exact solutions of the equations of motion of the sphere for the cases of stationary flows around a cylinder and around a body of revolution when the forces are noncentral. N. E. Zhukovskii [1] calculated the force acting on a fixed sphere in an arbitrary nonstationary flow. Kelvin [2] obtained the equations of motion of a sphere in a stationary flow of a liquid circulating through a hole in a solid. A formula for the force F, acting on a fixed small body of volume V in a stationary flow with speed v, was obtained by Taylor [3]: F = (T ₀ / v)Vv + 1/2_V v ² Here T₀ is the kinetic energy of an unbounded liquid in which a body moves with velocity v. This problem was solved in [3] through a direct integration of the pressure forces over the surface of the body in a flow defined by multipoles of the first and second orders at infinity.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 57–61, September–October, 1973.     

Poynting oscillations of a rigid disk supported by a neo-Hookean shaft

Simultaneous axial and torsional oscillations of a rigid disk attached to an elastomeric shaft are investigated. Five cases are solved exactly. The uncoupled, small amplitude axial and torsional oscillations of the disk are investigated for neo-Hookean and Mooney-Rivlin shafts with static stretch. The finite torsional vibration of the load superimposed on a static stretch of the shaft is studied for the Mooney-Rivlin model. Solutions for both small and finite amplitude, uniaxial vibrations of the body superimposed on a pretwisted neo-Hookean shaft with static stretch are derived. Simple bounds on the period for the finite motion are provided; and various universal frequency relations for neo-Hookean and Mooney-Rivlin materials are identified.Finally, the main problem of finite, uniaxial vibrations accompanied by a small twisting motion is studied for the neo-Hookean model. The exact periodic solution for the axial response is obtained; and the coupled, small torsional motion is then determined by Hill's equation. A stability criterion for the Mathieu-Hill equation is used to obtain stability maps in a physical parameter space. Geometrical conditions sufficient for universal stability of the motion are read from this graph. Instability of the torsional oscillation, the beating phenomenon and exchange of energies, and the relation of the stability diagram to amplitude bounds on the uncoupled, linearized motion sufficient to assure universal stability predicted for small amplitude vibrations, are discussed and described graphically with the aid of a numerical model. It is shown that an unstable configuration may be stabilized by increasing the diameter of the disk.     

A Generalized Formulation of Elastodynamics: Small on Rigid

The conventional theory of linearized elastodynamics addresses the case of motions that have small displacement gradients with respect to a reference configuration of the elastic body that is unstressed and at rest. Here, we develop a theory of much wider applicability in which the linearization is with respect to a reference configuration that is in rigid motion. More specifically, with an eye toward application of the theory to analysis of the motions of relatively inflexible machine parts, we view the motion as being composed of a rigid motion, which corresponds to the applied loads and initial conditions for the body under consideration, and an infinitesimal motion, in which the displacement from the rigid motion has a small gradient.     

The supersonic motion of a rigid body in an elastic medium with friction

The plane problem of supersonic steady motion of a body in an elastic medium is solved. Two possible cases of body motion are considered depending on its velocity. In the first case, the body moves at a velocity greater than the velocity of transverse waves but smaller than the velocity of longitudinal waves. In the second case, the body moves at a velocity greater than the velocity of longitudinal waves. An analytic solution of the problem under study is obtained and analyzed. It is shown that friction substantially influences the penetration process.     

The Method of Averaging for Euler's Equations of Rigid Body Motion

We formulate the method of averaging for perturbations of Euler's equations of rotational motion. Euler's equations are three strongly nonlinear coupled differential equations that can be viewed as a three dimensional oscillator. The method of averaging is used to determine the long-term influence of perturbation terms on the motion by averaging about the nominal rigid body motion. The treatment is applicable to a large class of motions including precession with large nutation – it is not restricted to small motions about simple spins or nearly axi-symmetric bodies. Three examples are shown that demonstrate the accuracy of the method's predictions.     

Quasi-Euler-Poinsot motion of a rigid body

The phase-plane method of nonlinear oscillation is used to discuss the influence of the small dissipation upon the Euler-Poinsot motion of a rigid body about a fixed point. The equations of phase coordinates are applied instead of Eulerian equations, and the global characteristics of the motion of rigid body are analysed according to the distribution and the type of the singular points. A Chaplygin's sphere on a rough plane, a rigid body in viscous medium and one with a cavity filled with viscous fluid are discussed as examples. It is shown that the motions of rigid bodies dissipated by various physical factors have a common qualitative character. The rigid body tends to make a permanent rotation about the principal axis of the largest moment of inertia. The transitive process can change from oscillatory to aperiodic with the decrease in dissipation.     

Impact of a body with a plane bottom on a thin liquid layer at a small angle

The problem of the impact of a body with a plane bottom (of the type of a box) on a thin liquid layer at a small angle is solved in the two-dimensional formulation. The nonlinear shallow water equations are used, together with the method of matched asymptotic expansions. It is found that at certain values of the input parameters of the problem the liquid pressure diminishes near the lower end of the body and becomes smaller than the atmospheric pressure, which results in liquid separation from the box bottom. The numerical results show that all input parameters of the problem have a considerable effect on the nature of body motion. The liquid separation effect on body motion is analyzed.     

External stability of resonances in the motion of an asymmetric rigid body with a strong magnet in the geomagnetic field

We consider a precession motion, close to the classical Lagrange case, of an asymmetric rigid body with a strong magnet in an orbit in the geomagnetic field. For the principal moment we take the restoring torque due to the interaction between the planet magnetic fields and the rigid body. The perturbing actions are due to small moments of the rigid body mass-inertial asymmetry and small constant moments. We show that these perturbations result in the realization of secondary resonance effects in the rotational motion of the rigid body caused by the influence of resonance denominators in higher-order approximations of the averaging method. These effects were discovered in the study of rotational motion of a satellite with a magnetic damper in the nearly Euler case. In the present paper, we analyze both the secondary resonance effects themselves and the external stability of resonances. We obtain conditions ensuring a decrease in the angular velocity of the rigid body rotation about its center of mass. We also discover several new laws of influence of resonances on the nonresonance evolution of slow variables, which is related to the appearance of stable resonances.     

Taking into account the influence of elastic elements of a free structure on the accuracy of its stabilization under a bang-bang control

We consider the problem of stabilization with respect to a prescribed position for the translational motion of a rigid body with interior material points connected with each other and with the exterior body by linear viscoelastic constraints. The motion occurs under the action of a constant exterior perturbation and a bang-bang control force that are directed along the line of motion. We assume that the bang-bang force control channel has a fixed delay, so that arbitrarily frequent switchings are impossible. We suggest a positional control ensuring the solution of this problem. We estimate the amplitude of the rigid body vibrations about the center of mass of the entire structure and the accuracy of stabilization of the prescribed position of the rigid body depending on the mechanical characteristics of the system and the control force magnitude. We also consider the problem of maximizing the stabilization accuracy depending on the control parameters. By way of example, we consider the controlled motion of a two-mass oscillatory system. This work is closely related to [1–3] and continues the studies of the guaranteed optimal bang-bang controllers with delay in the control channel [4–9]. The dynamics of a rigid body with elastic and dissipative elements was studied in [10] under the assumption that the period of natural vibrations and their decay time are small compared with the characteristic time of motion.