  We use numerical linked cluster expansions to study thermodynamic properties
of the two-dimensional spin-1/2 Ising, XY, and Heisenberg models with bimodal
random-bond disorder on the square and honeycomb lattices. In all cases, the
nearest-neighbor coupling between the spins takes values $\pm J$ with equal
probability. We obtain the disorder averaged (over all disorder configurations)
energy, entropy, specific heat, and uniform magnetic susceptibility in each
case. These results are compared with the corresponding ones in the clean
models. Analytic expressions are obtained for low orders in the expansion of
these thermodynamic quantities in inverse temperature.
