Quasi-momentum theorem in Riemann-Cartan space

The geometric formulation of motion of the first-order linear homogenous scleronomous nonholonomic system subjected to active forces is studied with the nonholonomic mapping theory. The quasi-Newton law, the quasi-momentum theorem, and the second kind Lagrange equation of dynamical systems are obtained in the Riemann-Cartan configuration spaces. By the nonholonomic mapping, a Euclidean configuration space or a Riemann configuration space of a dynamical system can be mapped into a Riemann-Cartan configuration space with torsion. The differential equations of motion of the dynamical system can be obtained in its Riemann-Cartan configuration space by the quasi-Newton law or the quasi-momentum theorem. For a constrained system, the differential equations of motion in its Riemann-Cartan configuration space may be simpler than the equations in its Euclidean configuration space or its Riemann configuration space. Therefore, the nonholonomic mapping theory can solve some constrained problems, which are difficult to be solved by the traditional analytical mechanics method. Three examples are given to illustrate the effectiveness of the method.     

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.     

Motion equations in redundant coordinates with application to inverse dynamics of constrained mechanical systems

The basis for any model-based control of dynamical systems is a numerically efficient formulation of the motion equations, preferably expressed in terms of a minimal set of independent coordinates. To this end the coordinates of a constrained system are commonly split into a set of dependent and independent ones. The drawback of such coordinate partitioning is that the splitting is not globally valid since an atlas of local charts is required to globally parameterize the configuration space. Therefore different formulations in redundant coordinates have been proposed. They usually involve the inverse of the mass matrix and are computationally rather complex. In this paper an efficient formulation of the motion equations in redundant coordinates is presented for general non-holonomic systems that is valid in any regular configuration. This gives rise to a globally valid system of redundant differential equations. It is tailored for solving the inverse dynamics problem, and an explicit inverse dynamics solution is presented for general full-actuated systems. Moreover, the proposed formulation gives rise to a non-redundant system of motion equations for non-redundantly full-actuated systems that do not exhibit input singularities.     

A constrained theory of a Cosserat point for the numerical solution of dynamic problems of non-linear elastic rods with rigid cross-sections

A constrained theory of a Cosserat point has been developed for the numerical solution of non-linear elastic rods. The cross-sections of the rod element are constrained to remain rigid but tangential shear deformations and axial extension are admitted. As opposed to the more general theory with deformable cross-sections, the kinetic coupling equations in the numerical formulation of the constrained theory are expressed in terms of the simple physical quantities of force and mechanical moment applied to the common ends of neighboring elements. Also, in contrast with standard finite element methods, the Cosserat element uses a direct approach to the development of constitutive equations. Specifically the kinetic quantities are determined by algebraic expressions which are obtained by derivatives of a strain energy function. Most importantly, no integration is needed over the element region. A number of example problems have been considered which indicate that the constrained Cosserat element can be used to model large deformation dynamic response of non-linear elastic rods.     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

Flexible Multibody Simulation and Choice of Shape Functions

The approach most widely used for the modelling of flexible bodies in multibody systems has been called the floating frame of reference formulation. In this methodology the flexible body motion is subdivided into a reference motion and deformation. The displacement field due to deformation is approximated by the Ritz method as a product of known shape functions and unknown coordinates depending on time only. The shape functions may be obtained using finite-element-models of flexible bodies in multibody systems, resulting in a detailed system representation and a high number of system equations. The number of system equations of such a nodal approach can be reduced considerably using a modal representation of deformation. This modal approach, however, leads to the fundamental problem of selecting the shape functions.The floating frame of reference formulation is reviewed here using a generic flexible body model, from which the various body models used in multibody simulations may be derived by formulation of specific constraint equations. Special attention is given in this investigation to the following subjects: The separation of flexible body motion into a reference motion and deformation requires the definition of a body reference frame, which in turn affects the choice of shape functions. Some alternatives will be outlined together with their advantages and disadvantages. Assuming the body deformation to be small, the system equations can be linearized. This may require considering geometric stiffening terms. The problem of how to compute these terms has been solved in literature on the instability of structures under critical loads. For finite element models the geometric stiffening terms are obtained from the tangential stiffness matrix. The generality of the flexible body model allows the definition of an object oriented data base to describe the system bodies. Such a data base includes a general interface between multibody- and finite-element-codes. By combining eigenfunctions and static deformation modes to represent body deformation one obtains a set of so-called quasi-comparison functions. When selected properly these functions can be shown to improve the representation of stresses significantly.     

多体动力学的几何积分方法研究进展

动力系统的几何积分研究是近20年来工程计算领域非常活跃的方向.多体动力学方程(微分方程, 微分代数方程)是一类典型的动力系统,将其从Lagrange体系向Hamilton系统过渡,目的在于从欧氏几何过渡到辛几何形态, 将对偶变量引入到力学研究中,然后利用辛几何的数学框架对多体系统动力学方程进行数值计算,可以预知多体动力学系统的一些定性信息,并在数值离散时能保持这些定性性质特征,尤其在表示关键的物理意义时需要强调保持这些几何性质.简要介绍多体系统(无约束多刚体系统、完整约束多刚体系统和柔性多体系统)的Hamilton正则方程的建立和几何积分方法的构造,着重介绍了在多体动力学计算中非常有应用前景的高阶辛算法(合成辛算法、分裂合成辛算法和辛精细积分法)、多辛算法,以及广义Hamilton 系统与Lie 群积分方法等计算几何力学方法, 并对Lie群积分的投影方法、流形局部坐标法等方法进行了阐述.      

Finite element modeling of contact conditions in multibody system dynamics

In this paper a new method is developed for the dynamic analysis of contact conditions in flexible multibody systems undergoing a rolling type of motion. The relative motion between the two contacting bodies is treated as a constraint condition describing their kinematic and geometric relations. Equations of motion of the system are presented in a matrix form making use of Kane's equations and finite element method. The method developed has been implemented in a general purpose program called DARS and applied to the simulation and analysis of a rotating wheel on a track. Both the bodies are assumed flexible and discretized using a three dimensional 8-noded isoparametric elements. The time variant constraint conditions are imposed on the nodal points located at the peripheral surfaces of the bodies under consideration. The simulation is carried out under two different boundary conditions describing the support of the track. The subsequent constraint forces associated with the generalized coordinates of the system are computed and plotted. The effects of friction are also discussed.     

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     

Dynamics of Articulated Structures. Part I. Theory

ABSTRACT

A finite element based method is developed for geometrically nonlinear dynamic analysis of spatial articulated structures; i.e., structures in which kinematic connections permit large relative displacement between components that undergo small elastic deformation. Vibration and static correction modes are used to account for linear elastic deformation of components. Kinematic constraints between components are used to define boundary conditions for vibration analysis and loads for static correction mode analysis. Constraint equations between flexible bodies are derived in a systematic way and a Lagrange multiplier formulation is used to generate the coupled large displacement-small deformation equations of motion. A lumped mass finite element structural analysis formulation is used to generate deformation modes. An intermediate-processor is used to calculate time-independent terms in the equations of motion and to generate input data for a large-scale dynamic analysis code that includes coupled effects of geometric nonlinearity and elastic deformation. Examples are presented and the effects of deformation mode selection on dynamic prediction are analyzed in Part II of the paper.     

Dynamic Analysis of a Flexible Vehicle Structure with Nonlinear Force Elements Using Substructure Technique∗

ABSTRACT

The substnicturing technique is used to analyze the motion of a flexible vehicle structure with nonlinear force elements. The number of effective degrees of freedom of the combined system model is greatly reduced by synthesizing the individually modeled substructures, using constraints of geometric compatibility at the interfaces. When combining the substructures, recursive-type nonlinear integral equations (NIE) are introduced instead of the conventional nonlinear differential equations (NDE). It is shown that the NIE formulation is computationally more efficient than the NDE formulation in simulating steady-state and transient responses of a flexible vehicle with nonlinear dampers. The NIE formulation is further applied to the nonlinear damping optimization problem, and it is found that this method is efficient for optimization.     

充液系统液体-多体耦合动力响应分析

提出了充液系统的液体-多体耦合力学模型，基于ALE有限元法和多体系统动力学理论，发展了液体-多体耦合动力响应分析的一种有效方法. 对于液体子系统，将其运动分解为随同贮箱的大位移运动和相对贮箱的大幅晃动，引入贮箱固连参考系中的任意拉格朗日-欧拉(ALE)运动学描述，建立了贮箱固连非惯性参考系中液体的ALE有限元方程，对液体有限元方程的缩聚大大减少了液体子系统的计算规模. 为了计及液体阻尼的影响，引入了液体修正的Rayleigh阻尼，避免了伪阻尼力的出现. 对于多体子系统，应用多体系统动力学理论建立动力学方程. 在此基础上详细导出了液体-多体耦合动力学方程，并采用预估-多重校正算法(PMA)和时间步长控制算法进行迭代求解，既保证了迭代收敛，又提高了计算效率. 所给算例成功求解了液体运输车辆系统的液体-多体耦合动力响应，深入分析了有关参数对系统动力响应的影响，获得了一些结论.     

非完整约束系统几何动力学研究进展:Lagrange理论及其它

近10年来, 非完整力学的发展主要集中在两个相互关联的方向上, 一个是非完整运动规划, 另一个则是非完整约束系统的几何动力学, 这两个研究方向都充分地利用了现代几何学, 如纤维丛理论、辛流形和Poisson流形结构等等.本文主要综述非完整约束系统几何动力学的外附型和内禀型Lagrange理论, 包括非定常力学系统所需要的射丛几何学的基本概念、射丛按约束的直和分解、约束流形上的水平分布、D'Alembert-Lagrange方程与Chaplygin方程的整体描述、以及Riemann-Cartan流形上的非完整力学, 文中对Chetaev条件和d-δ交换关系的几何意义作了深入讨论.除此之外, 简要评述非完整力学的Hamilton理论与赝Poisson结构、Noether对称性和Lie对称性、动量映射与对称约化、Vakonomic动力学等几个非常重要专题的研究进展.      

Finite element,stream function-vorticity solution of secondary flow in the channels of a rod bundle

A finite element method has been applied to predict the overall features of the fully developed turbulent flow in the non-circular channels of a rod bundle. The finite element discretization is based on the conventional Galerkin method using an isoparametric quadrilateral element with mixed interpolation. The primary axial flow and turbulent kinetic energy distributions have been predicted for fully developed turbulent flow conditions right up to the wall. The secondary velocity is represented by the stream function-vorticity formulation and the no-slip boundary conditions are explicitly introduced in the nonlinear equations by a boundary vorticity formula. The Newton-Raphson method is applied to the stream function-vorticity equations and solved simultaneously by the frontal solution technique. A one-equation eddy viscosity model of turbulence and an algebraic stress transport model have been used to predict primary axial velocity, secondary velocities and turbulent kinetic energy. The predictions obtained for a central subchannel of an equilateral-triangular rod array with p/d= 1.3 are in reasonable agreement with experimental data.     

Multi-Multiplier Ambient-Space Formulations of Constrained Dynamical Systems, with an Application to Elastodynamics

Various formulations of the equations of motion for both finite- and infinite-dimensional constrained Lagrangian dynamical systems are studied. The different formulations correspond to different ways of enforcing constraints through multiplier fields. All the formulations considered are posed on ambient spaces whose members are unrestricted by the need to satisfy constraint equations, but each formulation is shown to possess an invariant set on which the constraint equations and physical balance laws are satisfied. The stability properties of the invariant set within its ambient space are shown to be different in each case. We use the specific model problem of linearized incompressible elastodynamics to compare properties of three different ambient-space formulations. We establish the well-posedness of one formulation in the particular case of a homogeneous, isotropic body subject to specified tractions on its boundary. Accepted October 11, 2000?Published online April 23, 2001     

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.     

Numerical investigation of structural geometric nonlinearity effect in high-aspect-ratio wing using CFD/CSD coupled approach

An efficient and robust fluid–structure coupled methodology has been developed to investigate the linear and non-linear static aeroelastic behavior of flexible high-aspect-ratio wing. A three-dimensional open source finite element solver has been loosely coupled with an in-house Reynolds-averaged Navier–Stokes solver, designed for hybrid-unstructured meshes, to perform aero-structural coupled simulations. For volume mesh deformation and two-way data interpolation over non-matching grids interface, a radial basis function methodology combined with a data reduction algorithm has been used. This technique is efficient in handling large deflections and provides high-quality deformed meshes. Structural geometric nonlinearity has been considered to predict the deformations in the vertical and torsional directions caused by gravitational and aerodynamic loading. A multi-material finite element model has been generated to match the experimental configuration. Computational aeroelastic simulations were performed on an experimental high-aspect-ratio aeroelastic wing model with a slender body at the tip to get non-linear static deflections, twist and structure natural frequencies. The effect of the geometric nonlinearity is significant for large deformation analysis and has been highlighted in the predicted maximum tip deflection and twist. Good qualitative and quantitative agreement has been achieved between the predicted results and the available experimental data.     

Numerical Analysis of Coupled Stokes/Darcy Flows in Industrial Filtrations

Free flow channel confined by porous walls is a feature of many of the natural and industrial settings. Viscous flows adjacent to saturated porous medium occur in cross-flow and dead-end filtrations employed primarily in pharmaceutical and chemical industries for solid–liquid or gas–solid separations. Various mathematical models have been put forward to describe the conjugate flow dynamics based on theoretical grounds and experimental evidence. Despite this fact, there still exists a wide scope for extensive research in numerical solutions of these coupled models when applied to problems with industrial relevance. The present work aims towards the numerical analysis of coupled free/porous flow dynamics in the context of industrial filtration systems. The free flow dynamics has been expressed by the Stokes equations for the creeping, laminar flow regime whereas the flow behaviour in very low permeability porous media has been represented by the conventional Darcy equation. The combined free/porous fluid dynamical behaviour has been simulated using a mixed finite element formulation based on the standard Galerkin technique. A nodal replacement technique has been developed for the direct linking of Stokes and Darcy flow regimes which alleviates specification of any additional constraint at the free/porous interface. The simulated flow and pressure fields have been found for flow domains with different geometries which represent prototypes of actual industrial filtration equipment. Results have been obtained for varying values of permeability of the porous medium for generalised Newtonian fluids obeying the power law model. A series of numerical experiments has been performed in order to validate the coupled flow model. The developed model has been examined for its flexibility in dealing with complex geometrical domains and found to be generic in delivering convergent, stable and theoretically consistent results. The validity and accuracy of the simulated results has been affirmed by comparing with available experimental data.