This work studies scattering-induced elastic wave attenuation and phase velocity variation in three-dimensional untextured cubic polycrystals with statistically equiaxed grains using the theoretical second-order approximation (SOA) and Born approximation models and the grain-scale finite-element (FE) model, pushing the boundary towards strongly scattering materials. The results for materials with Zener anisotropy indices A > 1 show a good agreement between the theoretical and FE models in the transition and stochastic regions. In the Rayleigh regime, the agreement is reasonable for common structural materials with 1 < A <  3.2 but it deteriorates as A increases. The wavefields and signals from FE modelling show the emergence of very strong scattering at low frequencies for strongly scattering materials that cannot be fully accounted for by the theoretical models. To account for such strong scattering at A > 1, a semi-analytical model is proposed by iterating the far-field Born approximation and optimizing the iterative coefficient. The proposed model agrees remarkably well with the FE model across all studied materials with greatly differing microstructures; the model validity also extends to the quasi-static velocity limit. For polycrystals with A < 1, it is found that the agreement between the SOA and FE results is excellent for all studied materials and the correction of the model is not needed.