  Two-component (normal and degenerate stars) models are the simplest
realization of clusters with a mass spectrum because high mass stars evolve
quickly into degenerates, while low mass stars remain on the main-sequence for
the age of the universe. Here we examine the evolution of isolated globular
clusters using two-component Fokker-Planck (FP) models that include heating by
binaries formed in tidal capture and in three-body encounters. Three-body
binary heating dominates and the postcollapse expansion is self-similar, at
least in models with total mass M <= 3 x 10^5 M_\odot, initial half-mass radius
r_{h,i} >= 5 pc, component mass ratio m_2/m_1 <= 2, and number ratio N_1/N_2 <=
300 when m_2=1.4 M_\odot. We derive scaling laws for \rho_c, v_c, r_c, and r_h
as functions of m_1/m_2, N, M, and time t from simple energy-balance arguments,
and these agree well with the FP simulations. We have studied the conditions
under which gravothermal oscillations (GTOs) occur. If E_{tot} and E_c are the
energies of the cluster and of the core, respectively, and t_{rh} and t_c are
their relaxation times, then \epsilon \equiv (E_{tot}/t_{rh})/(E_c/t_{rc}) is a
good predictor of GTOs: all models with \epsilon>0.01 are stable, and all but
one with \epsilon < 0.01 oscillate. We derive a scaling law for \epsilon
against N and m_1/m_2 and compared with our numerical results. Clusters with
larger m_2/m_1 or smaller N are stabler.
