Resolving Controversies in the Application of the Method of Multiple Scales and the Generalized Method of Averaging

I compare application of the method of multiple scales with reconstitution and the generalized method of averaging for determining higher-order approximations of three single-degree-of-freedom systems and a two-degree-of-freedom system. Three implementations of the method of multiple scales are considered, namely, application of the method to the system equations expressed as second-order equations, as first-order equations, and in complex-variable form. I show that all of these methods produce the same modulation equations.I address the problem of determining higher-order approximate solutions of the Duffing equation in the case of primary resonance. I show that the conclusions of Rahman and Burton that the method of multiple scales, the generalized method of averaging, and Lie series and transforms might lead to incorrect results, in that spurious solutions occur and the obtained frequency–response curves bear little resemblance to the actual response, is the result of their using parameter values for which the neglected terms are the same order as the retained terms. I show also that spurious solutions cannot be avoided, in general, in any consistent expansion and their presence does not constitute a limitation of the methods. In particular, I show that, for the Duffing equation, the second-order frequency–response equation does not possess spurious solutions for the case of hardening nonlinearity, but possesses spurious solutions for the case of softening nonlinearity. For sufficiently small nonlinearity, the spurious solutions are far removed from the actual response. But as the strength of the nonlinearity increases, these solutions move closer to the backbone and eventually distort it. This is not a drawback of the perturbation methods but an indication of an application of the analysis for parameter values outside the range of validity of the expansion.Also, I address the problem of obtaining non-Hamiltonian modulation equations in the application of the method of multiple scales to multi-degree-of-freedom Hamiltonian systems written as second-order equations in time and how this problem can be overcome by attacking the state-space form of the governing equations. Moreover, I show that application of a variation of the method of Rahman and Burton to multi-degree-of-freedom systems leads to results that do not agree with those obtained with the generalized method of averaging.Contributed by Prof. R.A. Ibrahim.     

基于Gauss伪谱法的3D刚体摆姿态最优控制

研究Gauss伪谱法求解3D刚体摆姿态最优控制问题．针对其最优姿态控制问题,既要满足由任意位置运动到平衡位置姿态运动规划问题,又要满足系统含有动力学约束的力学模型问题,提出基于四元数来描述3D刚体摆的数学模型,建立3D刚体摆姿态的动力学和运动学方程,为了解决3D刚体摆在平衡位置处的姿态最优控制问题,设计基于Gauss伪谱算法的最优姿态开环控制器,得到了3D刚体摆的姿态最优控制轨迹,得到满足的可行解,通过仿真实验验证了其开环解在平衡位置的控制姿态最优性．     

Complex Inner Product Averaging Method for Calculating Normal Form of Ode

ＩｎｔｒｏｄｕｃｔｉｏｎＭｏｓｔｎｏｎｌｉｎｅａｒｄｙｎａｍｉｃｓｙｓｔｅｍｓａｒｅｍｕｌｔｉ＿ｄｅｇｒｅｅ＿ｏｆ＿ｆｒｅｅｄｏｍｓｙｓｔｅｍｓｉｎｔｈｅｆｉｅｌｄｏｆｅｎｇｉｎｅｅｒｉｎｇ[1,2 ],ｎｏｎｌｉｎｅａｒｄｙｎａｍｉｃｓｙｓｔｅｍｓｄｏｎｏｔｈａｖｅｅｘａｃｔｌｙｔｈｅｏｒｅｔｉｃｓｏｌｕｔｉｏｎｓ.Ｔｈｅｒｅｆｏｒｅ,ｉｔｉｓａｆｏｃａｌｐｏｉｎｔｃｏｎｃｅｒｎｅｄｂｙｓｃｉｅｎｔｉｓｔｓｔｏｓｅａｒｃｈｆｏｒａｐｐｒｏｘｉｍａｔｅｓｏｌｕｔｉｏｎｓｆｏｒａｌｏｎｇｔｉｍｅ ,ａｎｄｒ…     

The Equations of Motion of a Rigid Body with Constraints Explicitly Depending on Time

The general equations of motion of a body for an observer in S ₀ (an inertial frame) or for an observer in S ₁ (an accelerated frame) are derived. They allow us to determine, in any case, the inertial and gyroscopic forces and to find the difference between them. If the vector that determines the position of the body in S ₀ depends explicitly on time in S ₀, the work-energy principle yields a supplementary condition and these equations can be shown to be equivalent to the Painleve integrals in the Lagrange formulation. Since we are dealing with inertial frames, gyroscopic forces rather than inertial forces are taken into account. If another reference frame is used, we can choose it so that the position vector depends implicitly on time in S ₁ and another set of equations can be obtained for the motion of the body in S ₁. The work-energy principle yields a supplementary condition, but inertial forces should be added. Since an explicit time dependence does not exist in S ₁, gyroscopic forces do not exist as well and instead we have Coriolis forces that behave like gyroscopic forces     

Exact Averaging of Stochastic Equations for Transport in Random Velocity Field

We present new examples of exactly averaged multi-dimensional equation of transport of a conservative solute in a time-dependent random flow velocity field. The functional approach and a technique for decoupling the correlations are used. In general, the averaged equation is non-local. We study the special cases where the averaged equation can be localized and reduced to a differential equation of finite-order, where the problem of evolution of the initial plume (Cauchy problem) can be solved exactly. We present in detail the results of the analyses of two cases of exactly averaged problems for Gaussian and telegraph random velocity with an identical exponential correlation function, which are informative and convenient models for continuous and discontinuous random functions. The problems in which the field has sources of solute and boundaries are also examined. We study the behavior of different initial plumes for all times (evolutions and convergence) and show the manner in which they approach the same asymptotic limit for two stochastic distributions of flow-velocity. A comparison between exact solutions and solutions derived by the method of perturbation is also discussed.     

The Equations of Motion of a Plate with Residual Stress

We consider a three-dimensional body which is at rest in a cylindrical configuration of height 2ε. We assume that the material is residually stressed and that it responds elastically to deformations from the reference configuration. Under appropriate assumptions on the data, and using weak-convergence methods, we determine the limit, as ε goes to zero, of the elasto-dynamic problem. The plate problem obtained, as in the case without residual stress, splits into two problems: one governing the in-plane displacement and the other determining the out-of-plane motion.     

控制一类磁性刚体航天器的混沌姿态运动

本文研究一类磁性航天器的混沌姿态运动及其控制,建立了在近地球赤道面圆轨道上运动受万有引力矩、磁力矩作用磁性刚体航天器姿态运动的动力学方程。采用时间历程、Poincare截面、Lyapunov指数和功率谱对系统的动力学行为进行数值识别,结果表明随着磁场参数的增大系统动力学行为由准周期环面破裂而出现混沌。利用输入－输出反馈精确线性化的方法将航天器的混沌姿态控制运动控制为姿态静止和按给定的周期规律运动,数值结果表明该控制方法的有效性。     

Averaging principle for quasi-geostrophic motion under rapidly oscillating forcing

1　IntroductionandtheProblemPresentedClimateandgeophysicalflowsdefinetheenvironmentinwhichwelive ,andtheirunderstanding,beyondbeingafascinatingscientificchallenge,isoftremendouseconomicandsocialimportance.Aclassoflargescalegeophysicalfluidflows,suchastheGulfStreamandoceanicgyres,aremodelledbythequasi_geostrophicequation .Thequasi_geostrophic (QG)equationmodelsaclassoflargescalegeophysicalflows.ItisderivedasanapproximationoftherotatingNavier_StokesequationsbyanasymptoticexpansioninasmallRos…     

A Technique for Identification of The Principal Central Axis of Inertia in an Inhomogeneous Rigid Body

A technique is proposed to identify the principal central axis of inertia in a joint-suspended inhomogeneous rigid body from its vertical revolution. The problem is reduced to synchronizing the revolution of the body with the revolution of a fictitious mathematical pendulum     

Generalized Method of Separation of Variables on the Nonlinear Elastic-Plastic Wave Motion

The method of separation of variables is developed in addition to the Karman method and the method of characteristics for the wave motion of uniaxial stress in rods. Rate-independent theory is considered and it is shown that the plastic wave speed is independent of the constant n in the stress-strain law. Stresses are determined in explicit forms for two cases.     

流固偶合运动的边界元法

本文给出了流固偶合运动(包括物体散射辐射及偶合运动)的边界元法理论和应用.对于散射问题,求出了物体引起的散射势及入射波作用于物体的载荷.对于辐射问题,求出了辐射势及物体在流体中运动的附加质量和附加阻尼.偶合问题包括求其中包含的散射势和辐射势以及作用于物体之上的散射力、物体的附加质量、附加阻尼、物体在入射波作用下的运动.在偶合运动问题中,本文采取了边界积分方程与物体在流体中的运动方程联立求解的方法,并将其运用到边界元法的数值过程中.所编制的程序有较高的精度.最后给出了数值计算结果与理论解的比较.     

Geometrical Representation of the Motion of a Body Through a Medium

The paper presents a complete qualitative analysis of the model plane-parallel motion of a body (plate) through a resisting medium with jet or detached flow. A simplified system is analyzed on the phase plane. A geometrical interpretation is given to the motion of the plate     

空间刚架结构的广义逆力法

在结构计算中，根据算法中所采用的基本未知量的不同．结构分析方法可以分为力法、位移法和混合法。其中位移法由于适宜计算机处理而在结构计算领域得到了广泛的应用．经典力法相比之下应用就远不如位移法普遍，虽然力法本身在力学上有其独特的优势。广义逆矩阵做为一种较新的数学工具，自二十世纪五十年代诞生以来正日益表现出越来越旺盛的生命力。广义逆力法就是一种基于力法和广义逆矩阵理论的新的迭代算法。这种算法是一种完全适合计算机处理的力法方法。该算法的思路以及对于求解线弹性空间刚架结构问题的具体公式均在文中给出并给出了算例。从算例计算结果可以看到广义逆力法有着较好的计算效率和计算精度。该算法的提出为力法在计算机计算领域的应用开拓了新的发展空间。该算法在材料非线性问题和结构并行计算方面也有着较好的应用前景。     

弹性介质中几何和运动非线性问题的类孤波解

将弹性介质 的几何和运动非线性方程简化成具有电磁场中的Ｂｏｒｎ－Ｉｎｆｅｌｄ方程的形式，并证明了该方程的类孤波解的存在．     

基于卡尔曼滤波算法的汽车运动参数测试方法研究

提出了一种采用多维组合惯性测量元件进行汽车运动参数测试的方法,并利用卡尔曼滤波算法提高了测试系统的精度。仿真结果表明：卡尔曼滤波算法对汽车姿态、速度等的解算具有较好的效果。     

The control of a permanent rotational motion and equilibrium position of a spacecraft

This article has adopted an analytical method to obtain a non-linear control law to reach the exponential asymptotic stablity of the permanent rotational motion of a spacecraft. The control moments achieving this rotational motion are obtained. The control moments to establish exponential asymptotic stablity of the mentioned motion are obtained as non-linear functions of the phase coordinates of the spacecraft. The general solution of the equations of perturbed motion is derived. Furthermore, analysis and numerical simulation study of this solution are presented. For numerical examples the time needed for control is calculated. An equilibrium position of the spacecraft is proved to be exponentially asymptotically stable as a special case of the above-studied problem.     

The Dynamics of a Detaching Rigid Body

Consider a rigid block perfectly bonded to an elastic half-space. Suppose the block is impulsively struck so that it starts to oscillate with respect to its initial position. In the absence of dissipation, the motion will be harmonic. However, it may happen that the displacement becomes so large that the bond between the block and the half-space will begin to rupture. In this case, the subsequent motion of the rigid body is described by a second-order nonlinear ordinary differential equation, derived from the principle of conservation of mechanical energy. There is a critical bond length, with respect to a pre-assigned length scale. Below it the block will undergo harmonic motion with no rupture at low impact or eventual catastrophic rupture and instantaneous detachment at high impact. Above it the block will also undergo harmonic motion with no rupture at low impact. At very high impact, the block will undergo a period of controlled rupture ending with catastrophic rupture and instantaneous detachment. There is also an intermediate range of impacts which result in a period of controlled partial rupture after which rupture ceases and periodic motion ensues.     

中心有刚体质量的环形薄板的非线性强迫振动

本文研究中心带有刚体质量外部固定铰支或活动铰支的环形薄板的非线性强迫振动。考虑板的弯曲变形、面内位移和几何非线性，用哈米尔顿原理建立板的运动方程，用Ｋａｎｔｏｒｏｖｉｃｈ平均法消去时间变量，然后用数值积分求得非线性振动的振幅随激振力的大小及激振力的频率而变化的关系。求解过程中用打靶法逐步改进未知参量，以保证边界条件的满足。最后讨论薄膜力、激振力的分布、板的内外半径比等因素对响应的影响     

Equations of motion for mechanical systems: A unified approach

This paper presents a unified framework from which emerge the Lagrange equations, the Gibbs-Appell Equations and the Generalized Inverse Equations for describing the motion of constrained mechanical systems. The unified approach extends the applicability of the first two approaches to systems where the constraints are non-linear functions of the generalized velocities and are not necessarily independent. Furthermore, the approach leads to the Explicit Gibbs-Appell Equations.