An Internal Damping Model for the Absolute Nodal Coordinate Formulation

Introducing internal damping in multibody system simulations is important as real-life systems usually exhibit this type of energy dissipation mechanism. When using an inertial coordinate method such as the absolute nodal coordinate formulation, damping forces must be carefully formulated in order not to damp rigid body motion, as both this and deformation are described by the same set of absolute nodal coordinates. This paper presents an internal damping model based on linear viscoelasticity for the absolute nodal coordinate formulation. A practical procedure for estimating the parameters that govern the dissipation of energy is proposed. The absence of energy dissipation under rigid body motion is demonstrated both analytically and numerically. Geometric nonlinearity is accounted for as deformations and deformation rates are evaluated by using the Green–Lagrange strain–displacement relationship. In addition, the resulting damping forces are functions of some constant matrices that can be calculated in advance, thereby avoiding the integration over the element volume each time the damping force vector is evaluated.     

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.     

Modeling Method and Application of Rational Finite Element Based on Absolute Nodal Coordinate Formulation

In multibody system dynamics, the absolute nodal coordinate formulation(ANCF)uses power functions as interpolating polynomials to describe the displacement field. It can get accurate results for flexible bodies that undergo large deformation and large rotation. However, the power functions are irrational representation which cannot describe the complex shapes precisely, especially for circular and conic sections. Different from the ANCF representation,the rational absolute nodal coordinate formulation(RANCF) utilizes rational basis functions to describe geometric shapes, which allows the accurate representation of complicated displacement and deformation in dynamics modeling. In this paper, the relationships between the rational surface and volume and the RANCF finite element are provided, and the generalized transformation matrices are established correspondingly. Using these transformation matrices, a new four-node three-dimensional RANCF plate element and a new eight-node three-dimensional RANCF solid element are proposed based on the RANCF. Numerical examples are given to demonstrate the applicability of the proposed elements. It is shown that the proposed elements can depict the geometric characteristics and structure configurations precisely, and lead to better convergence in comparison with the ANCF finite elements for the dynamic analysis of flexible bodies.     

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.     

Use of Cholesky Coordinates and the Absolute Nodal Coordinate Formulation in the Computer Simulation of Flexible Multibody Systems

In a previous publication, procedures that can be used with the absolute nodal coordinate formulation to solve the dynamic problems of flexible multibody systems were proposed. One of these procedures is based on the Cholesky decomposition. By utilizing the fact that the absolute nodal coordinate formulation leads to a constant mass matrix, a Cholesky decomposition is used to obtain a constant velocity transformation matrix. This velocity transformation is used to express the absolute nodal coordinates in terms of the generalized Cholesky coordinates. The inertia matrix associated with the Cholesky coordinates is the identity matrix, and therefore, an optimum sparse matrix structure can be obtained for the augmented multibody equations of motion. The implementation of a computer procedure based on the absolute nodal coordinate formulation and Cholesky coordinates is discussed in this paper. Numerical examples are presented in order to demonstrate the use of Cholesky coordinates in the simulation of the large deformations in flexible multibody applications.     

Large Oscillations of a Thin Cantilever Beam: Physical Experiments and Simulation Using the Absolute Nodal Coordinate Formulation

Many papers have studied computer-aided simulations of elastic bodies undergoing large deflections and large deformations. But there have not been many attempts to check the validity of the numerical formulations used in these studies. The main aim of this paper is to demonstrate the validity of one such numerical formulation, the absolute nodal coordinate formulation (ANCF), by comparing the results it generates with the results of real experiments. Large oscillations of a thin cantilever beam are studied in this paper to numerically model the beam, which also accounts for the effects of an attached end-point weight and damping forces. The experiments were carried out using a high-speed camera and a data acquisition system.     

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.     

Formulation of Three-Dimensional Joint Constraints Using the Absolute Nodal Coordinates

A wide variety of mechanical and structural multibody systems consist ofvery flexible components subject to kinematic constraints. The widelyused floating frame of reference formulation that employs linear modelsto describe the local deformation leads to a highly nonlinear expressionfor the inertia forces and can be applied to only small deformationproblems. This paper is concerned with the formulation and computerimplementation of spatial joint constraints and forces using the largedeformation absolute nodal coordinate formulation. Unlike the floatingframe of reference formulation that employs a mixed set of absolutereference and local elastic coordinates, in the absolute nodalcoordinate formulation, global displacement and slope coordinates areused. The nonlinear kinematic constraint equations and generalized forceexpressions are expressed in terms of the absolute global displacementsand slopes. In particular, a new formulation for the sliding jointbetween two very flexible bodies is developed. A surface parameter isintroduced as an additional new variable in order to facilitate theformulation of this sliding joint. The constraint and force expressionsdeveloped in this paper are also expressed in terms of generalizedCholesky coordinates that lead to an identity inertia matrix. Severalexamples are presented in order to demonstrate the use of theformulations developed in the paper.     

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.     

Analysis of Thin Beams and Cables Using the Absolute Nodal Co-ordinate Formulation

The purpose of this paper is to present formulations for beam elements based on the absolute nodal co-ordinate formulation that can be effectively and efficiently used in the case of thin structural applications. The numerically stiff behaviour resulting from shear terms in existing absolute nodal co-ordinate formulation beam elements that employ the continuum mechanics approach to formulate the elastic forces and the resulting locking phenomenon make these elements less attractive for slender stiff structures. In this investigation, additional shape functions are introduced for an existing spatial absolute nodal co-ordinate formulation beam element in order to obtain higher accuracy when the continuum mechanics approach is used to formulate the elastic forces. For thin structures where bending stiffness can be important in some applications, a lower order cable element is introduced and the performance of this cable element is evaluated by comparing it with existing formulations using several examples. Cables that experience low tension or catenary systems where bending stiffness has an effect on the wave propagation are examples in which the low order cable element can be used. The cable element, which does not have torsional stiffness, can be effectively used in many problems such as in the formulation of the sliding joints in applications such as the spatial pantograph/catenary systems. The numerical study presented in this paper shows that the use of existing implicit time integration methods enables the simulation of multibody systems with a moderate number of thin and stiff finite elements in reasonable CPU time.     

基于压电曲壳单元的结构振动最优控制

曲壳结构已广泛应用于航空航天、航海等领域,结构微小的振动会极大地影响部件性能,抑制曲壳结构的振动就很必要。提出了压电曲壳结构的振动最优控制模型,利用空间曲壳单元理论及多点约束方程实现曲壳基结构壳元和压电壳元的耦合,推导了压电曲壳结构动力学有限元方程,并结合线性二次型最优控制理论,实现了压电作动器对曲壳结构的振动最优控制。数值算例表明,对曲壳结构进行动力分析和振动控制时,在达到同样精度及控制效果的情况下,与平壳相比,曲壳单元及作动器所需数目要少得多。     

高速碰撞下圆柱壳自由梁的穿孔特性

利用二级轻气炮加载,进行了球状2A12铝弹丸垂直撞击圆柱壳自由梁实验。并进行了弹丸速度、圆柱壳直径和壁厚等因素对穿孔直径影响的数值模拟,数值模拟结果和实验结果基本吻合。通过量纲分析和数值模拟结合,推导了穿孔直径与相关影响参数的经验关系式。研究结果表明：当圆柱壳直径和厚度不变时,高速撞击产生的穿孔直径在径向和轴向都随着弹丸速度增大而增大;当弹丸速度和圆柱壳厚度不变时,高速撞击产生的穿孔直径随着圆柱壳自由梁直径的增大而减小;当弹丸速度和圆柱壳直径不变时,穿孔直径随着圆柱壳厚度的增大而减小。

     

空间曲梁单元的扭转修正系数研究

三维空间曲梁有限单元模型是模拟曲梁结构的有效数值方法,可以考虑曲梁的弯扭耦合特性,最为符合曲梁的几何和受力特征.由于有限元法采用梁理论的平截面假定,空间曲梁单元上的扭转剪应力分布与实际曲梁截面上的扭转剪应力不同,从而会导致扭转刚度和扭转变形的计算失真.本文基于剪切应变能等效原理,推导了不同长宽比的矩形截面空间曲梁单元的...     

Parametric vibrations of a ribbed cylindrical shell with an attached beam

Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Prikladnaya Mekhanika, Vol. 24, No. 8, pp 47–55, August, 1988.     

Limit analysis of masonry vaults by means of curved shell finite elements and homogenization

The study of masonry vaults should take into account the essentials of the material “masonry” – i.e. heterogeneity, almost no resistance to tension combined with a good compressive strength and a high friction coefficient, as well as the overall importance of the geometry for achieving the equilibrium.In this paper, a new six-noded triangular curved element, specifically developed for the kinematic limit analysis of masonry shells, is presented. Plastic dissipation is allowed only at the interfaces (generalized cylindrical hinges) between adjoining elements for combined membrane actions, bending moment, torsion and out-of-plane shear, as it is required for the analysis of thick (Reissner–Mindlin) shells. An upper bound of the collapse load is so obtained, since, looking at the dual formulation, the admissibility of the stress state is imposed only at the element boundaries. Masonry strength domain at each interface between contiguous triangular elements is evaluated resorting to a suitable upper bound FE homogenization procedure. The model is assessed through several numerical simulations on a number of masonry shells experimentally tested until collapse. In particular, the dependence of the collapse multiplier on the mesh and on the material parameters (sensitivity analysis) is thoroughly discussed.