Recursive Linearization of Multibody Dynamics and Application to Control Design

1994 ◽  
Vol 116 (2) ◽  
pp. 445-451 ◽  
Author(s):  
Tsung-Chieh Lin ◽  
K. Harold Yae

The nonlinear equations of motion in multibody dynamics pose a difficult problem in linear control design. It is therefore desirable to have linearization capability in conjunction with a general-purpose multibody dynamics modeling technique. A new computational method for linearization is obtained by applying a series of first-order analytical approximations to the recursive kinematic relationships. The method has proved to be computationally more efficient. It has also turned out to be more accurate because the analytical perturbation requires matrix and vector operations by circumventing numerical differentiation and other associated numerical operations that may accumulate computational error.

Author(s):  
Tsung-Chieh Lin ◽  
K. Harold Yae

Abstract The non-linear equations of motion in multi-body dynamics pose a difficult problem in linear control design. It is therefore desirable to have linearization capability in conjunction with a general-purpose multibody dynamics modeling technique. A new computational method for linearization is obtained by applying a series of first-order analytical approximations to the recursive kinematic relationships. The method has proved to be computationally more efficient. It has also turned out to be more accurate because the analytical perturbation requires matrix and vector operations by circumventing numerical differentiation and other associated numerical operations that may accumulate computational error.


2020 ◽  
Vol 16 (3) ◽  
Author(s):  
Alejandro Cosimo ◽  
Federico J. Cavalieri ◽  
Javier Galvez ◽  
Alberto Cardona ◽  
Olivier Brüls

Abstract The aim of this work is to extend the finite element multibody dynamics approach to problems involving frictional contacts and impacts. The nonsmooth generalized-α (NSGA) scheme is adopted, which imposes bilateral and unilateral constraints both at position and velocity levels avoiding drift phenomena. This scheme can be implemented in a general purpose simulation code with limited modifications of pre-existing elements. The study of the woodpecker toy dynamics sets up a good example to show the capabilities of the NSGA scheme within the context of a general finite element framework. This example has already been studied by many authors who generally adopted a model with a minimal set of coordinates and small rotations. It is shown that good results are obtained using a general purpose finite element code for multibody dynamics, in which the equations of motion are assembled automatically and large rotations are easily taken into account. In addition, comparing results between different models of the woodpecker toy, the importance of modeling large rotations and the horizontal displacement of the woodpecker's sleeve is emphasized.


2010 ◽  
Vol 143-144 ◽  
pp. 592-597
Author(s):  
Qiu Sun Ye

An optional spatial 3D material solid reality picture may be showed on a 2D showing plane, its surface is commonly consisted of too many sub-surfaces of the 3D solid reality picture in both dispersal and continuation. It was a difficult problem to express all exact visibilities of solids sub-surfaces; because of some numerals computational error properties of visibilities are always unknown. So far, we cannot choose but select a proper dispersal degree of materials sub-surfaces, an approaching function of an optional material sub-surfaces B-Rep (Bound Representation), and numerals computation without any perturbation motion, etc. A novel numeral computational method of expressing exact visibilities of 3D spatial materials solids reality pictures with intelligent properties of VCR is given in this paper.


Author(s):  
Tamer M. Wasfy ◽  
Shahriar G. Ahmadi ◽  
Hatem M. Wasfy ◽  
Jeanne M. Peters

An explicit time integration multibody dynamics code is used to create a three-dimensional model of sand. Sand is modeled using discrete cubical particles with appropriate normal contact force and tangential friction force models. The model is used to predict the sand angle of repose and flow rate during discharge from a conical hopper. A penalty technique is used to impose normal contact constraints (including particle-particle, particle-hopper and particle-ground contact). An asperity-based friction model is used to model friction. A Cartesian Eulerian grid contact search algorithm is used to allow fast contact detection between particles. The governing equations of motion are solved along with contact constraint equations using a time-accurate explicit solution procedure. Parameter studies are performed in order to study the effects of the particle size and the orifice’s diameter of the hopper on the angle of repose and sand flow rate. The results of the simulations are validated using previously published experimental results.


Symmetry ◽  
2020 ◽  
Vol 12 (5) ◽  
pp. 773 ◽  
Author(s):  
Carmelo Militello ◽  
Leonardo Rundo ◽  
Luigi Minafra ◽  
Francesco Paolo Cammarata ◽  
Marco Calvaruso ◽  
...  

A clonogenic assay is a biological technique for calculating the Surviving Fraction (SF) that quantifies the anti-proliferative effect of treatments on cell cultures: this evaluation is often performed via manual counting of cell colony-forming units. Unfortunately, this procedure is error-prone and strongly affected by operator dependence. Besides, conventional assessment does not deal with the colony size, which is generally correlated with the delivered radiation dose or administered cytotoxic agent. Relying upon the direct proportional relationship between the Area Covered by Colony (ACC) and the colony count and size, along with the growth rate, we propose MF2C3, a novel computational method leveraging spatial Fuzzy C-Means clustering on multiple local features (i.e., entropy and standard deviation extracted from the input color images acquired by a general-purpose flat-bed scanner) for ACC-based SF quantification, by considering only the covering percentage. To evaluate the accuracy of the proposed fully automatic approach, we compared the SFs obtained by MF2C3 against the conventional counting procedure on four different cell lines. The achieved results revealed a high correlation with the ground-truth measurements based on colony counting, by outperforming our previously validated method using local thresholding on L*u*v* color well images. In conclusion, the proposed multi-feature approach, which inherently leverages the concept of symmetry in the pixel local distributions, might be reliably used in biological studies.


Author(s):  
Guang Dong ◽  
Zheng-Dong Ma ◽  
Gregory Hulbert ◽  
Noboru Kikuchi

The topology optimization method is extended for the optimization of geometrically nonlinear, time-dependent multibody dynamics systems undergoing nonlinear responses. In particular, this paper focuses on sensitivity analysis methods for topology optimization of general multibody dynamics systems, which include large displacements and rotations and dynamic loading. The generalized-α method is employed to solve the multibody dynamics system equations of motion. The developed time integration incorporated sensitivity analysis method is based on a linear approximation of two consecutive time steps, such that the generalized-α method is only applied once in the time integration of the equations of motion. This approach significantly reduces the computational costs associated with sensitivity analysis. To show the effectiveness of the developed procedures, topology optimization of a ground structure embedded in a planar multibody dynamics system under dynamic loading is presented.


Author(s):  
Francisco González ◽  
Pierangelo Masarati ◽  
Javier Cuadrado ◽  
Miguel A. Naya

Formulating the dynamics equations of a mechanical system following a multibody dynamics approach often leads to a set of highly nonlinear differential-algebraic equations (DAEs). While this form of the equations of motion is suitable for a wide range of practical applications, in some cases it is necessary to have access to the linearized system dynamics. This is the case when stability and modal analyses are to be carried out; the definition of plant and system models for certain control algorithms and state estimators also requires a linear expression of the dynamics. A number of methods for the linearization of multibody dynamics can be found in the literature. They differ in both the approach that they follow to handle the equations of motion and the way in which they deliver their results, which in turn are determined by the selection of the generalized coordinates used to describe the mechanical system. This selection is closely related to the way in which the kinematic constraints of the system are treated. Three major approaches can be distinguished and used to categorize most of the linearization methods published so far. In this work, we demonstrate the properties of each approach in the linearization of systems in static equilibrium, illustrating them with the study of two representative examples.


Author(s):  
Martin M. Tong

Numerical solution of the dynamics equations of a flexible multibody system as represented by Hamilton’s canonical equations requires that its generalized velocities q˙ be solved from the generalized momenta p. The relation between them is p = J(q)q˙, where J is the system mass matrix and q is the generalized coordinates. This paper presents the dynamics equations for a generic flexible multibody system as represented by p˙ and gives emphasis to a systematic way of constructing the matrix J for solving q˙. The mass matrix is shown to be separable into four submatrices Jrr, Jrf, Jfr and Jff relating the joint momenta and flexible body mementa to the joint coordinate rates and the flexible body deformation coordinate rates. Explicit formulas are given for these submatrices. The equations of motion presented here lend insight to the structure of the flexible multibody dynamics equations. They are also a versatile alternative to the acceleration-based dynamics equations for modeling mechanical systems.


1998 ◽  
Vol 120 (1) ◽  
pp. 8-14 ◽  
Author(s):  
Marco A. Arteaga

Control design of flexible robot manipulators can take advantage of the structural properties of the model used to describe the robot dynamics. Many of these properties are physical characteristics of mechanical systems whereas others arise from the method employed to model the flexible manipulator. In this paper, the modeling of flexible-link robot manipulators on the basis of the Lagrange’s equations of motion combined with the assumed modes method is briefly discussed. Several notable properties of the dynamic model are presented and their impact on control design is underlined.


Author(s):  
Martin Schulze ◽  
Stefan Dietz ◽  
Bernhard Burgermeister ◽  
Andrey Tuganov ◽  
Holger Lang ◽  
...  

Current challenges in industrial multibody system simulation are often beyond the classical range of application of existing industrial simulation tools. The present paper describes an extension of a recursive order-n multibody system (MBS) formulation to nonlinear models of flexible deformation that are of particular interest in the dynamical simulation of wind turbines. The floating frame of reference representation of flexible bodies is generalized to nonlinear structural models by a straightforward transformation of the equations of motion (EoM). The approach is discussed in detail for the integration of a recently developed discrete Cosserat rod model representing beamlike flexible structures into a general purpose MBS software package. For an efficient static and dynamic simulation, the solvers of the MBS software are adapted to the resulting class of MBS models that are characterized by a large number of degrees of freedom, stiffness, and high frequency components. As a practical example, the run-up of a simplified three-bladed wind turbine is studied where the dynamic deformations of the three blades are calculated by the Cosserat rod model.


Sign in / Sign up

Export Citation Format

Share Document