The purpose of this paper is to establish a new constitutive model for functionally graded concrete (FGC) with high fiber content, considering randomly distributed cylindrical fibers and spherical, ellipsoidal, and cubic aggregates. Considering the fiber spacing coefficient and adjusting for the reinforcing effect of the fibers, the new correction factors include the fiber orientation coefficient, the fiber inclination impact coefficient, and the fiber spacing coefficient. Moreover, the interaction between the upper and lower interfaces of the FGC and the interaction between the fibers and the matrix are also considered. The interface adopts a cohesive bond-slip model to define interface plastic damage. The interaction between the fibers and the matrix is based on the rigid body-spring concept, where the matrix block is first discretized into rigid body-spring elements, which then interact with the fibers to produce bond-slip forces. During model validation, the proposed constitutive model was introduced into Abaqus using Python code. Finite element models for FGC, including the conventional uniaxial tensile model, dog-bone model, and four-point bending model, were established. The results indicate that for the FGC numerical computation model, the critical fiber volume fraction ρs ranges between 6 % and 7 %. However, the exact value should be determined based on specific practical conditions. When ρs is within the range of 0 %–6 %, the mechanical properties of FGC, such as peak stress, tensile strength, and toughness, increase with the increase in ρs. However, when ρs exceeds 6 %–7 %, the growth trend begins to slow down. In some cases, peak stress and tensile strength may decrease with an increase in ρs. At this point, the influence of the fiber spacing coefficient must be considered. Finally, by comparing with the experimental results, it is proved that the proposed constitutive model is feasible for studying the mechanical properties of FGC.