6.1. Settlements
The main reason for reinforcing a road or railway with piles, apart from stabilising the embankment, is to avoid excessive settlements (downward vertical displacements). The numerically simulated settlements on the road in this study were analysed and compared with the criteria for serviceability limit states set by the Swedish Transport Administration in TK Geo 13 [
16]. The criteria of interest are the total settlements and the drainage gradient of the road pavement. The total settlement of the pavement is limited to 35 cm for the type of road studied in order to avoid seasonal flooding. The drainage gradient is defined as the inclination of the pavement between the crest and the side of the road. For the studied road, the design drainage gradient is 3.6%. The decrease in design drainage gradient due to settlements is limited to 1.1%, i.e., the final drainage gradient must be at least 2.5% to maintain sufficient water drainage of the pavement surface.
The settlements analysed were numerically computed on the crest of the road (point A), at the side of the road (point B), and at the embankment toe (point C), as shown in
Figure 1. The differential settlement between points A and B gives a measurement of the change in drainage gradient. Heaving at point C could indicate a potential problem with the stability of the embankment. Differential settlements between point C and points A and B might lead to damage to the geosynthetic reinforcement due to large strains.
In
Figure 9, the settlements
at points A, B and C were plotted against increasing centre-to-centre pile spacing (
) from 0.8 to 2.0 m for floating and semi-floating pile groups with square and triangular pile arrangements. No significant difference in
was found between the modelled square and triangular pile arrangements. Furthermore,
in the pavement (points A and B) remained within the serviceability limit state of 35 cm. The maximum values of
for floating and semi-floating pile groups were obtained on the crest (point A) of the embankment. As expected, the floating pile groups settled more than the semi-floating pile groups.
The magnitude of the differential settlements in
Figure 9 between the three points A, B and C did not vary much with the pile spacing
s. The largest changes in the differential settlements were observed when increasing the value of
s from 1.4 to 1.5 m and from 1.8 to 2.0 m. In both cases,
increased less at point C than at points A and B. The reason is an unchanged number of piles in the two cases at the same time as
s increased (see
Table 1), bringing the outermost pile closer to the embankment toe. The differential settlement between points A and B was approximately the same (about 1.5 cm) for both types of pile group and for all values of
. The decrease in drainage gradient was thus approximately 0.3% for all models, which is within the serviceability limit state of 1.1%. No local differential settlements (depressions) between the piles were observed in the pavement in the simulations. For the semi-floating pile group, the differential settlements between points B and C increased from 1.3 cm when
= 0.8 m to 2.0 cm when
= 2.0 m. For the floating pile group, it remained almost constant at 1.9 cm for the full range of
. The differential settlements observed in
Figure 9 are not of a magnitude that could damage the GR by plastic strains. No heaving was observed at point C, which is good from a stability point of view.
The results for 𝑢
𝑣 shown in
Figure 9 were compared with
for the case of no piles in order to quantify the reduction in
from pile-reinforcing the embankment. Two reference simulations were done without piles for the corresponding soil profiles of a floating and a semi-floating pile group;
is the average value of the settlement reduction in points A, B and C defined as a percentage
where
is
at point A on the crest of a piled embankment, and
is
at point A on the crest of an embankment without pile reinforcement. The quantities
,
,
, and
follow the same notation for points B and C. For an embankment with no pile group support,
at points A, B and C was 19, 18 and 12 cm, respectively, for the case of floating piles. The respective values were 17, 15 and 9 cm for the case of semi-floating piles. In
Figure 10a,
was plotted against
for floating and semi-floating pile groups with square and triangular pile arrangements. As expected, a semi-floating pile group reduced settlements more than a floating pile group did. However, increasing the value of
from 0.8 to 2.0 m had twice as large an influence on
for semi-floating piles as compared to floating piles. For
= 0.8 to 1.2 m,
remained nearly constant at 43% for the floating piles. Furthermore, for the floating piles, an increase in
s from 1.2 to 2.0 m reduced
almost linearly from 43% to 28%. For the semi-floating pile groups,
changed almost linearly over the range of modelled values of
, decreasing from 79% at
= 0.8 m to 45% at
= 2.0 m.
The resource efficiency of each installed pile can be quantified by dividing
with the number of piles per metre of road
(see
Table 1). The resulting ratio,
, gives an indication of how much each installed pile contributes to the settlement reduction for a given value of
. In
Figure 10b,
and
were plotted on separate vertical axes against
for floating and semi-floating pile groups with square and triangular pile arrangements.
increased with the increasing value of
, meaning that each individual pile contributed more to reducing the settlements as
was reduced. For both the floating and semi-floating piles,
increased almost linearly from
= 0.8 to 1.4 m. For
> 1.4 m,
still increased, although more slowly than between the simulated values of
. The gain in resource efficiency when increasing the value of
was less for
> 1.4 m compared to what was achieved in the interval 0.8 ≤
≤ 1.4 m. Thus, the pile spacing where
= 1.4 m is interesting as a design value. As previously mentioned, the pile spacing criterion in TK Geo 13 [
16] is 0.8 ≤
≤ 1.2 m. By increasing
above the maximum value, from 1.2 to 1.4 m, the simulated value of
increased by almost 1.3 times for both the floating and semi-floating piles. The value of
decreased from 9.2 to 6.4 when increasing
from 1.2 to 1.4 m, reducing the number of piles by almost one third.
Based on these results, it can be concluded that the pile arrangement, square or triangular, has no significant influence on the settlements in the road pavement or embankment toe of a lightly piled embankment. Furthermore, the serviceability limit state criteria in TK Geo 13 [
16] regarding total settlements and drainage gradient were not exceeded for the modelled range of
. It was also found that changes in
had greater influence on the magnitude of the settlements if the pile group is semi-floating instead of floating. The settlements with a floating pile group were approximately the same for any
within the pile spacing criterion 0.8 ≤
≤ 1.2 m in TK Geo 13, as seen in
Figure 9. Thus, it might only be beneficial to install the piles with
< 1.2 m if the piles can be driven down to a firm soil layer or bedrock.
6.2. Arching
In order to evaluate the arch formation in the modelled lightly piled embankment, the direction of the major principal effective stress (
) and the distribution of the vertical effective stress (
) were both analysed. The analysis specifically focused on comparing the arches formed for the square and triangular pile arrangements. As discussed in the literature study and shown in
Figure 5, the
-vectors align tangentially to the formed arch, and the value of
is lower underneath the arches and higher around the piles.
Figure 11 shows the distribution of
in the lower part of the embankment for a semi-floating pile group with
= 1.4 m. A detailed view of the distribution of
and the
-vectors, projected on a horizontal plane, are shown as red lines in
Figure 12. The length of a line corresponds to the relative magnitude of a
-vector projection. The greater the magnitude (i.e., the longer the line), the greater the stress rotation and arching.
Figure 11a,c, and
Figure 12a,c show the arch formation of the modelled piled road embankment, from now on referred to as the sloped embankment, for square and triangular pile arrangements, respectively.
Figure 11b,d, and
Figure 12b,d show the arch formation of an embankment with uniform embankment height for square and triangular pile arrangements, respectively. For all four cases, three horizontal planes were chosen at levels 1–3 marked in
Figure 11: plane 1 just above the pile heads, plane 3 near the top of the formed arches, and plane 2 in the middle of planes 1 and 3. The values of
presented in
Figure 11 are the average of the values of
on the three planes. Note that the legend in
Figure 11 has a scale of 0 ≤
≤ 150 kPa. The maximum value of
was approximately 1500 kPa on top of the pile heads below the embankment crest. Increasing the value of
above 1.4 m increased the value of
on top of and near the piles, and vice versa when decreasing the value of
below 1.4 m.
The two uniform embankments (
Figure 11b,d) were modelled to compare the numerical results of a sloped embankment with what is generally assumed in analytical calculations. When analytically calculating the ratio of load carried by the piles to the total load analytically, the combined load of the embankment fill weight and the traffic load is, in general, assumed to act vertically on the pile group as a uniform load. The height of the embankment is usually set equal to the height between the pile heads (or pile caps) and the embankment crest. Thus, the GRPSE design is on the safe side by assuming the largest possible load on the pile group and GR. The height of the two modelled uniform embankments was set to 2.50 m, equal to the average height of the embankment subjected to the traffic load
(see
Figure 11). The uniform embankment was modelled on top of seven pile rows (each row was seven piles wide) to simulate an infinitely vast piled embankment.
Figure 11a,b and
Figure 12a,b show the arch formation for the square pile arrangements. As expected, the arches formed symmetrically for the uniform embankment (
Figure 11b and
Figure 12b), with equal stresses in the transversal and the longitudinal road directions. The direction and magnitude of the
-vector projections show significant stress rotation between the piles in the uniform embankment. Zones with reduced value of
(destressed zones) formed between two adjacent piles, and a greater destressed zone formed between four piles, resulting in a square arch base (area under the arch). As seen in
Figure 11a and
Figure 12a, the resulting stress rotation between the piles in the longitudinal road direction for the sloped embankment was similar to the results for the uniform embankment. There was, however, a lack of stress rotation in the transversal road direction under the crest of the sloped embankment. Subsequently, the values of
were significantly greater along each pile row than between pile rows, especially below the embankment crest. Thus, the arches formed asymmetrically for the sloped embankment with square pile arrangement, with arches forming primarily in the longitudinal road direction.
Figure 11c,d and
Figure 12c,d show the arch formation for the triangular pile arrangements. As mentioned, the main reason the triangular pile arrangement was set as the design criterion for lightly piled embankments in TK Geo 13 [
16] was that a triangular pile arrangement would, in theory, create more stable arches than a square pile arrangement, in particular due to a shorter longest diagonal between the piles as the arch base is assumed to be triangular. As shown in
Figure 11d and
Figure 12d, the uniform embankment height resulted in near equal arch formation between three piles and equilateral triangular arch base. The slightly rhombus-shaped arch formation was due to an isosceles triangular pile arrangement. For the sloped embankment (
Figure 11c and
Figure 12c), the base of the arch was distinctly rhombus-shaped instead. The destressed zone was concentrated in the middle of the four piles. Thus, the longest span of the arch formed between four piles in a triangular pile arrangement was the long diagonal of the resulting rhombus arch base (
). This diagonal is longer than both the side of the expected triangular arch base (
) and the diagonal of the square arch base (
). Similar to the results of the square pile arrangement, the arching was asymmetrical with the arches forming primarily in the longitudinal road direction, as seen in
Figure 12c. Comparing
Figure 12a,c, the arching along the diagonal between two piles in the square pile arrangement is greater or equal to the longest diagonal between two piles in the triangular pile arrangements.
In order to compare the stability and shape of the arch for the square and triangular pile arrangement, cross-sections were chosen along the diagonal between two adjacent piles (marked A-A and B-B in
Figure 12).
Figure 13 and
Figure 14 show the distribution of
and
-vectors in the cross sections marked for the sloped and uniform embankments. As shown in
Figure 13 and
Figure 14, the shape of the arch was triangular for both the square and triangular pile arrangement, as well as for both sloped and uniform embankment. The shape of the arch in the uniform embankment was, for every modelled value of
, similar to the geotechnical centrifuge trapdoor test results by Rui et al. [
52] of an unreinforced embankment with
= 3.0. The
-vectors formed an enclosed triangular arch only for the triangular arrangement (
Figure 13b and
Figure 14b), translating to a more stable arch. However, the arching was similar for square and triangular pile arrangements in terms of the
distribution near the subsoil. The distribution of
in the embankment was more similar between the square and triangular pile arrangement for the uniform embankment (
Figure 14) than for the sloped embankment (
Figure 13). Furthermore, the results show that there was less overall arching in the sloped embankments than the uniform embankments. The value of
was greater between the piles on top of the GR and the subsoil for the sloped embankment (
Figure 13) than for the uniform embankment (
Figure 14). The direction of the
-vectors in
Figure 13 and
Figure 14 show that semi-circular arches formed underneath the GR in the uniform embankment, unlike in the sloped embankment.
To determine the effects of embankment model size on arch formation, the modelled piled road embankment was lengthened by seven pile rows and the road was widened from two to four lanes (adding six piles to each row). Lengthening the embankment made no difference to the arch formation. Widening the embankment still resulted in asymmetrical arch formation. The arches formed asymmetrically for all simulated values of
(
Table 1). Modelling the pile group as floating gave similar stress distribution and arch formation as modelling the pile group as semi-floating.
The main contributing factor to the asymmetrical arch formation in the sloped embankment is believed to be the spreading effect of the slopes. The effect is caused by the horizontal component of the earth pressure and the horizontal spreading of the traffic load over depth. This creates a horizontal load component acting on the slope as well. Fahmi Farag [
53] is one of many, and in particular one of the more detailed, studies of the spreading effect in geosynthetic-reinforced piled embankments (GRPSE). In the GRPSE design, the spreading effect is often taken into account when calculating the tensile force in the GR, but not when estimating the arch formation. This is also the case of the analytical methods presented in this paper. Thus, by assuming a uniform embankment, the spreading effect of the embankment slopes and subsequent asymmetrical arch formation is excluded in the GRPSE design.
6.3. Load Transfer
The piles in a GRPSE construction need to be able to carry the load that is transferred onto the piled heads by the arching and membrane effects (), as well as any load from negative shaft friction along the pile shaft. Thus, it is important to analyse at which value of s the total axial pile load () exceeds the structural bearing capacity of the timber piles; it is also of interest to analyse in order to quantify the arching and membrane effect.
The average distribution of
in the three modelled pile rows (
Figure 3) is summarised in
Figure 15 for both the floating and semi-floating pile groups, with 0.8 ≤
≤ 2.0 m. The difference in
between the square and triangular pile arrangements was insignificant (<3%) for the range of
. Thus, the values of
shown in
Figure 15 are the average results of the square and triangular pile arrangements. The difference in the value of
between the floating and semi-floating pile groups is mostly due to a difference in maximum pile toe resistance, set to 0 kN for the floating pile group and 14 kN for the semi-floating pile group. In addition, due to the slope of the embankment the maximum value of
in a pile row (
) was observed underneath the crest of the embankment, with the exception of
= 0.8 m and
= 1.0 m for the floating pile group. The distribution of
along a pile row was irregular, with a larger force acting on the outermost piles than on the second-outermost piles,
≤ 1.4 m for the floating pile group and
≤ 1.0 m for the semi-floating pile group. As previously shown, the resource efficiency in terms of settlement reduction (
) was almost linear (
Figure 10b) for
≤ 1.4 m, both for the floating and the semi-floating pile groups. Along with the distribution of
in
Figure 15, this suggests that the load bearing characteristics of the pile group were those of a combined unit for
≤ 1.4 m and those of single piles for
> 1.4 m.
The pile utilisation ratio as percentage (
), defined as
divided by the structural bearing capacity of 106 kN, is plotted against
in
Figure 16.
changed almost linearly for 0.8 ≤
≤ 2.0 m, increasing from approximately 33% at
= 0.8 m for all pile groups to approximately 114 and 120% at
= 2.0 m for the floating and semi-floating pile groups, respectively. The structural bearing capacity was reached (
≥ 100%) in one or more piles, for
> 1.6 m in the case of semi-floating piles and for
> 1.7 m in the case of floating piles. The maximum allowable value of
, in terms of structural bearing capacity, depends on the chosen factor of safety (
) defined as the ratio of the structural bearing capacity to
. Setting
equal to 1.2 or 1.3, for example, results in a maximum allowed
of roughly 88 kN or 81 kN (
= 83% or 76%), respectively. This corresponds to an
s equal to 1.5 m (
= 1.2) or 1.4 m (
= 1.3).
The average load transferred onto the pile heads from arching and membrane effects (
) is plotted against
in
Figure 17a. The proportion of embankment weight and traffic load carried by the piles is often referred to as pile efficacy in the literature (Van Eekelen and Han, 2020). Pile efficacy can be seen as an efficiency of the GRPSE design. The average pile efficacy (
) plotted against
is shown in
Figure 17b and calculated as
where
is the average load over the tributary area (
), as seen in
Figure 3. The difference in
(
Figure 17a) between the square and triangular pile arrangements was ≤3%, leading to the conclusion that the pile arrangement had no significant effect on pile head efficacy. For the modelled uniform embankment with a semi-floating pile group, used for comparison in the previous chapter, the difference in
was <2% between the two pile arrangements. In addition, modelling the pile group as floating instead of semi-floating resulted in a <3.5% difference in
(
Figure 17a), and thus had no significant effect on pile head efficacy. The difference in
is due to the larger relative displacement between the subsoil and pile head in the semi-floating pile group than in the floating pile group, which lead to a greater arching effect. This partly explains the observed difference in
and
between the floating and the semi-floating pile group in
Figure 15 and
Figure 16. From the results, it can be concluded that the stress distributions shown in
Figure 11 and
Figure 12 are comparable or deemed equal for the square and the triangular pile arrangements and for both the floating and semi-floating pile groups for a given value of
.
As seen in
Figure 17b,
decreased approximately linearly over the range of
for all pile groups by almost 25%, from about 60% at
= 0.8 m to above 35% at
= 2.0 m. When
= 100%, the sum of the embankment weight and traffic load is carried solely by the piles. When there is no arching or membrane effect,
< 5% for the present study. It should be noted that the lightly piled embankment design with no pile caps results in a small pile coverage ratio (
), defined as a percentage
where
is the equivalent square side length of a circular pile head with diameter
, and is calculated as
As previously mentioned,
was set to 200 mm for the timber piles, resulting in a range of
from 4.9% at
= 0.8 m to 0.8% at
= 2.0 m. Subsequently,
ranged from about 60% at
= 4.9% to above 35% at
= 0.8%. Similarly, Lai et al. [
35] performed DE simulations of a piled embankment with and without GR and for different values of
, and obtained a result of roughly 60% pile efficacy for the simulation with GR and
= 4.9% (
= 100 mm and
= 0.45 m). Briançon and Simon [
54] carried out a full-scale test on a piled embankment reinforced (
= 2.8%) with two layers of geogrid, similar to the load distribution layer in
Figure 4, and measured 74% pile efficacy. The simulated
in this study was about 55% at
= 3.1% (
= 1.0 m), which is significantly less in comparison. However, in the Kyoto Road project by Van Eekelen et al. [
2] on a monitored timber piled embankment (
= 4.4%), the pile efficacy was about 38%. Note that the Kyoto Road embankment fill material consisted of a sludge mixture with a friction angle lower than that of the embankment material simulated in this paper.
The results shown in
Figure 18 are the maximum and average values of the GR tensile load (
and
) above the pile heads, plotted against the modelled range of
. Modelling the piles in a triangular arrangement instead of a square arrangement increased
by 6% and
by 7% on average. The main reason for the difference is the choice of modelling the pile group as three pile rows. The offset in longitudinal road direction of every other pile (
Figure 3) in each pile row resulted in larger GR deflection and tensile load along the borders of the model for the triangular pile arrangement than for the square pile arrangements. Subsequently, adding more pile rows would decrease the influence of the tensile load along the model border on the average tensile load of the entire model. As seen in
Figure 18,
, was well below the GR tensile strength of 112 kN/m for the modelled cases. The largest resulting value of
for all modelled cases was almost 25 kN/m (1.1% strain), in the case of the triangularly arranged semi-floating pile group with
= 2.0 m. Briançon and Simon [
54] observed up to 0.7% strain in both geogrid layers at
= 2.8%, similar to the maximum simulated strain of roughly 0.6% at
= 3.1% (
= 1.0 m). The simulation results of
and
were larger on average by 20% and 23%, respectively, for the semi-floating pile groups than for the floating pile groups. The GR deflection and subsequent tensile load was greater for the semi-floating pile group than for the floating pile group due to the aforementioned difference in relative displacement between the two types of pile groups. Increased GR deflection leads to greater membrane effect and more load transferred onto the pile heads, as observed by Le Hello and Villard [
6], adding to the difference in
and
in
Figure 17 between the floating and semi-floating pile groups.
These results lead to the conclusion that the maximum allowable value of
in TK Geo 13 [
16] can be increased from 1.2 to 1.6 m without reaching the characteristic pile bearing capacity. In addition, the arching and membrane effect in a lightly piled embankment is not significantly influenced by the pile arrangement (square or triangular).
6.4. Comparison with Analytical Models
The numerical results were compared with results from analytical calculations, as no field measurements were conducted on the studied lightly piled embankment. Although the Extended Carlsson model [
23] is the analytical model recommended in TR Geo 13 [
16] for the design of lightly piled embankments, more advanced models were also included in the comparison in order to investigate their suitability for this type of GRPSE. Several analytical models for the design of GRPSE exist, and several of them represent national standards or recommendations. The analytical models considered most interesting for the present study are presented in
Table 3. The Swedish guidelines TR Geo 13 [
24] in combination with the Nordic guidelines [
23] were used for the calculations with the Extended Carlsson model [
26] and the SINTEF model [
27]. The British standards BS8006 [
22] were used for calculations with the H&R model [
29]. The French guidelines ASIRI [
30] were taken into consideration as a complement to BS8006, as the H&R model is also incorporated into ASIRI. The German recommendations EBGEO [
32] were used for the calculations with the Zaeske model [
32]. The calculations with the Concentric Arches (CA) model were done according to the Dutch guidelines CUR226 [
34]. The analytical models are ordered in
Table 3 according to their complexity, with the Extended Carlsson model being the simplest and the CA model the most advanced. The results of the analytical calculations were compared to the numerical results on pile head load (
), pile efficacy (
), and maximum GR tensile load (
).
Among the analytical models included in the comparison with the numerical model, the Zaeske model is the only one that supports a triangular pile arrangement. However, the only triangular pile arrangement supported in EBGEO is a square pile arrangement rotated 45 degrees. Thus, the formulas in Zaeske [
31] for a triangular pile arrangement were used, as they support the geometry of the TK Geo 13 triangular pile arrangement utilised in this study. Note that ASIRI presents geometrical recommendations but no analytical models for triangular pile arrangements. Furthermore, the Zaeske model is the only one of the included analytical models that takes into consideration an elevated GR, i.e., a vertical distance
between the pile head and the GR layer, when calculating the vertical stresses on the GR. The EBGEO [
32] contains recommendations for GR elevations; however, no elevation is included as a variable when calculating the vertical stress on the GR. Thus, the vertical stresses on the GR were calculated using the formulas derived by Zaeske [
31] for vertical stresses underneath the arch at
metres above the pile heads (where
= 0 m is in level with the pile heads), which are based on the same mechanical principles as the formulas prescribed in EBGEO. As previously mentioned,
= 0.45 m for the numerically modelled lightly piled embankment.
According to the analytical models included in the comparison and the corresponding guidelines (
Table 3), the calculated value of
is the sum of the tensile load in the GR due to the membrane effect (
) and the added tensile load in the GR due to the spreading effect (
) in the transversal road direction. The calculation of
differs between the included analytical models, and each method is briefly described further below.
is calculated similarly for all of the analytical models. In the Extended Carlsson, SINTEF and H&R model,
is calculated as
where
is the coefficient of active earth pressure,
is the embankment height of 2.5 m,
is the unit weight of the embankment material (calculated as 21.7 kN/m
3 for the partially saturated embankment), and
is the surcharge (traffic) load of 15 kPa, resulting in
= 18.1 kN/m. In the Zaeske and CA models,
is subtracted from
, which resulted to
= 13.1 kN/m.
The Extended Carlsson, SINTEF, and H&R models, and their corresponding standards or recommendations consider no subsoil support. A void is assumed to form underneath the GR due to soil consolidation and creep strains. That means that the entire embankment weight and surcharge load are transferred onto the pile heads, directly through the arching effect and indirectly through the membrane effect. Thus, = 100% for the Extended Carlsson, SINTEF and H&R models. is based on the geometry, load on the GR, GR stiffness, and maximum GR strain (). The value of is chosen based on the allowed short- and long-term strains. was set to 6% for the calculations with the Extended Carlsson, SINTEF, and H&R model, in accordance with the long-term criteria in TR Geo 13, the Nordic guidelines, BS8006, and ASIRI. For the H&R model calculations, the load on the GR was set to the minimum load prescribed in BS8006, which is equivalent to 15% of the embankment weight and the surcharge load.
Both the Zaeske and the CA models take subsoil support into consideration. , is calculated in a fashion similar to the Extended Carlsson, SINTEF, and H&R models, except that is calculated instead of selected. The value of is determined at vertical equilibrium between the loads acting on the GR and the sum of the GR reaction force and the subsoil support. The subsoil support is based on the modulus of subgrade reaction () of the subsoil, which can be set to zero if a void is expected to form below the GR. Based on oedometer tests with a constant rate of strain, was set to 304 kN/m3 for the Zaeske model and 151 kN/m3 for the CA model. was greater for the Zaeske model, as the model considers the stiffness and thickness of the embankment fill under the elevated GR.
Figure 19 shows the analytically and numerically determined values of
plotted against the pile centre-to-centre distance (
). The numerical results are the average values of
below the pavement for the square and triangular pile arrangements. The analytical results of the Zaeske model for the square and triangular pile arrangement are denoted as “Zaeske(Sq)” and “Zaeske(Tri)”, respectively. As previously mentioned, the analytical models included in this paper assume a uniform embankment height equal to the embankment crest height. Hence, the analytical models primarily estimate the load distribution below the pavement of the embankment. Thus, the numerical results shown in
Figure 19 are based on the average embankment height below the pavement. The weight of the soil beneath the arches (
) is excluded from the calculation results shown in
Figure 19 for the Extended Carlsson, SINTEF, and H&R models. Note that this underestimates the values of
for these three models, since
is normally set equal to the tributary load, defined as
, according to the respective guideline in
Table 3. For the Zaeske and CA models,
in
Figure 19 is the sum of the arched load (
) and the vertical component of
(
) along the pile head perimeter. Unlike in CUR226,
is not explicitly calculated in EBGEO. Thus, the Zaeske model results of
were calculated using the methods in CUR226, combined with Zaeske (2001). As shown in
Figure 19, the analytical results of
for the Extended Carlsson, SINTEF, and CA models were larger than the numerical results for all values of
. The results of the CA model almost followed the same trend as the numerical results, with an approximately 55% to 30% larger value of
at
= 0.8 to 2.0 m. The results of the Extended Carlsson and SINTEF models followed a similar trend as the numerical results, although the difference increased more than the CA results with increasing value of
. The results of the SINTEF model were constant for
≥ 1.6 m due to the calculated arch becoming taller than the embankment. The results of the H&R and Zaeske(Sq) models showed a similar trend, and were the closest to the numerical results for the full range of
. The results of the Zaeske(Tri) model of
were the closest to the numerical results for
≤ 1.3 m. Contrary to the numerical results, the results of the Zaeske model showed much smaller values of
in the triangular pile arrangement than in the square pile arrangement. In general, both of the numerical results of
were lower than the analytical results for
< 1.6 m. In the simulations, full arching or fully mobilised soil shear strength was not achieved, as the differential settlements between the pile heads and the subsoil were too small. In comparison, the analytical models assume full arching to be on the safe side with regard to the design axial load per pile. Thus, the numerical results were closer to the analytical results for
≥ 1.6 m, as the differential settlements between the pile heads and the subsoil became larger.
Figure 20 shows the resulting values of
for the modelled range of
. In the numerical results,
is taken as the average value of the quantity above those pile heads that are located below the pavement. As seen in
Figure 20, the Extended Carlsson, SINTEF, and H&R models give a much greater increase in
than the numerical results for an increasing value of
, which is due to the exclusion of subsoil support in these models. The results of the Zaeske(Sq) and CA models show a similar trend to the numerical results for the range of
, with the results of the CA model being closer in magnitude to the numerical results. The values of
were lower for Zaeske(Tri) than for Zaeske(Sq), since the sum of
and
in the transversal direction of the road was lower than
in the longitudinal direction with the triangular pile arrangement. The Zaeske(Tri) results of
decreased for
> 1.5 m as the long diagonal of the rhombus-shaped arch base (length of
) became too great in relation to the embankment height to support a stable arch. However, the numerical results showed no arch instability for any value of
s with the triangular pile arrangement. Overall, the numerical results of
were significantly lower than the analytical results, with the exception of Zaeske(Tri) for
= 2.0 m. This was firstly due to the active earth pressure not being fully mobilised in the numerical simulations, giving lower values of
in the analytical models. Secondly, the magnitudes of GR deflection were less in the simulation than in the analytical calculations, and resulted in lower values of
due to greater simulated subsoil support.
In conclusion, the results of the analytical calculations show that the numerical simulations are reliable in terms of load distribution in the embankment. The numerical and analytical results of were in good agreement overall. The numerical results of are deemed reasonable, though low compared with the analytical results. As expected, the more advanced analytical models (Zaeske and CA) were better at predicting the numerical results. The results of the Zaeske model came closest to the numerical results of . However, the difference in the results between the square and triangular pile arrangements was larger for the Zaeske model than for the numerical models. The CA model was the most consistent of the analytical models in the comparison, showing almost the same trend as the numerical results for not only but also . In contrast, the Extended Carlsson model greatly overestimated the numerical results, especially the results for . Thus, it would be more suitable to use the CA model instead of the Extended Carlsson model for the design of lightly piled embankments. However, due to its simplicity, the Extended Carlsson model could still be used for a first estimation if is substituted for .
The discussed guidelines (
Table 3) present geometrical limits for stable arch formation in terms of
, embankment height
and square pile head width (
) or pile/pile head diameter (
). These geometrical limits were compared with the hypothesis of an increased maximum allowable value of
(
) in TK Geo 13 [
16] for lightly piled embankments, currently at 1.2 m. EBGEO presents geometrical limits for both square and triangular pile arrangements.
Table 4 shows a summary of
for
= 2.5 m according to each of the studied guidelines;
is the diagonal centre-to-centre pile spacing, equal to
for the square pile arrangements and
for the triangular pile arrangements (longest diagonal between four piles). The results show that
= 1.6 m satisfies all the geometrical limits in
Table 4. For square and triangular pile arrangements,
equals to 2.2 or 1.6 m, respectively. This supports the hypothesis that
could be greater than 1.2 m for lightly piled embankments.