We calculate the correlation function of 79,091 galaxy clusters in the redshift region of \(0.05 \leq z \leq 0.5\) selected from the WH15 cluster catalog. With a weight of cluster mass, a significant baryon acoustic oscillation (BAO) peak is detected on the correlation function with a significance of \(3.9 \sigma\). By fitting the correlation function with a \(\Lambda\)CDM model curve, we find \(D_v(z = 0.331) r_d^{fid}/r_d = 1269.4 \pm 58\) Mpc which is consistent with the Planck 2015 cosmology. We find that the correlation functions of the higher mass sub-samples show a higher amplitude at small scales of \(r < 80~h^{-1}{\rm Mpc}\), which is consistent with our precious result. We find a clear signal of the `Finger-of-God' effect on the 2D correlation function of the whole sample, which indicates the random peculiar motion of central bright galaxies in the gravitation potential well of clusters.