Increasing demand in material and mechanical properties has led to production of complex composite structures. The composite structures, made of different materials, possess a variety of properties derived from each material. This has brought challenges in both analytical and numerical studies in thermal conduction which is of significant importance for thermoelastic problems. Therefore, a unified and effective approach would be desirable. The present study makes a first attempt to determining the analytical symplectic eigen solution for steady-state thermal conduction problem of multi-material crack. Based on the obtained symplectic eigen solution (including higher order expanding eigen solution terms), a new symplectic analytical singular element (SASE) for numerical modeling is constructed. It is concluded that composite structures composed of multi-material with complex geometric shapes can be modeled by the developed method, and the generalized flux intensity factors (GFIFs) can be solved accurately and efficiently. (C) 2016 Elsevier Ltd. All rights reserved.