A sub-domain method for free vibration analysis of arbitrarily shaped, concave membranes with a deep groove is proposed in the paper. The proposed method divides the concave membrane of interest into two convex regions. The vibration displacement(approximate solution) of each convex region is assumed by linearly superposing plane waves generated at edges of the region. A sub-system matrix for each convex region is extracted by applying a provisional boundary condition to the approximate solution. Finally, a system matrix, which of the determinant gives eigenvalues of the concave membrane, is made by considering the fixed boundary condition(displacement zero condition) at edges and the compatibility condition(the condition of continuity in displacement and slope) at the interface between the two regions. Case studies show that the proposed method is valid and accurate when the eigenvalues by the proposed are compared to those by NDIF method, FEM, or the exact method.