The occlusion effect (OE) refers to the phenomenon that more bone-conducted physiological sounds are transmitted into the earcanal when it is blocked and may cause discomfort on users of hearing protection devices. Models have been proposed to study the OE as they can help understand the physical mechanisms and can be used to evaluate the individual contribution on the OE of the factors that may affect it (i.e., occlusion device, ear anatomy, and stimulation). The existing finite element models developed to study the OE are limited by their truncated ear geometries. In order to progress in the understanding of the OE, the goal of this paper is to develop a finite element model of an entire head to predict the sound pressure field in its earcanals, open or occluded by earplugs. The model is evaluated by comparing the computed input mechanical impedances and OEs in various configurations with literature data. It is able to reproduce common behavior of the OE reported in the literature. In addition, the model is used to assess the effects on the simulated OEs of several parameters, including the modeling of the external air, the boundary condition at the head base and the material properties.