A method for detection of multiple open cracks in a slender Euler–Bernoulli beams is presented based on frequency measurements. The method is based on the approach given by Hu and Liang [J. Franklin Inst. 330 (5) (1993) 841], transverse vibration modelling through transfer matrix method and representation of a crack by rotational spring. The beam is virtually divided into a number of segments, which can be decided by the analyst, and each of them is considered to be associated with a damage parameter. The procedure gives a linear relationship explicitly between the changes in natural frequencies of the beam and the damage parameters. These parameters are determined from the knowledge of changes in the natural frequencies. After obtaining them, each is treated in turn to exactly pinpoint the crack location in the segment and determine its size. The forward, or natural frequency determination, problems are examined in the passing. The method is approximate, but it can handle segmented beams, any boundary conditions, intermediate spring or rigid supports, etc. It eliminates the need for any symbolic computation which is envisaged by Hu and Liang [J. Franklin Inst. 330 (5) (1993) 841] to obtain mode shapes of the corresponding uncracked beams. The proposed method gives a clear insight into the whole analysis. Case studies (numerical) are presented to demonstrate the method effectiveness for two simultaneous cracks of size 10% and more of section depth. The differences between the actual and predicted crack locations and sizes are less than 10% and 15% respectively. The numbers of segments into which the beam is virtually divided limits the maximum number of cracks that can be handled. The difference in the forward problem is less than 5%.