Abstract A study of existence and uniqueness of weak solutions to boundary value problems describing an elastic body with weakly nonlocal surface elasticity is presented. The chosen model incorporates the surface strain energy as a quadratic function of the surface strain tensor and the surface deformation gradients up to N th order. The virtual work principle, extended for higher‐order strain gradient media, serves as a basis for defining the weak solution. In order to characterize the smoothness of such solutions, certain energy functional spaces of Sobolev type are introduced. Compared with the solutions obtained in classical linear elasticity, weak solutions for solids with surface stresses are smoother on the boundary; more precisely, a weak solution belongs to where and .