Abstract: Using Kelvin's solutions and the concepts of finite-part integral,a set of hypersingular integral equations to solve the inclusion problems in two dimension elasticity is derived,and its numerical method is then proposed by combining the finite-part integral method with the boundary element method.Finally,several examples are carried out,and the numerical results of the stress intensity factors are obtained.