A numerical model that can simulate the nonlinear behavior of ultra high strength fiber-reinforced concrete (UHSFRC) structures subject to monotonic loadings is introduced. Since engineering material properties of UHSFRC are remarkably different from those of normal strength concrete and engineered cementitious composite, classification of the mechanical characteristics related to the biaxial behavior of UHSFRC, from the designation of the basic material properties such as the uniaxial stress-strain relationship of UHSFRC to consideration of the bond stress-slip between the reinforcement and surrounding concrete with fiber, is conducted in this paper in order to make possible accurate simulation of the cracking behavior in UHSFRC structures. Based on the concept of the equivalent uniaxial strain, constitutive relationships of UHSFRC are presented in the axes of orthotropy which coincide with the principal axes of the total strain and rotate according to the loading history. This paper introduces a criterion to simulate the tension-stiffening effect on the basis of the force equilibriums, compatibility conditions, and bond stress-slip relationship in an idealized axial member and its efficiency is validated by comparison with available experimental data. Finally, the applicability of the proposed numerical model is established through correlation studies between analytical and experimental results for idealized UHSFRC beams.