Integral bridge abutment with composite dowels: structural scheme and failure patterns

Zhihua Xiong, Di Di*, Hui Wang, Zhenhua PanCollege of Water Resources and Architectural Engineering, Northwest A&F University, Yangling, Chinaxiongzhihua_2013@126.com, http:// orcid.org/ 0000-0001-8796-10042022055844@nwafu.edu.cn, https://orcid.org/0009-0006-7465-966X18112096768@nwafu.edu.cn,2022050880@nwafu.edu.cn

Markus FeldmannInstitute of Steel Construction, RWTH Aachen University, Aachen, Germanyfeldmann@stb.rwth-aachen.de

INTRODUCTION

Integral bridges offer significant advantages over traditional bridges, including reduced construction costs by eliminating expansion joints and associated replacement. Additionally, they enhance bridge durability and safety, lower maintenance expenses, and significantly reduce lifecycle costs. However, the long-term shrinkage and creep of the concrete deck is inconsistent with the deformation of the steel girder, resulting in cumulative displacement at the abutment and even causing fatigue of the connector. At the meantime, the action of traffic load and temperature cycling also enhance the possibility of crack and fatigue of traditional connectors.
Integral bridges are generally categorized into three types: fully Integral Abutment Bridges (IABs), semi-integral abutment bridges, and extended deck slab bridges [1]. The economic and performance benefits of integral bridges are particularly notable [2]. In the United States, at least 27 states have recognized the cost-effectiveness of jointless bridges [3], which provide improved seismic resistance and superior performance under extreme conditions due to their enhanced structural redundancy. Seismic fragility curves and evaluations for the hybrid concrete-steel frame were reported in [4].
H-shaped steel piles are widely used in IABs due to their flexibility, high strength, and ability to accommodate horizontal deformation. Arsoy et al. [5] demonstrated that H-shaped steel piles are the optimal choice for IABs. In terms of soilstructure interactions of IABs under seismic action, parametric studies indicated that low friction and dilatancy angles in granular fill resulted in unfavorable settlement patterns [6]. Shaking table tests demonstrated that compressible inclusions between the abutment and backfill effectively reduce seismic accelerations transmitted to the bridge deck and mitigate settlements [7].
The girder-abutment joint is a critical component of composite girder integral bridges. In steel-concrete composite girders, shear connectors are essential for transmitting shear forces and preventing relative slip at the steel-concrete interface. This ensures effective composite action between the two materials. Kim et al. [8] investigated the bearing capacity of joints, demonstrating that installing welded nail connectors and stiffener connectors enhances joint stiffness, bearing capacity, and crack resistance.
Conventional shear connectors, such as stud connectors and perfobond strip connectors, are widely employed in steelconcrete composite girders. However, perfobond connectors encounter challenges related to rebar penetration. In contrast, composite dowels overcome these limitations by eliminating the need for specialized welding and delivering robust performance even with lower-strength materials [9], which stimulates this work using the composite dowels as connectors. With regarding to the fatigue performance of composite dowels, we [10-11] analyzed composite dowels crack propagation under flexural/pure shear, revealing the relationships between stress intensity factor and crack length.
The European Design Manual for Integral Composite Girder Bridges [12] proposes a joint configuration as shown in Fig.1(a) wherein steel girders are directly placed on steel piles. This design eliminates batch casting requirements, transmits shear forces via weld studs, and transfers compressive forces to the bridge abutment through a bottom flange bearing plate. Ashiduka et al. [13] introduced an alternative joint design as shown in Fig.1(b), utilizing direct holes in the girder's flange
and web to accommodate perforated steel bars. Riches et al. [14] proposed a configuration as shown in Fig.1(c) incorporating a bottom flange end plate to disperse concentrated loads, with tensile forces transmitted by steel reinforcement. The Japan Road Research Laboratory [15] conducted experimental and numerical studies on steel-concrete composite girder joints as shown in Fig.1(d), analyzing component-level load distribution. Briseghella et al. [16]] developed a design as shown in Fig.1(e) experimentally validated for hysteretic behavior, employing headless weld studs to transfer vertical shear forces, while tensile and compressive forces are transmitted through deck reinforcement and a bearing plate. Despite their merits, these configurations face practical challenges, including stud welding complexities, rebar penetration issues, concrete casting dead zones, and welding fatigue. To address these limitations, this study proposes a novel integral bridge abutment incorporating composite dowels girders as demonstrated in Fig 1(f), which is elaborated in the following sections.
Figure 1: Structural details of the joints.
In previous studies on typical steel girder-concrete abutment joints, two critical factors need attentions: the load-bearing capacity of the abutment joint and its ability to accommodate long-term thermal deformations. The joint must support vertical deck loads while maintaining sufficient flexibility to adapt to temperature-induced displacements [17-18]. Steelconcrete composite girders, when subjected to solar radiation and thermal cycling, experience uneven heat transfer, rendering the deck highly sensitive to temperature variations. Under extreme conditions, thermal effects may dominate over dead and live loads. Research indicates that temperature-induced elongation in jointless bridges can be predicted based on total temperature change [19], with bridge expansion length proportional to its total length, temperature variation, and coefficient of thermal expansion, typically distributed symmetrically at both ends [20].
However, in terms of steel-concrete composite dowel girder-abutment in this work, current integral bridge design codes don't cover. Furthermore, comprehensive concerns on the ultimate bending resistance, failure mechanisms, and load transfer modes of such systems-particularly those incorporating H-shaped steel pile flexible foundations under thermal effects remain addressing. Under this background, we propose an innovative integral abutment integrating composite dowels girder and connectors into the steel girder web, combined with H-shaped steel piles. This configuration enhances the interfacial bond between the concrete deck slab and steel girder while facilitating reinforcement bar penetration, thereby improving both construction efficiency and mechanical performance.

Integral bridge abutment with composite dowels

The novel integral abutment is consisted of Modified-Verbund-Fertigteil-Trager (MVFT) girder[10] and composite dowels connections, which is illustrated in Fig.1(f) and Fig.2. The open geometry of CLothoid (CL)-shaped connectors allows reinforcements to be directly threaded through predefined positions during construction, effectively resolving challenges such as rebar penetration difficulties and concrete casting dead zones. Composite dowels
are classified as ductile shear connectors [21], particularly suitable for steel-concrete composite sections without flanges. To enhance the bond between the steel pile and abutment, composite dowels are set in the tension zone of the H -shaped steel piles.
Figure 2: General configuration of the integral abutment.

Numerical study

The Finite Element (FE) model consists of five primary components: the concrete abutment, composite dowel girder, composite dowels, steel reinforcement, and H -shaped steel piles. The FE model, scaled at a ratio of 1:3, is shown in Fig. 2. The concrete abutment measures 800 mm in length and 500 mm in width along the bridge's longitudinal direction.
To account for contact interactions as well as material and geometric nonlinearities, the dynamic explicit method is adopted for the analysis. The quasi-static loading of the concrete abutment and steel girder is simulated using C3D8R elements, while the reinforcement is modeled with T3D2 elements. Given the complexity of the steel-concrete girder structure, localized mesh refinement is applied around the composite dowels to enhance accuracy, with a minimum element size of 10 mm .
For interaction definitions, the rebar is embedded within the abutment using the Embedded Region constraint. Since the composite dowels and the steel girder web can be fabricated from a single steel plate, they are modeled as a single integrated part. The interface between the steel reinforcement and surrounding concrete is also defined using the Embedded Region approach. The contact between the steel girder and concrete is modeled as hard contact in the normal direction, while tangential behavior is simulated using a penalty-based friction model. The 3D model of steel girder and abutment is demonstrated in Fig. 3.
Figure 3: Numeric model of integral abutment.

Verification model

To validate the accuracy of the FE model, a third-party experimental test was selected. Fig. 4 presents the elevation and side-view dimensions of the joint in the referenced test [21]. The front of the abutment includes a 30 mm thick pressure plate. The steel girder has a depth of 540 mm , with the upper and lower flanges and web each measuring 20 mm in thickness. The main girder is embedded 500 mm into the abutment. Near the front of the abutment, 50 mm diameter holes are provided in the flanges and web, spaced at 120 mm intervals. PBL connectors, consisting of 20 mm diameter steel bars, are inserted into these holes.
The verification model is developed based on the experimental setup, as shown in Fig. 4. Results indicate that the loaddisplacement response of the FE model closely matches the experimental data when the friction coefficient is set to 0.3 . In the model, all degrees of freedom at the constrained rigid reference point are restricted, except for rotation about the Zaxis. Displacement control is employed to prevent sudden failure, which could lead to excessive element deformation and affect ultimate load results.
To ensure consistency with the shear span ratio of the integral bridge and previous studies [22], a loading plate (modeled as an analytical rigid body) is positioned 1.75 times the deck depth from the abutment's front end and is rigidly connected to the deck. A displacement load is applied to the rigid body, with a prescribed forced displacement of 100 mm at the reference point of the loading plate to impose loading on the model. The ultimate load is defined as the maximum vertical reaction force of the loading plate during loading. A smooth amplitude curve is used for loading, and the dynamic internal energy ratio is controlled within 5 % 5 % 5%5 \%5% to ensure numerical stability.
For material modeling, the Concrete Damaged Plasticity model is used for concrete, with its constitutive behavior based on the General Design Code for Concrete Structures. Steel components are modeled using a bilinear strain-hardening constitutive model, following Eurocode 3 (2006). The material parameters are detailed in Tabs. 1 and 2.
Materials E / GPa E / GPa E//GPa\mathrm{E} / \mathrm{GPa}E/GPa ν / MPa ν / MPa nu//MPa\nu / \mathrm{MPa}ν/MPa σ s / MPa σ s / MPa sigma_(s)//MPa\sigma_{\mathrm{s}} / \mathrm{MPa}σs/MPa f t / MPa f t / MPa f_(t)//MPa\mathrm{f}_{\mathrm{t}} / \mathrm{MPa}ft/MPa f c / MPa f c / MPa f_(c)//MPa\mathrm{f}_{\mathrm{c}} / \mathrm{MPa}fc/MPa
C50 27.100 0.2 - - 50
Q345 206.000 0.3 345 500 -
Q460 206.000 0.3 460 600 -
HRB400 206.000 0.3 400 500 -
Materials E//GPa nu//MPa sigma_(s)//MPa f_(t)//MPa f_(c)//MPa C50 27.100 0.2 - - 50 Q345 206.000 0.3 345 500 - Q460 206.000 0.3 460 600 - HRB400 206.000 0.3 400 500 -| Materials | $\mathrm{E} / \mathrm{GPa}$ | $\nu / \mathrm{MPa}$ | $\sigma_{\mathrm{s}} / \mathrm{MPa}$ | $\mathrm{f}_{\mathrm{t}} / \mathrm{MPa}$ | $\mathrm{f}_{\mathrm{c}} / \mathrm{MPa}$ | | :---: | :---: | :---: | :---: | :---: | :---: | | C50 | 27.100 | 0.2 | - | - | 50 | | Q345 | 206.000 | 0.3 | 345 | 500 | - | | Q460 | 206.000 | 0.3 | 460 | 600 | - | | HRB400 | 206.000 | 0.3 | 400 | 500 | - |
Table 1: Elastic properties
Materials α / ( ) α / alpha//(^(@))\alpha /\left(^{\circ}\right)α/() Strain deviator f b 0 / f c 0 f b 0 / f c 0 f_(b0)//f_(c0)\mathrm{f}_{\mathrm{b} 0} / \mathrm{f}_{\mathrm{c} 0}fb0/fc0 K Coefficient of viscosity
Concrete 25 0.1 1.16 0.6667 0.0005
Materials alpha//(^(@)) Strain deviator f_(b0)//f_(c0) K Coefficient of viscosity Concrete 25 0.1 1.16 0.6667 0.0005| Materials | $\alpha /\left(^{\circ}\right)$ | Strain deviator | $\mathrm{f}_{\mathrm{b} 0} / \mathrm{f}_{\mathrm{c} 0}$ | K | Coefficient of viscosity | | :---: | :---: | :---: | :---: | :---: | :---: | | Concrete | 25 | 0.1 | 1.16 | 0.6667 | 0.0005 |
Table 2: Plastic properties.
Figure 4: Dimensions of girder-abutment test [21] (Unit: mm).
Concrete Elasticity modulus 31600
28-day cylinder f c f c f_(c)\mathrm{f}_{\mathrm{c}}fc 43.2
28-day cube f c f c f_(c)\mathrm{f}_{\mathrm{c}}fc 52.9
Elasticity modulus 198000
Yield strength 391
Ultimate strength 545
Concrete Elasticity modulus 31600 28-day cylinder f_(c) 43.2 28-day cube f_(c) 52.9 Elasticity modulus 198000 Yield strength 391 Ultimate strength 545| Concrete | Elasticity modulus | 31600 | | :---: | :---: | :---: | | | 28-day cylinder $\mathrm{f}_{\mathrm{c}}$ | 43.2 | | | 28-day cube $\mathrm{f}_{\mathrm{c}}$ | 52.9 | | | Elasticity modulus | 198000 | | | Yield strength | 391 | | | Ultimate strength | 545 |
HRB400 Elasticity modulus 202000
Yield strength 390
Ultimate strength 485
HRB400 Elasticity modulus 202000 Yield strength 390 Ultimate strength 485| HRB400 | Elasticity modulus | 202000 | | :---: | :---: | :---: | | | Yield strength | 390 | | | Ultimate strength | 485 |
Table 3: Material parameters (MPa).
Figure 5: Validation of FE results.
Figure 6: Comparison of failure modes in test and FE results.
As shown in Figs 5 and 6, the FE results closely match the experimental data from the third-party test. The numerical model exhibits slightly higher stiffness after entering the elastoplastic stage, which may be attributed to the absence of early-stage steel-concrete bonding in the ABAQUS model. In the model, the frictional force compensates for this bonding effect. In the experiment, the adhesive force disappears once the component enters the elastoplastic stage, while friction continues to influence the ABAQUS model, resulting in higher stiffness in the early and mid-stages of the FE analysis. Overall, the modeling approach is consistent with the experimental observations.
Fig. 8 illustrates the failure of the concrete in the numerical results, represented by the maximum plastic principal strain contour, which also aligns well with the experimental findings. The numerical model shows the steel girder's pull-out behavior, with cracks forming at the front of the platform and shear cracks along the side of the abutment, accurately simulating the observed failure modes.

Temperature effect

The minimum and maximum temperatures are key factors influencing the temperature effects on the structure. For the structure shown in Fig. 2, the temperature effect is specifically considered. Summer, being the hottest season of the year, exhibits the most significant temperature impact on bridges. To assess this, the monthly maximum and minimum temperature data for Xi'an over the past five years are analyzed, as shown in Fig. 7.
Correlation analysis and calculations were performed using MATLAB, with the Gevfit function applied to model the generalized extreme value (GEV) distribution. This approach analyzes the extreme temperature values to predict representative load values. The analysis assumes that the temperature samples are independent random variables, each following the same distribution. The generalized extreme value distribution is expressed by Eq. (1), where the representative values of the highest and lowest temperatures are 37.5 C 37.5 C 37.5^(@)C37.5^{\circ} \mathrm{C}37.5C and 15.7 C 15.7 C 15.7^(@)C15.7^{\circ} \mathrm{C}15.7C, respectively.
Figure 7: Monthly temperature data of Xi 'an in recent 5 years.
(1) H = exp { ( 1 + ζ x μ σ ) 1 / ζ } , 1 + ζ x μ σ > 0 (1) H = exp 1 + ζ x μ σ 1 / ζ , 1 + ζ x μ σ > 0 {:(1)H=exp{-(1+zeta(x-mu)/(sigma))^(-1//zeta)}","1+zeta(x-mu)/(sigma) > 0:}\begin{equation*} H=\exp \left\{-\left(1+\zeta \frac{x-\mu}{\sigma}\right)^{-1 / \zeta}\right\}, 1+\zeta \frac{x-\mu}{\sigma}>0 \tag{1} \end{equation*}(1)H=exp{(1+ζxμσ)1/ζ},1+ζxμσ>0
where x is the random variable, μ μ mu\muμ is the position parameter, σ σ sigma\sigmaσ is the scale parameter, ζ ζ zeta\zetaζ is the shape parameter, and the last three parameters are determined by the selected sample.
The maximum horizontal force of the temperature effect of the integral jointless bridge on the abutment and the following parameters are required according to [22]:
Abutment height H = 2.6 m H = 2.6 m H=2.6m\mathrm{H}=2.6 \mathrm{~m}H=2.6 m, deck width B = 1.6 m B = 1.6 m B=1.6m\mathrm{B}=1.6 \mathrm{~m}B=1.6 m, span length L = 20 m L = 20 m L=20m\mathrm{L}=20 \mathrm{~m}L=20 m, unit weight γ = 18 kN / m 3 , T max = 37.5 C , T min = 15.7 C γ = 18 kN / m 3 , T max = 37.5 C , T min = 15.7 C gamma=18kN//m^(3),T_(max)=37.5^(@)C,T_(min)=15.7^(@)C\gamma=18 \mathrm{kN} / \mathrm{m}^{3}, \mathrm{~T}_{\max }=37.5^{\circ} \mathrm{C}, \mathrm{T}_{\min }=15.7^{\circ} \mathrm{C}γ=18kN/m3, Tmax=37.5C,Tmin=15.7C, Coefficient of thermal expansion of concrete α = 10 × 10 6 / C α = 10 × 10 6 / C alpha=10 xx10^(-6)//^(@)C\alpha=10 \times 10^{-6} /{ }^{\circ} \mathrm{C}α=10×106/C. If calculated as beginning from zero, the horizontal spring for the whole abutment would be:
(2) k = p Δ L = H 2 γ 1.3 H B / 0.003 = 16224 kN m (3) Δ L k + F = 372.7 kN (2) k = p Δ L = H 2 γ 1.3 H B / 0.003 = 16224 kN m (3) Δ L k + F = 372.7 kN {:[(2)k=(p)/(Delta L)=(H)/(2)*gamma*1.3*H*B//0.003=16224kN*m],[(3)Delta L*k+F=372.7kN]:}\begin{align*} & k=\frac{p}{\Delta L}=\frac{H}{2} \cdot \gamma \cdot 1.3 \cdot H \cdot B / 0.003=16224 \mathrm{kN} \cdot m \tag{2}\\ & \Delta L \cdot k+F=372.7 \mathrm{kN} \tag{3} \end{align*}(2)k=pΔL=H2γ1.3HB/0.003=16224kNm(3)ΔLk+F=372.7kN

Parametric analysis

Aparametric analysis was conducted to examine the influence of various abutment parameters, including steel type, plate thickness, bridge deck reinforcement ratio and composite dowel spacing. Based on the General Design Code for Concrete Structures (GB 55008-2021), the parameters were set as follows:
Steel type: Q345, Q460; Steel plate thickness (t): 16-24 mm; Composite dowel spacing ( e x e x e_(x)\mathrm{e}_{\mathrm{x}}ex ): 200-800 mm;
Bridge deck reinforcement ratio ( @ 1 @ 1 @_(1)@_{1}@1 ): 1.01 % , 1.52 % , 2.11 % 1.01 % , 1.52 % , 2.11 % 1.01%,1.52%,2.11%1.01 \%, 1.52 \%, 2.11 \%1.01%,1.52%,2.11%; Deck width ( w s w s w_(s)w_{s}ws ): 400 600 mm 400 600 mm 400-600mm400-600 \mathrm{~mm}400600 mm; Abutment longitudinal width ( w c ) : 750 1500 mm ( w c ) : 750 1500 mm (wc):750-1500mm(\mathrm{w} c): 750-1500 \mathrm{~mm}(wc):7501500 mm; Abutment reinforcement ratio ( @ 2 @ 2 @_(2)@_{2}@2 ): 0.8 % 1.2 % 0.8 % 1.2 % 0.8%-1.2%0.8 \%-1.2 \%0.8%1.2%.
To validate the parametric study, a pilot integral abutment with composite dowels was designed and analyzed. The span of the integral abutment bridge with composite dowel girders is L = 20 m L = 20 m L=20m\mathrm{L}=20 \mathrm{~m}L=20 m. The deck consists of four girders, each B = 1.6 m B = 1.6 m B=1.6m\mathrm{B}=1.6 \mathrm{~m}B=1.6 m wide, as illustrated in Fig. 8. The total deck width is 8.5 m , while the concrete abutment has a longitudinal width of 0.9 m and a height of 2.6 m . The H-shaped steel piles are of type HW 400 × 400 × 20 × 20 mm 400 × 400 × 20 × 20 mm 400 xx400 xx20 xx20mm400 \times 400 \times 20 \times 20 \mathrm{~mm}400×400×20×20 mm.
A three-dimensional finite element model of the integral abutment joint, composite girder, and H -shaped steel piles was developed using ABAQUS. The load-bearing performance and failure mechanisms of the integrated joint were then systematically investigated. Here, H represents the composite girder height, while h denotes the steel girder height.
Figure 8: Cross-section of composite dowels girder (mm).

Influence of parameters of composite dowels girder

As shown in Fig. 9, increasing the steel strength of the composite dowel girder from Q345 to Q460, while keeping other parameters constant, results in a slight improvement in ultimate bearing capacity. For a girder height of 800 mm , the ultimate bearing capacity increases from 559 kN to 592 kN , representing a 5.9 % 5.9 % 5.9%5.9 \%5.9% increase. For a 1000 mm girder, the capacity rises from 652 kN to 675 kN , an increase of 3.5 % 3.5 % 3.5%3.5 \%3.5%. These results indicate that the ultimate bearing capacity of high-depth steel-concrete composite girders with composite dowels is minimally affected by steel strength. Therefore, a more economical steel grade can be selected without significantly compromising the load-bearing performance. As illustrated in Fig. 10, the bearing capacity of the steel girder base plate increases by 2.4 % 2.4 % 2.4%2.4 \%2.4% when its width expands from 400 mm to 500 mm , and by 12 % 12 % 12%12 \%12% when further increased to 600 mm . According to design specifications, the base plate width can be set at either 400 mm or 600 mm , depending on load-bearing requirements. A 400 mm width is suitable for lower design standards, whereas a 600 mm width is preferable for higher load-bearing demands.
As depicted in Fig. 11, for girder depths of 800 mm and 1000 mm , the ultimate bearing capacity does not increase with changes in web thickness or composite dowel spacing. The failure mode analysis reveals that the structure ultimately loses its load-bearing capacity due to concrete crushing beneath the floor, rather than limitations in the steel girder web.
Figure 9: Effect of the steel strength of composite dowels girder.
Figure 10: Effect of the width of steel girder bottom plate.
Figure 11: Effect of the center distance of composite dowels(a) and thickness of steel girder (b).
Three composite dowel arrangements were evaluated for embedding into the steel girder bottom plate at the bridge end: CL-1: Dowel spacing e x = 200 mm e x = 200 mm e_(x)=200mm\mathrm{e}_{\mathrm{x}}=200 \mathrm{~mm}ex=200 mm, with two rows spaced 260 mm apart at the girder end.
CL-2: Dowel spacing e x = 200 mm e x = 200 mm e_(x)=200mm\mathrm{e}_{\mathrm{x}}=200 \mathrm{~mm}ex=200 mm, with a horizontal arrangement of three rows, spaced 100 mm apart at the girder end.
CL-3: No composite dowels arranged.
Tab. 4 compares the bearing capacities of these configurations, showing that CL-2 outperforms the other arrangements. Compared to CL-1, CL-2 improves bearing capacity by 2.3 % 2.3 % 2.3%2.3 \%2.3%, while compared to CL-3, it increases by 6.1 % 6.1 % 6.1%6.1 \%6.1%. Despite featuring three rows of dowels, CL-2 does not compromise construction feasibility, making it the recommended configuration.
Table 4: Composite dowels combinations.
Figure 12: Effect of reinforcement ratio of abutment
Figure 13: Effect of reinforcement ratio of deck ϱ 2 ϱ 2 ϱ_(2)\varrho_{2}ϱ2

Influence of reinforcement ratio

Fig. 12 shows the effect of the reinforcement ratio of the abutment ( Q 1 Q 1 Q_(1)\boldsymbol{Q}_{\mathbf{1}}Q1 ) on the joint's load-bearing capacity. The results indicate a significant correlation between ϱ 1 ϱ 1 ϱ_(1)\varrho_{1}ϱ1 and the structural performance of the joint. Specifically, when ϱ 1 ϱ 1 ϱ_(1)\varrho_{1}ϱ1 increases from 0.7 % 0.7 % 0.7%0.7 \%0.7% to 1.5 % 1.5 % 1.5%1.5 \%1.5%, the ultimate bearing capacity of the joint improves by approximately 11 % 11 % 11%11 \%11%.
Similarly, Fig. 13 presents the influence of the reinforcement ratio of the deck ( ϱ 2 ϱ 2 ϱ_(2)\varrho_{2}ϱ2 ) on the joint's load-bearing capacity. The analysis reveals that when ϱ 2 ϱ 2 ϱ_(2)\varrho_{2}ϱ2 increases from 1 % 1 % 1%1 \%1% to 2 % 2 % 2%2 \%2%, the joint's bearing capacity enhances by approximately 13 % 13 % 13%13 \%13%. These findings highlight the crucial role of reinforcement in improving the structural integrity and performance of the integral abutment joint.

Influence of longitudinal bridge width on abutment W C W C W_(C)W_{C}WC

A s s sss is shown in Fig.14, the longitudinal bridge width ( 750 1500 mm 750 1500 mm 750-1500mm750-1500 \mathrm{~mm}7501500 mm ) exhibits a linear relationship with the overall bearing capacity of the joint. When the width increases from 750 mm to 1500 mm (a 100 % 100 % 100%100 \%100% increase), the bearing capacity improves by 35.8 % 35.8 % 35.8%35.8 \%35.8%. These findings suggest that the longitudinal bridge width of abutments can be optimized based on the specific requirements of the project to achieve the desired structural performance.
Figure 14: Effect of longitudinal bridge width on abutment.

FAILURE PATTERNS AND LOAD-TRANSITION MECHANISM - FAILURE PATTERNS

The failure modes of the joint are illustrated in Fig.s 15-18. In the finite element results, concrete failure is represented by the maximum plastic principal strain contour, where white areas indicate severe cracking. The steel girder failure is depicted using the longitudinal strain contour of the bridge.
As shown in Fig. 15, the stress distribution and 3 × 3 × 3xx3 \times3× deformation amplification of the steel girder at failure reveal that the web experiences significant strain, while the bottom plate and ends of the web undergo compression. Due to high stress, the bottom plate near the front of the abutment deforms significantly. During loading, the neutral axis of the steel girder section gradually moves downward, leading to rotation around the bottom plate near the front of the abutment. This results in cracks forming in the concrete web. The web plate of the embedded steel girder section is primarily subjected to shear forces, exhibiting significant shear stress. The relative angle between the inserted section and the extended section of the steel girder indicates strong interaction with the concrete abutment, demonstrating the high connection performance and tensile resistance of the composite dowels. Throughout the loading process, the integral abutment joints are subjected to large bending moments and shear forces, ultimately causing cracks in the concrete web due to excessive tensile and compressive stresses at the steel-concrete interface. As the external load increases, the number and width of cracks progressively expand.
Figure 15: Stress and deformation of steel main girder (MPa).
Figure 16: Plastic strain contour of steel girder.
Fig. 16 presents the plastic strain contours of steel girders for ex = 200 mm = 200 mm =200mm=200 \mathrm{~mm}=200 mm and ex = 600 mm = 600 mm =600mm=600 \mathrm{~mm}=600 mm when the joint reaches its ultimate bearing capacity, highlighting the failure modes. As ex increases, the plastic region at the interface between the embedded steel girder section and the extended section expands from the web to the bottom plate, creating a weakened section. Additionally, the tooth joints yield, and local crushing of concrete under the steel girder occurs once the weakened section reaches full plasticity.
According to the failure patterns of the concrete component shown in Fig. 17, tensile cracks initially develop at the top of the rear abutment due to its high stiffness and the rotation of the extended steel girder section. As loading continues, these cracks propagate, particularly in the concrete surrounding the composite dowels, caused by the rotation of the steel girder and concrete compression between the dowels. These cracks gradually extend toward the deck. With further loading, tensile stresses induce cracking in the deck, which spreads at an inclined angle toward the junction of the main girder and steel piles. The cracks on the abutment are shown in Fig. 18.
Figure 17: Failure mode.
Figure 18: Crack distribution.
During loading, the bending moment and shear force in the steel girder exceed its capacity, causing plastic deformation and yielding near the front of the abutment (Fig. 16). Eventually, concrete crushing beneath the steel girder base plate leads to structural failure. This phenomenon highlights the importance of considering steel-concrete interaction and their loadbearing properties under various stress states when designing composite girder structures.
Figure 19: Load-displacement curve.
To mitigate dynamic noise interference, the load-displacement curve was filtered (Fig. 19). After processing, the yield and ultimate points are more distinctly visible and align well with the original data. The yield point is determined as the tangent intersection between the line connecting the extreme point to the origin and the load-displacement curve.
In general, at the initial stage of loading, the abutment exhibits insignificant deformation, remaining within the elastic stage until reaching the yield point. As displacement increases, the bearing capacity gradually declines post-peak, forming a plateau region. This behavior occurs due to collapse of the concrete between composite dowels, surrounding concrete cracking, and shear resistance loss, while the penetrating steel bars between composite dowels continue to resist shear forces.
After reaching the peak load, the bearing capacity does not sharply decline, but deformation accelerates, demonstrating good ductility characteristics. The main girder concrete, abutment cracks, steel girder deflection, and reinforcement deformation progressively increase. Ultimately, the structure loses its load-bearing capacity due to: 1) Concrete crushing beneath the steel girder base plate; 2) Concrete web cracking; 3) Steel girder buckling.

LOAD-TRANSFERRING MECHANISM

As shown in Fig. 20, the stress on the composite dowels within the steel girder is minimal and does not reach the yield state. The stress beneath the composite dowels on the steel web is generally lower than that between the composite dowels. The composite dowels on the web do not directly contribute to the bearing capacity but enhance the collaborative load-bearing effect between the steel girder and the concrete. However, as the spacing between the dowels decreases, the stress on the composite dowels increases, although it remains significantly smaller than the stress on the steel girder web.
Fig. 20 and Fig. 21 demonstrate that the top section of the concrete abutment is damaged due to the action of the negative bending moment at the joint, where the steel girder acts as a lever. The concrete in the damaged area experiences a combination of tensile yield in the upper part and compression in the lower part. The external load is primarily resisted by three components of the joint: the concrete between the composite dowels, the steel bars in the concrete deck slab, and the concrete at the contact surface of the main girder. In the elastic stage, the bending moment from the external load is primarily borne by the steel girder's bottom plate. As the load increases, the steel web and joint enter the elastoplastic stage, the relative angle increases rapidly, and the load is gradually transferred to the steel bars in the deck.
When the plastic region of the steel girder web reaches its maximum, the joint enters the failure stage. The contact force between the steel girder bottom plate and the concrete decreases, and the load is predominantly borne by the concrete beneath the steel girder. With further load increase, the concrete beneath the steel girder is crushed, leading to the failure of the joint.
Figure 20: Local stress of composite dowels at failure
Figure 22: Stress of concrete and rebars between composite dowels.
Figure 21: Local strain of concrete at failure
Figure 23: Force transmission of nodes
Fig. 22 shows the maximum principal stress distribution in the concrete between the composite dowels under the ultimate load. It can be observed that the concrete in the hole experiences hydro-compression due to the constraints imposed by the surrounding concrete blocks and the presence of penetrating rebars. Therefore, during operation, the rebars within the composite dowels primarily experience shear stress, while the concrete between the dowels is subject to compressive stress. Based on Fig. 17 and Fig. 23, the primary stress transfer occurs at the contact point between the steel girder and the concrete beneath the plate. The shear force at the front surface of the platform, denoted as V , represents the joint's shear capacity. The forces F 1 F 1 F_(1)\mathrm{F}_{1}F1 and F 2 F 2 F_(2)\mathrm{F}_{2}F2 are the equivalent concentrated forces acting in the compression area of the steel girder, while M denotes the bending moment of the steel girder at the front surface of the platform, indicating the joint's bending capacity. C 1 C 1 C_(1)\mathrm{C}_{1}C1 and C 2 C 2 C_(2)\mathrm{C}_{2}C2 represent the positions of the equivalent concentrated forces in the compression area of the steel girder.

Theoretical analyses of the flexural bearing capacity of joints - Assumptions

Based on the theoretical analysis of the flexural bearing capacity of the joint of steel girder and concrete abutment, the following assumptions are adopted herein:
  1. Regardless of the adhesion and friction between the steel girder and the concrete abutment; 2) The assumption that the buried steel girder is regarded as a rigid body which is often used in the analysis of joints without connectors. When the connectors are set, the assumption is more realistic because the buried steel girder is shorter and has higher stiffness. Under the above assumptions, the compressive stress at the top and bottom of the steel girder is linearly distributed, and the peak strain ε c l ε c l epsi_(cl)\varepsilon_{c l}εcl and ultimate strain ε c u ε c u epsi_(cu)\varepsilon_{c u}εcu are 0.002 and 0.003 , respectively.
(4) ε b = 0.03 ( l e c ) / c (4) ε b = 0.03 ( l e c ) / c {:(4)epsi_(b)=0.03(le-c)//c:}\begin{equation*} \varepsilon_{b}=0.03(l e-c) / c \tag{4} \end{equation*}(4)εb=0.03(lec)/c
The force diagram of the joint in the limit state is demonstrated in Fig.24. The failure mode of the joint is assumed to be crushed concrete under the steel girder, and the other failure modes are not discussed this time. This requires ε b < 0.003 ε b < 0.003 epsi_(b) < 0.003\boldsymbol{\varepsilon}_{b}<0.003εb<0.003 in the limiting state. Under this condition, ε b ε b epsi_(b)\varepsilon_{b}εb of alternative Eq. (4) can obtain c / l e 0.5 c / l e 0.5 c//l_(e)◻0.5c / l_{e} \square 0.5c/le0.5.
Figure 24: Force transferring diagram in the limit state.
The stress-strain relationship of concrete is assumed to be:
(5) σ c = { 1000 f c ( ε 250 ε 2 ) , 0 < ε ε 1 f c , ε 1 < ε ε c u (5) σ c = 1000 f c ε 250 ε 2 , 0 < ε ε 1 f c , ε 1 < ε ε c u {:(5)sigma_(c)={[1000f_(c)(epsi-250epsi^(2))",",0 < epsi <= epsi_(1)],[f_(c),","epsi_(1) < epsi <= epsi_(cu)]:}:}\sigma_{c}= \begin{cases}1000 f_{c}\left(\varepsilon-250 \varepsilon^{2}\right), & 0<\varepsilon \leqslant \varepsilon_{1} \tag{5}\\ f_{c} & , \varepsilon_{1}<\varepsilon \leqslant \varepsilon_{c u}\end{cases}(5)σc={1000fc(ε250ε2),0<εε1fc,ε1<εεcu
The bending capacity of the steel girder-concrete joint is assumed to be the sum of the bending moment contribution of the bearing action and the bending moment contribution of the connector, namely:
(6) M u = P u ( a + c ) = M c + M s (6) M u = P u ( a + c ) = M c + M s {:(6)M_(u)=P_(u)(a+c)=M_(c)+M_(s):}\begin{equation*} M_{u}=P_{u}(a+c)=M_{c}+M_{s} \tag{6} \end{equation*}(6)Mu=Pu(a+c)=Mc+Ms

Determination of the rotation center

As shown in Fig. 24, the horizontal and vertical point lines are the parts where the vertical and horizontal relative displacements of steel girder and concrete abutment are zero respectively, and the center of deformation and rotation O is defined as the intersection point of these two lines. The balance equation combining vertical direction and horizontal direction in the ministry is as follows:
(7) C b + i = 1 n v i y C f + P u = 0 (7) C b + i = 1 n v i y C f + P u = 0 {:(7)C_(b)+sum_(i=1)^(n)v_(i)^(y)-C_(f)+P_(u)=0:}\begin{equation*} C_{b}+\sum_{i=1}^{n} v_{i}^{y}-C_{f}+P_{u}=0 \tag{7} \end{equation*}(7)Cb+i=1nviyCf+Pu=0
Make λ = c / l e , μ = a / l e λ = c / l e , μ = a / l e lambda=c//l_(e),mu=a//l_(e)\lambda=c / l_{e}, \mu=a / l_{e}λ=c/le,μ=a/le
Previous studies have concluded that:
(8) M C = C b ( 1 k 0 ) ( l e c ) + C f ( c 1 2 β 1 c ) (8) M C = C b 1 k 0 l e c + C f c 1 2 β 1 c {:(8)M_(C)=C_(b)(1-k_(0))(l_(e)-c)+C_(f)(c-(1)/(2)beta_(1)c):}\begin{equation*} M_{C}=C_{b}\left(1-k_{0}\right)\left(l_{e}-c\right)+C_{f}\left(c-\frac{1}{2} \beta_{1} c\right) \tag{8} \end{equation*}(8)MC=Cb(1k0)(lec)+Cf(c12β1c)
When λ 0.6 λ 0.6 lambda >= 0.6\lambda \geq 0.6λ0.6
(9) C b = 3 4 f c b l e ( 1 λ λ ) 2 ( 3 λ 1 ) (9) C b = 3 4 f c b l e 1 λ λ 2 ( 3 λ 1 ) {:(9)C_(b)=(3)/(4)f_(c)bl_(e)((1-lambda)/(lambda))^(2)(3lambda-1):}\begin{equation*} C_{b}=\frac{3}{4} f_{c} b l_{e}\left(\frac{1-\lambda}{\lambda}\right)^{2}(3 \lambda-1) \tag{9} \end{equation*}(9)Cb=34fcble(1λλ)2(3λ1)
(10) C f = 0.85 β 1 f c b l e λ (11) k = 2.75 λ 0.75 9 λ 3 (10) C f = 0.85 β 1 f c b l e λ (11) k = 2.75 λ 0.75 9 λ 3 {:[(10)C_(f)=0.85beta_(1)f_(c)bl_(e)lambda],[(11)k=(2.75 lambda-0.75)/(9lambda-3)]:}\begin{align*} C_{f} & =0.85 \beta_{1} f_{c} b l_{e} \lambda \tag{10}\\ k & =\frac{2.75 \lambda-0.75}{9 \lambda-3} \tag{11} \end{align*}(10)Cf=0.85β1fcbleλ(11)k=2.75λ0.759λ3
When 0.5 λ 0.6 0.5 λ 0.6 0.5◻lambda <= 0.6◻0.5 \square \lambda \leq 0.6 \square0.5λ0.6
(12) C b = f c b l e ( 1 11 9 λ ) (13) C f = 0.85 β 1 f c b l e λ (14) k = 13.5 33 λ + 20.5 λ 2 27 60 λ + 33 λ 2 (12) C b = f c b l e 1 11 9 λ (13) C f = 0.85 β 1 f c b l e λ (14) k = 13.5 33 λ + 20.5 λ 2 27 60 λ + 33 λ 2 {:[(12)C_(b)=f_(c)bl_(e)(1-(11)/(9)lambda)],[(13)C_(f)=0.85beta_(1)f_(c)bl_(e)lambda],[(14)k=(13.5-33 lambda+20.5lambda^(2))/(27-60 lambda+33lambda^(2))]:}\begin{align*} & C_{b}=f_{c} b l_{e}\left(1-\frac{11}{9} \lambda\right) \tag{12}\\ & C_{f}=0.85 \beta_{1} f_{c} b l_{e} \lambda \tag{13}\\ & k=\frac{13.5-33 \lambda+20.5 \lambda^{2}}{27-60 \lambda+33 \lambda^{2}} \tag{14} \end{align*}(12)Cb=fcble(1119λ)(13)Cf=0.85β1fcbleλ(14)k=13.533λ+20.5λ22760λ+33λ2
Subdivision point λ = 0.6 λ = 0.6 lambda=0.6\lambda=0.6λ=0.6, corresponding to Eq. (3) in ε = ε μ ε = ε μ epsi=epsi_(mu)\boldsymbol{\varepsilon}=\boldsymbol{\varepsilon}_{\mu}ε=εμ l.
It can be seen that λ λ lambda\lambdaλ is related to the bending moment contribution of the CL and the arrangement form of the CL. The solution for λ λ lambda\lambdaλ is very complex, but it can be inferred that the relative rotation of the center position O at the centroid of the composite dowels and the contribution to the bending moment of the composite dowels can be ignored without timing being obtained: When λ 0.6 λ 0.6 lambda >= 0.6\lambda \geq 0.6λ0.6,
( 55 34 β 1 2 ) λ 4 + ( 180 68 β 1 ) μ λ 3 + ( 420 μ 210 ) λ 2 + (15) ( 330 μ + 200 ) λ ( 60 μ + 45 ) = 0 55 34 β 1 2 λ 4 + 180 68 β 1 μ λ 3 + ( 420 μ 210 ) λ 2 + (15) ( 330 μ + 200 ) λ ( 60 μ + 45 ) = 0 {:[(55-34beta_(1)^(2))lambda^(4)+(180-68beta_(1))mulambda^(3)+(-420 mu-210)lambda^(2)+],[(15)(330 mu+200)lambda-(60 mu+45)=0]:}\begin{gather*} \left(55-34 \beta_{1}^{2}\right) \lambda^{4}+\left(180-68 \beta_{1}\right) \mu \lambda^{3}+(-420 \mu-210) \lambda^{2}+ \\ (330 \mu+200) \lambda-(60 \mu+45)=0 \tag{15} \end{gather*}(5534β12)λ4+(18068β1)μλ3+(420μ210)λ2+(15)(330μ+200)λ(60μ+45)=0
When 0.5 λ 0.6 0.5 λ 0.6 0.5◻lambda <= 0.6◻0.5 \square \lambda \leq 0.6 \square0.5λ0.6
(16) ( 17 40 β 1 2 + 41 54 ) λ 2 + ( 11 9 + 17 20 β 1 ) μ λ μ 0.5 = 0 (16) 17 40 β 1 2 + 41 54 λ 2 + 11 9 + 17 20 β 1 μ λ μ 0.5 = 0 {:(16)((17)/(40)beta_(1)^(2)+(41)/(54))lambda^(2)+((11)/(9)+(17)/(20)beta_(1))mu lambda-mu-0.5=0:}\begin{equation*} \left(\frac{17}{40} \beta_{1}^{2}+\frac{41}{54}\right) \lambda^{2}+\left(\frac{11}{9}+\frac{17}{20} \beta_{1}\right) \mu \lambda-\mu-0.5=0 \tag{16} \end{equation*}(16)(1740β12+4154)λ2+(119+1720β1)μλμ0.5=0
Mattock [23] assumed that μ μ mu\muμ is between 0.1-0.5, and calculated the relationship through Eq. (13). λ λ lambda\lambdaλ deviated little from its mean value of 0.66 , so λ λ lambda\lambdaλ was taken as 0.66 in the design. Liang found that for the joints of integral bridges, however, μ μ mu\muμ between 0.1-0.5 is lower and should be between 1-10, and the value of β β beta\betaβ ranges from 0.65-085. Corresponding to the compressive strength of concrete axis 28 56 MPa 28 56 MPa 28-56MPa28-56 \mathrm{MPa}2856MPa [21]. The mean is 0.56 as the value of λ λ lambda\lambdaλ, and the standard deviation is 4.4 % 4.4 % 4.4%4.4 \%4.4% Therefore, the coordinates relative to the center of rotation O are:
(17) x 0 = 0.44 l e (18) y 0 = 1 n i = 1 n y i (17) x 0 = 0.44 l e (18) y 0 = 1 n i = 1 n y i {:[(17)x_(0)=0.44l_(e)],[(18)y_(0)=(1)/(n)sum_(i=1)^(n)y_(i)]:}\begin{align*} x_{0} & =0.44 l_{e} \tag{17}\\ y_{0} & =\frac{1}{n} \sum_{i=1}^{n} y_{i} \tag{18} \end{align*}(17)x0=0.44le(18)y0=1ni=1nyi

Contribution from girder

(19) (19) {:(19)◻:}\begin{equation*} \square \tag{19} \end{equation*}(19)
ince λ = 0.56 λ = 0.56 lambda=0.56\lambda=0.56λ=0.56 and 0.5 λ 0.6 , C b , C f 0.5 λ 0.6 , C b , C f 0.5◻lambda <= 0.6,C_(b),C_(f)0.5 \square \lambda \leq 0.6, C_{b}, C_{f}0.5λ0.6,Cb,Cf and k k kkk in Eq.(6) are replaced by Eqn. (10)-(12)
M c = 0.85 f c b l e 2 ( 0.314 β 1 0.157 β 1 2 + 0.100 ) M c = 0.85 f c b l e 2 0.314 β 1 0.157 β 1 2 + 0.100 M_(c)=0.85f_(c)bl_(e)^(2)(0.314beta_(1)-0.157beta_(1)^(2)+0.100)M_{c}=0.85 f_{c} b l_{e}^{2}\left(0.314 \beta_{1}-0.157 \beta_{1}^{2}+0.100\right)Mc=0.85fcble2(0.314β10.157β12+0.100)
The maximum local stress f b f b f_(b)f_{b}fb of concrete under the limit state is assumed to be 0.85 f c 0.85 f c 0.85f_(c)0.85 f_{c}0.85fc in the above calculation, but the actual f b f b f_(b)f_{b}fb is related to the ratio of the width of the concrete abutment to the width of the steel girder ( w / b ) f b ( w / b ) f b (w//b)*f_(b)(w / b) \cdot f_{b}(w/b)fb can be expressed as:
(20) f b = K ( f c ) a ( w b ) c (20) f b = K f c a w b c {:(20)f_(b)=K(f_(c))^(a)((w)/(b))^(c):}\begin{equation*} f_{b}=K\left(f_{c}\right)^{a}\left(\frac{w}{b}\right)^{c} \tag{20} \end{equation*}(20)fb=K(fc)a(wb)c
Mattock[23] Tschemmernegg[24] Eurocode 2: -part 1-1 Hawkins[25]
K 4.5 4.5 1.0 1.0 ◻1.0\square 1.01.0 4.5
a 0.5 0.5 1.0 0.5
c 0.66 0.55/0.60 0.50 0.33
Range w / b 12 w / b 12 w//b <= 12w / b \leq 12w/b12 w / b 4.5 w / b 4.5 w//b <= 4.5w / b \leq 4.5w/b4.5 f b 3.0 f c f b 3.0 f c f_(b) <= 3.0f_(c)f_{b} \leq 3.0 f_{c}fb3.0fc -
Mattock[23] Tschemmernegg[24] Eurocode 2: -part 1-1 Hawkins[25] K 4.5 4.5 ◻1.0 4.5 a 0.5 0.5 1.0 0.5 c 0.66 0.55/0.60 0.50 0.33 Range w//b <= 12 w//b <= 4.5 f_(b) <= 3.0f_(c) -| | Mattock[23] | Tschemmernegg[24] | Eurocode 2: -part 1-1 | Hawkins[25] | | :--- | :--- | :--- | :--- | :--- | | K | 4.5 | 4.5 | $\square 1.0$ | 4.5 | | a | 0.5 | 0.5 | 1.0 | 0.5 | | c | 0.66 | 0.55/0.60 | 0.50 | 0.33 | | Range | $w / b \leq 12$ | $w / b \leq 4.5$ | $f_{b} \leq 3.0 f_{c}$ | - |
Table 5: Summarizes the coefficient values of f b f b f_(b)f_{b}fb calculated in the literature
In this paper, w / b = 4 w / b = 4 w//b=4w / b=4w/b=4 is selected based on the biased safety angle:
(21) f b = 4.5 f c ( w / b ) 0.66 (21) f b = 4.5 f c ( w / b ) 0.66 {:(21)f_(b)=4.5sqrt(f_(c))(w//b)^(0.66):}\begin{equation*} f_{b}=4.5 \sqrt{f_{c}}(w / b)^{0.66} \tag{21} \end{equation*}(21)fb=4.5fc(w/b)0.66
replacing 0.85 f c 0.85 f c 0.85f_(c)0.85 f_{c}0.85fc in Eq. (17) with Eq. (19), when β 1 β 1 beta_(1)\beta_{1}β1 varies between 0.65 0.85 , 2 β 1 β 1 2 + 0.637 0.65 0.85 , 2 β 1 β 1 2 + 0.637 0.65-0.85,2beta_(1)-beta_(1)^(2)+0.6370.65-0.85,2 \beta_{1}-\beta_{1}{ }^{2}+0.6370.650.85,2β1β12+0.637 is the mean value of 1.5645 , and the maximum error is 3.3 % 3.3 % 3.3%3.3 \%3.3%.
(22) M C = 1.106 l e 2 f c b w (22) M C = 1.106 l e 2 f c b w {:(22)M_(C)=1.106l_(e)^(2)sqrt(f_(c)bw):}\begin{equation*} M_{C}=1.106 l_{e}^{2} \sqrt{f_{c} b w} \tag{22} \end{equation*}(22)MC=1.106le2fcbw

Contribution from composite dowels

T he bearing capacity of a single composite dowel is determined according to [11], namely:
(23) F e = 2.6 η s k f c m c D , S t h d , e f f 3 e x 3 (23) F e = 2.6 η s k f c m c D , S t h d , e f f 3 e x 3 {:(23)F_(e)=2.6*eta_(s)*k*sqrt(f_(cm)c_(D,S))*root(3)(th_(d,eff))*root(3)(e_(x)):}\begin{equation*} F_{e}=2.6 \cdot \eta_{s} \cdot k \cdot \sqrt{f_{c m} c_{D, S}} \cdot \sqrt[3]{t h_{d, e f f}} \cdot \sqrt[3]{e_{x}} \tag{23} \end{equation*}(23)Fe=2.6ηskfcmcD,Sthd,eff3ex3
The symbols are described as:
η s η s eta_(s)\eta_{s}ηs reinforcement stirrup coefficient, 1 stirrups: η s = 0.9 , 2 η s = 0.9 , 2 eta_(s)=0.9,2\eta_{s}=0.9,2ηs=0.9,2 stirrups: η s = 1 , 3 η s = 1 , 3 eta_(s)=1,3\eta_{s}=1,3ηs=1,3 stirrups: η s = 1.1 ; η s = 1.1 ; eta_(s)=1.1;\eta_{s}=1.1 ;ηs=1.1;
k effect k effect  k_("effect ")k_{\text {effect }}keffect  size effect k = 1 + ( 200 / c D , S ) 0.5 k = 1 + 200 / c D , S 0.5 k^(-)=1+(200//c_(D,S))^(0.5)k^{-}=1+\left(200 / c_{D, S}\right)^{0.5}k=1+(200/cD,S)0.5;
f c m f c m f_(cm)f_{c m}fcm compressive strength of the concrete (MPa);
c D , S c D , S c_(D,S)c_{D, S}cD,S compressive strength of the concrete (m);
t t t quadt \quadt steel thickness (m);
e x e x e_(x)quade_{x} \quadex center distance of the composite dowels (m);
h d,eff h d,eff  h_("d,eff ")h_{\text {d,eff }}hd,eff  effective steel tooth height (m).
The bearing capacity of all composite dowels joints in the joint is assumed to be a polynomial relationship between F s F s F_(s)F_{s}Fs and F e F e F_(e)F_{e}Fe :
(24) F S = F N , F E M F C (24) F S = F N , F E M F C {:(24)F_(S)=F_(N,FEM)-F_(C):}\begin{equation*} F_{S}=F_{N, F E M}-F_{C} \tag{24} \end{equation*}(24)FS=FN,FEMFC

Contribution from reinforcement ratio

I n Fig. 25, the red part is the parameter sensitive region, and the blue part is the insensitive region. As mentioned above, the abutment reinforcement ratio ϱ 1 ϱ 1 ϱ_(1)\varrho_{1}ϱ1 increases from 0.7 % 0.7 % 0.7%0.7 \%0.7% to 1.3 % 1.3 % 1.3%1.3 \%1.3%, and the bearing capacity increases by about 3 % 3 % 3%3 \%3% with a change of 0.2 percentage points. The influence coefficient Z a Z a Z_(a)Z_{a}Za of the abutment reinforcement ratio can be defined
as: 1 1.11 ( 0.7 % 1 1.11 ( 0.7 % 1-1.11(0.7%1-1.11(0.7 \%11.11(0.7% to 1.3 % ) 1.3 % ) 1.3%)1.3 \%)1.3%). The abutment reinforcement ratio ϱ 2 ϱ 2 ϱ_(2)\varrho_{2}ϱ2 increases linearly by 14 % 14 % 14%14 \%14% from 1 % 1 % 1%1 \%1% to 2 % 2 % 2%2 \%2%, and an impact factor Z p Z p Z_(p)Z_{p}Zp 1-1.13 ( 1 % 1 % 1%1 \%1% to 2 % 2 % 2%2 \%2% ) can be defined.
Figure 25: Reinforcement Sensitivity Distribution.

BEARING CAPACITY FORMULA

Through the fitting of 180 sets of numeric results plus theoretic analysis, the ultimate bearing capacity formula of the integral abutment joint is constructed as follows. Given the concrete slab reinforcement ratio, the capacity formula is defined as:
(25) F N = Z a Z p F s (25) F N = Z a Z p F s {:(25)F_(N)=Z_(a)Z_(p)*F_(s):}\begin{equation*} F_{N}=Z_{a} Z_{p} \cdot F_{s} \tag{25} \end{equation*}(25)FN=ZaZpFs
The previous section found that the ultimate bearing capacity of 800 mm girder is about 74 % 74 % 74%74 \%74% of that of 1000 mm girder. Therefore, a girder depth influence coefficient h , h 800 = 0.74 , h 1000 = 1 h , h 800 = 0.74 , h 1000 = 1 h,h_(800)=0.74,h_(1000)=1h, \mathrm{~h}_{800}=0.74, \mathrm{~h}_{1000}=1h, h800=0.74, h1000=1 can be obtained.
(26) F N = Z a × Z p × ( 421 + 0.43 × ( 1.106 l e 2 f c b w / 2.9 + F e ) 0.00014 × ( 1.106 l e 2 f c b w / 2.9 ) 2 ) (26) F N = Z a × Z p × 421 + 0.43 × 1.106 l e 2 f c b w / 2.9 + F e 0.00014 × 1.106 l e 2 f c b w / 2.9 2 {:(26)F_(N)=Z_(a)xxZ_(p)xx(421+0.43 xx(1.106l_(e)^(2)sqrt(f_(c)bw)//2.9+F_(e))-0.00014 xx(1.106l_(e)^(2)sqrt(f_(c)bw)//2.9)^(2)):}\begin{equation*} F_{N}=Z_{a} \times Z_{p} \times\left(421+0.43 \times\left(1.106 l_{e}^{2} \sqrt{f_{c} b w} / 2.9+F_{e}\right)-0.00014 \times\left(1.106 l_{e}^{2} \sqrt{f_{c} b w} / 2.9\right)^{2}\right) \tag{26} \end{equation*}(26)FN=Za×Zp×(421+0.43×(1.106le2fcbw/2.9+Fe)0.00014×(1.106le2fcbw/2.9)2)
In Fig. 26 F 1 26 F 1 26F_(1)26 F_{1}26F1 is the finite element result and F 2 F 2 F_(2)F_{2}F2 is the formula value, it is clear that the ultimate bearing capacity obtained by the formula is very close to the finite element results, and the average error is within 3 % 3 % 3%3 \%3%.
Fig. 26: Comparison between formula results and FE values
Based on the developed formula for the integral abutment, a comparison of bearing capacity and steel consumption between the integral abutment proposed in this work and the traditional I-shaped steel girder with the same size is carried out. In Liang's test [21], the width of the bottom plate of the steel girder is 220 mm , the depth of the steel girder is 540 mm , the embedding depth of the steel girder is 360 mm , the reinforcement rate of the abutment is 0.72 % 0.72 % 0.72%0.72 \%0.72%, the width of the abutment is 500 mm , and the ultimate bearing capacity of the test is 488 kN . The size of the steel girder in this paper was adjusted to match that in the test. Employing the calculation formula proposed in this paper, the ultimate bearing capacity is calculated as 507 kN . It is found that the bearing capacity at the girder-joint mainly relies on the longitudinal tensile bearing capacity provided by the bridge deck reinforcement, while the concrete at the contact surface between the steel girder and the concrete abutment transfers the external load through compression. Compared with the traditional I-shaped steel girder in the experiment, under the same steel girder depth and bottom plate width, the steel consumption of the integral abutment in this paper is reduced by 23 % 23 % 23%23 \%23%. This saving arises in part because the longitudinal reinforcement was not included in Liang's test [21]. Once the steel girder yielded, it would lose bearing capacity. In contrast, for the integral abutment in this work, the presence of longitudinal reinforcement in the bridge deck and the concrete web allows the steel girder to retain load-bearing capacity even after yielding.

Conclusion

\checkmark The failure modes and load transfer mechanisms of the proposed integral bridge abutment with composite dowels are found. Failure initiates with cracking in the concrete surrounding the composite dowels, leading to a loss of shear resistance. However, the penetrating steel bars between the dowels do not yield and continue to carry shear forces. After reaching the peak load, the abutment maintains its bearing capacity without a sharp decrease, though deformation progresses rapidly, exhibiting good ductility. Cracking in the abutment, steel girder deflection, and steel bar deformation increase to varying degrees. Ultimately, failure occurs due to the local yielding of the steel girder.
\checkmark Typical girder dimensions and concrete reinforcement of the integral abutment are suggested. With respect to the girder depth range from 800 mm 1000 mm 800 mm 1000 mm 800mm-1000mm800 \mathrm{~mm}-1000 \mathrm{~mm}800 mm1000 mm, an optimal dowel spacing 300 350 mm 300 350 mm 300-350mm300-350 \mathrm{~mm}300350 mm and web thickness of 20 mm are recommended. Additionally, increasing the reinforcement ratio of the abutment from 1 % 1 % 1%1 \%1% to 2 % 2 % 2%2 \%2% enhances the joint's bearing capacity by approximately 13 % 13 % 13%13 \%13%.
\checkmark A formula for the ultimate bearing capacity of the integral abutment joint is proposed, which accounts for contributions from both the composite girder and connections. Validation against numerical results shows an average error below 3 % 3 % 3%3 \%3%, confirming its accuracy. This formula facilitates rapid design optimization and serves as a reliable tool for engineers to assess the performance of the abutment under multifactorial loading conditions, including thermal effects and seismic actions.

ACKNOWLEDGEMENTS

T he authors appreciate the support of Natural Science Foundation of Shaanxi Province (Grant No. 2023-JC-YB-360)

REFERENCES

[1] Pak, D., Bigelow, H., & Feldmann, M. (2017). Design of composite bridges with integral abutments. Steel Construction, 10(1), pp. 23-30.
[2] Mertz, D. (2013). Design Guide for Bridges for Service Life. Strategic Highway.
[3] Paraschos, A. (2016). Effects of wingwall configurations on the behavior of integral abutment bridges (Doctoral dissertation, University of Maryland, College Park).
[4] Pnevmatikos, N. G.,Papagiannopoulos, G. A.,Papavasileiou, G. S. (2019). Fragility curves for mixed concrete/steel frames subjected to seismic excitation, Soil Dynamics and Earthquake Engineering 116, pp. 709-713.
[5] Arsoy, S., Duncan, J. M., & Barker, R. M. (2002). Performance of piles supporting integral bridges. Transportation Research Record, 1808(1), pp. 162-167.
[6] Dreyer, T., Fronek, A., Fataar, H., Mpye, G. D., Skorpen, S., Kearsley, E. and van Zijl, G. (2025). Soil-Structure Interaction Analysis of Transition Zone Deformations due to Imposed Abutment Displacement Cycles in an Integral Bridge. Journal of Bridge Engineering, 30(4), 04025005.
[7] Fiorentino, G., Cengiz, C., De Luca, F., Mylonakis, G., Karamitros, D., Dietz, M. and Nuti, C. (2021). Integral abutment bridges: investigation of seismic soil- structure interaction effects by shaking table testing. Earthquake Engineering & Structural Dynamics, 50(6), pp. 1517-1538.
[8] Kim, S. H., Yoon, J. H., Kim, J. H., Choi, W. J. and Ahn, J. H. (2012). Structural details of steel girder-abutment joints in integral bridges: An experimental study. Journal of Constructional Steel Research, 70, pp. 190-212.
[9] Feldmann, M., Kopp, M. and Pak, D. (2016). Composite dowels as shear connectors for composite girders-background to the German technical approval. Steel Construction, 9(2), pp. 80-88.
[10] Xiong, Z., Pan, Z., Wolters, K., Liu, H., Cheng, L., Li, J. and Feldmann, M. (2024). Prefabricated Composite Dowels Girders for Small Bridges-Structural Design and Efficiency, Evaluation.ce/papers, 7(3-4), pp. 50-62.
[11] Xiong, Z., Li, J., Wolters, K., Mou, X., & Feldmann, M. (2025). Fatigue performance of composite dowels under combination of shear and tension - Part I - Steel component. Steel Construction,18(1), pp. 52-62.
[12] Feldmann, M., Naumes, J., Pak, D., Veljkovic, M., Nilsson, M., Eriksen, J. and Popa, N. (2010). Economic and durable design of composite bridges with integral abutments, European Commission, DOI: 10.2777/91014.
[13] Ashiduka, K., Miyata, H., Sakate, M., Kiso, S., Kurita, A., Shimotaniike, T. and Leistung, P. (2007). Design of hybrid frame bridge having spread foundation and proposal for rationalization of connection detail between steel girder and RC abutment. Kou Kouzou Rombunshuu. A, 53, pp. 936-945.
[14] Riches, O. J., Carstairs, N. A., & Jones, A. E. (2005, June). A simplified integral composite bridge connection. In Proceedings of the Institution of Civil Engineers-Bridge Engineering, 158(2), pp. 63-69).
[15] VTrans Integral Abutment Committee. (2008). Integral Abutment Bridge Design Guidelines, the State of Vermont. Agency of Transportation, Montpelier.
[16] Briseghella, B. and Zordan, T. (2015). An innovative steel-concrete joint for integral abutment bridges. Journal of Traffic and Transportation Engineering (English Edition), 2(4), pp. 209-222.
[17] Far, N. E., Maleki, S. and Barghian, M. (2015). Design of integral abutment bridges for combined thermal and seismic loads. Earthquakes and Structures, 9(2), pp. 415-430.
[18] LaFave, J. M., Riddle, J. K., Jarrett, M. W., Wright, B. A., Svatora, J. S., An, H. and Fahnestock, L. A. (2016). Numerical simulations of steel integral abutment bridges under thermal loading. Journal of Bridge Engineering, 21(10), 04016061.
[19] Kunin, J. and Alampalli, S. (2000). Integral abutment bridges: Current practice in United States and Canada. Journal of performance of constructed facilities, 14(3), pp. 104-111.
[20] Aktan, H., Attanayake, U. and Ulku, E. (2008). Combining link slab, deck sliding over backwall, and revising bearings (No. RC-1514).
[21] Liang, C., Liu, Y., Zhao, C., Lei, B. and Wu, J. (2018). Experimental and numerical study on an innovative girderabutment joint in composite bridges with integral abutments. Construction and Building Materials, 186, pp. 709-730.
[22] Veljkovic, M. (2007). Project: Economic and Durable Design of Composite Bridges with Integral Abutments.
[23] Mattock, A. H. and Gaafar, G. H. (1982, March). Strength of embedded steel sections as brackets. Journal Proceedings 79(2), pp. 83-93.
[24] Tschemmernegg, F. and Humer, C. (1988). The design of structural steel frames under consideration of the nonlinear behaviour of joints. Journal of Constructional Steel Research, 11(2), pp. 73-103.
[25] Hawkins, N. M. (1970). The bearing strength of concrete for strip loadings. Magazine of concrete research, 22(71), pp. 87-98.