Abstract: Thermoacoustic prime movers can steadily generate pressure oscillation without moving components on base of self-sustain effect. Here the details of simulation method for the thermoacoustic entire machine are presented, and two typical machines are calculated. A 4-port network method is used to build the transcendental equation with the frequency as a criterion to judge if the temperature distribution is correct. It is proved that the code can run robustly in wide range and output whatever results interested. Finally the calculated results are compared with the experiments of the real thermoacoustic prime mover. It shows that the error of frequency is no more than 1%, and the measurement of relative amplitude of pressure is about 70% of the calculated results. Finally some of influence factors are discussed.