This article examines the impact of inclined transverse cracking on the natural frequency of Euler-Bernoulli-model imperfect FG beams on an elastic foundation, taking into account the pinned-pinned and clamped-clamped boundary conditions. The equations are diffused using the classic finite element method (h-version). Material properties are considered to vary in both directions of the beam; width and thickness, using the power-law formula. The approximate porosity model is adopted with a uniform distribution. Cracked element stiffness is determined by reducing the cross-section of a bi-directional FG beam. The numerical results are compared with those of previous convergence studies for dimensionless fundamental frequencies. Case studies were carried out to evaluate the effect of the porosity values, the crack depth, the crack location, slope angle, and the Winkler-elastic foundation parameters on the first three natural frequencies of a beam with different support conditions and material gradient. In addition, the results demonstrated that the two-way distribution function plays an important role in the convergence of frequencies.