Three-dimensional thermo-mechanical stress analysis of composite joints with bi-adhesive bonding is presented using the full layerwise theory. Based on three-dimensional elasticity theory, sets of fully coupled governing differential equations are derived using the principle of minimum total potential energy and are simultaneously solved using the state space approach. Results show that bi-adhesive bonding substantially relieves the edge effects. Moreover, series of parametric studies reveal the nonlinear effects of bonding length ratio and the relative stiffness and coefficient of thermal expansion of the mid- and side-adhesives. It is also concluded that the optimum design of a bi-adhesive joint crucially depends on the loading/boundary conditions. (C) 2016 Elsevier Ltd. All rights reserved.