Determine the vertical deflection of the free end of the cantilever beam shown in Fig. 4.13(a).
Repeat using the Symbolic Math Toolbox in MATLAB.
Determine the vertical deflection of the free end of the cantilever beam shown in Fig. 4.13(a).
Repeat using the Symbolic Math Toolbox in MATLAB.
Let us suppose that the actual deflection of the cantilever at B produced by the uniformly distributed load is V_{ B } and that a vertically downward virtual unit load is applied at B before the actual deflection takes place. The external virtual work done by the unit load is, from Fig. 4.13(b), 1 V_{ B } The deflection, V_{ B }, is assumed to be caused by bending only; that is, we ignore any deflections due to shear. The internal virtual work is given by Eq. (4.21), which, since only one member is involved, becomes
W_{i, M}=\sum \int_{L} \frac{M_{ A } M_{v}}{E I} d x (4.21)
W_{i, M}=\int_{0}^{L} \frac{M_{ A } M_{v}}{E I} d x (i)
The virtual moments, M_{v}, are produced by a unit load, so we replace M_{v} by M_{1}. Then,
W_{i, M}=\int_{0}^{L} \frac{M_{ A } M_{1}}{E I} d x (ii)
At any section of the beam a distance x from the built-in end,
M_{ A }=-\frac{w}{2}(L-x)^{2}, \quad M_{1}=-1(L-x)
Substituting for M_{ A } and M_{ 1 } in Eq. (ii) and equating the external virtual work done by the unit load to the internal virtual work, we have
1 v_{ B }=\int_{0}^{L} \frac{w}{2 E I}(L-x)^{3} d x
which gives
1 v_{ B }=-\frac{w}{2 E I}\left[\frac{1}{4}(L-x)^{4}\right]_{0}^{L}
so that
v_{ B }=\frac{w L^{4}}{8 E I}
Note that V_{ B } is in fact negative, but the positive sign here indicates that it is in the same direction as the unit load.
Calculation of the vertical deflection of the free end of the cantilever shown in Fig. 4.13(a) is obtained through the following MATLAB file:
% Declare any needed variables
syms M_A M_V W_e W_i E I w x L v_B
% Define the moments due to the applied load (M_A) and the unit virtual load (M_V)
M_A = -0.5*w*(L-x)^2;
M_V = -1*(L-x);
% Define equations for the external (W_e) and internal (W_i) virtual work
W_e = 1*v_B;
W_i = int(M_A*M_V/(E*I),x,0,L); % From Eq. (4.21)
% Equate W_e and W_i, and solve for the free end vertical displacement (v_B)
v_B = solve(W_e-W_i,v_B);
% Output v_B to the Command Window
disp([‘v_B =’ char(v_B)])
The Command Window output resulting from this MATLAB file is as follows:
v_B = (L^4*w)/(8*E*I)