Dymamics of a rotating flexible manipulator with tip load

The dynamics of a flexible manipulator is investigated in this paper. From the point of view of dynamic blance, the motion equations of a rotating beam with tip load are established by using Hamilton's principle. By taking into account the effects of dynamic stiffening and dynamic softening, the stability of the system is proved by employing Lyapunov's approach. Furthermore, the method of power series is proposed to find the exact solution of the eigenvalue problem. The effects of rotating speed and tip load on the vibration behavior of the flexible manipulator are shown in numerical results. Supported by National Natural Science Foundation. of China     

Stability of curvilinear Euler-Bernoulli beams in temperature fields

In this work, stability of thin flexible Bernoulli-Euler beams is investigated taking into account the geometric non-linearity as well as a type and intensity of the temperature field. The applied temperature field T(x,z) is yielded by a solution to the 2D Laplace equation solved for five kinds of thermal boundary conditions, and there are no restrictions put on the temperature distribution along the beam thickness. Action of the temperature field on the beam dynamics is studied with the help of the Duhamel theory, whereas the motion of the beam subjected to the thermal load is yielded employing the variational principles.The heat transfer (Laplace equation) is solved with the use of the finite difference method (FDM) of the third-order accuracy, while the integrals along the beam thickness defining the thermal stress and moments are computed using Simpson's method. Partial differential equations governing the beam motion are reduced to the Cauchy problem by means of application of FDM of the second-order accuracy. The obtained ordinary differential equations are solved with the use of the fourth-order Runge-Kutta method.The problem of numerical results convergence versus a number of beam partitions is investigated. A static solution for a flexible Bernoulli-Euler beam is obtained using the dynamic approach based on employment of the relaxation/set-up method.Novel stability loss phenomena of a beam under the thermal field are reported for different beam geometric parameters, boundary conditions, and the temperature intensity. In particular, it has been shown that stability of the flexible beam during heating the beam surface essentially depends on the beam thickness.     

Control structure interaction in the nonlinear analysis of flexible mechanical systems

The effect of the control structure interaction on the feedforward control law as well as the dynamics of flexible mechanical systems is examined in this investigation. An inverse dynamics procedure is developed for the analysis of the dynamic motion of interconnected rigid and flexible bodies. This method is used to examine the effect of the elastic deformation on the driving forces in flexible mechanical systems. The driving forces are expressed in terms of the specified motion trajectories and the deformations of the elastic members. The system equations of motion are formulated using Lagrange's equation. A finite element discretization of the flexible bodies is used to define the deformation degrees of freedom. The algebraic constraint equations that describe the motion trajectories and joint constraints between adjacent bodies are adjoined to the system differential equations of motion using the vector of Lagrange multipliers. A unique displacement field is then identified by imposing an appropriate set of reference conditions. The effect of the nonlinear centrifugal and Coriolis forces that depend on the body displacements and velocities are taken into consideration. A direct numerical integration method coupled with a Newton-Raphson algorithm is used to solve the resulting nonlinear differential and algebraic equations of motion. The formulation obtained for the flexible mechanical system is compared with the rigid body dynamic formulation. The effect of the sampling time, number of vibration modes, the viscous damping, and the selection of the constrained modes are examined. The results presented in this numerical study demonstrate that the use of the driving forees obtained using the rigid body analysis can lead to a significant error when these forces are used as the feedforward control law for the flexible mechanical system. The analysis presented in this investigation differs significantly from previously published work in many ways. It includes the effect of the structural flexibility on the centrifugal and Coriolis forces, it accounts for all inertia nonlinearities resulting from the coupling between the rigid body and elastic displacements, it uses a precise definition of the equipollent systems of forces in flexible body dynamics, it demonstrates the use of general purpose multibody computer codes in the feedforward control of flexible mechanical systems, and it demonstrates numerically the effect of the selected set of constrained modes on the feedforward control law.     

Finite element modeling of contact conditions in multibody system dynamics

In this paper a new method is developed for the dynamic analysis of contact conditions in flexible multibody systems undergoing a rolling type of motion. The relative motion between the two contacting bodies is treated as a constraint condition describing their kinematic and geometric relations. Equations of motion of the system are presented in a matrix form making use of Kane's equations and finite element method. The method developed has been implemented in a general purpose program called DARS and applied to the simulation and analysis of a rotating wheel on a track. Both the bodies are assumed flexible and discretized using a three dimensional 8-noded isoparametric elements. The time variant constraint conditions are imposed on the nodal points located at the peripheral surfaces of the bodies under consideration. The simulation is carried out under two different boundary conditions describing the support of the track. The subsequent constraint forces associated with the generalized coordinates of the system are computed and plotted. The effects of friction are also discussed.     

A nonlinear finite-element approach for kineto-static analysis of multibody systems

Methods that treat rigid/flexible multibody systems undergoing large motion as well as deformations are often accompanied with inefficiencies and instabilities in the numerical solution due to the large number of state variables, differences in the magnitudes of the rigid and flexible body coordinates, and the time dependencies of the mass and stiffness matrices. The kineto-static methodology of this paper treats a multibody mechanical system to consist of two collections of bulky (rigid) bodies and relatively flexible ones. A mixed boundary condition nonlinear finite element problem is then formulated at each time step whose known quantities are the displacements of the nodes at the boundary of rigid and flexible bodies and its unknowns are the deformed shape of the entire structure and the loads (forces and moments) at the boundary. Partitioning techniques are used to solve the systems of equations for the unknowns, and the numerical solution of the rigid multibody system governing equations of motion is carried out. The methodology is very much suitable in modelling and predicting the impact responses of multibody system since both nonlinear and large gross motion as well as deformations are encountered. Therefore, it has been adopted for the studies of the dynamic responses of ground vehicle or aircraft occupants in different crash scenarios. The kineto-static methodology is used to determine the large motion of the rigid segments of the occupant such as the limbs and the small deformations of the flexible bodies such as the spinal column. One of the most dangerous modes of injury is the amount of compressive load that the spine experiences. Based on the developed method, a mathematical model of the occupant with a nonlinear finite element model of the lumbar spine is developed for a Hybrid II (Part 572) anthropomorphic test dummy. The lumbar spine model is then incorporated into a gross motion occupant model. The analytical results are correlated with the experimental results from the impact sled test of the dummy/seat/restraint system. With this extended occupant model containing the lumbar spine, the gross motion of occupant segments, including displacements, velocities and accelerations as well as spinal axial loads, bending moments, shear forces, internal forces, nodal forces, and deformation time histories are evaluated. This detailed information helps in assessing the level of spinal injury, determining mechanisms of spinal injury, and designing better occupant safety devices.     

Dynamic modeling of a revolute-prismatic flexible robot arm fabricated from advanced composite materials

A dynamic model for a two degree-of-freedom planar robot arm is derived in this study. The links of the arm, connected to prismatic and revolute joints, are considered to be flexible. They are assumed to be fabricated from either aluminum or laminated composite materials. The model is derived based on the Timoshenko beam theory in order to account for the rotary inertia and shear deformation. These effects are significant in modeling flexible links connected to prismatic joints. The deflections of the links are approximated by using a shear-deformable beam finite element. Hamilton's principle is implemented to derive the equations describing the combined rigid and flexible motions of the arm. The resulting equations are coupled and highly nonlinear. In view of the large number of equations involved and their geometric nonlinearity (topological and quadratic), the solution of the equations of motion is obtained numerically by using a stiff integrator.The digital simulation studies examine the interaction between the flexible and the rigid body motions of the robot arm, investigate the improvement in the accuracy of the model by considering the flexibility of all rather than some of the links of the arm, assess the significance of the rotary inertia and shear deformation, and illustrate the advantages of using advanced composites in the structural design of robotic manipulators.     

Systematic modeling of a chain of N-flexible link manipulators connected by revolute–prismatic joints using recursive Gibbs-Appell formulation

The main intent of this paper is to represent a systematic algorithm capable of deriving the equations of motion of N-flexible link manipulators with revolute–prismatic joints. The existence of the prismatic joints together with the revolute ones makes the derivation of governing equations difficult. Also, the variations of the flexible parts of the links, with respect to time cause the associated mode shapes of the links to vary instantaneously. So, to derive the kinematic and dynamic equations of motion for such a complex system, the recursive Gibbs-Appell formulation is applied. For a comprehensive and accurate modeling of the system, the coupling effects due to the simultaneous rotating and reciprocating motions of the flexible arms as well as the dynamic interactions between large movements and small deflections are also included. In this study, the links are modeled based on the Euler–Bernoulli beam theory and the assumed mode method. Also, the effects of gravity as well as the longitudinal, transversal and torsional vibrations have been considered in the formulations. Moreover, a recursive algorithm based on 3 ×  3 rotational matrices has been applied in order to derive the system’s dynamic equations of motion, symbolically and systematically. Finally, a numerical simulation has been performed by means of a developed computer program in order to demonstrate the ability of this algorithm in deriving and solving the equations of motion related to such systems.     

Non-linear dynamic response of a rotating radial Timoshenko beam with periodic pulse loading at the free-end

Consideration is given to the dynamic response of a Timoshenko beam under repeated pulse loading. Starting with the basic dynamical equations for a rotating radial cantilever Timoshenko beam clamped at the hub in a centrifugal force field, a system of equations are derived for coupled axial and lateral motions which includes the transverse shear and rotary inertia effects, as well. The hyperbolic wave equation governing the axial motion is coupled with the flexural wave equation governing the lateral motion of the beam through the velocity-dependent skew-symmetric Coriolis force terms. In the analytical formulation, Rayleigh-Ritz method with a set of sinusoidal displacement shape functions is used to determine stiffness, mass and gyroscopic matrices of the system. The tip of the rotating beam is subjected to a periodic pulse load due to local rubbing against the outer case introducing Coulomb friction in the system. Transient response of the beam with the tip deforming due to rub is discussed in terms of the frequency shift and non-linear dynamic response of the rotating beam. Numerical results are presented for this vibro-impact problem of hard rub with varying coefficients of friction and the contact-load time. The effects of beam tip rub forces transmitted through the system are considered to analyze the conditions for dynamic stability of a rotating blade with intermittent rub.     

大变形柔性多体系统非线性动力学数值仿真与实验研究

提出了一种作大范围运动柔性梁的非接触动态测试技术.在基于位移的柔性多体系统几何精确建模及非线性有限元分析技术的基础上,利用EAGLE-500运动分析系统及其相应的分析软件对作大范围运动钛合金柔性梁作了实验研究,并且利用之前提出的几何精确梁理论进行数值仿真.数值仿真结果与实验结果完全吻合,验证了作者所提的几何精确梁理论及...     

Nonlinear dynamics of an inclined beam subjected to a moving load

In this paper, the nonlinear dynamic response of an inclined pinned-pinned beam with a constant cross section, finite length subjected to a concentrated vertical force traveling with a constant velocity is investigated. The study is focused on the mode summation method and also on frequency analysis of the governing PDEs equations of motion. Furthermore, the steady-state response is studied by applying the multiple scales method. The nonlinear response of the beam is obtained by solving two coupled nonlinear PDEs governing equations of planar motion for both longitudinal and transverse oscillations of the beam. The dynamic magnification factor and normalized time histories of mid-pint of the beam are obtained for various load velocity ratios and the outcome results have been illustrated and compared to the results with those obtained from traditional linear solution. The appropriate parametric study considering the effects of the linear viscous damping, the velocity of the traveling load, beam inclination angle under zero or nonzero axial load are carried out to capture the influence of the effect of large deflections caused by stretching effects due to the beam’s immovable ends. It was seen that quadratic nonlinearity renders the softening effect on the dynamic response of the beam under the act of traveling load. Also in the case where the object leaves the inclined beam, its planar motion path is derived and the targeting accuracy is investigated and compared with those from the rigid solution assumption. Moreover, the stability analysis of steady-state response for the modes equations having quadratic nonlinearity was carried out and it was observed from the frequency response curves that for the considered parameters in the case of internal-external primary resonance, both saturation phenomenon and jump phenomenon can be predicted for the longitudinal excitation.     

Dynamic analysis of a flexible hub-beam system with tip mass

For a dynamic system of a rigid hub and a flexible cantilever beam, the traditional hybrid coordinate model assumes the small deformation in structural dynamics where axial and transverse displacements at any point in the beam are uncoupled. This traditional hybrid coordinate model is referred as the zeroth-order approximation coupling model in this paper, which may result in divergence to the dynamic problem of some rigid–flexible coupling systems with high rotational speed. In this paper, characteristics of a flexible hub-beam system with a tip mass is studied. Based on the Hamilton theory and the finite element discretization method, and in consideration of the second-order coupling quantity of the axial displacement caused by the transverse displacement of the beam, the rigid–flexible coupling dynamic model (referred as the first-order approximation coupling (FOAC) model in this paper) and the corresponding model in non-inertial system for the flexible hub-beam system with a tip mass are presented firstly, then the dynamic characteristics of the system are studied through numerical simulations under twos cases: the large motion of the system is known and is unknown. Simulation and comparison studies using both the traditional zeroth-order model and the proposed first-order model show that even small tip mass may affect dynamic characteristics of the system significantly, which may result in the largening of vibrating amplitude and the descending of vibrating frequency of the beam, and may affect end position of the hub-beam system as well. The effect of the tip mass becomes large along with the increasing of rotating speed of large motion of the system. When the large motion of the system is at low speed, the traditional ZOAC model may lead to a large error, whereas the proposed FOAC model is valid. When the large motion is at high speed, the ZOAC model may result in divergence to the dynamic problem of the flexible hub-beam system, while the proposed second model can still accurately describe the dynamic hub-beam system.     

Soft-impact dynamics of deformable bodies

Systems constituted by impacting beams and rods of non-negligible mass are often encountered in many applications of engineering practice. The impact between two rigid bodies is an intrinsically indeterminate problem due to the arbitrariness of the velocities after the instantaneous impact and implicates an infinite value of the contact force. The arbitrariness of after-impact velocities is solved by releasing the impenetrability condition as an internal constraint of the bodies and by allowing for elastic deformations at contact during an impact of finite duration. In this paper, the latter goal is achieved by interposing a concentrate spring between a beam and a rod at their contact point, simulating the deformability of impacting bodies at the interaction zones. A reliable and convenient method for determining impact forces is also presented. An example of engineering interest is carried out: a flexible beam that impacts on an axially deformable strut. The solution of motion under a harmonic excitation of the beam built-in base is found in terms of transverse and axial displacements of the beam and rod, respectively, by superimposition of a finite number of modal contributions. Numerical investigations are performed in order to examine the influence of the rigidity of the contact spring and of the ratio between the first natural frequencies of the beam and the rod, respectively, on the system response, namely impact velocity, maximum displacement, spring stretching and contact force. Impact velocity diagrams, nonlinear resonance curves and phase portraits are presented to determine regions of periodic motion with impacts and the appearance of chaotic solutions, and parameter ranges where the functionality of the non-structural element is at risk.     

Dynamic recursive decoupling method for closed-loop flexible mechanical systems

In this paper, a new method for the dynamic analysis of a closed-loop flexible kinematic mechanical system is presented. The kinematic and force models are developed using absolute reference, joint relative, and elastic coordinates as well as joint reaction forces. This recursive formulation leads to a system of loosely coupled equations of motion. In a closed-loop kinematic chain, cuts are made at selected auxiliary joints in order to form spanning tree structures. Compatibility conditions and reaction force relationships at the auxiliary joints are adjoined to the equations of open-loop mechanical systems in order to form closed-loop dynamic equations. Using the sparse matrix structure of these equations and the fact that the joint reaction forces associated with elastic degrees of freedom do not represent independent variables, a method for decoupling the joint and elastic accelerations is developed. Unlike existing recursive formulations, this method does not require inverse or factorization of large non-linear matrices. It leads to small systems of equations whose dimensions are independent of the number of elastic degrees of freedom. The application of dynamic decoupling method in dynamic analysis of closed-loop deformable multibody systems is also discussed in this paper. The use of the numerical algorithm developed in this investigation is illustrated by a closed-loop flexible four-bar mechanism.     

Computer Implementation of the Absolute Nodal Coordinate Formulation for Flexible Multibody Dynamics

Deformable components in multibody systems are subject to kinematic constraints that represent mechanical joints and specified motion trajectories. These constraints can, in general, be described using a set of nonlinear algebraic equations that depend on the system generalized coordinates and time. When the kinematic constraints are augmented to the differential equations of motion of the system, it is desirable to have a formulation that leads to a minimum number of non-zero coefficients for the unknown accelerations and constraint forces in order to be able to exploit efficient sparse matrix algorithms. This paper describes procedures for the computer implementation of the absolute nodal coordinate formulation' for flexible multibody applications. In the absolute nodal coordinate formulation, no infinitesimal or finite rotations are used as nodal coordinates. The configuration of the finite element is defined using global displacement coordinates and slopes. By using this mixed set of coordinates, beam and plate elements can be treated as isoparametric elements. As a consequence, the dynamic formulation of these widely used elements using the absolute nodal coordinate formulation leads to a constant mass matrix. It is the objective of this study to develop computational procedures that exploit this feature. In one of these procedures, an optimum sparse matrix structure is obtained for the deformable bodies using the QR decomposition. Using the fact that the element mass matrix is constant, a QR decomposition of a modified constant connectivity Jacobian matrix is obtained for the deformable body. A constant velocity transformation is used to obtain an identity generalized inertia matrix associated with the second derivatives of the generalized coordinates, thereby minimizing the number of non-zero entries of the coefficient matrix that appears in the augmented Lagrangian formulation of the equations of motion of the flexible multibody systems. An alternate computational procedure based on Cholesky decomposition is also presented in this paper. This alternate procedure, which has the same computational advantages as the one based on the QR decomposition, leads to a square velocity transformation matrix. The computational procedures proposed in this investigation can be used for the treatment of large deformation problems in flexible multibody systems. They have also the advantages of the algorithms based on the floating frame of reference formulations since they allow for easy addition of general nonlinear constraint and force functions.     

Finite deflection dynamics of elastic beams

Solutions are obtained for the problem of an infinite elastic beam subjected to essentially constant velocity boundary conditions at one point of the beam. The effects of finite deflections, normal force, rotatory inertia and shear deformation are included. The equations of the problem are converted into non-dimensional form and a perturbation approach is used to obtain a consistent approximation. Numerical solutions are obtained for the bending moment, shear force and the normal force for different velocities of impact. It is shown that the solution to the problem depends on a combined geometrical and material parameter which does not vary significantly for compact sections and a loading parameter which determines the amplitude of the response. Finally the linear Timoshenko beam theory is shown to predict the bending moment and shear force extremely well even when the deflections are large enough to cause appreciable stretching of the centroidal axis.     

Bending-torsional stability analysis of aerodynamically covered pipes with inclined terminal nozzle and concurrent internal and external flows

Stability analysis of a cantilevered pipe with an inclined terminal nozzle as well as simultaneous internal and external fluid flows is investigated in this study. The pipe is embedded in an aerodynamic cover with negligible mass and stiffness simply to streamline the external flow and avoid vortex induced vibrations. The structure of pipe is modeled as an Euler–Bernoulli beam and effects of internal fluid flow including flow-induced inertia, Coriolis and centrifugal forces and the follower force induced by the exhausting jet are taken into account. In addition, neglecting the compressibility effect and using the unsteady Wagner model, aerodynamic loading is determined as a distributed lateral load for any generic structural state. The integral form of coupled equations of motion are obtained using the Hamilton’s principle. Solution to the coupled flexural–torsional equations of motion is realized via the extended Galerkin method. After discretization of the equations of motion, an eigenvalue representation of the problem is obtained. Several parameter studies are then conducted to examine the effects of concurrent fluid flows and other related parameters on the stability margins of the system.     

多体系统动力学逆问题的一种处理方法

首先用控制力法建立多体系统动力学逆问题的运动微分方程，然后给出奇异值分解在求系统动力学响应中的应用及作用在系统中的控制力的解法，最后举了一个算例。     

On the static and dynamic stability of thin beam conveying fluid

In this paper, numerical investigation of the statical and dynamical stability of aligned and misaligned viscoelastic cantilevered beam is performed with a terminal nozzle in the presence of gravity in two cases: (1) effect of fluid velocity on the flutter boundary of beam conveying fluid and (2) effect of gravity on the buckling boundary of beam conveying fluid. The beam is assumed to have a large width-to-thickness ratio, so the out-of-plane bending rigidity is far higher than the in-plane bending and torsional rigidities. Gravity vector is considered in the vertical direction. Thus, deflection of the beam because of the gravity effect couples the in-plane bending and torsional equations. The beam is modeled by Euler–Bernoulli beam theory, with the flow-induced inertia, Coriolis and centrifugal forces along the beam considered as a distributed load along the beam. Furthermore, the end nozzle is regarded as a lumped mass and modeled as a follower axial force. The extended Hamilton’s principle and the Galerkin method are utilized to derive the bending–torsional equations of motion. The coupled equations of motion are solved as eigenvalue problems. Also, several cases are examined to study the impact of gravity, beam inclination angle, mass ratio, nozzle aspect ratio, bending-to-bending rigidity ratio and bending-to-torsional rigidity ratio on flutter and buckling margin of the system.

     

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.