The aim of this work is to set up a numerical model in order to simulate skull fracture, and to analyse the effects of the material discontinuity induced by this non-linearity on intracranial stress field. Mechanical parameters of cranial bone were found in the bibliography for the fmite element modelling of skull fracture. A failure criterion was validated with cranial bone samples and is applied to a new 3D anatomic skull model. The results comply well with the experimental bone response found in the bibliography. Finally, a simple spherical human head model was used to analyse the effects of skull fracture on head response in case of shock.