A variational formulation and the associated finite elemet equations have been derived for general three-dimensional deformation of a planar anisotropic rigid-plastic sheet metal which obeys the strain-rate potential proposed by BARLAT et al . By using the natural convected coordinate system, the effect of geometric change and the rotation of planar anisotropic axes are considered efficiently. In order to check the validity of present formulation, a cylindrical cup and a square cup deep drawing test was modeled. good agreement was found between the FE simulation and the experiment. The results have shown that the present formulation for planar anisotropic deformation can be efficiently applied to the analysis of sheet metal working processes for planar anisotropic nonferrous metals.