Based on the treatment that the rectangular function can be expanded into an approximate sum of complex Gaussian functions with finite numbers, the analytical expression for the focusing intensity distribution of a partially coherent Gaussian Schell-model (GSM) beam through an aperture lens was derived. The typical numerical examples were given and were compared to those obtained from numerically integral calculation. The results have shown that the method can significantly improve the numerical calculation efficiency.