A micromechanics approach along with a matrix method is employed to evaluate effective elastic, piezoelectric, piezomagnetic and thermoelasticity constants of generally anisotropic multilayer composites. We assume that multiferroic composites have perfect and nonperfect 2–2 connectivity amongst the phases. To derive effective constants, we propose several simplified approaches, such as: matrix expressions of the effective properties of the layered composite. The derived formulas are successfully applied to find the explicit expressions of the effective properties for some practical cases, such as: a multiferroic composite composed of an orthotropic pyroelectric phase and an orthotropic magnetostrictive phase; a multiferroic composite composed of an orthotropic piezoelectric phase, an orthotropic pyromagnetic phase and an orthotropic elastic substrate. Present method can efficiently handle the most general type of multilayers with an arbitrary number of general anisotropic layers. Analytical expressions for effective material properties are derived. The effects of crystallographic orientations and volume fractions of constituting layers on the thermo-magneto-electro-elastic (TMEE) coefficients can be investigated for practically any important multilayer composite. Note that the present model can be used for the material optimization of these TMEE multilayer’s to obtain a maximum product property such as the magnetoelectric, pyroelectric, and pyromagnetic coefficients. The same method can be used to predict the effective properties of poroelastic multilayers.