The paper presents novel boundary element technique for analysis of anisotropic thermomagnetoelectroelastic solids containing cracks and thin shell-like soft inclusions. Dual boundary integral equations of heat conduction and thermomagnetoelectroelasticity are derived, which do not contain volume integrals in the absence of distributed body heat and extended body forces. Models of 3D soft thermomagnetoelectroelastic thin inclusions are adopted. The issues on the boundary element solution of obtained equations are discussed. The efficient techniques for numerical evaluation of kernels and singular and hypersingular integrals are discussed. Nonlin-ear polynomial mappings are adopted for smoothing the integrand at the inclusion’s front, which is advantageous for accurate evaluation of field intensity factors. Special shape functions are introduced, which account for a square-root singularity of extended stress and heat flux at the inclusion’s front. Numerical example is presented.