In this paper, the tectonic stress field of Iran studied by a two-dimensional elastic finite element model (FEM) using inverse modeling technique. Using the named method, the relative magnitude of the tectonic forces and the regional parameters of the major blocks estimated. Geodetic and stress analysis results used for evaluating and validating the results of the developed inverse model. This developed model allows estimating the effect of different parameters of the stress pattern. The results denote that the direction of convergence plays an important role in distribution of the stress pattern, and the shape of the plate boundary is the main factor in controlling the shape of the distribution of the maximum compressional stress. The results also show that the magnitude of the maximum principal stress axes (δ1) reduce from the west and southwest to the north and northeast of Iran. The maximum stress magnitude is dominantly in the NW-Zagros, Azarbaijan, Qatar- Kazeron, Kope-Dagh, Zabol block and Lut block, while western part of Iran shows a thrusting mode stress state.