How to access normal vectors of each discretised tetrahedron element

Dear all,

I am trying to apply a variable compressive body force to a 3D cantilever model. The compressive force has no certain direction in a global scale.So I am trying to apply this compressive force on faces of each discretised tetrahedron element .Is there any way to access the normal vectors of each tetrahedron element?