A hybrid model for an anaerobic digestion process is proposed. The fermentation is assumed to be performed in two steps, acidogenesis and methanogenesis, by two bacterial populations. The model is based on mass balance equations, and the bacterial growth rates are represented by neural networks. In order to guarantee the biological meaning of the hybrid model (positivity of the concentrations, boundedness, saturation or inhibition of the growth rates) outside the training data set, a method that imposes constraints in the neural network is proposed. The method is applied to experimental data from a fixed bed reactor.