The size-dependent nonlinear thermomechanical vibration analysis of pre- and post-buckled tapered two-directional functionally graded (2D-FG) microbeams is presented in this study. In the context of the modified couple stress theory, the formulations are derived based on the parabolic shear deformation beam theory and von Karman nonlinear strains. Different thermomechanical material properties are assumed to be temperature-dependent and smoothly vary in both length and thickness directions using the power law and the physical neutral axis concept is employed. The nonlinear governing equations are derived using the Hamilton principle and the resulting variable coefficient equations of motion are solved using the differential quadrature method (DQM) and iterative Newton’s method for clamped-clamped and simply supported boundary conditions. Comparison studies are presented to validate the derived model and solution procedure. The impacts of induced thermal moments, temperature power index, two gradient indices, nonuniform cross-section, and microstructure length scale parameter on the frequency-temperature configurations are explored for both clamped and simply supported microbeams