Computational systems, methods, and articles of manufacture to predict at least one of residual stresses and distortion in quenched aluminum castings. Residual stresses and distortion may be predicted through incorporating thermal strains induced during quenching with the nonlinear constitutive behavior of quenched microstructures of a quenched aluminum casting, wherein thermal strains arise generally from non-uniform transient temperature distribution of the casting during quenching. The transient temperature distribution of the aluminum casting during quenching may be calculated based on heat transfer coefficients specific to one or more nodes, elements and/or zones on the surfaces of the aluminum casting. The nonlinear constitutive behavior of the quenched aluminum casting may be modeled as functions of temperatures, strain rates, and microstructure variations. A material constitutive model accounts for not only strain hardening and creep, but also precipitate hardening. Thereby, residual stresses and distortion in the quenched aluminum castings can be predicted using the embodiments.