Abstract: In this paper a numerical modeling of the stress and strain fields as well as the displacement rates of the Tanlu fault zone in the principal periods of its tectonic evolution is made using the Finit Element Method for a Newtonian liquid. The causative dynamic origin, the boundary forces and the effect of rock properties on the stress fields are discussed and analysed.This is aimed at further explaining that the repeated tectonic activity of the zone in its evolution is closely related to the motion of the North China platform.Lower Yangtze massif, the Pasific plate and the Philippine plate. Meanwhile, the geothermal anomaly accompaning the tectonic movement may also constitute one of the important dynamic factors to form the respective stress field, especially the tensile one of the fault zone.