A new mathematical model of thermoelectro viscoelasticity theory was constructed in the context of a new consideration of heat conduction with fractional order. The state space approach developed earlier by Ezzat was adopted for the solution of a one-dimensional problem in the presence of heat sources. The Laplace-transform technique was used. A numerical method was employed for the inversion of the Laplace transforms. According to the numerical results and their graphs, a conclusion about the new theory was constructed. Some comparisons are shown in figures to estimate the effect of the fractional order parameter on all of the studied fields. (C) 2011 Wiley Periodicals, Inc. J Appl Polym Sci, 2012