考虑刚弹耦合作用的柔性多体连续系统动力学建模

基于Ｈａｍｉｌｔｏｎ原理建立起一般柔性体连续系统的动力学建模方法,进而以水平面内作大范围回转运动的柔性梁为例,在Ｅｕｌｅｒ－Ｂｅｒｎｏｕｌｌｉ梁模型的假设前提下,根据轴向不可伸长的柔性梁的几何约束条件;推导出作大范围刚体运动的柔性梁连续系统的一致线性化振动微分方程．采用假设模态法对其离散化,导出考虑刚弹耦合作用的柔性梁有限维离散化动力学模型．最后给出仿真算例,验证了该方法的有效性．     

耦合变形对大范围运动柔性梁动力学建模的影响

柔性梁在作大范围空间运动时,产生弯曲和扭转变形,这些变形的相互耦合形成了梁在纵向以及横向位移的二次耦合变量。本文考虑了变形产生的几何非线性效应对运动柔性梁的影响,在其三个方向的变形中均考虑了二次耦合变量,利用弹性旋转矩阵建立了准确的几何非线性变形方程,通过Lagrange方程导出系统的动力学方程。仿真结果表明,在大范围运动情况下,仅在纵向变形中计及了变形二次耦合量的一次动力学模型,与考虑了完全几何非线性变形的模型具有一定的差异。     

基于初应力法的作大范围运动柔性梁动力学建模理论研究

研究了初应力法的作大范围运动柔性梁的建模理论.根据连续介质理论,考虑应变-位移中的非线性项,用一致质量有限元法对柔性梁进行离散,基于Jourdain速度变分原理导出定轴转动下大范围运动为自由的柔性梁刚-柔耦合动力学方程.从其刚柔耦合动力学方程出发,考虑在大范围运动已知情况下的结构动力学方程.通过引入准静态概念,把其结构动力学方程转化为准静态方程.对纵向和横向变形节点坐标进行坐标分离,解出与纵向变形相关的准静态方程,得到准静态时的纵向应力表达式,从而获得附加刚度项.并对此非惯性系下作大范围运动柔性梁的结构动力学方程进行数值仿真,对零次近似模型、一次近似模型、初应力法动力学模型的仿真结果进行分析,揭示三种模型的动力学性质的差异.     

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

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

柔性梁刚柔耦合碰撞动力学及变形传播研究

运用柔性多体系统刚柔耦合动力学理论,研究了作大范围回转运动柔性梁的碰撞动力学问题.考虑柔性梁的横向变形,以及横向变形引起的纵向缩短项即非线性耦合变形项.采用基于Hertz接触理论及非线性阻尼理论的非线性弹簧阻尼模型来求解碰撞过程中产生的碰撞力,运用第二类拉格朗日方程建立了系统的刚柔耦合碰撞动力学方程.编制仿真软件进行动力学仿真计算,得到了碰撞力和系统动力学响应,对比分析了不同动力学模型对系统动力学响应的影响.同时研究了碰撞导致的柔性梁横向变形传播的波动特性.     

作大范围空间运动柔性梁的刚-柔耦合动力学

研究带中心刚体的作大范围空间运动梁的刚-柔耦合动力学问题．从精确的应变-位移关系式出发，在动力学变分方程中，考虑了横截面转动的惯性力偶和与扭转变形有关的弹性力的虚功率，用速度变分原理建立了考虑几何非线性的空间梁的刚-柔耦合动力学方程，用有限元法进行离散．通过对空间梁系统的数值仿真研究扭转变形和截面转动惯量对系统动力学性态的影响．     

作大范围回转运动柔性梁斜碰撞动力学研究

为正确估计由于碰撞引起的多柔体系统动力学特性的变化,针对作大范围回转运动的柔性梁与一固定斜面发生斜碰撞的情况,在考虑刚柔耦合效应的多柔体系统动力学建模理论的基础上,利用假设模态法建立起重力场作用下的柔性梁一致线性化动力法向碰撞过程中系统的动力行为。基于Ｈｅｒｔｚ接触理论和非线性阻尼项建立法向碰撞接触模型,基于线性切向接触刚度建立柔性梁切向碰撞接触模型,提出的数值算法保证了计算结果的合理性,给出的仿     

两类基于局部标架梁单元的闭锁缓解方法

对于大转动、大变形柔性体的刚柔耦合动力学问题,基于李群SE(3)局部标架(local frame formulation, LFF)的建模方法能够规避刚体运动带来的几何非线性问题,离散数值模型中广义质量矩阵与切线刚度矩阵满足刚体变换的不变性,可明显地提高柔性多体系统动力学问题的计算效率. 有限元方法中,闭锁问题是导致单元收敛性能低下的主要原因, 例如梁单元的剪切以及泊松闭锁.多变量变分原理是缓解梁、板/壳单元闭锁的有效手段. 该方法不仅离散位移场,同时离散应力场或应变场, 可提高应力与应变的计算精度. 本文基于上述局部标架,研究几类梁单元的闭锁处理方法, 包括几何精确梁(geometrically exact beam formulation, GEBF)与绝对节点坐标(absolute nodal coordinate formulation, ANCF)梁单元. 其中, 采用Hu-Washizu三场变分原理缓解几何精确梁单元中的剪切闭锁,采用应变分解法缓解基于局部标架的ANCF全参数梁单元中的泊松闭锁. 数值算例表明,局部标架的梁单元在描述高转速或大变形柔性多体系统时,可消除刚体运动带来的几何非线性, 极大地减少系统质量矩阵和刚度矩阵的更新次数.缓解闭锁后的几类局部标架梁单元收敛性均得到了明显提升.      

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

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

考虑曲率纵向变形效应的大变形柔性梁刚柔耦合动力学建模与仿真

对在平面内做大范围转动的中心刚体柔性梁系统的动力学进行了研究,建立了考虑大变形效应的系统刚柔耦合动力学模型,并进行了动力学仿真.该动力学模型不但考虑了柔性梁横向弯曲变形和纵向变形(包含轴向拉伸变形和横向弯曲变形而引起的纵向缩短项),还考虑了纵向变形对曲率的影响,称为曲率纵向变形效应.在以往的研究中,柔性梁的横向弯曲变形能往往直接使用柔性梁横向弯曲变形来表达,并没有考虑纵向变形的影响.为了考虑柔性梁纵向变形对横向弯曲变形能的影响,在浮动坐标系下使用柔性梁参数方程形式的精确曲率公式来计算柔性梁的弯曲变形能.在此基础上建立了基于浮动坐标系的考虑曲率纵向变形效应的刚耦合动力学模型.论文给出了数值仿真算例,验证了本文所建的动力学模型既能适用于柔性梁的小变形问题,又能适用于大变形问题,且较现有高次刚柔耦合动力学模型更加适用于大变形问题的处理.论文还通过与能处理柔性梁大变形问题的绝对节点坐标法的比较,验证了模型的正确性.      

Three kinematic representations for modeling of highly flexible beams and their applications

Presented here are three kinematic representations of large rotations for accurate modeling of highly flexible beam-like structures undergoing arbitrarily large three-dimensional elastic deformation and/or rigid-body motion. Different methods of modeling torsional deformation result in different beam theories with different mathematical characteristics. Each of these three geometrically exact beam theories fully accounts for geometric nonlinearities and initial curvatures by using Jaumann strains, exact coordinate transformations, and orthogonal virtual rotations. The derivations are presented in detail, a finite element formulation is included, fully nonlinear governing equations and boundary conditions are presented, and the corresponding form for numerically exact analysis using multiple shooting methods is also derived. These theories are compared in terms of their appropriate application areas, possible singular problems, and easiness for use in modeling and analysis of multibody systems. Nonlinear finite element analysis of a rotating beam and nonlinear multiple shooting analysis of a torsional bar are performed to demonstrate the capability and accuracy of these beam theories.     

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

多柔体系统动力学主要研究由多个具有运动学约束、存在大范围相对运动的柔性部件构成的动力学系统的建模、计算和控制.多柔体系统不仅具有柔体大变形导致的几何非线性,更具有大范围刚体运动引起的几何非线性,其非线性程度远高于计算结构力学所研究的几何非线性问题.本文基于李群局部标架(local frame of Lie group, LFLG),讨论如何发展一套新的多柔体系统动力学建模和计算方法体系, 具体内容包括:基于局部标架的梁、板壳单元,适用于长时间历程计算的多柔体系统碰撞动力学积分算法,结合区域分解技术的大规模多柔体系统动力学并行求解器, 以及若干验证性算例.上述基于李群局部标架的方法体系可在计算中消除刚体运动带来的几何非线性问题,使柔体系统的广义惯性力、广义弹性力及其雅可比矩阵满足刚体运动的不变性,使多柔体系统动力学与大变形结构力学相互统一,有望推动新一代多柔体系统动力学建模和计算软件的发展.      

Nonlinear dynamics of three-dimensional belt drives using the finite-element method

In this paper, new nonlinear dynamic formulations for belt drives based on the three-dimensional absolute nodal coordinate formulation are developed. Two large deformation three-dimensional finite elements are used to develop two different belt-drive models that have different numbers of degrees of freedom and different modes of deformation. Both three-dimensional finite elements are based on a nonlinear elasticity theory that accounts for geometric nonlinearities due to large deformation and rotations. The first element is a thin-plate element that is based on the Kirchhoff plate assumptions and captures both membrane and bending stiffness effects. The other three-dimensional element used in this investigation is a cable element obtained from a more general three-dimensional beam element by eliminating degrees of freedom which are not significant in some cable and belt applications. Both finite elements used in this investigation allow for systematic inclusion or exclusion of the bending stiffness, thereby enabling systematic examination of the effect of bending on the nonlinear dynamics of belt drives. The finite-element formulations developed in this paper are implemented in a general purpose three-dimensional flexible multibody algorithm that allows for developing more detailed models of mechanical systems that include belt drives subject to general loading conditions, nonlinear algebraic constraints, and arbitrary large displacements. The use of the formulations developed in this investigation is demonstrated using two-roller belt-drive system. The results obtained using the two finite-element formulations are compared and the convergence of the two finite-element solutions is examined.     

Modeling of joints with clearance in flexible multibody systems

This paper is concerned with the modeling of joints with clearance within the framework of finite element based dynamic analysis of nonlinear, flexible multibody systems. For actual joints, clearance, lubrication and friction phenomena can significantly affect the dynamic response of the system. In this work, the effects of clearance and lubrication are studied for revolute and spherical joints. The formulation is developed within the framework of energy preserving and decaying time integration schemes that provide unconditional stability for nonlinear, flexible multibody systems. Numerical examples are presented that demonstrate the efficiency and accuracy of the proposed approach. The importance of modeling structural damping and limited driving power are discussed.     

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.     

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.     

Geometrically exact curved beam element using internal force field defined in deformed configuration

This paper presents a study on the development of high-performance finite elements for geometrically nonlinear analysis of frame structures with curved members. Based on the geometrically exact beam theory, a highly efficient and accurate mixed finite element is developed. A new approach is proposed for constructing the independent internal force field by including major terms satisfying equilibrium conditions in the deformed configuration. An element-level equilibrium iteration procedure is employed for the condensation of element internal degrees of freedom during the nonlinear solution. Numerical results are presented to demonstrate the excellent performance of the element developed, and it is shown that even when each structural member is modelled with just one element, accurate solutions can still be achieved.     

Dynamic Bifurcation of Multibody Systems

In this paper, the process of loss of stability of multibody systems and structures is analyzed. A novel approach is presented and applied to the statically loaded spatial systems for the analysis of a dynamic response of systems imposed on impact, high velocity compulsive motion, or percussive forces. The analysis is based on the solution of the dynamic equations and eigenvalue problem of systems, and of the resultant motion simulation. The flexible systems are discretized using the finite element method. The dynamic equations are derived with respect to the relative coordinates of the finite elements. Large flexible deflections due to a loss of stability are simulated. The initial forms of the possible deformations are defined by the computed eigenvectors solving the eigenvalue problem for the system stiffness matrix. The critical forces and system deflections are then analyzed. Examples of bifurcation of beam and beam structure imposed on compulsive motion, percussive forces, and impact are presented.