Many yield functions have been proposed in academia to describe the complicated shapes of yield surfaces of metals. However, many of the commercial FEM codes provide only classical yield functions. Moreover, it takes a long time for engineers to implement newly proposed yield functions to commercial FEM codes. The Japan Association for Nonlinear CAE (JANCAE), a non-profit organization, developed the Unified Material Model Driver for plasticity (UMMDp) subroutine suite with the cooperation of industry users of CAE and the engineers of software venders. This subroutine provides several anisotropic yield functions, and is applicable to most of the commercial FEM codes. The users can implement their own anisotropic yield functions easily using the UMMDp. This paper presents the basic framework of UMMDp and the development activities performed with many volunteers.