You are on the right track. Remember for the 2D HO, you will need to add the extra dimension to the kinetic energy part of the hamiltonian. For the energy expression, since you have two different force constants (an anisotropic oscillator), your energy will be a sum of two terms for each degree of freedom, namely: Enxny=hvx(nx+1/2)+hvy(ny+1/2). Then, you should plug in nx=ny=0 to get the energy of the lowest state.
So, you don't really need to solve the TISE for the 2D HO, you just need to be able to extend the 1D case by an extra dimension.
Hope this helps.