We　introduce　a　coupled　finite　and　boundary　element　method　for　elastic-plastic　analysis　over　multiscale　electronic　packaging　structures.　Based　on　the　finite　element-boundary　element　(FE-BE)　coupling　algorithm,　an　automatic　implementation　procedure　for　the　coupling　of　the　ABAQUS　with　a　self-written　elastic　BE　code　is　introduced　for　elastic　problems.　In　the　mixed　finite　element　method　(FEM)-boundary　element　method　(BEM)　model,　the　effective　stiffness　and　effective　forces　at　the　interfacial　boundary　are　evaluated　by　the　self-written　BE　code.　Then,　the　obtained　effective　stiffness　and　effective　forces　are　assembled　to　the　global　FE　formulations　by　using　the　user　subroutine　(UEL)　in　ABAQUS.　Numerical　simulation　of　structures　with　plastic　deformation,　stress　concentration,　etc.　is　carried　out　by　using　FEM　theory.　The　boundary　element　method　is　used　for　linear　elastic　domains　with　large-scale　structures.　The　proposed　method　offers　several　key　improvements　compared　with　current　analysis　methods　available　for　multiscale　electronic　packaging　structures.　The　benefits　are:　(i)　the　powerful　pre-　and　postprocessing　of　ABAQUS;　(ii)　the　higher　accuracy　of　the　solution;　(iii)　the　computational　cost　and　time　can　be　reduced　by　using　the　scheme;　and　(iv)　solving　systems　with　infinite　extension　by　using　the　BEM　as　a　supplement.　Furthermore,　we　demonstrate　the　ability　of　the　proposed　approach　to　handle　multiscale　structures　in　electronic　packaging　problems.　　Copyright　©　2022　by　ASME.