Next Article in Journal
Characteristics of Hydroxyapatite-Modified Coatings Based on TiO2 Obtained by Plasma Electrolytic Oxidation and Electrophoretic Deposition
Next Article in Special Issue
Cumulative Fatigue Damage of Composite Laminates: Engineering Rule and Life Prediction Aspect
Previous Article in Journal
Effects of Difluoro(oxalato)borate-Based Ionic Liquid as Electrolyte Additive for Li-Ion Batteries
Previous Article in Special Issue
Material and Process Tests of Heterogeneous Membranes Containing ZIF-8, SiO2 and POSS-Ph
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Haunch Geometry and Additional Steel Reinforcement on Load Carrying Capacity of Reinforced Concrete Box Culvert

1
Department of Civil Engineering, Ghulam Ishaq Khan Institute of Engineering Sciences and Technology, Topi, Swabi 23640, Pakistan
2
Faculty of Mechanical Engineering and Marine Technology, University of Rostock, Justus-von-Liebig-Weg 6, 18059 Rostock, Germany
3
Department of Life, Light and Matter, University of Rostock, Albert Einstein-Str. 25, 18059 Rostock, Germany
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(4), 1409; https://doi.org/10.3390/ma16041409
Submission received: 22 November 2022 / Revised: 26 January 2023 / Accepted: 6 February 2023 / Published: 8 February 2023
(This article belongs to the Special Issue Analysis and Design of Structures and Materials)

Abstract

:
The culverts are used to safely convey water under railways, highways, and overpasses. They are utilized in drainage areas or water channels and in areas where the bearing capacity of soil is low. The design and construction of this crucial infrastructure need to be improved to meet contemporary demands of reliability and affordability. Precast reinforced box culverts are popular alternatives as they ensure strength, durability, rigidity, and economy. This research seeks to develop an effective and affordable design improvement procedure for a precast box culvert using modern numerical tools. The Finite Element Method (FEM) based approach is used in studying the effects of haunch geometry and additional steel reinforcement on the load-bearing capacity of box culverts. A conventional box culvert is analyzed to create the numerical models in the Abaqus FEM code and to investigate the load-bearing capacity of culverts with an expanded span. The outcomes of the study reveal the critical places for stress concentration as well as the location of maximum damage. It is found that haunch geometry and additional reinforcement at these critical places significantly affect the load-carrying capacity of a culvert. From the comparison of capacity curves of models with and without haunches and diagonal reinforcement, it is found that a 25% increase in load-carrying capacity is achievable with the recommended changes. The proposed design improvement technique can be employed for the cost-effective and safe design of a concrete box culvert with larger span lengths and high water-flowing capacities. The findings of this study are expected to assist practitioners in strength enhancement tasks of box culverts for increased structural stability and drainage efficiency.

1. Introduction

A culvert is an essential part of a drainage system that helps to manage water flows and provides a dependable course for water to pass through embankments of highways and railways and helps prevent flooding [1]. There are two main types of culverts: rigid and flexible. Rigid culverts, made of concrete, are designed to withstand bending loads with minimal deformation. Flexible culverts, made of steel, are designed to interact with the soil structure and transfer loads differently [2]. The suitability of a culvert type is determined by its ease of construction and required drainage capacity. Box culvert structures are economical since they have immensely high rigidity due to their monolithic action, and these do not require isolated foundations. For small discharges, single-cell or box culverts are used; for large discharges, multi-cell or multi-bin culverts are normally utilized [3]. The top slab of box culverts is designed to support the live loads of a moving vehicle. Similarly, the walls and base slab resist the earth pressure and hydrostatic pressure from outside and from inside, respectively [4].
Climate change and meteorological variations were particularly considered during the original design of old hydraulic structures [5], and whenever there is heavy precipitation, well-planned culverts can drain the water both quickly and effectively. However, outdated or poorly designed culverts cause upstream floods and lead to critical losses and damage [6]. Therefore, the capacity of culverts should be enhanced to ensure adequate water discharging and load-carrying capabilities.
Efforts have been made in the past to develop cost-effective ways of enhancing the drainage and load-carrying capacity of culverts. One approach is the hydraulic modification of existing culverts to increase their capacity rather than completely rebuilding the entire structure. Contemporary design standards assess the performance of a culvert by computing the inflow and outflow of water based on established hydraulic theories while the least efficient flow control is chosen based on anticipated headwater levels. However, by implementing the better design improvement method, the flow and overall performance of the culvert could be improved [7].
In recent years, significant research has been conducted on experimental and numerical analysis and assessment of a box culvert. Crisfeld [8] analyzed the effects of cracking and the failure stress response of a box culvert. Marzouk et al. [9] devised a technique for the determination of the punching and deformation demand of a reinforced concrete (RC) box culvert. Polak [10] employed three-dimensional (3D) shell elements for shear and flexural strength assessment of RC box culverts.
Some researchers experimentally examined the effect of earth pressure on box culverts and compared the results with existing theoretical models to offer design improvements [11,12]. Several academics have studied box culverts to understand the impact of soil-structure interaction on the performance of box culverts [13,14]. Gong et al. carried out experimental and numerical evaluations to study the failure mechanism of reinforced concrete box culverts [15]. Anil and Ali [16] utilized FEM to validate the experimental findings of the box culvert. Moradi et al. [17] validated the experimental works of Maximos et al. [18] and Garg [19] through Finite Element Methods (FEM) and found that FEM tools can simulate the structural behavior of culverts with reasonable accuracy and lesser physical effort. Sinha and Sharma [20] analyzed the design of these small bridges with an emphasis on the specifications for the Indian road congress. Anderson et al. [21] proposed seismic design methods for the culvert. Esra [22] studied the response of long-span corrugated steel box culverts under ultimate loads and proposed a method to adjust the design force calculations. Kim et al. [23] investigated the failure mechanism of box culverts subjected to corrosion by employing FEM. Zenagebriel et al. [24] performed an experimental and numerical investigation on precast RC U-shaped box culverts to examine the effect of various loading conditions both on the slab and sides of the culvert.
Recently, a smart precast box culvert equipped with fiber optic sensors was designed. These sensors were able to detect and measure strain and stress on the culvert in order to monitor structural health via the stress–strain measurement [25]. Additionally, recent work has been focused on reducing the effects of load on box culverts under different loads of construction such as culverts near tunnels or next to train tracks [26,27,28].
Despite these efforts, the studies related to design improvements and load-carrying capacity enhancement of culverts with various strengthening methods have not been explored so far. Against this background, this paper uniquely presents a numerical assessment of a non-conventional box culvert with a larger span length through FEM. The purpose is to study the effects of various strengthening methods and haunch geometry with additional reinforcement on the strength and drainage capacity enhancement of the culvert structure. Furthermore, the FEM is used to simulate the structural behaviors of unstrengthened and strengthened culverts in terms of their deformation response and stress concentrations at critical places of the box culvert using Abaqus software [29].
FEM is an essential technique for predicting events that take place during the fabrication and operation of objects and structures. It has accelerated the advancements in numerous fields of research, such as fluid mechanics, mechanics of structures, fracture mechanics, heat conduction, and biomechanics. It has been proven that this approach, which is labor-intensive to build, calculate, and analyze models, is successful in resolving problems of design improvements [30,31,32]. The Abaqus FEM code has been effectively used in much culvert-related research for design and analysis [16,33,34,35,36,37].
The inner span of a box culvert is normally about 3.5 m [4]. For larger spans and discharges, the conventional size becomes ineffective and calls for the design of multiple sections thus leading to increased costs and construction efforts. A larger-span culvert is therefore an economical alternative. However, its structural integrity and efficiency need further exploration. This research provides design recommendations for larger-span culverts using the FEM approach. The load-bearing capacity of conventional culverts with an increased span length is evaluated, and the deformed shape, along with the stress distribution of the reference culvert, is utilized to indicate the crucial zones of failure. At these critical failure areas, additional reinforcement and haunch geometry are employed as capacity improvement measures. The efficiency of various strengthening techniques is evaluated on a larger span culvert to obtain the increased design capacity of a culvert for supporting heavier loads and conveying large quantities of water. The adopted analysis mechanism effectively utilizes the FEM techniques to analyze multiple scenarios and conveniently identify an efficient design improvement method for a culvert structure. The presented research work is expected to help practitioners in updating the existing design and analysis techniques and provide an effective methodology for the design enhancement of box culverts to increase both structural stability and drainage efficiency. The complete scheme of this study is depicted in Figure 1.

2. Materials and Methods

2.1. Model Description

The numerical model of the culvert was developed using the geometry and reinforcement details in the Abaqus FEM code for capacity enhancement. Moreover, the Concrete Damage Plasticity (CDP) model was used to simulate the constitutive behavior of concrete [38] as shown in Figure 2.
To simulate the inelastic behavior of concrete, the damaged plasticity model for concrete combines the ideas of isotropic damaged elasticity with isotropic tensile and compressive plasticity. The concrete damage model is a continuous, plasticity-based model. The model is intended for applications in which low confining pressures and monotonic, cyclic, and dynamic stresses are applied to the concrete. This model decomposes the total strain ε into elastic εel and plastic strains εpl as follows.
ε =   ε el + ε pl
The following equations give the stress–strain relations under uniaxial compression and tension:
σ t = 1 d t E o   ε t ε t pl
σ c = 1 d c E o ε c ε c pl
where Eo is the initial (undamaged) elastic stiffness of the material, ε t is total strain, ε c pl and ε t pl are compression and tensile equivalent plastic strains, and dt and dc are tension and compression damage parameters, respectively.
The two primary failure mechanisms inferred by the model are the compressive crushing and tensile cracking of the concrete material. In the plastic-damage model, a plastic-damage variable called “d” was defined. If plastic deformation occurs, the value of this variable increases. The non-dimensional variable “d” has a maximum value of one. When the plastic-damage variable attains its maximum value at a particular location in the solid, the total damage takes place in the form of a macroscopic crack. The two damage parameters for compressive and tensile behavior were used as functions of equivalent plastic strains. The relationship between the tension damage variable and cracking strain is shown in Figure 3. These parameters control the evolution of the yield surface and are related to failure mechanisms under tension and compression loading. The damage parameters can have values between zero and one, where one corresponds to a complete loss of strength and zero to the undamaged material.
The plastic flow of the material was simulated by flow parameters in the concrete damaged plasticity model including dilation angle (ψ), flow potential eccentricity (e), viscosity parameter (μ), and coefficient determining the shape of the deviatoric cross-section (K). The damaged plasticity model assumes non-associated potential flow. The flow potential G follows the Drucker-Prager hyperbolic function:
G = e σ t 0 tan φ 2 + q ¯ 2 p ¯ tan φ
where φ is the dilation angle measured in the p−q plane at high confining pressure, σ t 0   is the uniaxial tensile stress at failure, and e is referred to as the eccentricity, the rate at which the function approaches the asymptote. These continuous and smooth flow potential parameters ensure that the flow direction is defined uniquely. The function asymptotically approaches the linear Drucker-Prager flow potential at high confining pressure stress and intersects the hydrostatic pressure axis at 90.
The biaxial behavior of concrete material was considered with an initial equi-biaxial compressive yield stress (σbo) and an initial uniaxial compressive yield stress (σco) ratio (σboco). In the case of uniaxial loading, cracks move in a direction opposite to the stress direction and the available load-carrying area decreases as a result of the growth of cracks and rise in the effective stress. The effective load-carrying area is greatly diminished once a significant quantity of material fails. The effective uniaxial cohesion stress at any instance determines the size of the yield surface. The effective cohesion stress in compression and tension is provided below:
σ t ¯ = σ t 1 d t = E o ε t ε t pl
σ c ¯ = σ c 1 d c = E o ε c ε c pl
The traffic and embankment loads were modeled as equally distributed loads on the surface, but the soil and water loads were applied as depth-varying triangular loads. To derive load defection curves and stress distributions, the analytical equations were solved using the Abaqus solver of the Dassault Systèmes Simulia Corporation. The load–deflection curve was used to compute the load-carrying capacity, and the stress distribution of the structure was used to pinpoint the critical failure points. The geometry of the haunches and additional steel reinforcement in the optimized culvert structure was utilized to increase its capacity.

2.2. Geometric Details

A typical RC box culvert with an extended span of more than 3.5 m was selected as a reference model for this study. The box culvert’s structure is made up of a steel reinforcement cage and a concrete box, as seen in Figure 4. The main and distribution horizontal and vertical bars of the reinforcing cage had rebar sizes of 16 mm. The inner span of the culvert is 4.0 m, and the wall and slab thickness are both preserved at 0.4 m. The normal reinforcing spacing is 200 mm, and the concrete cover thickness is 40 mm.

2.3. Materials and Elements

The material properties of concrete for the CDP model utilized are based on past research [39]. These properties are listed in Table 1 for materials such as steel and concrete. Table 2 details the plastic flow parameters of the CDP model.
The elements of concrete and steel are discretized for FEM using 3D wire elements for reinforcement and three-dimensional solid elements for concrete components [24]. The components for reinforcement were 3D, 2 node (T3D2) elements. T3D2 elements are used to simulate slender, line-like structures that only enable axial loading along the element and do not support moments or forces that are perpendicular to the centerline. The concrete elements were meshed with continuum, 3D, 8 node elements with reduced integration (C3D8R). C3D8R elements are linear brick elements with fewer integration points and are standard elements for stress analysis. A courser mesh would yield faster solution convergence while a finer mesh would provide improved results with additional computational cost. The mesh convergence was undertaken by varying the mesh sizes of 200, 100, 50, and 20 mm for concrete, as well as reinforcement, and the results are provided in the section on modal validation.

2.4. Boundary Conditions, Interactions, and Loading

The bond between concrete and steel was taken into account using the embedded region method of the Abaqus FEM numerical package. By evaluating the stiffness of the reinforcement elements independently of the concrete elements, this method of modelling the steel reinforcement resolves the mesh restriction issue that arises in the discrete and smeared modelling of reinforcement. With this technique, the host element (concrete) and the slave element are perfectly joined (steel rebar). Additionally, this approach allows for the displacement of the surrounding concrete components to be compatible with the displacement of the steel bars. When applied to complex models, the embedded method is extremely helpful. However, this model expands the number of nodes and degrees of freedom; consequently, it needs more run time and costs more to compute.
The culvert’s bottom end was modelled as a fixed foundation using fixed boundary conditions. To simulate field conditions, the remaining surfaces received the necessary weight of soil and water to properly constrain the model with the required boundary conditions. With the appropriate loadings applied to the locations for lateral earth pressure and hydrostatic pressure, the top surface was loaded with a 75 mm layer of asphalt and a 0.3 m layer of soil [40]. The loading conditions and model boundary conditions are shown in Figure 5.
For the performance evaluation of various strengthening methods, the displacement-controlled loading of 120 mm was applied at the reference point to obtain the capacity curve. The motion of one group of nodes on the top slab was restricted by kinematic coupling constraints to rigid body motion defined by a reference node. These restrictions were helpful in controlling the deflection of one node, which was linked to a large number of nodes (the “coupling” nodes). The specification of a reference node, coupling nodes, and the constrained degrees of freedom at these nodes was provided for a kinematic coupling constraint. A kinematic link allowed the reference point (RP-1) and top surface to interact with one another. The load vs. displacement curve was then produced at the reference position by this method.
The critical stress locations in the reference model were sought from the analysis results. These crucial locations were subsequently retrofitted with various strengthening configurations to identify the efficient approach to strength enhancement.

3. Results and Discussion

3.1. Modal Validation

A numerical model of the culvert was prepared to indicate the effectiveness of the FEM technique and modal validation. The model was constructed based on the experimental work of Maximos et al. [18]. The experimental study used a full-scale precast box culvert for testing as per the American Society of Testing and Materials (ASTM) standards [41,42,43]. The specimen comprised dimensions of 2120 × 1220 × 203 mm where the first two dimensions correspond to inner short- and long-span lengths of culverts, and the third dimension represents the thickness of walls and slabs of a culvert. The tested segment was 1220 mm long, and the reinforcement was provided as per ASTM C1577. The specimen was made with a concrete compressive strength of 52.40 MPa and deformed bars of 470 MPa. The box culvert was subjected to a wheel load at the top of the slab using a footprint plate, and the structure was loaded until failure. The failure occurred when cracks developed in the negative moment sections of the side walls at a load of 244.65 kN (24.5 Tonnes).
For the numerical modeling of the test specimen, the procedures mentioned in Section 2.1, Section 2.3, and Section 2.4 were utilized to simulate the constitutive behavior of the materials and apply the required modeling parameters. The numerical model of the tested specimen for modal validation is shown in Figure 6. For modal calibration, the load vs. deflection response of the structure with different mesh sizes of steel and concrete was obtained at the point of application of loading. The results of this analysis are provided in Figure 7a, and Table 3 compares the ultimate load and deflection response of structures with different mesh sizes. An analysis of Figure 7a and Table 3 revealed that the error in the results was minimized when the mesh size was set at 50 mm. Moreover, as the mesh size was decreased below 50 mm, the results remained consistent and insignificant variations were observed. The results were not found to be mesh-sensitive due to the small size of the elements; however, the size of the element chosen was slightly larger to reduce the computational time. The corresponding element size of 100 mm was assigned to reinforcing steel and concrete as an optimum arrangement between computational time and modal performance [44]. The mesh size of 100 mm for concrete and steel components provided a comparable ultimate capacity with respect to experimental output. From the experiment, the test specimen exhibited an ultimate capacity of 24.5 tonnes while the FEM reported an ultimate capacity of 27 tonnes at failure. Furthermore, the crack patterns observed in the work of Maximos et al. and the FE model were also found to be similar as presented in Figure 7b. Therefore, the resulting numerical modeling parameters showed good agreement with the experimental results where an experiment/FEM ratio of 0.9 would be obtained and an error value of about 10% could be expected.

3.2. Parametric Analysis

The parametric analysis was conducted to evaluate the influence of various parameters on the load-carrying capacity of an RC box culvert. The reference model detailed in Section 2.2 was considered for the parametric study and a total of four parameters were considered in the analysis including the strength of concrete, strength of steel reinforcement, amount of steel reinforcement, and haunch geometry. The culvert size and shape, load applications, and boundary conditions remained unchanged for this parametric analysis. The details of the variables selected for parametric analysis are provided in Table 4.
The concrete strength was varied by using concrete mixes of different compressive strengths. Three different concrete strengths of 26.8, 40, and 52.4 MPa were considered. Three different steel strengths of 400, 470, and 550 MPa were included. Similarly, the influence of three different amounts of reinforcement was investigated. Similarly, the amount of reinforcement in the reference model was increased by 150% to 200%. Additionally, the haunch geometry was varied by changing the shape of the haunch. Three different haunch configurations of 300 mm size were studied i.e., a square haunch, a triangular haunch, and a semicircular haunch. The results of the parametric analysis are provided in Table 5 and Figure 8 to identify the influence of each parameter on the ultimate load-carrying capacity of a culvert structure.
A total of 12 analyses were undertaken for four parameters to evaluate their influence on the load-carrying capacity of the RC box culvert. The concrete strength showed a direct impact on the load-carrying capacity of the culvert. The high-strength concrete of double strength could improve the load-carrying capacity by about 25%. Similarly, the steel reinforcement strength and amount also showed a significant impact on the load-carrying capacity of the culvert where the maximum reinforcement amount provided the highest load-carrying capacity, followed by the moderate reinforcement, and the low reinforcement. The 200% increase in the reinforcement amount enhanced the load-carrying capacity by about 35%. Additionally, the haunch geometry showed a significant impact on the load-carrying capacity of the culvert. The provision of haunch geometry could increase the load-carrying capacity by about 20–25%. Overall, the amount of steel reinforcement and the haunch geometry showed a significant impact on the load-carrying capacity of the culvert, followed by the concrete and steel strength.
Since high-strength concrete and steel are more expensive than normal-strength reinforced concrete materials, the combined use of haunch geometry and reinforced concrete material should be encouraged for cost-reduction and effective utilization of resources. This can be realized by identifying the critically distressed areas in the structure where the material performance should be improved.
To identify the efficient approach to strength enhancement, the identified crucial locations were retrofitted with steel reinforcement configurations listed below.
  • Case 1: Near Surface Mounted (NSM) L-shaped reinforcement bars;
  • Case 2: NSM diagonal brackets;
  • Case 3: Concrete haunches and diagonal reinforcement.
The rebar size of 16 mm was used for these additional reinforcements. The NSM bars of cases 1 and 2 were embedded in concrete 20 mm from the surface. The spacing of all additional rebars was maintained at 200 mm, similar to the main reinforcement. For case 3, concrete haunches at a size of 0.3 m × 0.3 m were provided at the stress concentration locations. Figure 9a,b show the arrangements of bars for cases 1 and 2, respectively. The graphical representation of concrete haunches and diagonal bars of case 3 are presented in Figure 10. The load–displacement curves of these reinforcement configurations were compared, and their effectiveness was reported to suggest strength-enhancement strategies.

3.3. Stress Distribution and Ultimate Load-Carrying Capacity

The stress distribution in the reference model under the induced loading is shown in Figure 11. The corner intersection of the culvert’s top slab and walls revealed significant stress concentrations in addition to other typical places. Severe deformation at mid-span and in the joints of the walls and in the top slab was also observed. Excessive mid-span deflection is observed once the load-carrying capacity of the slab is exhausted. The deformation at joint locations further increased the mid-span deflection and reduced the load-carrying capacity. The ultimate load-carrying capacity of the reference culvert was calculated using the load–displacement curve shown in Figure 12. This figure depicts that the culvert structure could bear the maximum load of 210 tonnes with a small deformation of 22 mm. Once the load-resisting capacity of the structure is exhausted due to material damage, the structure experiences extremely large deflections with reduced load-bearing capabilities. The capacity of the reference structure was decreased to 100 tonnes until a 60 mm deformation. The increase in capacity again appeared up to 160 tonnes with a deformation of 80 mm because of the material hardening of steel rebars. No further increase in load-carrying capacity was observed subsequently, and the structure showed continuous large deformations. From Figure 11, it can be observed that the key stressed points are at the intersection of the top slab and walls. Moreover, from the stress distribution and deformed shape, it can be assumed that the strengthening of a culvert structure by means of additional reinforcement and a haunch geometry at the critical position would significantly improve the load-carrying capacity of a longer-span culvert. The strengthened culvert with an extended span would provide better water-flowing capability and would be able to withstand more loads.

3.4. Influence of Additional Steel Reinforcement and Concrete Haunches

In order to carefully follow the advancements in structural behavior and overall capacity, the influence of additional steel reinforcement and concrete haunches was investigated by comparing the load-displacement curves of the reference model and strengthened structures. The procedure described in the preceding section was utilized to determine the capacity of the retrofitted structure. The improved capacity curves of different cases are shown in Figure 12, next to the reference capacity curve. The results indicate that the installation of L-shaped NSM reinforcement would not be beneficial to improve the load-carrying capacity when compared to a diagonal bracket arrangement, which exhibited slight improvement in the capacity of the reference structure. The concrete haunch structure and additional diagonal reinforcements were found to be the most effective, providing a capacity of 270 tonnes. The culvert structure, strengthened with concrete haunches and diagonal bars, has the ability to increase the box culvert’s load-bearing capacity by around 25%. The difference between the peak load-carrying capacity of reference (208.3 tonnes) and the strengthened culvert (270 tonnes) was used to obtain this amount.
To summarize, this study evaluated the impact of different reinforcing methods on the load-carrying capacity of a culvert with an increased span length, which has not been undertaken to date. The NSM strengthening method has been adopted in past studies for various structural components such as beams and has shown some beneficial effects in enhancing structural performance [44]. However, in the case of culverts, this is a unique study and the first of its kind. The performance of NSM reinforcement in the capacity enhancement of a culvert was found to be insignificant compared to concrete haunch geometry and diagonal reinforcement. The current study uniquely identifies the performance of concrete haunch geometry and the diagonal reinforcement-based strengthening method and highlights the significance of the proposed design improvement method.
The haunch geometry and additional diagonal reinforcement played a crucial role in enhancing the structural resistance at the location of stress concentration and preventing the box culvert from excessive deformations. Figure 13 provides the comparison of induced deflections of 200 mm in the slabs of reference and haunched culverts. At the mid-span location, both reference and stiffened structures were forced with the same amount of deflection, and hence the deflections were almost similar; however, the deflection is reduced towards the wall and slab joint due to the provision of the haunch and diagonal reinforcement. The haunch and diagonal reinforcement at the joints provided additional stiffness at the joint location where the deflection in the stiffened structure is almost zero. The concrete haunch and diagonal bars also played a vital role in reducing the stress concentration at the critical joint location as shown in Figure 14. The quantity of elements experiencing a high number of stresses in reference culverts is reduced by adopting retrofitting measures. The quantity of elements with a smaller number of stresses is also increased in a retrofitted culvert which indicates that the stresses are well distributed in different elements of joints with lower absolute maximum values. The proposed methodological framework of the strength enhancement of culverts systematically identified the weaker zones in the structure and recommended the most effective strength enhancement method for culverts. The proposed framework presented in Figure 1 could be repeated to achieve the desired load-carrying capacity of a culvert structure with extended span lengths.

4. Conclusions

This study evaluated the impact of different strengthening methods on improving the load-carrying capacity of a culvert with increased span length and larger water-flowing capacity. The presented work demonstrated the effectiveness of the FEM for incorporating the required design improvements based on the output of the analysis. A typical culvert structure with an enlarged span was numerically modeled to act as a reference structure. NSM steel rebars and concrete haunches with diagonal strengthening methods were employed at critical failure locations of the structure for design improvement and strength enhancement of the reference model. The selected strength enhancement approaches had been commonly used in previous studies for various structural components, and the present work uniquely explored the application of these methods on culvert structures. The models with NSM reinforcement were found to be insignificant in the capacity enhancement of culverts compared to models where concrete haunch geometry and diagonal reinforcement were considered. The load–deflection curves of the reference model and the model strengthened with different strengthening methods were compared to obtain the ultimate load-carrying capacities of different models. The deformed shapes of the structure and stress distribution of a reference culvert and the culvert strengthened with the most effective design improvement methods were also considered to highlight the performance of a recommended approach of capacity enhancement. The following conclusions were drawn from the current study:
  • The link between the slab and the wall was found to have the highest likelihood of stress concentrations;
  • The NSM rebars of an L-shape and diagonal bracket slightly influence the capacity and could be useful for small repairs and restoration works;
  • The precast reinforced box culvert’s FE modelling indicates that adding diagonal reinforcement and a haunch construction could appreciably enhance its capacity. A 25% increase in the ultimate load-carrying capacity was found when the ultimate capacity of models with and without haunches was compared.
To fully understand the relevance of the suggested design improvement technique for cost-effective and safe design, the performance of box culverts with wider spans and under seismic loadings must also be examined. Furthermore, the influence of other retrofitting configurations, such as the Fiber Reinforced Polymer (FRP), to investigate strength enhancements should also be explored.

Author Contributions

Supervision, H.A.W., H.S., M.W. and A.R.; Conceptualization, methodology, visualization, H.A.W., M.I., M.N. and M.W.; Data curation, Software, H.A.W. and M.N.; writing—original draft, H.A.W., M.W., M.N. and M.I.; writing—review and editing, A.R. and H.S.; Validation, H.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. APC was funded by Deutsche Forschungsgemeinschaft (DFG) and the Open Access Publication Fund of the University of Rostock.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smart, P.; Graham Herbertson, J. (Eds.) Drainage Design; Blackie: Glasgow, UK, 1992. [Google Scholar]
  2. McGrath, T.J.; Moore, I.D.; Selig, E.T.; Webb, M.C.; Taleb, B. NCHRP Report 473: Recommended Specifications for Large-Span Culverts; Transportation Research Board: Washington, DC, USA, 2002. [Google Scholar]
  3. Mackie, K.; Stojadinovic, B. Seismic Demands for Performance-Based Design of Bridges, PEER Report 2003–16; Pacific Earthquake Engineering Research Center, University of California: Berkeley, CA, USA, 2003. [Google Scholar]
  4. AASHTO. LRFD Bridge Design Specifications, 9th ed.; American Association of State Highway and Transportation Officials: Washington, DC, USA, 2020. [Google Scholar]
  5. Stagrum, A.E.; Andenæs, E.; Kvande, T.; Lohne, J. Climate change adaptation measures for buildings—A scoping review. Sustainability 2020, 12, 1721. [Google Scholar] [CrossRef]
  6. Douglas, I.; Alam, K.; Maghenda, M.; Mcdonnell, Y.; McLean, L.; Campbell, J. Unjust waters: Climate change, flooding and the urban poor in Africa. Environ. Urban. 2008, 20, 187–205. [Google Scholar] [CrossRef]
  7. Jaeger, R.; Tondera, K.; Pather, S.; Porter, M.; Jacobs, C.; Tindale, N. Flow control in culverts: A performance comparison between inlet and outlet control. Water 2019, 11, 1408. [Google Scholar] [CrossRef]
  8. Crisfeld, M.A. Snap-through and snap- back response in concrete structures and the dangers of under-integration. Int. J. Numer. Meth. Eng. 1986, 22, 751–767. [Google Scholar] [CrossRef]
  9. Marzouk, H.; Osman, M.; Hussein, A. Punching shear of slabs: Crack size and size effects. Mag. Concr. Res. 2002, 54, 13–21. [Google Scholar] [CrossRef]
  10. Polak, M.A. Shell fnite element analysis of RC plates supported on columns for punching shear and fexure. Int. J. Comput.-Aided Eng. Softw. 2005, 22, 409–428. [Google Scholar] [CrossRef]
  11. Inyeop, C.; Woo, S.-K.; Woo, S.I.; Kim, J.; Lee, K. Analysis of Vertical Earth Pressure Acting on Box Culverts through Centrifuge Model Test. Appl. Sci. 2021, 12, 81. [Google Scholar]
  12. Lee, K.; Kim, J.; Woo, S.I. Analysis of Horizontal Earth Pressure Acting on Box Culverts through Centrifuge Model Test. Appl. Sci. 2022, 12, 1993. [Google Scholar] [CrossRef]
  13. Al-Zaidee, S.R.; Alsalmani, A.H.; Salam, A. Effects of Soil-structure Interaction and Different Simulation on the Response of Concrete Box Culverts. IOP Conf. Ser. Mater. Sci. Eng. 2020, 901, 012009. [Google Scholar] [CrossRef]
  14. Ozturk, K.F.; Cakir, T.; Araz, O. A comparative study to determine seismic response of the box culvert wing wall under influence of soil-structure interaction considering different ground motions. Soil Dyn. Earthq. Eng. 2022, 162, 107452. [Google Scholar] [CrossRef]
  15. Gong, Y.; Ma, Y.; Tan, G.; Bi, H.; Pang, Y.; Ma, C. Experimental Study and Numerical Simulation on Failure Process of Reinforced Concrete Box Culvert. Adv. Civ. Eng. 2020, 2020, 5423706. [Google Scholar] [CrossRef]
  16. Anil, G.; Ali, A. Finite-element modeling and analysis of reinforced concrete box culverts. J. Transp. Eng. 2009, 135, 121–128. [Google Scholar]
  17. Moradi, M.; Valipour, H.; Foster, S. Reserve of Strength in Inverted U-Shaped RC Culverts: Effect of Backfill on Ultimate Load Capacity and Fatigue Life. J. Bridge 2016, 21, 04015051. [Google Scholar] [CrossRef]
  18. Maximos, H.; Erdogmus, E.; Tadros, M.K. Fatigue evaluation for reinforced concrete box culverts. ACI Struct. J. 2010, 107, 13–20. [Google Scholar]
  19. Garg, A.K.; Abolmaali, A.; Fernandez, R. Experimental investigation of shear capacity of precast reinforced concrete box culverts. J. Bridge Eng. 2007, 12, 511–517. [Google Scholar] [CrossRef]
  20. Sinha, B.; Sharma, R. RCC box culvert—Methodology and design including computer method. J. Indian Roads Congr. 2009, 555, 198–219. [Google Scholar]
  21. Anderson, D.; Wang, J.; Martin, G.; Lam, I. Proposed seismic design methods for transportation culverts and drainage structures. Lifeline Earthq. Eng. Multihazard Environ. 2009, 15, 1187–1198. [Google Scholar]
  22. Esra, B.F. Response of long-span box type soil-steel composite structures during ultimate loading tests. J. Bridge Eng. 2009, 14, 496–506. [Google Scholar]
  23. Kim, D.; Bhowmik, S.; Willmer, J.; Leo, E. A case history and finite element modeling of a culvert failure. Geo-Frontiers 2011, 2011, 1943–5606. [Google Scholar]
  24. Zenagebriel, G.; Li, J.; Qiao, G. Experimental and Finite Element Analysis of Precast Reinforced Concrete U-Shaped Box Culvert. Iran. J. Sci. Technol. Trans. Civ. Eng. 2021, 45, 2133–2143. [Google Scholar] [CrossRef]
  25. Gebremedhn, Z.; Qiao, G.; Sun, L.; Li, H.; Bai, S.; Jin, H. Strain/Stress Self-Sensing Precast-Box Culvert: Smart Element Design, and Experimental Verification. Int. J. Civ. Eng. 2021, 19, 851–859. [Google Scholar] [CrossRef]
  26. Hussien, O.S. Reduction of box culvert stresses. Environ. Eng. Manag. J. 2020, 19, 1969–1974. [Google Scholar] [CrossRef]
  27. Lu, X.; Liu, Y.; Zhang, L. Research on Deformation Characteristics of Tunnel Box Culvert Structure under Different Construction Loads. IOP Conf. Ser. Earth Environ. Sci. 2021, 634, 012136. [Google Scholar] [CrossRef]
  28. Tarigan, J.; Nursyamsi, N.; Farhan, M. Precast box culvert analysis with direct load. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1122, 012020. [Google Scholar] [CrossRef]
  29. Smith, M. ABAQUS/Standard User’s Manual, Version 6.9; Dassault Systèmes Simulia Corp: Providence, RI, USA, 2009. [Google Scholar]
  30. Umravia, N.B.; Solanki, C.H. Numerical Analysis to Study Lateral Behavior of Cement Fly Ash Gravel Piles under the Soft Soil. Int. J. Eng. 2022, 35, 2111–2119. [Google Scholar] [CrossRef]
  31. Farahani, S.; Mohebkhah, A. Seismic Behavior of Direct Displacement-based Designed Eccentrically Braced Frames. Int. J. Eng. 2016, 29, 752–761. [Google Scholar]
  32. Mohsenzadeh, S.; Maleki, A.; Yaghin, M.A. Experimental and Numerical Study of Energy Absorption Capacity of Glass Reinforced SCC Beams. Int. J. Eng. 2019, 32, 1733–1744. [Google Scholar] [CrossRef]
  33. Wang Shuhong, A.J.; Jie, R.; Wang, P.; Liu, W. Simulation of force performance of precast rectangular box culvert and its potential failure mode. J. Northeast. Univ. 2018, 2, 260–265. [Google Scholar]
  34. Awwad, E.; Mabsout, M.; Sadek, S.; Tarhini, K. Finite element analysis of concrete box culverts. Comput. Civ. Build. Eng. 2000, 2000, 1051–1053. [Google Scholar]
  35. Gebremedhn, Z.; Qiao, G.; Li, J. Finite Element Modeling and Analysis of Precast Reinforced Concrete U-Shaped Box Culvert Using ABAQUS. Am. J. Civ. Eng. 2018, 6, 162. [Google Scholar] [CrossRef] [Green Version]
  36. Uddin, Z.; Waqas, H.A.; Husnain, M.; Ibrahim, M.; Zaib, M.N.H. Design Optimization of Precast Reinforced Box Culvert with the Application of Finite Element Modeling and Analysis. In Proceedings of the 1st International Conference on Advances in Civil & Environmental Engineering, University of Engineering & Technology Taxila, Rawalpindi, Pakistan, 22–23 February 2022. [Google Scholar]
  37. Ozturk, B. Seismic behavior of two monumental buildings in historical Cappadocia region of Turkey. Bull. Earthq. Eng. 2017, 15, 3103–3123. [Google Scholar] [CrossRef]
  38. Hibbitt, D.; Karlsson, B.; Sorensen, P. ABAQUS Analysis User’s Manual; Abaqus: Providence, RI, USA, 2022. [Google Scholar]
  39. Xiao, Y.; Chen, Z.; Zhou, J.; Leng, Y.; Xia, R. Concrete plastic-damage factor for finite element analysis: Concept, simulation, and experiment. Adv. Mech. Eng. 2017, 9, 1687814017719642. [Google Scholar] [CrossRef]
  40. Ali, A.M. Manual RC Box Culvert Analysis and Designing. Int. J. Innov. Sci. Res. Technol. 2020, 5, 91–115. [Google Scholar] [CrossRef]
  41. ASTM C1577; Standard Specification for Precast Reinforced Concrete Box Sections for Culverts, Storm Drains, and Sewers Designed According to AASHTO LRFD. ASTM International: West Conshohocken, PA, USA, 2005; p. 17.
  42. ASTM C1433; Standard Specifications for Precast Reinforced Concrete Box Sections for Culverts, Storm Drains, and Sewers. ASTM International: West Conshohocken, PA, USA, 2004; p. 25.
  43. Karalar, M.; Yeşil, M. Investigation on Seismic Behavior of Historical Tokatli Bridge under Near-Fault Earthquakes. Adv. Civ. Eng. 2021, 2021, 5596760. [Google Scholar] [CrossRef]
  44. Al-Ameedee, H.S.; Al-Khafaji, H.M. Improving the flexural behavior of RC beams strengthening by near-surface mounting. J. Mech. Behav. Mater. 2022, 31, 701–709. [Google Scholar] [CrossRef]
Figure 1. The methodological framework adopted in the current research study.
Figure 1. The methodological framework adopted in the current research study.
Materials 16 01409 g001
Figure 2. Schematic diagram of the CDP model (image used with the permission of publisher) [39].
Figure 2. Schematic diagram of the CDP model (image used with the permission of publisher) [39].
Materials 16 01409 g002
Figure 3. Concrete tension damage vs. cracking strain.
Figure 3. Concrete tension damage vs. cracking strain.
Materials 16 01409 g003
Figure 4. FEM of reinforcement cage and concrete box culvert.
Figure 4. FEM of reinforcement cage and concrete box culvert.
Materials 16 01409 g004
Figure 5. Details of model constrains and loading mechanism; (a) model constrains; (b) kinematic coupling for displacement-controlled loading.
Figure 5. Details of model constrains and loading mechanism; (a) model constrains; (b) kinematic coupling for displacement-controlled loading.
Materials 16 01409 g005
Figure 6. Numerical model based on the experimental work of Maximos et al. [18].
Figure 6. Numerical model based on the experimental work of Maximos et al. [18].
Materials 16 01409 g006
Figure 7. (a) Performance of numerical model with different mesh sizes and (b) crack patterns observed in the experiment of Maximos et al. [18] (image used with the permission of publisher).
Figure 7. (a) Performance of numerical model with different mesh sizes and (b) crack patterns observed in the experiment of Maximos et al. [18] (image used with the permission of publisher).
Materials 16 01409 g007
Figure 8. Summary of results of the parametric analysis.
Figure 8. Summary of results of the parametric analysis.
Materials 16 01409 g008
Figure 9. Arrangements of NSM reinforcement.
Figure 9. Arrangements of NSM reinforcement.
Materials 16 01409 g009
Figure 10. Introduction of concrete haunches and diagonal reinforcement for capacity enhancement.
Figure 10. Introduction of concrete haunches and diagonal reinforcement for capacity enhancement.
Materials 16 01409 g010
Figure 11. Stress distribution and location of stress concentration.
Figure 11. Stress distribution and location of stress concentration.
Materials 16 01409 g011
Figure 12. Load vs. displacement curves.
Figure 12. Load vs. displacement curves.
Materials 16 01409 g012
Figure 13. Comparison of deflections in the slabs of reference and haunched culverts.
Figure 13. Comparison of deflections in the slabs of reference and haunched culverts.
Materials 16 01409 g013
Figure 14. Stress distribution in the wall and slab joints of reference and haunched culverts.
Figure 14. Stress distribution in the wall and slab joints of reference and haunched culverts.
Materials 16 01409 g014
Table 1. Material information of the study culvert [39].
Table 1. Material information of the study culvert [39].
DetailsConcreteSteel
Density (Kg/m3)24007850
Compressive strength (MPa)26.8450
Tensile strength (MPa)2.39450
Poisson’s Ratio0.20.3
Modulus of elasticity (MPa)32,500200,000
Post yielding modulus of elasticity (MPa)-20,000
Table 2. Plasticity flow parameters for CDP model [39].
Table 2. Plasticity flow parameters for CDP model [39].
ParametersValues
Dilation angle (ψ) 40°
Eccentricity (e)0.1
Ratio of biaxial to uniaxial compressive yield stress (σboco)1.16
Coefficient determining the shape of the deviatoric cross-section (K) 0.66
Viscosity parameter (μ)0
Table 3. Results of model calibration.
Table 3. Results of model calibration.
Mesh SizeDisplacement Force
ConcreteSteelmmTonnes
200200535
1001003.4527
50504.524.0
20204.523.9
Selected mesh size100 mm for concrete and steel% Error in results of numerical analysis10
Table 4. Variables for parametric analysis.
Table 4. Variables for parametric analysis.
ParameterVariablesUnits
Strength of concrete26.8, 40, 52.4MPa
Strength of steel280, 450, 550MPa
Amount of steel100, 150, 200Percent
Haunch geometrysquare, triangular, and semicircular-
Table 5. Results of parametric analysis.
Table 5. Results of parametric analysis.
ParameterVariableUnits
Strength of concrete26.84052.4MPa
Ultimate load208.3233.3262.1Tonnes
Increase in capacity01225.8Percent
Strength of steel280450550MPa
Ultimate load192.8208.3222Tonnes
Increase in capacity−7.406.6Percent
Amount of steel100 150200Percent
Ultimate load208.3242.8280.1Tonnes
Increase in capacity016.634.5Percent
Haunch geometrysquare triangular semicircular-
Ultimate load261.5253.2258.1Tonnes
Increase in capacity25.621.623.9Percent
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Waqas, H.A.; Waseem, M.; Riaz, A.; Ilyas, M.; Naveed, M.; Seitz, H. Influence of Haunch Geometry and Additional Steel Reinforcement on Load Carrying Capacity of Reinforced Concrete Box Culvert. Materials 2023, 16, 1409. https://doi.org/10.3390/ma16041409

AMA Style

Waqas HA, Waseem M, Riaz A, Ilyas M, Naveed M, Seitz H. Influence of Haunch Geometry and Additional Steel Reinforcement on Load Carrying Capacity of Reinforced Concrete Box Culvert. Materials. 2023; 16(4):1409. https://doi.org/10.3390/ma16041409

Chicago/Turabian Style

Waqas, Hafiz Ahmed, Muhammad Waseem, Abdullah Riaz, Muhammad Ilyas, Muhammad Naveed, and Hermann Seitz. 2023. "Influence of Haunch Geometry and Additional Steel Reinforcement on Load Carrying Capacity of Reinforced Concrete Box Culvert" Materials 16, no. 4: 1409. https://doi.org/10.3390/ma16041409

APA Style

Waqas, H. A., Waseem, M., Riaz, A., Ilyas, M., Naveed, M., & Seitz, H. (2023). Influence of Haunch Geometry and Additional Steel Reinforcement on Load Carrying Capacity of Reinforced Concrete Box Culvert. Materials, 16(4), 1409. https://doi.org/10.3390/ma16041409

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop