Concrete is a widespread material all over the world. Due to this material’s heterogeneity and structural complexity, predicting the behavior of concrete structures under extreme environmental conditions is a very challenging task. High temperatures lead to microstructural changes which affect the macrostructural performance. In this context, computational tools that allow the simulation of structures may assist the analysis, by reproducing varied situations of thermal and mechanical loading and boundary conditions. In order to contribute to this scenario, this study proposes a numerical methodology to simulate the thermomechanical behavior of concrete under temperature gradients, through inverse analyses and a user subroutine implemented in Abaqus software. Thermal loading effects were considered as loading data for a damage model. Experimental data available in the literature was adopted for adjustment and validation purposes. The preliminary results presented herein encourage further improvements so as to allow realistic simulations of such an important aspect of concrete’s behavior.