ixed beam by two supports and concrete load calculate shear force diagram and bending moment diagram and deflection in matlab

main.m
% Fixed beam with two supports and concrete load

% Define constants
L = 10; % Length of the beam (in meters)
w = 5; % Uniformly distributed load (in kN/m)
E = 200 * 10^9; % Young's Modulus of the beam material (in Pa)
I = 10^-4; % Moment of inertia of the beam (in m^4)

% Create the symbolic variable x
syms x;

% Define the functions for load, shear force, bending moment, and deflection
w_x = @(x) w; % Uniformly distributed load

V_x = @(x) w*x - w*x^2/(2*L); % Shear force function
M_x = @(x) w*x^2/2 - w*x^3/(6*L); % Bending moment function

M_max = max(abs(M_x(x))); % Maximum bending moment

delta_x = w*x^4 / (24*E*I) * (L^3 - 2*L*x^2 + x^3); % Deflection function

% Plot the shear force diagram
figure;
fplot(V_x, [0, L], 'LineWidth', 1.5);
title('Shear Force Diagram');
xlabel('Position along the beam (m)');
ylabel('Shear Force (kN)');
grid on;

% Plot the bending moment diagram
figure;
fplot(M_x, [0, L], 'LineWidth', 1.5);
title('Bending Moment Diagram');
xlabel('Position along the beam (m)');
ylabel('Bending Moment (kN.m)');
grid on;

% Display the maximum bending moment
disp(['Maximum Bending Moment: ', num2str(M_max), ' kN.m']);

% Calculate and display the deflection at the midpoint of the beam
midpoint = L/2;
deflection = double(subs(delta_x, x, midpoint));
disp(['Deflection at midpoint: ', num2str(deflection), ' meters']);
1354 chars
45 lines