Based　on　a　consistent　couple　stress　elasticity　theory,　a　two-dimensional　(2D)　boundary　element　method　(BEM)　is　developed　to　solve　the　boundary　value　problems　in　the　couple　stress　elastic　materials.　The　displacement　boundary　integral　equations　(BIEs)　are　applied　to　compute　the　mechanical　quantities　on　the　boundary.　The　integral　representations　for　the　generalized　strains　and　stresses　at　internal　points　are　also　presented　together　with　the　explicit　expressions　for　the　high-order　fundamental　solutions.　The　collocation　method　is　adopted　for　the　spatial　discretization　of　the　BIEs.　Some　special　numerical　and　semi-analytical　techniques　are　provided　for　the　evaluation　of　weakly　and　strongly　singular　integrals.　The　sinh　transformation　is　used　to　deal　with　the　nearly　singular　integrals　for　near-boundary　internal　points　or　boundary　nodes　in　thin　structures.　An　error-controlled　integration　scheme　is　adopted　for　computing　the　non-singular　integrals　to　improve　the　accuracy.　Three　repre-sentative　2D　boundary　value　problems　in　the　consistent　couple　stress　elasticity　theory　are　studied　by　the　developed　BEM.　The　correctness　and　accuracy　of　the　present　BEM　results　are　verified　by　comparison　with　the　corresponding　analytical　results.　Then,　multiple　flat　elliptical　holes　in　a　plate　are　numerically　studied　to　show　the　interaction　effects　on　the　stress　concentration.