En el trabajo de investigación desarrollado en esta Tesis se propone una modelización completa y a la vez sencilla del transformador y trifásico de tres columnas para el funcionamiento en condiciones de saturación del núcleo. La modelización propuesta en este trabajo fue inicialmente desarrollada en el programa de simulación Pspice. Las simulaciones se han validado con un amplio experimental, en el cual se describe el proceso de obtención de los parámetros mediante ensayos de laboratorio fáciles de realizar. Uno de los principales objetivos de este trabajo es, que el modelo del transformador sea adecuado para integrarse a un programa de análisis de flujos de cargas con armónicos. Para esto, se ha desarrollado un algoritmo de resolución directa del régimen permanente del sistema de ecuaciones no lineal del transformador con saturación del núcleo. Las simulaciones del modelo desarrollado tanto en el programa de simulación Pspice como en el algoritmo de obtención del régimen permanente directo, se han validado contra las medidas de los ensayos de laboratorio. En ambos casos, las simulaciones han mostrado gran concordia con las medidas. Finalmente se hace una comparación entre los resultados de las simulaciones obtenidas con el programa de simulación y las obtenidas con el algoritmo. Los resultados son muy similares, lo cual indica que el algoritmo es adecuado.