The development of high-brightness X-ray free electron lasers (XFELs), such as hard X-ray self-seeding free electron lasers and XFEL oscillators (XFELOs), brings a severe challenge to the crystal monochromator due to a strong non-uniform thermal load. The distortion caused by spatial temperature gradients can severely affect the optical performance of crystals. Therefore, this paper presents a model to estimate the performance of non-uniform thermally distorted crystals. The model not only takes into account thermal strain, slope error and incident angle deviation, but also considers temperature-dependent factors such as the Debye–Waller factor and electric susceptibility. Our investigation indicates that the Debye–Waller factor reduces the height and bandwidth of rocking curves, and the impact of the electric susceptibility is tiny. The proposed model can describe the distortion of the reflectivity and transmissivity curves of non-uniform thermally loaded crystals and can be applied in the design of crystal monochromators, crystal splitters, crystal compressors and XFELOs.