This study presents a finite element investigation of steel fiber reinforced concrete beam-column joints under cyclic loading with the variation of axial load. The aim of the study is to investigate the influence of axial load variations on the seismic behavior of steel fiber-reinforced concrete (SFRC) beam-column joints. Nonlinear finite element analysis with a damaged plasticity model in ABAQUS/Standard is adopted. The finite element model is verified using experimental results conducted by other researchers. Six SFRC specimens with different column axial load ratios and a 2% volume fraction of steel fiber were simulated under reversed cyclic loading. The parameters investigated are maximum load-carrying capacity, stiffness degradation, energy dissipation and failure mode. The results indicated that an increase of column axial load has a valuable influence to delay the initiation of cracks and damage accumulation, slightly improvement of the joint stiffness and improves the energy dissipation of joints at the initial stage of loading. Moreover, when the axial load level increases up to 50% of the column capacity, no cracks observed in the joint area and no change in the maximum load-carrying capacity. However, when the axial load level of more than 50% of the column capacity, the cyclic stiffness decreased slightly due to the deterioration caused by crushing of concrete in column. Thus, the results revealed that the increase of column axial load improves the confinement of steel fiber reinforced concrete beam-column joints, however, a threshold limit could be required.