The mechanical properties of ultrahigh-strength steel (UHSS) produced through the gradient thermoforming process (GTP) are significantly influenced by residual stress and deformation. However, the precise distribution of these factors during GTP remains unclear. To address this issue, a detailed thermal model is developed, incorporating key heat transfer mechanisms, including coupled heat conduction and radiation, using the method of separation of variables. Furthermore, the displacement and stress fields of UHSS rectangular plates are derived using the differential quadrature method (DQM), accounting for phase transitions and external load. The combination of the separation of variables and DQM methods significantly improves computational efficiency. Furthermore, this study presents the first analysis of the combined effects of thermal load, phase transitions, and external forces on UHSS rectangular plates. Finally, through a combination of theoretical analysis, experimental validation, and case studies, the influence of material properties, geometric parameters, and external load on the temperature, stress, and displacement fields is examined. The findings reveal that heating rate, phase transitions, and external load play a critical role in the stress and displacement distribution in UHSS plates during GTP. This research provides a theoretical foundation for optimizing GTP process parameters and material design, ultimately enhancing the quality and performance of UHSS components.