In this work the nonlinear behavior of layered SAW resonators is studied with the help of Finite Element (FE) computations. The full calculations depend strongly on the availability of accurate tensor data. While there are accurate material data for linear computations, the complete sets of higher-order material constants, needed for nonlinear simulations, are still not available for relevant materials. To overcome this problem, scaling factors were used for each available nonlinear tensor. The approach here considers piezoelectricity, dielectricity, eletrostriction and elasticity constants up to fourth order. These factors act as a phenomenological estimate for incomplete tensor data. Since no set of fourth order material constants for LiTaO3 is available, an isotropic approximation for the fourth order elastic constants was applied. As a result, it was found that the fourth order elastic tensor is dominated by one fourth order Lamé constant. With the help of the FE model, derived in two different, but equivalent ways, we investigate the nonlinear behavior of a SAW resonator with a layered material stack. The focus was set to third order nonlinearity. Accordingly, the modeling approach is validated using measurements of third order effects in test resonators. In addition, the acoustic field distribution is analyzed.