Scribd
Upload
58 views2 pages

KDM Matlab Assignment

The document contains MATLAB code for calculating angles and lengths in a four-bar linkage mechanism. It computes the follower angle, coupler angle, and transmission angle based on given crank, coupler, and follower lengths. Additionally, it includes a section for calculating the coupler angle and slider length based on user input for the crank length and input angle.

Uploaded by

mayankmaroliya
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
58 views2 pages

KDM Matlab Assignment

The document contains MATLAB code for calculating angles and lengths in a four-bar linkage mechanism. It computes the follower angle, coupler angle, and transmission angle based on given crank, coupler, and follower lengths. Additionally, it includes a section for calculating the coupler angle and slider length based on user input for the crank length and input angle.

Uploaded by

mayankmaroliya
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd

clc;clear;close all;

%a = input ('Enter crank length : ');


%b = input ('Enter coupler length : ');
%c = input ('Enter follower length : ');
%d = input ('Enter fixed length : ');

a=20;
b=40;
c=45;
d=50;

k1 = d/a;
k2 = d/c;
k3 = (a^2-b^2+c^2+d^2)/(2*a*c);
k4 = d/b;
k5 = (c^2-d^2-a^2-b^2)/(2*a*b);

th3_o = zeros(13,1);
th3_c = zeros(13,1);
th4_o = zeros(13,1);
th4_c = zeros(13,1);
tr_o = zeros(13,1);
tr_c = zeros(13,1);
th2 = 0;

for i=1:13

%calculating follower angle (th4)


A = cosd(th2)-k1-(k2*cosd(th2))+k3;
B = -2*sind(th2);
C = k1-(k2+1)*(cosd(th2))+k3;
th4_o(i,1)=2*atand((-B - sqrt(B^2-(4*A*C)))/(2*A));
th4_c(i,1)=2*atand((-B + sqrt(B^2-(4*A*C)))/(2*A)) + 360;

%calculating coupler angle (th3)


D = cosd(th2)-k1+(k4*cosd(th2))+k5;
E = -2*sind(th2);
F = k1 + (k4-1)*cosd(th2) + k5;
th3_o(i,1)=2*atand((-E - sqrt(E^2-(4*D*F)))/(2*D));
th3_c(i,1)=2*atand((-E + sqrt(E^2-(4*D*F)))/(2*D)) + 360;

%calculating Transmission angle (tr)


tr_o(i,1) = acosd((b^2+c^2-a^2-d^2)/(2*b*c))+((a*d*cosd(th2)/(b*c)));
tr_c(i,1) = acosd((b^2+c^2-a^2-d^2)/(2*b*c))-((a*d*cosd(th2)/(b*c)));
%tr_c(i,1) = acosd((c^2+b^2-d^2-a^2)/(2*b*c))-((d*a/b*c)*cosd(th2));

th2=th2+30;
end
disp('Follower Angle');
disp([th4_o,th4_c]);
disp('Coupler Angle');
disp([th3_o,th3_c]);
disp('Transmission Angle');
%disp(tr_o);
disp([tr_o,tr_c]);

..................................................................................

clc;clear;close all;
a = input ('Enter crank length : ');
b = input ('Enter coupler length : ');
c = input ('Enter slider length : ');
th2 = input ('Enter Input angle : ');

th4 = 90;

%calculating coupler angle (th3)


th3 = asind((c-(a*sind(th2)))/b);
th3 = 360-th3;
%calculating the length of slider (d)
d = (a*cosd(th2)) + (b*cosd(th3));

disp('Coupler angle : ');


disp(th3);
disp('Length of slider : ');
disp(d);

You might also like