Zirconium tube is widely used in heat exchanger equipments in petrochemical process for significant corrosion resistance. The connection joint of tube-tubesheet is the weakest parts in a heat exchanger. The experiment and numerical analysis of different materials (zirconium tubes, titanium tubes and 16MnR tubesheets, 316L tubesheet) joints were performed in this paper. The expansion joints specimens were prepared at the pressure ranging from 28MPa to 38MPa. And pulling out test was performed from 20°C to 300°C. The finite element model of tube-tubesheet joint was established. The effect of expansion pressure, temperature and groove on the pulling out strength of joints was analyzed. Both the experiments and the finite element analysis show that the pull-out strength increases with the increasing expansion pressures. Working temperature also has a great effect on the connection strength of tube-to-tubesheet joints, especially for the zirconium and 316L joints, which have the most significant difference of thermal expansion coefficient between tube and tubesheet. The residual contacting pressure on the contact surface between tubes and the tubesheet is not uniformly distributed and two tightness bands are found near the surfaces of the tubesheet or at the two brinks of the groove on the tubesheet hole. Compared with the ungrooved joint, the residual contacting pressure on the tightness bands for the grooved joint is much higher, indicating a grooved joint has better tightness.