This manuscript develops for the first time a mathematical formulation of the dynamical behavior of bi-directional functionally graded porous plates (BDFGPP) resting on a Winkler–Pasternak foundation using unified higher-order plate theories (UHOPT). The kinematic displacement fields are exploited to fulfill the null shear strain/stress at the bottom and top surfaces of the plate without needing a shear factor correction. The bi-directional gradation of materials is proposed in the axial (x-axis) and transverse (z-axis) directions according to the power-law distribution function. The cosine function is employed to define the distribution of porosity through the transverse zdirection. Equations of motion in terms of displacements and associated boundary conditions are derived in detail using Hamilton’s principle. The two-dimensional differential integral quadrature method (2D-DIQM) is employed to transform partial differential equations of motion into a system of algebraic equations. Parametric analysis is performed to illustrate the effect of kinematic shear relations, gradation indices, porosity type, elastic foundations, geometrical dimensions, and boundary conditions (BCs) on natural frequencies and mode shapes of BDFGPP. The effect of the porosity coefficient on the natural frequency is dependent on the porosity type. The natural frequency is dependent on the coupling of gradation indices, boundary conditions, and shear distribution functions. The proposed model can be used in designing BDFGPP used in nuclear, marine, aerospace, and civil structures based on their topology and natural frequency constraints.