Hi Triambakeshwar Malladi,
It appears that you want to add animation for the coupler curve of crank slider mechanism. Here is a complete MATLAB code example for the same:
% Define the mechanism parameters
a = 5; % Length of the crank
b = 10; % Length of the connecting rod
e = 7; % Length of the coupler
e1 = 5; % Offset of the slider
ar = pi/6; % Angle of the coupler from the connecting rod
sx = 1; % Half-width of the slider
sy = 2; % Half-height of the slider
% Define the simulation parameters
numSteps = 100; % Number of simulation steps
T = linspace(0, 2*pi, numSteps); % Crank angle range from 0 to 2*pi
Bita = asin(a/b * sin(T)); % Angle of the connecting rod
% Initialize arrays to store coupler curve points
X = zeros(1, numSteps);
Y = zeros(1, numSteps);
% Calculate coupler curve points
for i = 1:numSteps
X(i) = a*cos(T(i)) + e*cos(ar - Bita(i));
Y(i) = a*sin(T(i)) + e*sin(ar - Bita(i));
end
% Create a figure for the animation
figure;
hold on;
title('Slider-Crank Mechanism with Coupler Curve');
xlim([-a-b-e, a+b+e]);
ylim([-a-b-e, a+b+e]);
xlabel('X');
ylabel('Y');
grid on;
% Animation loop
for i = 1:numSteps
% Coordinates of Slider Mechanism Link
Cx = [0, a*cos(T(i)), a*cos(T(i)) + b*cos(Bita(i)), X(i)];
Cy = [0, a*sin(T(i)), e1, Y(i)];
% Mechanism Plotting
crank = line([Cx(1), Cx(2)], [Cy(1), Cy(2)], 'color', 'r', 'linewidth', 5); % Crank
coupler = line([Cx(2), Cx(4)], [Cy(2), Cy(4)], 'color', 'g', 'linewidth', 5); % Coupler
connecting_rod = line([Cx(2), Cx(3)], [Cy(2), Cy(3)], 'color', 'b', 'linewidth', 5); % Connecting rod
% Slider representation
Slider1 = line([Cx(3)+sx, Cx(3)+sx], [Cy(3)+sy, Cy(3)-sy], 'color', 'k', 'linewidth', 1);
Slider2 = line([Cx(3)-sx, Cx(3)-sx], [Cy(3)+sy, Cy(3)-sy], 'color', 'k', 'linewidth', 1);
% Plotting the coupler curve
plot(X(1:i), Y(1:i), 'm', 'Linewidth', 1);
% Animation of Plotting by Deleting the previous plotting
pause(0.05);
if i < numSteps
delete(crank);
delete(coupler);
delete(connecting_rod);
delete(Slider1);
delete(Slider2);
end
end
% Display final position (if needed)
crank = line([Cx(1), Cx(2)], [Cy(1), Cy(2)], 'color', 'r', 'linewidth', 5); % Crank
coupler = line([Cx(2), Cx(4)], [Cy(2), Cy(4)], 'color', 'g', 'linewidth', 5); % Coupler
connecting_rod = line([Cx(2), Cx(3)], [Cy(2), Cy(3)], 'color', 'b', 'linewidth', 5); % Connecting rod
% Display crank radius and slide distance
Crank_Radius = a;
Slide_Distance = max(X) - min(X);
disp(['Crank Radius: ', num2str(Crank_Radius)]);
disp(['Slide Distance: ', num2str(Slide_Distance)]);
In this example missing parameters are assigned with some value for the demonstration.
Best Regards,
Abhishek Chakram
