In this solution some basic concepts and formulas of Mechanisms Kinematics are used. For more information, refer to any standard textbook or drop a comment below. Please give a Thumbs Up, if solution is helpful.
Solution :
The matlab code to create the plots is follows,
r = 5/12; % Crank radius (= r in
ft)
l = 30/12; % Length of Connecting Rod in ft
n = l/r;
w = 20; % Angular Velocity
t = 0:0.1:(pi); % Theta angle
b = asin(r*sin(t)/l); % Beta angle
m = 5+10; % Total mass
g = 32.3; % Acc. due to gravity
Fg = m*g; % Gravity force
Fi = m*r*(w^2)*(cos(t) + (cos(t)/n));
F = Fi+Fg; % Total force of piston
Fc = F./cos(b); % Thurst along connecting rod
Fn = Fc.*sin(b); % Thrust on cylinder sides
Ft = Fc./sin(t+b); % Crank effort
Fr = Fc.*cos(t+b); % Thrust on bearings
aD = r*(w^2)*(cos(t) + (cos(t)/n)); %
Acceleration of D
aC = -(w^2).*sin(t).*((n^2 -1 )./((n^2 - (sin(t)).^2)).^(3/2)); %
Angulr acceleration of C
subplot(1,2,1)
plot(t,F)
title('Various Forces')
ylabel('Force (in Lbs)')
xlabel('Theta angle, t (in radians)')
grid on
hold on
plot(t,Fc,'cyan')
plot(t,Fn,'black')
plot(t,Ft,'red')
plot(t,Fr,'green')
legend('F','Fc','Fn','Ft','Fr')
hold off
subplot(1,2,2)
plot(t,aD)
title('Various Accelertions')
ylabel('Acceleration in ft/s^2 or rad/s^2')
xlabel('Theta angle, t (in radians)')
grid on
hold on
plot(t,aC,'red')
legend('aD','aC')
NOTE : Feel free to ask further queries. Your positive rating would be appreciated and motivate me !
1 of 2 MAE 371 Spring 2019 Project For the slider crank configuration, shown below, derive the expression for the velocity and acceleration of the piston, A (taken positive to the right), as a f...