REVIEW OF DYNAMIC MODELING AND STABILITY ANALYSIS OF A BICYCLE$^{\bf 1)}$

自行车发明于两个多世纪前. 这一看似古老的交通工具在为人们提供出行便利的同时,其独特的运动特性及动力学性质 也吸引了来自数学、物理及力学等多个学科相关学者的兴趣. 大体上,自行车可以描述为具有 7 个自由度和 4 个非完整约束的多刚体系统. 但由于前后车轮之间复杂的运动耦合关系,使得自行车的约束方程和动力学模型变得异常复杂, 导致对自行车的稳定性存在一些模糊认识. 本文针对经典的 Carvallo-Whipple 自行车构型,系统回顾了历史上自行车动力学研究中的相关问题,这些问题包括：(1) 自行车在复杂曲面上的几何约束和非完整约束的数学描述;(2) 自行车系统内在的对称性及守恒量; (3) 自行车动力学的各类建模方法; (4) 自行车运动的相对平衡点及稳定性分析,包括水平面上的匀速直线运动及旋转对称曲面上的匀速圆周运动;(5) 影响自行车自稳定性的结构参数等. 本文最后对自行车动力学实验和控制方面的研究工作进行了回顾,并对自行车今后的研究给出了展望.     

非完整约束系统几何动力学研究进展:Lagrange理论及其它

近10年来, 非完整力学的发展主要集中在两个相互关联的方向上, 一个是非完整运动规划, 另一个则是非完整约束系统的几何动力学, 这两个研究方向都充分地利用了现代几何学, 如纤维丛理论、辛流形和Poisson流形结构等等.本文主要综述非完整约束系统几何动力学的外附型和内禀型Lagrange理论, 包括非定常力学系统所需要的射丛几何学的基本概念、射丛按约束的直和分解、约束流形上的水平分布、D'Alembert-Lagrange方程与Chaplygin方程的整体描述、以及Riemann-Cartan流形上的非完整力学, 文中对Chetaev条件和d-δ交换关系的几何意义作了深入讨论.除此之外, 简要评述非完整力学的Hamilton理论与赝Poisson结构、Noether对称性和Lie对称性、动量映射与对称约化、Vakonomic动力学等几个非常重要专题的研究进展.      

动力学与控制论力学系统的自由度

力学系统的自由度定义源自描述系统位形的独立坐标数.在分析力学发展过程中,人们通过对非完整约束的研究,将其拓展为独立的坐标 变分数.本文指出,对于含非完整约束的力学系统,该定义存在不妥之处,给出的自由度会过度限制系统的力学行为.文中研究力学系统在状态空间中的可达流形,指出可达流形维数与描述系统动力学的一阶常微分方程组的最少未知函数个数一致,例如Gibbs-Appell方程与广义速度方程联立的未知函数个数,进而将可达流形维数的一半定义为系统自由度.通过含黏弹性支承的振动系统、在倾斜平面上运动的冰橇等案例,讨论了单个非完整约束导致的半自由度概念,指出其力学意义和与相邻整数自由度的关系.此外,文中还给出两个非完整约束导致系统减少一个自由度的案例,讨论了系统的切丛和余切丛维数.      

具有单面非线性非完整约束的动力学系统的运动

本文提出并研究了具有单面非线性非完整约束的动力学系统的运动问题,给出了描述系统运动的Ｂｏｌｔｚｍａｎｎ－Ｈａｍｅｌ型方程,并举例说明结果的应用。     

论力学系统的自由度

力学系统的自由度定义源自描述系统位形的独立坐标数.在分析力学发展过程中,人们通过对非完整约束的研究,将其拓展为独立的坐标变分数.本文指出,对于含非完整约束的力学系统,该定义存在不妥之处,给出的自由度会过度限制系统的力学行为.文中研究力学系统在状态空间中的可达流形,指出可达流形维数与描述系统动力学的一阶常微分方程组的最少未知函数个数一致,例如Gibbs-Appell方程与广义速度方程联立的未知函数个数,进而将可达流形维数的一半定义为系统自由度.通过含黏弹性支承的振动系统、在倾斜平面上运动的冰橇等案例,讨论了单个非完整约束导致的半自由度概念,指出其力学意义和与相邻整数自由度的关系.此外,文中还给出两个非完整约束导致系统减少一个自由度的案例,讨论了系统的切丛和余切丛维数.     

一类分段光滑非线性平面运动系统响应特性研究

针对由一个线性子系统和一个非线性子系统构成的两自由度非自治分段光滑非线性平面运动系统的响应特性开展了研究. 该分段光滑非线性模型可用来确定对称转子/定子系统的主要碰摩响应, 且在反映非光滑系统典型特性上具有明显的特征:(1) 切换分界面是由两自由度坐标共同决定的一个幅值曲面;(2) 子系统周期解与分界面的擦碰, 不是发生在一个点上, 而是同时发生在解的所有点上;(3) 完整系统未发现由两子系统共同作用而产生的周期解. 因此, 对于该非光滑系统响应特性的研究, 很难直接利用目前有关非光滑系统平衡点和周期解分岔分析的方法. 为此, 尝试了根据子系统的响应特性, 划分出完整系统响应对分界面处切换的敏感区和非敏感区, 并针对非敏感区可由子系统解的特性求得完整系统的响应, 而针对敏感区通过子系统动力学特征的分析有助于解释完整系统响应的生成机制.     

一类多体系统非完整运动最优规划的拟牛顿算法

研究带有非完整约束的一类多体系统运动规划问题。多体系统中的非完整约束通常是由不可积的速度约束或不可积的守恒律引起。在系统动量和动量矩守恒情况下,动力学方程降阶为非完整形式约束方程,系统的控制问题可转化为无漂移系统的非完整运动规划问题。文中首先导出具有多体开链系统的非完整运动模型。利用最优控制理论和最优化技术,采用输入参数化的方法将连续的最优控制问题转化为离散的最优控制问题,提出一种非完整多体系统运动规划的拟牛顿算法。最后将该方法用于自由漂浮的空间三连杆机构,仿真结果验证了该方法的有效性。     

变质量可控力学系统的Gauss原理和Appell方程

力学系统的运动依赖于力和约束,人们可以借助于力来控制运动(称为动力学控制),也可以借助于约束来控制运动(称为运动学控制)。我们研究一类力学系统,它的约束依赖于某些控制参数。得到可控力学系统的Lagrange方程、Hamilton方程和AppelI方程,用Gauss原理导出一阶非线性非完整系统广义坐标下的Appell 方程。本文考虑非线性非完整系统,导出了变质量可控力学系统在广义坐标和准坐标下的Gauss原理及Appell方程,最后举例说明其应用。本文主要结果是(1.5),(1.6),(1.20),(2.1),(3.13)及(3.14)。     

非完整系统稳定性的若干进展

介绍非完整系统稳定性理论的某些近代进展,包括非完整系统平衡位置关于全部变量的稳定性,平衡状态流形的稳定性,平衡位置关于部分变量稳定性及其与关于全部变量的稳定性的关系,平稳运动的稳定性,以及非完整控制系统的镇定．同时,讨论非完整系统稳定性的几个主要应用,并给出几个未来研究方向的建议      

THE GENERALITY OF THE GENERAL EQUATIONS OF DYNAMICS 1)

理论力学中动力学普遍方程,在分析力学中称为d'Alembert--Lagrange原理。动力学普遍方程之普遍在于,由它不仅可导出动力学普遍定理,可导出完整约束系统和非完整约束系统的运动微分方程,还可导出积分变分原理。     

An approximate decoupled dynamics and kinematics analysis of legless locomotion

We present a novel analysis technique to understand the dynamics of a recently described locomotion mode called legless locomotion. Legless locomotion is a locomotion mode available to a legged robot when it becomes high-centered, that is, when its legs do not touch the ground. Under these conditions, the robot may still locomote in the plane by swinging its legs in the air, rocking on its body, and taking advantage of the nonholonomic contact constraints. Legless locomotion is unique from all previously studied locomotion modes, since it combines the effect of oscillations due to controls and gravity, nonholonomic contact constraints, and a configuration-dependent inertia. This complex interaction of phenomena makes dynamics analysis and motion planning difficult, and our proposed analysis technique simplifies the problem by decoupling the robot’s oscillatory rotational dynamics from its contact kinematics and also decoupling the dynamics along each axis. We show that the decoupled dynamics models are significantly simpler, provide a good approximation of the motion, and offer insight into the robot’s dynamics. Finally, we show how the decoupled models help in motion planning for legless locomotion.     

Forward and inverse dynamics of nonholonomic mechanical systems

This paper deals with the forward and the inverse dynamic problems of mechanical systems subjected to nonholonomic constraints. The intrinsically dual nature of these two problems is identified and utilised to develop a systematic approach to formulate and solve them according to an unified framework. The proposed methodology is based on the fundamental equations of constrained motion which derive from Gauss’s principle of least constraint. The main advantage arising from using the fundamental equations of constrained motion is that they represent an effective method capable to derive the generalised acceleration of a mechanical system, constrained in general by a set of nonholonomic constraints, together with the generalized constraint forces (forward dynamics). When the constraint equations are used to represent the desired behaviour of the mechanical system under study, the generalised constraint forces deriving from the fundamental equations of constrained motion provide the control actions which reproduce the specified motion for the system (inverse dynamics). This approach is systematically extended to underactuated mechanical systems introducing a new method named underactuation equivalence principle. The underactuation equivalence principle is founded on the key idea that the underactuation property of a mechanical system can be mathematically represented using a particular set of nonholonomic constraint equations. Two simple case-studies are reported to exemplify the proposed methodology. In the first case-study the computation of the generalised constraint forces relative to the revolute joint constraints of a physical pendulum is illustrated. In the second case-study the calculation of the control action which solves the swing-up problem for an inverted pendulum is described.     

爬壁机器人系统的Noether对称性和守恒量

爬壁机器人的运动是一种模仿壁虎爬行的运动, 爬壁机器人的运动可分解为四肢带动身体的运动, 先前的研究都是基于牛顿力学的方法. 本文采用Lagrange 力学的方法建立爬壁机器人系统的运动方程, 并运用Lie群分析方法建立该系统的Noether对称性理论, 得出爬壁机器人的运动规律. 首先, 给出非完整爬壁机器人系统的动能、势能和Lagrange函数以及所受的非完整约束, 从而建立了非完整爬壁机器人系统的Lagrange方程; 其次, 引入关于时间和广义坐标的无限小变换, 提出了非完整爬壁机器人系统的Hamilton作用量和Hamilton作用量的基本变分公式; 第三, 给出爬壁机器人系统 Noether对称性变换和广义准对称变换的定义, 判据和存在的Noether守恒量, 并提出了非保守完整系统和非保守非完整爬壁机器人系统的Noether定理; 最后, 以圆锥面上爬壁机器人为例, 对给出的守恒量直接进行积分给出圆锥面上爬壁机器人整体运动的精确解和四肢运动的数值解, 发现了该爬壁机器人的运动规律, 很好地验证了非完整爬壁机器人系统的Noether对称性理论. 本文的研究为Lie群分析方法应用于其他复杂的机器人系统以及柔性机器人系统的对称性求解提出了一种新的对称性求解方法.      

Trajectory planning and tracking control for autonomous bicycle robot

Control of the autonomous bicycle robot offers considerable challenges to the field of robotics due to its nonholonomic, underactuated, and nonminimum-phase properties. Furthermore, instability and complex dynamic coupling make the trajectory planning of the bicycle robot even more challenging. In this paper, we consider both trajectory planning and tracking control of the autonomous bicycle robot. The desired motion trajectory of the contact point of the bicycle’s rear wheel is constructed using the parameterized polynomial curve that can connect two given endpoints with associated tangent angles. The parameters of the polynomial curve are determined by minimizing the maximum of the desired roll angle’s equilibrium of the bicycle, and this optimization problem is solved by the particle swarm optimization algorithm. Then, a control scheme that can achieve full-state trajectory tracking while maintaining the bicycle’s balance is proposed by combining a planar trajectory tracking controller with a roll angle balance controller. Simulation results are presented to demonstrate the effectiveness of the proposed method.     

Direct sensitivity analysis of multibody systems with holonomic and nonholonomic constraints via an index-3 augmented Lagrangian formulation with projections

Optimizing the dynamic response of mechanical systems is often a necessary step during the early stages of product development cycle. This is a complex problem that requires to carry out the sensitivity analysis of the system dynamics equations if gradient-based optimization tools are used. These dynamics equations are often expressed as a highly nonlinear system of ordinary differential equations or differential-algebraic equations, if a dependent set of generalized coordinates with its corresponding kinematic constraints is used to describe the motion. Two main techniques are currently available to perform the sensitivity analysis of a multibody system, namely the direct differentiation and the adjoint variable methods. In this paper, we derive the equations that correspond to the direct sensitivity analysis of the index-3 augmented Lagrangian formulation with velocity and acceleration projections. Mechanical systems with both holonomic and nonholonomic constraints are considered. The evaluation of the system sensitivities requires the solution of a tangent linear model that corresponds to the Newton–Raphson iterative solution of the dynamics at configuration level, plus two additional nonlinear systems of equations for the velocity and acceleration projections. The method was validated in the sensitivity analysis of a set of examples, including a five-bar linkage with spring elements, which had been used in the literature as benchmark problem for similar multibody dynamics formulations, a point-mass system subjected to nonholonomic constraints, and a full-scale vehicle model.     

轮式移动机器人非完整运动规划的最优控制方法

本文讨论轮式动机器人非完整运动的最优规划问题,利用约束与最优控制理论建立数学模型,考虑系统的非完整约束特性,提出轮式移动机器人运动规划的最优控制算法。通过数值仿真,表明该方法的有效性。     

Dynamics of a small bubble in a compressible fluid

An efficient (simplified) method for solving problems of spherically symmetric dynamics of a small gas bubble in a compressible fluid is proposed. The method is based on the joint use of the full problem statement (the gas dynamics equations for the gas and the fluid) and its relevant simplifications. Some approximate statements are discussed. In the proposed method, the rarefaction and compression of the gas during the slow motion of the bubble surface is assumed to be uniform over the bubble volume. At the same time the fluid in the thin zone adjacent to the bubble is considered to be slightly compressible. Otherwise the gas dynamics equations are used for the gas and the fluid. The dynamics of the fluid in the thick external zone are described by the linear acoustics only. The proposed simplified method and two others used in literature are estimated by comparison of their numerical results with those obtained in full statement. Copyright © 2000 John Wiley & Sons, Ltd.     

HYBRID OPTIMIZATION STRATEGY FOR MOTION PLANNING OF SPACE ROBOT SYSTEM WITH PRISMATIC JOINT

研究了自由漂浮带滑移铰空间机器人非完整运动规划的最优控制问题,提出一种由高斯伪谱法求解可行解与直接打靶法求解最优解相结合的混合优化策略.首先,根据多体系统动力学理论建立空间机器人的动力学模型,给定系统的初始和目标位形,将空间机器人运动规划问题描述成博尔察（Bolza）型最优控制问题;然后,利用高斯伪谱法将最优控制问题离散为非线性规划问题,求解在较少勒让德-高斯（Legendre-Gauss,LG）点时状态变量和控制变量对应的可行解;最后,在LG点处离散控制变量,作为直接打靶法的初值,利用序列二次规划算法求解空间机器人系统的优化运动轨迹和最优控制输入.通过数值仿真,系统优化运动轨迹光滑平稳,最优控制输入也能很好地满足各种约束条件,仿真结果验证了该混合优化策略的鲁棒性和有效性.