I had a fairly useless comment posted earlier which I apologize for. It turns out much is known though. Below everything is restricted to zero applied field. For a finite square grid of size $m$ by $n$ with periodic boundary conditions, the expression for the partition function $Z$ (from which the entropy per site can be worked out) was given already by [Kaufmann](https://doi.org/10.1103/PhysRev.76.1232) (note that you can find freely accessible copies of this paper through google) very shortly after Onsager's solution: $Z=\frac{1}{2}(2\sinh(2H))^{mn/2}\left\{\prod_{r=1}^n\left(2\cosh(\frac{m}{2}\gamma_{2r-1})\right)+\prod_{r=1}^n\left(2\sinh(\frac{m}{2}\gamma_{2r-1})\right)+\prod_{r=1}^n\left(2\cosh(\frac{m}{2}\gamma_{2r})\right)+\prod_{r=1}^n\left(2\sinh(\frac{m}{2}\gamma_{2r})\right)\right\}$ where, $H=J/k_BT$, $H'=J'/k_BT$, $\cosh\gamma_j=\cosh(2H)\cosh(2H')-\sinh(2H)\sinh(2H')\cos(\pi j/n)$, with $J,J'$ being the coupling constants in the horizontal and vertical directions, respectively. A fuller explanation is in chapter IV of McCoy and Wu's book "The Two-dimensional Ising model". Another way to derive the partition function which has proved useful for [generalizations to antiperiodic boundary conditions](https://doi.org/10.1088/0305-4470/35/25/304) was worked out in [this paper by V.N. Plechko](https://doi.org/10.1007/BF01017042). On an infinitely long cylinder with circumference $L$ I was able to find in a [2004 paper by Huang et al](https://arxiv.org/abs/cond-mat/0407731v3) an expression for the free energy (with $J=J'=1$): $f=-\frac{1}{2}\ln(4z)-\frac{1}{2L}\sum_p\int_0^{2\pi}\frac{d\phi}{2\pi}\ln\left[z+z^{-1}-\Phi_p(\phi)\right]$ where now $z=\sinh(2\beta)$ and $\Phi_p(\phi)=\cos\phi+\cos\frac{2\pi p}{L}$.