大范围运动柔性机械臂斜碰撞动力学分析

对作大范围运动柔性机械臂系统,进行斜碰撞动力学分析.基于柔性多体系统刚柔耦合动力学理论,计入耦合变形项,全面考虑大范围刚体运动与弹性小变形运动的耦合,建立系统连续动力学方程.引入斜碰撞力学模型,将法向和切向碰撞力以广义力的形式加入动力学方程中,对系统进行斜碰撞动力学建模分析.法向碰撞模型选取基于连续接触力法的非线性弹簧阻尼模型,切向碰撞模型选取一种修正Coulomb摩擦模型,对切向摩擦力进行统一描述.给出接触、分离判据,实现不同状态的动力学模型转换与求解.对斜碰撞全局动力学进行了仿真验证,分析了柔性机械臂全局过程的动力学特性变化以及碰撞对大范围运动和小变形运动的作用,并对比了不同碰撞方向对大范围运动、变形、机械能、碰撞力等动力学参数的影响.     

多体系统中规则柔性板的离散建模方法研究

把柔性梁的离散坐标法有限段法扩展到规则柔性板中,视柔性板为带关节柔性(刚度、阻尼)的多刚体系统,详细阐述了离散坐标法的基本思想、理论依据,采用牛顿欧拉方法建立了动力学方程.借助通用有限元软件和动力学仿真程序验证了离散坐标法可以解决具有几何非线性变形的规则柔性板构件的多体系统动力学问题.     

考虑刚-柔-热耦合的板结构多体系统的动力学建模

对热载荷作用下中心刚体与大变形薄板多体系统的动力学建模问题进行研究.基于Kirchhoff假设,从格林应变和曲率与绝对位移的非线性关系式出发,推导了非线性广义弹性力阵,用绝对节点坐标法建立了大变形矩形薄板的有限元离散的动力学变分方程.为了考虑刚体姿态运动、弹性变形和温度变化的相互耦合作用,推导了热流密度与绝对节点坐标之间的关系式.引入系统的运动学约束方程,建立了中心刚体-矩形板多体系统的考虑刚-柔-热耦合的热传导方程和带拉格朗日乘子的第一类拉格朗日动力学方程.为了有效地提高计算效率,将改进的中心差分法和广义-α法相结合,求解热传导方程和动力学方程,差分后的方程通过牛顿迭代法耦合求解.对刚-柔耦合和刚-柔-热三者耦合两种模型的仿真结果进行比较表明,刚体运动对温度梯度和热变形的影响显著.此外,本文建模方法考虑了几何非线性项,因此也考虑了热膨胀引起的轴向变形对横向变形的影响.     

充气梁展开仿真研究

充气展开结构的动力学仿真主要涉及柔性结构非线性接触,基于有限变形的柔性多体动力学等关键问题,以充气梁为模型给出了适合仿真的材料模型、单元算法、接触算法和充气展开算法.利用软件LS-Dyna成功模拟了折叠充气梁的展开过程,推导出充气梁的弯曲刚度与压力的关系并给出计算结果,同时报告了充气展开仿真分析中,压力过高会导致充气梁失稳的现象,接着考察了压力对充气结构固有频率的影响,最后比较了不同压力下的冲击响应.仿真结果表明该方法研究充气结构动力学是可行的,有效的.     

基于等几何分析的柔性多体系统建模方法研究

近三十年来,柔性多体系统动力学取得长足进步,尤其是以绝对节点坐标方法(Absolute Nodal Coordinate Formulation, ANCF)为代表的非线性有限元已被用来处理复杂的柔性多体系统动力学问题.但绝对节点坐标方法采用斜率矢量作为广义坐标,导致系统自由度多,计算效率低.针对柔性多体系统,基于非均匀有理B样条(Non-Uniform Rational B-Splines, NURBS)曲线和曲面分别提出了Euler-Bernoulli细长梁单元和Kirchhoff-Love薄壳单元,在完全拉格朗日格式下,根据Green应变张量对单元变形进行描述,结合第二类Piola-Kirchhoff应力张量给出单元应变能公式,推导了单元的弹性力和弹性力雅可比矩阵表达式,最后通过静力学及动力学数值算例对提出的两类单元的性能进行对比和验证,为柔性多体系统建模提供了一种精确高效的新单元.     

空间充气展开绳网系统捕获目标自抗扰控制研究

空间充气展开绳网系统是依靠充气梁展开绳网进行目标捕获的航天器系统,具有更好的稳定性和可操控性.然而由于充气梁和绳网的大柔性变形以及捕获失稳自旋目标后的未知碰撞,使得捕获后的航天器姿态稳定控制困难.本文主要基于自抗扰控制解决了空间充气展开绳网系统捕获目标后的姿态稳定和消旋难题.首先,基于理想薄膜充压失效理论和绝对节点坐标方法建立充气展开绳网系统动力学模型,而后设计了航天器姿态稳定自抗扰控制器,用于实时估计并补偿系统捕获过程中未知惯量目标与捕获机构的碰撞干扰.仿真结果表明,动力学模型能够模拟捕获过程中充气梁的屈曲失效及碰撞特性,自抗扰控制器能够有效抑制碰撞带来的干扰,实现空间充气展开绳网系统捕获后的高精度姿态稳定控制,同时能够在有限时间内对自旋目标实现消旋.     

一种在轨可展开天线的多柔体动力学建模与计算

研究了一种考虑重力梯度的大型空间可展开天线桁架的动力学建模与计算问题.使用绝对节点坐标法和绝对节点坐标参考节点法建立了可展开结构的刚柔耦合动力学模型.利用离散方向导数构建了动力学方程的保能量动量时间积分算法,通过计算获得了含有多个模块的大型空间可展开天线在轨展开以及展开后轨道机动过程轨道-姿态-变形耦合动响应.通过数值计算结果与经典算例结果、商业软件ADAMS计算结果的对比分析验证了提出方法的正确性.     

基于Bezier插值方法的作大范围旋转运动锥形悬臂梁动力学研究

研究了作大范围旋转运动高度和宽度均沿着梁长度方向变化的锥形悬臂梁动力学问题. 采用 Bezier插值方法对柔性梁的变形场进行描述,考虑柔性梁的纵向拉伸变形和横向弯曲变形,计入由于横向弯曲变形引起的纵向缩短,即非线性耦合项. 运用第二类拉格朗日方程推导出作旋转运动锥形梁的动力学方程,并编制了动力学仿真软件,对作旋转运动锥形梁的频率和动力学响应进行研究. 结果表明: 不同锥形梁截面的动力学响应和系统频率将有明显差异,因此对实际系统合理建模,将能得到更为精确的结果.     

变拓扑柔性多体系统接触碰撞动力学研究

在实际工程领域中存在着大量接触碰撞等非连续动力学问题,现有的解决柔性多体系统连续动力学过程的建模理论与方法,已经无法解决或无法很好解决这些问题.本文基于变拓扑思想,提出了附加接触约束的柔性多体系统碰撞动力学建模理论；通过设计柔性圆柱杆接触碰撞实验,验证了所提出附加约束接触碰撞模型的有效性；针对柔性多体系统全局动力学仿真面临时间和空间的多尺度问题,提出多变量的离散方法,从而提高了柔性多体系统非连续动力学的仿真效率.     

多体系统动力学仿真中的接触碰撞模型分析

为了明确多体系统动力学仿真时接触碰撞模型的选择依据,对恢复系数模型,IMPACT函数模型,迟滞阻尼模型三种模型进行了分析讨论,并指出了模型中各参数的物理意义.最后以小球与平面碰撞问题为例,通过理论推导与试验相结合的方法确定了模型的参数,在ADAMS软件中进行了仿真计算,得到了三种模型的计算结果.计算结果表明,恢复系数方法只能得到碰撞后的速度;IMPACT函数模型可以得到较好的仿真结果,但是参数确定困难;修正的迟滞阻尼模型可以解决恢复系数接近1的接触问题.     

Experimental validation of rigid-flexible coupling dynamic formulation for hub–beam system

The aim of this paper is to compare the accuracy of the absolute nodal coordinate formulation and the floating frame of reference formulation for the rigid-flexible coupling dynamics of a three-dimensional Euler–Bernoulli beam by numerical and experimental validation. In the absolute nodal coordinate formulation, based on geometrically exact beam theory and considering the torsion effect, the material curvature of the beam is derived, and then variational equations of motion of a three-dimensional beam are obtained, which consist of three position coordinates, two slope coordinates, and one rotational coordinate. In the floating frame of reference formulation, the displacement of an arbitrary point on the beam is described by the rigid-body motion and a small superimposed deformation displacement. Based on linear elastic theory, the quadratic terms of the axial strain are neglected, and the curvatures are simplified to the first order. Considering both the linear damping and the quadratic air resistance damping, the equations of motion of the multibody system composed of air-bearing test bed and a cantilevered three-dimensional beam are derived based on the principle of virtual work. In order to verify the results of the computer simulation, two experiments are carried out: an experiment of hub–beam system with large deformation and a dynamic stiffening experiment. The comparison of the simulation and experiment results shows that in case of large deformation, the frequency result obtained by the floating frame of reference formulation is lower than that obtained by the experiment. On the contrary, the result obtained by the absolute nodal coordinate formulation agrees well with that obtained by the experiment. It is also shown that the floating frame of reference formulation based on linear elastic theory cannot reveal the dynamic stiffening effect. Finally, the applicability of the floating frame of reference formulation is clarified.     

A geometrically exact beam element based on the absolute nodal coordinate formulation

In this study, Reissner’s classical nonlinear rod formulation, as implemented by Simo and Vu-Quoc by means of the large rotation vector approach, is implemented into the framework of the absolute nodal coordinate formulation. The implementation is accomplished in the planar case accounting for coupled axial, bending, and shear deformation. By employing the virtual work of elastic forces similarly to Simo and Vu-Quoc in the absolute nodal coordinate formulation, the numerical results of the formulation are identical to those of the large rotation vector formulation. It is noteworthy, however, that the material definition in the absolute nodal coordinate formulation can differ from the material definition used in Reissner’s beam formulation. Based on an analytical eigenvalue analysis, it turns out that the high frequencies of cross section deformation modes in the absolute nodal coordinate formulation are only slightly higher than frequencies of common shear modes, which are present in the classical large rotation vector formulation of Simo and Vu-Quoc, as well. Thus, previous claims that the absolute nodal coordinate formulation is inefficient or would lead to ill-conditioned finite element matrices, as compared to classical approaches, could be refuted. In the introduced beam element, locking is prevented by means of reduced integration of certain parts of the elastic forces. Several classical large deformation static and dynamic examples as well as an eigenvalue analysis document the equivalence of classical nonlinear rod theories and the absolute nodal coordinate formulation for the case of appropriate material definitions. The results also agree highly with those computed in commercial finite element codes.     

Effect of the centrifugal forces on the finite element eigenvalue solution of a rotating blade: a comparative study

In this study, the effect of the centrifugal forces on the eigenvalue solution obtained using two different nonlinear finite element formulations is examined. Both formulations can correctly describe arbitrary rigid body displacements and can be used in the large deformation analysis. The first formulation is based on the geometrically exact beam theory, which assumes that the cross section does not deform in its own plane and remains plane after deformation. The second formulation, the absolute nodal coordinate formulation (ANCF), relaxes this assumption and introduces modes that couple the deformation of the cross section and the axial and bending deformations. In the absolute nodal coordinate formulation, four different models are developed; a beam model based on a general continuum mechanics approach, a beam model based on an elastic line approach, a beam model based on an elastic line approach combined with the Hellinger–Reissner principle, and a plate model based on a general continuum mechanics approach. The use of the general continuum mechanics approach leads to a model that includes the ANCF coupled deformation modes. Because of these modes, the continuum mechanics model differs from the models based on the elastic line approach. In both the geometrically exact beam and the absolute nodal coordinate formulations, the centrifugal forces are formulated in terms of the element nodal coordinates. The effect of the centrifugal forces on the flap and lag modes of the rotating beam is examined, and the results obtained using the two formulations are compared for different values of the beam angular velocity. The numerical comparative study presented in this investigation shows that when the effect of some ANCF coupled deformation modes is neglected, the eigenvalue solutions obtained using the geometrically exact beam and the absolute nodal coordinate formulations are in a good agreement. The results also show that as the effect of the centrifugal forces, which tend to increase the beam stiffness, increases, the effect of the ANCF coupled deformation modes on the computed eigenvalues becomes less significant. It is shown in this paper that when the effect of the Poisson ration is neglected, the eigenvalue solution obtained using the absolute nodal coordinate formulation based on a general continuum mechanics approach is in a good agreement with the solution obtained using the geometrically exact beam model.     

Rigid-flexible coupling dynamics of three-dimensional hub-beams system

In the previous research of the coupling dynamics of a hub-beam system, coupling between the rotational motion of hub and the torsion deformation of beam is not taken into account since the system undergoes planar motion. Due to the small longitudinal deformation, coupling between the rotational motion of hub and the longitudinal deformation of beam is also neglected. In this paper, rigid-flexible coupling dynamics is extended to a hub-beams system with three-dimensional large overall motion. Not only coupling between the large overall motion and the bending deformation, but also coupling between the large overall motion and the torsional deformation are taken into account. In case of temperature increase, the longitudinal deformation caused by the thermal expansion is significant, such that coupling between the large overall motion and the longitudinal deformation is also investigated. Combining the characteristics of the hybrid coordinate formulation and the absolute nodal coordinate formulation, the system generalized coordinates include the relative nodal displacement and the slope of each beam element with respect to the body-fixed frame of the hub, and the variables related to the spatial large overall motion of the hub and beams. Based on precise strain-displacement relation, the geometric stiffening effect is taken into account, and the rigid-flexible coupling dynamic equations are derived using velocity variational principle. Finite element method is employed for discretization. Simulation of a hub-beams system is used to show the coupling effect between the large overall motion and the torsional deformation as well as the longitudinal deformation. Furthermore, conservation of energy in case of free motion is shown to verify the formulation.     

Coupled thermo-structural analysis of a bimetallic strip using the absolute nodal coordinate formulation

A bimetallic strip consists of two different metal pieces that are bonded together. Due to the different coefficients of thermal expansion, exposing the strip to temperature induces thermal stresses that cause the structure to bend. Most often, incremental finite-element methods that introduce element nodal coordinates have been successfully applied to analyze the thermally induced vibrations in such systems. The exposure of these bimetallic strips to high temperatures results in large deflections and deformations, where the effects of the rigid-body motion and large rotations must be taken into account. For classic, non-isoparametric elements such as beams and plates the incremental methods do not result in zero strains under arbitrary, rigid-body motion. Therefore, in this paper a new model of a bimetallic strip is proposed based on a coupled thermo-structural analysis using the absolute nodal coordinate formulation. The applied, non-incremental, absolute nodal coordinate formulation uses a set of global displacements and slopes so that the beam and the plate elements can be treated as isoparametric elements. In order to simulate the bimetallic strip’s dynamic response, the formulation of the shear-deformable beam element had to be extended with thermally induced stresses. This made it possible to model the coupled thermo-structural problem and to represent the connectivity constraints at the interface between the two strips of metal. The proposed formulation was verified by comparing the responses using a general-purpose finite-element software.     

Analysis of Large Flexible Body Deformation in Multibody Systems Using Absolute Coordinates

To consider large deformation problems in multibody system simulations afinite element approach, called absolute nodal coordinate.formulation,has been proposed. In this formulation absolute nodal coordinates andtheir material derivatives are applied to represent both deformation andrigid body motion. The choice of nodal variables allows a fullynonlinear representation of rigid body motion and can provide the exactrigid body inertia in the case of large rotations. The methodology isespecially suited for but not limited to modeling of beams, cables andshells in multibody dynamics.This paper summarizes the absolute nodal coordinate formulation for a 3D Euler–Bernoulli beam model, in particular the definition of nodal variables, corresponding generalized elastic and inertia forces and equations of motion. The element stiffness matrix is a nonlinear function of the nodal variables even in the case of linearized strain/displacement relations. Nonlinear strain/displacement relations can be calculated from the global displacements using quadrature formulae.Computational examples are given which demonstrate the capabilities of the applied methodology. Consequences of the choice of shape.functions on the representation of internal forces are discussed. Linearized strain/displacement modeling is compared to the nonlinear approach and significant advantages of the latter, when using the absolute nodal coordinate formulation, are outlined.     

Large Deflection Analysis of a Thin Plate: Computer Simulations and Experiments

Many previous studies have conducted computer-aided simulations ofelastic bodies undergoing large deflections and deformations, but therehave not been many attempts to validate their numerical results. Thesubject of this paper is a thin clamped plate undergone large vibrationdue to attached end-point weight. The main aim of this paper is to showthe validity of the absolute nodal coordinate formulation (ANCF) bycomparing to the real experiments. Large oscillations of thin plates arestudied in the paper with taking into account effects of an attachedend-point weight and aerodynamic damping forces. The physicalexperiments are carried out using a high-speed camera and dataacquisition system. For numerical modeling of the plate, the absolutenodal coordinate formulation is used.     

Panel flutter analysis of plate element based on the absolute nodal coordinate formulation

In this study, aeroelastic analysis of a plate subjected to the external supersonic airflow is carried out. A 3-D rectangular plate element of variable thickness based on absolute nodal coordinate formulation (ANCF) has been developed for the structural model. In the approach to the problem, a continuum mechanics approach for the definition of the elastic forces within the finite element is considered. Both shear strain and transverse normal strain are taken into account. Linearized first-order potential (piston) theory is coupled with the structural model to account for pressure loading. Aeroelastic equations using ANCF are derived and solved numerically. Values of critical dynamic pressure are obtained by a modal approach, in which the mode shapes are obtained by ANCF. All the formulations and the computations are built up in a FORTRAN 90 computer program after it was confirmed by Mathematica^?, ver. 5. The results of free vibration analysis and flutter are compared with the available references and reasonable good agreement has been found. However, some results indicate that the known problem of locking (ANCF with uniform thickness) still persist in the current developed formulation.     

Poisson modes and general nonlinear constitutive models in the large displacement analysis of beams

Most existing formulations for structural elements such as beams, plates and shells do not allow for the use of general nonlinear constitutive models in a straightforward manner. Furthermore, such structural element models, due to the nature of the generalized coordinates used, do not capture some Poisson modes such as the ones that couple the deformation of the cross section of the structural element and stretch and bending. In this paper, beam models that employ general nonlinear constitutive equations are presented using finite elements based on the nonlinear absolute nodal coordinate formulation. This formulation relaxes the assumptions of the Euler–Bernoulli and Timoshenko beam theories, and allows for the use of general nonlinear constitutive models. The finite elements based on the absolute nodal coordinate formulation also allow for the rotation as well as the deformation of the cross section, thereby capturing Poisson modes which can not be captured using other beam models. In this investigation, three different nonlinear constitutive models based on the hyper-elasticity theory are considered. These three models are based on the Neo–Hookean constitutive law for compressible materials, the Neo–Hookean constitutive law for incompressible materials, and the Mooney–Rivlin constitutive law in which the material is assumed to be incompressible. These models, which allow capturing Poisson modes, are suitable for many materials and applications, including rubber-like materials and biological tissues which are governed by nonlinear elastic behavior. Numerical examples that demonstrate the implementation of these nonlinear constitutive models in the absolute nodal coordinate formulation are presented. The results obtained using the nonlinear and linear constitutive models are compared in this study. These results show that the use of nonlinear constitutive models can significantly enhance the performance and improve the computational efficiency of the finite element models based on the absolute nodal coordinate formulation. The results also show that when linear constitutive models are used in the large deformation analysis, singular configurations are encountered and basic formulas such as Nanson’s formula are no longer valid. These singular deformation configurations are not encountered when the nonlinear constitutive models are used.