Improvement in accuracy of the finite element method in analysis of stability of Euler–Bernoulli and Timoshenko beams

The problem of buckling of the Euler–Bernoulli and Timoshenko beams is analyzed by the finite element method. Significant improvement in accuracy of the method is obtained by replacing the discontinuous function of the bending moment related to the approximation of the eigenfunction obtained by FEM by a “smoothed” function in the Rayleigh quotient. The smoothed function is obtained by fitting the discontinuous one using the least square technique.