1 Introduction

Many hollow slab bridges are designed to be skewed bridges to satisfy the traffic requirements. The live-load response of this kind of bridge could be affected by span length, bridge deck thickness and skewed angle, etc. To describe and simplify the effect of these parameters, the distribution factor (DF) is proposed by assigning fractions of transmitted load amongst the primary structural members (Harris 2010; Song et al. 2003).

Most hollow slab bridges in the design plan were designed as simply supported. However, it was reported that even slight changes in boundary conditions have a considerable effect on the results (Schulz et al. 1995). Bridge test in Ontario had shown that slab-on-girder bridges are usually significantly stiffer in flexure because of horizontal restraint provided by the girder bearings (Bakht and Jaeger 1988). Emon and Nowak (2001) proposed that the hinge-roller supports can be partially fixed (frozen) due to many factors and they found that the code-specified DF values without considering the effect of possible support fixity can be too conservative. The actual boundary condition for the bridge was quite important for analyzing the bridge structural behavior. Many field tests were conducted on the existing bridges and the field test responses were suggested to calibrate the FEM for further research (Seo et al. 2017; Hodson et al. 2011, Bechtel et al. 2010). Bridge response under loading test was the interaction of various parameters. However, field test data on the existing bridge couldn’t calibrate all of these parameters. As a result, parameters had few effect on the bridge behavior were considered to be the design value in the FEM. For hollow slab bridges, adjacent girders were connected transversely by the hinge joint to attract load together. For accurate analysis, transverse connection between adjacent girders and support condition are two important parameters that need to be calibrated in the FEM for this type of bridge. Multi-stage loading tests on the constructing bridge were an effective way to calibrate these two parameters, respectively.

Many studies investigated the effect of secondary elements such as parapets, sidewalk and diaphragm on the DFs and concluded that AASHTO LRFD specification was conservative without considering these secondary elements (Mabsout et al. 1997; Conner and Huo 2006; Namy et al. 2015). Geometric dimensions of the bridge remained unchanged in these studies. However, geometric dimensions of the bridge were the variables of the DF calculation formula in AASHTO LRFD specification. For a certain kind of bridge, geometric dimensions could be varied in an appropriate range. As a result, the designed DFs for different geometric dimensions were also changed in a certain range. So the effect of geometric dimensions, which were varied in a rational range, on the DFs should be figured out first before discussing the effect of secondary elements.

The aim of this study is to identify the supporting condition and transverse connection of the hollow slab bridge in the field tests using multi-stage construction of full-scale test in situ. First, a concentrated loading test was conducted on the prefabricated girder that settled on the supports before the girders were connected transversely. The response of the prefabricated girder could be used to calibrate the support condition of the hollow slab in the FEM. Second, a vehicle loading test was performed on the bridge after the prefabricated girders were transversely connected. Transverse connection between adjacent girders could be calibrated with the bridge response under the vehicle loading test. In addition, a parametric study was conducted to investigate the influence of various parameters on the DFs of this hollow slab bridge including span length, skew and bridge deck thickness. Finally, the effect of these parameters on the DFs was quantified compared with the AASHTO LRFD specification.

2 Bridge Description

A constructing hollow slab bridge, which is located on WeiXu expressway in Henan province, China, was selected for this study. The superstructure was simply supported prefabricated concrete hollow slab girders, which was shown in Fig. 1. The span of bridge was 13 m in length with a 12.412 m clear span. The width of the bridge was 13.5 m and skewed at \(25^\circ\). The hollow slab girders had a height of 0.7 m and a width of 0.99 m, and covered 0.1 m thick layer of cast-in-place reinforced deck. The bridge girders were prefabricated and the deck was cast. The connection of structures in the bridge superstructure is shown in Fig. 2. The concrete of hinge joints and bridge deck was cast at the same time. The concrete used in the slab girders, hinge joints and bridge deck had a designed compressive strength of 50 MPa and modulus of elasticity of \(3.45 \times 10^{4} {\text{MPa}}\). The superstructure was composed of 13 girders, in which the leftmost girder was denominated G1 while the rightmost girder denominated G13.

Fig. 1
figure 1

Field tested bridge (units: mm).

Fig. 2
figure 2

Connection configuration of the bridge superstructure.

2.1 Concentrated Force Loading Test

2.1.1 General

This loading phase was performed after the prefabricated girders were installed and the adjacent girders had not been connected transversely. Before this loading test, a preloading test was carried out on the tested girders to ensure the accuracy of the instruments. G1 (exterior girder), G4 (interior girder) and G7 (medium girder) girders were selected to be the test girders. The loading girder, which weighed 160.44 kN, was also a prefabricated hollow slab girder. One end of the loading girder was simply supported on the bearings in the mid-span of the two tested girders and the other end of the loading girder was simply supported outside the abutment (Fig. 3). As a result, the analysis of the tested girder could be simplified to be a prefabricated girder with a concentrated force on the mid-span. So the concentrated force applied on the mid-span of each tested girder was 80.22 kN. The concentrated force loading tests were performed on the G1 and G4, G4 and G7, respectively. The dimension of the girder and material properties were measured in this loading phrase.

Fig. 3
figure 3

Concentrated force loading test.

2.1.2 Instrument Arrangement

During the loading test, four strain gauges were placed on the middle section of each tested girders to measure strains, while one dial indicator was placed to measure the mid-span deflection. Where, two dial indicators were installed on the tested girders near the supports to measure the support settlement. The deflection caused by the concentrated force on the mid-span of the tested girder could be eventually determined. Figures 4 and 5 show the gauge configurations. Totally, for the three slabs, the gauges for the strain (D1–D12) were located on the top and bottom flanges of the girder, while the dial gauges for the mid-span deflection (DN1–DN3) were located on the bottom flange of the girder. The strain gauges were sealed by insulating tape as moisture-proof measurement. The lead wires from these gauges went through the hinge joint from the top to the bottom, so that the gauges could be used in the next loading phase after the bridge deck was cast. A strain gauge was placed on a stress free concrete as the common temperature measurement. The test data acquired in this loading phase was expected to study the actual support conditions of the tested bridge.

Fig. 4
figure 4

Strain and deflection gauges.

Fig. 5
figure 5

Measurement instrument position (units: mm). Note: Di: strain gauges; DNi: deflection gauges.

3 Vehicle Loading Test

3.1 General

This loading phase was performed after the slabs were connected transversely, the bridge deck was cast and the concrete strength met the design requirements. In order to excite the normal response of the bridge, Heavy truck loaded 360 kN was used as the loading vehicle. The truck was weighed and measured before its arrival on site. Figure 6 shows the parameters of the loading vehicle.

Fig. 6
figure 6

Loading vehicle measurement (units: mm).

During the loading test, the loading truck was positioned at each of the five transverse loading cases (LC) to produce different response (Fig. 7). When the loading truck was driven to the predetermined position and the engine was stalled, there was an interval of 5 min to collect the strain and deflection data until the effect of the moving truck on the bridge disappeared. Once the collection of the data was complete in one location, the loading truck was driven off the bridge and waited for 5 min so that the deformation of the tested girder could be restored, then, the vehicle loading test could be performed again on the next transverse location. Loading and unloading tests were performed repeatedly to acquire test data.

Fig. 7
figure 7

Vehicle loading test configurations (units: mm).

3.2 Instrument Arrangement

Strain gauges and dial indicators were all installed on the bottom of the mid-span section of the girder (Fig. 8), In which, six gauges that were mounted in the top flange of the three girders in the former loading test were also used in this loading phase to measure the strain of the top flange. Again, two indicators were placed at the supports and one indicator was placed in the mid-section for each girder. Consequently, there were 19 strain gauges and 39 dial indicators in this loading phase.

Fig. 8
figure 8

Instrument arrangement in the vehicle loading test. Note: Pi: strain point, PNi: deflection point.

4 Finite-Element Model

In order to assure the correction of the finite element analysis (FEA), two different FEA programs—MIDAS CIVIL and ANSYS were conducted. In MIDAS CIVIL, the main girder and virtual transverse beam were modeled by the beam element while the bridge deck was modeled by the plate elements. In ANSYS program, beam188 elements were used to model the main girder and virtual transverse beam, shell63 element were used to model bridge deck. For the beam element in MIDAS CIVIL and ANSYS, each node for these elements had three translational degrees of freedom and three rotational degrees of freedom. As for the plate element and shell63 element, there were six degrees of freedom at each node were used to model bridge deck. Material and structural properties in the FEM were based on the bridge design plan as well as the collected and measured information about the bridge supplemented with engineering judgment(Seo et al. 2015). The concrete grade using for the girder, hinge joint and bridge deck is C50, The modulus of elasticity of the concrete is \(3.45 \times 10^{4} {\text{MPa}}\). As specified in the AASHTO LRFD specification, live load deflection is a service issue, not a strength issue. So in the calculation of the FEM, the nonlinearity of the concrete materials was not considered. Figure 9 shows the FE models of the bridge.

Fig. 9
figure 9

Finite element models.

4.1 Boundary Conditions

Elastic spring elements were used to simulate actual behavior of supports as shown in Fig. 10 (Eom and Nowak 2001). To model this kind of connection, general connection options in elastic connection were selected in MIDAS CIVIL, while the combin14 elements were modeled in ANSYS model. In this calibrating process, main girders were not connected transversely. Elastic spring elements were added to restrict the horizontal movement of the bridge on the base of the simply support condition according to the design plan. The stiffness of the spring elements was represented by K values. The suitable K values were found by comparing the deflections at the mid-span from FEM analysis with those from the concentrated force loading test by trial and error analysis. The K values used in the analysis are shown in Table 1, in which, girder 1 and 13 had different K values with the other girders because of their different section from the others. The K values that acquired from different FE models were the same.

Fig. 10
figure 10

Boundary condition used in finite elements.

Table 1 Spring coefficients used in FEM.

4.2 Transverse Connection

After the boundary conditions of the FEM were determined, transverse connections between adjacent girders were calibrated with the field test data acquired in vehicle loading tests. 100 mm thick plate elements and shell63 elements were included in this modeling phrase to simulate the effect of bridge deck. Beam elements and beam188 elements used for transverse connection between adjacent girders were defined as virtual transverse beam elements, which had no quality but only transverse stiffness between adjacent girders. The cross sections of virtual transverse beams were rectangles with a width of 1000 mm. The height of the virtual transverse beam was iterated to fit the test results. The truck loads were applied as concentrated loads and were positioned on the FEM as presented in the vehicle loading test. By trial and error, the suitable virtual transverse beam height was determined by comparing correlation and error between FEM and field test responses. Percent error,\(\delta_{\text{p}}\), and correlation coefficient, \(\rho\) (Seo et al. 2017) were used in the comparison:

$$\delta_{p} = 100 \cdot \left[ {\frac{{\sum {\left( {\omega_{f} - \omega_{a} } \right)^{2} } }}{{\sum {\left( {\omega_{f} } \right)^{2} } }}} \right]$$
(1)
$$\rho = \left[ {\frac{{\sum {\left( {\omega_{f} - \overline{{\omega_{f} }} } \right)\left( {\omega_{a} - \overline{{\omega_{a} }} } \right)} }}{{\sum {\sqrt {\left( {\omega_{f} - \overline{{\omega_{f} }} } \right)^{2} \left( {\omega_{a} - \overline{{\omega_{a} }} } \right)^{2} } } }}} \right]$$
(2)

where \(\omega_{f}\), \(\omega_{a}\), \(\overline{{\omega_{f} }}\) and \(\overline{{\omega_{a} }}\) = field deflection response, FEM deflection responses, sample mean values of \(\omega_{f}\), and sample mean values of \(\omega_{a}\), respectively.

Virtual transverse beam height with highest correlation coefficient and the lowest percent error was determined to be the ultimate height. Table 2 presents a summary of the percent errors and correlation coefficients among different virtual transverse beam heights. The correlation coefficients of different transverse beam heights are 1, which means using virtual transverse beam method could simulate transverse connection of adjacent girders pretty well. Virtual transverse beam with 300 mm height has the lowest percent error in both FEM analyses. So the virtual transverse beam with a width of 1000 mm and a height of 300 mm was found to forecast most closely actual transverse connections between adjacent girders.

Table 2 Model accuracy for different virtual transverse beam heights.

5 Live Load Distribution Factors

The actual DF of each girder can be calculated using the experimental deflection data. The effective deflection under vehicle loading was calculated by the following Eq. 3 as shown in Fig. 11.

Fig. 11
figure 11

The deflection of the test girder.

$$\omega_{e} { = }\omega_{t} - \omega_{r} - \omega_{c}$$
(3)

where \(\omega_{e}\) is the effective measured deflection of girders under loads, \(\omega_{t}\) is the total measured deflection of girders, \(\omega_{r}\) is the measured residual deflection of girders, \(\omega_{c}\) is the measured deflection at both ends of girders.

The formula used here to calculate DF for the ith girder is Eq. 4 (Harris et al. 2008), which has been widely used in many live-load tests (Civitillo et al. 2014; Harris et al. 2016; Waldron et al. 2005). These DFs are representative of a single-truck loading resisted by each girder. DFs in five different loading conditions were shown in Fig. 12.

Fig. 12
figure 12

DFs in different loading conditions.

$$DF_{i} = \left( {{\raise0.7ex\hbox{${\omega_{\text{ i}} }$} \!\mathord{\left/ {\vphantom {{\omega_{\text{ i}} } {\sum\limits_{{{\text{i}} = 1}}^{13} {\omega_{\text{ i}} } }}}\right.\kern-0pt} \!\lower0.7ex\hbox{${\sum\limits_{{{\text{i}} = 1}}^{13} {\omega_{\text{ i}} } }$}}} \right) \times N_{vehicles}$$
(4)

where, \(DF_{i}\) is load distribution factor of ith girder; \(\omega_{i}\) is maximum deflection of ith girder; \(N_{vehicles}\) is the number of vehicles on the bridge.

It can be concluded that the girder which suffered the traffic load directly experiences the greatest deflection, whereas girders far away from the loading position experience less. DFs for exterior girder were 52% greater than the DFs for the most heavily loaded interior girder, the reason is that in this test there are no railings which usually offered a significant stiffening contribution to the exterior girders and the exterior girder is unable to distribute load to an interior member.

Figure 13 shows the deflections at mid-span on each girder under different loading conditions. It can be seen from Fig. 13 that the deflection curves were nearly overlapped for the five loading conditions. It can be concluded that the FEM is reliable to conduct further researches.

Fig. 13
figure 13

Deflection distribution in different loading conditions.

Distribution factors for single and multilane loaded are different in AASHTO specification. Only one lane loaded was conducted in the field test. The number of lane was changed in the calibrated FEM to observe the deflection of girders. The loading vehicle used in the FEM was the same as it used in field test and the vehicle wheel loads were applied as concentrated loads in the FEM. The bridge was divided into three 3.75 m-wide-lanes. FEM with concentrated load was shown in Fig. 14. Loading vehicles were positioned on the middle of each lane laterally. The longitudinal position of trucks was calculated as the position producing the maximum bending moment at midspan. The number and position of the load vehicles were systematically varied to consider all possible load combinations. In all cases the vehicles were placed following the \(25^\circ\) skew.

Fig. 14
figure 14

FEM with three loaded lanes.

Interior and exterior distribution factors are calculated separately in the AASHTO specification. As such, the maximum distribution factor for the one, two and three loaded lanes for both the interior and exterior girders were calculated using the maximum deflection of different girders. Calculated moment DFs for the interior and exterior girders for each load cases are listed in Table 3. DFs calculated by the AASHTO LRFD were also listed in Table 3.

Table 3 Distribution factors for different number of lanes.

Table 3 shows that using AASHTO LRFD specifications resulted in a conservative DFs estimate. The AASHTO DFs were much higher than those calculated using the calibrated FEM for exterior girders. For one lane loaded, the AASHTO LRFD DFs were 1.5 times larger than that of the FEM DFs for exterior girders and 1.2 times larger for interior girders. For multilane loaded, the AASHTO LRFD DFs were 1.4 times larger than that of the FEM DFs for exterior girders and were very close to the FEM DFs for interior girders. DFs for exterior girder were controlling DFs in AASHTO LRFD specification, while the controlling DFs were interior girder DFs in the FEM. This is probably caused by the lack of parapets in the FEM, which could provide stiffness for the exterior girders. DFs for interior girder in AASHTO LRFD specification was 1% to 18% higher than the DFs in FEM. For exterior girder, DFs of AASHTO LRFD specification was 37% to 48% higher than the DFs in FEM.

As the loading lane increased from one to two in the FEM, DFs increased 17% for exterior girders and 34% for interior girder. While as the loading lane increased from two to three, the DFs increased by 2% and 9% for exterior and interior girders, respectively. DFs could be different for one lane loaded and two lanes loaded, but less difference for two lanes loaded and three lanes loaded. In multilane loaded cases, DFs of three loaded lanes were more closed to the AASHTO LRFD specification compared with the DFs of two loaded lanes. When the bridge was filled with vehicles on the lanes that it can be divided into, the DFs were more closed to the AASHTO LRFD specification for this multilane loaded case.

6 Parametric Study

A parametric study was performed for hollow slab bridge to investigate the parameters affecting the DFs of vehicle loads. The effect of parameters on the DFs were quantified and compared with AASHTO LRFD specifications. Span length, bridge skew and bridge deck thickness were three major parameters in this parametric study. A baseline FEM was determined using the dimensions and supporting conditions of the field tested bridge. The baseline model was adjusted to study the effect of one parameter while the other parameters were kept in the baseline value. While changes in certain parameters (i.e., span length) do affect other parameters (i.e., girder stiffness), the effect was minimal and they could be evaluated independently (Zokaie 2000). Therefore, span length was varied from 8 to 20 m, which was the specified length of the hollow slab bridge, in the FEM. The effect of skew was investigated by adjusting the skew angle from 0 to \(60^\circ\). The bridge deck thickness was varied from 0.1 m to 0.5 m by adjusting the thickness of plate element. DFs calculated from the FEM were compared with those calculated ones according to the AASHTO LRFD specifications, this ratio of which could reflect the difference between FEM and AASHTO LRFD specifications.

DFs for one loaded lane and three loaded lanes were calculated in the Table 4 to investigate the effect of span length. It can be seen from Table 4 that DFs for exterior girder increased 10% for one loaded lane and 42% for three loaded lanes, while for interior girder, DFs decreased 30% for one loaded lane and 9% for three loaded lanes. As the DFs calculated by MIDAS CIVIL and ANSYS were really close, so the DFs in Fig. 15 were the average values of the DFs in two different FEM analyses. In Fig. 15, DF ratios increased and were closed to 1 as the span length increased for both one loaded lane and three loaded lanes. This trend showed that AASHTO LRFD specification became less conservative as the span length increased. The ratio of the FEM to AASHTO distribution factor for exterior girder ranged from 0.45 to 0.81 for one loaded lane and from 0.45 to 0.76 for three loaded lanes. While for interior girder, the ratio of the FEM to AASHTO LRFD specification distribution factor was 0.76 to 0.91 for one lane and 0.84 to 0.95 for three lanes. The lack of parapets in the FEM, which could provide stiffness for the exterior girder, made the DF of exterior girder lower than that of interior girder.

Table 4 DFs for span parameter in the FEM.
Fig. 15
figure 15

FEM/AASHTO ratio for length parameter. Note: REXT-i: DF ratio of exterior girder in i loaded lane; RINT-i: DF ratio of interior girder in i loaded lane.

DFs in different bridge skew were shown in Table 5. Compared with the DFs of the straight bridge (of which the skew is \(0^\circ\)), the DFs decreased 5% for exterior girder and 4% for interior girder when the bridge skew was \(5^\circ\). For exterior girder, DFs were not changed for the skew of \(10^\circ\), \(15^\circ\) and \(25^\circ\), while the DFs decreased 3% to 11% as the skew increased from 35 to \(60^\circ\). DFs of interior girder were not changed as the skew increased from \(10\) to \(45^\circ\) and the DFs increased 3% when the skew was \(60^\circ\). DFs in Fig. 16 were also the average DF values of two different FEM analyses. It can be seen in the Fig. 16 that AASHTO became less conservative as the skew increased from \(10^\circ\). The ratio of FEM and AASHTO was greater than 1 for interior girder when the skew was \(60^\circ\), while the DFs of exterior girder were less than design value for each of skew angle.

Table 5 DFs for skew parameter in the FEM.
Fig. 16
figure 16

FEM/AASHTO ratio for skew parameter.

The calculation results from MIDA CIVIL and ANSYS were pretty close. DFs in Fig. 17 were the average values of the two FEM analyses. DFs for exterior girder increased 4% and 5% for one and three loaded lanes, respectively, while for interior girder, DFs decreased 3% and 2% for one and three loaded lanes, respectively. Neither the DF of exterior nor the interior girder was influenced by the deck thickness significantly in Fig. 17.

Fig. 17
figure 17

DFs for bridge deck thickness parameter.

7 Conclusions

A hollow slab bridge was subjected to a live-load test and detailed FE analyses. DFs for AASHTO LRFD specification and actual bridge were compared. In addition, the effects of span length, skew and bridge deck thickness on the DFs were also investigated in the parametric study. From the analyses, the following conclusions were obtained:

  1. 1.

    Compared with the DFs in the FEM of this hollow slab bridge, AASHTO LRFD DFs for interior girder (1–18%) were less conservative than exterior girder (37–48%).

  2. 2.

    The AASHTO LRFD specification became less conservative (8%) when the loading vehicles were filled in the lane that the bridge can be divided to for multilane loaded cases.

  3. 3.

    AASHTO LRFD specification became less conservative as the span length increased for this hollow slab bridge. DFs of exterior girder increased 10–42% while of interior girder decreased 9–30% as the span length increased.

  4. 4.

    AASHTO LRFD specification became less conservative as the skew was greater than \(10^\circ\). DFs of FEM for interior girder were 1.4 times larger than the DFs in AASHTO LRFD specification when the skew was \(60^\circ\), while the DFs of exterior girder were less than design value with an increase in the skew.

  5. 5.

    Deck thickness had little effects on the DFs of hollow slab bridge.

8 Discussion

The parapets have been proved to provide more stiffness to the exterior girder compared with the interior girders (Barker 2001; Zhou et al. 2015). As the vehicle loading test was conducted on the bridge before the parapet was constructed, the parapets were not modeled in the FEM. So the effect of the parapets on DFs was not discussed in this paper.