The low swelling property of magnesium oxide concrete is a significant feature that can be used to control the cracking of mass concrete. Based on the characteristics of the chemical reaction, this work proposes a coupled hydro-thermo-mechanical model that can be implemented with the finite element method for predicting the autogenous volumetric deformation of magnesium concrete. By introducing the degree of the hydration reaction of magnesia and the degree of the hydration reaction of cementitious materials as intermediate variables of the chemical reaction system, a prediction model of the concrete temperature and chemical fields is established, and using this model, the effect of the temperature on the reaction rate can be considered in real time. In addition, by combining the relationship between the degree of the hydration reaction of magnesium oxide and the comprehensive expansion of concrete, a mathematical model for calculating the expansion stress of magnesia concrete was established. The algorithms were derived by mathematical equations, and the simulation results were compared to the experimental temperature and autogenous volumetric strain curves, which showed that the hydration model provides a relatively high accuracy. The model was also applied to an arch dam, and the coupled thermo-chemical-mechanical responses of mass concrete during construction were investigated. Simulation results show that the increase in temperature (hydration of cementitious material) and expansion volumetric deformation (hydration of MgO) of the concrete on the upstream and downstream surfaces lags obviously behind that of the inner regions. Quantitative analysis for differences of internal and external expansion is worthy of further attention and study on a basis of further experimental data as well as monitored data.