Abstract: A complex variable method is presented of stress and displacement problems for non-circular cold region tunnels. The complex variable method considers the non-circle tunnel support, the frost heaving circle and the unfrozen surrounding rock at the same time, making the problem from the single connected domain to multi connected domain problem, and can not directly apply the classical complex function theory to solve the problem of the stress and displacement of the non-circular tunnel. The analytical formula in Zeta plane orthogonal curvilinear coordinate system of non-circular tunnel lining to frozen surrounding rock to unfrozen surrounding rock system is derived by using classical complex variable method and continuity condition with power series and conformal transformation. Then, the frost heaving force and frost heaving deformation of the non-circular cold region tunnel are obtained by conformal transformation. This complex variable method is applied to Zhegushan tunnel research, and the analytical solution of Zhegushan tunnel openning frost heave stress and frost heave displacement is obtained. Comparing the analytical solution with the numerical solution to verify the correctness of the analytical solution. It can be seen from the results that frost heaving force has an obvious influence on the lining, and the additional circumferential stresses caused by frost heaving at the vault, arch foot and arch bottom are significantly increased. The additional normal stresses caused by frost heaving at the vault and both sides of arch shoulder are significantly increased. Because of the geometric configurations of lining, the frost heaving deformation is uneven, which leads to uneven distribution of frost heaving forces. Compared with previous studies, which consider the frost heaving ring as a single ring, the complex variable method is more practical. And the research results were expected to provide a referenced basis for elasto-plastic analysis of cold-region tunnel.