NUMERICAL INTEGRATION OF CONSTRAINED MULTI-BODY DYNAMICAL SYSTEMS USING 5 T H ORDER EXACT ANALYTIC CONTINUATION ALGORITHM
|
|
- Randolf Kelley Jackson
- 5 years ago
- Views:
Transcription
1 (Preprint) AAS NUMERICAL INTEGRATION OF CONSTRAINED MULTI-BODY DYNAMICAL SYSTEMS USING 5 T H ORDER EXACT ANALYTIC CONTINUATION ALGORITHM Ahmad Bani Younes, and James Turner Many numerical integration methods have been developed for predicting the evolution of the response of dynamical systems. Standard algorithms approach approximate the solution at a future time by introducing using a truncated power series representation that attempts to recover an n-th order Taylor series approximation, while only numerically sampling a single derivative model. An exact fifthorder analytic continuation method is presented for integrating constrained multibody vector-valued systems of equations, where the Jacobi form of the Routh-Voss equations of motion simultaneously generates the acceleration and Lagrange multiplier solution. The constraint drift problem is addressed by introducing an analytic continuation method that rigorously enforces the kinematic constraints through five time derivatives. The proposed approach is expected to be particularly useful for stiff dynamical systems, as well as systems where implicit integration formulations are introduced. Numerical examples are presented that demonstrate the effectiveness of the proposed methodology. INTRODUCTION Analytic Continuation Method An exact fifth-order analytic continuation method is presented for integrating vector-valued constrained multi-body dynamical systems of equations. The major innovation of this paper is the introduction of computational differentiation (CD) techniques which enable the development of a rigorous n-th order Taylor series analytic continuation model for vector systems of differential equations. Classical numerical methods sample a single ordinary differential equation motion several times to develop a weighted approximation for the integrated solution. The proposed analytic continuation method computes exact derivatives for all of the terms retained in the series approximation. Two classes of dynamical systems are considered: (1) unconstrained differential equations, and (2) constrained multi-body dynamical systems. Special attention is given to handling matrix-vector equations were the matrix is time-varying and is too large to be symbolically inverted. Both constant and time-varying kinematic constraints are considered for multi-body systems. Numerical results are presented that demonstrate the proposed analytical/numerical modeling approach. For general applications in science and engineering, the system equations of motion are described PhD Candidate, Aerospace Engineering, Texas A&M University, College Station, Texas, olalahmad@gmail.com. Research Professor, Aerospace Engineering, Texas A&M University, College Station, Texas, turner@aero.tamu. edu. 1
2 by the following first-order vector differential equation: ẋ = f(t, x); given x(t 0 ) = a (1) where x R n denotes the state vector, t denotes time, and ẋ = dx dt. A rigorous Taylor series expansion is developed for Eq.(1), that defines the future response of the system by summing contributions form n-time derivatives of f( ), where time is evaluated at the point t = a. Analytically continuing x(t) leads to the following Taylor series expansion x(t) = x(a) + x (a) 1! (t a) + x (a) 2! (t a) 2 + x (a) (t a) 3 + (2) 3! Theoretically, the solution is complete once a numerical method is specified for evaluating the time derivatives for: dx dt, d2 x, d3 x, unfortunately these time derivatives are generally not available because of the analyst time required for deriving and coding and verifying the derivative calculations. dt 2 dt 3 Classical methods overcome this problem by a repeated sampling and averaging strategy. The main contribution of this paper is to demonstrate that CD tools can be brought to bear on the problem of numerically integrating constrained-complex nonlinear vector systems of ordinary differential equations by returning to the fundamental Taylor expansion that defines the evolution of the system response. To this end, CD is used to numerically evaluate derivative approximations for: dx dt, d2 x dt 2, d3 x dt 3, which are then used in the truncated series x(t+h) x(t)+ 5 1 d i x i! h i. This ex- dt i trapolation equation replaces the weighted averaging approach of conventional solution algorithms. Numerical time derivative values are computed from the following cascade of tensor contraction operations: ẋ = f ẍ = f.ẋ x (3) = 2 f.ẋ.ẋ + f.ẍ x (4) = 3 f.ẋ.ẋ.ẋ f.ẍ.ẋ + 2 f.ẋ.ẍ + f.x (3) (3) x (5) = 4 f.ẋ.ẋ.ẋ.ẋ f.ẍ.ẋ.ẋ f.ẋ.ẍ.ẋ f.x (3).ẋ + 3 f.ẋ.ẋ.ẍ f.ẍ.ẍ + 2 f.ẋ.x (3) + f.x (4) These equations assume that time does not appear explicitly in the differential equation. Higherorder time derivative terms are easily generated. Turner s computational differentiation-based Object- Oriented Coordinate Embedding Algorithm (OCEA) provides numerical values the n-th order Jacobians. The analyst only codes f(x, t); OCEA uses hidden operator-overloaded software-based tools for automatically deriving, coding, and evaluating the Jacobian tensors required in Eq.(2). Without the availability of CD Eq.(3) only represents a theoretical solution for an exact series approximation. COMPUTATIONAL DIFFERENTIATION Computational differentiation is a specialized topic in applied mathematics and computer science for developing and fielding software tools for numerically evaluating partial derivative models. 1, 2 This paper makes use of Turner s Object-Oriented Coordinate Embedding Algorithm (OCEA) program for computing 1 st -4 th -order mixed partial derivative models. 3, 4, 5, 6 OCEA is a CD tool for computing arbitrarily complex partial derivative models. At compile time, OCEA uses the programmer s math model as a template for deriving, coding, and generating an executable for simultaneously compiling the simulation and sensitivity models. No symbolic or finite difference tools 2
3 are used for any of the gradient tensor calculations; all results are immediately processed to produce numerical results for all partial derivative orders. Each time a program is complied OCEA derives, assemblies, and codes the partial derivative model in the program executable, yielding numerical results that are accurate to the working level of precision for the machine. This approach benefits the user in four ways: (1) only the basic math model is programmed and checked out; (2) no analyst effort is required to either derive or code sensitivity models; (3) the user recaptures the development time normally required for developing sensitivity math models, coding, and validating nonlinear and high-order models; and (4) most importantly, math model changes are automatically handled each time the code is complied. Hand derived models are very vulnerable to model changes that can potentially force a new derivation and coding for the sensitivity model. OCEA transforms all math and intrinsic functions to embed multiple levels of the chain rule of calculus for building partial derivative models. OCEA consists of a suite of programs/modules that process the user supplied software math models as a template for generating high-order sensitivity models. Minimal modifications of the User s existing software models are required for enabling derivative-enhanced calculations. Users identify independent variables for the derivative calculations and define the variables for which partial derivatives are required (e.g., typically one is required to convert real variable(s) to an OCEA variable TYPE(EB), which instructs the complier that partial derivatives are computed for the TYPE(EB) variable(s). The software accepts user math models coded in standard FORTRAN 95/2003 language constructs. At compile time the partial derivative solution is assembled by identifying variable data types, intrinsic/lib functions, and inserting Function/Subroutine calls in the complier generated executable for computing the associated partial derivatives. From the user perspective, OCEA behaves as a language extension for FORTRAN 95/2003: its hidden operations exploit operator overloading and user-defined data types for handling all memory management details and derivative enhanced operations. It should not be surprising that OCEA introduces a computational overhead when compared to highly optimized hand coded partial derivative models. The trade off is this: if the basic math model requires X man-months to develop, then the derivation, coding, and validation of the sensitivity model can add 5X-10X man-months to the project development effort. Since the computer time required for deriving and compiling the sensitivity solution is measured in seconds vs. 5X-10X manmonths for an analyst, OCEA s impact is both clear and unambiguous for developing and solving real-world projects subject to man-power resource constraints. OCEA is particularly valuable in the normally fluid engineering design environment, where frequent design changes and -what ifexperiments must be carried out to fully explore the opportunities available in the notional system design space. The analyst always has the most up to date model. This is in stark contrast the case of a hand-derived model where even seemingly simple model changes can devastate all previous derivation efforts; thereby, forcing a restart for the derivation and coding effort from scratch for each new model is a daunting unwelcome task. OCEA s derivative enhanced variables are defined as abstract compound data objects, where objects such as the variable F, defined below, consist of the following list of concatenated data in computer memory: F := { F F 2 F n F } where m F denotes the m th order tensor gradient operation. The only variable visible to the analyst is F. Numerical values for the sub-object component values of the tensor gradient operators are obtained by using structure constructor designators (e.g., m F = F %T m, where T m denotes the tensor order). A detailed description of all of OCEA s capabilities are found in the software user manual. 7 3
4 CONSTRAINED MULTI-BODY DYNAMICAL SYSTEMS A problem of long standing in the science and engineering community is concerned with providing accurate integrations of math models that consist of both differential and algebraic equations. All numerical techniques produce drift in the constraint values that must be addressed for the simulations to remain physically meaningful. Many strategies have been proposed for addressing this problem, ranging from differential correction to defining artificial constraint differential equations. 8, 9, 10 Our goal is to exactly model the constraint through several time derivatives so that the constraint drift is minimized. Other researchers have considered developing CD-based integration methods; however, these efforts have not dealt with the time-varying system matrix derivative handled in this formulation. The new capability is critically important for enabling CD methods to handle real-world applications in engineering and science. To this end, the following constrained multi-body dynamics system is considered M q = f + [B] T λ subject to [B] q + b = 0 (4) where M denotes the mass matrix, q is the generalized coordinates of the system, ( q = dq dt, q = d 2 q ), f denotes all forces acting on the system, [B] denotes the velocity constraint matrix, λ dt 2 ( R Nc 1 ) denotes the constraint Lagrange multiplier, N c is the number of constraints. This system of equations is re-cast in the Jacobi form of the Routh-Voss equations as [ M B T B 0 ] [ ] [ ] q f = λ ḃ Ḃ q MẎ = F (5) Equation (5) completely describes the constrained system dynamics and recovers both the system acceleration and the negative of the constraint Lagrange multiplier. Equation (5) is known as the mass descriptor form of the constrained equations of motion. An exact series expansion for the system position and velocity is obtained by repeatedly differentiating Eq.(5) w.r.t. time, yielding the cascade of equations M Ẏ = F Ẏ = M 1 F MŸ + MẎ = F M Y... Ÿ = M 1 ( F MẎ) + 2 MŸ + MẎ = F... Y = M 1 ( F 2 MŸ MẎ) (6).. where the time derivatives of the mass descriptor matrix M and the force vector F are calculated using this equation Ψ = Ψ(q, q) Ψ = Ψ = N q N q ( Ψ q i ( Ψ q i + Ψ ) q i q i q i + Ψ... q i q i N q q i + j=1 + 2 Ψ q i q i q q j + 2 Ψ )) q i q j j q i q j ( 2 Ψ q i q i q j q j + 2 Ψ q i q j q i q j (7) 4
5 where Ψ can be replaced by M or F, and N q denotes to the number of generalized coordinates. The numerical solutions account for the curvature of the constraints through four derivatives. These equations permit the following exact Taylor expansions for the constrained response ( 1 R(t + h) R(t) + V(t)h + S Yh Yh 5) 2!Ẏh2 3!Ÿh3 4! 5! ( 1 V(t + h) V(t) + S 1!Ẏh Yh Yh 4) (8) 2!Ÿh2 3! 4! S = [I 0] Selecting Operator Yielding 6 th order prediction for the position and 5 th order prediction for the velocity. Several examples are presented for both unconstrained and constrained systems. Constraint drift during integration is expected to be minimized during the simulation of complex mechanical systems. Numerical experiments will be performed to determine the impact of the proposed high-order constraint tracking algorithm. NUMERICAL EXAMPLES Acceleration Constraints This section considers presents one of the common constrained problems in engineering. Figure (1) represents slider crank mechanism connected to a block mass that is allowed to slide on the horizontal plane. Therefore, we have one constraint equation that forces the M 3 mass moves horizontally. The system equations of motion can be derived using the Lagrangian approach assuming the bars are too thin and the mass center locates at the middle point of the body. The generalized coordinates of the system are given by (q 1 = θ, q 2 = θ 2 ). To start the derivation, one needs to write the position vector for each mass: Y M 1, L 1 M 2, L 2 M3 θ 1 θ 2 X Figure 1. Slider Crank Mechanism. r 1 = l 1 2 ( ) cos(q1 ) ; r sin(q 1 ) 2 = 2r 1 + l ( ) ( ) 2 cos(q2 ) cos(q2 ) ; r 2 sin(q 2 ) 3 = 2r 1 + l 2 sin(q 2 ) (9) The Lagrangian equations of motion can be written in the form: d ( L ) L dt q q = Q nc + [B] T λ (10) 5
6 where L(q, q, t) = T (q, q, t) V (q, t). Q nc is the non-conservative force vector; assumed zero. T and V are the kinetic and potential energy of the system, respectively. T = M i (ṙ i.ṙ i ); V = g M i l i sin(q i ) + M 3 g(l 1 sin(q 1 ) l 2 sin(q 2 )) (11) where g is the earth gravity constant (= 9.81m/s 2 ). And the system Pfaffian constraint is given by: ( ) q1 [B] q + b = 0 [l 1 cos(q 1 ) l 2 cos(q 2 )] = 0 (12) q 2 Following a few mathematical manipulations, the mass coefficient matrix and the force vector can be written as: (0.25M 1 + M 2 + M 3 )l 2 1 (0.5M 2 + M 3 )l 1 l 2 cos(q 1 + q 2 ) l 1 cos(q 1 ) M = (0.5M 2 + M 3 )l 1 l 2 cos(q 1 + q 2 ) (0.25M 2 + M 3 )l2 2 l 2 cos(q 2 ) l 1 cos(q 1 ) l 2 cos(q 2 ) 0 (0.5M 2 + M 3 )l 1 l 2 q 2 2 F sin(q 1 + q 2 ) (0.5M 1 + M 3 )l 1 gcos(q 1 ) = (0.5M 2 + M 3 )l 1 l 2 q 1 2sin(q 1 + q 2 ) (0.5M 2 M 3 )l 2 gcos(q 2 ) (13) l 1 q 1 2sin(q 1) l 2 q 2 2sin(q 2) We do not claim that this is the best way to integrate the multi-body dynamics. However, the advantages are noticeable if it is compared with the most common numerical integrator; Runge- Kutta (RK4). The solution is compared with Runge-Kutta (RK4) to demonstrate these advantages. Figure 2 shows the numerical integration results obtained by the two cases. Figure (a) shows the time histories of the state variables as integrated by the Analytic Integration Method. Figure (b) shows the % error between the solution obtained by Runge-Kutta and the analytic integration methods. Figures (c) and (d) show xy-trajectories as integrated by the Runge-Kutta method and Analytic Integration Method, respectively. One can notice the drift is less in the Analytic Integration case. Velocity and Acceleration Constraints The previous developments only constrain the solution at the acceleration level. The solution accuracy is further improved by constraining both the velocity and acceleration simultaneously. This is accomplished by generalizing the previous developments. Assembling the equation of motion and both the velocity constraint, and acceleration constraint, one obtains: M q B T λ = F M 0 B T q F B q = b B q + Ḃ q = ḃ 0 B 0 Λ = b (14) B Ḃ 0 λ ḃ The mass matrix is thus extended to include the velocity constraints coefficients. This produces a non-square non-invertable matrix. To fix this problem, however, recall the velocity correction term q + = q + M 1 B T Λ, which when introduced above yields M q B T λ = F M 0 B T q F B( q + M 1 B T Λ) = b B q + Ḃ( q + M 1 B T Λ) = ḃ 0 BM 1 B T 0 Λ = b B q B ḂM 1 B T 0 λ ḃ Ḃ q (15) 6
7 (a) (b) (c) (d) Figure 2. Acceleration Constraints Results. This system of equations is implemented and compared with the previous results, as well as with Runge-Kutta method. The Numerical Results are presented in Figure 3; (a) shows the constrained force and constrained momentum histories with time. (b) shows the % error between the solution obtained by Runge-Kutta and the analytic integration methods (both methods with velocity constraints). (e) and (f) show xy-trajectories as integrated by the Runge-Kutta method and Analytic Integration Method (both methods with velocity constraints), respectively. One can notice the drift is slightly reduced in the Analytic Integration case, as it is compared with the solution without velocity constraints. Adding the constrained velocity constraint is seen to further improve the solution accuracy. SUMMERY AND CONCLUSION A computational differentiation method is presented for analytically computing a Taylor series representation of the response nonlinear vector systems. Formulations are presented for constrained multi-body systems. Planer slider crank mechanism problem is presented that demonstrates the proposed analytic Taylor series integration of the constrained nonlinear ordinary vector differential equations produce very accurate results. Further research is required for establishing the utility of these new approaches for addressing engineering level of fidelity simulations. Constrained multibody systems showed very rigorous effect on enforcing the kinematic constraints and decrease the 7
8 (a) (b) (c) (d) Figure 3. Velocity and Acceleration Constraints Results. constraint drift level. REFERENCES [1] C. H. Bischof, A. Carle, P. D. Hovland, P. Khademi, and A. Mauer, ADIFOR 2.0 User s Guide (Revision D), tech. rep., Mathematics and Computer Science Division Technical Memorandum no. 192 and Center for Research on Parallel Computation Technical Report CRPC S, [2] A. Griewank, On Automatic Differentiation, Mathematical Programming (M. Iri and K. Tanabe, eds.), pp , Dordrecht: Kluwer Academic Publishers, [3] J. L. Junkins, J. D. Turner, and M. Majji, Generalizations and Applications of the Lagrange Implicit Function Theorem, Special Issue: The F. Landis Markley Astronautics Symposium, The Journal of the Astronautical Sciences, Vol. 57, January-June 2009, pp [4] Macsyma, Inc, Macsyma, Symbolic/numeric/graphical mathematics software: Mathematics and System Reference Manual, 16th ed., [5] J. Turner, Automated Generation of High-Order Partial Derivative Models, Vol. 41, AIAA, August 2003, pp [6] J. D. Turner, M. Majji, and J. L. Junkins, Keynote Paper: Fifth-Order Exact Analytic Continuation Numerical Integration Algorithm, Nanjing, China, Presented to International Conference on Computational and Experimental Engineering and Sciences, April [7] J. Turner, OCEA User Manual. Amdyn System, [8] E. Routh, Dynamics of a System of Rigid Bodies, Advanced Part. London:Macmillan/St. Martin: Reprinted 1955 by Dovr, New York, 6th ed., [9] E. Routh, Dynamics of a System of Rigid Bodies, Elementry Part. London:Macmillan/St. Martin: Reprinted 1960 by Dovr, New York, 7th ed.,
9 [10] H. Schaub and J. Junkins, Analytical mechanics of space systems. AIAA education series, American Institute of Aeronautics and Astronautics,
HIGH-ORDER STATE FEEDBACK GAIN SENSITIVITY CALCULATIONS USING COMPUTATIONAL DIFFERENTIATION
(Preprint) AAS 12-637 HIGH-ORDER STATE FEEDBACK GAIN SENSITIVITY CALCULATIONS USING COMPUTATIONAL DIFFERENTIATION Ahmad Bani Younes, James Turner, Manoranjan Majji, and John Junkins INTRODUCTION A nonlinear
More informationAdvanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich
Advanced Dynamics - Lecture 4 Lagrange Equations Paolo Tiso Spring Semester 2017 ETH Zürich LECTURE OBJECTIVES 1. Derive the Lagrange equations of a system of particles; 2. Show that the equation of motion
More informationRADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION
AAS 15-440 RADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION Austin B. Probe, * Brent Macomber,* Julie I. Read,* Robyn M. Woollands,* and John L. Junkins
More informationAN EMBEDDED FUNCTION TOOL FOR MODELING AND SIMULATING ESTIMATION PROBLEMS IN AEROSPACE ENGINEERING
AAS -8 AN EMBEDDED FUNCTION TOOL FOR MODELING AND SIMULATING ESTIMATION PROBLEMS IN AEROSPACE ENGINEERING INTRODUCTION D. Todd Griffith *, James D. Turner, and John L. Junkins An automatic differentiation-based
More informationLinear Feedback Control Using Quasi Velocities
Linear Feedback Control Using Quasi Velocities Andrew J Sinclair Auburn University, Auburn, Alabama 36849 John E Hurtado and John L Junkins Texas A&M University, College Station, Texas 77843 A novel approach
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 informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationAnalytical Mechanics. of Space Systems. tfa AA. Hanspeter Schaub. College Station, Texas. University of Colorado Boulder, Colorado.
Analytical Mechanics of Space Systems Third Edition Hanspeter Schaub University of Colorado Boulder, Colorado John L. Junkins Texas A&M University College Station, Texas AIM EDUCATION SERIES Joseph A.
More informationUSING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS
(Preprint) AAS 12-629 USING CARLEMAN EMBEDDING TO DISCOVER A SYSTEM S MOTION CONSTANTS John E. Hurtado and Andrew J. Sinclair INTRODUCTION Although the solutions with respect to time are commonly sought
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 informationHigh-Order Representation of Poincaré Maps
High-Order Representation of Poincaré Maps Johannes Grote, Martin Berz, and Kyoko Makino Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA. {grotejoh,berz,makino}@msu.edu
More informationESM 3124 Intermediate Dynamics 2012, HW6 Solutions. (1 + f (x) 2 ) We can first write the constraint y = f(x) in the form of a constraint
ESM 314 Intermediate Dynamics 01, HW6 Solutions Roller coaster. A bead of mass m can slide without friction, under the action of gravity, on a smooth rigid wire which has the form y = f(x). (a) Find the
More informationInfinite series, improper integrals, and Taylor series
Chapter 2 Infinite series, improper integrals, and Taylor series 2. Introduction to series In studying calculus, we have explored a variety of functions. Among the most basic are polynomials, i.e. functions
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 informationPHY321 Homework Set 10
PHY321 Homework Set 10 1. [5 pts] A small block of mass m slides without friction down a wedge-shaped block of mass M and of opening angle α. Thetriangular block itself slides along a horizontal floor,
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationLecture 21. MORE PLANAR KINEMATIC EXAMPLES
Lecture 21. MORE PLANAR KINEMATIC EXAMPLES 4.5c Another Slider-Crank Mechanism Figure 4.24 Alternative slider-crank mechanism. Engineering-analysis task: For = ω = constant, determine φ and S and their
More informationChapter 4 Statics and dynamics of rigid bodies
Chapter 4 Statics and dynamics of rigid bodies Bachelor Program in AUTOMATION ENGINEERING Prof. Rong-yong Zhao (zhaorongyong@tongji.edu.cn) First Semester,2014-2015 Content of chapter 4 4.1 Static equilibrium
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Numerics Michael Bader SCCS Technical University of Munich Summer 018 Recall: Molecular Dynamics System of ODEs resulting force acting on a molecule: F i = j
More informationSwing-Up Problem of an Inverted Pendulum Energy Space Approach
Mechanics and Mechanical Engineering Vol. 22, No. 1 (2018) 33 40 c Lodz University of Technology Swing-Up Problem of an Inverted Pendulum Energy Space Approach Marek Balcerzak Division of Dynamics Lodz
More informationFigure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody
Lecture 27. THE COMPOUND PENDULUM Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody diagram The term compound is used to distinguish the present
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
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 informationClassical Mechanics Review (Louisiana State University Qualifier Exam)
Review Louisiana State University Qualifier Exam Jeff Kissel October 22, 2006 A particle of mass m. at rest initially, slides without friction on a wedge of angle θ and and mass M that can move without
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 informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationRobotics 2 Robot Interaction with the Environment
Robotics 2 Robot Interaction with the Environment Prof. Alessandro De Luca Robot-environment interaction a robot (end-effector) may interact with the environment! modifying the state of the environment
More informationROBOTICS Laboratory Problem 02
ROBOTICS 2015-2016 Laboratory Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with or without friction along the horizontal
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
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 informationLarge Steps in Cloth Simulation. Safeer C Sushil Kumar Meena Guide: Prof. Parag Chaudhuri
Large Steps in Cloth Simulation Safeer C Sushil Kumar Meena Guide: Prof. Parag Chaudhuri Introduction Cloth modeling a challenging problem. Phenomena to be considered Stretch, shear, bend of cloth Interaction
More informationDynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18
Dynamics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Dynamics Semester 1, 2016-17 1 / 18 Dynamics Dynamics studies the relations between the 3D space generalized forces
More informationIYGB Mathematical Methods 1
IYGB Mathematical Methods Practice Paper B Time: 3 hours Candidates may use any non programmable, non graphical calculator which does not have the capability of storing data or manipulating algebraic expressions
More informationLecture 31. EXAMPLES: EQUATIONS OF MOTION USING NEWTON AND ENERGY APPROACHES
Lecture 31. EXAMPLES: EQUATIONS OF MOTION USING NEWTON AND ENERGY APPROACHES Figure 5.29 (a) Uniform beam moving in frictionless slots and attached to ground via springs at A and B. The vertical force
More informationIntroduction to Robotics
J. Zhang, L. Einig 277 / 307 MIN Faculty Department of Informatics Lecture 8 Jianwei Zhang, Lasse Einig [zhang, einig]@informatik.uni-hamburg.de University of Hamburg Faculty of Mathematics, Informatics
More informationUsing Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver
Using Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver Jeffrey P. Thomas, Earl H. Dowell, and Kenneth C. Hall Duke University, Durham, NC 27708 0300 A novel nonlinear
More informationMathematics for Intelligent Systems Lecture 5 Homework Solutions
Mathematics for Intelligent Systems Lecture 5 Homework Solutions Advanced Calculus I: Derivatives and local geometry) Nathan Ratliff Nov 25, 204 Problem : Gradient and Hessian Calculations We ve seen that
More informationRobotics. Dynamics. University of Stuttgart Winter 2018/19
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler, joint space control, reference trajectory following, optimal operational
More informationThe Implicit Function Theorem with Applications in Dynamics and Control
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 2010, Orlando, Florida AIAA 2010-174 The Implicit Function Theorem with Applications in Dynamics
More informationCOMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.
BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationA Sliding Mode Controller Using Neural Networks for Robot Manipulator
ESANN'4 proceedings - European Symposium on Artificial Neural Networks Bruges (Belgium), 8-3 April 4, d-side publi., ISBN -9337-4-8, pp. 93-98 A Sliding Mode Controller Using Neural Networks for Robot
More informationReview of Engineering Dynamics
Review of Engineering Dynamics Part 1: Kinematics of Particles and Rigid Bodies by James Doane, PhD, PE Contents 1.0 Course Overview... 4.0 Basic Introductory Concepts... 4.1 Introduction... 4.1.1 Vectors
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 informationChapter 11. Special Relativity
Chapter 11 Special Relativity Note: Please also consult the fifth) problem list associated with this chapter In this chapter, Latin indices are used for space coordinates only eg, i = 1,2,3, etc), while
More informationConstraint Based Control Method For Precision Formation Flight of Spacecraft AAS
Constraint Based Control Method For Precision Formation Flight of Spacecraft AAS 06-122 Try Lam Jet Propulsion Laboratory California Institute of Technology Aaron Schutte Aerospace Corporation Firdaus
More informationChapter 1 Mathematical Preliminaries and Error Analysis
Numerical Analysis (Math 3313) 2019-2018 Chapter 1 Mathematical Preliminaries and Error Analysis Intended learning outcomes: Upon successful completion of this chapter, a student will be able to (1) list
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationTaylor series. Chapter Introduction From geometric series to Taylor polynomials
Chapter 2 Taylor series 2. Introduction The topic of this chapter is find approximations of functions in terms of power series, also called Taylor series. Such series can be described informally as infinite
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... Introduce the topic of numerical methods Consider the Error analysis and sources of errors Introduction A numerical method which
More informationGradient Descent. Dr. Xiaowei Huang
Gradient Descent Dr. Xiaowei Huang https://cgi.csc.liv.ac.uk/~xiaowei/ Up to now, Three machine learning algorithms: decision tree learning k-nn linear regression only optimization objectives are discussed,
More informationSlovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS
Slovak University of Technology in Bratislava Institute of Information Engineering, Automation, and Mathematics PROCEEDINGS of the 18 th International Conference on Process Control Hotel Titris, Tatranská
More informationCovariant Formulation of Electrodynamics
Chapter 7. Covariant Formulation of Electrodynamics Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 11, and Rybicki and Lightman, Chap. 4. Starting with this chapter,
More informationORBIT 14 The propagator. A. Milani et al. Dipmat, Pisa, Italy
ORBIT 14 The propagator A. Milani et al. Dipmat, Pisa, Italy Orbit 14 Propagator manual 1 Contents 1 Propagator - Interpolator 2 1.1 Propagation and interpolation............................ 2 1.1.1 Introduction..................................
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
More informationSolving scalar IVP s : Runge-Kutta Methods
Solving scalar IVP s : Runge-Kutta Methods Josh Engwer Texas Tech University March 7, NOTATION: h step size x n xt) t n+ t + h x n+ xt n+ ) xt + h) dx = ft, x) SCALAR IVP ASSUMED THROUGHOUT: dt xt ) =
More informationUsing Spreadsheets to Teach Engineering Problem Solving: Differential and Integral Equations
Session 50 Using Spreadsheets to Teach Engineering Problem Solving: Differential and Integral Equations James P. Blanchard University of Wisconsin - Madison ABSTRACT Spreadsheets offer significant advantages
More informationDecentralized PD Control for Non-uniform Motion of a Hamiltonian Hybrid System
International Journal of Automation and Computing 05(2), April 2008, 9-24 DOI: 0.007/s633-008-09-7 Decentralized PD Control for Non-uniform Motion of a Hamiltonian Hybrid System Mingcong Deng, Hongnian
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More informationTwo-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra
Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Politecnico di Milano Department of Aerospace Engineering Milan, Italy Taylor Methods and Computer Assisted Proofs
More informationDefinition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.
5 Vector fields Last updated: March 12, 2012. 5.1 Definition and general properties We first need to define what a vector field is. Definition 5.1. A vector field v on a manifold M is map M T M such that
More informationDirect Optimal Control and Costate Estimation Using Least Square Method
21 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July 2, 21 WeB22.1 Direct Optimal Control and Costate Estimation Using Least Square Method Baljeet Singh and Raktim Bhattacharya
More informationNONLINEAR EQUATIONS AND TAYLOR S THEOREM
APPENDIX C NONLINEAR EQUATIONS AND TAYLOR S THEOREM C.1 INTRODUCTION In adjustment computations it is frequently necessary to deal with nonlinear equations. For example, some observation equations relate
More informationOPTIMAL SPACECRAF1 ROTATIONAL MANEUVERS
STUDIES IN ASTRONAUTICS 3 OPTIMAL SPACECRAF1 ROTATIONAL MANEUVERS JOHNL.JUNKINS Texas A&M University, College Station, Texas, U.S.A. and JAMES D.TURNER Cambridge Research, Division of PRA, Inc., Cambridge,
More informationAdaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties
Australian Journal of Basic and Applied Sciences, 3(1): 308-322, 2009 ISSN 1991-8178 Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties M.R.Soltanpour, M.M.Fateh
More informationAN INTRODUCTION TO LAGRANGE EQUATIONS. Professor J. Kim Vandiver October 28, 2016
AN INTRODUCTION TO LAGRANGE EQUATIONS Professor J. Kim Vandiver October 28, 2016 kimv@mit.edu 1.0 INTRODUCTION This paper is intended as a minimal introduction to the application of Lagrange equations
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationTrajectory Planning from Multibody System Dynamics
Trajectory Planning from Multibody System Dynamics Pierangelo Masarati Politecnico di Milano Dipartimento di Ingegneria Aerospaziale Manipulators 2 Manipulator: chain of
More information7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
More informationPart of the advantage : Constraint forces do no virtual. work under a set of virtual displacements compatible
FORCES OF CONSTRAINT Lagrangian formalism : Generalized coordinate Minimum set of Eqns Part of the advantage : Constraint forces do no virtual work under a set of virtual displacements compatible with
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationMultibody simulation
Multibody simulation Dynamics of a multibody system (Euler-Lagrange formulation) Dimitar Dimitrov Örebro University June 16, 2012 Main points covered Euler-Lagrange formulation manipulator inertia matrix
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationAttitude Determination for NPS Three-Axis Spacecraft Simulator
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 6-9 August 4, Providence, Rhode Island AIAA 4-5386 Attitude Determination for NPS Three-Axis Spacecraft Simulator Jong-Woo Kim, Roberto Cristi and
More informationEN Nonlinear Control and Planning in Robotics Lecture 2: System Models January 28, 2015
EN53.678 Nonlinear Control and Planning in Robotics Lecture 2: System Models January 28, 25 Prof: Marin Kobilarov. Constraints The configuration space of a mechanical sysetm is denoted by Q and is assumed
More informationSECTION C: CONTINUOUS OPTIMISATION LECTURE 9: FIRST ORDER OPTIMALITY CONDITIONS FOR CONSTRAINED NONLINEAR PROGRAMMING
Nf SECTION C: CONTINUOUS OPTIMISATION LECTURE 9: FIRST ORDER OPTIMALITY CONDITIONS FOR CONSTRAINED NONLINEAR PROGRAMMING f(x R m g HONOUR SCHOOL OF MATHEMATICS, OXFORD UNIVERSITY HILARY TERM 5, DR RAPHAEL
More informationMilne s Spreadsheet Calculator Using VBA with Excel Programming for solving Ordinary Differential Equations
Milne s Spreadsheet Calculator Using VBA with Excel Programming for solving Ordinary Differential Equations Miss. S. P. Hingmire 1, Dr. M.R. Gosavi 2, Dr. N. A.Patil 3 1 Research Scholar, Asst.Professor,
More informationIn most robotic applications the goal is to find a multi-body dynamics description formulated
Chapter 3 Dynamics Mathematical models of a robot s dynamics provide a description of why things move when forces are generated in and applied on the system. They play an important role for both simulation
More informationCHAPTER 1 INTRODUCTION TO NUMERICAL METHOD
CHAPTER 1 INTRODUCTION TO NUMERICAL METHOD Presenter: Dr. Zalilah Sharer 2018 School of Chemical and Energy Engineering Universiti Teknologi Malaysia 16 September 2018 Chemical Engineering, Computer &
More informationKinematics. Chapter Multi-Body Systems
Chapter 2 Kinematics This chapter first introduces multi-body systems in conceptual terms. It then describes the concept of a Euclidean frame in the material world, following the concept of a Euclidean
More informationIn this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:
Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationCOSY INFINITY Version 7
COSY INFINITY Version 7 Kyoko Makino and Martin Berz Department of Physics and Astronomy and National Superconducting Cyclotron Laboratory Michigan State University, East Lansing, MI 48824 Abstract An
More informationNumerical Data Fitting in Dynamical Systems
Numerical Data Fitting in Dynamical Systems A Practical Introduction with Applications and Software by Klaus Schittkowski Department of Mathematics, University of Bayreuth, Bayreuth, Germany * * KLUWER
More informationLINEAR AND NONLINEAR PROGRAMMING
LINEAR AND NONLINEAR PROGRAMMING Stephen G. Nash and Ariela Sofer George Mason University The McGraw-Hill Companies, Inc. New York St. Louis San Francisco Auckland Bogota Caracas Lisbon London Madrid Mexico
More informationAN EVALUATION SCHEME FOR THE UNCERTAINTY ANALYSIS OF A CAPTIVE TRAJECTORY SYSTEM
24th INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES AN EVALUATION SCHEME FOR THE UNCERTAINTY ANALYSIS OF A CAPTIVE TRAJECTORY SYSTEM Sean Tuling Defencetek, CSIR Keywords: Uncertainty Analysis, CTS
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationOptimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints
Optimization using Calculus Optimization of Functions of Multiple Variables subject to Equality Constraints 1 Objectives Optimization of functions of multiple variables subjected to equality constraints
More informationConvergence of a Gauss Pseudospectral Method for Optimal Control
Convergence of a Gauss Pseudospectral Method for Optimal Control Hongyan Hou William W. Hager Anil V. Rao A convergence theory is presented for approximations of continuous-time optimal control problems
More informationCase Study: The Pelican Prototype Robot
5 Case Study: The Pelican Prototype Robot The purpose of this chapter is twofold: first, to present in detail the model of the experimental robot arm of the Robotics lab. from the CICESE Research Center,
More informationLeast-squares Solutions of Linear Differential Equations
1 Least-squares Solutions of Linear Differential Equations Daniele Mortari dedicated to John Lee Junkins arxiv:1700837v1 [mathca] 5 Feb 017 Abstract This stu shows how to obtain least-squares solutions
More informationModule III: Partial differential equations and optimization
Module III: Partial differential equations and optimization Martin Berggren Department of Information Technology Uppsala University Optimization for differential equations Content Martin Berggren (UU)
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationDynamics and Controls of a Generalized Frequency Domain Model Flexible Rotating Spacecraft
SpaceOps Conferences 5-9 May 24, Pasadena, CA SpaceOps 24 Conference.254/6.24-797 Dynamics and Controls of a Generalized Frequency Domain Model Flexible Rotating Spacecraft Tarek A. Elgohary, James D.
More informationIntroduction to gradient descent
6-1: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction to gradient descent Derivation and intuitions Hessian 6-2: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction Our
More informationA Study on Linear and Nonlinear Stiff Problems. Using Single-Term Haar Wavelet Series Technique
Int. Journal of Math. Analysis, Vol. 7, 3, no. 53, 65-636 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ijma.3.3894 A Study on Linear and Nonlinear Stiff Problems Using Single-Term Haar Wavelet Series
More information