A combined surface and bulk TCAD damage model for the analysis of radiation detectors operating at HL-LHC fluences