Simulation of planar flexible multibody systems with clearance and lubricated revolute joints

Modeling of clearance joints plays an important role in the analysis and design of multibody mechanical systems. Based on the absolute nodal coordinate formulation (ANCF), a new computational methodology for modeling and analysis of planar flexible multibody systems with clearance and lubricated revolute joints is presented. A planar absolute nodal coordinate formulation based on the locking-free shear deformable beam element is implemented to discretize the flexible bodies. A continuous contact-impact model is used to evaluate the contact force, in which energy dissipation in the form of hysteresis damping is considered. A force transition model from hydrodynamic lubrication forces to dry contact forces is introduced to ensure continuity in the joint reaction force. A comprehensive study with different lubrication force models has also been carried out. The generalized-α method is used to solve the equations of motion and several efficient methods are incorporated in the proposed model. Finally, the methodology is validated by two numerical examples.     

A new model for dry and lubricated cylindrical joints with?clearance in spatial flexible multibody systems

A new approach to model and analyze flexible spatial multibody systems with clearance of cylindrical joints is presented in this work. The flexible parts are modeled by using absolute nodal coordinate formulation (ANCF)-based elements, while the rigid parts are described by employing the natural coordinate formulation (NCF), which can lead to a constant system mass matrix for the derived system equations of motion. In a simple way, a cylindrical joint with clearance is composed of two main elements, that is, a journal inside a bearing. Additionally, a lubricant fluid can exist between these two mechanical elements to reduce the friction and wear and increase the system??s life. For the case in which the joint is modeled as a dry contact pair, a technique using a continuous approach for the evaluation of the contact force is applied, where the energy dissipation in the form of hysteresis damping is considered. Furthermore, the frictional forces developed in those contacts are evaluated by using a modified Coulomb??s friction law. For the lubricated case, the hydrodynamic theory for dynamically loaded journal bearings is used to compute the forces generated by lubrication actions. The lubricated model is based on the Reynolds equation developed for the case of journal bearings with length-to-diameter ratios up to 1. Using this approach, the misalignment of the journal inside the bearing can be studied. Finally, two demonstrative examples of application are used to provide results that support the discussion and show the validity of the proposed methodologies.     

计及热应变的空间曲梁的刚-柔耦合动力学

研究带中心刚体的作大范围运动的空间曲梁的刚-柔耦合动力学.结合混合坐标法和绝对坐标法的特点,取与中心刚体大范围运动有关的变量和柔性梁各单元节点相对中心刚体连体基的位移和斜率作为广义坐标,建立了一种新的柔性梁的刚柔耦合模型.基于精确的应变和位移的关系式,根据Jourdian速度变分原理,建立了带中心刚体柔性曲梁的有限元离散的动力学方程.数值对比了空间曲梁系统和空间直梁系统的刚柔耦合动力学性质,用能量守恒规律验证了文中曲梁模型的合理性.在此基础上,在应变能中计及热应变,研究温度增高引起的曲梁的热膨胀对系统的动力学性态的影响.     

Efficient Evaluation of the Elastic Forces and the Jacobian in the Absolute Nodal Coordinate Formulation

This paper develops a new procedure for evaluating the elastic forces, the elastic energy and the jacobian of the elastic forces in the absolute nodal coordinate formulation. For this procedure, it is fundamental to use some invariant sparse matrices that are integrated in advance and have the property of transforming the evaluation of the elastic forces in a matrix multiplication process. The use of the invariant matrices avoids the integration over the volume of the element for every evaluation of the elastic forces. Great advantages can be achieved from these invariant matrices when evaluating the elastic energy and calculating the jacobian of the elastic forces as well. The exact expression of the jacobian of the differential system of equations of motion is obtained, and some advantages of using the absolute nodal coordinate formulation are pointed out. Numerical results show that there is important time saving as a result of the use of the invariant matrices.     

求解大柔度梁/绳索的ANCF单元及隐式迭代格式

为模拟大柔度梁/绳索结构的变形和大范围运动,基于绝对节点坐标方法ANCF(Absolute nodal coordi-nate formulation)和HHT(Hilber-Hughes-Taylor)积分方法,建立了ANCF单元的隐式动力学迭代格式.得到了简洁的节点等效力向量,且进一步导出了切线刚度矩阵的全部公式,...     

An exact nonlinear hybrid-coordinate formulation for flexible multibody systems

The previous low-order approximate nonlinear formulations succeeded in capturing the stiffening terms, but failed in simulation of mechanical systems with large deformation due to the neglect of the high-order deformation terms. In this paper, a new hybrid-coordinate formulation is proposed, which is suitable for flexible multibody systems with large deformation. On the basis of exact strain–displacement relation, equations of motion for flexible multibody system are derived by using virtual work principle. A matrix separation method is put forward to improve the efficiency of the calculation. Agreement of the present results with those obtained by absolute nodal coordinate formulation (ANCF) verifies the correctness of the proposed formulation. Furthermore, the present results are compared with those obtained by use of the linear model and the low-order approximate nonlinear model to show the suitability of the proposed models. The project supported by the National Natural Science Foundation of China (10472066, 50475021).     

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.     

Investigation on the choice of boundary conditions and shape functions for flexible multi-body system

The objective of this investigation is to examine the correctness and efficiency of the choice of boundary conditions when using assumed mode approach to simulate flexible multi-body systems. The displacement field due to deformation is approximated by the Rayleigh-Ritz assumed modes in floating frame of reference (FFR) formulation. The deformations obtained by the absolute nodal coordinate (ANC) formulation which are transformed by two sets of reference coordinates are introduced as a criterion to verify the accuracy of the simulation results by using the FFR formulation. The relationship between the deformations obtained from different boundary conditions is revealed. Numerical simulation examples demonstrate that the assumed modes with cantilevered-free, simply-supported and freefree boundary conditions without inclusion of rigid body modes are suitable for simulation of flexible multi-body system with large over all motion, and the same physical deformation can be obtained using those mode functions, differ only by a coordinate transformation. It is also shown that when using mode shapes with statically indeterminate boundary conditions, significant error may occur. Furthermore, the slider crank mechanism with rigid crank is accurate enough for investigating boundary condition problem of flexible multi-body system, which cost significant less simulating time.     

Implicit and explicit integration in the solution of the absolute nodal coordinate differential/algebraic equations

This investigation is concerned with the use of an implicit integration method with adjustable numerical damping properties in the simulation of flexible multibody systems. The flexible bodies in the system are modeled using the finite element absolute nodal coordinate formulation (ANCF), which can be used in the simulation of large deformations and rotations of flexible bodies. This formulation, when used with the general continuum mechanics theory, leads to displacement modes, such as Poisson modes, that couple the cross section deformations, and bending and extension of structural elements such as beams. While these modes can be significant in the case of large deformations, and they have no significant effect on the CPU time for very flexible bodies; in the case of thin and stiff structures, the ANCF coupled deformation modes can be associated with very high frequencies that can be a source of numerical problems when explicit integration methods are used. The implicit integration method used in this investigation is the Hilber–Hughes–Taylor method applied in the context of Index 3 differential-algebraic equations (HHT-I3). The results obtained using this integration method are compared with the results obtained using an explicit Adams-predictor-corrector method, which has no adjustable numerical damping. Numerical examples that include bodies with different degrees of flexibility are solved in order to examine the performance of the HHT-I3 implicit integration method when the finite element absolute nodal coordinate formulation is used. The results obtained in this study show that for very flexible structures there is no significant difference in accuracy and CPU time between the solutions obtained using the implicit and explicit integrators. As the stiffness increases, the effect of some ANCF coupled deformation modes becomes more significant, leading to a stiff system of equations. The resulting high frequencies are filtered out when the HHT-I3 integrator is used due to its numerical damping properties. The results of this study also show that the CPU time associated with the HHT-I3 integrator does not change significantly when the stiffness of the bodies increases, while in the case of the explicit Adams method the CPU time increases exponentially. The fundamental differences between the solution procedures used with the implicit and explicit integrations are also discussed in this paper.     

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.     

On the Use of Implicit Integration Methods and the Absolute Nodal Coordinate Formulation in the Analysis of Elasto-Plastic Deformation Problems

In the general theory of continuum mechanics, the state of rotation and deformation of material points can be uniquely defined from the displacement field by using the nine independent components of the displacement gradients. For this reason, the use of the absolute rotation parameters as nodal coordinates, without relating them to the displacement gradients, leads to coordinate redundancy that leads to numerical and fundamental problems in many existing large rotation finite element formulations. Because of this fundamental problem, special measures that require modifications of the numerical integration methods were proposed in the literature in order to satisfy the principle of work and energy. As demonstrated in this paper, no such measures need to be taken when the finite element absolute nodal coordinate formulation is used since the principle of work and energy are automatically satisfied. This formulation does not suffer from the problem of coordinate redundancy and ensures the continuity of stresses and strains at the nodal points. In this study, the use of the implicit integration methods with the consistent Lagrangian elasto-plastic tangent moduli is examined when the absolute nodal coordinate formulation is used. The performance of different numerical integration methods in the dynamic analysis of large elasto-plastic deformation problems is investigated. It is shown that all these methods, in the case of convergence, yield a solution that satisfies the principle of work and energy without the need of taking any special measures. Semi-implicit integration methods, however, can lead to numerical difficulties in the case of very stiff problems due to the linearization made in these methods in order to avoid the iterative Newton--Raphson procedure. It is also demonstrated that the use of the consistent Lagrangian-plastic tangent moduli derived in this investigation using the absolute nodal coordinate formulation leads to better convergence of the iterative Newton--Raphson procedure used in the implicit integration methods.     

Nonlinear formulation for flexible multibody system with large deformation

In this paper, nonlinear modeling for flexible multibody system with large deformation is investigated. Absolute nodal coordinates are employed to describe the displacement, and variational motion equations of a flexible body are derived on the basis of the geometric nonlinear theory, in which both the shear strain and the transverse normal strain are taken into account. By separating the inner and the boundary nodal coordinates, the motion equations of a flexible multibody system are assembled. The advantage of such formulation is that the constraint equations and the forward recursive equations become linear because the absolute nodal coordinates are used. A spatial double pendulum connected to the ground with a spherical joint is simulated to investigate the dynamic performance of flexible beams with large deformation. Finally, the resultant constant total energy validates the present formulation. The project supported by the National Natural Science Foundation of China (10472066, 10372057). The English text was polished by Yunming Chen.     

Coupled vibrations of a partially fluid-filled cylindrical container with an internal body including the effect of free surface waves

Internal bodies (baffles) are used as damping devices to suppress the fluid sloshing motion in fluid-structure interaction systems. An analytical method is developed in the present article to investigate the effects of a rigid internal body on bulging and sloshing frequencies and modes of a cylindrical container partially filled with a fluid. The internal body is a thin-walled and open-ended cylindrical shell that is coaxially and partially submerged inside the container. The interaction between the fluid and the structure is taken into account to calculate the sloshing and bulging frequencies and modes of the coupled system using the Rayleigh quotient, Ritz expansion and Galerkin method. It is shown that the present formulation is an appropriate and new approach to tackle the problem with good accuracy. The effects of fluid level, number of nodal diameters, internal body radius and submergence ratio on the dynamic characteristics of the coupled system are also investigated.     

On a Fluid-Structure Model in Which the Dynamics of the Structure Involves the Shear Stress Due to the Fluid

In this paper we consider a model for fluid-structure interaction. The hybrid system describes the interaction between an incompressible fluid in a three-dimensional container with interior a fixed domain and a thin elastic plate, the interface, which coincides with a flexible flat part of the surface of the vessel containing the fluid. The motion of the fluid is described by the linearized Navier–Stokes equations and the deformation of the plate by the classical plate equations for in-plane motions, modified to include the viscous shear stress which the fluid exerts on the plate as well as damping of Kelvin–Voigt type. We establish the existence of a unique weak solution of the interactive system of partial differential equations by considering an appropriate variational formulation. Uniform stability of the energy associated with the model is shown under the assumption that the potential plate energy is dominated by the dissipation induced by the viscosity of the fluid. The retention of the physical parameters in the problem is an a priori requirement in this physical condition.      

Development of flexible telescopic boom model using absolute nodal coordinate formulation sliding joint constraints with LuGre friction

In this investigation, a modeling procedure of a telescopic boom of cranes is developed using the absolute nodal coordinate formulation together with the sliding joint constraints. Since telescopic booms are extracted and retracted under various operating conditions, the overall length of the boom changes dynamically, leading to the time-variant vibration characteristics. For modeling the telescopic structure of booms, a special care needs to be exercised since the location of the sliding contact point moves along the deformable axis of the flexible boom and the solution to a moving boundary problem is required. This issue indeed makes the modeling of the telescopic boom difficult, despite the significant needs for the analysis. It is, therefore, the objective of this investigation to develop a modeling procedure for the flexible telescopic boom by considering the sliding contact condition with the dynamic frictional effect. To this end, the sliding joint constraint developed for the absolute nodal coordinate formulation is employed for describing relative sliding motion between flexible booms, while flexible booms are modeled using the beam element of the absolute nodal coordinate formulation, which allows for modeling the large rotation and deformation of the structure.     

A New Flexible Multibody Beam Element Based on the Absolute Nodal Coordinate Formulation Using the Global Shape Function and the Analytical Mode Shape Function

Several techniques for the reduced dimensionality of finite elementformulations were considered as component mode reduction methods in themiddle sixties. These techniques are widely used in flexiblemultibody simulations for solving small deformation problems. Theabsolute nodal coordinate formulation for solving large rotation anddeformation problems has been established as a full finite elementmethod instead of using similar kinds of reduction techniques. In thispaper, a reduced order absolute nodal coordinate formulation is newlyestablished by introducing the global beam shape function and theanalytical deformation modes as a full finite element. This formulationleads to a constant and symmetric mass matrix as the conventionalabsolute nodal coordinate formulation, and makes it possible to reducethe number of elements and system coordinates of the beam structurewhich undergoes large rotations and large deformations. Numericalexamples show that the excellent agreements between thepresent formulation and the conventional absolute nodal coordinateformulation using a large number of elements are examined. These results demonstratethat the present formulation has high accuracy in the sense that thepresent solutions are similar to the conventional ones with fewersystem coordinates, and high efficiency in computation.