The post-construction settlement of roadbed is directly related to the safety performance of highway. The prediction of creep behavior is of great significance to ensure the long-term safety of roadbed engineering. In finite element numerical simulation of creep behavior, the meshing of complex structural models requires a large amount of computational resources, as well as a large amount of computational time. In this paper, a new arbitrary polygonal hybrid stress element (PHSEM) considering creep is proposed for roadbed settlement problems in engineering. Based on the hybrid stress element method and the creep behavior of roadbed soil, the basic element format of PHSEM is derived. The analysis shows that this element introduces a higher-order stress field, which can better improve the calculation accuracy. The numerical calculation model of roadbed creep is further established, and the comparative analysis is carried out with MARC. The study shows that: the numerical simulation curves of PHSEM are well fitted, and the deviation is within a reasonable range, which verifies the validity of arbitrary polygonal stress hybrid element considering creep; the element can be a polygon with variable number of edges when meshing, which can well simulate the roadbed model with multiple media and realistically calculate the stress distribution in different parts; the establishment of PHSEM provides a new idea for creep research in practical engineering.