Performance Prediction of GFRP-Reinforced Concrete Deep Beams Containing a Web Opening in the Shear Span

: This study aimed to investigate the nonlinear structural behavior of concrete deep beams internally reinforced with glass fiber-reinforced polymer (GFRP) reinforcing bars and containing a web opening of various sizes and locations within the shear span. Three-dimensional (3D) numerical simulation models were developed for large-scale GFRP-reinforced concrete deep beams (300 mm × 1200 mm × 5000 mm) with a shear span-to-depth ratio ( a / h ) of 1.04. Predictions of the numerical models were validated against published experimental data. A parametric study was conducted to examine the effect of varying the opening size and location on the shear response. Results of the numerical analysis indicated that the strength of the deep beam models with an opening in the middle of the shear span decreased with an increase in either the opening width or height. The rate of the strength reduction caused by increasing the opening height was, however, more significant than that produced by increasing the opening width. Placing a web opening in the compression zone close to the load plate was very detrimental to the beam strength. Conversely, a negligible strength reduction was recorded when the web opening was placed in the tension side above the flexural reinforcement and away from the natural load path. Data of the parametric study were utilized to introduce simplified analytical formulas capable of predicting the shear capacity of GFRP-reinforced concrete deep beams with a web opening in the shear span.


Introduction
Discontinuity regions (D-regions) are formed in reinforced concrete (RC) beams due to statical or geometric discontinuities [1,2].The former are regions near concentrated loads and support reactions, whereas the latter are regions adjacent to web openings or abrupt changes in cross-section [1,2].Concrete deep beams, i.e., a shear span-to-depth ratio (a/h) ≤ 2, with and without cutouts, are influenced by both statical and geometric discontinuities [1,2].The discontinuity in statical loading or geometry causes a complex flow of internal stresses and nonlinear distribution of longitudinal strains within the cross-section.As such, the traditional beam theory (Bernoulli hypothesis) applied in the analysis of conventional RC structural members (B-regions), is not valid for the analysis of D-regions [1,2].Traditionally, D-regions in structural concrete members reinforced with steel reinforcing bars have been designed using dissimilar empirical equations that are not universally applicable [3,4].The problem becomes more challenging when conventional steel reinforcing bars are replaced by nonmetallic reinforcement such as fiber-reinforced polymer (FRP) composites.
The load-deflection response of FRP-reinforced deep beams is typically bilinear [11][12][13][14][15][16][17][18][19][20][21][22].This behavior is comparable to that reported for steel-reinforced deep beams [27][28][29].GFRPreinforced deep beams typically fail in a brittle manner because of the crushing in the diagonal compression strut, which is consistent with the failure mode of steel-reinforced deep beams [21,22,27].Deep beams reinforced with BFRP/GFRP reinforcing bars exhibited a greater midspan deflection and a lower shear capacity relative to those of their steelreinforced counterparts due to the lower elastic modulus, lower transverse strength, and lower dowel action of the BFRP/GFRP reinforcing bars [12,22,25].Increasing the section height reduced the normalized shear stress at ultimate load of BFRP/GFRP-reinforced deep beams [14,22].The reduction in the normalized shear stress at ultimate load of BFRPreinforced deep beams due to an increase in the section size was greater than that of their steel-reinforced counterparts [20,22].The bond performance of embedded FRP rebars in concrete was inferior to that of steel rebars [26].The bond quality of FRP rebars can be enhanced by using high-strength concrete materials to improve the load-bearing resistance of structural members such as RC deep beams [26].
Numerous studies have examined the shear behavior of FRP-reinforced concrete deep beams without web reinforcement [11][12][13][14][15][16][17][18][19][20].The shear capacity of solid deep beams reinforced with FRP bars improved with an increase in the value of a/h, the reinforcement ratio/modulus of elasticity of longitudinal FRP bars, and concrete compressive strength [11][12][13][14][15][16][17][18][19][20].An inverse linear correlation between the shear capacity and the cubic root of a/d was reported for solid deep beams reinforced with BFRP bars, where d is the beam effective depth [17,18].The increase in the shear capacity of FRP-reinforced deep beams due to decreasing a/h can be ascribed to an increase in the angle of inclination and/or width of the inclined concrete strut that governed the beam failure [15].Increasing the GFRP reinforcement ratio in concrete deep beams controlled the widening of shear cracks, and hence increased the shear strength [12,13].The enhancement in the shear strength of GFRP-reinforced concrete deep beams due to an increase in the concrete compressive strength may have a threshold, after which a minor additional increase in the shear capacity may be recorded [12].It is noteworthy that a significant increase in the longitudinal reinforcement ratio together with the concrete strength could be detrimental to the shear strength of GFRP-reinforced short beams without web reinforcement due to a change in the mode of failure from strut crushing to diagonal splitting [19].
Little information is available in the literature on the shear behavior of FRP-reinforced deep beams with web reinforcement and/or cutouts [21][22][23][24].The use of vertical GFRP stirrups improved the shear capacity of solid GFRP-reinforced large-scale deep beams by 20%, whereas the use of horizontal web reinforcement solely was detrimental to the shear strength because of the high tensile strains in the horizontal bars that induced deterioration and softening of the concrete in the diagonal strut [21].An increase in the value of a/h for deep beams with BFRP web reinforcement reduced the normalized ultimate shear strength but had an ignorable influence on the size effect [22].The presence of a web opening, with an opening width-to-shear span ratio of w o /a = 0.27 and an opening height-to-beam depth ratio of h o /h = 0.25, in the middle of the shear span of a large-scale GFRP-reinforced deep beam, with a/h = 1 and minimum shear reinforcement, reduced the shear capacity by 54%, relative to that of its solid counterpart [23].Providing additional GFRP shear reinforcement around the web opening increased the shear capacity by 20-56% compared to that of the deep beam without additional reinforcement around the openings [23].The use of diagonal GFRP reinforcement in the upper and lower chords of concrete deep beams with a web opening improved the shear capacity [24].Placing the diagonal GFRP reinforcement in the lower chord rather than in the upper chord was more effective in improving the shear capacity of GFRP-reinforced deep beams with a web opening [24].Despite the interesting findings and useful information offered in the aforementioned study [24], the beams did not include minimum GFRP shear reinforcement, which is impractical and not in compliance with the requirements of international design guides and standards (e.g., CSA S806 [5] and ACI 440.1R [6]).
Previous studies highlighted the lack of knowledge on the shear behavior of GFRPreinforced concrete deep beams with a web opening in the shear span [30].The interaction between the opening size and opening location in GFRP-reinforced deep beams was not investigated.This study aims to fill this gap through numerical analysis.Threedimensional (3D) simulation models capable of predicting the structural shear behavior of GFRP-reinforced concrete beams with D-regions were developed using ATENA ® software (v.5.9.0k) [31] and validated against experimental data from [23].A parametric study was conducted to investigate the effect of key parameters on the shear response of GFRP-reinforced deep beams.Refined analytical formulas were then developed to assist practitioners in predicting the shear capacity of GFRP-reinforced deep beams with a web opening.
The findings of the present study are anticipated to assist practitioners and researchers in designing concrete D-regions reinforced with GFRP bars.The outcomes of the study are anticipated to advance the development of design guidelines and standards on reinforcing concrete structures with GFRP bars.The widespread use of GFRP reinforcing bars instead of the conventional steel reinforcement in construction would reduce repair cycles and operational costs, thus leading to long-lasting sustainable structures.

Model Development 2.1. Geometry of the Beam Models
Details of the reinforcement of the deep beam models used for verification against the experimental data from [23], DB-S, DB-O1, DB-O2, and DB-O3, are shown in Figure 1a-d, respectively.The deep beam models were 5000 mm long, 300 mm wide, and 1200 mm deep.The effective span of the beams was 3000 mm, the shear span was 1250 mm, and the shear span-to-depth ratio (a/h) was 1.04.One deep beam model was solid (DB-S), whereas the other three models (DB-O1, DB-O2, and DB-O3) had a web opening in the middle of the shear span with a width (w o ) of 340 mm and height (h o ) of 304 mm, which corresponded to an opening width-to-shear span ratio (w o /a) of 0.27 and an opening height-to-beam depth ratio (h o /h) of 0.25.Figure 1a shows the location of monitoring points on the GFRP reinforcing bars of the solid deep beam model, whereas Figure 1b shows the location of monitoring points on the GFRP reinforcing bars of a deep beam model with an opening.For clarity purposes the monitoring points are not shown in Figure 1c,d

Constitutive Laws of Materials
Figure 2a,b show the concrete compressive hardening and softening, respectively, whereas the tensile softening law is shown in Figure 2c.The uniaxial constitutive laws of concrete in compression and tension begin with a linear relationship, with a slope equal to the concrete modulus of elasticity (Ec).The compressive hardening phase starts at a stress value of f'co = 2.1ft, where ft is the uniaxial tensile strength, and ends at a peak stress, f'c and a corresponding plastic strain, εcp.The values of ft and εcp are generated automatically by the software based on the concrete compressive strength [31].The post peak compressive stress decreases linearly with an increase in the displacement (wc) throughout the length scale (Lc).The value of wc can be calculated at any plastic strain εp through the relationship wc = (εp − εcp)Lc.The default value of the displacement at zero stress (wd) is 0.5 mm [31].The tensile softening law decreases exponentially from the peak tensile stress, ft, function of the crack opening displacement (wt) that is calculated by the product of the fracturing strain (εcf) and the length scale (Lt) [31].The crack opening at the complete release of stress (wtc) is determined automatically by the software based on ft and the specific fracture energy of the material (Gf).The cylindrical compressive strength (fc') of the concrete was 37 MPa, except for DB-O3, which had an fc' value of 45 MPa.Key input parameters of the concrete material used in the analysis are provided in Table 1.

Constitutive Laws of Materials
Figure 2a,b show the concrete compressive hardening and softening, respectively, whereas the tensile softening law is shown in Figure 2c.The uniaxial constitutive laws of concrete in compression and tension begin with a linear relationship, with a slope equal to the concrete modulus of elasticity (E c ).The compressive hardening phase starts at a stress value of f ′ co = 2.1f t , where f t is the uniaxial tensile strength, and ends at a peak stress, f ′ c and a corresponding plastic strain, ε cp .The values of f t and ε cp are generated automatically by the software based on the concrete compressive strength [31].The post peak compressive stress decreases linearly with an increase in the displacement (w c ) throughout the length scale (L c ).The value of w c can be calculated at any plastic strain ε p through the relationship w c = (ε p − ε cp )L c .The default value of the displacement at zero stress (w d ) is 0.5 mm [31].The tensile softening law decreases exponentially from the peak tensile stress, f t , function of the crack opening displacement (w t ) that is calculated by the product of the fracturing strain (ε cf ) and the length scale (L t ) [31].The crack opening at the complete release of stress (w tc ) is determined automatically by the software based on f t and the specific fracture energy of the material (G f ).The cylindrical compressive strength (f c ′ ) of the concrete was 37 MPa, except for DB-O3, which had an f c ′ value of 45 MPa.Key input parameters of the concrete material used in the analysis are provided in Table 1.
The behavior of GFRP reinforcing bars was modeled as linear elastic until failure.The vertical GFRP stirrups having a diameter of 12 mm (No. 12) had a cross-sectional area of 127 mm 2 , an elastic modulus of 50.0 GPa, and a tensile strength of 1019 MPa at the straight portion [23].It is noteworthy that the tensile strength of the bent portions of the vertical stirrups (459 MPa) was 45% of that of straight portions (1019 MPa), as reported in Reference [23].The strength reduction that typically occurs in the bent portions of FRP bars is attributed to the bending of fibers and the stress concentration that happens at the beginning of the bent portion, inducing premature failure at this location [5,6,32,33].As such, the vertical GFRP stirrups were divided into segments so that a reduced tensile strength can be assigned to the bent portions of the stirrups.The horizontal GFRP bars with a diameter of 16 mm (No. 16) had a cross-sectional area of 198 mm 2 , a tensile strength of 1184 MPa, and an elastic modulus of 62.6 GPa [14].The flexural GFRP bars in the tension side having a diameter of 25 mm (No. 25) had a cross-sectional area of 507 mm 2 , a tensile strength of 1000 MPa, and an elastic modulus of 66.4 GPa [14].A perfect bond was assumed between the GFRP reinforcing bars and the surrounding concrete.The validity of the perfect bond assumption between the GFRP bars and the concrete and its adequacy to produce a reasonable prediction for the load capacity of GFRP-reinforced concrete beams was reported in previous studies [34,35].The comparative analysis between predictions of the numerical models developed in the present study and the corresponding published experimental data, presented in Section 3 of this manuscript, further verifies the validity of this assumption.The numerical models developed and verified in the present study can be used in future research to compare the behavior of GFRP-reinforced concrete beams with a web opening to other similar beams reinforced with conventional steel bars with various yielding strengths.It is noteworthy that the contribution of the GFRP reinforcing bars around the opening is mostly provided by the vertical stirrups around the opening and the continuous horizontal bars in the top and bottom chords located above and below the opening, respectively.The horizontal bars intercepted by the opening are stopped near the surface of the opening.In most cases, these bars are not intercepted by the critical shear crack causing failure.As such, the contribution of the horizontal bars cut by the opening to the shear resistance provided by the reinforcement may be considered not critical relative to that of other GFRP reinforcement crossed by the critical load path in the upper and bottom chords.Future research should further investigate the effect of the anchorage length and bond condition of the horizontal bars cut by the opening on the response of the GFRP-reinforced deep beams.The behavior of GFRP reinforcing bars was modeled as linear elastic until failure.The vertical GFRP stirrups having a diameter of 12 mm (No. 12) had a cross-sectional area of 127 mm 2 , an elastic modulus of 50.0 GPa, and a tensile strength of 1019 MPa at the straight portion [23].It is noteworthy that the tensile strength of the bent portions of the vertical stirrups (459 MPa) was 45% of that of straight portions (1019 MPa), as reported in Reference [23].The strength reduction that typically occurs in the bent portions of FRP bars is attributed to the bending of fibers and the stress concentration that happens at the beginning of the bent portion, inducing premature failure at this location [5,6,32,33].As

Types of Element and Boundary Conditions
Solid 3D macro-elements were used to model the concrete and steel plates.The GFRP bars were modeled as one-dimensional elements embedded into the macro-elements of the concrete.One half of the beam was modeled to optimize the processing time.A mesh sensitivity analysis was carried out to determine the optimum mesh size, which was found to be 50 mm.The support plate was restricted from movement in the transverse and vertical directions.Surface support was also used to prevent the surface at the plane of symmetry from moving in the direction normal to the other symmetrical part of the beam.The beams were subjected to two concentrated loads placed symmetrically around the midspan.Concentrated loads represent loads of columns planted on deep beams, which is commonly used in practical settings in ground and/or transfer floors of high-rise buildings.It is noteworthy that practical applications may also include deep beams loaded through a slab.In such a case, the beams are subjected to uniform loads.Deep beams with uniform loads may fail at a higher shear force than those of similar beams with concentrated loads [27].The loading scheme adopted in the analysis was displacement-controlled at a prescribed displacement of 0.1 mm/step applied at the middle of the top surface of the loading plate.Monitoring points were used to measure the load and midspan deflection.The top and bottom steel plates were connected to the beam through fixed contacts.The standard Newton-Raphson iterative solution implemented in ATENA ® software (v.5.9.0k) was adopted [31].Figure 3a shows a typical numerical model for a solid deep beam, while Figure 3b shows a typical numerical model for a beam with a web opening in the middle of the shear span.General 3D views showing the flexural and web reinforcements of the numerical models are provided in Figure 4.
the reinforcement may be considered not critical relative to that of other GFRP reinforcement crossed by the critical load path in the upper and bottom chords.Future research should further investigate the effect of the anchorage length and bond condition of the horizontal bars cut by the opening on the response of the GFRP-reinforced deep beams.

Types of Element and Boundary Conditions
Solid 3D macro-elements were used to model the concrete and steel plates.The GFRP bars were modeled as one-dimensional elements embedded into the macro-elements of the concrete.One half of the beam was modeled to optimize the processing time.A mesh sensitivity analysis was carried out to determine the optimum mesh size, which was found to be 50 mm.The support plate was restricted from movement in the transverse and vertical directions.Surface support was also used to prevent the surface at the plane of symmetry from moving in the direction normal to the other symmetrical part of the beam.The beams were subjected to two concentrated loads placed symmetrically around the midspan.Concentrated loads represent loads of columns planted on deep beams, which is commonly used in practical settings in ground and/or transfer floors of high-rise buildings.It is noteworthy that practical applications may also include deep beams loaded through a slab.In such a case, the beams are subjected to uniform loads.Deep beams with uniform loads may fail at a higher shear force than those of similar beams with concentrated loads [27].The loading scheme adopted in the analysis was displacementcontrolled at a prescribed displacement of 0.1 mm/step applied at the middle of the top surface of the loading plate.Monitoring points were used to measure the load and midspan deflection.The top and bottom steel plates were connected to the beam through fixed contacts.The standard Newton-Raphson iterative solution implemented in ATENA ® software (v.5.9.0k) was adopted [31].Figure 3a shows a typical numerical model for a solid deep beam, while Figure 3b shows a typical numerical model for a beam with a web opening in the middle of the shear span.General 3D views showing the flexural and web reinforcements of the numerical models are provided in Figure 4.

Model Verification
The numerical models were verified against experimental data of four large-scale RC deep beams published in one of the top-tier journals [23].One of these beams was solid (DB-S), whereas the other three beams (DB-O1, DB-O2, and DB-O3) had a web opening in the shear span with different amounts of GFRP reinforcement around the opening.Model verification against experimental data of such large-scale beams may be considered adequate given the scope and complexity of the topic.It is noteworthy that no studies on large-scale GFRP-reinforced deep beams with a web opening are currently available in the literature except for that used in the model verification of the present work.Future research should consider studying the effect of using different types of reinforcing bars (e.g., carbon, glass, steel with different yield strengths) on the behavior of concrete deep beams with a web opening in the shear span.

Load-Deflection Response
The numerically predicted responses of DB-S, DB-O1, DB-O2, and DB-O3 are compared to experimental data from [23] in Figure 5.The numerically predicted responses started with a linear branch until flexural cracks initiated at approximately 500 kN.There was a deviation between the pre-cracking stiffness predicted numerically and that measured experimentally.This behavior is expected, particularly for such a large-scale deep beam.The actual deep beams could have had shrinkage cracks due to the large

Model Verification
The numerical models were verified against experimental data of four large-scale RC deep beams published in one of the top-tier journals [23].One of these beams was solid (DB-S), whereas the other three beams (DB-O1, DB-O2, and DB-O3) had a web opening in the shear span with different amounts of GFRP reinforcement around the opening.Model verification against experimental data of such large-scale beams may be considered adequate given the scope and complexity of the topic.It is noteworthy that no studies on large-scale GFRP-reinforced deep beams with a web opening are currently available in the literature except for that used in the model verification of the present work.Future research should consider studying the effect of using different types of reinforcing bars (e.g., carbon, glass, steel with different yield strengths) on the behavior of concrete deep beams with a web opening in the shear span.

Load-Deflection Response
The numerically predicted responses of DB-S, DB-O1, DB-O2, and DB-O3 are compared to experimental data from [23] in Figure 5.The numerically predicted responses started with a linear branch until flexural cracks initiated at approximately 500 kN.There was a deviation between the pre-cracking stiffness predicted numerically and that measured experimentally.This behavior is expected, particularly for such a large-scale deep beam.The actual deep beams could have had shrinkage cracks due to the large surface area prior to testing, which might have reduced its initial stiffness in the pre-cracking stage.The initiation of the flexural cracks caused a change in the load-deflection response.The beam models exhibited a quasilinear response in the post-cracking phase due to a progressive development of cracks during loading, which was in alignment with the experimental results.The post-cracking stiffness of the deep beams with openings was lower than that of the solid beam.Another change in slope was observed prior to reaching the ultimate load in DB-S, probably because of the development of a new major crack and/or localized rupture at the bent portion of a vertical stirrup.It is noteworthy that the post-cracking stiffness of all beams predicted numerically almost coincided with that measured experimentally.Table 2 presents the numerical and experimental ultimate loads along with the corresponding deflection capacities.The difference between the numerical and experimental ultimate loads was within a 12% error band.The deviation between the deflection capacity predicted numerically and that obtained from the tests did not exceed 15%, except for DB-O3 which showed a deviation of 28% in the deflection capacity.Generally, the difference between the numerical and experimental results is within the acceptable margin of error.and/or localized rupture at the bent portion of a vertical stirrup.It is noteworthy that the post-cracking stiffness of all beams predicted numerically almost coincided with that measured experimentally.Table 2 presents the numerical and experimental ultimate loads along with the corresponding deflection capacities.The difference between the numerical and experimental ultimate loads was within a 12% error band.The deviation between the deflection capacity predicted numerically and that obtained from the tests did not exceed 15%, except for DB-O3 which showed a deviation of 28% in the deflection capacity.Generally, the difference between the numerical and experimental results is within the acceptable margin of error.

Crack Pattern and Failure Mode
The crack patterns predicted numerically at different stages of loading for DB-S are shown in Figure 6, whereas Figure 7 shows those for a typical deep beam model with a web opening (DB-O2).Schematics of the crack patterns obtained from the experiments based on data obtained from [23] for DB-S and a typical deep beam with a web opening are shown in Figure 8a,b respectively.From Figure 6, it can be seen that the numerical prediction of crack development for DB-S indicated the initiation of flexural cracks prior to shear cracks.The experimental data from [23] also indicated the initiation of flexural cracks early at 18% of the ultimate load before the development of any shear cracks.Diagonal cracks were then formed in the shear span as the load progressed.Additional shear cracks were developed with an increase in the applied load defining the direction of a concrete diagonal strut between the load and support points.Eventually, DB-S failed by crushing of the diagonal concrete strut formed in the shear span, which was in alignment with the mode of failure reported in the literature [23], as shown in Figure 8a.cracks early at 18% of the ultimate load before the development of any shear cracks.Diagonal cracks were then formed in the shear span as the load progressed.Additional shear cracks were developed with an increase in the applied load defining the direction of a concrete diagonal strut between the load and support points.Eventually, DB-S failed by crushing of the diagonal concrete strut formed in the shear span, which was in alignment with the mode of failure reported in the literature [23], as shown in Figure 8a. to shear cracks.The experimental data from [23] also indicated the initiation of flexural cracks early at 18% of the ultimate load before the development of any shear cracks.Diagonal cracks were then formed in the shear span as the load progressed.Additional shear cracks were developed with an increase in the applied load defining the direction of a concrete diagonal strut between the load and support points.Eventually, DB-S failed by crushing of the diagonal concrete strut formed in the shear span, which was in alignment with the mode of failure reported in the literature [23], as shown in Figure 8a. Figure 7 shows the crack patterns predicted numerically at different stages of loading for a typical deep beam model with a web opening (DB-O2).Shear cracks initiated first at the opposite corners of the opening closer to the support and load plates.This behavior was in agreement with the experimental data from [23], which indicated the early formation of shear cracks at the opening corner at 17 to 24% of the ultimate load.A further increase in load resulted in the development of flexural cracks and propagation of these Figure 7 shows the crack patterns predicted numerically at different stages of loading for a typical deep beam model with a web opening (DB-O2).Shear cracks initiated first at the opposite corners of the opening closer to the support and load plates.This behavior was in agreement with the experimental data from [23], which indicated the early formation of shear cracks at the opening corner at 17 to 24% of the ultimate load.A further increase in load resulted in the development of flexural cracks and propagation of these shear cracks toward the load and support plates along with the development of flexural cracks.The shear cracks that developed earlier at the opening corners were then stabilized and other shear cracks then developed in the top and bottom chords.Eventually, the beam model with a web opening failed due to the formation of two independent major shear cracks that developed in the top and bottom chords.The major shear crack causing failure in the top chord connected the edge of the load plate to the top corner of the opening closer to the support plate.The major shear crack causing failure in the bottom chord connected the edge of the support plate and the bottom corner of the opening closer to the load point.
The crack patterns at the different stages of loading and the failure mechanism predicted numerically were in good agreement with the sequence of crack propagation and failure mode observed during the experimental tests reported in the literature [23], as shown in Figure 8b.The stress in this location was almost equal to that of the tensile strength of the bent portion of GFRP bars, indicating localized rupture of GFRP at this location.These numerical findings are consistent with the published experimental data [23], which indicated that crushing of the diagonal concrete strut in DB-S was accompanied by a localized rupture at the bent portion of the vertical stirrups.For DB-O1, the stress in the bent portions of the vertical GFRP stirrups at the sides of the opening almost reached the tensile strength of the bent portion of GFRP bars (Figure 9b).The experimental findings also indicated that rupture of the bent portions of the vertical GFRP stirrups near the load and support plates occurred at the ultimate load along with the crushing of concrete along the diagonal cracks that developed in the upper and lower chords [23].In agreement with the published data [14], the predicted stresses in the vertical stirrups for DB-O2 and DB-O3, having extra GFRP reinforcement around the opening, were below the tensile strength of the bent portion of GFRP bars (Figure 9c).These results verify the effectiveness of the extra GFRP bars in relieving the web reinforcement, which allowed the beam models to develop a higher load capacity.The stresses in the horizontal web reinforcement in DB-O2 and DB-O3 at peak load were also lower than those of DB-O1 (Figure 10).From Figure 11, it can be seen that the maximum stress in the flexural reinforcement at peak load for DB-S was 390 MPa (i.e., 39% of the tensile strength of straight GFRP bars).The models with a web opening failed at a lower load than that of DB-S, and hence, the maximum stress in the flexural reinforcement was lower at an average value of 265 MPa (i.e., 27% of the tensile strength of straight GFRP bars).Published experimental data [23] verified the nonoccurrence of the rupture of the flexural GFRP reinforcement in the tested deep beams.It is noteworthy that the negative stress in the top GFRP longitudinal reinforcement in DB-S at ultimate load was higher than that of their counterparts with a web opening because DB-S failed at a higher load.stress in the flexural reinforcement was lower at an average value of 265 MPa (i.e., 27% of the tensile strength of straight GFRP bars).Published experimental data [23] verified the nonoccurrence of the rupture of the flexural GFRP reinforcement in the tested deep beams.It is noteworthy that the negative stress in the top GFRP longitudinal reinforcement in DB-S at ultimate load was higher than that of their counterparts with a web opening because DB-S failed at a higher load.

Parametric Study
Following the validation of the developed numerical models, a parametric study was conducted to study the effect of different opening sizes and locations on the shear behavior of GFRP-reinforced deep beams with a web opening.The parametric study was divided into two phases.The first phase focused on the analysis of deep beams with a web opening in the middle of the shear span having different sizes.The second phase focused on the analysis of deep beams with a web opening at different locations with respect to the natural load path within the shear span, noting that the natural load path is the line with a web opening failed at a lower load than that of DB-S, and hence, the maximum stress in the flexural reinforcement was lower at an average value of 265 MPa (i.e., 27% of the tensile strength of straight GFRP bars).Published experimental data [23] verified the nonoccurrence of the rupture of the flexural GFRP reinforcement in the tested deep beams.It is noteworthy that the negative stress in the top GFRP longitudinal reinforcement in DB-S at ultimate load was higher than that of their counterparts with a web opening because DB-S failed at a higher load.

Parametric Study
Following the validation of the developed numerical models, a parametric study was conducted to study the effect of different opening sizes and locations on the shear behavior of GFRP-reinforced deep beams with a web opening.The parametric study was divided into two phases.The first phase focused on the analysis of deep beams with a web opening in the middle of the shear span having different sizes.The second phase focused on the analysis of deep beams with a web opening at different locations with respect to the natural load path within the shear span, noting that the natural load path is the line with a web opening failed at a lower load than that of DB-S, and hence, the maximum stress in the flexural reinforcement was lower at an average value of 265 MPa (i.e., 27% of the tensile strength of straight GFRP bars).Published experimental data [23] verified the nonoccurrence of the rupture of the flexural GFRP reinforcement in the tested deep beams.It is noteworthy that the negative stress in the top GFRP longitudinal reinforcement in DB-S at ultimate load was higher than that of their counterparts with a web opening because DB-S failed at a higher load.

Parametric Study
Following the validation of the developed numerical models, a parametric study was conducted to study the effect of different opening sizes and locations on the shear behavior of GFRP-reinforced deep beams with a web opening.The parametric study was divided into two phases.The first phase focused on the analysis of deep beams with a web opening in the middle of the shear span having different sizes.The second phase focused on the analysis of deep beams with a web opening at different locations with respect to the natural load path within the shear span, noting that the natural load path is the line

Parametric Study
Following the validation of the developed numerical models, a parametric study was conducted to study the effect of different opening sizes and locations on the shear behavior of GFRP-reinforced deep beams with a web opening.The parametric study was divided into two phases.The first phase focused on the analysis of deep beams with a web opening in the middle of the shear span having different sizes.The second phase focused on the analysis of deep beams with a web opening at different locations with respect to the natural load path within the shear span, noting that the natural load path is the line connecting the edges of the support plate and the load plate.The results of the parametric study include the load-deflection response, crack propagation, failure mechanism, and stresses in the GFRP reinforcement.

Deep Beam Models with Different Opening Sizes
Parameters of the deep beam models with different opening sizes are summarized in Table 3.The variables included values of w o /a (0.16, 0.27, and 0.32) and values of h o /h (0.17, 0.25, and 0.33).The deep beam models are labeled DB-WX-HY, where X denotes the value of w o /a and Y refers to h o /h.For instance, DB-W0.16     The installation of the web opening in the middle of the shear span reduced the shear cracking load and the post-cracking stiffness relative to those of the solid deep beam model.An increase in the opening size further compromised the response of the beam models.The reductions in the cracking load and post-cracking stiffness caused by increasing the opening height were more pronounced for the models with the greater w o /a of 0.32.The beam models with the web opening failed at a deflection significantly lower than that of the solid beam model DB-S.Table 4 presents the ultimate load and the deflection at ultimate load for the beam models.At the same w o /a of 0.16, strength reductions of 31, 37, and 49% were recorded for the deep beam models with h o /h values of 0.17, 0.25, and 0.33, respectively.The strength reduction caused by increasing the opening height was more pronounced at the greater w o /a values.At a w o /a of 0.27, strength reductions of 35, 42, and 52% were recorded for the deep beam models with h o /h values of 0.17, 0.25, and 0.33, respectively.At a w o /a of 0.32, strength reductions of 39, 47, and 61% were recorded for the deep beam models with h o /h values of 0.17, 0.25, and 0.33, respectively.Although the deflection capacity of the beam models with a web opening was significantly lower than that of the solid deep beam model, changing the opening size had a minor impact on the deflection capacity.The deflection capacities of the beam models with different opening heights were insignificantly different.The deflection capacities of the beam models with w o /a values of 0.16, 0.27, and 0.32 were 56, 64, and 68% of those of the solid beam model DB-S.The strength of the beam models with web openings was normalized to that of the solid beam model DB-S, and then plotted against the opening size in Figure 16.It is evident that the strength of the beam models decreased with an increase in either the opening width or height.The rate of the strength reduction tended to be significant when the value of w o /a increased from 0.27 and 0.32 (Figure 16a).Similarly, the strength degraded at a higher rate when the value of h o /h increased from 0.25 and 0.33 as shown in Figure 16b.The rate of the strength reduction caused by increasing the opening height (Figure 16b) was, however, more significant than that produced by increasing the opening width (Figure 16a).It is noteworthy that the strength of deep beams with a web opening in the middle of the shear span typically decreases when reducing the angle of inclination of the lower load path connecting the edge of the support plate and the bottom corner of the opening nearest to the load plate (angle of inclination of crack D in Figure 8b) [36,37].The increase in the opening height caused a greater reduction in the angle of inclination of the lower load path than that caused by increasing the opening width, which justifies the higher rate of strength reduction exhibited by the beam models with an increased opening height.wo/a values of 0.16 and 0.27, extensive shear cracks also developed along the upper load path in the top chord connecting the edge of the load plate to the opposite corner of the opening, causing failure along the upper load path simultaneously with that which occurred along the lower load path.It is noteworthy that the lower load path is vulnerable to failure before the upper load path because of the transverse strains caused by the stresses in the flexural tensile reinforcement at the bottom of the beam model.In contrast, the deep beam models with the opening width of 400 mm, having a respective wo/a of 0.32, failed due to the excessive widening of shear cracks at the opening corners closer to the load and support plates along with excessive rotation that occurred in the top and bottom chords.

GFRP Stresses
General 3D views of the stresses in the vertical GFRP stirrups predicted numerically for the beam models of this group are shown in Figures 20-22.For the beam models with the smaller wo/a value of 0.16, the maximum stresses in the vertical GFRP stirrups were in the straight portions of the stirrup near the corner of the web opening closer to the support plate (Figure 20).The maximum GFRP stress at peak load for DB-W0.16-H0.17with wo/a = 0.16 and ho/h = 0.17 was 48% of the tensile strength of straight GFRP bars.The respective value for the other two models with the greater ho/h values of 0.25 and 0.33 was, on average, 57% of the tensile strength of straight GFRP bars.For the beam models with wo/a = 0.27, the maximum stresses in the vertical GFRP stirrups were in the straight portions of the stirrup near both corners of the web opening closer to the support and load plates (Figure 21).The maximum GFRP stresses at peak load for the beam models with wo/a = 0.27 and ho/h values of 0.17, 0.25, and 0.33 were 42, 47, and 54% of the tensile strength of straight GFRP bars, respectively.The vertical GFRP bars at the sides of the opening exhibited higher stresses in the beam models with the greater opening height because increasing the opening height reduced the contribution of the concrete to the shear resistance, and hence, higher stresses were transferred/resisted by the GFRP stirrups.For the beam models with the greatest wo/a value of 0.32, the maximum stresses in the vertical

GFRP Stresses
General 3D views of the stresses in the vertical GFRP stirrups predicted numerically for the beam models of this group are shown in Figures 20-22.For the beam models with the smaller w o /a value of 0.16, the maximum stresses in the vertical GFRP stirrups were in the straight portions of the stirrup near the corner of the web opening closer to the support plate (Figure 20).The maximum GFRP stress at peak load for DB-W0.16-H0.17with w o /a = 0.16 and h o /h = 0.17 was 48% of the tensile strength of straight GFRP bars.The respective value for the other two models with the greater h o /h values of 0.25 and 0.33 was, on average, 57% of the tensile strength of straight GFRP bars.For the beam models with w o /a = 0.27, the maximum stresses in the vertical GFRP stirrups were in the straight portions of the stirrup near both corners of the web opening closer to the support and load plates (Figure 21).The maximum GFRP stresses at peak load for the beam models with w o /a = 0.27 and h o /h values of 0.17, 0.25, and 0.33 were 42, 47, and 54% of the tensile strength of straight GFRP bars, respectively.The vertical GFRP bars at the sides of the opening exhibited higher stresses in the beam models with the greater opening height because increasing the opening height reduced the contribution of the concrete to the shear resistance, and hence, higher stresses were transferred/resisted by the GFRP stirrups.For the beam models with the greatest w o /a value of 0.32, the maximum stresses in the vertical GFRP stirrups were in the straight portions of the stirrup near the corner of the web opening closer to the support plate, except for DB-W0.32-H0.17,which exhibited maximum GFRP stresses at both corners of the web opening (Figure 22).The maximum GFRP stresses at peak load for the beam models with a w o /a of 0.32 were in the range of 34 to 42% of the tensile strength of straight GFRP bars, and they also tended to increase with an increase in the opening height.In contrast, for the beam models with the same opening height and different opening widths, it can be seen that the stresses in the vertical GFRP stirrups at the sides of the opening decreases with an increase in the opening width.This behavior could be due to an increased transfer of shear stresses through the upper and lower chords for the beam models with the greater opening width.It is noteworthy that the stresses in the bottom portion of the vertical legs of the GFRP stirrup placed at the side of the opening closer to the support plate in most of the beam models with w o /a values of 0.16 and 0.27 were in the range of 90 to 100% of the tensile strength of the bent portion of GFRP bars for the beam models.These results implied the possible rupture of these vertical GFRP For the beam models with the greater w o /a of 0.32, the stresses in the vertical GFRP bars were below the rupture strength of the bent portions of the GFRP bars.
stirrups at the sides of the opening decreases with an increase in the opening width.This behavior could be due to an increased transfer of shear stresses through the upper and lower chords for the beam models with the greater opening width.It is noteworthy that the stresses in the bottom portion of the vertical legs of the GFRP stirrup placed at the side of the opening closer to the support plate in most of the beam models with wo/a values of 0.16 and 0.27 were in the range of 90 to 100% of the tensile strength of the bent portion of GFRP bars for the beam models.These results implied the possible rupture of these vertical GFRP stirrups at the bent portions in most of the beam models with wo/a values of 0.16 and 0.27.For the beam models with the greater wo/a of 0.32, the stresses in the vertical GFRP bars were below the rupture strength of the bent portions of the GFRP bars.height and different opening widths, it can be seen that the stresses in the vertical GFRP stirrups at the sides of the opening decreases with an increase in the opening width.This behavior could be due to an increased transfer of shear stresses through the upper and lower chords for the beam models with the greater opening width.It is noteworthy that the stresses in the bottom portion of the vertical legs of the GFRP stirrup placed at the side of the opening closer to the support plate in most of the beam models with wo/a values of 0.16 and 0.27 were in the range of 90 to 100% of the tensile strength of the bent portion of GFRP bars for the beam models.These results implied the possible rupture of these vertical GFRP stirrups at the bent portions in most of the beam models with wo/a values of 0.16 and 0.27.For the beam models with the greater wo/a of 0.32, the stresses in the vertical GFRP bars were below the rupture strength of the bent portions of the GFRP bars.None of the horizontal web reinforcing bars or the flexural reinforcement at the tension side in the deep beam models of this group reached their tensile strength.The stresses in the horizontal web reinforcement and tensile flexural reinforcement predicted numerically for the beam models of this group are shown in Figures 23-28.For the beam models with the smaller wo/a value of 0.16, the maximum stresses in the horizontal web reinforcement were on average 258 MPa, which corresponded to 22% of the tensile strength of the GFRP bars (Figure 23).The maximum stress in their tensile flexural reinforcement ranged from 190 to 249 MPa, which corresponds to 19 to 25% of the tensile strength of the GFRP bars (Figure 24).The beam models with the wo/a value of 0.27 exhibited maximum stresses of 240 to 348 MPa (i.27 and 28, respectively.These results indicate that for the beam models with the web opening in the middle of the shear span, varying the opening size had an insignificant effect on the maximum stresses of the horizontal web reinforcement and tensile flexural reinforcement recorded at peak load.Generally, the maximum stresses recorded in the horizontal web reinforcement were in the range of 21 to 38% of the GFRP tensile strength, whereas those of the tensile flexural reinforcement were in the range of 19 to 25% of the GFRP tensile strength. reinforcement ranged from 190 to 249 MPa, which corresponds to 19 to 25% of the tensile strength of the GFRP bars (Figure 24).The beam models with the wo/a value of 0.27 exhibited maximum stresses of 240 to 348 MPa (i.e., 21 to 30% of the tensile strength of the GFRP bars) in the horizontal web reinforcement and 200 to 236 MPa (i.e., 20 to 24% of the tensile strength of the GFRP bars) in the tensile flexural reinforcement as shown in Figures 25 and 26, respectively.The beam models with the wo/a value of 0.32 exhibited maximum stresses of 324 to 382 MPa (i.e., 27 to 38% of the tensile strength of the bars) in the horizontal web reinforcement and 169 to 226 MPa (i.e., 17 to 23% of the tensile strength of the GFRP bars) in the tensile GFRP flexural reinforcement as shown in Figures 27 and 28, respectively.These results indicate that for the beam models with the web opening in the middle of the shear span, varying the opening size had an insignificant effect on the maximum stresses of the horizontal web reinforcement and tensile flexural reinforcement recorded at peak load.Generally, the maximum stresses recorded in the horizontal web reinforcement were in the range of 21 to 38% of the GFRP tensile strength, whereas those of the tensile flexural reinforcement were in the range of 19 to 25% of the GFRP tensile strength.

Deep Beam Models with Different Opening Locations
The parameters of the deep beam models with different opening locations are summarized in Table 5.The deep beam models of this group had the same web opening size of wo = 340 mm and ho = 304 mm (i.e., wo/a = 0.27 and ho/h = 0.25).The primary variable of this group was the location of the web opening.The center of the web opening was located at a distance xo from the face of the support plate and yo from the bottom soffit of the beam model.The locations of the center of the opening in the horizontal and vertical directions were normalized to the clear shear span (Xc) and beam depth (h), respectively.Values of the term xo/Xc were 0.25, 0.50, and 0.75, whereas those of the term yo/h were 0.

Deep Beam Models with Different Opening Locations
The parameters of the deep beam models with different opening locations are summarized in Table 5

Model Designation
Opening Location (mm)

Load-Deflection Response
Figure 32a-c show the load-deflection responses of the deep beam models with a web opening having an xo/Xc of 0.25, 0.50, and 0.75 respectively.The response of the solid deep beam model DB-S is included in the figures for the purpose of comparison.Figure 32a shows that installation of the web opening near the support at an xo/Xc of 0.25 seriously impaired the response of the beam models.The detrimental effect of the web opening was more pronounced for the beam models DB-X0.25-Y0.33 and DB-X0.25-Y0.50 with the lower yo/h values of 0.33 and 0.50, respectively, because the opening in these two models fully interrupted the natural load path.The web opening in the beam model DB-X0.25-Y0.75 did not interrupt the natural load path, and hence, its post-cracking stiffness was almost the same as that of DB-S.The beam model DB-X0.25-Y0.75failed, however, at lower load and deflection capacities than those of DB-S.Although the web opening in DB-X0.25-Y0.75 did not interrupt the natural load path, it could have affected the width of the middle portion of a bo?le-shaped strut that could have formed in the shear span.In addition, placing the web opening in the compression zone above the natural load path hindered the transfer of the load to the support through other struts formed above the

Load-Deflection Response
Figure 32a-c show the load-deflection responses of the deep beam models with a web opening having an x o /X c of 0.25, 0.50, and 0.75 respectively.The response of the solid deep beam model DB-S is included in the figures for the purpose of comparison.Figure 32a shows that installation of the web opening near the support at an x o /X c of 0.25 seriously impaired the response of the beam models.The detrimental effect of the web opening was more pronounced for the beam models DB-X0.25-Y0.33 and DB-X0.25-Y0.50 with the lower y o /h values of 0.33 and 0.50, respectively, because the opening in these two models fully interrupted the natural load path.The web opening in the beam model DB-X0.25-Y0.75 did not interrupt the natural load path, and hence, its post-cracking stiffness was almost the same as that of DB-S.The beam model DB-X0.25-Y0.75failed, however, at lower load and deflection capacities than those of DB-S.Although the web opening in DB-X0.25-Y0.75 did not interrupt the natural load path, it could have affected the width of the middle portion of a bottle-shaped strut that could have formed in the shear span.In addition, placing the web opening in the compression zone above the natural load path hindered the transfer of the load to the support through other struts formed above the natural load path.Figure 32b shows that the behavior of the beam models with the web opening installed at an x o /X c of 0.50 was affected by the location of the center of the opening with respect to the bottom face of the beam (y o ).Although the behavior of all models was inferior to that of the solid deep beam model DB-S, the degradation in the behavior intensified with an increase in the distance between the center of the opening and the bottom face of the beam.These results indicate that placing the web opening above the natural load path was more detrimental to the structural behavior than placing it below the natural load path.The presence of a web opening in the compression zone above the natural load path interrupted the formation of struts in this region, and hence, compromised the structural response.Figure 32c shows that the behavior of the beam models with the web opening installed at x o /X c of 0.   Table 6 presents the ultimate load and the deflection capacity for the beam models of this group.For the beam models having a web opening with an x o /X c of 0.25, strength reductions of 43 and 49% were recorded for DB-X0.25-Y0.33 and DB-X0.25-Y0.50having y o /h values of 0.33 and 0.50, respectively.It is noteworthy that the web opening in DB-X0.25-Y0.33 and DB-X0.25-Y0.50fully interrupted the natural load path.When the opening was pushed upward away from the natural load path in DB-X0.25-Y0.75, a lower strength reduction of 34% was recorded.The deflection capacity of the beam models containing a web opening with an x o /X c of 0.25 was, on average, 57% of that of DB-S.For the beam models having a web opening with an x o /X c of 0.50, strength reductions of 24, 42, and 52% were recorded for DB-X0.50-Y0.33,DB-X0.50-Y0.50, and DB-X0.50-Y0.75having y o /h values of 0.33, 0.50, and 0.75, respectively.These results indicate that the strength reduction for the beam models with an x o /X c = 0.5 was intensified with an increase in the distance between the center of the opening and the bottom face of the beam, noting that the opening in these three models interrupted the natural load path.The reduction in the ultimate load was accompanied by a reduction in the corresponding deflection capacity of 35 to 61%.For the beam models having a web opening with an x o /X c of 0.75, strength reductions of 7, 28, and 56% were recorded for DB-X0.75-Y0.33,DB-X0.75-Y0.50, and DB-X0.75-Y0.75having y o /h values of 0.33, 0.50, and 0.75, respectively.The beam model DB-X0.75-Y0.33exhibited a negligible strength reduction of 7% because the web opening was in the tension side and did not interrupt the natural load path.In contrast, the beam model DB-X0.75-Y0.75exhibited a significant strength reduction of 56% because the web opening was in the compression zone close to the load plate and fully interrupted the natural load path.The deflection capacity of DB-X0.75-Y0.33 with a y o /h of 0.33 was almost the same as that of DB-S.Conversely, the beam model DB-X0.75-Y0.75 with a y o /h of 0.75 exhibited a 45% reduction in the deflection capacity relative to that of DB-S.The relationships between the strength reduction caused by the web opening and the opening location are provided in Figure 33.For the beam models with a y o /h of 0.33 and 0.50, the strength increased with an increase in the horizontal distance measured from the face of the support within the shear span.The strength of the beam models with a y o /h of 0.33 tended to be higher than that of their counterparts with a y o /h of 0.50, and this behavior was more evident with an increase in the distance from the face of the support plate.The opening in the beam models with a y o /h = 0.33 was just above the longitudinal reinforcement in the tension side, which is already cracked due to the tensile stresses in the flexural reinforcement that act as a tie.Moving the opening in the beam models with y o /h = 0.33 in the horizontal direction away from the support plate reduced its interruption for the natural load path, and hence, increased the beam strength.The beam models with a y o /h of 0.75 exhibited an opposite trend, where the strength decreased with an increase in the distance measured from the face of the support within the shear span.Moving the opening in the beam models with y o /h = 0.75 in the horizontal direction away from the support plate increased its interruption for the natural load path, and hence, reduced the beam strength.These results verified the detrimental effect of placing a web opening in the compression zone closer to the load plate (x o /X c = 0.75 and y o /h = 0.75).The results also implied that a negligible strength reduction could be obtained when the web opening was placed in the tension side above the flexural reinforcement but away from the natural load path and the support plate (x o /X c = 0.75 and y o /h = 0.33).
reinforcement in the tension side, which is already cracked due to the ten the flexural reinforcement that act as a tie.Moving the opening in the beam yo/h = 0.33 in the horizontal direction away from the support plate reduced i for the natural load path, and hence, increased the beam strength.The beam a yo/h of 0.75 exhibited an opposite trend, where the strength decreased w in the distance measured from the face of the support within the shear spa opening in the beam models with yo/h = 0.75 in the horizontal direction a support plate increased its interruption for the natural load path, and henc beam strength.These results verified the detrimental effect of placing a w the compression zone closer to the load plate (xo/Xc = 0.75 and yo/h = 0.75).T implied that a negligible strength reduction could be obtained when the we placed in the tension side above the flexural reinforcement but away from th path and the support plate (xo/Xc = 0.75 and yo/h = 0.33).

GFRP Stresses
General 3D views of the stresses in the vertical GFRP stirrups predicted numerically for the beam models of this group are shown in Figures 38-40.The maximum stress in the vertical GFRP stirrups did not exceed the tensile strength of the bent portion (459 MPa), except in model DB-X0.50-Y0.50,which exhibited a maximum stress of 477 MPa in a location close to the bent portion.For the beam models with xo/Xc = 0.25 (Figure 38), the

GFRP Stresses
General 3D views of the stresses in the vertical GFRP stirrups predicted numerically for the beam models of this group are shown in Figures 38-40.The maximum stress in the vertical GFRP stirrups did not exceed the tensile strength of the bent portion (459 MPa), except in model DB-X0.50-Y0.50,which exhibited a maximum stress of 477 MPa in a location close to the bent portion.For the beam models with xo/Xc = 0.25 (Figure 38), the maximum stress in the vertical GFRP stirrups ranged from 335 to 431 MPa.The location of the maximum stress in the vertical GFRP stirrups varied based on the value of yo/h.When the opening was close to the tension face (DB-X0.25-Y0.33), the maximum stress was in the vertical GFRP stirrup at the side of the opening.When the opening was pushed upward (DB-X0.25-Y0.50 and DB-X0.25-Y0.75), the maximum stress was in the short vertical GFRP stirrups located in the lower chord below the opening.Figure 39

Deep Beam Models with Different Opening Locations
Kong and Sharp [38] developed Equation ( 1) to determine the ultimate shear strength (Vu) of steel-reinforced deep beams with a web opening, noting that the ultimate load (Pu)

Deep Beam Models with Different Opening Locations
Kong and Sharp [38] developed Equation (1) to determine the ultimate shear strength (V u ) of steel-reinforced deep beams with a web opening, noting that the ultimate load (P u ) of the deep beams of the present study equals 2V u .In these equations, X c = clear shear span, f t = tensile strength of the concrete, b = width of the beam, h = total depth of the beam, A = area of an individual web bar or a main reinforcing bar, y 1 = depth at which a typical bar intersects a potential critical shear crack, α 1 = angle of inclination between a typical bar and the critical shear crack, C 1 = 1.4 for normal weight concrete, and C 2 = 300 N/mm 2 for deformed steel bars.The coefficients a 1 = w o /2X c and a 2 = h o /2h define the opening size.The coefficients K 1 = X o /X c and K 2 = Y o /h define the position of the opening, where X o = horizontal distance between the center of the opening and the inner face of the support plate and Y o = vertical distance between the center of the opening and the bottom face of the beam.The coefficient λ = 1.0 for the main longitudinal bars, whereas for web reinforcing bars, λ = 1.5.Table 7 compares predictions of Kong's Equation (1) [38] for the GFRP-reinforced deep beams with a web opening in the middle of the shear span included in the parametric study to the strengths predicted by the numerical analysis.It is evident that Kong's Equation (1) consistently overestimated the ultimate loads of the deep beam models with a web opening in the middle of the shear span by up to 20%.The unconservative predictions provided by Kong's Equation (1) [38] could be attributed to the reduced dowel action and Young's modulus of the GFRP longitudinal reinforcement.As such, Equation ( 2) is proposed in the present study to represent a refined formula to estimate V u of GFRP-reinforced concrete deep beams with a web opening in the middle of the shear span, where A f = individual area of a main reinforcing bar, A w = individual area of a web reinforcing bar, E f = elastic modulus of the main GFRP reinforcing bars (66.4 GPa), and E s = elastic modulus of steel bars (200 GPa).As shown in Table 7, predictions of the refined analytical formula provided a conservative prediction for the ultimate loads of the deep beam models with an opening in the middle of the shear span.[38] for the GFRP-reinforced deep beams with a web opening at different locations within the shear span.It is evident that Kong's Equation (1) [38] provided inconsistent results.In many cases, the prediction of Equation ( 1) was significantly higher than the strength obtained from the numerical analysis.In some other cases, the prediction of Equation (1) was lower.It is noteworthy that Equation (2) proposed in the current study did not provide satisfactory results when the center of the web opening was shifted from the middle of the shear span.These results indicate that for GFRP-reinforced beams with a web opening shifted from the middle of the shear span, a modified Equation should be adopted.As such, the refined formula given in Equation ( 3) is proposed for the prediction of the shear strength of GFRP-reinforced deep beams with a web opening shifted from the midpoint of the shear span.Two coefficients were introduced in this formula, namely λ 1 and λ 2 , to account for the opening size and location, respectively.These two coefficients were proposed based on careful examination of the ultimate loads of the deep beam models with a web opening shifted from the middle of the shear span.The results of the parametric study showed an increase in the ultimate load with an increase in the distance from the support for the beam models with K 2 values of 0.33 and 0.5.The results of the parametric study also indicated that when the center of the opening was located at a K 2 value of 0.75, a reduction in the ultimate load was recorded with an increase in the distance from the support.The coefficient λ 2 reflects the trend of the results of the parametric study.The predictions of Equation ( 3) are in good agreement with the results of the numerical analysis.It is noteworthy that Equation (3) also provided reasonable predictions (deviation ≤ 12%) for the shear strength of the GFRP-reinforced concrete deep beams with an opening in the middle of the shear span having h o /h ≤ 0.25.

Conclusions
This study provided new knowledge on the structural behavior of GFRP-reinforced concrete deep beams with a web opening through a numerical analysis.The simulation models developed in the current study served as a numerical platform for performance prediction of GFRP-reinforced deep beams with a web opening.The size and location of the web opening within the shear span played a primary role in the behavior of the deep beam models.Placing a web opening in the compression zone close to the support or the load plate was very detrimental to the beam strength.In contrast, a negligible strength reduction was recorded when the web opening was placed in the tension side above the flexural reinforcement and away from the natural load path and the support plate.Specific conclusions of the work are summarized below.

•
For the beam models with y o /h values of 0.33 and 0.50, the strength increased with an increase in the distance measured from the face of the support within the shear span.The strength of the beam models with a y o /h of 0.33 tended to be higher than that of their counterparts with a y o /h of 0.50, and this behavior was more evident with an increase in the distance from the face of the support plate.The beam models with a y o /h of 0.75 exhibited an opposite trend, where the strength decreased with an increase in the distance measured from the face of the support within the shear span because such a movement resulted in an opening closer the load plate.

•
For the beam models having a web opening closer to the support plate (x o /X c = 0.25), strength reductions of 43 and 49% were recorded at y o /h values of 0.33 and 0.50, respectively.When the opening was pushed upward away from the natural load path (x o /X c = 0.25 and y o /h = 0.75), a lower strength reduction of 34% was recorded.

•
For the beam models having a web opening with an x o /X c of 0.50, strength reductions of 24, 42, and 52% were recorded at y o /h values of 0.33, 0.50, and 0.75, respectively.• The beam model with x o /X c = 0.75 and y o /h = 0.33 exhibited a negligible strength reduction of 7% because the web opening was in the tension side and did not interrupt the natural load path.In contrast, the beam model with a web opening closer to the load plate (x o /X c = 0.75 and y o /h = 0.75) exhibited a significant strength reduction of 56% because the web opening was in the compression zone close to the load plate and fully interrupted the natural load path.

•
For the beam models with a web opening in the middle of the shear span, the strength decreased with an increase in either the opening width or height.The rate of the strength reduction caused by increasing the opening height was, however, more significant than that produced by increasing the opening width.At the same w o /a of 0.16, strength reductions of 31-49% were recorded for the deep beam models having an opening in the midpoint of the shear span with h o /h values of 0.17-0.33.More pronounced respective strength reductions of 39-61% were recorded for the deep beam models with the greater w o /a of 0.32.• The existing empirical equation for concrete deep beams reinforced with conventional steel bars with a web opening provided unconservative and/or inconsistent predictions for the ultimate load of the beam models reinforced with GFRP bars.

•
Refined empirical equations were introduced for shear strength prediction of GFRPreinforced concrete deep beams with a web opening of different sizes and locations within the shear span.The refined analytical formulas tended to provide conservative/reasonable predictions for the shear capacity of the GFRP-reinforced concrete deep beams considered in the present study.• The simulation models developed and verified in the present study can be used as a numerical platform in future research to study the effect of using different types of reinforcing bars (e.g., carbon, glass, steel with different yield strengths) on the behavior of concrete deep beams with and without a web opening in the shear span.Future research should further investigate the effect of the anchorage length and bond condition of GFRP reinforcing bars on the response of the GFRP-reinforced deep beams.

Fibers
. The longitudinal reinforcement in the tension side consisted of eight GFRP bars with a diameter of 25 mm (No. 25) placed at a distance d = 1100 mm from the compression face.The compression reinforcement consisted of two GFRP bars with a diameter of 16 mm (No. 16).The internal shear reinforcement consisted of vertical and horizontal GFRP bars with diameters of 12 mm (No. 12) and 16 mm (No. 16), respectively, placed at a spacing of 200 mm in both directions.The deep beam model DB-O2 had 2 extra double-leg GFRP stirrups in the vertical direction (one at each side of the opening), 2 extra horizontal GFRP bars above the opening, and 2 extra horizontal GFRP bars below the opening.The deep beam model DB-O3 had 4 extra double-leg GFRP stirrups in the vertical direction (one at each side of the opening and two crossed by the opening), 4 extra horizontal GFRP bars above the opening, and 4 extra horizontal GFRP bars below the opening.All the extra GFRP reinforcing bars in deep beam models DB-O2 and DB-O3 had a diameter of 12 mm (No. 12).Steel plates (200 mm × 300 mm × 30 mm) were placed at the load and support points.

Fibers 2024 ,Figure 1 .
Figure 1.Details of the deep beam models, noting that bars with red, green, and black colors represent the main flexural reinforcement, the horizontal web reinforcement, and the vertical web reinforcement, respectively (dimensions are in mm): (a) DB-S; (b) DB-O1; (c) DB-O2; (d) DB-O3.

Figure 1 .
Figure 1.Details of the deep beam models, noting that bars with red, green, and black colors represent the main flexural reinforcement, the horizontal web reinforcement, and the vertical web reinforcement, respectively (dimensions are in mm): (a) DB-S; (b) DB-O1; (c) DB-O2; (d) DB-O3.

Figure 4 .
Figure 4. Flexural and web reinforcements: (a) DB-S (Bars with red, green, and brown colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively); (b) DB-O1 (Bars with cyan, red, and green colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively); (c) DB-O2 (Bars with cyan, red, and green colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively); (d) DB-O3 (Bars with cyan, red, and green colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively).

Figure 4 .
Figure 4. Flexural and web reinforcements: (a) DB-S (Bars with red, green, and brown colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively); (b) DB-O1 (Bars with cyan, red, and green colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively); (c) DB-O2 (Bars with cyan, red, and green colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively); (d) DB-O3 (Bars with cyan, red, and green colors represent flexural tensile reinforcement, horizontal web reinforcement, and vertical web reinforcement, respectively).

Figure 8 .
Figure 8. Schematics showing the crack patterns obtained from the experiments based on data from [23]: (a) solid beam; (b) beam with a web opening.

Figures 9 -
show the stresses in the vertical stirrups, horizontal web reinforcement, and flexural reinforcement of the deep beam models at peak load, respectively.Portions of the web GFRP reinforcement crossing the diagonal strut in DB-S exhibited the highest stresses, as shown in Figure9a.The vertical concentrated load caused indirect tensile stresses in the top horizontal portion of the vertical stirrup under the load point in DB-S.The stress in this location was almost equal to that of the tensile strength of the bent portion of GFRP bars, indicating localized rupture of GFRP at this location.These numerical findings are consistent with the published experimental data[23], which indicated that crushing of the diagonal concrete strut in DB-S was accompanied by a localized rupture at the bent portion of the vertical stirrups.For DB-O1, the stress in the bent portions of the vertical GFRP stirrups at the sides of the opening almost reached the tensile strength of the bent portion of GFRP bars (Figure9b).The experimental findings also indicated that rupture of the bent portions of the vertical GFRP stirrups near the load and support plates occurred at the ultimate load along with the crushing of concrete along the diagonal cracks that developed in the upper and lower chords[23].In agreement with the published data[14], the predicted stresses in the vertical stirrups for DB-O2 and DB-O3, having extra GFRP reinforcement around the opening, were below the tensile strength of the bent portion of GFRP bars (Figure9c).These results verify the effectiveness of the extra GFRP bars in relieving the web reinforcement, which allowed the beam models to develop a higher load capacity.The stresses in the horizontal web reinforcement in DB-O2 and DB-O3 at peak load were also lower than those of DB-O1 (Figure10).From Figure11, it can be seen that the maximum stress in the flexural reinforcement at peak load for DB-S was 390 MPa (i.e., 39% of the tensile strength of straight GFRP bars).The models with a web opening failed at a lower load than that of DB-S, and hence, the maximum stress in the flexural reinforcement was lower at an average value of 265 MPa (i.e., 27% of the tensile strength of straight GFRP bars).Published experimental data[23] verified the nonoccurrence of the rupture of the flexural GFRP reinforcement in the tested deep beams.It is noteworthy that the negative stress in the top GFRP longitudinal reinforcement in DB-S at ultimate load was higher than that of their counterparts with a web opening because DB-S failed at a higher load.

4. 1 . 1 .
Figure 15a-c show the load-deflection responses of the deep beam models with a web opening of different sizes placed in the middle of the shear span.The response of the solid deep beam model was included in the figures for the purpose of comparison.The installation of the web opening in the middle of the shear span reduced the shear cracking load and the post-cracking stiffness relative to those of the solid deep beam model.An increase in the opening size further compromised the response of the beam models.The reductions in the cracking load and post-cracking stiffness caused by increasing the opening height were more pronounced for the models with the greater w o /a of 0.32.The beam models with the web opening failed at a deflection significantly lower than that of the solid beam model DB-S.

Figure 19 .
Figure 19.Crack pattern for a typical model failed due to an excessive widening of shear cracks at the opening corners (DB-W0.32-H0.25):(a) at 25% of peak load; (b) at 50% of peak load; (c) at 75% of peak load; (d) at 100% of peak load.

Fibers 2024 ,
12, 66 18 of 35 stirrups at the bent portions in most of the beam models with w o /a values of 0.16 and 0.27.

Figure 21 .FibersFigure 22 .
Figure 21.Three-dimensional view of stresses in GFRP bars for models with w o /a = 0.27 (MPa): (a) DB-W0.27-H0.17;(b) DB-W0.27-H0.25;(c) DB-W0.27-H0.33.Fibers 2024, 12, x FOR PEER REVIEW 19 of 35 e., 21 to 30% of the tensile strength of the GFRP bars) in the horizontal web reinforcement and 200 to 236 MPa (i.e., 20 to 24% of the tensile strength of the GFRP bars) in the tensile flexural reinforcement as shown in Figures 25 and 26, respectively.The beam models with the wo/a value of 0.32 exhibited maximum stresses of 324 to 382 MPa (i.e., 27 to 38% of the tensile strength of the bars) in the horizontal web reinforcement and 169 to 226 MPa (i.e., 17 to 23% of the tensile strength of the GFRP bars) in the tensile GFRP flexural reinforcement as shown in Figures 27 and 28, respectively.These results indicate that for the beam models with the web opening in the

Figure 22 .
Figure 22.Three-dimensional view of stresses in GFRP bars for models with w o /a = 0.32 (MPa): (a) DB-W0.32-H0.17;(b) DB-W0.32-H0.25;(c) DB-W0.32-H0.33.None of the horizontal web reinforcing bars or the flexural reinforcement at the tension side in the deep beam models of this group reached their tensile strength.The stresses in the horizontal web reinforcement and tensile flexural reinforcement predicted numerically for the beam models of this group are shown in Figures23-28.For the beam models with the smaller w o /a value of 0.16, the maximum stresses in the horizontal web reinforcement were on average 258 MPa, which corresponded to 22% of the tensile strength of the GFRP bars (Figure23).The maximum stress in their tensile flexural reinforcement ranged from 190 to 249 MPa, which corresponds to 19 to 25% of the tensile strength of the GFRP bars (Figure24).The beam models with the w o /a value of 0.27 exhibited maximum stresses of 240 to 348 MPa (i.e., 21 to 30% of the tensile strength of the GFRP bars) in the horizontal web reinforcement and 200 to 236 MPa (i.e., 20 to 24% of the tensile strength of the GFRP
. The deep beam models of this group had the same web opening size of w o = 340 mm and h o = 304 mm (i.e., w o /a = 0.27 and h o /h = 0.25).The primary variable of this group was the location of the web opening.The center of the web opening was located at a distance x o from the face of the support plate and y o from the bottom soffit of the beam model.The locations of the center of the opening in the horizontal and vertical directions were normalized to the clear shear span (X c ) and beam depth (h), respectively.Values of the term x o /X c were 0.25, 0.50, and 0.75, whereas those of the term y o /h were 0.33, 0.50, and 0.75.The deep beam models are labeled DB-WM-HN, where M denotes the value of x o /X c and N refers to y o /h.For instance, DB-X0.25-Y0.33 is the designation of the deep beam model having x o /X c of 0.25 and y o /h 0.33.Figures 29-31 show details of the reinforcement for the beam models of this group.

Figure 30 .Figure 30 .
Figure 30.Details of deep beam models with xo/Xc = 0.50, noting that bars with red, green, and black colors represent the main flexural reinforcement, the horizontal web reinforcement, and the vertical web reinforcement, respectively (dimensions are in mm): (a) yo/h = 0.33; (b) yo/h = 0.50; (c) yo/h = 0.75.

Figure 31 .
Figure 31.Details of deep beam models with xo/Xc = 0.75 and different opening vertical locations, noting that bars with red, green, and black colors represent the main flexural reinforcement, the horizontal web reinforcement, and the vertical web reinforcement, respectively (dimensions are in mm): (a) yo/h = 0.33; (b) yo/h = 0.50; (c) yo/h = 0.75.

Figure 31 .
Figure 31.Details of deep beam models with x o /X c = 0.75 and different opening vertical locations, noting that bars with red, green, and black colors represent the main flexural reinforcement, the horizontal web reinforcement, and the vertical web reinforcement, respectively (dimensions are in mm): (a) y o /h = 0.33; (b) y o /h = 0.50; (c) y o /h = 0.75.
75 was seriously affected by the location of the center of the opening in the vertical direction (y o ).The behavior of the beam model DB-X0.75-Y0.33 was not seriously affected by the opening since it was provided at the tension side (y o /h = 0.33) away from the natural load path.In contrast, the behavior of the beam model DB-X0.75-Y0.75 with y o /h = 0.75 was seriously compromised since the opening was very close to the load plate and fully interrupted the natural load path.Initially, the behavior of the beam model DB-X0.75-Y0.50 with y o /h = 0.50 almost coincided with that of DB-X0.75-Y0.33 with y o /h = 0.33 until the load reached a value of approximately 1250 kN.Next, a degradation in the behavior of the beam model DB-X0.75-Y0.50occurred until it failed at load and deflection capacities lower than those of DB-X0.75-Y0.33.It is noteworthy that a change in the slope of the load-deflection response was observed prior to reaching the ultimate load, probably because of the development of new major cracks.

FibersFigure 32 .
Figure 32.Effect of opening location on the load-deflection response: (a) xo/Xc = 0.25 (b) xo/Xc = 0.50; (c) xo/Xc = 0.75.Table 6 presents the ultimate load and the deflection capacity for the beam models of this group.For the beam models having a web opening with an xo/Xc of 0.25, strength reductions of 43 and 49% were recorded for DB-X0.25-Y0.33 and DB-X0.25-Y0.50having yo/h values of 0.33 and 0.50, respectively.It is noteworthy that the web opening in DB-X0.25-Y0.33 and DB-X0.25-Y0.50fully interrupted the natural load path.When the opening was pushed upward away from the natural load path in DB-X0.25-Y0.75, a lower strength reduction of 34% was recorded.The deflection capacity of the beam models

Figure 32 .
Figure 32.Effect of opening location on the load-deflection response: (a) x o /X c = 0.25; (b) x o /X c = 0.50; (c) x o /X c = 0.75.

Figure 33 .
Figure 33.Effect of the opening location on the strength.

Figure 33 .
Figure 33.Effect of the opening location on the strength.4.2.2.Crack Pattern and Failure MechanismFigures[34][35][36][37] show the crack development and propagation at different stages of loading along with the final pattern for sample deep beam models of this group.The beam models of this group, except those with the web opening close to the compression face of the beam (y o /h = 0.75), exhibited first cracking at the opposite corner of the opening closer to the support and load points at a load value ≤ 25% of their strength.The beam models with y o /h = 0.75 exhibited very few flexural cracks at 25% of the strength, whereas shear cracks were initiated at the opening corners at a load value higher than 25% of the strength.The shear cracks initiated earlier at the opening corners propagated diagonally toward the support and load plates as the load progressed.Flexural cracks were also initiated at the midspan, within the shear span, and in the bottom chord below the opening.Additional shear cracks were developed with an increase in the applied load creating a band of shear cracks.The final crack pattern and failure mode was dependent on the location of the opening.The beam models with the web opening away from the natural load path (DB-X0.25-Y0.75 and DB-X0.75-Y0.33)exhibited failure of the diagonal strut connecting the load and support plates.The beam model DB-X0.75-Y0.50failed due to extensive diagonal shear cracking crossing the opening center.Other beam models failed due to the development of extensive diagonal cracks that caused the failure of the concrete along the upper and lower load paths or failure of the concrete along one of them combined with extensive diagonal cracks along the other load path as well as excessive widening of shear cracks at the opening corners.

Figure 34 .
Figure 34.Crack pattern for a typical model that failed along the strut connecting the load and support plates (DB-X0.25-Y0.75):(a) at 25% of peak load; (b) at 50% of peak load; (c) at 75% of peak load; (d) at 100% of peak load.

Figure 35 .
Figure 35.Crack pattern for a typical model that failed due to extensive diagonal shear cracking at the opening corners (DB-X0.75-Y0.50):(a) at 25% of peak load; (b) at 50% of peak load; (c) at 75% of peak load; (d) at 100% of peak load.

Figure 37 .
Figure 37. Crack pattern for a typical model that failed due to failure of the concrete along the lower load path and excessive widening of shear cracks at the opening corner (DB-X0.25-Y0.50):(a) at 25% of peak load; (b) at 50% of peak load; (c) at 75% of peak load; (d) at 100% of peak load.

4. 2 Figure 37 .
Figure 37. Crack pattern for a typical model that failed due to failure of the concrete along the lower load path and excessive widening of shear cracks at the opening corner (DB-X0.25-Y0.50):(a) at 25% of peak load; (b) at 50% of peak load; (c) at 75% of peak load; (d) at 100% of peak load.

Figure 40 .
Figure 40.Three-dimensional view of stresses in GFRP bars for models with x o /X c = 0.75 (MPa): (a) DB-X0.75-Y0.33;(b) DB-X0.75-Y0.50;(c) DB-X0.75-Y0.75.Figures 41-46 show stresses in the horizontal web reinforcement and tensile flexural reinforcement predicted numerically for the beam models of this group.None of the horizontal web reinforcing bars or the flexural reinforcement at the tension side in the deep beam models of this group reached their tensile strength.Generally, the maximum stress in the horizontal GFRP bars was on average 275 MPa (23% of the tensile strength of straight GFRP bars), with a minimum of 188 MPa (16% of the tensile strength of straight GFRP bars) and a maximum of 325 MPa (27% of the tensile strength of straight GFRP bars).The stresses in the flexural tensile GFRP reinforcement within the shear span were almost constant, except at the region near the support which exhibited reduced values of GFRP stresses.The near-uniform stress distribution of the GFRP stress in the flexural GFRP reinforcing bars within the shear span indicated that they acted as a tie, which is in alignment with the behavior of D-regions.The maximum stress in the flexural GFRP bars was on average 240 MPa (24% of the tensile strength of straight GFRP bars), with a minimum of 176 MPa (18% of the tensile strength of straight GFRP bars) and a maximum of 386 MPa (39% of the tensile strength of straight GFRP bars).The value of the maximum GFRP stress was dependent on the value of the ultimate load.The beam models with the higher ultimate load typically exhibited higher GFRP stresses in the flexural reinforcement at peak load.For instance, the smallest longitudinal GFRP stresses of 176, 184, and 190 MPa were recorded for the beam models DB-X0.25-Y0.50,DB-X0.50-Y0.75, and DB-X0.75-Y0.75having the lowest ultimate loads.The greatest longitudinal GFRP stress of 386 MPa was recorded for the beam model DB-0.75-0.33 having the greatest ultimate load.

Table 1 .
Key input parameters of the concrete.
* Values between parenthesis are for DB-O3.

Table 1 .
Key input parameters of the concrete.
* Values between parenthesis are for DB-O3.

Models Ultimate Load (kN) Deflection Capacity (mm) Experimental 1 Numerical Error (%) 2 Experimental 1 Numerical Error (%) 2
-H0.25 is the designation of the deep beam model having a w o /a of 0.16 and h o /h of 0.25.Figures12-14show details of the reinforcement of the numerical models with different opening sizes in the middle of the shear span.wo values of 200, 340, and 400 mm correspond to wo/a values of 0.16, 0.27, and 0.32, respectively.** ho values of 200, 304, and 400 mm correspond to ho/h values of 0.17, 0.25, and 0.33, respectively. *

Table 3 .
Parameters of deep beam models with different opening sizes.w o values of 200, 340, and 400 mm correspond to w o /a values of 0.16, 0.27, and 0.32, respectively.** h o values of 200, 304, and 400 mm correspond to h o /h values of 0.17, 0.25, and 0.33, respectively. *

Table 4 .
Numerical results of the deep beam models with different opening sizes.

Table 4 .
Numerical results of the deep beam models with different opening sizes.

Table 5 .
Parameters of deep beam models with different opening locations.

Table 5 .
Parameters of deep beam models with different opening locations.

Table 6 .
Numerical results of the deep beam models with different opening locations.x o values of 262.5, 525, and 787.5 mm correspond to x o /X c values of 0.25, 0.50, and 0.75, respectively.** y o values of 400, 590, and 900 mm correspond to y o /h values of 0.33, 0.50, and 0.75, respectively. *

Table 7 .
(1)parison between predictions of the analytical formulas and numerical results of deep beams with an opening in the middle of the shear span.Deep Beams with a Web Opening Shifted from the Midpoint of the Shear Span Table8compares predictions of Kong's Equation(1) 1Error (%) = Formula−Numerical Numerical ×100.

Table 8 .
Comparing predictions of the analytical formulas and numerical results of deep beams with an opening shifted from the midpoint of the shear span.