[Display omitted] • A Quasi-Static FE based CDP model for analysis of GFRP-RC deep beams is proposed. • The impacts of main design variables on the behavior of GFRP-RC deep beams are assessed. • The suitability of five international design codes is evaluated. • A new design formula for the ultimate concrete shear contribution of the GFRP-RC deep beams is proposed. The main objective of this article is to propose a new, more robust, and comprehensive design equation for Glass-Fiber-Reinforced-Polymer (GFRP) reinforced concrete (RC) Deep Beams (DBs). The investigation is based on a numerical model using Quasi-Static Two-Dimensional FE simulation of GFRP-RC DBs. The numerical model is validated by calibrating four parameters (concrete dilatation angle ψ , concrete tensile strength f ct , fracture energy G f , mesh size) against data from ten full-scale experiments on GFRP-RC DBs available in the literature. Additionally, a numerical parametric study based on a full factorial design with 5 factors (and between 2 and 5 modalities) is conducted using the validated quasi-static FE model, which includes 360 large-scale GFRP-RC DBs. The impacts of primary design variables on the concrete contribution V c to the shear capacity of GFRP-RC DBs are assessed. The effects of various geometrical and mechanical variables, such as cross-sectional dimensions (width b × height H), concrete compressive strength grade (f c), GFRP properties (modulus E f , and strength f fu), GFRP reinforcement ratio (ρ f), and shear span to depth ratios (a / d f), are examined to understand the shear behavior of GFRP-RC DBs. The numerical simulations demonstrate that the proposed FE model accurately captures the behavior of GFRP-RC DBs. Building upon these findings, a new shear design equation is proposed for DBs, and its accuracy and robustness are compared with five established international design codes (ACI-440.1R-15, ISIS-Canada, CAN/CSA-S806-12, JSCE-1997, and BS-8110). The existing design codes significantly underestimate the concrete contribution to shear capacity, leading to uneconomical designs. In contrast, the new design equation provides a fair prediction of V c , with a standard error of 7.734 % and 20.84 % compared to the FE and experimental datasets, respectively. This represents a substantial improvement over the existing design guides. [ABSTRACT FROM AUTHOR]