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

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

基于SE(3)群局部标架的5/6 Dofs CB壳单元

基于李群局部标架(local frame of Lie group,LFLG)的多柔体系统动力学建模方法可自然消除刚体运动带来的几何非线性,使系统的广义弹性力、广义惯性力及其雅可比矩阵满足刚体运动的不变性.本文融合李群局部标架思想和基于连续体(continuum-based,CB)的壳理论,提出基于SE(3)群局部标架...     

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

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

考虑几何非线性和热效应的刚-柔耦合动力学

温度增高和温度梯度会引起梁的纵向、横向变形位移,在一定程度上对刚-柔耦合规律产生影响.该文考虑热应变,从平面梁的非线性的应变与位移关系式出发,建立了刚体运动、弹性变形和温度相互耦合的有限元离散的热传导方程和动力学方程.研究热流作用下的中心刚体-简支梁系统的刚-柔耦合动力学性质,揭示了几何非线性项和热应变对弹性变形和刚体运动影响.     

基于高斯原理的多体系统动力学建模

基于高斯最小拘束原理,以释放中的绳系卫星为背景,建立地球引力场内变长度大变形柔索联系的多体系统动力学模型. 利用基尔霍夫动力学比拟方法将柔索的变形转化为刚性截面沿中心线的转动,使包含刚性分体与变形体的刚柔耦合系统转化为由柔索的刚性截面与刚性分体组成的广义多刚体系统. 由于刚性截面的局部小变形沿弧坐标的积累不受限制,适合描述柔索的超大变形. 文中对此刚柔耦合多体系统导出其在地球引力场中的拘束函数,考虑各分体在空间中相对位置的几何约束条件,利用拉格朗日乘子构成以条件极值问题为特征的数学模型. 将高斯原理用于多体系统动力学的建模,其特点是以寻求函数极值的变分方法代替微分方程,通过数值计算直接得出运动规律. 其形式统一,不随系统拓扑结构的变化而改变,也无需区分树系统或非树系统.对于带控制的多体系统,动力学分析还可根据技术需要与系统的优化结合进行.      

DYNAMIC MODELING OF MULTI-BODY SYSTEM BASED ON GAUSS'S PRINCIPLE

基于高斯最小拘束原理,以释放中的绳系卫星为背景,建立地球引力场内变长度大变形柔索联系的多体系统动力学模型. 利用基尔霍夫动力学比拟方法将柔索的变形转化为刚性截面沿中心线的转动,使包含刚性分体与变形体的刚柔耦合系统转化为由柔索的刚性截面与刚性分体组成的广义多刚体系统. 由于刚性截面的局部小变形沿弧坐标的积累不受限制,适合描述柔索的超大变形. 文中对此刚柔耦合多体系统导出其在地球引力场中的拘束函数,考虑各分体在空间中相对位置的几何约束条件,利用拉格朗日乘子构成以条件极值问题为特征的数学模型. 将高斯原理用于多体系统动力学的建模,其特点是以寻求函数极值的变分方法代替微分方程,通过数值计算直接得出运动规律. 其形式统一,不随系统拓扑结构的变化而改变,也无需区分树系统或非树系统.对于带控制的多体系统,动力学分析还可根据技术需要与系统的优化结合进行.     

作大运动弹性薄板中的几何非线性 与耦合变形1)

导出作大范围刚体运动弹性薄板包括了几何非线性和中面变形之间的相互耦合(耦合变形)的动力学控制方程.分析了几何非线性和耦合变形各自对系统动力学性质的影响,得到了在传统方法上只考虑几何非线性,系统将通过同宿轨分岔过渡到混沌运动;若在传统方法上考虑耦合变形,系统稳定且数值解收敛,与实际情形相符.     

作大运动弹性薄板中的几何非线性与耦合变形

导出作大范围刚体运动弹性薄板包括了几何非线性和中面变形之间的相互耦合（耦合变形）的动力学控制方程．分析了几何非线性和耦合变形各自对系统动力学性质的影响,得到了在传统方法上只考虑几何非线性,系统将通过同宿轨分岔过渡到混沌运动;若在传统方法上考虑耦合变形,系统稳定且数值解收敛,与实际情形相符．     

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

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

含芯拧绞绳非线性弯曲动力学特性分析与研究

建立细长缆索大柔性多体动力学模型时,现实存在的复杂捻制几何构型多不予考虑,而是将柔索简化为材料均匀梁进行描述,致使运动仿真模型与物理实际存在一定差距. 为此,研究一种典型非线性拧绞绳股的大变形等效动力学建模方法,考虑准静态与大范围运动情况下绳股内的线接触,计算了受摩擦力及弯曲曲率影响的绳股可变弯曲刚度,通过等效梁模型避免了绳股精细建模时的大规模计算消耗. 基于连续介质力学与绝对节点坐标方法,建立了拧绞绳惯性广义坐标下的多柔体动力学模型. 为了验证等效模型的可行性,与基于有限段方法建立的精细模型进行对比仿真分析,通过位形验证了等效模型的精度. 进一步地,根据力载作用下的准静态构型,研究了特定构型绳股弯曲刚度沿轴向的分布规律;通过自重力下一端固定柔性绳摆自由运动仿真并与传统均匀梁模型相比,研究了模型弯曲特性的差异. 最后,根据能量守恒原理分析了摩擦耗散系统内各种能量间的相互转化. 拧绞绳大变形等效动力学模型能够提高绳索动力系统运动预测的仿真计算效率,还能为钢丝绳参数与构型设计提供依据.      

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

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

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

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

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.     

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 Variational Approach to Dynamics of Flexible Multibody Systems

Abstract

This paper presents a variational formulation of constrained dynamics of flexible multibody systems, using a vector-variational calculus approach. Body reference frames are used to define global position and orientation of individual bodies in the system, located and oriented by position of its origin and Euler parameters, respectively. Small strain linear elastic deformation of individual components, relative to their body reference frames, is defined by linear combinations of deformation modes that are induced by constraint reaction forces and normal modes of vibration. A library of kinematic couplings between flexible and/or rigid bodies is defined and analyzed. Variational equations of motion for multibody systems are obtained and reduced to mixed differential-algebraic equations of motion. A space structure that must deform during deployment is analyzed, to illustrate use of the methods developed