We develop an enrichment wall-model within the spectral element method (SEM) framework for large-eddy simulations of wall-bounded flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions. The performance of the SEM enrichment wall-modeling approach is assessed in turbulent channel flow LES for a range of Reynolds numbers. The results show that the proposed method improves solution accuracy on under-resolved near-wall grids as compared to the traditional shear stress wall-models.