This work presents a comparative analysis of thermal models for photovoltaic modules using Grammatical Evolution (GE) and Differential Evolution (DE) across four photovoltaic (PV) technologies: crystalline silicon (c-Si), amorphous silicon (a-Si), cadmium telluride (CdTe), and organic (OPV), under three sky conditions: sunny, cloudy, and diffuse. Temperature data were collected through a monitoring system in a photovoltaic cube, measuring temperatures on the horizontal face (top PV module) as well as environmental parameters (irradiance, ambient temperature, wind speed and direction, relative humidity). Three empirical models from the literature (Sandia, Faiman, and Obiwulu) were compared with 10 models generated by GE+DE using a Global Performance Index (GPI) to evaluate the accuracy of the models, considering five statistical metrics. The results show that in 11 out of 12 scenarios, the generated models outperform the empirical models, highlighting the importance of relative humidity in model accuracy. This work extends previous research, providing more accurate predictive models for the operating temperature of photovoltaic modules.