Nemerical algorithms were developed to analyze the behavior of the double lining as well as ground mass separately or simultaneously. A lining interface element was especially developed, verified and applied to the study on the coupled interaction of shotcrete and the concrete lining. It could be known fro parameter studys on double lining support systems that as the contact surface between shotcrete and concrete lining was rougher, the more decreased bearing capacity against the cracking of the system. If the thickness of the shotcrete increased, the bearing capacity of the double lining also increased linearly with the thickness. If the thickness of the concrete lining increased, the bearing capacity of the double lining had the relationship of the characteristic S-shape of a sigmoid function with the thickness. When the thickness increased over a given value, it was not useful to increase more the thickness because bearing capacity had no remarkable change. It could be concluded that the behavior of the shotcrete and concrete lining was generally reversed before and after the ratio of horizontal to vertical earth preassure of 1.0 and 0.5 respectively. Therefore, we could guess that the movement which two shotcrete and concrete lining deflect toward each other around the crown caused a friction between two linings and thus this disadvantageous effect could contribute to reducing the bearing capacity against the cracking.