El presente documento propone un desarrollo matemático de los enfoques de plasticidad y daño para simular procesos de conformado de chapa metálica. Se enfoca en la predicción numérica de la deformación de la chapa metálica durante el proceso de embutición profunda cuando aparece una grieta. Se proponen ecuaciones constitutivas de plasticidad anisotrópica. Se utiliza una integración totalmente implícita de las ecuaciones constitutivas de acoplamiento y conduce a dos ecuaciones escalares locales no lineales que se resuelven mediante el método de Newton. El modelo desarrollado permite predecir el inicio de grietas en chapas metálicas durante operaciones de conformado en frío. El modelo numérico se implementa en el software ABAQUS utilizando subrutinas definidas por el usuario, que son VUMAT y UMAT. La precisión del modelo elastoplástico anisotrópico completamente acoplado con daño dúctil se evalúa mediante ejemplos numéricos.