Zhou Int J Concr Struct Mater (2021) 15:51

https://doi.org/10.1186/s40069-021-00490-9 International Journal of Concrete

Structures and Materials

RESEARCH Open Access

Investigation of Bursting Stress and Spalling

Stress in Post‑tensioned Anchorage Zones
Linyun Zhou*

Abstract
Over the past decades, considerable efforts have been made to quantify the bursting forces in the post-tensioned
anchorage zones based on the simplified model or fitting formulas, however reproducing the transverse stress
distribution is still a challenging topic, which is also important to detail the reinforcing details in the anchorage zones,
especially for cracking control. To address this issue, this paper is devoted to seeking an elasticity solution for trans-
verse stresses in the anchorage zones, and providing a more rational equation for transverse distribution in anchorage
zones. The sum function of normal stresses is employed to solve the stresses filed in the anchorage zones with con-
centric load and two eccentric loads. The bursting stresses in the concentric anchorage zones and spalling stresses
in the eccentric anchorage zones are verified by the photoelastic tests. The transverse stresses along the symmetry
axis of the eccentric anchorage zones can be handled as a concentric single anchorage zone with equivalent bear-
ing plate width. Moreover, according to the concept of force stream tube, the profiles of isostatic line of compression
(ILCs) are determined and validated, which confirms the existence of ILCs.
Keywords: prestressed concrete, anchorage zone, bursting stress, spalling stress, photoelastic test, isostatic line of
compression

1 Introduction finite element analysis, experimental studies and the

The anchorage zone is defined as the portion of struc- strut-and-tie models.
ture in which the concentrated post-tensioning forces are Based on the theory of elasticity, Guyon (1953) studied
transferred from the anchorage device to the structure. In the behavior of a concentric load applied over a rectangu-
post-tensioned anchorage zones, the spread of the con- lar body, and then provided a formulation to determine
centrated forces will produce transverse tensile stresses, the bursting stress distribution ahead of a concentric end
usually referred to as bursting stresses, along the tendon anchor for different bearing plate ratios. This relation was
path. The resultant of the bursting stresses is called burst- very influential in the design of anchorage zones. Yet-
ing forces. Practically, both the bursting forces and the tram and Robbins (1970, 1971) reported that the location
bursting stresses are used to detail the reinforcing details of uniform stress and the effect of Poisson’s ratio were
in the anchorage zone. insignificant to the stress distribution based on elastic
To date, a large number of studies of the prediction of finite element analysis. Moreover, the extension of burst-
bursting stresses and bursting forces in anchorage zone ing stresses of I-sections was further than in rectangu-
have been conducted, including linear elastic studies, lar sections and the flanges played an important role in
reducing the spalling stresses. Sanders and Breen (1997)
conducted experimental studies of anchorage zones
with 36 specimens tested. By assuming that the result-
*Correspondence: [email protected] ing force at the end of the specimen would be shifted
School of Transportation, Southeast University, Sipailou 2, Xuanwu in the direction of the crack, the dispersion angle of the
district, Nanjing 210096, People’s Republic of China compressive strut should be reduced, then a modified
Journal information: ISSN 1976-0485/eISSN 2234-1315

© The Author(s) 2021. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which
permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the
original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or
other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line
to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory
regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this
licence, visit http://​creat​iveco​mmons.​org/​licen​ses/​by/4.​0/.
Zhou Int J Concr Struct Mater (2021) 15:51 Page 2 of 13

strut-and-tie model was proposed to predict the ulti- Over the past decades, considerable efforts have been made
mate load of the anchorage zones. Wollmann (1991) con- to quantify the bursting forces in the post-tensioned anchor-
ducted three experimental tests of concentric anchorage age zones based on the simplified model or fitting formulas,
zones to study the influence of reaction forces on burst- however reproducing the transverse stress distribution is still
ing stresses. For small tendon eccentricity, the influence a challenging topic, which is also important to detail the rein-
of reaction forces on linear-elastic bursting forces and forcing details in the anchorage zone. To address this issue, this
the resulting bursting force is conservative and can be paper is devoted to seeking an analytical solution for transverse
neglected. Besides, the reaction forces have little influ- stresses in the anchorage zones, and providing a more rational
ence on the failure mode of anchorage zones. Based on equation for transverse distribution in anchorage zones.
the results of those tests and finite element analyses, a
fitting formula was proposed for calculating the burst- 2 Stress Solution of Anchorage Zones
ing forces for typical anchorage zones (Breen et al., 1994), 2.1 Sum Function Method
which had been adopted in the AASHTO-LRFD Bridge According to the theory of elasticity (Timoshenko &
Design Specifications (2014) since 1994. Foster and Goodier, 1951), the equilibrium equation of the plane
Rogowsky (1997) studied the bursting stress distributions stress problem can be presented by
in the case of service load behavior, and distributions in
∂σx ∂τxy
the case of ultimate load behavior of anchorage zones by + = 0, (1a)
using rotating crack finite element model. They stated ∂x ∂y
that the bursting stress distribution was flatter than that
of the linear analysis due to the stress redistribution after ∂σy ∂τxy
cracking, which led to a conservative estimation of the
+ = 0, (1b)
∂y ∂x
bursting forces, and then a modified equation was pro-
posed to estimate the load capacity of the anchorage and then, the equilibrium differential equation can be
zones. After analyzing the stress magnitudes and distri- given by
butions of rectangular post-tensioned anchorage zones
∂ 2I
using ultra-high performance concrete (UHPC), Kim and ∇ 2 τxy = − , (2)
Kim (2017) concluded that the use of UHPC gives signifi- ∂x∂y
cant reduction of anchorage zone size and no reinforce- 2 2
in which ∇ 2 = ∂x∂
2 + ∂y2 ; I = σx + σy is the sum function
∂
ments are required. Yun (2005) proposed the nonlinear
strut-and-tie model approach to predict the behavior of normal stresses.
of structural concrete, and then extended this study to The stress solution should satisfy the harmonious equa-
the ultimate strengths of post-tensioned anchorages by tion, that is
checking the occurrence of a nodal zone failure mecha- ∇ 2 I = 0. (3)
nism, geometric compatibility condition, and the struc-
tural instability of the model struts and ties. Following As can be seen, the sum function, I , can be figured out
the concept of isostatic line of compression (ILCs) ini- from harmonious equation; and then the shear stress
tialed by Guyon (1953), Sahoo et al. (2009) presented an can be obtained by substituting I into Eq. (2); finally, the
analytical equation for estimating the bursting stress and stress components will be solved based on the equilib-
bursting forces by introducing certain boundary condi- rium equations and boundary conditions.
tions of ILCs. However, this model was imperfect due to
some questionable assumptions on the boundary condi- 2.2 Solution of Stress Components
tions (Windisch, 2010). He and Liu (2011) proposed the For the problem as illustrated in Fig. 1, prescribed forces,
compression dispersion model (CDM) by excluding the q(y), are symmetrically loaded on the plate. According to
unreasonable assumptions in Sahoo’s model. Neverthe- symmetries, the sum function, I , can be assumed as
less, the CDM failed to reproduce the bursting stresses 
along the tendon path due to insufficient boundary con- I= f (x) cos βy, (4)
ditions. Based on the equilibrium condition at the far-end
in which β = nπ/h, f (x) is the undetermined functions,
of the anchorage zones, Zhou et al. (2015) updated the
and h is the height of the plate.
CDM by incorporating two more boundary conditions,
Substituting Eq. (4) into Eq. (3) yields
and verified by finite element analysis. However, in this
updated model, the location of interface section of com-

[f ′′ (x) − β 2 f (x)] cos βy = 0. (5)
pressive and tensile stress was assumed to be constant,
what is in contradiction with reality. Then, the sum function, I , can be solved by
Zhou Int J Concr Struct Mater (2021) 15:51 Page 3 of 13

y y

q(y) h x q(y)
a h x
P P

l
Fig. 1 Plate with symmetrically prescribed forces. l=2h
Fig. 2 Concentric anchorage zone.


I= (an chβx + bn shβx) cos βy, (6)  an
σy = (cn shβx + dn chβx + chβx
in which coefficients an and bn can be determined by 2
boundary conditions. an bn
+ βxshβx + βxchβx
Substituting sum function into Eq. (2), the shear stress 2 2
can be derived by bn
an + shβx) cos βy + D2 (x),
(9b)

τxy = (cn chβx + dn shβx + βxchβx 2
2 (7)
bn an bn
+ βxshβx) sin βy , τxy =

(cn chβx + dn shβx + βxchβx + βxshβx) sin βy ,
2 2 2 (9c)
where cn and dn are constraints determined by boundary
in which function D1 (y) and D2 (x) are determined by the
conditions.
stress boundary conditions.
According to the relationship of normal stress and
shear stress illustrated in equilibrium equation, we have
∂σx ∂τxy 2.3 Transverse Stress in Concentric Single Anchorage Zone
=−
∂x ∂y According to Huang (2012), the anchorage shape and
 an duct hole have little effect on bursting stresses and
=− β(cn chβx + dn shβx + βxchβx bursting forces, thus the concentric anchorage zone
2
bn can be modeled as a plane stress problem as depicted
+ βxshβx) cos βy , in Fig. 2. The concentric anchor force, P , loaded on the
2 (8a)
bearing plate can be expressed by the Fourier series,
∂σy ∂τxy an bn 
(10)

=− =− β(cn shβx + dn chβx + chβx + βshβx q(y) = A0 − An cos βy,
∂y ∂x 2 2
an bn
+ βxshβx + βxchβx) sin βy. in which
2 2
(8b) aq
A0 = , (11a)
Finally, the stress components can be derived by inte- h
grating the above equations:
2q
 an An = sin aβ, (11b)
σx = − (cn shβx + dn chβx + βxshβx hβ
2
an bn where q = P/(at); a is the bearing plate width; and t is
− chβx + βxchβx
2 2 the thickness of the anchorage zone.
bn From the above equation, the anchor force can be
− shβx) cos βy + D1 (y),
2 (9a) treated as two parts, one of which is a uniform load,
Zhou Int J Concr Struct Mater (2021) 15:51 Page 4 of 13

0, and the other a polynomial distributed load,

q11 (y) = A aq  an an
σx = − − (dn chβx + βxshβx − chβx) cos βy,
q12 (y) = An cos βy. h 2 2
For the case of uniform load, the stress components in (16a)
the anchorage zone can be given by
 an an
σy = (dn chβx + chβx + βxshβx) cos βy,
2 2
σx = −q11 = −aq/h, (12a) (16b)
 an
τxy = (dn shβx + βxchβx) sin βy . (16c)
σy = 0, (12b) 2
According to Eq. (16b), the transverse stress along the
τxy = 0. (12c) tendon path, i.e., y = 0, can be obtained by
While for the case of distributed load, the stress bound-  an an
σy = (dn chβx + chβx + βxshβx)
ary conditions can be presented by 2 2
 an (17)
 ≈ (1 + βx − βl)chβx.
τxy x=±l = 0, (13a) 2
 Representing the above equation by the fourth-order
σx |x=±l = An cos βy, (13b) Taylor series expansion, that is,
 An πx π 2 x2 π 4 x4
 l σy = ( − 5.28)(1 + + ).
2 24h4
σy dx = 0, (13c) sh2nπ h 2h
−l (18)
Cartesian coordinate system is established in the anchor-
h age zone as shown in Fig. 3, the direction parallel to the

σx dy = P. (13d) load is set to be the x-axis. Therefore, the bursting stress in
−h
the concentric anchorage zone can be approximated by the
Imposing the stress boundary equations into Eq. (9), we first two terms:
have  
12P sin(π a/2h)(h − a) h sin(π a/2h)
2shβl 2An σTc = − (x − )
an = An ≈ , (14a) h6 t 4
βl + shβlchβl shβl (x − h)3 .
(19)
bn = cn = 0, (14b)
It can be seen that the interface section
√ of compressive
h sin(π a/2h)
βlchβl An βl and tensile stress is located at x0 = 6 , where the
dn = −An ≈− , (14c) transverse stress is equal to zero. And the transverse stress
βl + shβlchβl shβl
reaches its maximum value σT ,max at x = (h + 3x0 )/4
D1 (y) = D2 (x) = 0. (14d) 0.57P(h − a)
σTc,max ≈ . (20)
Therefore, the stress induced by the distributed load can h2 t
be presented by Thus, the bursting force along the tendon axis Tb can be
 an an written as
σx = − (dn chβx + βxshβx − chβx) cos βy,
2 2  h
a
(15a) Tb = σTc tdx ≈ 0.24P(1 − ). (21)
 an an x0 h
σy = (dn chβx + chβx + βxshβx) cos βy,
2 2 It can be seen that the bursting forces in concentric
(15b)
anchorage zones agree well with the formulas given by cur-
an
rent design codes (ACI 318-08 2008, AASHTO, 2014, CEB-

τxy = (dn shβx + βxchβx) sin βy . (15c)
2 FIP 1990)
Finally, the stress components in the concentric anchor- a
age zone can be derived by Tb = 0.25P(1 − ). (22)
h
Zhou Int J Concr Struct Mater (2021) 15:51 Page 5 of 13

P y P
a a
e h x e

a a
P P

l=2h
Fig. 4 Anchorage zone with two anchor forces.

2shβl 2Bn
an = Bn ≈ , (25a)
βl + shβlchβl shβl

Fig.3 Concentric anchorage zone with single load. βlchβl Bn βl

dn = −Bn ≈− . (25b)
βl + shβlchβl shβl
2.4 Transverse Stress in Anchorage Zones with Two
Anchors Then, the transverse stress, σy , in the Cartesian coordi-
As for the problem illustrated in Fig. 4, the two anchor nate system as shown in Fig. 5, can be approximated by
forces, having their centers at a distance e in the x-direc-
tion, can be expressed by 4q eπ aπ π π π
σy ≈ cos sin (1 + x)ch( x − 2π ) cos y.
π sh2π h h h h h
(26)

q(y) = B0 + Bn cos βy, (23a)

As can be seen, the magnitude of transverse stress is

in which mainly dependent on the eccentricity e. When the anchor
devices are placed at the quarter points, that is e = h/2,
2aq
B0 = , (23b) no transverse stress will be produced in the anchorage
h zone. While for the case of the anchor forces acting over
half of the depth, that is, e > h/2, strain compatibility will
4q produce spalling stress between the anchors, and the dis-
Bn = cos eβ sin aβ. (23c)
βh persion of concentric loads will generate bursting stress.
Similarly, the spalling stress along the x-axis, σTs, and the
Similarly, the anchor forces can be divided into two
transverse bursting stress along the tendon path, σTb, can
parts, that is, a uniform load, q
21 (y) = B0, and a polyno-
be approximated by
mial distributed load, q22 (y) = Bn cos βy.
Following the similar method as mentioned before, the

26P sin πas /2h(h − as )
stress solution of the anchorage zone with two eccentric σTs =
h6 t
anchor loads can be given by   
h sin πas /2h
2aq  an an x− (x − h)3 ,
6
σx = −
h
− (dn chβx + βxshβx − chβx) cos βy,
2 2 (27a)
(24a)

 an an sin πa/2h′ (h′ − a)
6P
σy = (dn chβx + chβx + βxshβx) cos βy, σTb =−
2 2 h′6 t
(24b) 
′
 
h sin πa/2h′
 an x− (x − h′ )3 ,
τxy = (dn shβx + βxchβx) sin βy , (24c) 6
(27b)
2
in which in which as = 2e − h; and h′ = h − e.
Zhou Int J Concr Struct Mater (2021) 15:51 Page 6 of 13

the quarter points, the spalling stress will be produced

between the anchors, and its distribution can be approxi-
mated by the negative value of the bursting stress in the
concentric anchorage zone with bearing plate width
as = 2e − h. Thus, a unified equation can be presented
to express the transverse stress along the centerline of
anchorage zone

12P sin πae /2h(h − ae )
σTb = ∓
h6 t

h sin πae /2h
[x − ](x − h)3 ,
6 (31)
in which ae is the equivalent bearing plate width and can
be defined by

e e < h/2
ae = . (32)
Fig. 5 Anchorage zone with two anchor forces. e − h/2 e > h/2

3 Validation
Equation (27) shows that the bursting stress along the 3.1 Test Program
tendon path in the anchorage zone with two eccentric In order to verify the accuracy of the proposed equa-
anchor loads is equal to the bursting stress in the concen- tions, 10 photoelastic tests (Timoshenko & Goodier,
tric anchorage zone with the depth h′ = h − e, which is 1951) have been carried out. All the test specimens are
the principle of symmetrical prism proposed by Guyon made of polycarbonate, and have the rectangular cross
(1953). section b × t = 80 × 8 mm. Young’s modulus of poly-
The corresponding spalling forces, Ts, and bursting carbonate is taken to be E = 3200MPa, Poisson’s ratio
forces, Tb, can be given by as ν = 0.37, and photoelastic material fringe value as
P(h − as ) f = 308N /order/m. Details of the test specimens are
Ts ≈ 0.12 , (28a) listed in Table 1. All of the test specimens are monotoni-
h
cally loaded over the full thickness so that the specimen
a would behave as two-dimensional as possible. Specimens
Tb ≈ 0.12P(1 − ). (28b) C1–C5 were tested to investigate the effect of the bearing
h−e
plate ratio to the transverse bursting stress in concentric
While for the case of e < h/2, the transverse stress anchorage zones as shown in Fig. 6; while M1–M5 were
along the x-axis, σTb, can be rearranged by tested to study the influence of eccentricities on trans-
 verse stress.
12P sin π ab /2h(h − ab ) Fig. 7 illustrates the experimental setup to obtain
σT = −
 h6 t  photoelasticity images, which consists of light source,
polarizer, quarter wave plate, analyzer and camera. Two

h sin π ab /2h
x− (x − h)3 , setups were used, one for employing white light and the
6 (29)
other for monochromatic light, the use of which will be
where ab = 2e. explained later.
Then, the bursting forces, Tb, can be obtained by
3.2 Test Procedure
e
Tb ≈ 0.24P(1 − ). (30) During the test, the light will be linearly polarized by
h the polarizer and circularly polarized by the first quar-
It can be seen that when the multi-anchors are placed ter wave plate. After passing through the test specimen,
inside the quarter points, that is, the anchor spacing is the light is linearly polarized by the second quarter
less than 0.5h, the multi anchorage zone can be handled wave plate and passes through the analyzer, which cre-
as a concentric anchorage zone with the bearing plate ates dark or bright regions, i.e., the fringe pattern. The
width ab = 2e; while the two anchors are placed outside
Zhou Int J Concr Struct Mater (2021) 15:51 Page 7 of 13

Table 1 Details of the test specimens. allel to the ground and that of the analyzer is vertical
Specimen height Bearing Anchor Applied load
to the ground. Capture the isoclines patterns.
h (mm) plate ratio spacing P (N) 3. Rotate polarizer and analyzer by 5 degrees simulta-
a/h e/h neously. Capture the isoclines patterns.
C1 80 0.1 0 33
4. Repeat step 3 until the a complete rotation of 90
C2 80 0.2 0 58
degrees. Capture the isoclines patterns (see Figs. 8b
C3 80 0.3 0 60
and 9b).
C4 80 0.4 0 60
C5 80 0.5 0 66
M1 80 0.1 0.3 32
M2 80 0.1 0.4 38
3.3 Stress Solution
M3 80 0.1 0.5 40
According to the photoelasticity images captured by
M4 80 0.1 0.6 40
camera, the isochromatic and isocline patterns can be
M5 80 0.1 0.8 40
sketched as shown in Fig. 10a. Following the shear differ-
ence method (Frocht 1941), the stress components in the
desired integration path MN (see Fig. 10b) can be given
by
Polycarbonate specimen i �x
(σx )i = (σx )i−1 − (�τxy )i−1 , (33a)
�y

(σy )i = (σx )i ± |m cos 2α|, (33b)

h = 0 .08
m
τxy = sin 2α, (33c)
2

l = 0.16 in which
Fig. 6 Test specimen (unit: m).

Quarter wave Quarter wave

Polarizer plate plate Analyzer

Light source Camera

Specimen

Fig. 7 Setup for analysis by photoelasticity.

fringe pattern is then captured by the camera. The brief

test can be illustrated as following:

1. Apply load gradually to the test specimen in the cir-

cular polariscope system obtained from the white
light source. Capture the isochromatic patterns (see
Figs. 8a and 9a). Fig. 8 Photoelasticity images for the concentric anchorage zones
2. In the dark-field orthogonal linear polarization sys- with a/h = 0.1; a isochromatic pattern; and b isocline patterns at 30
tem obtained from the monochromatic light, make degrees.
sure that the polarization axis of the polarizer is par-
Zhou Int J Concr Struct Mater (2021) 15:51 Page 8 of 13

3.4 Concentric Anchorage Zones

Fig. 12 shows the comparisons of the bursting stresses
in the concentric anchorage zones calculated by differ-
ent methods. It can be seen that the proposed equation
can predict the transverse tensile stress accurately, and
the maximum relative error is only 3.2%. It also should
be noted that the interface section of compressive and
tensile stress, where transverse stress σb = 0, varies with
the bearing plate ratio. In recent previous analytical stud-
ies of the bursting stresses (He et al., 2011; Zhou et al.,
2015), the interface section was assumed to be constant.
Actually, as early as 1950s, Guyon (1953) presented a
fitting formula to capture the location of this section as
shown in Fig. 13. It is clear that both Eq. (19) and Guyon’s
formula (1953) produce a good agreement with the test
results, and the predicted values of Eq. (19) are on aver-
age 96.4% of the test values, with a standard deviation of
Fig. 9 Photoelasticity images for the multi concentric anchorage 2.7%; these two numbers are 104.8% and 3.7%, respec-
zones with anchor spacing e/h = 0.6; a isochromatic pattern; and b
tively, for Guyon equation (1953).
isocline patterns at 10 degrees.
Fig. 14 illustrates the comparison of bursting force
with different bearing plate ratios. As can be seen, both
Eq. (21) and AASHTO (2014) equation agree well with
�τxy = (τxy )AB − (τxy )CD , (34a) the test results. The average relative deviations are only
5.7 and 3.8%, respectively.
i (�τxy )i−1 + (�τxy )i
(�τxy )i−1 = , (34b) 3.5 Two Anchorage Zones
2 The transverse stress along the centerline of specimen for
where m is the order isochromatic, α is the inclination of the anchorage zone with two eccentric anchor loads with
the isocline with respect to x-axis. different anchor spacing is illustrated in Fig. 15. It can be
Table 2 lists the stress solution of the section MN, at seen that the proposed Eq. (30) give very close results to
a depth of 0.5 h below the loading plate, in the concen- the test transverse stress. The average relative deviation is
tric anchorage zone with a/h = 0.1. Fig. 11a, b show the only 4.8%. Fig. 16 shows the influence of the anchor spac-
test normal stress components in anchorage zone with ing on the transverse tensile forces in anchorage zones
single concentric load and two eccentric anchor loads. (the size of anchor plate is kept constant at a = 0.1 h). It

5 80°
y
0
40°
1
2
4
30° A B
3 20°
∆y N x
M N M 0 1 2 3 19 20
10°
C D
∆x
0°

(a) (b)
Fig. 10 Photoelasticity images for concentric anchorage zones with a/h = 0.1: a isochromatic and isocline sketches; and b path of integration.
Zhou Int J Concr Struct Mater (2021) 15:51 Page 9 of 13

P=33N P=20N P=20N

-225.6 59.3
24.7

-120.3 x x
16.9

-89.4
23.1 -11.1

-74.6 -11.1
19.3
x y -46.6
14.0 -66.9
-7.1

x 9.4 x y y -21.3
y x

(a) (b)
Fig. 11 Test stress isoline of the anchorage zones (unit kPa): a concentric anchorage zone with a/h = 0.1; and b multi anchorage zone with
e/h = 0.8.

Table 2 Stress components calculated by shear difference method.

Node Section AB Section CD Section MN
m α τAB(kPa/order) m α(°) τAB(kPa/order) σx σy τxy
(order) (°) (order) (kPa) (kPa) (kPa)

0 0.26 3.8 0.02 0.76 1.88 0.02 0.0 − 24.1 0.0

1 0.58 10.9 0.11 0.88 8.0 0.12 0.3 − 29.1 − 0.2
2 0.92 16.5 0.25 1.12 9.9 0.19 1.2 − 35.3 − 0.6
3 1.38 16.7 0.38 1.75 12.9 0.38 2.3 − 46.4 − 1.3
4 1.61 16.7 0.44 1.83 14.2 0.43 4.2 − 49.4 − 2.1
5 1.78 15.3 0.45 2.01 12.2 0.41 6.9 − 58.4 − 2.7
6 1.98 13.7 0.45 2.24 9.8 0.37 9.2 − 66.1 − 3.2
7 2.28 10.2 0.40 2.49 7.4 0.32 12.3 − 74.5 − 3.3
8 2.51 7.2 0.31 2.66 4.6 0.21 15.8 − 81.0 − 2.7
9 2.69 4.1 0.19 2.75 3.2 0.15 18.4 − 87.1 − 1.6
10 2.76 0.0 0.00 2.83 0.0 0.00 19.0 − 87.8 0.0
Note: the values are for “at the location of 0.5 h below the loading plate in the concentric anchorage zone with a/h = 0.1”

Fig. 12 Transverse stress in the concentric anchorage zones with

different bearing plate ratio. Fig. 13 The interface section of compressive and tensile stress.
Zhou Int J Concr Struct Mater (2021) 15:51 Page 10 of 13

Fig. 16 Transverse tensile forces in the multi anchorage zones.

Fig. 14 Bursting forces in concentric anchorage zones.

(1953), there must exist an ILC along the centerline of the

indicates that the proposed equation produces a good anchorage zone. Referring to Fig. 18a, the force stream
agreement with the test results, the average relative devi- tube (Kelly & Elsley, 1995) bounded by two ILCs is con-
ation is only 4.0%. sidered. The force over flow tube should satisfy

4 Towards the Existence of ILCs

Fj = Fj−1 = P�a/a, (35)
As early as 1953, Guyon (1953) used the ILCs to visualize in which P is the anchor force; and a is width of anchor-
the flow of forces in post-tensioned anchorage zones as age zone; a is the vertical distance between the center-
shown in Fig. 17. Each individual ILC transfers an equal line line and ILC; Fj is the force in the stream tube in the
share of the force, which induces a certain amount of arbitrary section and can be defined by
transverse stress depending on its curvature. However,
this is just a conceptual model, and the existence of ILCs (σx,Aj + σx,Bj )t
Fj = yj , (36)
is presently unproven. 2 cos
θj
2
Although the longitudinal compressive stress in
anchorage zones will be effected by the duct hole sig- in which σy,Aj and σy,Bj are the x-direction normal stresses
nificantly (Park et al., 2020), it has little effect on the of the centerline and ILC at x = Aj , respectively; t is the
distribution of transverse bursting stress (Huang, 2012). thickness of anchorage zone; yj is the vertical distance
Thus, the anchorage zone can be idealized as a concrete between the centerline line and ILC at x = Aj ; and θj is
block without hole. According to the concept of Guyon the local angular orientation of the force tube wall.

Fig.15 Transverse stress in multi anchorage zone.

Zhou Int J Concr Struct Mater (2021) 15:51 Page 11 of 13

According to the geometrical condition as shown in

P/6
P/6 Fig. 19, the slope of the ILCs at the inflection point can be
approximated by
P/6
P/6

dy 
P P/6
P/6 = tan θ ≈ 2m/(1 − m2 ), (40)
h a dx x=x0
P/6
P/6
in which m = (�h − �a)/l ; h and a are the verti-
P/6
P/6 cal value of ILCs at starting point and ending point,
ILCs
P/6
P/6 respectively.
Table 3 lists the bursting forces in the concentric
Fig. 17 ILCs in anchorage zones.
anchorage zone predicted by the proposed equation and
ILCs. As can be seen, the predictions of ILCs agree well
with the test results. The average discrepancy between
Substituting Eq. (35) into Eq. (36) leads to the ILCs predictions and test results is only 13%, which is
mainly induced by iterative process. Therefore, the ILCs
2P �a θj
yj = cos . (37) found by proposed method can reveal the load transfer
(σx,Aj + σx,Bj )t a 2 mechanism in anchorage zones. That is, ILCs are real
existence in the anchorage zones and transmit the anchor
According to the geometrical condition of the ILC as
loads to the far-end section.
illustrated in Fig. 18b, we have
yj = yj−1 + �x tan θj−1 . (38) 5 Conclusions
To clear the stress field in the anchorage zones, an elas-
As can be seen from the above equation, if y1 = a ticity solution of the transverse stress is derived based
and x is specified, the vertical location of the ILCs on the sum function method. The sum function of nor-
will be determined by an iterative solution. Fig. 19 illus- mal stresses can be figured out from harmonious equa-
trates the shapes of ILCs in concentric anchorage zone tion, and then the stress field in the anchorage zone can
with different bearing plate ratios, in which only a sin- be determined by equilibrium equations and boundary
gle ILC is plotted for the sake of simplification. conditions. According to the photoelastic test results and
Thus, the bursting forces can be derived by their comparison with the proposed equations, the fol-
lowing conclusions can be drawn:
P d2y

dy 
 
T = σtran tdx = − tdx = −P .
t dx2 dx x=x0
(1) The proposed equations can reproduce the trans-
(39) verse stress in concentric anchorage zones, and is
verified to have good accuracy by comparisons with
test results. On the other hand, the other existing

(a) (b)

Fig. 18 Determination the profile of ILC in concentric anchorage zone: a force stream tube; and b location of ILC.
Zhou Int J Concr Struct Mater (2021) 15:51 Page 12 of 13

Abbreviations
a: Width of the bearing plate along the direction of the beam height; h
: Height of the beam; t : Thickness of the anchorage zone; P : Anchor force; e
: Eccentricity of the anchoring point to the centroid anchorage zone; σx , σy ,
τxy: Normal stresses and shear stress; I = σx + σy: Sum function of normal
stresses; σTc: Bursting stresses distribution of concentric anchorage zone;
σT ,max: Maximum bursting stresses in concentric anchorage zone; Tb: Burst-
ing force in concentric anchorage zone; σTb: Bursting stresses distribution of
eccentric anchorage zone; σTs: Spalling stresses of eccentric anchorage zone;
Tb: Bursting force in eccentric anchorage zone; Ts: Spalling force in eccentric
anchorage zone; σT : Transverse stress along the centerline of anchorage
zone; ae: Equivalent bearing plate width; Fj: Force in stream tube at arbitrary
section; yj: Vertical distance between the centerline line and ILC; θj: Local
angular orientation of the force tube wall; h, a: Vertical value of ILCs at
starting point and ending point, respectively.

Acknowledgements
The research team for new long-span bridges in Southeast University, Nanjing,
China, is gratefully acknowledged for the test program in this paper.

Authors’ contributions
Fig. 19 ILCs in the concentric anchorage zones with different
The author has contributed substantially to all aspects of this article. The
bearing plate ratio. author read and approved the final manuscript, and has the appropriate
permission and rights to the reported data. All authors read and approved the
final manuscript.

analytical models cannot predict the transverse Authors’ information

stress distribution accurately, due to the incorrect Linyun Zhou is a post-doc in Southeast University, School of Transportation.
His research interest included strut-and-tie model, bridge engineering, and
location of the interface section of compressive and UHPC.
tensile stress.
(2) The transverse stress along the symmetry axis of Funding
This study was supported by the National Natural Science Foundation of
the anchorage zone with two eccentric anchor loads China (Grant No. 51808208) and the fundamental research funds for the
can be handled as a concentric single anchorage central universities (Grant No. 2242021R20011).
zone with equivalent bearing plate width.
Availability of data and materials
(3) It is very interesting that the bursting stress along All data, models, or code that support the findings of this study are available
the tendon path in the anchorage zone with two from the corresponding author upon reasonable request. All data, models,
anchors is equal to the bursting stress in the con- and code generated or used during the study appear in the submitted article.
centric anchorage zone with the depth h′ = h − e,
which is the principle of symmetrical prism. Declarations
(4) The major significance of this investigation not only Competing interests
lies in reproducing the transverse stress in anchor- The authors declare that there have no competing interests in this paper.
age zones, but also in confirming the existence of
Received: 13 March 2021 Accepted: 9 December 2021
the ILCs.

References
AASHTO. (2014). AASHTO LRFD bridge design specifications. American Associa-
Table 3 Bursting forces calculated by different methods. tion of State Highway and Transportation Officials.
ACI committee 318, (2008). Building Code Requirements for Structural Concrete
a/h h (mm ) m Bursting forces T (ACI 318–08) and Commentary (ACI 318R-08). Farmington Hills (MI), Ameri-
can Concrete Institute.
Test Equation (21)/ ILCs/Test Breen, J. E., Burdet, O., Roberts, C., Sanders, D., & Wollmann G. (1994). Anchor-
Test age zone reinforcement for post-tensioned concrete girders, National
Cooperative Highway Research Program Report 356, Washington, D. C.,
0.1 0.08 0.104 0.23P 0.93 0.81 p. 8–11.
0.2 0.08 0.105 0.20P 0.96 0.84 Comité Euro-International du Béton, Fédération International de la Précon-
trainte (CEB-FIP), (1990). Practical design of structural concrete. London, UK.
0.3 0.08 0.106 0.17P 0.97 0.86
Foster, J., & Rogowsky, M. (1997). Bursting forces in concrete panels resulting
0.4 0.08 0.107 0.15P 0.98 0.89 from in-plane concentrated loads. Magazine of Concrete Research, 180,
0.5 0.08 0.109 0.12P 1.02 0.93 231–240.
Frocht, M.M. (1941). Photoelasticity. John Wiley & Sons, Inc.
Mean 0.97 0.87
Guyon, Y. (1953). Prestressed concrete. John Wiley and Sons Inc.
Standard deviation 0.03 0.05
Zhou Int J Concr Struct Mater (2021) 15:51 Page 13 of 13

He, Z., & Liu, Z. (2011). Investigation of bursting forces in anchorage zones:
Compression-dispersion models and unified design equation. ASCE
Journal of Bridge Engineering, 16(6), 820–827.
Huang, H. Q. (2012). Reinforcement design for anchorage zone in post tensioned
end block based on strut and tie model. Master thesis, Southeast University.
Kelly, D., & Elsley, A. (1995). Procedure for determining load paths in elastic
continua. Engineering Computations, 5, 415–424.
Kim, J. S., & Kim, T. H. (2017). A stress analysis of the post-tensioned anchorage
zones using UHPC. Key Engineering Materials, 737, 500–504.
Park, Y. H., Kim, M. Y., Park, J. M., & Jeon, S. J. (2020). An improved equation
for predicting compressive stress in posttensioned anchorage zone.
Advances in Civil Engineering, 2020, 5635060.
Sahoo, D. K., Singh, B., & Bhargava, P. (2009). Investigation of dispersion of com-
pression in bottle-shaped struts. ACI Structural Journal, 106, 178–186.
Sanders, D., & Breen, J. E. (1997). Post-tensioned anchorage zones with single
straight concentric anchorages. ACI Structural Journal, 2, 146–158.
Timoshenko, S., & Goodier, J. (1951). Theory of elasticity. Mcgraw-Hill book
company.
Windisch, A. (2010). Investigation of dispersion of compression in bottle-
shaped struts-discussion. ACI Structural Journal, 107, 124–126.
Wollmann, G. (1991). Design of anchorage zones in post-tensioned concrete
structures. PhD dissertation, University of Texas at Austin.
Yettram, A. L., & Robbins, K. (1970). Anchorage zone stresses in post-tensioned
uniform members with eccentric and multiple anchorages. Magazine of
Concrete Research, 73, 209–218.
Yettram, A. L., & Robbins, K. (1971). Anchorage zone stresses in axially post-
tensioned I-section members with end-blocks. Magazine of Concrete
Research, 74, 37–42.
Yun, Y. M. (2005). Evaluation of ultimate strength of post-tensioned anchorage
zones. Journal of Advanced Concrete Technology, 3, 149–159.
Zhou, L., Liu, Z., & He, Z. (2015). Further investigation of transverse stresses and
bursting forces in post-tensioned anchorage zones. Structural Concrete,
16, 84–92.

Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in pub-
lished maps and institutional affiliations.