The paper presents general boundary element approach for analysis of thermomagnetoelectroelastic solids containing shell‐like electrically conducting inclusions with high magnetic permittivity. The latter are modeled with opened surfaces with certain boundary conditions on their faces. Rigid displacement and rotation, along with constant electric and magnetic potentials of inclusions are accounted for in these boundary conditions. Formulated boundary value problem is reduced to a system of singular boundary integral equations, which is solved numerically by the boundary element method. Special attention is paid to the field singularity at the front line of a shell‐like inclusion. Special shape functions are introduced, which account for this square‐root singularity and allow accurate determination of field intensity factors. Approaches for fast and accurate numerical evaluation of anisotropic thermomagnetoelectroelastic kernels are discussed, which are crucial in derivation of fast and precise boundary element approach. Numerical examples are presented.