This　paper　presents　a　finite　element　method　and　boundary　element　method　(FEM-BEM)　coupling　scheme　to　study　the　steady-state　2D　heat　transfer　problem　for　structures　with　multi-scale　features.　Based　on　the　geometric　information　or　material　properties,　the　considered　model　is　divided　into　FE　sub-domain　and　BE　sub-domain.　The　closed　BE　sub-domain　in　the　mixed　FEM-BEM　model　is　defined　as　a　finite　super-element　for　which　the　effective　stiffness　matrix　and　effective　fluxes　are　obtained　by　the　self　written　BE　code　and　hence　can　be　assembled　to　the　FE　systems　by　the　user　subroutine　(UEL)　provided　in　ABAQUS.　This　method　offers　a　number　of　key　improvements　compared　with　current　analysis　methods　available　for　multi-scale　structures.　These　improvements　are:　(i)　the　flexibility　and　convenience　by　using　the　powerful　pre/post-processing　in　ABAQUS;　(ii)　the　higher　accuracy　of　the　solution　over　the　classical　FEM;　and　(iii)　the　computational　cost　and　time　can　be　highly　reduced　by　using　the　coupling　method.　Numerical　results　are　compared　with　reference　solutions　and　demonstrate　the　accuracy　and　effectiveness　of　the　proposed　approach.　©　2022　Elsevier　Ltd