Development and numerical implementation for an elastoplastic constitutive model for anisotropic and asymmetric materials are presented in this paper. The Coulomb-Mohr yield criterion was modified to consider both the anisotropic and asymmetric properties. The modified yield criterion is an isotropic function of the principal values of a symmetric matrix which is linearly transformed from the Cauchy stress space. In addition to the constitutive equation, the numerical treatment for the singularity in the vertex region of yield surface and stress integration algorithm based on elastoplasticity were presented. In order to assess the accuracy of numerical algorithm, isoerror maps were considered. Also, extension of a strip with a circular hole was simulated and results compared with those obtained using the (smooth) Mises yield criterion to validate stress output for a complex stress state.