These are the tasks that I am trying to accomplish. What I do not understand is how to find a way to extract a certain value for u. Is it even possible to solve for a certain value of w?
I also don't understand how to extract values from my solved matrix "x" and plot them vs theta. If Nb and Torque are being solve in the loop, how is it possible to extract the data and plot them.
I'm NOT looking for someone to write the code. I just don't know where to start.
--> IN YOUR COMMENTS INDICATE The value for w (to nearest tenth of a rad/s) at which your box loses contact WHEN u = 0.8 and ALSO when u = 0.2.
--> PLOT as an output of your Code NB vs THETA and Drive Torque T vs THETA (two separate plots--figure 1 & figure 2). The plots should be produced when your program runs and be well labeled with units etc.
theta=30
sign=(cosd (theta)/abs (cosd (theta)))
u=.8
w=(2*pi/5)
R=1
c=4.05/12
H=18/12
y=H-(R*sind(theta/2))
Md=50/32.2
Id=(.5*Md*R^2)
Mp=75/32.2
Mb=30/32.2
alpha=0
u=.8
L=(6/12*R)
g=32.2
ax=(-R*w^2)*cosd(theta)
ab=(-R*w^2*cosd(theta))-(R*w^2*sind(theta))
ay=(-R*w^2)*sind(theta)
a=[1 0 1 0 0 0 0 0 0 0 0; 0 1 0 1 0 0 0 0 0 0 0; 0 0 -R*sind(theta) R*cosd(theta) 0 0 0 0 0 0 1; 0 0 0 0 1 0 1 0 0 0 0; 0 0 0 0 0 1 0 1 0 0 0; 0 0 0 0 -R*sind(theta) R*cosd(theta) 0 0 0 0 0; 0 0 -1 0 -1 0 0 0 (4*abs(sign)) 0 0; 0 0 0 -1 0 -1 0 0 (-4*u*sign) 0 0; 0 0 -c (3*R) -c (-3*R) 0 0 (-4*(y-c)) 0 0; 0 0 0 0 0 0 0 0 -sign u*sign 0; 0 0 0 0 0 0 0 0 (-u*sign) 1 0]
C=[0 Md*g Id*alpha 0 Md*g Id*alpha Mp*ax Mp*ay+Mp*g 0 Mb*ab Mb*g]
x=C/inv(a)
w = 0;
dw = 0.01;
rownum = 9;
colnum = 9;
while a(rownum,colnum) >= 0
w = w + dw;
end
