We introduce a suitable notion of generalized Hessian and show that it can be used
to construct approximations by means of piecewise linear functions to the solutions of variational
problems of second order. An important guideline of our argument is taken from the theory of the
Gamma-convergence. The convergence of the method is proved for integral functionals whose integrand is
convex in the Hessian and satisfies standard growth conditions