We present a computational model of amygdala neural networks. It is used to simulate neuronal activation in amygdala nuclei at different stages of aversive conditioning experiments with rats. Our model is based on neurobiological data. Simple formal neurons and an adaptive Hebbian rule are the key elements of the model. The results are compatible with neuronal activation maps obtained with C-Fos markers. The model also enables interesting predictions.