The aim of this work is to present a computationally efficient algorithm to simulate thedeformations suffered by a viscoplastic body in a solidification process. This type ofproblems involves a nonlinearity due to the considered thermo-elastic-viscoplastic law. Inour previous papers, this difficulty has been solved by means of a duality method, knownas Bermúdez–Moreno algorithm, involving a multiplier which was computed with a fixed pointalgorithm or a Newton method. In this paper, we will improve the former algorithms bymeans of a generalized duality method with variable parameters and we will presentnumerical results showing the applicability of the resultant algorithm to solidificationprocesses. Furthermore, we will describe a numerical procedure to choose a constantparameter for the Bermúdez–Moreno algorithm which gives good results when it is applied tosolidification processes.
