  We compute the relaxation of the total energy related to a variational model
for nematic elastomers, involving a nonlinear elastic mechanical energy
depending on the orientation of the molecules of the nematic elastomer, and a
nematic Oseen--Frank energy in the deformed configuration. The main assumptions
are that the quasiconvexification of the mechanical term is polyconvex and that
the deformation belongs to an Orlicz-Sobolev space with an integrability just
above the space dimension minus one, and does not present cavitation. We
benefit from the fine properties of orientation-preserving maps satisfying that
regularity requirement proven in \cite{HS} and extend the result of \cite{MCOl}
to Orlicz spaces with a suitable growth condition at infinity.
