Simulation, Transfer Function
|
|
- Ilene Chase
- 6 years ago
- Views:
Transcription
1 Max Force (lb) Displacement (in) 1 ME313 Homework #12 Simulation, Transfer Function Last Updated November 17, 214. Repeat the car-crash problem from HW#6. Use the Matlab function lsim with ABCD format to perform the simulation. Note, that, for the car crash problem, the input u=. Before the simulation, compute the eigenvalues of [A], and the time constants and frequencies of oscillation. Set the time range for the simulation to be 3 times the longest time constant. Again consider the carcrash problem. Only now, write a Matlab script that will perform a simulation for a sequence of seat-belt stiffnesses k 2 ranging from k 2 =5 lb/ft to k 2 =5 klb/ft and seatbelt damping coefficients c 2 ranging from c 2 =2 lb s/ft to c 2 =2 lb s/ft. Use the function lsim to perform each simulation. After each simulation, extract the maximum force exerted by the seatbelt on the occupant, and the maximum relative displacement between the occupant and car. Plot both of these quantities versus the seat-belt stiffness k 2 and seatbelt damping ratio c 2 as a surface Your plots should look like those shown above Damping Coefficient (lbs/ft) Damping Coefficient (lbs/ft) Maximum Occupant Displacement Maximum Force on Occupant Stiffness (klb/ft 2 3 Stiffness (klb/ft
2 2 Esfandiari and Lu text, problem set 4.3, problems 1 and 2. Only compute the transfer function for the specified output. I.e., for problem 1, the specified output is x 2 (t). And, o for problem 1, determine T(j ), T(j ), arg[t(j )], and the steady-state output x 2 (t) if f(t)=5sin(2t). o For problem 2, determine T(j ), T(j ), arg[t(j )], and the steady-state output x 1 (t) if f(t)=.5sin(1.5t). Esfandiari and Lu text, problem set 4.3, problem 7. Defining q 1 (t) as the output o Determine the transfer function T(s) between Q 1 (s) as output and V(s) as input. o Plot the amplitude and phase response of q 1 (t) assuming that v(t)=3sin( t), and the driving frequency varies over the range < <1 rad/sec. o Use the amplitude and phase plot to determine the steady state response v(t) when =.5 rad/sec. It will be useful to you to use the command grid on; with your plot.
3 Displacement (m) 3 Handout problem Solve the problem by plotting the frequency response of the vertical motion of the car versus frequency, and identify the natural frequencies by locating the frequencies at which peaks in amplitude of car motion appear. Hint Suppose the road surface vertical variation is y(x) =.25sin ( 2π x) ft, where x (in ft) is the horizontal displacement of the vehicle along the road. This would model the road variation with a period of 2 ft as described in the problem statement. Then suppose the vehicle is travelling at a constant velocity v. Then, since v=x/t, the road surface variation as a function of time will be y(t) =.25sin ( 2π vt). 2 Again handout problem Use lsim and tf to simulate the motion of the car when driven over a speed bump. Perform the simulation over a time range of t 1 sec. Simulate the bump as an input of y(t)=.1m of t 1 sec, and y(t)=m of 1 t 1 sec. Repeat the simulation when a damper c=2 lb s/ft is added between the wheel and car (you will need to derive a new transfer function to include the damper). Plot both simulations on the same graph. Your result should look like: Car Without Damper Car With Damper Bump Time (sec)
4 Code for Car Crash Simulation clear all;close all;clc; HW 7 Linear Simulation Problem ICS x1= Vehicle position x2= Vehicle velocity x3= Occupant position (inertial) x4= Occupant velocity (inertial) xo(1)=; xo(2)=-(528/36)*5; xo(3)=; xo(4)=-(528/36)*5; k1=1e3; c1=5; m1=3/32.2; k2=5e3; c2=5 ; m2=15/32.2; A=[ 1 ; -(k1+k2)/m1 -(c1+c2)/m1 k2/m1 c2/m1 ; 1 ; k2/m2 c2/m2 -k2/m2 -c2/m2]; B=[ ; ; ; ]; C=[1-1 ]; D=; [E,L]=eig(A) wn1=abs(l(1,1));z1=-real(l(1,1))/wn1;wd1=wn1*sqrt(1-z1^2);td1=(2*pi)/wd1; tau1=1/(z1*wn1); wn2=abs(l(3,3));z2=-real(l(3,3))/wn2;wd2=wn2*sqrt(1-z2^2);td2=(2*pi)/wd2; tau2=1/(z2*wn2); Set up time range tend=5*max([tau1 tau2]); N=5; t=linspace(,tend,n); dt=t(2)-t(1); Perform simulation u=zeros(1,length(t)); sys=ss(a,b,c,d); [y,t,x]=lsim(sys,u,t,xo); y=x(:,1)-x(:,3); f=(-k2*(x(:,3)-x(:,1))-c2*(x(:,4)-x(:,2))); subplot(3,1,1); plot(t,x(:,1)*12,t,x(:,3)*12) title('inertial Coordinates of Car and Occupant') axis([ tend -1*12 1*12]); legend('x1','x3'); ylabel('displacement (in)');
5 Force (lb) Displacement (in) Displacement (in) subplot(3,1,2); plot(t,y*12) title('relative Displacement Between Car and Occupant') axis([ tend 1.1*min(y)*12 1.1*max(y)*12 ]); ylabel('displacement (in)'); subplot(3,1,3); plot(t,f) title('force of Seat-Belt on Occupant') axis([ tend 1.1*min(f) 1.1*max(f)]); ylabel('force (lb)'); xlabel('time'); Result of Simulation 1 Inertial Coordinates of Car and Occupant x1 x Relative Displacement Between Car and Occupant Force of Seat-Belt on Occupant Time
6 Maximum force and relative displacement for varying seatbelt stiffness and damping ratio. Code below to generate data shown in plots contained in assignment (although plot axis limits are included in the code). clear all; close all; Compute maximum force exerted by seatbelt on occupant for a range of seatbelt stiffnesses k1=1e3; c1=5; m1=3/32.2; k2=5e3; c2=5 ; m2=15/32.2; Nk2=2; k2t=linspace(5,5e3,nk2); Nc2=2; c2t=linspace(2,2,nc2); ICS x1= Vehicle position x2= Vehicle velocity x3= Occupant position (inertial) x4= Occupant velocity (inertial) xo(1)=; xo(2)=-(528/36)*5; xo(3)=; xo(4)=-(528/36)*5; for i1=1:nk2 for j1=1:nc2 k1=1e3; c1=5; m1=3/32.2; k2=5e3; c2=5 ; m2=15/32.2; k2=k2t(i1); c2=c2t(j1); A=[ 1 ; -(k1+k2)/m1 -(c1+c2)/m1 k2/m1 c2/m1 ; 1 ; k2/m2 c2/m2 -k2/m2 -c2/m2]; B=[ ; ; ; ]; C=[1-1 ]; D=; [E,L]=eig(A) wn1=abs(l(1,1));z1=-real(l(1,1))/wn1; tau1=1/(z1*wn1); wn2=abs(l(3,3));z2=-real(l(3,3))/wn2; tau2=1/(z2*wn2); tend=5*max([tau1 tau2]); N=5; t=linspace(,tend,n); u=zeros(1,length(t)); sys=ss(a,b,c,d); [y,t,x]=lsim(sys,u,t,xo); y1=x(:,1)-x(:,3); ymax(j1,i1)=max(abs(y1));
7 f=(-k2*(x(:,3)-x(:,1))-c2*(x(:,4)-x(:,2))); fmax(j1,i1)=max(abs(f)); end end subplot(2,1,1); surf(k2t/1,c2t,ymax*12) title('maximum Occupant Displacement') xlabel('stiffness (klb/ft'); ylabel('damping Coefficient (lbs/ft)'); zlabel('displacement (in)'); subplot(2,1,2); surf(k2t/1,c2t,fmax) title('maximum Force on Occupant') xlabel('stiffness (klb/ft'); ylabel('damping Coefficient (lbs/ft)'); zlabel('max Force (lb)');
8
9 Code for Esfandiari and Lu text, set 4.3, problem 7 and accompanying instructions V=3; Nw=5; w=linspace(,1,nw); for i1=1:nw s=1i*w(i1); T(i1)=(2*s+4)/(2*s^3+6*s^2+9*s+2); end subplot(2,1,1); plot(w,v*abs(t)) ylabel('amplitude Q_1'); xlabel('frequency (rad/sec)'); grid on; subplot(2,1,2); plot(w,angle(t)*(18/pi)) ylabel('phase {\phi} (deg)'); xlabel('frequency (rad/sec)'); grid on; Resulting Output 6 Amplitude Q Frequency (rad/sec) Phase (deg) Frequency (rad/sec) At =5 rad/sec, the amplitude of q1(t) is approximately 3, and the phase is approximately -65 o. So, the steady state solution q 1(t) is q 1(t)=3sin(.5t-65o).
10
11
12 Car Displacement Amplitude x3 (ft) Car Displacement Amplitude x3 (ft) Code for frequency response clear all;close all;clc ms=625/32.2; mu=19/32.2; ks=1e3; kt=8e3; bs=; Nv=5; v=linspace(1,1,nv); w=((2*pi)/2)*v; for i1=1:nv s=i*w(i1); N=bs*kt*s+ks*kt; D=ms*mu*s^4+(bs*ms+bs*mu)*s^3+(ks*ms+kt*ms+ks*mu)*s^2+bs*kt*s+ks*kt; T(i1)=N/D; end plot(v,abs(t)*.15) ylabel('car Displacement Amplitude x3 (ft)'); xlabel('velocity (ft/sec)'); title(['b_s=',num2str(bs)]); axis([ 1 1]); Frequency response for bs= and bs=4 1 b s = Velocity (ft/sec) 1 b s = Velocity (ft/sec)
13 Code for simulation of car over speed bump clear all; close all; clc Car Suspension Frequency Response ms=625/32.2; mu=19/32.2; ks=1e3; kt=8e3; bs=4; b=[bs*kt ks*kt]; a=[ms*mu bs*ms+bs*mu ks*ms+kt*ms+ks*mu bs*kt ks*kt]; TFsys=tf(b,a); Nt=5; t=linspace(,1,nt); x1=zeros(1,nt); for i1=1:nt if (t(i1) <= 1) x1(i1)=.1; else x1(i1)=; end end x3(1)=; x3(2)=; [x3,t]=lsim(tfsys,x1,t,x3); plot(t,x3,t,x1) ylabel('car Displacement x3 (ft)'); title(['shock Absorber b_s=',num2str(bs),' lbs/ft']); xlabel('time (sec)');
14 Car Displacement x3 (ft) Car Displacement x3 (ft).25 Shock Absorber b s = lbs/ft Time (sec).2 Shock Absorber b s =4 lbs/ft Time (sec)
Analysis of Linear State Space Models
Max f (lb) Max y (in) 1 ME313 Homework #7 Analysis of Linear State Space Models Last Updated November 6, 213 Repeat the car-crash problem from HW#6. Use the Matlab function lsim to perform the simulation.
More informationState Space Representation
ME Homework #6 State Space Representation Last Updated September 6 6. From the homework problems on the following pages 5. 5. 5.6 5.7. 5.6 Chapter 5 Homework Problems 5.6. Simulation of Linear and Nonlinear
More informationLesson 14: Van der Pol Circuit and ode23s
Lesson 4: Van der Pol Circuit and ode3s 4. Applied Problem. A series LRC circuit when coupled via mutual inductance with a triode circuit can generate a sequence of pulsing currents that have very rapid
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 informationProblem 1 Matlab code
Problem 1 Matlab code x=[.5377 1.8339-2.2588.8622.3188-1.377 -.4336.3426 3.5784 2.7694-1.3499 3.349.7254 -.631.7147 -.25]; X=fft(x) fs=1;n=16;dt=1/fs;t=n*dt; f=[ 1/T 2/T 3/T 4/T 5/T 6/T 7/T 8/T]' N=16
More informationME 3143 System Dynamics & Modeling
ME 3143 System Dynamics & Modeling Project II Car Suspension System Analysis Using Bond-Graph Approach Group # 1 Group members Cooperation (%) Lauren Baxter 14.2857/100 Michael Bryant 14.2857/100 Ross
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
R. M. Murray Fall 2004 CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 101/110 Homework Set #2 Issued: 4 Oct 04 Due: 11 Oct 04 Note: In the upper left hand corner of the first page
More informationProb. 1 SDOF Structure subjected to Ground Shaking
Prob. 1 SDOF Structure subjected to Ground Shaking What is the maximum relative displacement and the amplitude of the total displacement of a SDOF structure subjected to ground shaking? magnitude of ground
More informationMath Lab 9 23 March 2017 unid:
Math 2250-004 Lab 9 23 March 2017 unid: Name: Instructions and due date: Due: 30 March 2017 at the start of lab. If extra paper is necessary, please staple it at the end of the packet. For full credit:
More informationApplication 7.4 Earthquake-Induced Vibrations of Multistory Buildings
Application 7.4 Earthquake-Induced Vibrations of Multistory Buildings In this application you are to investigate the response to transverse earthquake ground oscillations of a seven-story building like
More informationParticle Motion. Typically, if a particle is moving along the x-axis at any time, t, x()
Typically, if a particle is moving along the x-axis at any time, t, x() t represents the position of the particle; along the y-axis, yt () is often used; along another straight line, st () is often used.
More informationMatlab Sheet 3 - Solution. Plotting in Matlab
f(x) Faculty of Engineering Spring 217 Mechanical Engineering Department Matlab Sheet 3 - Solution Plotting in Matlab 1. a. Estimate the roots of the following equation by plotting the equation. x 3 3x
More informationSection 3.7: Mechanical and Electrical Vibrations
Section 3.7: Mechanical and Electrical Vibrations Second order linear equations with constant coefficients serve as mathematical models for mechanical and electrical oscillations. For example, the motion
More informationFeedback Control Systems
ME Homework #0 Feedback Control Systems Last Updated November 06 Text problem 67 (Revised Chapter 6 Homework Problems- attached) 65 Chapter 6 Homework Problems 65 Transient Response of a Second Order Model
More informationCDS 101: Lecture 2.1 System Modeling
CDS 101: Lecture 2.1 System Modeling Richard M. Murray 4 October 2004 Goals: Define what a model is and its use in answering questions about a system Introduce the concepts of state, dynamics, inputs and
More informationa. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules.
Lab #1 - Free Vibration Name: Date: Section / Group: Procedure Steps (from lab manual): a. Follow the Start-Up Procedure in the laboratory manual. Note the safety rules. b. Locate the various springs and
More informationCh 3.7: Mechanical & Electrical Vibrations
Ch 3.7: Mechanical & Electrical Vibrations Two important areas of application for second order linear equations with constant coefficients are in modeling mechanical and electrical oscillations. We will
More informationParticle Motion. Typically, if a particle is moving along the x-axis at any time, t, x()
Typically, if a particle is moving along the x-axis at any time, t, x() t represents the position of the particle; along the y-axis, yt () is often used; along another straight line, st () is often used.
More informationChapter 14 Oscillations. Copyright 2009 Pearson Education, Inc.
Chapter 14 Oscillations 14-1 Oscillations of a Spring If an object vibrates or oscillates back and forth over the same path, each cycle taking the same amount of time, the motion is called periodic. The
More informationChapter 14 Oscillations
Chapter 14 Oscillations If an object vibrates or oscillates back and forth over the same path, each cycle taking the same amount of time, the motion is called periodic. The mass and spring system is a
More information1st Simulation: Convert the quarter-car model to state-space form and run it using "lsim" built-in function
Lab 3 Chris Rini 7/16/2015 Contents Prepare Matlab to handle this script (Wipe the slate) Define variables and assumptions Set initial conditions (ICs) Establish the road profile using the old roadgenerator2
More informationMath 308 Exam II Practice Problems
Math 38 Exam II Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More information( ) for t 0. Rectilinear motion CW. ( ) = t sin t ( Calculator)
Rectilinear motion CW 1997 ( Calculator) 1) A particle moves along the x-axis so that its velocity at any time t is given by v(t) = 3t 2 2t 1. The position x(t) is 5 for t = 2. a) Write a polynomial expression
More informationMATH2351 Introduction to Ordinary Differential Equations, Fall Hints to Week 07 Worksheet: Mechanical Vibrations
MATH351 Introduction to Ordinary Differential Equations, Fall 11-1 Hints to Week 7 Worksheet: Mechanical Vibrations 1. (Demonstration) ( 3.8, page 3, Q. 5) A mass weighing lb streches a spring by 6 in.
More informationThe... of a particle is defined as its change in position in some time interval.
Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle
More informationLesson 11: Mass-Spring, Resonance and ode45
Lesson 11: Mass-Spring, Resonance and ode45 11.1 Applied Problem. Trucks and cars have springs and shock absorbers to make a comfortable and safe ride. Without good shock absorbers, the truck or car will
More informationSecond Order Systems
Second Order Systems independent energy storage elements => Resonance: inertance & capacitance trade energy, kinetic to potential Example: Automobile Suspension x z vertical motions suspension spring shock
More information6.003 Homework #6. Problems. Due at the beginning of recitation on October 19, 2011.
6.003 Homework #6 Due at the beginning of recitation on October 19, 2011. Problems 1. Maximum gain For each of the following systems, find the frequency ω m for which the magnitude of the gain is greatest.
More informationMath Lab 9 3 November 2016 unid:
Math 2250-004 Lab 9 3 November 2016 unid: Name: Instructions and due date: Due: 10 November 2016 at the start of lab. If extra paper is necessary, please staple it at the end of the packet. For full credit:
More informationCE 513: STATISTICAL METHODS
28-8-217/CE 68 CE 513: STATISTICAL METHODS IN CIVIL ENGINEERING Lecture: Introduction to Fourier transforms Dr. Budhaditya Hazra Room: N-37 Department of Civil Engineering 1 Fourier Analysis Fourier Series
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 informationLab Experiment 2: Performance of First order and second order systems
Lab Experiment 2: Performance of First order and second order systems Objective: The objective of this exercise will be to study the performance characteristics of first and second order systems using
More informationUnit 5 PreCalculus Review
Class: Date: Unit 5 PreCalculus Review Multiple Choice Identify the choice that best completes the statement or answers the question. 1. Find the terminal point P (x, y) on the unit circle determined by
More informationMATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam
MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html
More information2.003 Engineering Dynamics Problem Set 10 with answer to the concept questions
.003 Engineering Dynamics Problem Set 10 with answer to the concept questions Problem 1 Figure 1. Cart with a slender rod A slender rod of length l (m) and mass m (0.5kg)is attached by a frictionless pivot
More informationCHAPTER 12 TIME DOMAIN: MODAL STATE SPACE FORM
CHAPTER 1 TIME DOMAIN: MODAL STATE SPACE FORM 1.1 Introduction In Chapter 7 we derived the equations of motion in modal form for the system in Figure 1.1. In this chapter we will convert the modal form
More informationLAB #6 The Swaying Building
LAB #6 The Swaying Building Goal: Determine a model of the swaying of a skyscraper; estimating parameters Required tools: Matlab routines pplane, ode45, plot; M-files; systems of differential equations.
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 informationAP Calculus. Particle Motion. Student Handout
AP Calculus Particle Motion Student Handout 016-017 EDITION Use the following link or scan the QR code to complete the evaluation for the Study Session https://www.surveymonkey.com/r/s_sss Copyright 016
More informationMotion, Speed, Velocity & Acceleration. Physical Science Bella Vista Middle School
Motion, Speed, Velocity & Acceleration Physical Science Bella Vista Middle School What Is Motion? Motion is when an object changes place or position. To properly describe motion, you need to use the following:
More informationHOMEWORK ANSWERS. Lesson 4.1: Simple Harmonic Motion
DEVIL PHYSICS HOMEWORK ANSWERS Tsokos, Chapter 3 Test Lesson 4.1: Simple Harmonic Motion 1. Objectives. By the end of this class you should be able to: a) Understand that in simple harmonic motion there
More informationSession 2: MDOF systems
BRUFACE Vibrations and Acoustics MA Academic year 7-8 Cédric Dumoulin (cedumoul@ulb.ac.be) Arnaud Deraemaeker (aderaema@ulb.ac.be) Session : MDOF systems Exercise : Multiple DOFs System Consider the following
More informationChapter 14 Oscillations. Copyright 2009 Pearson Education, Inc.
Chapter 14 Oscillations Oscillations of a Spring Simple Harmonic Motion Energy in the Simple Harmonic Oscillator Simple Harmonic Motion Related to Uniform Circular Motion The Simple Pendulum The Physical
More informationa) Find the equation of motion of the system and write it in matrix form.
.003 Engineering Dynamics Problem Set Problem : Torsional Oscillator Two disks of radius r and r and mass m and m are mounted in series with steel shafts. The shaft between the base and m has length L
More informationTo convert a speed to a velocity. V = Velocity in feet per seconds (ft/sec) S = Speed in miles per hour (mph) = Mathematical Constant
To convert a speed to a velocity V S ( 1.466) V Velocity in feet per seconds (ft/sec) S Speed in miles per hour (mph) 1.466 Mathematical Constant Example Your driver just had a rear-end accident and says
More informationHonors Algebra 2 Chapter 14 Page 1
Section. (Introduction) Graphs of Trig Functions Objectives:. To graph basic trig functions using t-bar method. A. Sine and Cosecant. y = sinθ y y y y 0 --- --- 80 --- --- 30 0 0 300 5 35 5 35 60 50 0
More informationCar Dynamics using Quarter Model and Passive Suspension; Part V: Frequency Response Considering Driver-seat
357 Car Dynamics using Quarter Model and Passive Suspension; Part V: Frequency Response Considering Driver-seat Galal Ali Hassaan Emeritus Professor, Department of Mechanical Design & Production, Faculty
More informationSection , #5. Let Q be the amount of salt in oz in the tank. The scenario can be modeled by a differential equation.
Section.3.5.3, #5. Let Q be the amount of salt in oz in the tank. The scenario can be modeled by a differential equation dq = 1 4 (1 + sin(t) ) + Q, Q(0) = 50. (1) 100 (a) The differential equation given
More informationSection 6.1 Angles and Radian Measure Review If you measured the distance around a circle in terms of its radius, what is the unit of measure?
Section 6.1 Angles and Radian Measure Review If you measured the distance around a circle in terms of its radius, what is the unit of measure? In relationship to a circle, if I go half way around the edge
More informationME 563 HOMEWORK # 5 SOLUTIONS Fall 2010
ME 563 HOMEWORK # 5 SOLUTIONS Fall 2010 PROBLEM 1: You are given the lumped parameter dynamic differential equations of motion for a two degree-offreedom model of an automobile suspension system for small
More informationLectures Chapter 10 (Cutnell & Johnson, Physics 7 th edition)
PH 201-4A spring 2007 Simple Harmonic Motion Lectures 24-25 Chapter 10 (Cutnell & Johnson, Physics 7 th edition) 1 The Ideal Spring Springs are objects that exhibit elastic behavior. It will return back
More informationChapter 4: The State Space and Numerical Simulation
Chapter 4: The State Space and Numerical Simulation Samantha Ramirez Preview Questions What mathematical formulation might we apply to the models we derived in the previous chapter in order to facilitate
More informationFigure XP3.1 (a) Mass in equilibrium, (b) Freebody diagram, (c) Kinematic constraint relation Example Problem 3.1 Figure XP3.1 illustrates a mass m
LECTURE 7. MORE VIBRATIONS ` Figure XP3.1 (a) Mass in equilibrium, (b) Freebody diagram, (c) Kinematic constraint relation Example Problem 3.1 Figure XP3.1 illustrates a mass m that is in equilibrium and
More informationEE 4343/ Control System Design Project LECTURE 10
Copyright S. Ikenaga 998 All rights reserved EE 4343/5329 - Control System Design Project LECTURE EE 4343/5329 Homepage EE 4343/5329 Course Outline Design of Phase-lead and Phase-lag compensators using
More informationA-level Physics (7407/7408)
A-level Physics (7407/7408) Further Mechanics Test Name: Class: Date: September 2016 Time: 55 Marks: 47 Page 1 Q1.The diagram shows a strobe photograph of a mark on a trolley X, moving from right to left,
More informationECE382/ME482 Spring 2005 Homework 7 Solution April 17, K(s + 0.2) s 2 (s + 2)(s + 5) G(s) =
ECE382/ME482 Spring 25 Homework 7 Solution April 17, 25 1 Solution to HW7 AP9.5 We are given a system with open loop transfer function G(s) = K(s +.2) s 2 (s + 2)(s + 5) (1) and unity negative feedback.
More informationProject 2: Using linear systems for numerical solution of boundary value problems
LINEAR ALGEBRA, MATH 124 Instructor: Dr. T.I. Lakoba Project 2: Using linear systems for numerical solution of boundary value problems Goal Introduce one of the most important applications of Linear Algebra
More informationCDS 101: Lecture 2.1 System Modeling
CDS 0: Lecture. System Modeling Richard M. Murray 4 Octoer 004 Goals: Define what a model is and its use in answering questions aout a system Introduce the concepts of state, dynamics, inputs and outputs
More informationSession 1 : Fundamental concepts
BRUFACE Vibrations and Acoustics MA1 Academic year 17-18 Cédric Dumoulin (cedumoul@ulb.ac.be) Arnaud Deraemaeker (aderaema@ulb.ac.be) Exercise 1 Session 1 : Fundamental concepts Consider the following
More informationEx. 1. Find the general solution for each of the following differential equations:
MATH 261.007 Instr. K. Ciesielski Spring 2010 NAME (print): SAMPLE TEST # 2 Solve the following exercises. Show your work. (No credit will be given for an answer with no supporting work shown.) Ex. 1.
More informationME 328 Machine Design Vibration handout (vibrations is not covered in text)
ME 38 Machine Design Vibration handout (vibrations is not covered in text) The ollowing are two good textbooks or vibrations (any edition). There are numerous other texts o equal quality. M. L. James,
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 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 informationSchool of Mechanical Engineering Purdue University
Case Study ME375 Frequency Response - 1 Case Study SUPPORT POWER WIRE DROPPERS Electric train derives power through a pantograph, which contacts the power wire, which is suspended from a catenary. During
More informationx 3x 1 if x 3 On problems 8 9, use the definition of continuity to find the values of k and/or m that will make the function continuous everywhere.
CALCULUS AB WORKSHEET ON CONTINUITY AND INTERMEDIATE VALUE THEOREM Work the following on notebook paper. On problems 1 4, sketch the graph of a function f that satisfies the stated conditions. 1. f has
More informatione jωt = cos(ωt) + jsin(ωt),
This chapter introduces you to the most useful mechanical oscillator model, a mass-spring system with a single degree of freedom. Basic understanding of this system is the gateway to the understanding
More informationLecture 6 mechanical system modeling equivalent mass gears
M2794.25 Mechanical System Analysis 기계시스템해석 lecture 6,7,8 Dongjun Lee ( 이동준 ) Department of Mechanical & Aerospace Engineering Seoul National University Dongjun Lee Lecture 6 mechanical system modeling
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 informationContents. Dynamics and control of mechanical systems. Focus on
Dynamics and control of mechanical systems Date Day 1 (01/08) Day 2 (03/08) Day 3 (05/08) Day 4 (07/08) Day 5 (09/08) Day 6 (11/08) Content Review of the basics of mechanics. Kinematics of rigid bodies
More informationFinal Examination Thursday May Please initial the statement below to show that you have read it
EN40: Dynamics and Vibrations Final Examination Thursday May 0 010 Division of Engineering rown University NME: General Instructions No collaboration of any kind is permitted on this examination. You may
More informationMAE 143B - Homework 7
MAE 143B - Homework 7 6.7 Multiplying the first ODE by m u and subtracting the product of the second ODE with m s, we get m s m u (ẍ s ẍ i ) + m u b s (ẋ s ẋ u ) + m u k s (x s x u ) + m s b s (ẋ s ẋ u
More informationChapter 14 Preview Looking Ahead
Chapter 14 Preview Looking Ahead Text: p. 438 Slide 14-1 Chapter 14 Preview Looking Back: Springs and Restoring Forces In Chapter 8, you learned that a stretched spring exerts a restoring force proportional
More informationEECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total
More informationWheel and Axle. Author: Joseph Harrison. Research Ans Aerospace Engineering 1 Expert, Monash University
Wheel and Axle Author: Joseph Harrison British Middle-East Center for studies & Research info@bmcsr.com http:// bmcsr.com Research Ans Aerospace Engineering 1 Expert, Monash University Introduction A solid
More informationChapter 5. Section 5.1. Section ( x ) ( y ) 7. ( x ) ( y ) (0, 3 + 5) and (0, 3 5)
9 Chapter Section.. 0. ( x ) ( y ). ( x 7 ) + ( y+ ) = 9 7. ( x ) ( y ) 8 + + 0 = 8 + 8 = 9.. (0, + ) and (0, ). (.60786, 7.6887). (-.07,.8) 7. 9.87 miles Section. 70 0 -.. 00. 0 7. 9.. 8 9.. miles 7.
More informationEnergy and the Simple Harmonic Oscillator
OpenStax-CNX module: m42244 Energy and the Simple Harmonic Oscillator OpenStax College his wor is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 3.0 Abstract Determine
More informationInternational Journal of Multidisciplinary and Current Research
International Journal of Multidisciplinary and Current Research Research Article ISSN: 2321-3124 Available at: http://ijmcr.com Theoretical and Numerical Analysis of Half Car Vehicle Dynamic Model Subjected
More informationThe area under the velocity/time curve is equal to the total change in displacement
Mousetrap.car.notes.problems Topics that will be studied with mousetrap cars are: motion in one dimension under constant acceleration torque and its relationship to angular and linear acceleration angular
More informationAPPLICATIONS FOR ROBOTICS
Version: 1 CONTROL APPLICATIONS FOR ROBOTICS TEX d: Feb. 17, 214 PREVIEW We show that the transfer function and conditions of stability for linear systems can be studied using Laplace transforms. Table
More informationMA FINAL EXAM Green May 5, You must use a #2 pencil on the mark sense sheet (answer sheet).
MA 600 FINAL EXAM Green May 5, 06 NAME STUDENT ID # YOUR TA S NAME RECITATION TIME. You must use a # pencil on the mark sense sheet (answer sheet).. Be sure the paper you are looking at right now is GREEN!
More informationC. points X and Y only. D. points O, X and Y only. (Total 1 mark)
Grade 11 Physics -- Homework 16 -- Answers on a separate sheet of paper, please 1. A cart, connected to two identical springs, is oscillating with simple harmonic motion between two points X and Y that
More informationAircraft Dynamics First order and Second order system
Aircraft Dynamics First order and Second order system Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Aircraft dynamic
More information4-5 Graphing Other Trigonometric Functions
Locate the vertical asymptotes, and sketch the graph of each function. 1. y = 2 tan x 4. y = 3 tan 2. 5. 3. 6. y = tan 3x esolutions Manual - Powered by Cognero Page 1 7. y = 2 tan (6x π) 10. 8. 11. y
More informationPHY 3221 Fall Homework Problems. Instructor: Yoonseok Lee. Submit only HW s. EX s are additional problems that I encourage you to work on.
PHY 3221 Fall 2012 Homework Problems Instructor: Yoonseok Lee Submit only HW s. EX s are additional problems that I encourage you to work on. Week 1: August 22-24, Due August 27 (nothing to submit) EX:
More information11. Some applications of second order differential
October 3, 2011 11-1 11. Some applications of second order differential equations The first application we consider is the motion of a mass on a spring. Consider an object of mass m on a spring suspended
More informationFuzzy Logic Control for Half Car Suspension System Using Matlab
Fuzzy Logic Control for Half Car Suspension System Using Matlab Mirji Sairaj Gururaj 1, Arockia Selvakumar A 2 1,2 School of Mechanical and Building Sciences, VIT Chennai, Tamilnadu, India Abstract- To
More informationECE2111 Signals and Systems UMD, Spring 2013 Experiment 1: Representation and manipulation of basic signals in MATLAB
ECE2111 Signals and Systems UMD, Spring 2013 Experiment 1: Representation and manipulation of basic signals in MATLAB MATLAB is a tool for doing numerical computations with matrices and vectors. It can
More informationChapter 5 Design. D. J. Inman 1/51 Mechanical Engineering at Virginia Tech
Chapter 5 Design Acceptable vibration levels (ISO) Vibration isolation Vibration absorbers Effects of damping in absorbers Optimization Viscoelastic damping treatments Critical Speeds Design for vibration
More information(a) Torsional spring-mass system. (b) Spring element.
m v s T s v a (a) T a (b) T a FIGURE 2.1 (a) Torsional spring-mass system. (b) Spring element. by ky Wall friction, b Mass M k y M y r(t) Force r(t) (a) (b) FIGURE 2.2 (a) Spring-mass-damper system. (b)
More informationM A : Ordinary Differential Equations
M A 2 0 5 1: Ordinary Differential Equations Essential Class Notes & Graphics D 19 * 2018-2019 Sections D07 D11 & D14 1 1. INTRODUCTION CLASS 1 ODE: Course s Overarching Functions An introduction to the
More informationMATH 100 Introduction to the Profession
MATH 100 Introduction to the Profession Differential Equations in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2012 fasshauer@iit.edu MATH 100 ITP 1 What
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 informationAnswers and Hints to Review Questions for Test (a) Find the general solution to the linear system of differential equations Y = 2 ± 3i.
Answers and Hints to Review Questions for Test 3 (a) Find the general solution to the linear system of differential equations [ dy 3 Y 3 [ (b) Find the specific solution that satisfies Y (0) = (c) What
More information6.003 Homework #6 Solutions
6.3 Homework #6 Solutions Problems. Maximum gain For each of the following systems, find the frequency ω m for which the magnitude of the gain is greatest. a. + s + s ω m = w This system has poles at s
More informationProblem Weight Score Total 100
EE 350 EXAM IV 15 December 2010 Last Name (Print): First Name (Print): ID number (Last 4 digits): Section: DO NOT TURN THIS PAGE UNTIL YOU ARE TOLD TO DO SO Problem Weight Score 1 25 2 25 3 25 4 25 Total
More informationHomework 2: Kinematics and Dynamics of Particles Due Friday Feb 12, Pressurized tank
EN40: Dynamics and Vibrations Homework 2: Kinematics and Dynamics of Particles Due Friday Feb 12, 2016 School of Engineering Brown University 1. The figure illustrates an idealized model of a gas gun (used,
More informationHomework Assignment 3
ECE382/ME482 Fall 2008 Homework 3 Solution October 20, 2008 1 Homework Assignment 3 Assigned September 30, 2008. Due in lecture October 7, 2008. Note that you must include all of your work to obtain full
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationEntering Matrices
MATLAB Basic Entering Matrices Transpose Subscripts The Colon Operator Operators Generating Matrices Element-by-element operation Creating a Plot Multiple Data Sets in One Graph Line Styles and Colors
More informationLecture 4&5 MATLAB applications In Signal Processing. Dr. Bedir Yousif
Lecture 4&5 MATLAB applications In Signal Processing Dr. Bedir Yousif Signal Analysis Fourier Series and Fourier Transform Trigonometric Fourier Series Where w0=2*pi/tp and the Fourier coefficients a n
More information