Determine the rotation, that is, the slope, of the beam ABC shown in Fig. 4.14(a) at A.
Repeat using the Symbolic Math Toolbox in MATLAB.
Determine the rotation, that is, the slope, of the beam ABC shown in Fig. 4.14(a) at A.
Repeat using the Symbolic Math Toolbox in MATLAB.
The actual rotation of the beam at A produced by the actual concentrated load, W, is \theta_{ A }. Let us suppose that a virtual unit moment is applied at A before the actual rotation takes place, as shown in Fig. 4.14(b). The virtual unit moment induces virtual support reactions of R_{v, A }(=1 / L) acting downward and R_{v, C }(=1 / L) acting upward. The actual internal bending moments are
M_{ A }=+\frac{W}{2} x, \quad 0 \leq x \leq L / 2M_{ A }=+\frac{W}{2}(L-x), \quad L / 2 \leq x \leq L
The internal virtual bending moment is
M_{v}=1-\frac{1}{L} x, \quad 0 \leq x \leq L
The external virtual work done is 1 \theta_{ A } (the virtual support reactions do no work, as there is no vertical displacement of the beam at the supports), and the internal virtual work done is given by Eq. (4.21). Hence,
W_{i, M}=\sum \int_{L} \frac{M_{ A } M_{v}}{E I} d x (4.21)
1 \theta_{ A }=\frac{1}{E I}\left[\int_{0}^{L / 2} \frac{W}{2} x\left(1-\frac{x}{L}\right) d x+\int_{L / 2}^{L} \frac{W}{2}(L-x)\left(1-\frac{x}{L}\right) d x\right] (i)
Simplifying Eq. (i), we have
\theta_{ A }=\frac{W}{2 E I L}\left[\int_{0}^{L / 2}\left(L x-x^{2}\right) d x+\int_{L / 2}^{L}(L-x)^{2} d x\right] (ii)
Hence,
\theta_{ A }=\frac{W}{2 E I L}\left\{\left[L \frac{x^{2}}{2}-\frac{x^{3}}{3}\right]_{0}^{L / 2}-\frac{1}{3}\left[(L-x)^{3}\right]_{L / 2}^{L}\right\}
from which
\theta_{ A }=\frac{W L^{2}}{16 E I}
Calculation of the slope of the beam shown in Fig. 4.14(a) at A is obtained through the following MATLAB file:
% Declare any needed variables
syms M_A M_V W x L theta_A E I
% Define the moments due to the applied load (M_A) and the unit virtual load (M_V)
M_A = [W/2*x; % 0 >= x >= L/2
W/2*(L-x)]; % L/2 >= x >= L
M_V = 1 – x/L; % 0 >= x >= L
% Define equations for the external (W_e) and internal (W_i) virtual work
W_e = 1*theta_A;
W_i = (int(M_A(1)*M_V,x,0,L/2) + int(M_A(2)*M_V,x,L/2,L))/(E*I); % From Eq. (4.21)
% Equate W_e and W_i, and solve for the slope of the beam (theta_A)
theta_A = solve(W_e-W_i,theta_A);
% Output theta_A to the Command Window
disp([‘theta_A =’ char(theta_A)])
The Command Window output resulting from this MATLAB file is as follows:
theta_A=(L^2*W)/(16*E*I)