▎ 摘　　要

Due to the inevitable presence of random defects, unpredictable grain boundaries in macroscopic samples, stress concentration at clamping points, and unknown load distribution in the investigation of graphene sheets, uncertainties are crucial and challenging issues that require more exploration. The application of the Kriging surrogate model in vibration analysis of graphene sheets is proposed in this study. The Latin hypercube sampling method effectively propagates the uncertainties in geometrical and material properties of the finite element model. The accuracy and convergence of the Kriging surrogate model are confirmed by a comparison with the reported references. The uncertainty analysis for both Zigzag and Armchair graphene sheets are compared and discussed.