<p>This paper presents a stochastic analysis of composite spherical shells (CSS) featuring cutouts, examining how the size of the cutouts, and the orientation angle of the plies influence the natural frequency (NF). Using finite element (FE) methods, we model laminated composite spherical shells with square cutouts to capture the first three stochastic NF, influenced by material and geometric properties variations. The analysis incorporates brute force Monte Carlo simulation to address the inherent randomness in these parameters. The first three natural frequencies is validated through FE analysis and mesh convergence studies. The research looks at the overall effects of the unknown NF on several input parameters, including thickness, ply orientation angle, Poisson’s ratio, modulus of elasticity, shear modulus, and mass density. The results are validate from previously published results, and found accuracy up to 97.7 %. Sensitivity analysis and highlights the influence of variability in these parameters, providing insights into optimizing design for improved performance and reliability, and found that twist angle, ply orientation angle and plate thickness has more impact as compare to other parameters. The findings offer valuable contributions to understanding the dynamic behaviour and reliability of composite spherical shells with cutouts under stochastic conditions.</p>