Glass fiber–reinforced polymer (GFRP) reinforcement provides an effective alternative to conventional steel in concrete structures due to its corrosion resistance. Nevertheless, the lower elastic modulus of GFRP necessitates careful consideration of serviceability behavior in GFRP-reinforced concrete members. This study presents a numerical sectional analysis model for predicting the flexural response and ultimate capacity of hybrid reinforced concrete beams incorporating embedded GFRP profiles in combination with either mild steel or GFRP reinforcement bars under monotonic static loading. The proposed model employs realistic nonlinear stress–strain relationships for concrete and steel, together with secant moduli of elasticity evaluated at different loading stages. Particular emphasis is placed on detailed stress distribution in flexural sections, including the contribution of tension stiffening in the post-cracking regime. The formulation integrates nonlinear constitutive material behavior with theoretical sectional equilibrium to evaluate the effective flexural secant stiffness. For practical serviceability assessment and to reduce dependence on complex analytical procedures, strain vectors and stiffness matrix components are derived using elasticity coefficients that reflect modulus degradation obtained from numerical analysis. The accuracy of the model is verified through comparison with experimental results, including ultimate flexural capacity and moment–deflection responses. Many crucial parameters were studied, such as the longitudinal reinforcement ratio, type of reinforcement, concrete compressive strength, position of the I-GFRP profile, and rotation of the I-GFRP profile. The results of this study demonstrated that both the longitudinal reinforcement ratio and the rotation of the I-GFRP profile have a significant influence on the ultimate load capacity and deflection behavior. The close agreement between numerical predictions and experimental observations demonstrates the reliability and applicability of the proposed model for structural engineering analysis and design.
