La aproximación automatizada de soluciones a ecuaciones diferenciales que involucran discontinuidades a través de superficies en evolución se aborda. La tecnología de elementos finitos se ha desarrollado hasta el punto en que ahora es posible modelar discontinuidades en evolución de forma independiente a la malla subyacente, lo cual es particularmente útil para simular fallas en sólidos. Sin embargo, el enfoque sigue siendo tedioso de programar, especialmente en el caso de problemas acoplados donde se emplean una variedad de bases de elementos finitos y donde se puede utilizar una mezcla de campos continuos y discontinuos. Abordamos este punto explorando la posibilidad de emplear técnicas de generación de código automatizado para modelar discontinuidades. Los espacios de funciones y formas variacionales se definen en un lenguaje que se asemeja a la notación matemática, y el código informático para modelar discontinuidades se genera automáticamente. Se esclarecen los principios subyacentes al enfoque y se presentan varios ejemplos bidimensionales y tridimensionales para diferentes ecuaciones.