In this paper, numerical methods for wrinkle prediction of thin membrane were studied by finite element analysis. Techniques using membrane and shell elements were applied for triangular membrane. In case of membrane element method, the wrinkling was accounted for by the wrinkle algorithm of property modification, which was implemented to ABAQUS as a user subroutine. In case of shell method, geometrically nonlinear post-buckling analysis was performed to obtain the wrinkle deformation explicitly. The wrinkling deformation was induced by seeding the mesh with a random geometric imperfection. The results were investigated focusing on the mesh convergence and the solution accuracy.