In this paper, analysis of linear and nonlinear buckling of relatively thick orthotropic graphene sheets is carried out under mechanical load based on elasticity theories. With the help of nonlocal elasticity theory, the principle of virtual work, first order shear theory and Von-Karman nonlinear strain, the dominant relationship in terms of obtained displacements has been obtained, and the method of differential quadrature (DQ) with non-uniform distribution points (Chebyshev -Gauss-Lobato) has been used. To check the validity, the obtained results have been compared with other resources, and the effects of nonlocal coefficient, thickness, radius and elastic base on the dimensionless buckling loads were calculated and investigated. Moreover, the results of analysis using nonlocal and local theory were compared together. It can be noticed that the dimensionless buckling loads on graphene sheets increased more with a decrease in flexibility as far boundary condition is concerned. Additionally, with an increase in the sheet radius, the variation between nonlocal and local analysis results will be more.