In　this　paper,　the　nodal-based　continuous-discontinuous　deformation　analysis　method　(NCDDAM)　is　enriched　to　simulate　the　thermo-mechanical　coupling　effects　in　rock　masses.　A　distance-based　contact　potential　algorithm　is　first　incorporated　into　NCDDAM　to　avoid　the　dependency　of　element　shape　and　size　on　the　calculation　of　contact　force　between　different　blocks.　Then,　three　types　of　heat　conduction　models,　which　can　deal　with　heat　conduction　in　continuum　and　discontinuum,　are　incorporated　into　NCDDAM　to　simulate　the　heat　conduction　effects　of　rock　masses.　Finally,　a　two-way　staggered　algorithm　is　adopted　in　the　context　of　NCDDAM　to　simulate　the　thermo-mechanical　coupling　effects　in　rock　masses.　Several　benchmark　examples　are　used　to　verify　the　correctness　of　the　enriched　NCDDAM　in　handling　contact　problems,　heat　conduction　problems　and　thermo-mechanical　coupling　problems.　The　effects　of　time　step　size　on　the　accuracy　of　NCDDAM　for　both　thermal　simulation　and　mechanical　simulation　are　investigated　detailly.　In　addition,　the　thermal　cracking　processes　of　rock　masses　for　both　experiment　and　engineering　scales　are　simulated　by　the　enriched　NCDDAM.　The　numerical　results　indicate　that　the　enriched　NCDDAM　is　a　powerful　tool　to　simulate　the　thermo-mechanical　coupling　processes　of　rock　masses.　©　2024　Elsevier　B.V.