Construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body

We consider the problem of construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body so as to ensure the transition of the rigid body from an arbitrary initial angular position to the required final angular position. For the functionals to be minimized, we use combined performance functionals, one of which characterizes the expenditure of time and of the squared modulus of the angular momentum vector in a given proportion, while the other characterizes the expenditure of time and momentum of the modulus of the angular momentum vector necessary to change the rigid body orientation. The control (the vector of the rigid body angular momentum) is assumed to be bounded in the modulus. The problem is solved by using Pontryagin’s maximum principle and the quaternion differential equation [1, 2] relating the vector of the dynamically symmetric rigid body angular momentum to the quaternion of orientation of the coordinate system rotating with respect to the rigid body about its dynamical symmetry axis at an angular velocity proportional to the angular momentum vector projection on the axis. The use of such a model of rotational motion leads to the problem of optimal control with the moving right end of the trajectory and significantly simplifies the analytic study of the problem of construction of optimal laws of variation in the angular momentum vector, because this model explicitly exploits the body angular momentum quaternion (control) instead of the rigid body absolute angular velocity quaternion. We construct general analytic solutions of the differential equations for the boundary-value problems which form systems of nine nonlinear differential equations. It is shown that the process of solving the differential boundary-value problems is reduced to solving two scalar algebraic transcendental equations.     

An Internal Damping Model for the Absolute Nodal Coordinate Formulation

Introducing internal damping in multibody system simulations is important as real-life systems usually exhibit this type of energy dissipation mechanism. When using an inertial coordinate method such as the absolute nodal coordinate formulation, damping forces must be carefully formulated in order not to damp rigid body motion, as both this and deformation are described by the same set of absolute nodal coordinates. This paper presents an internal damping model based on linear viscoelasticity for the absolute nodal coordinate formulation. A practical procedure for estimating the parameters that govern the dissipation of energy is proposed. The absence of energy dissipation under rigid body motion is demonstrated both analytically and numerically. Geometric nonlinearity is accounted for as deformations and deformation rates are evaluated by using the Green–Lagrange strain–displacement relationship. In addition, the resulting damping forces are functions of some constant matrices that can be calculated in advance, thereby avoiding the integration over the element volume each time the damping force vector is evaluated.     

Biquaternion solution of the kinematic control problem for the motion of a rigid body and its application to the solution of inverse problems of robot-manipulator kinematics

The problem of reducing the body-attached coordinate system to the reference (programmed) coordinate system moving relative to the fixed coordinate system with a given instantaneous velocity screw along a given trajectory is considered in the kinematic statement. The biquaternion kinematic equations of motion of a rigid body in normalized and unnormalized finite displacement biquaternions are used as the mathematical model of motion, and the dual orthogonal projections of the instantaneous velocity screw of the body motion onto the body coordinate axes are used as the control. Various types of correction (stabilization), which are biquaternion analogs of position and integral corrections, are proposed. It is shown that the linear (obtained without linearization) and stationary biquaternion error equations that are invariant under any chosen programmed motion of the reference coordinate system can be obtained for the proposed types of correction and the use of unnormalized finite displacement biquaternions and four-dimensional dual controls allows one to construct globally regular control laws. The general solution of the error equation is constructed, and conditions for asymptotic stability of the programmed motion are obtained. The constructed theory of kinematic control of motion is used to solve inverse problems of robot-manipulator kinematics. The control problem under study is a generalization of the kinematic problem [1, 2] of reducing the body-attached coordinate system to the reference coordinate system rotating at a given (programmed) absolute angular velocity, and the presentedmethod for solving inverse problems of robotmanipulator kinematics is a development of the method proposed in [3–5].     

Geometric methods in determining rigid-body dynamics

This research develops a measurement system using linear accelerometers to determine the three-dimensional, six degrees of freedom, impact response of an anthropomorphic test device (dummy). A procedure using spherical geometric analysis (SGA) was developed. It uses three triaxial accelerometer clusters for determining angular velocity, angular acceleration, and linear acceleration. SGA differs in its calculation of angular velocity from other procedures which determine rigid-body motion. Unlike procedures which use linear accelerometers to determine angular velocity by integration of angular acceleration, SGA uses the topology of the sphere to obtain both angular acceleration and angular velocity through algebraic manipulation of the output from the linear accelerations. The validation of SGA is accomplished by the use of hypothetical as well as experimental data.     

Passenger air-bag study using geometric analysis of rigid-body motion

The goal of this research was to develop a system for measuring the force on the chest of an anthropomorphic test device (dummy) as it penetrates into an inflated-air-bag system during an impact. A spherical geometric analysis (SGA) method was developed to track dummy motion when film documentation was not reasonable, such as during air-bag deployment. The method uses three linear accelero-meter triaxial clusters for assessment of angular velocity and linear acceleration in a three-dimensional sense. From these measurements, angular acceleration and position, and linear velocity and position are derived.     

A Simple Method for the Determination of Angular Velocity and Acceleration of a Spherical Motion Through Quaternions

Given a quaternion representation of a spherical motion of a rigid body with respect to another body, acting as a reference frame, this contribution presents a simple and straightforward method for determining both the angular velocity and angular acceleration of the moving body with respect to the reference frame. Instead of employing orthogonal matrices or their linear invariants, this contribution makes use of quaternions avoiding, in this way, the series of matrix identities or theorems that are required in a pair of previous approaches.     

On the formulation of the planar ANCF triangular finite elements

In this paper, new planar isoparametric triangular finite elements (FE) based on the absolute nodal coordinate formulation (ANCF) are developed. The proposed ANCF elements have six coordinates per node: two position coordinates that define the absolute position vector of the node and four gradient coordinates that define vectors tangent to coordinate lines (parameters) at the same node. To shed light on the importance of the element geometry and to facilitate the development of some of the new elements presented in this paper, two different parametric definitions of the gradient vectors are used. The first parametrization, called area parameterization, is based on coordinate lines along the sides of the element in the reference configuration, while the second parameterization, called Cartesian parameterization, employs coordinate lines defined along the axes of the structure (body) coordinate system. The fundamental differences between the ANCF parameterizations used in this investigation and the parametrizations used for conventional finite elements are highlighted. The Cartesian parameterization serves as a unique standard for the triangular FE assembly. To this end, a transformation matrix that defines the relationship between the area and the Cartesian parameterizations is introduced for each element in order to allow for the use of standard FE assembly procedure and define the structure (body) inertia and elastic forces. Using Bezier geometry and a linear mapping, cubic displacement fields of the new ANCF triangular elements are systematically developed. Specifically, two new ANCF triangular finite elements are developed in this investigation, namely four-node mixed-coordinate and three-node ANCF triangles. The performance of the proposed new ANCF elements is evaluated by comparison with the conventional linear and quadratic triangular elements as well as previously developed ANCF rectangular and triangular elements. The results obtained in this investigation show that in the case of small and large deformations as well as finite rotations, all the elements considered can produce correct results, which are in a good agreement if appropriate mesh sizes are used.     

Kinematic problem of optimal nonlinear stabilization of angular motion of a rigid body

The problem of optimal transfer of a rigid body to a prescribed trajectory of preset angular motion is considered in the nonlinear statement. (The control is the vector of absolute angular velocity of the rigid body.) The functional to be minimized is a mixed integral quadratic performance criterion characterizing the general energy expenditure on the control and deviations in the state coordinates.Pontryagin’s maximum principle is used to construct the general analytic solution of the problem in question which satisfies the necessary optimality condition and ensures the asymptotically stable transfer of the rigid body to any chosen trajectory of preset angular motion. It is shown that the obtained solution also satisfies Krasovskii’s optimal stabilization theorem.     

Geometry and differentiability requirements in multibody railroad vehicle dynamic formulations

The dynamic equations of multibody railroad vehicle systems can be formulated using different sets of generalized coordinates; examples of these sets of coordinates are the absolute Cartesian and trajectory coordinates. The absolute coordinate based formulations do not require introducing an intermediate track coordinate system since all the absolute coordinates are defined in the global system. On the other hand, when the trajectory coordinates are used, a track coordinate system that follows the motion of a body in the railroad vehicle system is introduced. This track coordinate system is defined by the track geometry and the distance traveled by the body along the track centerline. The configuration of the body with respect to the track coordinate system is defined using five coordinates; two translations and three Euler angles. In this paper, the formulations based on the absolute and trajectory coordinates are compared. It is shown that these two sets of coordinates require different degrees of differentiability and smoothness. When an elastic contact formulation is used to study the wheel/rail dynamic interaction, there are significant differences in the order of the derivatives required in both formulations. In fact, as demonstrated in this study, in the absence of a contact constraint formulation, higher order derivatives with respect to geometric parameters are still required when the equations are formulated using the trajectory coordinates. The formulation of the constraints used in the analysis of the wheel/rail contact is discussed and it is shown that when the absolute coordinates are used, only third order derivatives need to be evaluated. The relationship between the track frame used in railroad vehicle dynamics and the Frenet frame used in the theory of curves to describe the curve geometry is also discussed in this paper. Based on the analysis presented in this paper, the advantages and drawbacks of a hybrid method which employs both the absolute and trajectory coordinates and planar contact conditions in order to reduce the number of contact constraints and relax the differentiability requirements are discussed. In this method, the absolute coordinates are used to formulate the equations of motion of the railroad vehicle system. The absolute coordinate solution can be used to determine the trajectory coordinates and their time derivatives. Using the trajectory coordinates, the motion of the body in the vehicle with respect to the track coordinate system can be predicted and used in the formulation of the planar contact model.     

两种坐标系下惯导传递对准效果比较

分别建立了发射点惯性坐标系和当地地理坐标系下的惯导误差传播模型,以及典型匹配模式下的传递对准基本方程。基于运载火箭主动段飞行,分析了姿态机动运动对两种坐标系下传递对准的不同影响。数学模型和仿真分析表明,两种坐标系下的传递对准在原理和系统特性上均有较大差异。典型匹配方案的对准效果对比表明,发射点惯性坐标系内的"角速度+加速度"匹配是一种能够实现子惯导快速、准确初始化的有效方法。     

Rigid-body pose and twist estimation using an accelerometer array

Summary A novel technique for the determination of the pose and the twist of rigid bodies using point-acceleration data is proposed. These data are collected from an accelerometer array, which is a kinematically redundant set of triaxial accelerometers. Because orientational error in the installation of the accelerometers can be fatal to the accuracy of the results, a calibration procedure based on the consistency of the point accelerations is outlined. The formulation developed is then utilized in the simulation analysis of two sample motions. The relations required to estimate the pose and the twist are derived in a body-fixed frame. The body angular acceleration and angular velocity, in this order, are determined directly from the acceleration data; the body attitude is then computed through integration.This research work was supported by the Natural Sciences and Engineering Research Council of Canada, under Research Grants OGP0004532 and OGP0000967.     

On Conservation and Balance Laws in Micromorphic Elastodynamics

In this paper, following Noether’s theorem we investigate the Lie point symmetries of linear micromorphic elastodynamics (linear elastodynamics with microstructure). Conservation and balance laws of linear, micromorphic elastodynamics are derived. We generalize the J, L and M integrals for this theory. In addition, we give the Eshelby stress tensor, pseudomomentum vector, field intensity vector, Hamiltonian, angular momentum tensor and scaling flux generalized to micromorphic elastodynamics.      

Attitude maneuver of liquid-filled spacecraft with a flexible appendage by momentum wheel

Attitude maneuver of liquid-filled spacecraft with an appendage as a cantilever beam by momentum wheel is studied.The dynamic equations are derived by conservation of angular momentum and force equilibrium principle.A feedback control strategy of the momentum wheel is applied for the attitude maneuver.The residual nutation of the spacecraft in maneuver process changes with some chosen parameters,such as steady state time,locations of the liquid container and the appendage,and appendage parameters.The results indicate that locations in the second and fourth quadrants of the body-fixed coordinate system and the second quadrant of the wall of the main body are better choices for placing the liquid containers and the appendage than other locations if they can be placed randomly.Higher density and thicker cross section are better for lowering the residual nutation if they can be changed.Light appendage can be modeled as a rigid body,which results in a larger residual nutation than a flexible model though.The residual nutation decreases with increasing absolute value of the initial sloshing angular height.     

Construction of optimal laws of variation of the angular momentum vector of a rigid body

We consider the problem of constructing optimal preset laws of variation of the angular momentum vector of a rigid body taking the body from an arbitrary initial angular position to the required terminal angular position in a given time. We minimize an integral quadratic performance functional whose integrand is a weighted sum of squared projections of the angular momentum vector of the rigid body. We use the Pontryagin maximum principle to derive necessary optimality conditions. In the case of a spherically symmetric rigid body, the problem has a well-known analytic solution. In the case where the body has a dynamic symmetry axis, the obtained boundary value optimization problem is reduced to a system of two nonlinear algebraic equations. For a rigid body with an arbitrarymass distribution, optimal control laws are obtained in the form of elliptic functions. We discuss the laws of controlled motion and applications of the constructed preset laws in systems of attitude control by external control torques or rotating flywheels.     

Lumped parameter analytic modeling and behavioral simulation of a 3-DOF MEMS gyro-accelerometer

A new analytical model of a 3-degree-of-freedom(3-DOF) gyro-accelerometer system consisting of a 1-DOF drive and 2-DOF sense modes is presented. The model constructs lumped differential equations associated with each DOF of the system by vector analysis. The coupled differential equations thus established are solved analytically for their responses in both the time and frequency domains. Considering these frequency response equations, novel device design concepts are derived by forcing the sense phase to zero, which leads to a certain relationship between the structural frequencies, thereby causing minimization of the damping effect on the performance of the system. Furthermore, the feasibility of the present gyro-accelerometer structure is studied using a unique discriminatory scheme for the detection of both gyro action and linear acceleration at their events. This scheme combines the formulated settled transient solution of the gyro-accelerometer with the processes of synchronous demodulation and filtration, which leads to the in-phase and quadrature components of the system's output signal. These two components can be utilized in the detection of angular motion and linear acceleration. The obtained analytical results are validated by simulation in a MATLAB/Simulink environment, and it is found that the results are in excellent agreement with each other.     

落体法测刚体转动惯量的测量方法比较

用JIJG—I刚体转动惯量实验仪的同组数据,分别用测角加速度法和测时间法计算待测圆盘的转动惯量,结果表明,测角加速度法比测时间法的测量精度可提高1到2个数量级,具有明显优点;同时分析了测时间法产生误差的主要原因,也指出了使用测角加速度法的注意事项．     

动态环境下光纤陀螺误差辨识与补偿技术

为了提高光纤陀螺在高动态环境下的测量精度,需要精确地辨识角加速度信息以便有效地补偿。针对直接对陀螺的角速度信息微分处理后得到角加速度的方法误差较大的问题,提出了将微分后的角加速度信息分为线性和非线性两个部分,其中线性部分采用Savitzky-golay最小二乘拟合,而非线性部分则采用RBF神经网络技术进行拟合。上述处理方法能更真实地反映实际物理过程,具有较强的自适应性和较好的拟合效果。通过试验验证,证明了该方法的有效性和准确性,提高了角加速度辨识精度,比直接微分的方法测量精度提高二个数量级,有效地补偿了陀螺仪在高动态环境下的测量精度。     

GENERALIZATION OF KIRCHHOFF KINETIC ANALOGY TO THIN ELASTIC SHELLS 1)

将弹性细杆的"Kirchhoff动力学比拟"方法推广到弹性薄壳,使弹性薄壳的变形在物理概念上和刚体的运动对应, 在数学表述上等同,从而可以用刚体动力学的理论和方法研究弹性薄壳的变形,为连续的弹性薄壳提供新的离散化方法. 在直法线假设下,在弹性中面上构筑空间正交轴系, 此轴系沿坐标线"运动"的角速度构成两自变量的弯扭度. 沿两个坐标线的弯扭度表达了弹性薄壳的变形和位形,证明了弯扭度之间以及弯扭度与中面切矢间的相容关系. 用Euler角和Lam$\acute{e}$系数表达了非完整约束和中面位形的微分方程,用弯扭度和Lam$\acute{e}$系数表达了应变和应力以及内力及其本构方程.导出了用分布内力集度表达的弹性薄壳在变形后位形上的平衡偏微分方程组,方程的形式与刚体动力学的Euler方程和弹性细杆的Kirchhoff方程具有相似性,实现了Kirchhoff动力学比拟对弹性薄壳的推广.总结了弹性薄壳静力学和刚体动力学以及弹性细杆静力学在概念上的比拟关系.最后给出了一个算例. 为研究弹性薄壳的变形和运动提供新的建模方法和研究思路.也可进一步推广到弹性薄壳动力学.