We present large-eddy simulations (LES) of separation and reattachment of a flat-plate turbulent boundary-layer flow. Instead of resolving the near wall region, we develop a two-dimensional virtual wall model which can calculate the timeand space-dependent skin-friction vector field at the wall, at the resolved scale. By combining the virtual-wall model with the stretched-vortex subgrid-scale (SGS) model, we construct a self-consistent framework for the LES of separating and reattaching turbulent wall-bounded flows at large Reynolds numbers. The present LES methodology is applied to two different experimental flows designed to produce separation/reattachment of a flat-plate turbulent boundary layer at medium Reynolds number Re θ based on the momentum boundary-layer thickness θ. Comparison with data from the first case at Re θ = 2000 demonstrates the present capability for accurate calculation of the variation, with the streamwise co-ordinate up to separation, of the skin friction coefficient, Re θ , the boundary-layer shape factor and a non-dimensional pressure-gradient parameter. Additionally the main large-scale features of the separation bubble, including the mean streamwise velocity profiles, show good agreement with experiment. At the larger Re θ = 11 000 of the second case, the LES provides good postdiction of the measured skin-friction variation along the whole streamwise extent of the experiment, consisting of a very strong adverse pressure gradient leading to separation within the separation bubble itself, and in the recovering or reattachment region of strongly-favourable pressure gradient. Overall, the present two-dimensional wall model used in LES appears to be capable of capturing the quantitative features of a separation-reattachment turbulent boundary-layer flow at low to moderately large Reynolds numbers.