In engineering, it is quite often to have applications of the heat transfer of conduction having domain heat generation present inside. The paper aims to present boundary element formulations for analyzing the three-dimensional heat-conduction in dissimilarly bonded anisotropic media involving quadratic volume heat sources. In this paper, the additional volume integral present in the boundary integral equation is analytically transformed to the boundary surface for the volume heat sources modeled by quadratic functions. The technique of domain-mapping is employed to treat the three-dimensional anisotropic heat conduction in multiply adjoined media with proper interfacial conditions provided. For showing our successful implementation, a few example cases are analyzed with verification of independent analyses by the finite element method.