Subsurface ocean tides act as a mechanism to dissipate tidal energy in icy satellite interiors. We numerically model the effect of an ice shell on ocean tides using non-linear bottom drag for the first time. We demonstrate that subsurface oceans experience tidal pressurization due to the confining nature of the ice shell, and find that Enceladus' eccentricity forcing can generate up to 2.2 kPa of pressure excess at the ocean surface. Existing free surface oceanic energy dissipation scaling laws are extended to subsurface oceans, and are benchmarked against our numerical results to within 10 %. We show that for the large bodies Ganymede, Europa and Titan, an ice shell increases eccentricity tidal heating due to self-gravity, whereas the shel...