  We present a minimalistic approach to simulations of force transmission
through granular systems. We start from a configuration containing cohesive
(tensile) contact forces and use an adaptive procedure to find the stable
configuration with no tensile contact forces. The procedure works by
sequentially removing and adding individual contacts between adjacent beads,
while the bead positions are not modified. In a series of two-dimensional
realizations, the resulting force networks are shown to satisfy a linear
  constraint among the three components of average stress, as anticipated by
recent theories. The coefficients in the linear constraint remain nearly
  constant for a range of shear loadings up to about .6 of the normal loading.
The spatial distribution of contact forces shows strong concentration along
``force chains". The probability of contact forces of magnitude f shows an
exponential falloff with f. The response to a local perturbing force is
concentrated along two characteristic rays directed downward and laterally.
