Efficient thermo-mechanical model for solidification processes and its applications in steel continuous casting

A 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 thermoelastic-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 iteration of the constitutive law. Then, the global finite element equations describing the boundary value problem are solved using full Newton-Raphson iteration. The procedure has been implemented into the commercial package Abaqus [1] 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. 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. The model is validated both with a semi-analytical solution from Weiner and Boley [2] as well as with an in-house finite element code CON2D [3,4,7,8] specialized in thermomechanical modeling of continuous casting. Both finite element codes are then applied to simulate temperature and stress development of a slice through the solidifying steel shell