Thermal aging significantly deteriorates the mechanical properties and service performance of epoxy resins used for the high-voltage bushing. Current studies on the thermal aging behavior of epoxy resins mainly focus on experimental observations. However, an in-depth understanding of the mechanism of thermal aging of epoxy resins requires the monitoring of structural evolution of epoxy resins during thermal aging at the molecular level. To thoroughly analyze the intrinsic factors affecting structural evolution and the effect of thermal aging on the mechanical properties of epoxy resin for high-voltage bushing, epoxy resin models with different crosslinking degrees were established and thermal aging treatments at various temperatures and time periods were carried out by molecular dynamics simulation. It was found that the tensile strength of the epoxy resin was enhanced with the increase of the crosslinking degree, which was related to the elevation of the proportion of C-N and O-H bonds in its structure. With the increase of thermal aging temperature, the tensile strength of the epoxy resin decreased, which was related to the formation of weak bonds. At the early stage of thermal aging and after a period of high-temperature thermal aging, the strength of epoxy resin significantly decreases. The thermal aging of the epoxy resin is accelerated under external loading. In addition, the crosslinking degree and curing agent also affect the thermal aging resistance of epoxy resins. The results of this study can provide guidance for predicting and improving the thermal aging resistance of epoxy resins. Materials Studio was used to construct molecular models and complete crosslinking reactions. DGEBA and 44DDS (or 33DDS) were mixed at a ratio of 2:1, followed by crosslinking reaction. During this process, the Nose method was used to control temperature, the Berendsen method was used to control pressure, and the polymer consistent force field (PCFF) was used to control the motion and force of atoms. Isobaric-isothermal ensemble (NPT ensemble) was used to heat up epoxy resin models to various thermal aging temperatures of 400K, 500K, 600K and 700K. The models were maintained at these temperatures for different thermal aging times of 100ps, 200ps, 300ps, 400ps, 500ps, 600ps, 700ps and 800ps. Afterwards, the models were cooled down to 300K and subjected to uniaxial tensile testing at this temperature with a strain rate of 1 × 109s-1. The structural configurations and stress-strain data during the tensile process were recorded. The flow stress of the material was derived by counting the average stress in the 20-50% strain interval.