Fatigue cracks nucleate at the surface of the shaft operating under cyclic load. These cracks then grow into semi-elliptical shapes while extending into the shaft. The paper presents numerical estimation of stress intensity factor (SIF) of such cracks in circular solid shaft when subjected to bending and twisting loads in elastic conditions. The finite element model is first validated with the help of results reported elsewhere for a single surface crack. Later on, multiple coplanar surface cracks are modeled and the effect of the number of cracks and the distance between them on SIF of the main crack is investigated. The results are in line with the expectations.Nomenclature