Curing of epoxy resins poses a particular challenge in terms of modeling, experimental investigation, and numerical implementation, as it is a thermo-chemo-mechanical process. Several constitutive relations are required to model these processes, yielding numerous material parameters. The calibration of the constitutive relations must be performed using multiple steps, wherein uncertainties unavoidably propagate. In this study, we investigate the propagation of uncertainties during both the multi-step calibration procedure and the numerical simulation of curing processes with the identified parameters. For both, we employ the first-order second-moment method, which is carefully evaluated through coverage tests and by comparing it to the Monte Carlo method as a reference. It is demonstrated that the first-order second-moment method efficiently yields reasonable results, although providing only a first-order approximation of the highly nonlinear stochastic model response.