Accurate magnetostatic field computation for complex 3D ferromagnetic objects can be widely used to solve various engineering problems, while mesh quality is a key factor that can greatly affect the calculation precision. For complex shaped 3D objects, their mesh quality can be improved greatly by using hybrid mesh strategy with the advantages of hexahedral elements and tetrahedral elements. In this paper, a hybrid mesh integral equation method (IEM) for 3D magnetostatic field computation is realized based on mixed mesh with hexahedrons and tetrahedrons. Furthermore, computational accuracy of the proposed method is demonstrated by simulation and experimental examples with 2.10% and 3.20%, respectively. Thus, foundations for the mixed mesh strategy for the complex shaped 3D objects and fast computational methods with hybrid mesh IEM are established.