Abstract: The existing methods are difficult to construct hexahedral meshes. To address this problem, an interactive construction method of complex hexahedral mesh model based on volumetric subdivision is presented. The user firstly needs to interactively construct the model skeleton, by placing cubes at the nodes of the skeleton structure, and performing interactive operations such as rotation, translation, and scaling of the node cubes. Then, through the connection between the nodes and the topological split operation, the initial control mesh can be constructed. Further, interpolatory Catmull-Clark volumetric subdivision method is used to generate hexahedral meshes with different resolutions. Finally, the padding operation is used to eliminate the degraded elements at the boundary and improve the quality of the hexahedral mesh elements to obtain the final hexahedral mesh. The results of numerical example show that the method can easily and efficiently generate hexahedral meshes interactively, and eliminates the intermediate steps of generating volume meshes from surface meshes. It has some practical applications in finite element analysis, isogeometric analysis and geometric modeling in animation.