Results are reported from wall-modelled large eddy simulations (WMLES) of a zero pressure gradient flat-plate turbulent boundary layer (TBL) flow performed using unstructured computational meshes. In particular, two meshes are considered: a hex-dominant and a polyhedral. The resolution of the meshes is kept constant with respect to the local thickness of the TBL. The WMLES predictions are evaluated by comparison with reference data from direct numerical simulation (DNS) and semi-empirical expressions for the development of integral quantities along the TBL. Good agreement is observed for the skin friction coefficient, mean streamwise velocity and the Reynolds stresses. Also, the influence of the location of the sampling (matching) point of the employed algebraic wall-stress model is investigated. It is found that moving the sampling point to the third consecutive off-the-wall cell centre leads to a significant improvement in the prediction of the mean wall shear stress, as opposed to sampling for the wall-adjacent cell.