This　paper　studied　the　transient　response　of　a　semi-permeable　crack　between　two　dissimilar　anisotropic　piezoelectric　layers　by　a　time-domain　BEM　with　a　sub-domain　technique　and　an　iterative　process　for　the　non-linear　crack-face　boundary　conditions.　The　present　time-domain　BEM　uses　a　quadrature　formula　for　the　temporal　discretization　to　approximate　the　convolution　integrals　and　a　collocation　method　for　the　spatial　discretization.　A　universal　matrix-form　displacement　extrapolation　formula　and　its　explicit　formula　are　used　to　determine　the　dynamic　intensity　factors.　Several　examples　are　presented　and　discussed　to　show　the　effects　of　the　electrical　crack-face　boundary　condition　on　dynamic　intensity　factors.　(C)　2013　Elsevier　Ltd.　All　rights　reserved.