In the analysis a 2 dimension element is combined with the boundary punition element(spring element) to calculate stress distribution at the periphery of a tunnel.