Abstract: Until now a lot of investigations have been done on the thermoelasto-plastic analysis. Its genaral approach is to add up all the external forces on the structure to form a so-called load vector and to divide the element states into elastic, plastic and elastplastic elements, but the variation of the element states during the loading procedure and the loading path are always ingnored. In fact these two factors do affect the reliability of the analytic results. In order to improve the accuracy of the thermoelasto-plastic analysis, two algorithms are proposed which imitate the loading path and take the shifting of the element states into account. The external forces are expressed as the functions of time respectivley. Thereby the analysis procedure corresponds to the loading procedure. The element states are classified into six kinds according to the different stress-strain states. Thus a program system is developed. Two examples are given to show the difference between the analysis results under different loading path.