This paper proposes a method to jointly match multiple 3D meshes by maximizing pairwise feature affinities and cycle consistency across models. It formulates the matching problem as a low-rank matrix recovery problem and uses nuclear norm relaxation for rank minimization. An alternating minimization algorithm is used to efficiently solve the optimization problem. Experimental results show the method provides an order of magnitude speed-up compared to state-of-the-art algorithms based on semi-definite programming, while achieving competitive performance. It also introduces a distortion term to the pairwise matching to help match reflexive sub-parts of models distinctly.