Matlab Examples. Dayalbagh Educational Institute. From the SelectedWorks of D. K. Chaturvedi Dr.
|
|
- Doris Morton
- 6 years ago
- Views:
Transcription
1 Dayalbagh Educational Institute From the SelectedWorks of D. K. Chaturvedi Dr. Winter July 10, 2012 Matlab Examples D. K. Chaturvedi, Dr., Dayalbagh Educational Institute Available at:
2 EXERCISES ON MATLAB Prof. D.K. Chaturvedi Dayalbagh Educational Institute Dayalbagh, Agra (U.P.) Phone: DKC-31
3 Problem-1 The math behind bungee jumping This demonstration shows how to use MATLAB to model a simple physics problem faced by a college student. During spring break, John Smith wants to go bungee jumping. John has to determine which elastic cord is best suited for his weight. Elastic Cord A B C Spring Constant =5 N/m =40 N/m =500 N/m The air resistance that the bungee jumper faces is R= a1*v - a2* v *v Where A1=1 A2=1 The length of the unstretched cord is 30m. The bungee jumpers is 80m above the ground. DKC-32
4 To solve this physics problem, we need to: 1. Determine all the forces acting ON the body. 2. Draw a free body diagram. 3. Apply Newton's second law. 4. Solve the equation. 1. Determine the forces acting ON the body. Weight (W): W = m*g m = 90 kg g = 10 m/s^2 Air Resistance (R): R=a1*v+a2* v *v a1=1 a2=1 v=dx/dt Force from the elastic cord (Fe): Fe= k*x if x>0 0 if x<0 DKC-33
5 2. Draw the free body diagram. Please note that we have selected downwards as the positive axis. 3. Apply Newton's second law. Net forces=m*a W-R-Fe=m*a mg-fe-a1*v+a2* v *v=m*a where v=dx/dt and a=dv/dt Solution: The math behind bungee jumping We need to use the ODE solvers. The ODE solvers solve initial value problems for ordinary differential equations (ODEs). In this example we are using ODE45. The syntax for ODE45 is: [T,Y] = ode45(odefun,tspan,y0,options,p1,p2...) odefun A function that evaluates the right-hand side of the differential equations. All solvers solve systems of equations in the form y = f(t,y) tspan A vector specifying the interval of integration, [t0,tf]. To obtain solutions at specific times (all increasing or all decreasing), use tspan = [t0,t1,...,tf].y0a vector of initial conditions. Options Optional integration argument created using the DKC-34
6 odeset function. See odeset for details. p1,p2.. Optional parameters that the solver passes to odefun and all the functions specified in options First we need to create the odefun function: We need to rewrite our equation so that it will be in the form of y'=f(t,y) Assume : X1=x X2=dx/dt Therefore X1dot = x2; X2dot = g - Fe/m - a1/m*x2 - a2* x2 *x2 function dxdt = bungeeode (t,x,k) m=90; g=10; a1=1; a2=1; W=m*g; R= a1*x(2)+a2*abs(x(2))*x(2); if x(1)>0 Fe = k*x(1); else Fe = 0; end dxdt= [ x(2) ; (W-Fe-R)/m]; DKC-35
7 The time span can be from 0 to 50 s. To use the default options, we assign the option as [ ]. In a script file: figure [t,xsol]=ode45(@bungeeode,[0 50],[-30 0],[],5); plot(t,50-xsol(:,1)) figure [t,xsol]=ode45(@bungeeode,[0 50],[-30 0],[],40); plot(t,50-xsol(:,1)) DKC-36
8 DKC-37
9 Result: The best choice is ELASTIC CORD B. DKC-38
10 Problem 2 The physics of baseball A pitcher throws a baseball at 50m/s at an angle of 36.7 degrees in the air. How far will the ball go before it hits the ground? How high will the ball go? How long does it take to hit the ground? There is more than one way to solve this problem: 1. Use Calculus 2. Use MATLAB ODE solvers 3. Solve it symbolically Using MATLAB ODE solvers: function dxdt = trajectoryode(t, x) dxdt(1)= x(2); dxdt(2,1)= 0; dxdt(3,1) = x(4); dxdt(4,1) = -9.8; MATLAB script: v0= 20; theta = deg2rad(36.7); v0x = v0*cos(theta); v0y = v0*sin(theta); x0=0; y0=0; [t,xsol]=ode45(@trajectoryode,[0 2.5],... [x0 v0x y0 v0y],[]); x=xsol(:,1); vx=xsol(:,2); y=xsol(:,3); vy=xsol(:,4); plot(x,y) axis([ ]) hold on for i= 1:size(x) plot(x,y,'+') pause(0.2) end axis([ ]) plot(t,y) plot(t,vy) DKC-39
11 Solving the question symbolically: x0=0; y0=0; v0= 20; theta = deg2rad(36.7); v0x = v0*cos(theta); v0y = v0*sin(theta); syms t ay=-9.8; vy=int(ay,t)+v0y; y= int(vy,t)+y0; ax=0; vx=int(ax,t)+v0x; x= int(vx,t)+x0; To reduce number of digits use the function VPA: Y=vpa(y,3); X=vpa(x,3); Pretty(X) Pretty(Y) 16.0 t t t As expected: Ezplot(X,[0 2.5]) Ezplot(Y,[0 2.5]) Time to reach back to the ground: solve(y) ans = [0] [ ] Range=subs(X,t, ) For maximum height: Solve(vy) DKC-40
12 ans = MaxHeight=subs(Y,t, ) Problem 7 Trajectory of Parachute % Parachute Trajectory clear all; Vx=350; Vy=30; dt=0.01; t=0; tsim=10; rhow= ; Cd=0.7; S=19.625; W=50; G=9.81; x=0; H=250; Drag=0; n=round((tsim-t)/dt); for i=1:n X1(i,:)=[t, Vx, Vy, x, H, Drag]; V=sqrt(Vx^2+Vy^2); Q=(rhow*V^2)/2; dvx=-(q*cd*s*g*vx)/(w*v); dvy=-(q*cd*s*g*vy)/(w*v)-g; dx=(vx+dt*dvx/2); dh=(vy+dt*dvy/2); Vx=Vx+dt*dVx; Vy=Vy+dt*dVy; x=x+dt*dx; H=H+dt*dh; Drag=-(W/G)*(dVx/Vx)*V; t=t+dt; end subplot(2,2,1) plot(x1(:,1),x1(:,2)); title('vx') Problem: Energy Absorber System % Energy Absorber clear all; M=40000; J1=131; J2=131; K1= /17; K2= /17; K3=330*10^3; C1=16.427; C2=16.427; DKC-41
13 R1=0.75; R2=0.75; DTHETA1=0; DTHETA2=0; T=0.0087; X=0.0; DX=87.3; L1=17 L2=17; B=0; NUM=0; t=0; TSIM=6; DT=0.01; ENG=0; THETA1=0;THETA2=0;A=0; N1=1;K=1;F=1;EXTN=0; EXTN1=0; X3=0; THETA3=0; P=0; D=0; DTHETA3=0; N=0; EC=0; TR=0; for i=t:dt:tsim DDX=-((K1*(X-sin(A)*EXTN))+(K2*(X-sin(A)*EXTN))*sin(A))/M; DDTHETA1=-((C1*DTHETA1^2+K1*((sin(A)*EXTN-X)*R1))+D)/J1; DDTHETA2=-((C2*DTHETA2^2+K1*((sin(A)*EXTN-X)*R2))+D)/J2; % K*((sin(A)*extn-X*R2))J2; DTHETA1=DTHETA1+DT*DDTHETA1; DTHETA2=DTHETA2+DT*DDTHETA2; if X<=6.0 DTHETA1=0; DTHETA2=0; end DTHETA3=DT*DDTHETA1; THETA1=THETA1+DT*DTHETA1; THETA2=THETA2+DT*DTHETA2; THETA3=DT*DTHETA1; DX=DX+DT*DDX; X=X+DT*DX; X3=DT*DX; EXTN=EXTN+R1*THETA3; EXTN1=(X-sin(A)*EXTN); P=(C1*DTHETA1.^2-K1*((sin(A)*EXTN-X)*R1))/K3; D=(THETA3-P)*K3; DRTR=(C1*DTHETA1.^2); dral=(ddtheta2*j1)/r1; A=atan(X/30.25); TAPELOAD1=(K1*EXTN1); GLOAD=DDX/9.8; B=sqrt(X^ ^2); L1=B+30; L2=B+30; J1=J1/1.000;%J1-(L1*0.68)*(R1^ )); %J1= J2=J2/1.000;%J2-(L2*0.68)*(R2^ )); %J2/1.0008; F=(M*DDX); %((K1*(X-SIN(A)*SXTN))+(K2*(X-SIN(A)*EXTN)))*SIN(A); %abs(m*ddx) ENG=(0.5*F*X); V=R1*DTHETA1; N=DTHETA1*30/pi; EC=EC+2*(C1*DTHETA1^2+12^3*EXTN1*R1)*DTHETA1*DT; TR=-F*sin(A)*R1; if (-THETA1/(pi/2)+N1)<= & R1>=0.190 DKC-42
14 R1=R1-T/4; %(T*THETA3/2*pi); R2=R2-T/4; %(T*THETA3/2*pi); N1=N1+1; end dddx=-ddx/9.8; time1(k,:)=[i]; X1(K,:)=[X,DX,dddx, K1, DTHETA1,EC, TAPELOAD1]; K=K+1; if (X<2) TAPELOAD1=0; K1=0; end if (X>=2 & X<=5) K1=0.3*stif1(TAPELOAD1); %elseif X>2&X<5 K1=125695; %elseif X>5&X<20 K1=114695; %elseif X>20&X<42 K1=129308; elseif (X>5 & X<=15) K1=0.15*stif1(TAPELOAD1); end K2=K1; else K1= /L1; % if DTHETA1<60, C1=C1/1.1; % elseif DTHETA1>60 & DTHETA1<90 C1=C1*1.01; else C1=9.18; %end end X1(:,1)=time1; subplot(2,2,1) plot(x1(:,1),x1(:,2)) title('vel vs dis') xlabel('distance, m') ylabel('velocity') subplot(2,2,2) plot(x1(:,1),x1(:,3)) title('decelration') xlabel('distance, m') ylabel('g-load') subplot(2,2,3) plot(x1(:,1),x1(:,5)) title('angular Velocity') xlabel('distance, m') ylabel('rad/sec') subplot(2,2,4) plot(x1(:,1),x1(:,7)) title('tape load') xlabel('distance, m') ylabel('joules') DKC-43
15 Problem - 3 Model a basic cruise control system Let us model a basic cruise control system. We first assume that the Inertia of the vehicle wheels are negligible, and it is understood that the viscous friction created from the car's speed is completely in the direction of the opposing motion of the automobile. Therefore, our cruise control problem is now reduced to a single mass and a constant damper system as shown in Figure 1. Figure 1: Simple Cruise Control System. Given: b = 55 N*sec/m (damper constant), m = 1220 kg (vehicle mass), u = 1000 N (engine force), x (position of the vehicle in m), v (velocity of the vehicle in m/s) Objective: Once you determine the modeling equations of this control system, find the response of the car velocity when the input (u) is 1000 N. Use MATLAB to determine the following control characteristics of this cruise control system: 1. Transfer Function 2. Determine Poles: Are the poles stable? 3. Step Response: At what amplitude does it settle to? 4. Root-Locus: Explain the behavior. Calculated: State-Space variable form DKC-44
16 Extra Credit: If you were to vary the mass of the automobile from 1220 kg to either 500 kg or 2000 kg, what effect does this do to the response of the cruise control system? Solution: Model a basic cruise control system Figure 1: Simple Cruise Control System. %// Given values m = 1220; % mass of the vehicle b = 55; % damper constant u = 1000; % Input (engine force) % Setup the state-space matrices A = [0 1; 0 -b/m]; B = [0; 1/m]; C = [0 1]; D = [0]; DKC-45
17 % 1. Calculate the Transfer Function sys = ss(a,b*u,c,d); sys1 = tf(sys); (Answer) Transfer function: s % 2. Determine the poles and zeros and plot them num = [0.8197]; % zeros den = [ ]; % poles [z p k] = tf2zp(num,den); figure; pzmap(p,z); axis([ ]); (Answer) From the Pole-Zero plot, we can determine that the poles are stable since they are in the left-half plane of the s-plane. % 3. Do a step response of the system figure; sys = ss(a,b*u,c,d); sys1 = tf(sys); step(sys1); DKC-46
18 (Answer) Based on the step response plot, the system settles out at a magnitude value of K = % 4. Plot the Root Locus figure; rlocus(num,den); axis([ ]); (Answer) Looking at the root locus plot, we can tell that the system is definitely stable with one pole at p = With an increase in gain (K) the pole moves to the left in the s-plane towards infinity. DKC-47
ET3-7: Modelling I(V) Introduction and Objectives. Electrical, Mechanical and Thermal Systems
ET3-7: Modelling I(V) Introduction and Objectives Electrical, Mechanical and Thermal Systems Objectives analyse and model basic linear dynamic systems -Electrical -Mechanical -Thermal Recognise the analogies
More information2 Solving Ordinary Differential Equations Using MATLAB
Penn State Erie, The Behrend College School of Engineering E E 383 Signals and Control Lab Spring 2008 Lab 3 System Responses January 31, 2008 Due: February 7, 2008 Number of Lab Periods: 1 1 Objective
More information1.1 Write down the acceleration vector, in polar coordinates, in terms of time derivatives of r and
EN4: Dynamics and Vibrations Homework 3: Solving equations of motion for particles Solutions School of Engineering Brown University 1. The figure shows a simple idealization of a type of MEMS gyroscope
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationEXAMPLE: MODELING THE PT326 PROCESS TRAINER
CHAPTER 1 By Radu Muresan University of Guelph Page 1 EXAMPLE: MODELING THE PT326 PROCESS TRAINER The PT326 apparatus models common industrial situations in which temperature control is required in the
More informationUNIVERSIDADE FEDERAL DO AMAZONAS (UFAM) MATERIAL DE AULA PRÁTICA 01 MODELAGEM E SIMULAÇÃO MATEMÁTICA DE SISTEMAS FÍSICOS
UNIVERSIDADE FEDERAL DO AMAZONAS (UFAM) IDENTIFICAÇÃO DA DISCIPLINA Professor Pierre Vilar Dantas Disciplina Laboratório de Sistemas de Controle Horário Sextas, 16h00min às 18h00min Data 17/08/2018 Turma(s)
More informationChapter 5. The Laws of Motion
Chapter 5 The Laws of Motion The Laws of Motion The description of an object in There was no consideration of what might influence that motion. Two main factors need to be addressed to answer questions
More informationPhysics 1A, Week 2 Quiz Solutions
Vector _ A points north and vector _ B points east. If _ C = _ B _ A, then vector _C points: a. north of east. b. south of east. c. north of west. d. south of west. Find the resultant of the following
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
PH 105 Exam 2 VERSION A Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) Is it possible for a system to have negative potential energy? A)
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
PH 105 Exam 2 VERSION B Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) A boy throws a rock with an initial velocity of 2.15 m/s at 30.0 above
More informationChapter 5. The Laws of Motion
Chapter 5 The Laws of Motion The Laws of Motion The description of an object in motion included its position, velocity, and acceleration. There was no consideration of what might influence that motion.
More informationChapter 5. The Laws of Motion
Chapter 5 The Laws of Motion Sir Isaac Newton 1642 1727 Formulated basic laws of mechanics Discovered Law of Universal Gravitation Invented form of calculus Many observations dealing with light and optics
More informationMITOCW free_body_diagrams
MITOCW free_body_diagrams This is a bungee jumper at the bottom of his trajectory. This is a pack of dogs pulling a sled. And this is a golf ball about to be struck. All of these scenarios can be represented
More informationChapter 5 Force and Motion
Chapter 5 Force and Motion Chapter Goal: To establish a connection between force and motion. Slide 5-2 Chapter 5 Preview Slide 5-3 Chapter 5 Preview Slide 5-4 Chapter 5 Preview Slide 5-5 Chapter 5 Preview
More informationOrdinary Differential Equations (ODE)
Ordinary Differential Equations (ODE) Why study Differential Equations? Many physical phenomena are best formulated mathematically in terms of their rate of change. Motion of a swinging pendulum Bungee-jumper
More informationPhysics 1401V October 28, 2016 Prof. James Kakalios Quiz No. 2
This is a closed book, closed notes, quiz. Only simple (non-programmable, nongraphing) calculators are permitted. Define all symbols and justify all mathematical expressions used. Make sure to state all
More informationChapter 4. The Laws of Motion
Chapter 4 The Laws of Motion Classical Mechanics Describes the relationship between the motion of objects in our everyday world and the forces acting on them Conditions when Classical Mechanics does not
More informationSECTION 1.2. DYNAMIC MODELS
CHAPTER 1 BY RADU MURESAN Page 1 ENGG4420 LECTURE 5 September 16 10 6:47 PM SECTION 1.2. DYNAMIC MODELS A dynamic model is a mathematical description of the process to be controlled. Specifically, a set
More informationExam 2 Fall 2013
95.141 Exam 2 Fall 2013 Section number Section instructor Last/First name Last 3 Digits of Student ID Number: Answer all questions, beginning each new question in the space provided. Show all work. Show
More informationIn this activity, we explore the application of differential equations to the real world as applied to projectile motion.
Applications of Calculus: Projectile Motion ID: XXXX Name Class In this activity, we explore the application of differential equations to the real world as applied to projectile motion. Open the file CalcActXX_Projectile_Motion_EN.tns
More informationAnnouncements. Principle of Work and Energy - Sections Engr222 Spring 2004 Chapter Test Wednesday
Announcements Test Wednesday Closed book 3 page sheet sheet (on web) Calculator Chap 12.6-10, 13.1-6 Principle of Work and Energy - Sections 14.1-3 Today s Objectives: Students will be able to: a) Calculate
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
PH105-007 Exam 2 VERSION A Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) A 1.0-kg block and a 2.0-kg block are pressed together on a horizontal
More informationPhysics 11 Chapter 10: Energy and Work
Physics 11 Chapter 10: Energy and Work It is good to have an end to journey toward; but it is the journey that matters, in the end. Ursula K. Le Guin Nobody made a greater mistake than he who did nothing
More informationi. Indicate on the figure the point P at which the maximum speed of the car is attained. ii. Calculate the value vmax of this maximum speed.
1. A 0.20 kg object moves along a straight line. The net force acting on the object varies with the object's displacement as shown in the graph above. The object starts from rest at displacement x = 0
More informationAppendix 3B MATLAB Functions for Modeling and Time-domain analysis
Appendix 3B MATLAB Functions for Modeling and Time-domain analysis MATLAB control system Toolbox contain the following functions for the time-domain response step impulse initial lsim gensig damp ltiview
More informationWhat are Numerical Methods? (1/3)
What are Numerical Methods? (1/3) Numerical methods are techniques by which mathematical problems are formulated so that they can be solved by arithmetic and logic operations Because computers excel at
More informationClosed-form Method to Evaluate Bike Braking Performance
Human Power ejournal, April 4, 13 Closed-form Method to Evaluate Bike Braking Performance Junghsen Lieh, PhD Professor, Mechanical & Materials Engineering Wright State University, Dayton Ohio 45435 USA
More informationAmanda Stowers. Identifiers: Quantity Identifier Type Value Density of links Constant 1 kg/m 3. Angle to in direction Variable Initially 0
MIPSI: Effect of changing link lengths on a double pulum Question: How does changing the length of a single link in a double pulum affect the time it takes for the pulum s trajectory to change from its
More informationObject Impact on the Free Surface and Added Mass Effect Laboratory Fall 2005 Prof. A. Techet
Object Impact on the Free Surface and Added Mass Effect.016 Laboratory Fall 005 Prof. A. Techet Introduction to Free Surface Impact Free surface impact of objects has applications to ocean engineering
More informationFoundations of Physical Science. Unit One: Forces and Motion
Foundations of Physical Science Unit One: Forces and Motion Chapter 3: Forces and Motion 3.1 Force, Mass and Acceleration 3.2 Weight, Gravity and Friction 3.3 Equilibrium, Action and Reaction Learning
More informationChapter 5 Force and Motion
Force F Chapter 5 Force and Motion is the interaction between objects is a vector causes acceleration Net force: vector sum of all the forces on an object. v v N v v v v v Ftotal Fnet = Fi = F1 + F2 +
More informationEQUATIONS OF MOTION: RECTANGULAR COORDINATES
EQUATIONS OF MOTION: RECTANGULAR COORDINATES Today s Objectives: Students will be able to: 1. Apply Newton s second law to determine forces and accelerations for particles in rectilinear motion. In-Class
More informationThe roots are found with the following two statements. We have denoted the polynomial as p1, and the roots as roots_ p1.
Part II Lesson 10 Numerical Analysis Finding roots of a polynomial In MATLAB, a polynomial is expressed as a row vector of the form [an an 1 a2 a1 a0]. The elements ai of this vector are the coefficients
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationFALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Thursday, 11 December 2014, 6 PM to 9 PM, Field House Gym
FALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Thursday, 11 December 2014, 6 PM to 9 PM, Field House Gym NAME: STUDENT ID: INSTRUCTION 1. This exam booklet has 13 pages. Make sure none are missing 2.
More informationPotential energy functions used in Chapter 7
Potential energy functions used in Chapter 7 CHAPTER 7 CONSERVATION OF ENERGY Conservation of mechanical energy Conservation of total energy of a system Examples Origin of friction Gravitational potential
More informationPHYSICS. Chapter 5 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 5 Lecture RANDALL D. KNIGHT Chapter 5 Force and Motion IN THIS CHAPTER, you will learn about the connection between force and motion.
More informationChapter 6: Work and Kinetic Energy
Chapter 6: Work and Kinetic Energy Suppose you want to find the final velocity of an object being acted on by a variable force. Newton s 2 nd law gives the differential equation (for 1D motion) dv dt =
More informationPRELAB IMPULSE AND MOMENTUM
Impulse Momentum and Jump PRELAB IMPULSE AND MOMENTUM. In a car collision, the driver s body must change speed from a high value to zero. This is true whether or not an airbag is used, so why use an airbag?
More informationDynamics 4600:203 Homework 09 Due: April 04, 2008 Name:
Dynamics 4600:03 Homework 09 Due: April 04, 008 Name: Please denote your answers clearly, i.e., box in, star, etc., and write neatly. There are no points for small, messy, unreadable work... please use
More informationIntroduction to MATLAB
Introduction to MATLAB Violeta Ivanova, Ph.D. Educational Technology Consultant MIT Academic Computing violeta@mit.edu http://web.mit.edu/violeta/www/iap2006 Topics MATLAB Interface and Basics Linear Algebra
More informationCHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN
CHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN 6. Introduction Frequency Response This chapter will begin with the state space form of the equations of motion. We will use Laplace transforms to
More information6.1 Simulation of Simple Systems
LESSON 6.1: SIMULATION OF SIMPLE SYSTEMS 273 6.1 Simulation of Simple Systems In the previous lessons of Part II we looked at many simplified, yet realistic examples of dynamic systems. We have learned
More informationComputer lab for MAN460
Computer lab for MAN460 (version 20th April 2006, corrected 20 May) Prerequisites Matlab is rather user friendly, and to do the first exercises, it is enough to write an m-file consisting of two lines,
More informationPhysics 101. Hour Exam I Fall Last Name: First Name Network-ID Discussion Section: Discussion TA Name:
Last Name: First Name Network-ID Discussion Section: Discussion TA Name: Instructions Turn off your cell phone and put it away. Keep your calculator on your own desk. Calculators cannot be shared. This
More informationFall 2009 Sample Problems Exam 2
Sample Problems Exam 2 1. (24 points) In the table below is shown three physical situations in which two objects move and interact. The assumptions you are to make about the objects and the system are
More informationTime Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors
Time Response of Dynamic Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Multi-dimensional trajectories Numerical integration Linear and nonlinear systems Linearization
More informationCEE 271: Applied Mechanics II, Dynamics Lecture 17: Ch.15, Sec.2 4
1 / 38 CEE 271: Applied Mechanics II, Dynamics Lecture 17: Ch.15, Sec.2 4 Prof. Albert S. Kim Civil and Environmental Engineering, University of Hawaii at Manoa Tuesday, October 16, 2012 2 / 38 PRINCIPLE
More informationChapter 6. Circular Motion and Other Applications of Newton s Laws
Chapter 6 Circular Motion and Other Applications of Newton s Laws Circular Motion Two analysis models using Newton s Laws of Motion have been developed. The models have been applied to linear motion. Newton
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More information*ANSWER KEY * ANSWER KEY* ANSWER KEY* Newton's First Law of Motion Study Guide
*ANSWER KEY * ANSWER KEY* ANSWER KEY* Newton's First Law of Motion Study Guide Newton's First Law of Motion (Law of Inertia) An object at rest will remain at rest unless acted on by an unbalanced force.
More informationProf. Rupak Mahapatra. Physics 218, Chapter 7 & 8 1
Chapter 7, 8 & 9 Work and Eergy Prof. Rupak Mahapatra Physics 218, Chapter 7 & 8 1 Checklist for Today EOC Exercises from Chap 7 due on Monday Reading of Ch 8 due on Monday Physics 218, Chapter 7 & 8 2
More informationSimulink Modeling Tutorial
Simulink Modeling Tutorial Train system Free body diagram and Newton's law Model Construction Running the Model Obtaining MATLAB Model In Simulink, it is very straightforward to represent a physical system
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Spring Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs)
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs) is
More informationMath 56 Homework 1 Michael Downs. ne n 10 + ne n (1)
. Problem (a) Yes. The following equation: ne n + ne n () holds for all n R but, since we re only concerned with the asymptotic behavior as n, let us only consider n >. Dividing both sides by n( + ne n
More informationMULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. 1) A baseball is thrown vertically upward and feels no air resistance. As it is rising A) both
More informationHealy/DiMurro. Vibrations 2016
Name Vibrations 2016 Healy/DiMurro 1. In the diagram below, an ideal pendulum released from point A swings freely through point B. 4. As the pendulum swings freely from A to B as shown in the diagram to
More informationSOLUTION T 1 + U 1-2 = T C(31.5)(2.5)A10 6 B(0.2)D = 1 2 (7)(v 2) 2. v 2 = 2121 m>s = 2.12 km>s. Ans. (approx.
4 5. When a 7-kg projectile is fired from a cannon barrel that has a length of 2 m, the explosive force exerted on the projectile, while it is in the barrel, varies in the manner shown. Determine the approximate
More informationInaugural University of Michigan Science Olympiad Invitational Tournament. Hovercraft
Inaugural University of Michigan Science Olympiad Invitational Tournament Test length: 50 Minutes Hovercraft Team number: Team name: Student names: Instructions: Do not open this test until told to do
More informationRead textbook CHAPTER 1.4, Apps B&D
Lecture 2 Read textbook CHAPTER 1.4, Apps B&D Today: Derive EOMs & Linearization undamental equation of motion for mass-springdamper system (1DO). Linear and nonlinear system. Examples of derivation of
More informationCEE 271: Applied Mechanics II, Dynamics Lecture 9: Ch.13, Sec.4-5
1 / 40 CEE 271: Applied Mechanics II, Dynamics Lecture 9: Ch.13, Sec.4-5 Prof. Albert S. Kim Civil and Environmental Engineering, University of Hawaii at Manoa 2 / 40 EQUATIONS OF MOTION:RECTANGULAR COORDINATES
More informationGyroscopes and statics
Gyroscopes and statics Announcements: Welcome back from Spring Break! CAPA due Friday at 10pm We will finish Chapter 11 in H+R on angular momentum and start Chapter 12 on stability. Friday we will begin
More informationDifferential Equations and the Parachute Problem
Differential Equations and the Parachute Problem Ronald Phoebus and Cole Reilly May 10, 2004 Abstract The parachute problem is a classical first semester differential equations problem often introduced
More informationDr. Ian R. Manchester
Dr Ian R. Manchester Week Content Notes 1 Introduction 2 Frequency Domain Modelling 3 Transient Performance and the s-plane 4 Block Diagrams 5 Feedback System Characteristics Assign 1 Due 6 Root Locus
More informationThe dimensions of an object tend to change when forces are
L A B 8 STRETCHING A SPRING Hooke s Law The dimensions of an object tend to change when forces are applied to the object. For example, when opposite forces are applied to both ends of a spring, the spring
More informationGeneral Physics I Spring Applying Newton s Laws
General Physics I Spring 2011 Applying Newton s Laws 1 Equilibrium An object is in equilibrium if the net force acting on it is zero. According to Newton s first law, such an object will remain at rest
More informationPhysics I (Navitas) EXAM #2 Spring 2015
95.141 Physics I (Navitas) EXAM #2 Spring 2015 Name, Last Name First Name Student Identification Number: Write your name at the top of each page in the space provided. Answer all questions, beginning each
More informationANSWER'SHEET' 'STAPLE'TO'FRONT'OF'EXAM'! Name:!!!CWID:!!! Lab'section'(circle'one):' 6!(W!3pm)! 8!(W!7pm)!!!
ANSWER'SHEET' 'STAPLE'TO'FRONT'OF'EXAM' Name: CWID: Lab'section'(circle'one):' 6(W3pm) 8(W7pm) Multiplechoice: 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. Shortanswer: 16. 17. 18. 19. 20. 5(R7pm)
More informationANSWER'SHEET' 'STAPLE'TO'FRONT'OF'EXAM'! Name:!!!CWID:!!! Lab'section'(circle'one):' 6!(W!3pm)! 8!(W!7pm)!!!
ANSWER'SHEET' 'STAPLE'TO'FRONT'OF'EXAM' Name: CWID: Lab'section'(circle'one):' 6(W3pm) 8(W7pm) Multiplechoice: 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. Shortanswer: 16. 17. 18. 19. 20. 5(R7pm)
More informationExtra Credit Final Practice
Velocity (m/s) Elements of Physics I Spring 218 Extra Credit Final Practice 1. Use the figure to answer the following questions. Which of the lines indicate a) Constant velocity b) Constant acceleration
More informationAMS 27L LAB #6 Winter 2009
AMS 27L LAB #6 Winter 2009 Symbolically Solving Differential Equations Objectives: 1. To learn about the MATLAB Symbolic Solver 2. To expand knowledge of solutions to Diff-EQs 1 Symbolically Solving Differential
More informationSection /07/2013. PHY131H1F University of Toronto Class 9 Preclass Video by Jason Harlow. Based on Knight 3 rd edition Ch. 5, pgs.
PHY131H1F University of Toronto Class 9 Preclass Video by Jason Harlow Based on Knight 3 rd edition Ch. 5, pgs. 116-133 Section 5.1 A force is a push or a pull What is a force? What is a force? A force
More information2. To study circular motion, two students use the hand-held device shown above, which consists of a rod on which a spring scale is attached.
1. A ball of mass M attached to a string of length L moves in a circle in a vertical plane as shown above. At the top of the circular path, the tension in the string is twice the weight of the ball. At
More informationNEWTON S LAWS OF MOTION
NAME SCHOOL INDEX NUMBER DATE NEWTON S LAWS OF MOTION 1. 1995 Q21 P1 State Newton s first law of motion (1 mark) 2. 1998 Q22 P1 A body of mass M is allowed to slide down an inclined plane. State two factors
More informationChapters 5-6. Dynamics: Forces and Newton s Laws of Motion. Applications
Chapters 5-6 Dynamics: orces and Newton s Laws of Motion. Applications That is, describing why objects move orces Newton s 1 st Law Newton s 2 nd Law Newton s 3 rd Law Examples of orces: Weight, Normal,
More informationForces. Prof. Yury Kolomensky Feb 9/12, 2007
Forces Prof. Yury Kolomensky Feb 9/12, 2007 - Hooke s law - String tension - Gravity and Weight - Normal force - Friction - Drag -Review of Newton s laws Today s Plan Catalog common forces around us What
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 informationLaboratory 1. Solving differential equations with nonzero initial conditions
Laboratory 1 Solving differential equations with nonzero initial conditions 1. Purpose of the exercise: - learning symbolic and numerical methods of differential equations solving with MATLAB - using Simulink
More informationChapter 4 Forces Newton s Laws of Motion
Chapter 4 Forces Newton s Laws of Motion Forces Force A vector quantity that changes the velocity vector of an object. When you hit a baseball, the velocity of the ball changes. Can be a push or a pull
More informationPHYSICS 220 Lecture 04 Forces and Motion in 1 D Textbook Sections
PHYSICS 220 Lecture 04 Forces and Motion in 1 D Textbook Sections 3.2 3.6 Lecture 4 Purdue University, Physics 220 1 Last Lecture Constant Acceleration x = x 0 + v 0 t + ½ at 2 v = v 0 + at Overview v
More informationMidterm α, Physics 1P21/1P91
Midterm α, Physics 1P21/1P91 Prof. D. Crandles March 1, 2013 Last Name First Name Student ID Circle your course number above No examination aids other than those specified on this examination script are
More informationEE/ME/AE324: Dynamical Systems. Chapter 4: Block Diagrams
EE/ME/AE324: Dynamical Systems Chapter 4: Block Diagrams and Computer Simulation Block Diagrams A block diagram is an interconnection of: Blocks representing math operations Wires representing signals
More informationPHYSICS I RESOURCE SHEET
PHYSICS I RESOURCE SHEET Cautions and Notes Kinematic Equations These are to be used in regions with constant acceleration only You must keep regions with different accelerations separate (for example,
More informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
More informationMain Ideas in Class Today
Main Ideas in Class Today After today s class, you should be able to: Understand what a system is and determine whether it is isolated Apply Conservation of Linear Momentum to determine a final velocity
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE201: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE21: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationPhysics 125, Spring 2006 Monday, May 15, 8:00-10:30am, Old Chem 116. R01 Mon. 12:50 R02 Wed. 12:50 R03 Mon. 3:50. Final Exam
Monday, May 15, 8:00-10:30am, Old Chem 116 Name: Recitation section (circle one) R01 Mon. 12:50 R02 Wed. 12:50 R03 Mon. 3:50 Closed book. No notes allowed. Any calculators are permitted. There are no trick
More informationDynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras
Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 09 Characteristics of Single Degree - of -
More informationPHYSICS. Chapter 5 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 5 Lecture RANDALL D. KNIGHT Chapter 5 Force and Motion IN THIS CHAPTER, you will learn about the connection between force and motion.
More informationUnit 06 Examples. Stuff you asked about:
Unit 06 Examples Today s Concepts: 1. Force due to gravity 2. Force due to strings 3. Force due to springs (just a little bit) Mechanics Lecture 5, Slide 1 Stuff you asked about: Can you go over all of
More informationSolutions PHYS 250 Exam 2 Practice Test
Solutions PHYS 250 Exam 2 Practice Test 1C This questions really deals with conservation of energy, ΔK +ΔU = 0. The main problem is determining the initial height, h, of the man just as he starts to swing.
More informationEF 151 Final Exam - Spring, 2016 Page 1 Copy 1
EF 151 Final Exam - Spring, 016 Page 1 Copy 1 Name: Section: Instructions: Sit in assigned seat; failure to sit in assigned seat results in a 0 for the exam. Put name and section on your exam. Put seating
More informationDynamics: Forces and Newton s Laws of Motion
Lecture 7 Chapter 5 Physics I Dynamics: Forces and Newton s Laws of Motion Course website: http://faculty.uml.edu/andriy_danylov/teaching/physicsi Today we are going to discuss: Chapter 5: Force, Mass:
More information9. Introduction and Chapter Objectives
Real Analog - Circuits 1 Chapter 9: Introduction to State Variable Models 9. Introduction and Chapter Objectives In our analysis approach of dynamic systems so far, we have defined variables which describe
More informationNon-textbook problem #I: The kinetic energy of a body depends on its mass and speed as. K = 1 2 mv2 (1) m 1 v 2 1 = 1 2 m 2v 2 2 (2)
PHY 309 K. Solutions for Problem set # 7. Non-textbook problem #I: The kinetic energy of a body depends on its mass and speed as K = 1 mv (1) Therefore, two bodies of respective masses m 1 and m and speeds
More informationPHYSICS GUIDESHEET UNIT 5. - ENERGY SUBUNIT - ENERGY CONVERSIONS POTENTIAL AND KINETIC ENERGY ACTIVITY LESSON DESCRIPTION SCORE/POINTS
1 NAME PERIOD PHYSICS GUIDESHEET UNIT 5. - ENERGY SUBUNIT - ENERGY CONVERSIONS POTENTIAL AND KINETIC ENERGY ACTIVITY LESSON DESCRIPTION SCORE/POINTS 1. NT CLASS OVERHEAD NOTES (5 pts/page) /20 (Plus 5
More informationChapter 8 Solutions. The change in potential energy as it moves from A to B is. The change in potential energy in going from A to B is
Chapter 8 Solutions *8. (a) With our choice for the zero level for potential energy at point B, U B = 0. At point A, the potential energy is given by U A = mgy where y is the vertical height above zero
More information