Steel Fiber Reinforced Concrete (SFRC) allows overcoming brittleness and weakness in tension, the main drawbacks of plain concrete. The goal of the present presentation is to present an ad-hoc numerical strategy to account for the contribution of the fibers in the simulation of the mechanical response of SFRC. In the model presented, the individual fibers immersed in the concrete bulk are accounted for in their actual location and orientation. The selected approach is based on the ideas introduced in the Immersed Boundary (IB) methods. These methods were developed to account for 1D (or 2D) solids immersed in 2D (or 3D) fluids. Here, the concrete bulk is playing the role of the fluid and the cloud of steel fibers is acting as the immerse boundary (that is a 1D structure in a 2D or 3D continuous). Thus, the philosophy of the IB methodology is used to couple the behavior of the two systems, the concrete bulk and the fiber cloud, precluding the need of matching finite element meshes. In the proposed approach, the meshes of the concrete bulk and fiber cloud are independent and the models are coupled imposing displacement compatibility and equilibrium of the two systems. The concrete bulk is modeled using any nonlinear model. The constitutive model for the fibers is designed to account for the complex interaction between fibers and concrete. The fiber models are based on previous investigations describing the concrete-fiber interaction and its dependence on the factors identified to be relevant: shape of the fiber (straight or hooked) and angle between the fiber and the crack plane. 3D examples with fibers distributed and oriented randomly are reproduced using the proposed approach.