Reinforced concrete beam-column joint (BCJ) has a great role to maintain the safety of buildings under influence of seismic forces. During seismic actions, overturning moments are produced by lateral loads that are translated as axial loads in the columns. BCJ behavior changes with frame column axial load level. This paper presents BCJ behavior under varying column axial load from the tensile capacity to the compression capacity of column using finite element model (FE Model). A two dimensional (2D) FE Model taking into account material non-linearity was proposed. The proposed FE Model is verified with experimental results available in literature with varying column load level and failed in different modes. The comparison between experimental and numerical results indicates that the FE Model is able to simulate the performance of BCJs and is able to capture the different failure modes with acceptable accuracy. Moreover, the proposed FE Model is used to conduct a parametric study to investigate behavior of two specimens designed according to Eurocode recommendations under varying column load levels. An interaction diagram was introduced for each specimen expressed the behavior under varying column axial load using FE Model and equations of Eurocode. Difference in results appeared between equations of Eurocode and FE Model results in joint shear strength capacity.