Strain-driven (ɛD) and stress-driven (σD) two-phase local/nonlocal integral models (TPNIM) are applied to study the size-dependent nonlinear post-buckling behaviors of functionally graded porous Timoshenko microbeams. The differential governing equations and boundary conditions of motion are derived on the basis of the principle of minimum potential energy and expressed in nominal form. The integral relation between local and nonlocal stresses is transformed unitedly into equivalent differential form with constitutive constraints for ɛD- and σD-TPNIMs. Bending deflection and cross-sectional rotation are derived explicitly with Laplace transformation for linear buckling. Taking into account boundary conditions and constitutive constraints, a general eigenvalue problem is obtained to determine linear buckling mode shape (LBMS) and nominal buckling load (NBL). Local and nonlocal LBMS based Ritz–Galerkin methods and general differential quadrature method (GDQM) based Newton–Raphson’s method are utilized to obtain the numerical solutions for linear and nonlinear NBLs. Numerical results show that nonlocal LBMS based Ritz–Galerkin method and GDQM would lead to same prediction for linear NBLs as analytical method, and nonlocal LBMSs based Ritz–Galerkin and GDQM based Newton–Raphson’s method would lead to exact same prediction for nonlinear post-buckling loads. However, local LBMS based Ritz–Galerkin method is failed to provide accurate prediction for both linear NBLs and nonlinear post-buckling loads, through the difference between local and nonlocal LBMSs is not significant. The influence of nonlocal parameters, porous distribution patterns and boundary conditions on the linear buckling and nonlinear post-buckling behaviors are investigated numerically.