The cracking-bridging capacity provided by steel fibers improves both the toughness and durability of concrete. several parameters such as geometrical and mechanical features of fibers and concrete affect the behavior of FRC, therefore, studing behaviors is necessary. In previous research, analysis of FRCs have been done by experimental investigations and also section analysis have been done by experimental and numerical investigations. In numerical investigations, the modeling have been proceed with continuum element. According to past researches, the effectiveness of distribution and orientation in concrete matrix play the main roles in FRC behavior. The objective of present thesis was to develop a new model for simulating fibers in real conditions and at the end of the thesis the FRC behavior have been studied under static loads. For modeling, the nonflexible fibers have been chosen and the proposed distributions in previous studies applied to the models and finally a governed distribution was proposed. In present thesis, a distribution algorithm is written in a numerical solver software (ABAQUS) which uses finite element method (FEM). This algorithm defined by a python code. For concrete element, damage concrete plasticity model is chosen and for steel element elasto-plastic behavior is considered. The influence of correlation between fiber and concrete is one of the most important facts on FRC behavior. So this fact is investigated related to whole previous suggested models and it has been shown that cohesive behavior contact is the most appropriate model for expressing continuity between fiber and concrete. Because application of mentioned model is very time consuming and needs high cupable hardware so another numerical method which is called embedded model is used. In sequence, a non-linear solution method is used to analysis the FRC models. Finally, two cylindrical and cubic specimens are chosen for numerical tests. The diameter and height of cylindrical specimen were 150 mm and 300 mm respectively, with 0.5 percent fiber volume fraction and also the cubic specimen has the length of 400 mm, height of 100 mm and width of 100 mm containing 1 percent fiber volume fraction. These two models have been investigated under static loads and their behaviors under bending and tensile were studied. Numerical predictions are compared to observations from experiments to demonstrate the ability to predict such FRC meaningfully. There is a good agreement between numerical and experimental results. By using the proposed numerical model in this thesis, the behavior of FRC could be investigated under various static and dynamic loads.