Virtual Node (VN) method is widely adopted to handle satellite network topological dynamics. However, conventional VN method is insufficient when earth rotation and inter-plane phase difference are considered. An improved VN method based on Celestial Sphere Division is proposed to overcome the defects of the conventional method. An optimized inter-satellite link connecting mode is derived to achieve maximal available links. The optimal VN division solution and addressing scheme are designed to generate a nearly static virtual network and solve the asynchronous switches caused by inter-plane phase difference. Comparison results demonstrate the advantages of proposed method.