Effect of damping and wave parameters on offshore structure under random excitation

This paper describes the linearized and nonlinear dynamic response of a tension leg platform (TLP) to random waves and current forces. The forcing term of the equation of motion is inherently nonlinear due to the nonlinear drag force. Two analysis procedures are used: nonlinear time domain analysis and linear frequency domain analysis. For the nonlinear analysis, the random wave particle velocities and accelerations are simulated for a given wave spectrum. The nonlinear equation of motion is then integrated directly to obtain the system response statistics. For the linear frequency domain analysis, the nonlinear drag force is linearized through an introduction of linearization coefficients. The main objective of this paper is to investigate the effect of the structural damping and wave parameters on both nonlinear and linear dynamic response of the TLP by parametric studies. The results of stochastic nonlinear and linear dynamic response of the TLP, with and without the presence of current, are presented and compared.     

Effect of the deformation in the inertia forces on the inverse dynamics of planar flexible mechanical systems

The inverse dynamics problem for articulated structural systems such as robotic manipulators is the problem of the determination of the joint actuator forces and motor torques such that the system components follow specified motion trajectories. In many of the previous investigations, the open loop control law was established using an inverse dynamics procedure in which the centrifugal and Coriolis inertia forces are linearized such that these forces in the flexible model are the same as those in the rigid body model. In some other investigations, the effect of the nonlinear centrifugal and Coriolis forces is neglected in the analysis and control system design of articulated structural systems. It is the objective of this investigation to study the effect of the linearization of the centrifugal and Coriolis forces on the nonlinear dynamics of constrained flexible mechanical systems. The virtual work of the inertia forces is used to define the complete nonlinear centrifugal and Coriolis force model. This nonlinear model that depends on the rate of the finite rotation and the elastic deformation of the deformable bodies is used to obtain the solution of the inverse dynamics problem, thus defining the joint torques that produce the desired motion trajectories. The effect of the linearization of the mass matrix as well as the centrifugal and Coriolis forces on the obtained feedforward control law is examined numerically. The results presented in this investigation are obtained using a slider crank mechanism with a flexible connecting rod.     

Multi-step transversal and tangential linearization methods applied to a class of nonlinear beam equations

A novel and continuously parameterized form of multi-step transversal linearization (MTrL) method is developed and numerically explored for solving nonlinear ordinary differential equations governing a class of boundary value problems (BVPs) of relevance in structural mechanics. A similar family of multi-step tangential linearization (MTnL) methods is also developed and applied to such BVP-s. Within the framework of MTrL and MTnL, a BVP is treated as a constrained dynamical system, i.e. a constrained initial value problem (IVP). While the MTrL requires the linearized solution manifold to transversally intersect the nonlinear solution manifold at a chosen set of points across the axis of the independent variable, the essential difference of the present MTrL method from its previous version [Roy, D., Kumar, R., 2005. A multistep transversal linearization (MTL) method in nonlinear structural dynamics. J. Sound Vib. 17, 829–852.] is that it has the flexibility of treating nonlinear damping and stiffness terms as time-variant damping and stiffness terms in the linearized system. The resulting time-variant linearized system is then solved using Magnus’ characterization [Magnus, W., 1954. On the exponential solution of differential equations for a linear operator. Commun. Pure Appl. Math., 7, 649–673.]. Towards numerical illustrations, response of a tip loaded cantilever beam (Elastica) is first obtained. Next, the response of a simply supported nonlinear Timoshenko beam is obtained using a variationally correct (VC) model for the beam [Marur, S., Prathap, G., 2005. Nonlinear beam vibration problems and simplification in finite element model. Comput. Mech. 35(5), 352–360.]. The new model does not involve any simplifications commonly employed in the finite element formulations in order to ease the computation of nonlinear stiffness terms from nonlinear strain energy terms. A comparison of results through MTrL and MTnL techniques consistently indicate a superior quality of approximations via the transversal linearization technique. While the usage of tangential system matrices is common in nonlinear finite element practices, it is demonstrated that the transversal version of linearization offers an easier and more general implementation, requires no computations of directional derivatives and leads to a consistently higher level of numerical accuracy. It is also observed that higher order versions of MTrL/MTnL with Lagrangian interpolations may not work satisfactorily and hence spline interpolations are suggested to overcome this problem.     

Generalization of the method of full approximation and its applications

This paper presents a generalized form of the method of full approximation.By usingthe concept of asymptotic linearization and making the coordinate transformationsincluding the nonlinear functionals of dependent variables,the original nonlinear problemsare linearized and their higher-order solutions are given in terms of the first-termasymptotic solutions and corresponding transformations.The analysis of a model equationand some problems of weakly nonlinear oscillations and waves with the generalized methodshows that it is effective and straightforward.     

Robust trajectory tracking control for a laboratory helicopter

This paper presents a robust nonlinear control strategy to deal with the trajectory tracking control problem for a laboratory helicopter. The helicopter model is considered as a nominal one with uncertainties such as unmodeled nonlinear dynamics, parametric uncertainties, and external disturbances. The proposed control approach incorporates the feedback linearization technique (FLT) and the signal compensation technique. The FLT is first applied to achieve the linearization of the nominal nonlinear model for reducing the conservation of the robust compensator design. A nominal controller based on the linear quadratic regulation method is designed for the linearized nominal system, whereas a robust compensator is introduced to restrain the influences of the uncertainties. It is shown that the trajectory tracking errors of the closed-loop system are ultimately bounded, and the boundaries can be specified by choosing the controller parameters. Simulation and experimental results on the lab helicopter verify the effectiveness of the proposed method.     

Development of the “binary interaction” theory for entangled polydisperse linear polymers

Based on the theoretical principles previously described in the literature, the development of the “naïve” binary interaction model is detailed in this paper. The new theory is effectively a sweeping generalization of the “Double Reptation” model. The “switch function” has been shown to be an essential feature of any constraint release model for Doi–Edwards type molecular models that invoke the concept of a discrete slip-link tube and is used in our formulation. Using the assumption of a constant entanglement density, a slip link linear density evolution equation is derived to rigorously count matrix entanglements. This function has no counterpart in the conventional Doi–Edwards theory, or its derivatives, and is absolutely required to properly generalize the “Double Reptation” model so that nonlinear flows can be modeled. The binary interaction polydispersity model is complex mathematically but can be rigorously and justifiably simplified by suppressing the tube coordinate dependence using a boundary layer analysis. The simplification process can be continued to the continuum level to create a hierarchy of approximate binary interaction models, thereby making large-scale numerical simulations of complex flows viable, indeed straightforward.     

Vibration absorbers for a rotating flexible structure with cyclic symmetry: nonlinear path design

This paper analytically investigates the nonlinear dynamics of order-tuned vibration absorbers applied to cyclic rotating flexible structures under traveling wave (TW) engine-order excitation. The primary cyclic structure is assumed to be governed by linear vibrations and the nonlinear absorber response arises from large amplitude kinematic effects. These dynamics are captured by a lumped-parameter model that consists of N blades with one blade mode and one absorber per blade, which are arranged with cyclic symmetry on a rotating disk. The governing equations of motion are formulated for arbitrary absorber paths to allow investigation of the absorber path design for nonlinear response. This paper extends previous work by the authors, which considered the linearized blade and absorber dynamics of a similar system. Several intriguing features of the dynamics were uncovered, most notably the existence of an absorber tuning range that avoids resonance at any rotation speed. Of particular interest is the existence and stability of the steady-state TW response to TW excitation, as experienced in turbomachinery, and how these are affected by selection of the absorber paths, which fix the linear and nonlinear tuning characteristics. It is shown that the TW response, which is unique for the linearized system, also exists for the weakly nonlinear model and can be captured by an equivalent two degree of freedom model obtained using the symmetry of the excitation and system response. The forced response exhibits the usual characteristics of a weakly nonlinear system, specifically, bistability and the attendant hysteresis near resonance. More significantly, it does not experience any additional instabilities associated with the symmetry. That is, the desired TW response is robust to nonlinear effects in the absorber, which allows use of the simple equivalent model for selection of absorber tuning parameters. For good performance and robustness, the linear absorber tuning should be in the “no-resonance zone” described by the linear theory and the absorber paths should have a slightly softening nonlinear characteristic.     

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].     

Error estimate of the modified point-mass trajectory model of an artillery shell

The article deals with the motion of an axially symmetric spinning artillery shell in the gravity field under the action of the system of aerodynamic forces and moments adopted in ballistics. As the starting point, the system of differential equations of motion of the shell is taken, which is obtained from the original “accurate” system by its linearization in the variables describing the angular motion of the symmetry axis and by additional linearization in the angle between the velocity vector of the center of mass and the vertical plane (l-system). This article examines the system of differential equations of the translational motion and axial rotation of the shell which describes its modified point-mass trajectory model as applied to l-system (m-system). By small parameter methods, an estimate is obtained for the difference of the solution of l-system with given initial data and the solution of m-system with the same initial data for the variables of translational motion and axial rotation. This analytical evaluation is built in such a way that it corresponds with certain numerical estimates for components of the translational motion and axial rotation. It is observed that, under accepted assumptions, m-system and l-system determine the translational motion of the shell with the same order of the error as compared to the original “accurate” nonlinear system of equations of motion of the shell. But m-system does not contain rapidly oscillating variables describing the angular motion of the symmetry axis, and so its numerical integration requires tens of times less computational resources than the numerical integration of l-system. Numerical simulation data are represented.     

Stochastic optimal active control of a half-car nonlinear suspension under random road excitation

The random response analysis and the stochastic optimal active control of a half-car model with nonlinear suspension stiffness and damping traversing a rough road are studied in this paper. The road roughness height is modeled as the output of a first-order linear filter to Gaussian white noise. Considering the hysteretic nonlinear stiffness and the square damping of the vehicle model, the response statistics of the nonlinear suspension with active control are obtained by using the equivalent linearization method. The performance indexes of the active suspension are evaluated and compared with those of the corresponding passive suspension. It is found that the nonlinear active suspension gives a better vehicle performances like ride comfort, suspension stroke and overall performance. Finally, the theoretical results are verified through Monte Carlo simulation.     

TWO-DIMENSIONAL VORTEX-INDUCED VIBRATION OF CABLE SUSPENSIONS

Resonant responses of suspended elastic cables driven by a steady current are investigated. Phenomenological fluid force models for alternate vortex-shedding are coupled with the nonlinear partial differential equations of cable motion. Decoupled cross-flow and in-line vortex-induced vibrations (VIV) are examined first using linearized and nonlinear cable models. The linearized cable model predicts well the basic characteristics of VIV and the nonlinear cable model captures the hysteresis often observed in experiments. Next, coupled cross-flow and in-line vibrations are evaluated by considering two principal coupling mechanisms: (i) cable structural nonlinearities, and (ii) coupled fluid lift and drag. Attention is focused on a “worst-case” resonant response where the natural frequencies for cable modes in the cross-flow and in-line directions are in the same 1:2 ratio as the excitation frequencies associated with lift and drag. The inclusion of cable structural nonlinearities alone leads to coupled responses that differ qualitatively (i.e., in number and stability of periodic motions) when compared to those of the decoupled model. The inclusion of coupled fluid lift and drag produces non-planar “figure eight” motions of the cable cross-section that exhibit similar characteristics to those previously measured on spring supported cylinders.     

纤维增强陶瓷基复合材料的D失效判据

经典唯象强度理论适用于正交各向异性线弹性体。对于非线性纤维增强复合材料,通过加卸载试验和损伤力学的分析方法,可以得到一种虚拟的线性化应力-应变关系;依据损伤等效假设,针对线性损伤和非线性损伤,对基于应力的经典二次失效准则进行变换,建立了一种基于损伤的强度理论,即“D失效判据”,这一强度理论可以作为经典判据的补充和扩展。针对平纹编织C/SiC复合材料的拉/剪组合试验,进行了实例计算,结果表明：利用D失效判据预测的失效包络线比蔡-希尔准则的预测曲线低,而且,失效曲线的形式与材料的损伤演化规律相关。     

中心刚体-柔性梁系统的最优跟踪控制

对考虑阻尼影响的中心刚体-柔性梁系统的动力特性和主动控制进行研究. 研究 中考虑了3种动力学模型：一次近似耦合模型、一次近似简化模型和线性化模型. 一次近 似模型中同时考虑了柔性梁的轴向变形和横向变形. 若在一次近似耦合模型中忽略轴向变 形的影响,则可得出一次近似简化模型. 线性化模型是对一次近似简化模型的线性化处理. 另外研究中考虑了3种阻尼因素：结构阻尼、风阻、中心刚体轴承处的阻尼. 控制设计采 用最优跟踪控制方法. 给出了从物理测量中提取模态坐标的滤波器方法. 研究结果显 示,一次近似简化模型能够有效地对系统的动力学行为进行描述;阻尼对系统的动力学特 性有着重要影响;当系统大范围运动为低速时,模态滤波器能够较好地提取出控制律所需 的模态坐标,最优跟踪控制方法能够使得系统跟踪所期望的运动轨迹,并且柔性梁的弹性 振动可得到抑制.     

考虑间隙反馈控制时滞的磁浮车辆稳定性研究

常导磁吸型(EMS)磁悬浮列车在悬浮控制中的每个环节,时滞是不可避免的,当时滞超过一定程度后,系统有可能失稳.本文针对EMS磁浮列车控制环节的临界时滞与车辆参数(如运行速度、反馈控制增益、导轨参数和悬挂参数)的关系开展研究.建立了磁浮车辆/导轨耦合动力学模型,车辆包含1节车辆和4个磁浮架,考虑车辆的10个自由度,每个磁浮架上包含4个悬浮电磁铁.导轨模拟为一系列简支Bernoulli-Euler梁,采用模态叠加法对导轨振动方程进行求解.采用传统线性电磁力模型实现车辆和轨道的耦合.采用比例-微分控制算法对电磁铁电流进行反馈控制,实现车辆稳定悬浮,并假设时滞均发生在控制环节,且只考虑间隙反馈控制环节的时滞.采用四阶龙格库塔法对耦合系统动力学方程进行求解,编写了数值仿真程序,计算得到车辆导轨耦合系统在考虑间隙反馈控制时滞时的响应.将系统运动发散时的时滞大小视为临界时滞,开展了参数规律影响分析.通过分析,给出了提高时滞条件下车辆稳定性的方法,包括增大导轨的弯曲刚度和阻尼比,减小间隙反馈控制增益并增大速度反馈控制增益,以及增大二系悬挂阻尼.      

Breakdown of elasticity theory for jammed hard-particle packings: conical nonlinear constitutive theory

Hard-particle packings have provided a rich source of outstanding theoretical problems and served as useful starting points to model the structure of granular media, liquids, living cells, glasses, and random media. The nature of “jammed” hard-particle packings is a current subject of keen interest. We demonstrate that the response of jammed hard-particle packings to global deformations cannot be described by linear elasticity (even for small particle displacements) but involves a “conical” nonlinear constitutive theory. It is the singular nature of the hard-particle potential that leads to the breakdown of linear elasticity. Interestingly, a nonlinear theory arises because the feasible particle displacements (leading to unjamming) depend critically on the local spatial arrangement of the particles, implying a directionality in the feasible strains that is absent in particle systems with soft potentials. Mathematically, the set of feasible strains has a conical structure, i.e., components of the imposed strain tensor generally obey linear inequalities. The nature of the nonlinear behavior is illustrated by analyzing several specific packings. Finally, we examine the conditions under which a packing can be considered to “incompressible” in the traditional sense.     

Large displacement analysis of shear deformable composite beams with interlayer slips

This paper presents a novel geometric non-linear finite element formulation for the analysis of shear deformable two-layer beams with interlayer slips. We adopt the co-rotational approach where the motion of the element is decomposed into two parts: a rigid body motion which defines a local coordinate system and a small deformational motion of the element relative to this local coordinate system. The main advantage of this approach is that the transformation matrices relating local and global quantities are independent to the choice of the geometrical linear local element. The effect of transverse shear deformation of the layers is taken into account by assuming that each layer behaves as a Timoshenko beam element. The layers are assumed to be continuously connected and partial interaction is considered by considering a continuous relationship between the interface shear flow and the corresponding slip. In order to avoid curvature and shear locking phenomena, the local linear element is formulated using “exact” displacement shape functions derived from the closed-form solution of the governing equations of a two-layer beam element. Finally, three numerical applications are presented in order to assess the performance of the proposed formulation.     

Adaptive nonlinear model predictive control design of a flexible-link manipulator with uncertain parameters

This paper presents a novel adaptive nonlinear model predictive control design for trajectory tracking of flexible-link manipulators consisting of feedback lineariza-tion, linear model predictive control, and unscented Kalman filtering. Reducing the nonlinear system to a linear system by feedback linearization simplifies the optimization prob-lem of the model predictive controller significantly, which, however, is no longer linear in the presence of parame-ter uncertainties and can potentially lead to an undesired dynamical behaviour. An unscented Kalman filter is used to approximate the dynamics of the prediction model by an online parameter estimation, which leads to an adaptation of the optimization problem in each time step and thus to a better prediction and an improved input action. Finally, a detailed fuzzy-arithmetic analysis is performed in order to quantify the effect of the uncertainties on the control structure and to derive robustness assessments. The control structure is applied to a serial manipulator with two flexible links containing uncertain model parameters and acting in three-dimensional space.     

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

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