柔性多体系统产生动力刚化原因的研究

传统的柔性多体系统建模理论由于对柔性体的变形及其与大范围运动产生惯性力之间的耦合处理得过于简单,所以在分析存在高速大范围运动柔性多体系统的动力学性态时会得到完全错误的结论。本文将通过对作大范围运动弹性薄板的讨论来揭示产生这种错误的及探讨对传统性多体系统建模理论作出改进的对策     

作大范围运动弹性梁刚—柔耦合动力学建模

利用弹性梁的变形理论和 Hamilton力学原理对作大范围运动弹性梁的刚 -柔耦合动力学建模理论进行了研究。分析了大范围运动对弹性梁的横向振动和纵向振动的影响 ,得到了大范围运动与弹性梁的中线耦合变形之间的耦合作用对该系统动力学性质有显著的影响 ,从而提出了作大范围运动弹性梁的刚柔耦合动力学模型     

Multibody Formulation with Large Bending Finite Elements (LBFE)

The goal of this paper is to present a flexible multibody formulation for Euler-Bernoulli beams involving large displacements. This method is based on a discretisation of internal and kinetic energies. The beam is represented by its line of centroids and each section is oriented by a frame defined by three Euler angles. We apply a finite element formulation to describe the evolution of these angles along the neutral fibre and describe the internal energy. The kinetic energy is approximated as the one of two rigid bars tangent to the neutral fibre at the ends of the beam. We derive the equations of motion from a Lagrange formulation. These equations are solved using the Newmark method or/and the Newton-Raphson technique. We solve some very classic problems taken from the literature as the curved beam presented by Simo [Simo, J. C., ‘A three-dimensional finite-strain rod model. the three-dimensional dynamic problem. Part I’, Comput. Meths. Appl. Mech. Engrg. 49, 1985, 55–70; Simo, J. C. and Vu-Quoc, L., ‘A three-dimensional finite-strain rod model, Part II: Computationals aspects’, Comput. Meths. Appl. Mech. Engrg. 58, 1988, 79–116] and Lee [Lee, Kisu, ‘Analysis of large displacements and large rotations of three-dimensional beams by using small strains and unit vectors’, Commun. Numer. Meth. Engrg. 13, 1997, 987–997] or the rotational rod presented by Avello [Avello, A., Garcia de Jalon, J., and Bayo, E., ‘Dynamics of flexible multibody systems using cartesian co-ordinates and large displacement theory’, Int. J. Num. Methods in Engineering 32, 1991, 1543–1563] and Simo [Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part I’ Jour. of Appl. Mech. 53, 1986, 849–854; Simo, J. C. and Vu-Quoc, L., ‘On the dynamics of flexible beams under large overall motions – the planar case. Part II’, Jour. of Appl. Mech. 53, 1986, 855–863].     

旋转中心刚体-FGM梁刚柔热耦合动力学特性研究

对旋转中心刚体-功能梯度材料(functionally graded material,FGM)梁刚柔热耦合动力学特性进行研究.FGM梁为物理性能参数沿厚度方向呈幂律分布的欧拉伯努利梁.考虑柔性梁的横向弯曲变形和轴向拉伸变形, 并计入横向弯曲变形引起的纵向缩短,即非线性耦合变形量.考虑变截面空心梁在外部高温、内冷通道冷却情况下的热力耦合对系统动力学特性的影响,求解得到FGM梁沿厚度方向分布的温度场, 进而在本构关系中计入热应变.采用假设模态法对柔性梁变形场进行离散,运用第二类拉格朗日方程推导得到系统的刚柔热耦合动力学方程,并编制动力学仿真软件, 然后通过仿真算例对系统的动力学问题进行研究.结果表明:不同截面梁动力学响应差异较大, 因此需对实际系统合理建模;大范围运动已知时, 考虑热冲击载荷的FGM梁将有效抑制横向弯曲变形,而大范围运动恒定时随热冲击的叠加会出现高频振荡; 大范围运动未知时,外力矩和热冲击载荷相互作用产生热力耦合效应, 导致系统呈现高频振荡,同时与中心刚体大范围旋转运动产生刚柔热耦合效应.      

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.     

Geometric nonlinear formulation for thermal-rigid-flexible coupling system

This paper develops geometric nonlinear hybrid formulation for flexible multibody system with large deformation considering thermal efect. Diferent from the conventional formulation, the heat flux is the function of the rotational angle and the elastic deformation, therefore, the coupling among the temperature, the large overall motion and the elastic deformation should be taken into account. Firstly,based on nonlinear strain–displacement relationship, variational dynamic equations and heat conduction equations for a flexible beam are derived by using virtual work approach,and then, Lagrange dynamics equations and heat conduction equations of the first kind of the flexible multibody system are obtained by leading into the vectors of Lagrange multiplier associated with kinematic and temperature constraint equations. This formulation is used to simulate the thermal included hub-beam system. Comparison of the response between the coupled system and the uncoupled system has revealed the thermal chattering phenomenon. Then, the key parameters for stability, including the moment of inertia of the central body, the incident angle, the damping ratio and the response time ratio, are analyzed. This formulation is also used to simulate a three-link system applied with heat flux. Comparison of the results obtained by the proposed formulation with those obtained by the approximate nonlinear model and the linear model shows the significance of considering all the nonlinear terms in the strain in case of large deformation. At last, applicability of the approximate nonlinear model and the linear model are clarified in detail.     

Comparison of Various Techniques for Modelling Flexible Beams in Multibody Dynamics

The modelling of flexible elements in mechanical systems has been widely investigated through several methods issuing from both the area of structural mechanics and the field of multibody dynamics. As regards the latter discipline, beside the problem of the generation of the multibody equations of motion, the choice of a spatial discretization method for modelling flexible elements has always been considered as a critical phase of the modelling. Although this subject is abundantly tackled in the open-literature, the latter probably lacks an objective comparison between the most commonly used approaches.This contribution presents an extensive investigation of several discretization techniques of flexible beams, in a pure multibody context. In particular, it is shown that shape functions based on power series monomials are very suitable and versatile to model beams being part of a multibody system and thus constitutes an interesting alternative to finite element analysis. For this purpose, a symbolic multibody program, in which various discretization techniques were implemented, was generalized to compute the equations of motion of a general multibody system containing flexible beams.     

基于李群局部标架的多柔体系统动力学建模与计算

多柔体系统动力学主要研究由多个具有运动学约束、存在大范围相对运动的柔性部件构成的动力学系统的建模、计算和控制.多柔体系统不仅具有柔体大变形导致的几何非线性,更具有大范围刚体运动引起的几何非线性,其非线性程度远高于计算结构力学所研究的几何非线性问题.本文基于李群局部标架(local frame of Lie group,...     

Modeling and analysis of a coupled rigid-flexible system

ＩｎｔｒｏｄｕｃｔｉｏｎＭａｎｙｍｅｃｈａｎｉｃａｌｓｙｓｔｅｍｓｓｕｃｈａｓｒｏｂｏｔｍａｎｉｐｕｌａｔｏｒ,ｔｕｒｂｉｎｅｂｌａｄｅｓ,ｈｅｌｉｃｏｐｔｅｒｐｒｏｐｅｌｌｅｒａｎｄｍａｎ＿ｍａｄｅｓａｔｅｌｌｉｔｅｗｉｔｈｆｌｅｘｉｂｌｅａｔｔａｃｈｍｅｎｔｓａｒｅｃｏｍｐｏｓｅｄｏｆｒｉｇｉｄａｎｄｆｌｅｘｉｂｌｅｂｏｄｉｅｓ．Ｔｈｅｙａｒｅｃａｌｌｅｄｃｏｕｐｌｅｄｒｉｇｉｄ＿ｆｌｅｘｉｂｌｅｓｙｓｔｅｍｓ．Ｃｏｒｒｅｃｔｐｒｅｄｉｃｔｉｏｎｓｏｆｔｈｅｂｅｈａｖｉｏｒｏｆｆｌｅｘｉｂｌ…     

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.     

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.     

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.     

考虑剪切效应的旋转FGM楔形梁刚柔耦合动力学建模与仿真

本文对一类中心刚体-柔性梁系统在大范围转动下的刚柔耦合动力学问题进行了研究. 柔性梁为功能梯度材料(functionally graded materials, FGM)楔形变截面梁,材料体积分数在梁轴向呈幂律分布变化. 以弧长坐标来描述柔性FGM梁的几何位移关系,分别使用倾角和拉伸应变变量描述柔性梁的横向弯曲和纵向拉伸变形,并计及剪切效应. 采用假设模态法离散变形场,运用第二类拉格朗日方程进行方程推导,得到系统考虑剪切效应的刚柔耦合动力学模型. 基于全新的刚柔耦合动力学建模理论,研究不同轴向材料梯度分布的FGM楔形梁,通过数值仿真计算,分析讨论不同的转速、梯度分布规律以及变截面参数对系统动力学特性的影响. 结果表明,剪切效应对大高跨比的FGM楔形梁的变形影响较为明显,不容忽略;材料梯度分布规律和截面参数的选取均会对旋转FGM楔形梁的动力学响应和频率产生较大影响. 本文提出的考虑剪切效应的倾角刚柔耦合动力学模型是对以往非剪切模型的进一步完善,可应用于工程中的 Timoshenko梁结构的动力学问题求解.      

Nonlinear dynamic analysis on rigid-flexible coupling system of an elastic beam

Previous work examined the effect of the attached stiffness matrix terms on stability of an elastic beam undergoing prescribed large overall motion. The aim of the present work is to extend the nonlinear formulations to an elastic beam with free large overall motion. Based on initial stress method, the nonlinear coupling equations of elastic beams are obtained with free large overall motion and the attached stiffness matrix is derived by solving sub-static formulation. The angular velocity and the tip deformation of the elastic pendulum are calculated. The analytical results show that the simulation results of the present model are tabled and coincide with the one-order approximate model. It is shown that the simulation results accord with energy conservation principle.     

刚柔耦合建模理论的实验验证

传统的混合坐标建模理论是零次近似方法，对刚柔耦合问题的描述存在缺陷，研究以一个由中心刚体，柔性梁及末端质量组成的刚柔耦合系统对象，建立了精确的一次近似的刚柔耦合动力学方程，在该模型中计及了结构阻尼及风阻的影响，利用单轴气浮台动力学实验平台，通过与实验数据的比较，说明了传统零次近似方法在某些条件下已不能下确描述柔性的刚性正确性和可靠性。     

刚-柔耦合动力学系统的建模理论研究

刚－柔耦合动力学系统的传统的混合坐标方法是零次近似方法,在建模过程中,直接套用的结构动力学的小变形假设,忽略了变形位移的高次耦合变形量．本文对柔性梁建立较零次近似更精确的高次耦合动力学模型,从连续介质力学理论出发,在变形位移中,计及横向位移引起的轴向缩短,导出变形位移的二次耦合量．用一致质量有限元方法对梁进行离散,基于Jourdain速度变分原理导出大范围运动为自由的柔性梁的刚－柔耦合动力学方程．计算了柔性重力摆的角速度和摆端点的横向变形,揭示零次近似模型和耦合模型的刚－柔耦合动力学性质的根本差异．     

链状柔性多体机器人系统动力学研究

本文基于Ｊｏｕｒｄａｉｎ变分原理建立了具有链状拓扑结构柔性多体机器人系统动力学通用模型，用在一致质量有限单元法及正则模态分析基础上引入的模态坐标描述构件的弹性形，用独立坐标描述相邻板件间的大位移运动，每个铰容许１－６个自由度，组强非线性惯性耦合的封闭形式的系统动力学微分方程组，文末对单弹性臂和双弹性臂机器人操作手进行动力学仿真。     

平动弹性梁的刚-柔耦合动力学

本文建立了作大范围平动弹性梁的刚-柔耦合动力学控制方程。分析了大范围平动对弹性梁变形运动动力学性质的影响，发现了大范围平动与变形运动之间的耦合动力学与大范围转动与变形运动之间的耦合动力学存在显著的差异。大范围平动使弹性梁的刚度降低，同时使系统阻尼增加；而大范围转动使弹性梁的刚度增加，同时使系统产生了能量转换的陀螺效应。因此，柔性多体系统刚-柔耦合动力建模中必须包括大范围平动与柔性体变形运动之间的耦合动力学效应。     

大范围运动细长柔性空间结构动力学特性分析

自由-自由边界无约束状态的细长柔性空间结构大范围运动时的动力学特性对整体结构运动分析和运动控制系统设计具有极其重要的作用。通过浮动坐标系建立结构的运动学关系;借助假设模态法对结构变形进行变量分离;利用Lagrange’s方程建立了结构的刚柔耦合振动方程;再通过Rayleigh-Ritz法,以无大范围运动时的振型函数作为基本解组,得到了大范围运动影响下的结构振动特征方程,求解该方程得到了结构频率和振型。通过几组数值算例的对比分析,指出了非耦合模型和耦合模型下结构频率及振型之间的差异。     

刚-柔耦合多体系统动力学建模与数值仿真

柔性多体系统动力学传统的混合坐标建模方法忽略了变形位移的高次耦合变形量，是一种零次近似方法，其适用范围受到限制。本文以中心刚体、柔性粱及末端质量组成的刚柔耦合系统为对象，考虑了有粘性阻尼及风阻的情况。在柔性粱的纵向变形位移中计及了横向位移引起的轴向变形，并采用有限元方法和Hamilton变分原理导出了系统的刚柔耦合一次近似的动力学方程。该方程充分计及了中心刚体的大范围运动与柔性粱的弹性变形运动的相互耦合，并采用一致的方法引入了阻尼因素。文中最后提供了一个比“动力刚化”问题更具有一般性的仿真计算反例，进一步说明了零次近似方法在处理某些刚一柔耦合动力学问题时的缺陷，同时表明了由一次近似模型可得到正确合理的结论。