A direct elimination scheme is presented for the solution of large scale voxel vibration analysis. The regularity of the stiffness matrix of the voxel models is used for decision of the order of elimination. By this process, elimination of degrees of freedoms can be performed without storing all non-zero profile of the global coefficient matrix. Sample solutions verify that large-scale voxel models can be analyzed with the smaller amount of in-core memory and out-core device storage.