The batch comminution equation for multi-component mineral systems is a multidimensional integrodifferential equation [GRAPHICS] that cannot be solved analytically, The phase space (g, D) was discretized and the equation formulated as [GRAPHICS] Solutions to this discretized form of the equation are easy to generate provided that models for the selection function S-ij and the multi-component breakage function b(ijkl) are known, A general solution to the batch comminution equation that does not rely on the random fracture assumption is derived in this paper and is given by [GRAPHICS] with the coefficients given by a(ijmj) = 0 if i not equal m [GRAPHICS] This solution includes the following non-random fracture processes: selective phase breakage, differential breakage, preferential breakage, liberation by detachment and boundary-region fracture. This solution was compared with experimental data obtained from batch tests on two ores. The model was found to be reliable and the derived Andrews-Mika diagram is shown. The derived models for the selection function and the Andrews-Mika diagram can be used to simulate the effects of the processes that produce non random fracture, and thus to assess their relative importance, and also to simulate the behavior of the ore in continuous milling and concentrating circuits. An assumption of random fracture provides an important simplification because, under this assumption, the selection function is independent of the particle composition. The random fracture assumption has been used by researchers who developed models for the prediction of mineral liberation by comminution. The solution of the batch comminution equation under this assumption was found to be [GRAPHICS] with the coefficients given by [GRAPHICS] and [GRAPHICS]