This paper deals with modelling the effect of local fibre volume fraction variability, fibre misalignment and fibre strength variability on the longitudinal tensile strength of unidirectional plies with finite element analysis. Variability is accounted for by generating spatially-correlated fields of fibre misalignment and volume fraction. This information is then translated into local mechanical properties and orientations in finite element models of the ply, which are virtually tested in longitudinal tension. Monte Carlo simulations were performed to evaluate the effect of different sources of material variability, i.e. local fibre strength, fibre volume fraction and misalignment. Ply strength predictions lowered when including the variability of local volume fraction and fibre misalignment in the modelling, showing a better agreement with experiments for the carbon/epoxy system investigated.