A dimensionless discrete state-space mathematical model is proposed for stress analysis of an axially-symmetric saturated poroelastic circular plate being loaded by temperature rise of one of the flat faces. Computation of stress by means of this model is free from a need for differentiation. The curves of thermal postbuckling, thermal nonlinear bending, radial stress, circumferential stress and resultant radial force have been achieved by means of six state-space variables and differential quadrature procedure. The effects of loading temperature, boundary conditions, porosity, pore distribution and fluid pore pressure on the thermal post buckling and nonlinear bending phenomena have been investigated. The critical postbuckling thermal load has been determined by means of a numerical algorithm that does not converge to trivial solution because it starts from a relatively great numerical thermal load. Contrary to the governing equation in usual form, the governing equations in state space together with the state equations equalize the number of the governing assembled difference equations and that of the unknowns, and eventually do not require to apply the so-called auxiliary extra boundary conditions at the nodes adjacent to the boundary nodes. As another advantage, in state-spatial viewpoint, there is not any necessity for computing the coefficients of derivatives with order higher than unity.