A General Purpose Formulation for Nonsmooth Dynamics With Finite Rotations: Application to the Woodpecker Toy

2020 ◽  
Vol 16 (3) ◽  
Author(s):  
Alejandro Cosimo ◽  
Federico J. Cavalieri ◽  
Javier Galvez ◽  
Alberto Cardona ◽  
Olivier Brüls
Keyword(s):  
Finite Element ◽  
Multibody Dynamics ◽  
Equations Of Motion ◽  
Horizontal Displacement ◽  
General Purpose ◽  
Nonsmooth Dynamics ◽  
Unilateral Constraints ◽  
Simulation Code ◽  
Frictional Contacts ◽  
Large Rotations

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.

A General Purpose Formulation for Nonsmooth Dynamics Including Large Rotations: Application to the Woodpecker Toy

2019 ◽  
Author(s):  
Javier Galvez ◽  
Alejandro Cosimo ◽  
Federico J. Cavalieri ◽  
Alberto Cardona ◽  
Olivier Brüls
Keyword(s):  
Finite Element ◽  
Equations Of Motion ◽  
General Purpose ◽  
Rigid Bodies ◽  
Nonsmooth Dynamics ◽  
Minimal Set ◽  
Frictional Contacts ◽  
Bilateral Constraints ◽  
Large Rotations ◽  
Element Approach

Abstract The aim of this work is to extend the finite element multibody dynamics approach to problems involving frictional contacts and impacts. Since rigid bodies and joints involve bilateral constraints, it is important to avoid any drift phenomenon. Therefore, the nonsmooth generalized-α method is used, which imposes the constraints both at position and at velocity levels. Its low intrusiveness allows one to reuse an existing library of elements without major modifications. The study of the woodpecker toy dynamics sets up a good example to show the capabilities of the nonsmooth generalized-α within the context of a general finite element framework. This example has already been studied by many authors who generally adopt a model with a minimal set of coordinates and small rotations. We show that using a finite element approach, the equations of motion can be assembled automatically, and large rotations can be easily considered.

Adding Non-Smooth Analysis Capabilities to General-Purpose Multibody Dynamics by Co-Simulation

2013 ◽  
Author(s):  
Matteo Fancello ◽  
Pierangelo Masarati ◽  
Marco Morandini
Keyword(s):  
Differential Algebraic Equations ◽  
General Purpose ◽  
Rigid Body Dynamics ◽  
Nonsmooth Dynamics ◽  
Algebraic Equations ◽  
Unilateral Constraints ◽  
Frictional Contacts ◽  
Continuous Contact ◽  
Hard Constraints ◽  
Dynamics Problems

Multi-rigid-body dynamics problems with unilateral constraints, like frictionless and frictional contacts, are characterized by nonsmooth dynamics. The issue of nonsmoothness can be addressed with methods that apply a mathematical regularization, called continuous contact methods; alternatively, hard constraints with complementarity approaches can be proficiently used. This work presents an attempt at integrating consistently modeled unilateral constraints in a general purpose multibody formulation and implementation originally designed to address intrinsically smooth problems. The focus is on the analysis of generally smooth problems, characterized by significant multidisciplinarity, with the need to selectively include nonsmooth events localized in time and in specific components of the model. A co-simulation approach between the smooth Differential-Algebraic Equations solver and the classic Moreau-Jean timestepping approach is devised as an alternative to entirely redesigning a monolithic nonsmooth solver, in order to provide elements subject to frictionless and frictional contact in the general-purpose, free multibody solver MBDyn. The implementation uses components from the INRIA’s Siconos library for the solution of Complementarity Problems. The proposed approach is applied to several problems of increasing complexity to empirically evaluate its properties and versatility. The applicability of the family of second-order accurate, A/L stable multistep integration algorithms used by MBDyn to nonsmooth dynamics is also discussed and assessed.

Kinematics for unilateral constraints in multibody dynamics

2016 ◽  
Vol 22 (8) ◽  
pp. 1654-1687
Author(s):  
P Lidström
Keyword(s):  
Multibody Dynamics ◽  
Distance Function ◽  
Contact Surface ◽  
Equations Of Motion ◽  
Time Derivative ◽  
Unilateral Constraint ◽  
Contact Forces ◽  
Unilateral Constraints ◽  
Constraint Forces ◽  
Friction Forces

This paper is concerned with the kinematics of unilateral constraints in multibody dynamics. These constraints are related to the contact between parts and the principle of impenetrability of matter and have the property that they may be active, in which case they give rise to constraint forces, or passive, in which case they do not give rise to constraint forces. In order to check whether the constraint is active or passive a distance function between parts of the multibody is required. The paper gives a rigorous definition of the distance function and derives certain of its properties. The unilateral constraint may then be expressed in terms of this distance function. The paper analyses the transitions from passive constraints to active and vice versa. Sufficient regularity of the transplacements of the parts and their boundary surfaces will lead to specific properties of the time derivative of the distance function. When the unilateral constraint is active then the parts are geometrically in contact and there is a certain contact surface that, in specific cases, may degenerate into a point. If the parts are in mechanical contact over the contact surface then there will be an interaction between the parts given by contact forces, such as normal and friction forces. Parts in contact may be at rest relative to one another, over the contact surface, or they may be in relative sliding motion. The transition from non-sliding contact to sliding and from sliding to non-sliding is discussed and necessary conditions on the relative velocity and the traction vector are derived. Appropriate complementary conditions are then formulated. These are instrumental when the technique of linear complementarity is used in order to find solutions to the equations of motion.

Finite Element Method of Flexible Multibody Dynamics Based on a New Kind of Floating Frame

1991 ◽  
Author(s):  
You-Fang Lu ◽  
Zhao-Hui Qi ◽  
Bin Wang ◽  
Guan-Min Feng
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Numerical Analysis ◽  
Multibody Dynamics ◽  
Equations Of Motion ◽  
Flexible Multibody Dynamics ◽  
Flexible Multibody ◽  
Floating Frame ◽  
Element Method

Abstract A new kind of floating frame whose parameters do not appear in equations of motion as additional unknowns is defined. Numerical analysis of flexible multibody dynamics is much facilitated by using finite-element iteration of the corresponding equations based on this concept.

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
Keyword(s):  
Multibody Dynamics ◽  
Equations Of Motion ◽  
Control Design ◽  
Numerical Differentiation ◽  
Difficult Problem ◽  
General Purpose ◽  
Computational Method ◽  
Computational Error ◽  
Dynamics Modeling ◽  
Analytical Approximations

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.

Modeling Rotorcraft Dynamics With Finite Element Multibody Procedures

2009 ◽  
Author(s):  
Olivier A. Bauchau
Keyword(s):  
Finite Element ◽  
Multibody Dynamics ◽  
Equations Of Motion ◽  
Maximum Efficiency ◽  
Flexible Bodies ◽  
Modal Reduction ◽  
Nonlinear Finite Element Methods ◽  
Deformable Bodies ◽  
Joint Element ◽  
Key Aspects

This paper describes a multibody dynamics approach to the modeling of rotorcraft systems and reviews the key aspects of the simulation procedure. The multibody dynamics analysis is cast within the framework of nonlinear finite element methods, and the element library includes rigid and deformable bodies as well as joint elements. No modal reduction is performed for the modeling of flexible bodies. The structural and joint element library is briefly described. The algorithms used to integrate the resulting equations of motion with maximum efficiency and robustness are discussed. Various solution procedures, static, dynamic, stability, and trim analysis, are presented. Post-processing and visualization issues are also addressed. Finally, the paper concludes with selected rotorcraft applications.

scholarly journals A Symbolic Approach to the Multibody Modeling of Road Vehicles

2017 ◽  
Vol 09 (05) ◽  
pp. 1750068 ◽  
Author(s):  
Roberto Lot ◽  
Matteo Massaro
Keyword(s):  
Equations Of Motion ◽  
Automatic Generation ◽  
General Purpose ◽  
Rigid Bodies ◽  
Simulation Code ◽  
Symbolic Form ◽  
Multibody Modeling ◽  
Mixed Coordinates ◽  
Computer Algebra Software ◽  
Symbolic Approach

This paper introduces MBSymba, an object-oriented language for the modeling of multibody systems and the automatic generation of equations of motion in symbolic form. MBSymba has built upon the general-purpose computer algebra software Maple and it is freely available for teaching and research purposes. With MBSymba, objects such as points, vectors, rigid bodies, forces and torques, and the relationships among them may be defined and manipulated both at high and low levels. Absolute, relative or mixed coordinates may be used, as well as combination of infinitesimal and noninfinitesimal variables. Once the system has been modeled, Lagrange’s and/or Newton’s equations can be derived in a quasi-automatic way, either in an inertial or noninertial reference frame. Equations can be automatically converted into Matlab, C/C++ or Fortan code to produce stand alone, numerically optimized simulation code. MBSymba is particularly suited for the modeling of ground, water or air vehicles; therefore, the mathematical model of a passenger car with trailer is illustrated as a case study. Time domain simulations, steady state analysis and stability results are also presented.

HOTINT: A Script Language Based Framework for the Simulation of Multibody Dynamics Systems

2013 ◽  
Author(s):  
Johannes Gerstmayr ◽  
Alexander Dorninger ◽  
Rafael Eder ◽  
Peter Gruber ◽  
Daniel Reischl ◽  
...  
Keyword(s):  
Finite Element ◽  
Multibody Dynamics ◽  
Industrial Applications ◽  
Rigid Bodies ◽  
Interface Element ◽  
Script Language ◽  
Simulation Code ◽  
Research Education ◽  
Element Simulation ◽  
Graphical Visualization

The multibody dynamics and finite element simulation code has been developed since 1997. In the past years, more than 10 researchers have contributed to certain parts of HOTINT, such as solver, graphical user interface, element library, joint library, finite element functionality and port blocks. Currently, a script-language based version of HOTINT is freely available for download, intended for research, education and industrial applications. The main features of the current available version include objects like point mass, rigid bodies, complex point-based joints, classical mechanical joints, flexible (nonlinear) beams, port-blocks for mechatronics applications and many other features such as loads, sensors and graphical objects. HOTINT includes a 3D graphical visualization showing the results immediately during simulation, which helps to reduce modelling errors. In the present paper, we show the current state and the structure of the code. Examples should demonstrate the easiness of use of HOTINT.

The Influence of Loading Position in A Priori High Stress Detection using Mode Superposition

2020 ◽  
Vol 1 (1) ◽  
pp. 93-102
Author(s):  
Carsten Strzalka ◽  
◽  
Manfred Zehn ◽  
Keyword(s):  
Finite Element ◽  
Degrees Of Freedom ◽  
Equations Of Motion ◽  
A Priori ◽  
High Stress ◽  
Von Mises Stress ◽  
Detailed Knowledge ◽  
Variable Loading ◽  
Numerical Effort ◽  
Von Mises

For the analysis of structural components, the finite element method (FEM) has become the most widely applied tool for numerical stress- and subsequent durability analyses. In industrial application advanced FE-models result in high numbers of degrees of freedom, making dynamic analyses time-consuming and expensive. As detailed finite element models are necessary for accurate stress results, the resulting data and connected numerical effort from dynamic stress analysis can be high. For the reduction of that effort, sophisticated methods have been developed to limit numerical calculations and processing of data to only small fractions of the global model. Therefore, detailed knowledge of the position of a component’s highly stressed areas is of great advantage for any present or subsequent analysis steps. In this paper an efficient method for the a priori detection of highly stressed areas of force-excited components is presented, based on modal stress superposition. As the component’s dynamic response and corresponding stress is always a function of its excitation, special attention is paid to the influence of the loading position. Based on the frequency domain solution of the modally decoupled equations of motion, a coefficient for a priori weighted superposition of modal von Mises stress fields is developed and validated on a simply supported cantilever beam structure with variable loading positions. The proposed approach is then applied to a simplified industrial model of a twist beam rear axle.

Sign in / Sign up

Export Citation Format

Share Document