A variational model is developed to exactly determine both stress and displacement fields at free edges of general symmetric composite laminate strips under thermomechanical loads. By partitioning the total stresses in a composite with free edges into unperturbed (without free edge) and perturbation stresses and using the minimum complementary energy principle, the optimal stress and displacement fields are derived that exactly satisfy equilibrium, compatibility, boundary and continuity conditions. The paper extends the theory of stress-based variational stress transfer so that the effects of both applied traction and displacement loads can be considered. It has been also shown how the displacement field can be determined for a stress-based variational approach.The results are compared to refined finite element results. The superiority of the developed method over finite element method, both in terms of accuracy and computational efficiency, is discussed.The method is also applicable to thin-ply laminates and is computationally efficient.