We establish the Euler-Lagrange inclusion of a nonsmooth integral functional defined on Orlicz-Sobolev spaces. This result is achieved through variational techniques in nonsmooth analysis and an integral representation formula for the Clarke generalized gradient of locally Lipschitz integral functionals defined on Orlicz spaces.