Element Stresses 

for the beam element, shouldn't the derivative of phi have (12*x)/L^3 for V_1 and same for V_2? this would change the values when evaluating the derivate at x=L