This paper presents an isogeometric analysis (IGA) based on Bézier extraction to investigate the stability of a circular tunnel in cohesive-frictional soils subjected to uniform surcharge loading. The basic concept of IGA is to use non-uniform rational B-splines (NURBS) basis functions for both the parametrisation of the geometry in computer-aided design (CAD) and the finite approximation for analysis. In this study, IGA based on Bézier extraction is applied to model a circular tunnel's geometry discretisation. The soil is modelled to be a homogeneous material obeying the Mohr-Coulomb failure criterion. The ultimate surcharge loading and failure mechanisms of the circular tunnel are obtained via upper bound limit analysis. Based on the results of isogeometric analysis, design tables and charts are generated for various soil properties, internal friction angles and ratios of tunnel cover to tunnel diameter. The numerical results for plane strain circular tunnels are compared with those of other available solutions to prove the reliability and efficiency of the present method.