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

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

A new locking-free formulation for planar,shear deformable,linear and quadratic beam finite elements based on the absolute nodal coordinate formulation

Many widely used beam finite element formulations are based either on Reissner’s classical nonlinear rod theory or the absolute nodal coordinate formulation (ANCF). Advantages of the second method have been pointed out by several authors; among the benefits are the constant mass matrix of ANCF elements, the isoparametric approach and the existence of a consistent displacement field along the whole cross section. Consistency of the displacement field allows simpler, alternative formulations for contact problems or inelastic materials. Despite conceptional differences of the two formulations, the two models are unified in the present paper.     

Generalization of Plate Finite Elements for Absolute Nodal Coordinate Formulation

We propose a way to generate new finite elements in the absolute nodalcoordinate formulation (ANCF) and use a generalization of displacementfields and degrees of freedom (d.o.f.) of ordinary finite elements usedin structural mechanics. Application of this approach to 16- and12-d.o.f. rectangle plate elements as well as to 9-d.o.f. triangleelement gives, accordingly, 48-, 36- and 27-d.o.f. ANCF plate elements.We perform a thorough study of a 48-d.o.f. Hermitian element. Its shapefunction set is a Cartesian product of sets of one-dimensional shapefunctions for beam elements. Arguments of the shape functions aredecoupled, that is why an explicit calculation of terms of equations ofmotion leads to single integration only. We develop several models ofelastic forces of different complexity with their Jacobian matrices.Convergence and accuracy of the finite element is demonstrated ingeometrically nonlinear static and dynamic test problems, as well as inlinear analysis of natural frequencies.     

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.     

A formal procedure and invariants of a transition from conventional finite elements to the absolute nodal coordinate formulation

In this paper, finite elements based on the absolute nodal coordinate formulation (ANCF) are studied. The formulation has been developed by various authors for the dynamical simulation of large-displacement and large-rotation problems in flexible multibody dynamics. This study introduces a procedure to track the general geometrical properties of ANCF elements back to their prototypes in the conventional finite-element method (FEM), which deals with small-displacement problems. In this study, it is shown that each known ANCF element can be derived from a conventional FEM using a universal transform. Moreover, some important static and dynamic properties of the elements in small-displacement problems are automatically preserved. In the past, the authors of each newly proposed ANCF element have made unnecessary efforts to show the consistency of the above mentioned properties.     

A study of moderately thick quadrilateral plate elements based on the absolute nodal coordinate formulation

Finite element analysis using plate elements based on the absolute nodal coordinate formulation (ANCF) can predict the behaviors of moderately thick plates subject to large deformation. However, the formulation is subject to numerical locking, which compromises results. This study was designed to investigate and develop techniques to prevent or mitigate numerical locking phenomena. Three different ANCF plate element types were examined. The first is the original fully parameterized quadrilateral ANCF plate element. The second is an update to this element that linearly interpolates transverse shear strains to overcome slow convergence due to transverse shear locking. Finally, the third is based on a new higher order ANCF plate element that is being introduced here. The higher order plate element makes it possible to describe a higher than first-order transverse displacement field to prevent Poisson thickness locking. The term “higher order” is used, because some nodal coordinates of the new plate element are defined by higher order derivatives. The performance of each plate element type was tested by (1) solving a comprehensive set of small deformation static problems, (2) carrying out eigenfrequency analyses, and (3) analyzing a typical dynamic scenario. The numerical calculations were made using MATLAB. The results of the static and eigenfrequency analyses were benchmarked to reference solutions provided by the commercially available finite element software ANSYS. The results show that shear locking is strongly dependent on material thickness. Poisson thickness locking is independent of thickness, but strongly depends on the Poisson effect. Poisson thickness locking becomes a problem for both of the fully parameterized element types implemented with full 3-D elasticity. Their converged results differ by about 18 % from the ANSYS results. Corresponding results for the new higher order ANCF plate element agree with the benchmark. ANCF plate elements can describe the trapezoidal mode; therefore, they do not suffer from Poisson locking, a reported problem for fully parameterized ANCF beam elements. For cases with shear deformation loading, shear locking slows solution convergence for models based on either the original fully parameterized plate element or the newly introduced higher order plate element.     

Curve-induced distortion of polynomial space curves,flat-mapped extension modeling,and their impact on ANCF thin-plate finite elements

The Absolute Nodal Coordinate Formulation (ANCF) is a relatively new nonlinear finite element type that uses Hermite splines for shape functions. In this investigation, the ANCF is examined as a possible tool for use in modeling the media in flexible media transport systems, such as printers, copy machines, and roll-to-roll systems. However, it is demonstrated using an example of a thin plate-type ANCF finite element that these elements can suffer from significant membrane locking, which can be problematic for paper or paper-like media. One source of this locking is identified to be a property of all parametric curves that are composed of polynomials. The property is that for parametric polynomial curves, changes in the state of curvature of the curve cause changes in the distribution of points along the curve. This property is labeled Curve-Induced Distortion (CID) by the authors of this paper. CID can cause axial and membrane strain distortion in elements, causing them to be overly stiff. A new solution method is proposed to directly counteract CID in finite elements that use cubic Hermite curves for shape functions, specifically for modeling problems in which bending occurs primarily around one axis, such as paper in printing and media transport machinery. This method is labeled Flat-Mapped Extension Modeling (FMEM). FMEM is a mixed field method that uses a 1D Hermite polynomial kinematically linked to the 3D Hermite curve to represent the axial displacement field. FMEM significantly reduces the effect of CID in the ANCF element tested here. This investigation demonstrates using a single ANCF plate element type that the ANCF’s accuracy can be significantly improved by FMEM with only a small increase in computational cost. It is shown with this plate-element example that without correcting CID, the ANCF element tested is computationally much slower than contemporary methods like the co-rotational formulation for similar accuracy. But with FMEM, the ANCF is significantly faster than the co-rotational formulation for similar accuracy.     

A simple method to impose rotations and concentrated moments on ANC beams

Recently introduced ANC beam elements furnish a simple formulation that allows to solve nonlinear problems of beams, including those with large displacements and strains, as well as complex nonlinear (inelastic) materials. The success and simplicity of these finite elements is mainly due to the fact that the only nodal degrees of freedom that they employ are displacements, and rotations are thus completely avoided. This in turn makes it very difficult to apply concentrated moments or to impose rotations at specific nodes of a finite element mesh. In this article, we present a simple enhancement to this beam formulation that allows to apply these two types of boundary conditions in a simple manner, making ANC beam elements more versatile for both multibody and structural applications.     

Finite element analysis of lateral buckling for beam structures

The application of finite element analysis to lateral buckling problems, locating the critical points and tracing the postbifurcation path, is treated on the basis of a geometrically nonlinear formulation for a beam with small elastic strain but with possibly large rotations. The existing finite element formulations for thin beams are examined in the aspect of application to bifurcation problems, such as lateral buckling, and the choice of an appropriate rotation parameter for representing incremental or variational rotations in finite element formulations is discussed in relation to locating bifurcation points. This is illustrated through several numerical examples and followed by appropriate discussion.     

Finite element formulation for dynamics of planar flexible multi-beam system

In some previous geometric nonlinear finite element formulations, due to the use of axial displacement, the contribution of all the elements lying between the reference node of zero axial displacement and the element to the foreshortening effect should be taken into account. In this paper, a finite element formulation is proposed based on geometric nonlinear elastic theory and finite element technique. The coupling deformation terms of an arbitrary point only relate to the nodal coordinates of the element at which the point is located. Based on Hamilton principle, dynamic equations of elastic beams undergoing large overall motions are derived. To investigate the effect of coupling deformation terms on system dynamic characters and reduce the dynamic equations, a complete dynamic model and three reduced models of hub-beam are prospected. When the Cartesian deformation coordinates are adopted, the results indicate that the terms related to the coupling deformation in the inertia forces of dynamic equations have small effect on system dynamic behavior and may be neglected, whereas the terms related to coupling deformation in the elastic forces are important for system dynamic behavior and should be considered in dynamic equation. Numerical examples of the rotating beam and flexible beam system are carried out to demonstrate the accuracy and validity of this dynamic model. Furthermore, it is shown that a small number of finite elements are needed to obtain a stable solution using the present coupling finite element formulation.     

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.     

A nonlinear two-node superelement for use in flexible multibody systems

A nonlinear two-node superelement is proposed for the modeling of flexible complex-shaped links for use in multibody simulations. Assuming that the elastic deformations with respect to a corotational reference frame remain small, substructuring methods may be used to obtain reduced mass and stiffness matrices from a linear finite element model. These matrices are used in the derivation of potential and kinetic energy expressions of the nonlinear two-node superelement. By evaluating Lagrange’s equations, expressions for the internal and external forces acting on the superelement can be obtained. The inertia forces of the superelement are derived in terms of absolute nodal velocities and accelerations, which greatly simplifies the dynamic formulation. Three examples are included. The first two examples are used to validate the method by comparing the results with those obtained from nonlinear beam element solutions. We consider a benchmark simulation of the spin-up motion of a flexible beam with uniform cross-section and a similar simulation in which the beam is simultaneously excited in the out-of-plane direction. Results from both examples show good agreement with simulation results obtained using nonlinear finite beam elements. In a third example, the method is applied to an unbalanced rotating shaft, illustrating the potential of the proposed methodology for a more complex geometry.     

Slope discontinuities in the finite element absolute nodal coordinate formulation: gradient deficient elements

In this paper, the treatment of the slope discontinuities in the finite element absolute nodal coordinate formulation (ANCF) is discussed. The paper explains the fundamental problems associated with developing a constant transformation that accounts for the slope discontinuities in the case of gradient deficient ANCF finite elements. A procedure that allows for the treatment of slope discontinuities in the case of gradient deficient finite elements which do not employ full parameterization is proposed for the special case of commutative rotations. The use of the proposed procedure leads to a constant orthogonal element transformation that describes the element initial configuration. As a consequence, one obtains in the case of large deformation and commutative rotations, a constant mass matrix for the structures. In order to achieve this goal, the concept of the intermediate finite element coordinate system is invoked. The intermediate finite element coordinate system used in this investigation serves to define the element reference configuration, follows the rotation of the structure, and maintains a fixed orientation relative to the structure coordinate system. Since planar rotations are always commutative, the procedure proposed in this investigation is applicable to all planar gradient deficient ANCF finite elements.     

Flexible Multibody Systems Models Using Composite Materials Components

The use of a multibody methodology to describe the large motion of complex systems that experience structural deformations enables to represent the complete system motion, the relative kinematics between the components involved, the deformation of the structural members and the inertia coupling between the large rigid body motion and the system elastodynamics. In this work, the flexible multibody dynamics formulations of complex models are extended to include elastic components made of composite materials, which may be laminated and anisotropic. The deformation of any structural member must be elastic and linear, when described in a coordinate frame fixed to one or more material points of its domain, regardless of the complexity of its geometry. To achieve the proposed flexible multibody formulation, a finite element model for each flexible body is used. For the beam composite material elements, the sections properties are found using an asymptotic procedure that involves a two-dimensional finite element analysis of their cross-section. The equations of motion of the flexible multibody system are solved using an augmented Lagrangian formulation and the accelerations and velocities are integrated in time using a multi-step multi-order integration algorithm based on the Gear method.     

Absolute nodal coordinate plane beam formulation for multibody systems dynamics

A new plane beam dynamic formulation for constrained multibody system dynamics is developed. Flexible multibody system dynamics includes rigid body dynamics and superimposed vibratory motions. The complexity of mechanical system dynamics originates from rotational kinematics, but the natural coordinate formulation does not use rotational coordinates, so that simple dynamic formulation is possible. These methods use only translational coordinates and simple algebraic constraints. A new formulation for plane flexible multibody systems are developed utilizing the curvature of a beam and point masses. Using absolute nodal coordinates, a constant mass matrix is obtained and the elastic force becomes a nonlinear function of the nodal coordinates. In this formulation, no infinitesimal or finite rotation assumptions are used and no assumption on the magnitude of the element rotations is made. The distributed body mass and applied forces are lumped to the point masses. Closed loop mechanical systems consisting of elastic beams can be modeled without constraints since the loop closure constraints can be substituted as beam longitudinal elasticity. A curved beam is modeled automatically. Several numerical examples are presented to show the effectiveness of this method.     

First order sensitivity analysis of flexible multibody systems using absolute nodal coordinate formulation

Design sensitivity analysis of flexible multibody systems is important in optimizing the performance of mechanical systems. The choice of coordinates to describe the motion of multibody systems has a great influence on the efficiency and accuracy of both the dynamic and sensitivity analysis. In the flexible multibody system dynamics, both the floating frame of reference formulation (FFRF) and absolute nodal coordinate formulation (ANCF) are frequently utilized to describe flexibility, however, only the former has been used in design sensitivity analysis. In this article, ANCF, which has been recently developed and focuses on modeling of beams and plates in large deformation problems, is extended into design sensitivity analysis of flexible multibody systems. The Motion equations of a constrained flexible multibody system are expressed as a set of index-3 differential algebraic equations (DAEs), in which the element elastic forces are defined using nonlinear strain-displacement relations. Both the direct differentiation method and adjoint variable method are performed to do sensitivity analysis and the related dynamic and sensitivity equations are integrated with HHT-I3 algorithm. In this paper, a new method to deduce system sensitivity equations is proposed. With this approach, the system sensitivity equations are constructed by assembling the element sensitivity equations with the help of invariant matrices, which results in the advantage that the complex symbolic differentiation of the dynamic equations is avoided when the flexible multibody system model is changed. Besides that, the dynamic and sensitivity equations formed with the proposed method can be efficiently integrated using HHT-I3 method, which makes the efficiency of the direct differentiation method comparable to that of the adjoint variable method when the number of design variables is not extremely large. All these improvements greatly enhance the application value of the direct differentiation method in the engineering optimization of the ANCF-based flexible multibody systems.     

Design of a planar multibody dynamic system with ANCF beam elements based on an element-wise stiffness evaluation procedure

The absolute nodal coordinate formulation (ANCF) has been widely applied for large deformation analysis in flexible multibody dynamics. Although the formulation led to stable solutions for time integration under large rotations and deformations, excessive time consumption was recorded. The nonlinear relationship between the deformation and the internal force accounted for repeated adjustment to the force equilibrium state as the structure deformed. In this research, an equivalent model of the ANCF beam structure was constructed. The stiffness evaluation method was applied in an element-wise manner. In this model, the irrelevant parts were separated from those that relate to the displacements and design parameters enabling efficient updates of internal forces to achieve force equilibrium. Therefore, by using this model, optimization problems, in which displacements as well as design parameters keep changing can be efficiently approached. To verify the proposed method, two examples of optimization problems related to a free-falling pendulum and a slider-crank mechanism are demonstrated.     

A general three-dimensional L-section beam finite element for elastoplastic large deformation analysis

In this paper, the development of a general three-dimensional L-section beam finite element for elastoplastic large deformation analysis is presented. We propose the generalized interpolation scheme for the isoparametric formulation of three-dimensional beam finite elements and the numerical procedure is developed for elastoplastic large deformation analysis. The formulation is general and effective for other thin-walled section beam finite elements. To show the validity of the formulation proposed, a 2-node three-dimensional L-section beam finite element is implemented in an analysis code. As numerical examples, we first perform elastic small and large deformation analyses of a cantilever beam structure subjected to various tip loadings, and elastoplastic large deformation analysis of the same structure under reversed cyclic tip loading. We then analyze the failures of simply supported beam structures of different lengths and slenderness ratios under elastoplastic large deformation. The same problems are solved using refined shell finite element models of the structures. The numerical results of the L-section beam finite element developed here are compared with the solutions obtained using shell finite element analyses. We also discuss the numerical solutions in detail.