An Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation
|
|
- Peregrine Logan
- 5 years ago
- Views:
Transcription
1 An Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation Johannes Gerstmayr 7. WORKSHOP ÜBER DESKRIPTORSYSTEME March 2005, Liborianum, Paderborn, Germany Austrian Academy of Sciences, Institute for Technical Mechanics Johannes Kepler University of Linz AUSTRIA Support of the author by a grant of the Austrian Academy of Sciences (APART-scholarship) is gratefully acknowledged.
2 Overview Goal: Flexible multibody simulation Time integration code HOTINT Nonlinear equations in K-form Large scale systems: Structure and improvement of Jacobian Examples: Large number of equations Example: High order of convergence
3 Goal Develop formalisms for the modeling of flexible multibody systems: new flexible bodies and special joints o Multibody system simulation (3D industrial applications) Rigid bodies (non-constant mass matrix) Flexible bodies (non-linear stiffness matrix) Special joints (e.g. sliding joint, including ODEs) Hydraulics, controlled actuators (very stiff) o Needed: Time-integration code Mixed first order / second order ODE / DAE large number of equations (>1000) Integrated multibody and time-integration code
4 Flexible Multibody System Formulations Two approaches for flexible MBS we are interested: Floating Frame of Reference Formulation (FFRF) Absolute Nodal Coordinate Formulation (ANCF) Very stiff equations ( condition numbers as high as possible ) Numerical solution for MBS with small deformations: ANCF: + constant mass matrix + nonlinear stiffness matrix can be linearized + approximate Jacobian, no factorization for all time steps (except for nonl. algebra. equ.) FFRF: - non-constant mass matrix + If coupling of deformations and rotations is week Jacobian can be re-used for larger number of time steps (typically 10s-1000s) + Modal analysis can be used in the framework of FFRF
5 High order Implicit Runge Kutta (IRK) schemes HOTINT High Order Time INTegrator Mixed first order + second order ODE + algebraic equations solver Generally implemented for implicit Runge Kutta schemes Properties: Object-oriented multibody formulation (C++) Most operations done on element level, object oriented Gauss, LobattoIIIA, LobattoIIIC, RadauIA, RadauIIA; up to s=20 Specially adapted to (flexible) multibody system simulation Adapted to redundant coordinates multibody formulation, sparse system matrices Object oriented structure (Classes, Inheritance): Element (Class): Rigid body, flexible body, joint, contact, control, actuator, sensor, Element provides first/second order ODE and algebraic equations, knows what to do Forces: point load, torque, body load,
6 Main structure of HOTINT Equations of motion for mixed 1 st, 2 nd order and algebraic equations, semi-implicit form: (, u&, x, z) u&& = F2 ( u, u&, x, z) F1 ( u, u&, x, z) G( u, u&, x, z) M u x& = 0 = with u...second order ODE variables x... first order ODE variables z...algebraic variable M (non-constant) Mass matrix F 2 includes damping, gyroscopic, elastic and external forces; constraint forces F 1 right-hand-side of first order equations (e.g. hydraulics, control, ) G G G G G algebraic constraints,,,, is available u u& x z
7 Structure of HOTINT (contd.) General Implicit Runge Kutta methods are used to solve the DAE Index 2 formulation mostly used (numerical drift off in many cases neglectable, or stabilized) Runge-Kutta methods written in K-form unknowns: K, K z, nonlinear equations for every stage i: iv ix, i ( Kiu ) ( vi ) ( ) i iv F = 2 ui, vi, xi, zi Kix F1 ( ui, vi, xi, zi ) 0 G( u, v, x, z ) M K i i i i vi = v0 + τ j= 1 AijK jv n, with ui = u0 + τ j= 1 AijK ju, n xi = x0 + τ j= 1 AijK jx n K iu = v i is eliminated The Jacobian is computed for every element w.r.t. ( u, v, x, z ) i i i i The evaluation step does not require the factorization of the mass matrix For first order equations, g-form, adaptive stepsize RadauIIA-methods see E. Hairer and G. Wanner. Stiff differential equations solved by Radau methods or the RADAU5 code. Available via WWW at URL, ftp://ftp.unige.ch/pub/doc/math/stiff/radau5.f (1996).
8 Sketch of the Jacobian for the K-form IRK method
9 Transformation of the Jacobian Split up system into one part with dominating band-structure ( ; mostly terms of mass and stiffness matrix) and remaining part (mostly algebraic equations). The Jacobian is written as J bb Jbb Jbz qb Rb = J J q R zb zz z z Requirement: J bb must be regular for rigid bodies with Euler parameters, the algebraic equations are part of the bb-system and are sparse matrices (Lagrange multipliers; depends on joints) J zb J J bz Size of is much larger than size of bb The solution is computed by means of the transformed equation: ( ) 1 with the matrix Q = J J J J. zz zz zb bb bz J zz 1 1 qz = Qzz ( Rz JzbJbbRb) ( ) q = J R J q 1 b bb b bz z Alternatives: Sparse matrix solver (SuperLU, PARDISO)
10 Advantages (compared to standard ODE/DAE integrators and multistep integrators) Smaller number of unknowns for time-integration compared to implicit-dae approaches (DASSL) or compared to RadauIIA (constant mass matrix) No symbolic inversion of possibly non-linear and non-constant mass matrix High order of convergence; orders up to 20 (10-stage-Gauss) possible for real-life multibody systems No restarting is necessary after a discontinuous step Factorization and back-substitution is not expensive compared to evaluation (~1K-10K nonlinear equations). RadauIIA (s>1) methods are useful for the development/verification of new joints (Index 3) Numerical damping is low for most methods (Gauss, LobattoIIIA; RadauIIA, s>1 ok) Disadvantages Increasing size of equations for nonlinear system, larger Jacobian than in multistep methods Condition numbers become worse with larger number of stages Error estimator based on half-steps (similar to idea of Richardson-extrapolation)
11 Example: Sliding 3D Beam with Eccentricity Beam modeled by absolute nodal coordinate formulation, example model from Sugiyama et al. sliding joint flexible cable sliding beam
12 Computational times for dense and transformed sparse factorisation Example: Sliding flexible beam along flexible cable (Sugiyama et al. 2003) 0.2 seconds simulation, 2000 steps, trapezoidal rule sparse + transformed dense factorization linear quadratic computational time unknowns
13 Main parts of the computational effort in the multibody simulation Computational effort versus unknowns: 100% 90% 80% 70% 60% 50% 40% 30% Rest Factorize Apply Evaluation 20% 10% 0% Evaluation: Computation of mass matrix, elastic and external forces, constraints, Apply: Explicit computations (multiplications and back-substitution)
14 Example: 3D pantograph/catenary system Integrated simulation: Pantograph (spring-damper system with, rigid bodies) and Catenary (flexible, large deformation, ANCF) as one multibody system 279 flexible bodies and 4 rigid bodies, ca second order ODEs, ~1 hour for 8seconds
15 Contact force at pan-head for 32 m/s (left) and 45 m/s (right).
16 2-Arm crane with slender pipe, contact and friction (with M. Stangl) Actuators (represented by moments M1 and M2) placed in first and second hinge dry friction in third hinge Degrees of freedom: 4 Rigid body DOF / 4 elastic DOF Discontinuity: contact with surface and dry friction in hinge 3 nonlinear control of actuators in hinges 1 and 2 with limiter Control: PD-Control 3y MR3(t) 4y -q3(t) 3x 4x w4(x,t) q4(t) w2(x,t) w3(x,t) 2y q2(t) 2x 1y Iy w1(x,t) M2 M1 q1(t) 1x Ix h FK
17 Animation: Control of Crane-mechanism with contact and friction moment
18 Optimal stages for computational time? Test case: Absolute accuracy 1e-4, adaptive timestep, simulate 3 seconds discontinuous system, 2-arm crane with pipe, contact and friction Computational time [s] stages Computational time, Radau IIA, stages 2 to 8, ODE-Order 3 to 15
19 Number of time steps and number of computed Jacobians for different number of stages Radau IIA, absolute accuracy =10-4 Number of timesteps 7 x stages Number of Jacobians stages
20 ODE - Accuracy Convergence of Continuous system (ODE-variables) (no friction/contact, const. control parameters) RadauIIA s=3 RadauIIA s=6 Gauss s=3 Gauss s=6 investigated: RadauIIA RadauIA Gauss LobattoIIIA/C RadauIIA and LobattoIIIA perform best Timestepsize [s] Good agreement with theory, see Hairer and Wanner, Solving ordinary differential equations I (II), Springer Verlag Berlin Heidelberg, 1991
21 Convergence of Continuous system (no friction/contact, const. Control) Computational time versus accuracy with 2 arm-crane with pipe, adaptive timestep stages 3 stages 4 stages 5 stages 6 stages 7 stages 8 stages Computation Time [s] ODE-Accuracy
22 Conclusion Implicit Runge Kutta methods show good performance for large flexible multibody systems Nearly order-n formalism due to numerical methods Good convergence properties for small multibody systems and higher order Integrated time integration and object oriented multibody code is favorable Outlook / Question: Automatic stepsize + order + method? Error estimation for stiff equations?
Numerical Integration of Equations of Motion
GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.
More informationA MULTI-BODY ALGORITHM FOR WAVE ENERGY CONVERTERS EMPLOYING NONLINEAR JOINT REPRESENTATION
Proceedings of the ASME 2014 33rd International Conference on Ocean, Offshore and Arctic Engineering OMAE2014 June 8-13, 2014, San Francisco, California, USA OMAE2014-23864 A MULTI-BODY ALGORITHM FOR WAVE
More informationINTEGRATION OF THE EQUATIONS OF MOTION OF MULTIBODY SYSTEMS USING ABSOLUTE NODAL COORDINATE FORMULATION
INTEGRATION OF THE EQUATIONS OF MOTION OF MULTIBODY SYSTEMS USING ABSOLUTE NODAL COORDINATE FORMULATION Grzegorz ORZECHOWSKI *, Janusz FRĄCZEK * * The Institute of Aeronautics and Applied Mechanics, The
More informationA 3D finite element method for flexible multibody systems
Multibody Syst Dyn (2006) 15:309 324 DOI 10.1007/s11044-006-9009-3 A 3D finite element method for flexible multibody systems Johannes Gerstmayr Joachim Schöberl Received: 22 July 2004 / Accepted: 4 March
More informationNumerical integration of DAE s
Numerical integration of DAE s seminar Sandra Allaart-Bruin sbruin@win.tue.nl seminar p.1 Seminar overview February 18 Arie Verhoeven Introduction to DAE s seminar p.2 Seminar overview February 18 Arie
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationLecture «Robot Dynamics»: Dynamics and Control
Lecture «Robot Dynamics»: Dynamics and Control 151-0851-00 V lecture: CAB G11 Tuesday 10:15 12:00, every week exercise: HG E1.2 Wednesday 8:15 10:00, according to schedule (about every 2nd week) Marco
More informationAN ALGORITHM FOR TOPOLOGY OPTIMIZATION
AN ALGORITHM FOR TOPOLOGY OPTIMIZATION OF MULTIBODY SYSTEMS TOHEED GHANDRIZ Master s thesis 2014:E52 Faculty of Science Centre for Mathematical Sciences Numerical Analysis CENTRUM SCIENTIARUM MATHEMATICARUM
More informationInvestigation on the Most Efficient Ways to Solve the Implicit Equations for Gauss Methods in the Constant Stepsize Setting
Applied Mathematical Sciences, Vol. 12, 2018, no. 2, 93-103 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.711340 Investigation on the Most Efficient Ways to Solve the Implicit Equations
More informationTime integration. DVI and HHT time stepping methods in Chrono
Time integration DVI and HHT time stepping methods in Chrono Time Integration in Chrono Two classes of time stepping methods in Chrono Time steppers for smooth dynamics Classical multibody dynamics rigid
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More informationOn the Diagonal Approximation of Full Matrices
On the Diagonal Approximation of Full Matrices Walter M. Lioen CWI P.O. Box 94079, 090 GB Amsterdam, The Netherlands ABSTRACT In this paper the construction of diagonal matrices, in some
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More information19.2 Mathematical description of the problem. = f(p; _p; q; _q) G(p; q) T ; (II.19.1) g(p; q) + r(t) _p _q. f(p; v. a p ; q; v q ) + G(p; q) T ; a q
II-9-9 Slider rank 9. General Information This problem was contributed by Bernd Simeon, March 998. The slider crank shows some typical properties of simulation problems in exible multibody systems, i.e.,
More informationSoftware Verification
EXAMPLE 6-6 LINK SUNY BUFFALO DAMPER WITH LINEAR VELOCITY EXPONENT PROBLEM DESCRIPTION This example comes from Section 5 of Scheller and Constantinou 1999 ( the SUNY Buffalo report ). It is a two-dimensional,
More informationNumerical Methods for Rigid Multibody Dynamics
Numerical Methods for Rigid Multibody Dynamics Claus Führer Centre for Mathematical Sciences Lund University Lappenranta 2012 Unit 0: Preface These notes serve as a skeleton for the compact course. They
More informationA New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations
A New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations Rei-Wei Song and Ming-Gong Lee* d09440@chu.edu.tw, mglee@chu.edu.tw * Department of Applied Mathematics/
More informationApproach based on Cartesian coordinates
GraSMech course 2005-2006 Computer-aided analysis of rigid and flexible multibody systems Approach based on Cartesian coordinates Prof. O. Verlinden Faculté polytechnique de Mons Olivier.Verlinden@fpms.ac.be
More informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More informationModelling Physical Phenomena
Modelling Physical Phenomena Limitations and Challenges of the Differential Algebraic Equations Approach Olaf Trygve Berglihn Department of Chemical Engineering 30. June 2010 2 Outline Background Classification
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationFunctional Mockup Interface (FMI)
Functional Mockup Interface (FMI) A framework for coarse-grained parallel time integration in nonlinear system dynamics PinT 2015 4th Workshop on Parallel-in-Time Integration May 2015, Dresden, Germany
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationq 1 F m d p q 2 Figure 1: An automated crane with the relevant kinematic and dynamic definitions.
Robotics II March 7, 018 Exercise 1 An automated crane can be seen as a mechanical system with two degrees of freedom that moves along a horizontal rail subject to the actuation force F, and that transports
More informationMethods of Analysis. Force or Flexibility Method
INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses
More informationCable-Pulley Interaction with Dynamic Wrap Angle Using the Absolute Nodal Coordinate Formulation
Proceedings of the 4 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'17) Toronto, Canada August 21 23, 2017 Paper No. 133 DOI: 10.11159/cdsr17.133 Cable-Pulley Interaction with
More informationSYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationTechnische Universität Berlin
Technische Universität Berlin Institut für Mathematik M7 - A Skateboard(v1.) Andreas Steinbrecher Preprint 216/8 Preprint-Reihe des Instituts für Mathematik Technische Universität Berlin http://www.math.tu-berlin.de/preprints
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationME751 Advanced Computational Multibody Dynamics
ME751 Advanced Computational Multibody Dynamics Inverse Dynamics Equilibrium Analysis Various Odd Ends March 18, 2010 Dan Negrut, 2010 ME751, UW-Madison Action speaks louder than words but not nearly as
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationTime-adaptive methods for the incompressible Navier-Stokes equations
Time-adaptive methods for the incompressible Navier-Stokes equations Joachim Rang, Thorsten Grahs, Justin Wiegmann, 29.09.2016 Contents Introduction Diagonally implicit Runge Kutta methods Results with
More informationDynamics and control of mechanical systems
Dynamics and control of mechanical systems Date Day 1 (03/05) - 05/05 Day 2 (07/05) Day 3 (09/05) Day 4 (11/05) Day 5 (14/05) Day 6 (16/05) Content Review of the basics of mechanics. Kinematics of rigid
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationLecture «Robot Dynamics»: Dynamics 2
Lecture «Robot Dynamics»: Dynamics 2 151-0851-00 V lecture: CAB G11 Tuesday 10:15 12:00, every week exercise: HG E1.2 Wednesday 8:15 10:00, according to schedule (about every 2nd week) office hour: LEE
More informationSolving Orthogonal Matrix Differential Systems in Mathematica
Solving Orthogonal Matrix Differential Systems in Mathematica Mark Sofroniou 1 and Giulia Spaletta 2 1 Wolfram Research, Champaign, Illinois, USA. marks@wolfram.com 2 Mathematics Department, Bologna University,
More informationNumerical Methods for Engineers
Numerical Methods for Engineers SEVENTH EDITION Steven C Chopra Berger Chair in Computing and Engineering Tufts University Raymond P. Canal Professor Emeritus of Civil Engineering of Michiaan University
More informationExample 37 - Analytical Beam
Example 37 - Analytical Beam Summary This example deals with the use of RADIOSS linear and nonlinear solvers. A beam submitted to a concentrated load on one extremity and fixed on the other hand is studied.
More informationControl of Earthquake Induced Vibrations in Asymmetric Buildings Using Passive Damping
Control of Earthquake Induced Vibrations in Asymmetric Buildings Using Passive Damping Rakesh K. Goel, California Polytechnic State University, San Luis Obispo Abstract This paper summarizes the results
More informationIMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 8 (1), 17-30 (2002) IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM MAŁGORZATA JANKOWSKA 1, ANDRZEJ MARCINIAK 1,2 1 Poznań University
More informationA Rosenbrock Nystrom state space implicit approach for the dynamic analysis of mechanical systems: II method and numerical examples
73 A Rosenbrock Nystrom state space implicit approach for the dynamic analysis of mechanical systems: II method and numerical examples D Negrut *, A Sandu, E J Haug 3, F A Potra 4 and C Sandu 5 MSCsoftware,
More informationLesson Rigid Body Dynamics
Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationMulti Linear Elastic and Plastic Link in SAP2000
26/01/2016 Marco Donà Multi Linear Elastic and Plastic Link in SAP2000 1 General principles Link object connects two joints, i and j, separated by length L, such that specialized structural behaviour may
More informationContents. Dynamics and control of mechanical systems. Focus on
Dynamics and control of mechanical systems Date Day 1 (01/08) Day 2 (03/08) Day 3 (05/08) Day 4 (07/08) Day 5 (09/08) Day 6 (11/08) Content Review of the basics of mechanics. Kinematics of rigid bodies
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationBoundary Nonlinear Dynamic Analysis
Boundary Nonlinear Dynamic Analysis Damper type Nonlinear Link Base Isolator type Nonlinear Link Modal Nonlinear Analysis : Equivalent Dynamic Load Damper type Nonlinear Link Visco-Elastic Damper (VED)
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More informationDynamics. describe the relationship between the joint actuator torques and the motion of the structure important role for
Dynamics describe the relationship between the joint actuator torques and the motion of the structure important role for simulation of motion (test control strategies) analysis of manipulator structures
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More informationStable Adaptive Co-simulation: A Switched Systems Approach. Cláudio Gomes, Benoît Legat, Raphaël M. Jungers, Hans Vangheluwe
Stable Adaptive Co-simulation: A Switched Systems Approach Cláudio Gomes, Benoît Legat, Raphaël M. Jungers, Hans Vangheluwe Agenda 1. Background 2. Adaptive Orchestration 3. Contribution 4. Conclusion
More informationEuropean Consortium for Mathematics in Industry E. Eich-Soellner and C. FUhrer Numerical Methods in Multibody Dynamics
European Consortium for Mathematics in Industry E. Eich-Soellner and C. FUhrer Numerical Methods in Multibody Dynamics European Consortium for Mathematics in Industry Edited by Leif Arkeryd, Goteborg Heinz
More informationA multibody dynamics model of bacterial
A multibody dynamics model of bacterial Martin Servin1 1 UMT Research Lab - Department of Physics Umea University August 26, 2015 (1 : 17) Nonsmooth multidomain dynamics What has this to do with bacterias?
More information1.053J/2.003J Dynamics and Control I Fall Final Exam 18 th December, 2007
1.053J/2.003J Dynamics and Control I Fall 2007 Final Exam 18 th December, 2007 Important Notes: 1. You are allowed to use three letter-size sheets (two-sides each) of notes. 2. There are five (5) problems
More informationNumerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.
Numerical Methods for Engineers An Introduction with and Scientists Applications using MATLAB Third Edition Amos Gilat Vish- Subramaniam Department of Mechanical Engineering The Ohio State University Wiley
More informationA complete strategy for efficient and accurate multibody dynamics of flexible structures with large lap joints considering contact and friction
Multibody Syst Dyn (2017) 40:407 436 DOI 10.1007/s11044-016-9555-2 A complete strategy for efficient and accurate multibody dynamics of flexible structures with large lap joints considering contact and
More informationVibration Dynamics and Control
Giancarlo Genta Vibration Dynamics and Control Spri ringer Contents Series Preface Preface Symbols vii ix xxi Introduction 1 I Dynamics of Linear, Time Invariant, Systems 23 1 Conservative Discrete Vibrating
More informationInvestigation of falling control rods in deformed guiding tubes in nuclear reactors using multibody approaches
The th Joint International Conference on Multibody System Dynamics June 8, 8, Lisboa, Portugal Investigation of falling control rods in deformed guiding tubes in nuclear reactors using multibody approaches
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationThe Absolute Nodal Coordinate Formulation
The Absolute Nodal Coordinate Formulation ANCF Antonio Recuero Dan Negrut May 27, 2016 Abstract This white paper describes the fundamentals of the nonlinear nite element theory used to implement ANCF nite
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationInverse Dynamics of Flexible Manipulators
Technical report from Automatic Control at Linköpings universitet Inverse Dynamics of Flexible Manipulators Stig Moberg, Sven Hanssen Division of Automatic Control E-mail: stig@isy.liu.se, soha@kth.se
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More informationME8230 Nonlinear Dynamics
ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass
More informationDynamics of Multibody Systems: Conventional and Graph-Theoretic Approaches
Dynamics of Multibody Systems: Conventional and Graph-Theoretic Approaches by Professor John J. McPhee, P.Eng. Systems Design Engineering (Cross-appointed to Mechanical Engineering) University of Waterloo,
More informationTheory & Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati
Theory & Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 5 Torsional Vibrations Lecture - 4 Transfer Matrix Approach
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationModule 4: Numerical Methods for ODE. Michael Bader. Winter 2007/2008
Outlines Module 4: for ODE Part I: Basic Part II: Advanced Lehrstuhl Informatik V Winter 2007/2008 Part I: Basic 1 Direction Fields 2 Euler s Method Outlines Part I: Basic Part II: Advanced 3 Discretized
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. I - Numerical Analysis and Methods for Ordinary Differential Equations - N.N. Kalitkin, S.S.
NUMERICAL ANALYSIS AND METHODS FOR ORDINARY DIFFERENTIAL EQUATIONS N.N. Kalitkin Institute for Mathematical Modeling, Russian Academy of Sciences, Moscow, Russia S.S. Filippov Keldysh Institute of Applied
More informationRobotics. Dynamics. Marc Toussaint U Stuttgart
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Rigid body physics Particle system Most simple instance of a physics system Each object (body) is a particle Each particle
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationNumerical Methods for Differential Equations
CHAPTER 5 Numerical Methods for Differential Equations In this chapter we will discuss a few of the many numerical methods which can be used to solve initial value problems and one-dimensional boundary
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationDynamic Simulation of the EMS Maglev Vehicle-Guideway-Controller Coupling System
Dynamic Simulation of the EMS Maglev Vehicle-Guideway-Controller Coupling System Wanming ZHAI*, Chunfa ZHAO**, Chengbiao CAI** * National Power Traction Lab, Southwest Jiaotong University, Chengdu 610031,
More informationIntroduction to SAMCEF MECANO
Introduction to SAMCEF MECANO 1 Outline Introduction Generalized coordinates Kinematic constraints Time integration Description and paramerization of finite rotations Acknowledgements Michel Géradin (ULg)
More informationRicharson Extrapolation for Runge-Kutta Methods
Richarson Extrapolation for Runge-Kutta Methods Zahari Zlatevᵃ, Ivan Dimovᵇ and Krassimir Georgievᵇ ᵃ Department of Environmental Science, Aarhus University, Frederiksborgvej 399, P. O. 358, 4000 Roskilde,
More informationRobotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007
Robotics & Automation Lecture 25 Dynamics of Constrained Systems, Dynamic Control John T. Wen April 26, 2007 Last Time Order N Forward Dynamics (3-sweep algorithm) Factorization perspective: causal-anticausal
More informationEffect of Angular movement of Lifting Arm on Natural Frequency of Container Lifting Mechanism using Finite Element Modal Analysis
Effect of Angular movement of Lifting Arm on Natural Frequency of Container Lifting Mechanism using Finite Element Modal Analysis Khodu M Dhameliya, 2 Utpal V Shah, 3 Dhaval Makwana, 4 Mansi Yadav, 5 Ishankumar
More informationBindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 18. HW 7 is posted, and will be due in class on 4/25.
Logistics Week 12: Monday, Apr 18 HW 6 is due at 11:59 tonight. HW 7 is posted, and will be due in class on 4/25. The prelim is graded. An analysis and rubric are on CMS. Problem du jour For implicit methods
More informationDynamic Model of a Badminton Stroke
ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: +45 994 9317 / Fax: +45 9815
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationAdvanced methods for ODEs and DAEs
Advanced methods for ODEs and DAEs Lecture 13: Overview Bojana Rosic, 13. Juli 2016 Dynamical system is a system that evolves over time possibly under external excitations: river flow, car, etc. The dynamics
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationModule SimBeam with Element of Linear Variable Cross Sections. SIMPACK USER Meeting 2011
Module SimBeam with Element of Linear Variable Cross Sections by Prof. Dr. Oskar Wallrapp, Wessling, Germany Dr. Stefan Dietz, SIMPACK AG Byunghoon Song (Trainee Student) SIMPACK AG SIMPACK USER Meeting
More informationCable installation simulation by using a multibody dynamic model
Multibody Syst Dyn (2013) 30:433 447 DOI 10.1007/s11044-013-9364-9 Cable installation simulation by using a multibody dynamic model Cai Jin Yang Di Feng Hong Ge Xue Ren Zhi Hua Zhao Received: 20 February
More informationQuadratic SDIRK pair for treating chemical reaction problems.
Quadratic SDIRK pair for treating chemical reaction problems. Ch. Tsitouras TEI of Chalkis, Dept. of Applied Sciences, GR 34400 Psahna, Greece. I. Th. Famelis TEI of Athens, Dept. of Mathematics, GR 12210
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationControl of constrained spatial three-link flexible manipulators
Control of constrained spatial three-link flexible manipulators Sinan Kilicaslan, M. Kemal Ozgoren and S. Kemal Ider Gazi University/Mechanical Engineering Department, Ankara, Turkey Middle East Technical
More informationFinal Exam Solution Dynamics :45 12:15. Problem 1 Bateau
Final Exam Solution Dynamics 2 191157140 31-01-2013 8:45 12:15 Problem 1 Bateau Bateau is a trapeze act by Cirque du Soleil in which artists perform aerial maneuvers on a boat shaped structure. The boat
More informationChapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs
Chapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs In the previous chapter, we investigated stiffness in ODEs. Recall that an ODE is stiff if it exhibits behavior on widelyvarying timescales.
More information