Abstarct: In many countries, tunnels in soils and rocks are constructed by using the New Austrian Tunneling Method NATM. Hakim highway twin tunnels with a length of about 1000 meters, has been constructed by the traditional method. Useful width of the way has been considered 14.5 m and allowed longitudinal slope has been taken 5%. The tunnel has been implemented with approximate width of 17.9 meters and a height of 13.9 meters by means of the NATM method in seven stages with side drifts (SD) methods. The aim of this thesis is the stability analysis of the highway tunnels, temporary support design, and control of structures displacement. In this study, stability analysis of tunnels was performed by using the Sakurai hazard warning level, numerical method and instrumentation results. Finally combination support system was proposed by using the equivalent sections method for the studied tunnel. Induced displacements are empirically controlled by adjusting the speed of excavation, distance between tunnel face and support, partial-face excavation and closure of invert. In this thesis, by using the finite difference method, the impact of round length, optimum excavation sequence, excavation patterns and interaction of two tunnels were investigated. According to the results, whatever the difference between the round length of the roof and bench in each drift be lesser then the difference between the displacements and trends will be lesser. Results show that the rapid closure of the supporting ring and the excavation stage of the central gallery (middle drift) are the most important factors in controlling tunnel deformations and surface settlements in soft ground tunneling. The rapid ring closure by adopting a larger excavation volume is more effective than adopting a smaller excavation volume for limiting the tunnel crown deformations and surface settlements in soft ground tunneling at shallow depth if the face stability is adequately maintained. Four different patterns were analyzed with different sequences. According to the results, reducing the free span distance greatly reduces maximum displacements. In this case remaining bench of the final displacement decreased of more than 70%. Closer full support activation distance reduces final displacements more than 50%. If two adjacent tunnels are constructed far apart such that their influence zones do not overlap, then the individual tunnels can be considered separately as single tunnels and analyzed as such. special attention was paid to the influence of lagging distance between the twin tunnel excavated faces (LT) and the load-transfer mechanism between the two tunnels. As well as the interaction between two tunnels becomes negligible for values of the pillar width to tunnel diameter ratio of W/D = 1.0 or greater. Finally, to control and bear the displacement around the tunnel by considering the different elements of the FLAC3D software and modeling them in the Situation of area, suitable support system was determined for tunnels, and Safety factors of support system was obtained by PCACOL software.