The present article describes an analytical approach to investigate the nonlinear free vibrational behavior of a two-directional functionally graded porous (TDFGP) cone-shaped shell resting on elastic substrates. A mixture model with the power-law distribution functions is utilized to obtain the TDFGP shell’s mechanical properties. In addition, it is assumed that the internal porosities in the matrix materials can be dispersed into two independent patterns, either even or uneven porosity distribution. The field displacement in the TDFGP shell is approximated using first-order shear deformation theory. Hamilton’s principle and von Karman’s large deformation assumptions give a dynamic model for the TDFGP shell. The Galerkin decomposition method is employed to discretize the governing partial differential equations into time-dependent ordinary differential equations. The harmonic balance method is used to analytically solve the final time-dependent equations and acquire the nonlinear frequency response of the TDFGP shell. Research has been conducted on porosity, gradient indexes, elastic substrates, boundary conditions, and geometric ratios through parametric studies to investigate the outcomes and their applicability to real-world issues.