My suggestion: implement a cutting algorithm for polyhedra. I mean split a polyhedron in two by means of a plane.
For every vertex you can compute the algebraic distance to the plane.
Consider every face of the polyhedron in turn and cut it with the plane. If all vertices are on the same side of the plane, there is no intersection. You will keep the face as such or completely discard it. If there are vertices on either sides, then you will keep some of the edges intact, split those that cross and discard others. You will reconstruct the cut face by joining the vertices in the proper order (after sorting the piercing points along the intersection line).
After doing that, you will have a new set of faces forming a polyhedron, with a cover face missing. You will reconstruct it by joining the piercing points, using the same edges that were used to close the faces. (Actually, you may end up forming several faces because the cross section can be made of several pieces.)
When you are able to cut a polyhedron with a plane, then you are able to find its intersection with an arbitrary convex polyhedron, such as a cube.
What I just described is a generalization to 3D of the Sutherland–Hodgman clipping algorithm.
There is a happy case when the polyhedron vertices are all on the same size of some of the planes. You may start the work by doing this test. But for the other cases, there is no real shortcut.
You can implement bounding box tests on the volumes, faces and edges, hoping to get speedup, but the more you put such tests, the less they are efficient.
The case of intersection with a convex polyhedron is easier to implement as all faces remain convex and the topological changes are simpler.