A new analytical solution based on the Ritz method is presented in this paper for analyzing the free vibration and buckling behavior of porous bi-directional functionally graded (2D-FG) beams under various boundary conditions. The solution is based on first-order shear deformation theory (FSDT). The selection of solution functions used in Ritz methods distinguishes the methods from each other and determines the accuracy of the analytical solution. To accurately capture the system's behavior and achieve the desired results, these functions have been carefully selected as a combination of polynomial and trigonometric expressions tailored as mixed series functions for each boundary condition. The study considers three types of porosity, namely PFG-1, PFG-2, and PFG-3. The equations of motion are derived using Lagrange's principle, taking into account the power-law variation of the beam material components throughout the volume. The non-dimensional fundamental frequencies and critical buckling loads are calculated for different boundary conditions, gradation exponents in the x and z directions (px, pz), slenderness (L/h), porosity coefficient (e), and porosity types. Initially, the accuracy of the mixed series functions is investigated for non-porous bi-directional functionally graded beams, and the numerical results are compared with existing literature to validate the proposed solution. Subsequently, the paper focuses on analyzing the influence of porosity on the free vibration and buckling behavior of bi-directional functionally graded beams using the developed solution method.