This article proposes an arbitrary polygonal hybrid stress element considering gravity. It derives an arbitrary polygonal hybrid stress element considering gravity alone for slope stability related engineering analysis. In the stability analysis of slopes, slope disasters caused by gravity erosion have recently become an urgent problem to be solved through engineering. The traditional finite element analysis of slope stability faces problems such as a large number of divided elements and slow calculation efficiency. By introducing high-order stress fields through stress hybridization elements, accurate results can be simulated using a small number of elements. When dividing the mesh, most of the cell shapes are asymmetric, and the shape of the cell can be any polygon, which can simulate the geometric shape of complex slopes well and more accurately calculate the stress distribution in different parts, thus accurately simulating the stability situation in engineering. By comparing with the corresponding commercial software MARC 2020, the effectiveness and efficiency of the element were verified. It has been proven that any polygonal hybrid stress element has the advantage of flexible mesh division, which can obtain high-order stress fields and stress concentration phenomena with fewer elements. Applying this element to practical problems of slopes in engineering has also achieved good calculation results.