SUMMARYA new, computationally efficient algorithm has been implemented to solve for thermal stresses, strains, and displacements in realistic solidification processes which involve highly nonlinear constitutive relations. A general form of the transient heat equation including latent-heat from phase transformations such as solidification and other temperature-dependent properties is solved numerically for the temperature field history. The resulting thermal stresses are solved by integrating the highly nonlinear thermo-elastic-viscoplastic constitutive equations using a two-level method. First, an estimate of the stress and inelastic strain is obtained at each local integration point by implicit integration followed by a bounded Newton-Raphson (NR) iteration of the constitutive law. Then, the global finite element equations describing the boundary value problem are solved using full NR iteration. The procedure has been implemented into the commercial package Abaqus (Abaqus Standard Users Manuals, v6.4, Abaqus Inc., 2004) using a user-defined subroutine (UMAT) to integrate the constitutive equations at the local level. Two special treatments for treating the liquid/mushy zone with a fixed grid approach are presented and compared. ISS, vol. 76, 1993) specialized in thermo-mechanical modelling of continuous casting. Both finite element codes are then applied to simulate temperature and stress development of a slice through the solidifying steel shell in a continuous casting mold under realistic operating conditions including a stress state of generalized plane strain and with actual temperaturedependent properties. Other local integration methods as well as the explicit initial strain method used in CON2D for solving this problem are also briefly reviewed and compared.