In this work analytical solutions are developed to study the free vibration of composite beams under axial loading. The beam with a single delamination is modeled as four interconnected Euler-Bernoulli beams using the delamination as their boundary. The continuity and the equilibrium conditions are satisfied between the adjoining beams. The studies show that the sizes and the locations of the delaminations significantly influence the natural frequencies and mode shapes of the beam. A monotonic relation between the natural frequency and the axial load is predicted.