1. Isometric muscle force and the surface electromyogram (EMG) were simulated from a model that predicted recruitment and firing times in a pool of 120 motor units under different levels of excitatory drive. The EMG-force relationships that emerged from simulations using various schedules of recruitment and rate coding were compared with those observed experimentally to determine which of the modeled schemes were plausible representations of the actual organization in motor-unit pools. 2. The model was comprised of three elements: a motoneuron model, a motor-unit force model, and a model of the surface EMG. Input to the neuron model was an excitatory drive function representing the net synaptic input to motoneurons during voluntary muscle contractions. Recruitment thresholds were assigned such that many motoneurons had low thresholds and relatively few neurons had high thresholds. Motoneuron firing rate increased as a linear function of excitatory drive between recruitment threshold and peak firing rate levels. The sequence of discharge times for each motoneuron was simulated as a random renewal process. 3. Motor-unit twitch force was estimated as an impulse response of a critically damped, second-order system. Twitch amplitudes were assigned according to rank in the recruitment order, and twitch contraction times were inversely related to twitch amplitude. Nonlinear force-firing rate behavior was simulated by varying motor-unit force gain as a function of the instantaneous firing rate and the contraction time of the unit. The total force exerted by the muscle was computed as the sum of the motor-unit forces. 4. Motor-unit action potentials were simulated on the basis of estimates of the number and location of motor-unit muscle fibers and the propagation velocity of the fiber action potentials. The number of fibers innervated by each unit was assumed to be directly proportional to the twitch force. The area of muscle encompassing unit fibers was proportional to the number of fibers innervated, and the location of motor-unit territories were randomly assigned within the muscle cross section. Action-potential propagation velocities were estimated from an inverse function of contraction time. The train of discharge times predicted from the motoneuron model determined the occurrence of each motor-unit action potential. The surface EMG was synthesized as the sum of all motor-unit action-potential trains. 5. Two recruitment conditions were tested: narrow (limit of recruitment < 50% maximum excitation) and broad recruitment range conditions (limit of recruitment > 70% maximum excitation).(ABSTRACT TRUNCATED AT 400 WORDS)