calculate shear force diagram and bending moment diagram and deflection for beam fixed by two supports in matlab

main.m
% Define properties of the beam
L = 5; % Beam length in meters
E = 200e9; % Young's modulus in Pa
I = 2e-4; % Area moment of inertia in m^4
w = 10; % Distributed load in N/m

% Calculate reactions at the supports
syms R1 R2;
eq1 = R1 + R2 - w*L == 0;
eq2 = R1*L/2 - w*L^2/2 - R2*L == 0;
sol = solve([eq1, eq2], [R1, R2]);

R1 = sol.R1;
R2 = sol.R2;

% Shear force and bending moment equations
x = 0:0.01:L;
V = R1 - w*x;
M = R1*x - w*x.^2/2;

% Deflection equation (assuming a point load at the center of the beam)
syms x1;
delta = (w/(24*E*I)) * (x1^3 - 4*L*x1^2 + 6*L^2*x1);

% Plotting the diagrams and deflection curve
figure;
subplot(3,1,1);
plot(x, V, 'r', 'LineWidth', 2);
xlabel('Position (m)');
ylabel('Shear Force (N)');
title('Shear Force Diagram');

subplot(3,1,2);
plot(x, M, 'b', 'LineWidth', 2);
xlabel('Position (m)');
ylabel('Bending Moment (Nm)');
title('Bending Moment Diagram');

subplot(3,1,3);
x1_val = 0:0.01:L;
delta_val = subs(delta, x1, x1_val);
plot(x1_val, delta_val, 'g', 'LineWidth', 2);
xlabel('Position (m)');
ylabel('Deflection (m)');
title('Deflection Curve');
1095 chars
46 lines