The　dynamic　crack　propagation　in　composites　with　inclusions　is　analyzed　by　a　time-domain　boundary　element　method　(BEM)　in　conjunction　with　the　sub-region　technique.　The　crack-growth　direction　and　the　crack-tip　instantaneous　velocity　are　determined　by　the　maximum　circumferential　stress　criterion.　The　instantaneous　velocity　is　well　smoothed　by　a　bisection　technique.　New　crack-tip　elements　of　inconstant　length　are　added　to　the　active　crack-tip　to　simulate　the　fast　crack　growth.　Square　root　shape　functions　are　adopted　as　to　describe　the　proper　asymptotic　behavior　in　the　vicinity　of　the　crack-tips.　The　computation　time　for　the　dynamic　fracture　problems　in　composites　with　multiple　inclusions　is　reduced　by　a　numerical　method.　The　influences　of　the　inclusions　on　the　shielding　ratio,　the　crack　growing　path　and　the　crack-tip　instantaneous　speed　are　well　investigated.　Crown　Copyright　(C)　2011　Published　by　Elsevier　Ltd.　All　rights　reserved.