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

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

On the use of absolute interface coordinates in the floating frame of reference formulation for flexible multibody dynamics

In this work a new formulation for flexible multibody systems is presented based on the floating frame formulation. In this method, the absolute interface coordinates are used as degrees of freedom. To this end, a coordinate transformation is established from the absolute floating frame coordinates and the local interface coordinates to the absolute interface coordinates. This is done by assuming linear theory of elasticity for a body’s local elastic deformation and by using the Craig–Bampton interface modes as local shape functions. In order to put this new method into perspective, relevant relations between inertial frame, corotational frame and floating frame formulations are explained. As such, this work provides a clear overview of how these three well-known and apparently different flexible multibody methods are related. An advantage of the method presented in this work is that the resulting equations of motion are of the differential rather than the differential-algebraic type. At the same time, it is possible to use well-developed model order reduction techniques on the flexible bodies locally. Hence, the method can be employed to construct superelements from arbitrarily shaped three dimensional elastic bodies, which can be used in a flexible multibody dynamics simulation. The method is validated by simulating the static and dynamic behavior of a number of flexible systems.     

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.     

System reduction in multibody dynamics of wind turbines

A system reduction scheme is devised related to a multibody formulation from which the dynamic response of a wind turbine is determined. In this formulation, each substructure is described in its own frame of reference, which is moving freely in the vicinity of the moving substructure. The Ritz bases spanning the reduced system comprises of rigid body modes and some dynamic low-frequency elastic eigenmodes compatible to the kinematic constraints of the related substructure. The high-frequency elastic modes are presumed to cause merely quasi-static displacements, and thus are included in the expansion via a quasi-static correction. The results show that by using the derived reduction scheme it is only necessary with 2 dynamical modes for the blade substructure when the remaining modes are treated as quasi-static. Moreover, it is shown that it has little to none effect if the gyroscopic stiffness matrix during a stopped situation or under nominal operational conditions is used to derive the functional basis of the modal expansion.     

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

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

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.     

Topology optimization based on level set for a flexible multibody system modeled via ANCF

A topology optimization methodology is proposed for the flexible multibody system undergoing both large overall motion and large deformation. The system of concern is modeled via the absolute nodal coordinate formulation. The equivalent static load method is employed to transform the topology optimization of the nonlinear dynamic response of the system into a static one, and evaluated to adapt to the absolute nodal coordinate formulation by splitting the elastic deformations of the flexible components from the overall motions of those components. During the static topology optimization, the material interface is implicitly described as the zero level set of a higher-dimensional scalar function. Then, the semi-implicit level set method with the additive operator splitting algorithm is employed to solve the corresponding Hamilton-Jacobi partial differential equation. In addition, the expert evaluation method of weights based on the grey theory is utilized to define the objective function, and a modified augmented Lagrange multiplier method is proposed to treat the inequality volume constraint so as to avoid the oscillation and drift of the volume. Finally, two numerical examples are provided to validate the proposed methodology.     

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.     

Structural sensitivity analysis of flexible multibody systems modeled with the floating frame of reference approach using the adjoint variable method

For the efficient analysis and optimization of flexible multibody systems, gradient information is often required. Next to simple and easy-to-implement finite difference approaches, analytical methods, such as the adjoint variable method, have been developed and are now well established for the sensitivity analysis in multibody dynamics. They allow the computation of exact gradients and require normally less computational effort for large-scale problems. In the current work, we apply the adjoint variable method to flexible multibody systems with kinematic loops, which are modeled using the floating frame of reference formulation. Thereby, in order to solve ordinary differential equations only, the equations of motion are brought into minimal form using coordinate partitioning, and the constraint equations at position and velocity level are incorporated in the adjoint dynamics. For testing and illustrative purposes, the procedure is applied to compute the structural gradient for a flexible piston rod of a slider–crank mechanism.     

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.     

Flexible Multibody Dynamics: Review of Past and Recent Developments

In this paper, a review of past and recent developments in the dynamics of flexible multibody systems is presented. The objective is to review some of the basic approaches used in the computer aided kinematic and dynamic analysis of flexible mechanical systems, and to identify future directions in this research area. Among the formulations reviewed in this paper are the floating frame of reference formulation, the finite element incremental methods, large rotation vector formulations, the finite segment method, and the linear theory of elastodynamics. Linearization of the flexible multibody equations that results from the use of the incremental finite element formulations is discussed. Because of space limitations, it is impossible to list all the contributions made in this important area. The reader, however, can find more references by consulting the list of articles and books cited at the end of the paper. Furthermore, the numerical procedures used for solving the differential and algebraic equations of flexible multibody systems are not discussed in this paper since these procedures are similar to the techniques used in rigid body dynamics. More details about these numerical procedures as well as the roots and perspectives of multibody system dynamics are discussed in a companion review by Schiehlen [79]. Future research areas in flexible multibody dynamics are identified as establishing the relationship between different formulations, contact and impact dynamics, control-structure interaction, use of modal identification and experimental methods in flexible multibody simulations, application of flexible multibody techniques to computer graphics, numerical issues, and large deformation problem. Establishing the relationship between different flexible multibody formulations is an important issue since there is a need to clearly define the assumptions and approximations underlying each formulation. This will allow us to establish guidelines and criteria that define the limitations of each approach used in flexible multibody dynamics. This task can now be accomplished by using the absolute nodal coordinate formulation which was recently introduced for the large deformation analysis of flexible multibody systems.     

Geometric nonlinear effects on the planar dynamics of a pivoted flexible beam encountering a point-surface impact

Flexible-body modeling with geometric nonlinearities remains a hot topic of research by applications in multibody system dynamics undergoing large overall motions. However, the geometric nonlinear effects on the impact dynamics of flexible multibody systems have attracted significantly less attention. In this paper, a point-surface impact problem between a rigid ball and a pivoted flexible beam is investigated. The Hertzian contact law is used to describe the impact process, and the dynamic equations are formulated in the floating frame of reference using the assumed mode method. The two important geometric nonlinear effects of the flexible beam are taken into account, i.e., the longitudinal foreshortening effect due to the transverse deformation, and the stress stiffness effect due to the axial force. The simulation results show that good consistency can be obtained with the nonlinear finite element program ABAQUS/Explicit if proper geometric nonlinearities are included in the floating frame formulation. Specifically, only the foreshortening effect should be considered in a pure transverse impact for efficiency, while the stress stiffness effect should be further considered in an oblique case with much more computational effort. It also implies that the geometric nonlinear effects should be considered properly in the impact dynamic analysis of more general flexible multibody systems.     

Validation of flexible multibody dynamics beam formulations using benchmark problems

As the need to model flexibility arose in multibody dynamics, the floating frame of reference formulation was developed, but this approach can yield inaccurate results when elastic displacements becomes large. While the use of three-dimensional finite element formulations overcomes this problem, the associated computational cost is overwhelming. Consequently, beam models, which are one-dimensional approximations of three-dimensional elasticity, have become the workhorse of many flexible multibody dynamics codes. Numerous beam formulations have been proposed, such as the geometrically exact beam formulation or the absolute nodal coordinate formulation, to name just two. New solution strategies have been investigated as well, including the intrinsic beam formulation or the DAE approach. This paper provides a systematic comparison of these various approaches, which will be assessed by comparing their predictions for four benchmark problems. The first problem is the Princeton beam experiment, a study of the static large displacement and rotation behavior of a simple cantilevered beam under a gravity tip load. The second problem, the four-bar mechanism, focuses on a flexible mechanism involving beams and revolute joints. The third problem investigates the behavior of a beam bent in its plane of greatest flexural rigidity, resulting in lateral buckling when a critical value of the transverse load is reached. The last problem investigates the dynamic stability of a rotating shaft. The predictions of eight independent codes are compared for these four benchmark problems and are found to be in close agreement with each other and with experimental measurements, when available.     

On the influence of model reduction techniques in topology optimization of flexible multibody systems using the floating frame of reference approach

Employing the floating frame of reference formulation in the topology optimization of dynamically loaded components of flexible multibody systems seems to be a natural choice. In this formulation the deformation of flexible bodies is approximated by global shape functions, which are commonly obtained from finite element models using model reduction techniques. For topology optimization these finite element models can be parameterized using the solid isotropic material with penalization (SIMP) approach. However, little is known about the interplay of model reduction and SIMP parameterization. Also securing the model reduction quality despite major changes of the design during the optimization has not been addressed yet. Thus, using the examples of a flexible frame and a slider-crank mechanism this work discusses the proper choice of the model reduction technique in the topology optimization of flexible multibody systems.     

A method for improving dynamic solutions in flexible multibody dynamics

This paper presents a method for improving dynamic solutions that are obtained from the dynamic simulation of flexible multibody systems. The mode-acceleration concept in linear structural dynamics is utilized in the proposed method for improving accuracy in the postprocessing stage. A theoretical explanation is made on why the proposed method improves the dynamic solutions in the context of the mode-acceleration method. A mode-acceleration equation for each flexible body is defined and the load term in the right hand side of the equation is represented as a combination of space-dependent and time-dependent terms so that efficient computation of dynamic solutions can be achieved. The load term is obtained from dynamic simulation of a flexible multibody system and a finite element method is used to compute dynamic solutions by quasi-static analyses. Numerical examples show the effectiveness of the proposed method.     

Optimization of flexible components of multibody systems via equivalent static loads

An optimization methodology that iteratively links the results of multibody dynamics and structural analysis software to an optimization method is presented to design flexible multibody systems under dynamic loading conditions. In particular, rigid multibody dynamic analysis is utilized to calculate dynamic loads of a multibody system and a structural optimization algorithm using equivalent static loads transformed from the dynamic loads are used to design the flexible components in the multibody dynamic system. The equivalent static loads, which are derived from equations of motion, are used as multiple loading conditions of linear structural optimization. A simple example is solved to verify the proposed methodology and the pelvis part of the biped humanoid, a complex multibody system which consists of many bodies and joints, is redesigned using the proposed methodology.     

Weakly and fully coupled methods for structural optimization of flexible mechanisms

The paper concerns a detailed comparison between two optimization methods that are used to perform the structural optimization of flexible components within a multibody system (MBS) simulation. The dynamic analysis of flexible MBS is based on a nonlinear finite element formulation. The first method is a weakly coupled method, which reformulates the dynamic response optimization problem in a two-level approach. First, a rigid or flexible MBS simulation is performed, and second, each component is optimized independently using a quasi-static approach in which a series of equivalent static load (ESL) cases obtained from the MBS simulation are applied to the respective components. The second method, the fully coupled method, performs the dynamic response optimization using the time response obtained directly from the flexible MBS simulation. Here, an original procedure is proposed to evaluate the ESL from a nonlinear finite element simulation, contrasting with the floating reference frame formulation exploited in the standard ESL method. Several numerical examples are provided to support our position. It is shown that the fully coupled method is more general and accommodates all types of constraints at the price of a more complex optimization process.     

Analysis of multibody systems with flexible plates using variational graph-theoretic methods

This work treats the problem of modelling multibody systems with structural flexibility. By combining linear graph theory with the principle of virtual work and finite elements, a dynamic formulation is obtained that extends graph-theoretic (GT) modelling methods to the analysis of thin flexible plates for multibody systems. The system is represented by a linear graph, in which nodes represent reference frames on flexible plates, and edges represent components that connect these frames. To generate the equations of motion with elastic deformations, the flexible plates are discretized using a triangular thin shell finite element based on the discrete Kirchhoff criterion and can be used to discretize bidirectional bodies such as satellite panels, flatbed trailers, and mechanisms with plates. Three flexible systems with plates are analyzed to illustrate the performance of this new variational graph-theoretic formulation and its ability to generate directly a set of motion equations for flexible multibody systems (FMS) without additional user input.     

Cable installation simulation by using a multibody dynamic model

A major concern when installing the cables into the underground conduit is minimizing the tensile forces exerted on the cables as they are pulled. This knowledge makes it possible to avoid over conservative design practices and to achieve substantial saving during construction. A general computing algorithm of predicting the tensile force of the cable pulled through the underground conduit with an arbitrary configuration is presented in this paper, which is based on multibody system dynamic formulation. The presented multibody dynamic model for this problem consists of the cable, the underground conduit, and the interaction between the cable and the conduit. In this paper, the cable is modeled by the finite cable element based on an absolute nodal coordinate formulation. The interaction between the cable and the underground conduit is described by the Hertz contact theory. Numerical examples are presented to illustrate the effectiveness and efficiency of the proposed method for estimating the cable tension.     

A systematic method of dynamics for flexible robot manipulators

In this article, a systematic method to derive dynamic equations of motion for flexible robot manipulators is developed by using the Lagrangian assumed modes method. The proposed method can be applied to dynamic simulation and control system design for flexible robot manipulators. In the proposed method, the link deflection is described by a truncated modal expansion. The operations of only 3x3 matrices and/or 3 × 1 vectors exist in the method. All the dynamics computations are performed in the link coordinate systems, where the kinematics informations are computed with the forward recursion from the base to the hand tip and the dynamics informations are computed with the return recursion. As generally compared with other existing methods, the method proposed in this article is, computationally, more simple, systematic, and efficient. A computational simulation for a single-link flexible robot manipulator is presented to verify the proposed method. © 1992 John Wiley & Sons, Inc.