Suppose your 3-form is tau; to extract the integrand, simply use tau[1,2,3].expr(), or tau[frame,1,2,3,chart].expr() if the vector frame and the chart are not the default ones on the manifold. Note that expr() returns a symbolic expression, which can be passed directly to integrate. For instance, for the example of question 41228, this would be
sage: tau.display(spherical.frame(), spherical)