Computing time-ordered products via Wick's theorem is fairly straightforward, schematically,
$$\mathcal{T} \left\{ \phi_1 \phi_2 \dots \phi_n\right\} = \; : \sum \mathrm{all \;possible \; contractions}:$$
where colons denote normal-ordering, and for simplicity we have chosen a real scalar field, and the notation $\phi_n$ denotes a field evaluated at the point $x_n$. A contraction of a pair of scalar fields at different points is the propagator evaluated at the difference of the respective space-time points, i.e. the Feynman propagator $\Delta_{ij} \equiv \Delta(x_i-x_j)$. Consider as a first example the product,
$$\mathcal{T} \left\{ \phi_1 \phi_2 \phi_3 \phi_4 \right\} = \; :\phi_1 \phi_2 \phi_3 \phi_4 : + \Delta_{12} :\phi_3\phi_4: + \Delta_{34}:\phi_1 \phi_2:$$
$$+ \Delta_{13}:\phi_2 \phi_4: + \Delta_{24}:\phi_1 \phi_3: + \Delta_{14}:\phi_2 \phi_3: + \Delta_{23}:\phi_1 \phi_4:$$
If we consider, for example, a Dirac field then we must take into account that they anti-commute when expanding the time-ordered product. For every swap, we introduce a minus sign. In addition, the propagator is only constructed from the contraction of a Dirac conjugate field with another field; any other type of contraction vanishes. In your example,
$$\mathcal{T} \left\{ \bar{\psi}_x \psi_x \bar{\psi}_y \psi_y \right\} = \; :\bar{\psi}_x \psi_x \bar{\psi}_y \psi_y: + \Delta_{xx}:\bar{\psi}_y \psi_y:$$
$$+ \Delta_{yy}:\bar{\psi}_x \psi_x: + \Delta_{xy}:\psi_x \bar{\psi}_y: + \Delta_{yx} :\bar{\psi}_x \psi_y:$$
Notice in this case there are no overall minus signs because there were an even number of required operator swaps. For additional resources on Wick's theorem, I recommend:
- Tong's lecture notes, they provide a proof of Wick's theorem, followed by a few examples, as well as applications of the theorem in computing scattering amplitudes. It is also shown which terms in the expansion correspond to certain processes, and how one selects which yield contributions to a particular Feynman diagram.
- Srednicki's Quantum Field Theory which includes a treatment of Wick's theorem, with examples, as well as problems to work through.
- Tutorial on Wick's Theorem which features a long computation using Wick's theorem, in complete explicit detail, and shows how to handle certain other insertions in time-ordered products, such as gamma matrices.
Unfortunately, you will not be able to find dozens of examples in a single resource because the computations are quite mechanical and tedious. In addition, the approach to scattering amplitudes involving Wick's theorem is mostly obsolete, and the preferred approach uses Feynman rules, and Feynman diagrams, as well as methods built upon them.