hw6_sol - MEEN 364 Fall 2004 Homework Set 6 Due October 15,...

Info iconThis preview shows pages 1–3. Sign up to view the full content.

View Full Document Right Arrow Icon
MEEN 364 Homework Set 6 Fall 2004 October 8, 2004 Homework Set 6 – Due October 15, 2004 @ 5:00 PM 1. (a) We have (1.1) 0 sin ) ( 2 = + + + θ & & & R g R l To obtain the equilibrium point, the derivative terms are equated to zero. ( ) 0 , 0 = = & & & 0 ) sin( 0 = g π n = 0 , n = 0, ± 1, ± 2 The only meaningful value, however is 0 0 = (b) Expanding each term in (1.1) by Taylor’s series about the equilibrium point and neglecting the higher order terms, we get () 0 0 ) 0 ( sin sin ) 0 ( ) 0 ( ) 0 ( 0 2 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 = + + + + + + + + = = = = = = = = = = = = = = = = = = & & & & & & & & & & & & & & & & & & & & & & R R g g R R R l l 0 0 0 0 0 0 0 = + + + + + + + g l & & 0 lg θθ ⇒+ && = 0 g l = (c) The simulation results for the linear and non-linear cases are given below: 0 1 2 3 4 5 6 7 8 9 10 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 time ( s ) ( t ) ( rad ) θ vs tim e Linear Non-linear 0 1 2 3 4 5 6 7 8 9 10 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 time ( s ) ω ω vs time Linear Non-linear 1
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
MEEN 364 Homework Set 6 Fall 2004 October 8, 2004 1) M-file ---------------------------------------------------------------------------------------- % single-pendulum function xp=sing(t,x) % linear state-equations for a single-pendulum. % x=[theta thetadot] % xp is the time derivative of the state-vector, x. R=0.2;l=1;g=9.8; x1=x(1); xp(1,1)=x(2); xp(2,1)=-g/l*x1; ---------------------------------------------------------------------------------------- function xp=sing(t,x) % nonlinear state-equations for a single-pendulum. R=0.2;l=1;g=9.8; x1=x(1); % the angle xp(1,1)=x(2); % angular derivative is the angular velocity xp(2,1)=-(g*sin(x1)+R*x(2)*x(2))/(l+R*x1); % Newton's 2nd law ---------------------------------------------------------------------------------------- % plot [t1,x1]=ode45(@Lsing,[0 10],[pi/6 0]'); %[0 10]: time span, [pi/6 0]: initial conditions [t2,x2]=ode45(@sing,[0 10],[pi/6 0]'); figure(1),plot(t1,x1(:,1),':',t2,x2(:,1)); xlabel('time ({\its})');ylabel('\theta({\itt}) ({\itrad})'); title('\theta vs time');legend('Linear','Non-linear');
Background image of page 2
Image of page 3
This is the end of the preview. Sign up to access the rest of the document.

This note was uploaded on 04/18/2010 for the course MEEN 364 taught by Professor Staff during the Spring '08 term at Texas A&M.

Page1 / 5

hw6_sol - MEEN 364 Fall 2004 Homework Set 6 Due October 15,...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online