In this paper we propose a systematic methodology for improving computational efficiency of models describing the dynamics of materials with memory as part of multilayered structures, in particular in thermoelectric shape memory alloys actuators. The approach, based on a combination of the centre manifold technique and computer algebra, is systematic in a sense that it allows us to derive computational models with arbitrary order of accuracy with respect to certain small parameters. Computational results demonstrating the efficiency of the proposed methodology in reproducing the dynamics of austeniticmartensitic phase transformations upon thermoelectric cooling are presented.