AbstractA 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. © 2011 Wiley Periodicals, Inc. J Appl Polym Sci, 2012