We consider a constitutive model for the behavior of elastic flexoelectric materials including strain gradient fields and polarization gradient fields. This model is based on a stored elastic energy density function which depends on four independent variables: the polarization field and the polarization field gradient as well as the strain field and the strain field gradient. A generalized Toupin variational approach is utilized to find the governing equations (constitutive relations, equilibrium equations and boundary conditions) of the material. The present model is then applied to the problem of a thick walled cylindrical tube of elastic isotropic flexoelectric material, subjected to axisymmetric loading. The resulting radial displacement field noticeably differs from the elastic and strain gradient elastic cases.