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 - coordinate transformation, angular velocity vector, description of velocity and acceleration in relatively moving frames. Euler angles, Review of methods of momentum and angular momentum of system of particles, inertia tensor of rigid body. Dynamics of rigid bodies - Euler's equation, application to motion of symmetric tops and gyroscopes and problems of system of bodies. Kinetic energy of a rigid body, virtual displacement and classification of constraints. D Alembert s principle. Introduction to generalized coordinates, derivation of Lagrange's equation from D Alembert s principle. Small oscillations, matrix formulation, Eigen value problem and numerical solutions. Modelling mechanical systems, Introduction to MATLAB, computer generation and solution of equations of motion. Introduction to complex analytic functions, Laplace and Fourier transform. PID controllers, Phase lag and Phase lead compensation. Analysis of Control systems in state space, pole placement, computer simulation through MATLAB. 1 Contents Focus on 4 Introduction to modeling of mechanical systems 4 Modeling fundamentals of elements of mechanical systems Inertia element (mass) Energy storage elements (spring) Energy dissipation elements (dampers) 4 Examples and exercises 2
Introduction 4 A model is an abstraction of the physical world Used for analysis and design, commonly before realization of the physical system Can be obtained from first principles or experimentally The modeling purpose determines level of abstraction Can be complex enough, but no more than the physical object 4 Model derivation from first principles Input Using physical laws to derive the models Models based on first principles provides better undertanding Such models can also use empirical data to determine parameters, and validate the model System Output Give known input (empirical data) and observe output à fit model to your data - Applicable for complex systems and blackbox systems 3 Introduction 4 Derivation of Differential Equations (Des) is an important part of modeling of mechanical systems. 4 Sources of such DEs Newton s 2 nd Law F = ma Euler s 2 nd Law M = Iα Hooke s Law F = kx Lagrangian Hamiltonian d dt (L (q j - (L (q j = 0 dx dt - (H (q j = 0 4
Introduction Modeling of static and dynamic systems Static systems Relationship is static (not changing with time) Dynamic systems Relationship changes with time, past inputs and initial conditions influence current output (à it is dynamic) Only current input determines the output, which reacts immediately Relationship is represented by algebraic equations Output takes time to respond to inputs Relationship is represented by differential equations Input Output (system response) System Volt, 5v (Ang. Speed = 200 rpm) El-motor Force and motion sources cause elements to respond 5 1) Energy storage elements q Inertia elements q Stiffness (spring) elements 2) Energy dissipation elements q q Dampers Friction Mechanical systems can be either translational or rotational, for instance, for energy storage elements. Translational system - Inertia element: mass, m - Stiffness element: k For rotational system: - Inertia element: Mass moment of inertia, I - Stiffness element: k t 6 These are passive (nonenergy producing) devices
1) Inertia elements q Inertia elements store energy in the form of KE q Commonly components are not added to a mech. system for the purpose of adding inertia q Since all materials have mass, however, the mass or inertia element often may represents an undesirable effect in the system. q There are some applications in which mass itself serves a useful function, e.g., flywheels. Response of engine with flywheel Flywheels are used as energy-storage devices or as a means of smoothing out speed fluctuations in in mechanical systems such as IC engines or other machines 7 Response of engine without flywheel 1) Inertia elements... modelling q Inertia elements are modelled as follows Particle (point mass), m Rid body of mass, m Symbol for translational element T, θ F, x Torque θ = T I F, x m Moment of inertia x = F m Mathematical model from Newton s 2 nd Law Note that, for a rigid body in translation, every particle in the body has the same velocity and acceleration, while in flexible bodies the velocities and/or accelerations can vary Symbol for rotational element Mathematical model from Newton s 2 nd Law Note: No real rigid body exists when being accelerated, and the inertia element is a model, not a real/physical object. 8
1) Inertia elements... Mass moment of inertia for typical geometries q Particle of mass m, rotating at end of a massless rod of length l: I = ml 2 q A cylinder (r, m), rotating about center axis: I = 1 2 mr2 q Solid sphere (r, m), rotating about center axis: I = 2 3 mr2 I = m r 2 dm q Uniform rod (l, m), rotating about its center: q Uniform rod (l, m), rotating about its end: I = 1 3 ml2 I = 1 12 ml2 Note: moment of inertia of rigid bodies with complex geometry are determined experimentally 9 Inertia elements... Stiffness elements (springs) Real springs - Are neither pure nor ideal - Can have inertia and friction Pure spring is a mathematical model, - It is not a real device - It is ideal: massless, no damping and linear Force: F = k(x 1 x 2 ) = kx For ideal spring element Torque: T = k t (θ 1 θ 2 ) = ktθ Note: The inverse of stiffness is called compliance (c) / flexibility F x x F k Energy: E = 1 2 kx2 10 c
Inertia elements... Stiffness elements (springs) Real springs are not purely linear - Linearization of a nonlinear spring - Taylor s expansion F(x) = F(x A ) + (x x 0 ). EF EG 0 + G HG A I! Ignoring higher order derivatives è F(x) - F(x A ) = (x x 0 ). è F = kx EF EG 0 I E K F EG K 0 + L M! F(x A ) E N F EG N 0 +... Nonlinear F(x) 11 x A Inertia elements... Stiffness elements (springs) Real springs are not purely linear 4 Coil spring 4 Cantilever beam spring 4 Tension rod spring 4 Torsion bar spring 4 clamped-end beam spring 4 ring spring 4 rubber spring (shock mount) 4 hydraulic (oil) spring 4 pneumatic (air) spring k = Gd 5 8nD 3 k = 3EI L 3 k t = 12 k = AE L EI L
Inertia elements... Stiffness elements (springs) Key steps to modelling of mechanical systems 1) Establish an inertial (fixed) coordinate system 2) Identify discrete mechanical system elements (masses, springs, dampers) and isolate them 3) Determine min. number of variables needed to uniquely define the configuration of system and subtract constraints 4) Draw the free body diagram for each inertia element 5) Write equations that relate each element s loading to deformation 6) Apply Newton s 2nd Law to each element: Translational motion: Newton s 2 nd Law: F = ma Rotational motion: Euler s 2 nd Law: 13 T = Iα Inertia elements... Stiffness elements (springs) Example: 2 mass and 3 spring system* 14
Energy dissipation elements: Damping elements 4 Dampers converts the mechanical energy to thermal energy. 4 Two main forms Friction damping (coulomb dry friction damping) Viscoelastic damping or viscous damping (most common type) F Static friction F Kinetic friction F x x x Linear damper Coulomb friction Drag damper F = Cx F = C Sin(x ) 15 F = Cx 2 4 Applications: A damper element is used to model a device designed into a system (e.g., automotive shock absorbers) Sprang mass B Unsprang mass 16
Energy dissipation elements: Damping elements... 4 Damper force or torque is directly proportional to the relative velocity of its two ends. 4 F D = B d x 1 x dt 2 = B dx = Bx for translational motion dt 4 T D = B d (θ 1 θ dt 2) = B dθ dt = Bθ - for rotational motion 4 Where C is the damping coefficient In some textbooks B is used as a damping coeff. C Symbol of a damper Power dissipated Fx = C(x ) 2 (translational motion) 17 Example 2*: mass spring and damper system C 1 C 2 C 3 18
Exercise: Modelling the suspension system of a vehicle k 2 x 2 m 2 C x 1 m 1 m 1 = Unsprung mass m 2 = Spring mass, a quarter of the car s mass C = damping coeff. of shock absorber k 1 = Tire stiffness (damping in tire is neglected) k 2 = Stiffness of suspension coil spring x 0 = motion of the ground stiffness x 1 = coordinate for motion of unsprung mass x 2 = coordinate for motion of sprung mass k 1 x 0 Derive the DEs that describe the motion of the suspension system and fomulate it in matrix form 19 This lecture has focused on Summary 4 Explanation of the building blocks of modeling mechanical systems: - Inertia elements, - energy storage (spring) elements and - energy dissipation (damper) elements 4 Use of Newton s 2 nd Law and Euler s 2 nd Law to develop mathematical models of the modeling elements Next: Modelling mechanical systems, Introduction to MATLAB, computer generation and solution of equations of motion.. 20