0Centrifuge tests on strip footings on sand with a weak layer
Ziccarelli, Valore, Muscolino and Fioravante
Geotechnical Research, 2017, 4(1), 47-64 http://dx.doi.org/10.1680/jgere.16.00021 Paper 16.00021
Received 18/12/2016; accepted 10/02/2017 Published online 15/03/2017
Keywords: geotechnical engineering/models (physical)/shallow foundations
Published with permission by the ICE under the CC-BY license. (http://creativecommons.org/licenses/by/4.0/)
Institution of Civil Engineers
publishing
Centrifuge tests on strip footings on sand with a weak layer
^ Maurizio Ziccarelli PhD
Assistant Professor, Dipartimento di Ingegneria Civile, Ambientale, Aerospaziale, dei Materiali, Universita degli Studi di Palermo, Palermo, Italy
^ Calogero Valore
Full Professor, Dipartimento di Ingegneria Civile, Ambientale, Aerospaziale, dei Materiali, Universita degli Studi di Palermo, Palermo, Italy
^ Sandro Rino Muscolino PhD
Research Fellow, Dipartimento di Ingegneria Civile, Ambientale, Aerospaziale, dei Materiali, Universita degli Studi di Palermo, Palermo, Italy
J Vincenzo Fioravante PhD
Full Professor, Dipartimento di Ingegneria, Universita degli Studi di Ferrara, Ferrara, Italy
Tests on small-scale physical models of a strip footing resting on a dense sand bed containing a thin horizontal weak soil layer were carried out at normal gravity (1g). The results, reported in a companion paper, point out that the weak layer plays an important role in the failure mechanism and the ultimate bearing capacity of the footing if it falls within the ground volume relevant to the behaviour of the sand-footing system. The same problem was also investigated by means of centrifuge tests on reduced-scale models at 25g and 40g. The results of these tests, reported and discussed in this paper, confirm that failure mechanisms are governed substantially by the presence of the weak layer if its depth does not exceed a critical value and highlight marked scale effects involving the ultimate bearing capacity related essentially to the mean equivalent stress level in the soil beneath and around the footing. Equivalent bearing capacity factors, N*r, for footings on a dense sand bed containing a thin weak layer are derived from experimental results and are proposed in the paper.
Notation qlim ultimate (or limit) bearing pressure (at peak) or ultimate
B footing width bearing capacity
Bm base width of model qlim,0 ultimate bearing pressure (at peak) of footing on
BP base width of prototype homogeneous sand bed
Cu uniformity coefficient to thickness of weak layer
c L1p cohesion intercept of sand Zi depth from the ground surface of the weak layer
A relative density zm depth from the ground surface of the deepest point of
d particle diameter the failure mechanism
d50 mean particle size gdi dry weight of sand
E Young's modulus gd2 dry unit weight of weak layer
e void ratio gsl specific weight of sand
eo initial void ratio 8 angle of shearing resistance of the footing-sand interface
emax maximum void ratio dl angle of friction of the glass-sand interface
emin minimum void ratio e emersion angle of the failure surface; QL and Sr on the
Gsi specific gravity of sand left and right sides of the footing, respectively
g gravity acceleration n Poisson's ratio
Ko coefficient of Earth pressure at rest p settlement of the footing
L footing length p* density (or volumic mass)
^m lateral extent of failure mechanism plim settlement of the footing in correspondence of qlim
N ratio between the centrifuge test and gravity s ' normal effective stress
accelerations sv vertical effective stress
n porosity t shear stress
no initial porosity fl angle of shearing resistance of sand
Q vertical load applied to the footing flcv angle of shearing resistance of sand at constant volume
q mean vertical pressure acting on the footing base (or at critical porosity)
Volume 4 Issue GR1 with a weak layer
Ziccarelli, Valore, Muscolino and Fioravante
f1p peak angle of shearing resistance of sand
mean equivalent angle of shearing strength f2p angle of shearing resistance of the weak layer y1p peak dilation angle of sand
Introduction
Minor geological and geotechnical details can have great relevance for seepage and consolidation processes as well as for the movements and stability of natural and manmade geotechnical systems (Leonards, 1982; Rowe, 1972; Terzaghi, 1929). The simplest of such details is probably exemplified by a thin horizontal weak soil layer interbedded in a mass of stiffer soil. This problem was recently investigated with reference to the ultimate bearing capacity of strip footings by 1g small-scale model tests discussed in a companion paper (Valore et al., 2017).
These tests highlighted that the weak layer, despite its thinness, can markedly affect the failure mechanism and significantly reduce the ultimate bearing capacity. However, it is well known that 1g tests on reduced-scale models suffer from severe limitations due to scale effects associated, first of all, with the very low stresses in the granular soil of the model (e.g. de Beer (1965), Vesic (1975), Kimura et al. (1985), Bolton and Lau (1989), Kusakabe et al. (1991), Ueno et al. (2001), Zhu et al. (2001), Lau and Bolton (2011a)).
To investigate scale effects on failure mechanisms and on the bearing capacity, centrifuge experiments at enhanced gravity of 25g and 40g were carried out. It is well known that these kinds of tests also serve the important purpose of providing reliable physical data to verify numerical methods as clearly pointed out by Ng (2014). Ten centrifuge tests were performed.
The problem dealt with can be schematised with reference to the front view of the reduced-scale model shown in Figure 1. The strip footing is a rigid punch and rests on the surface of a dry
Q Footing
Sand 1 z-.
1 Weak layer 3:
1-L--I
Figure 1. Scheme for the formulation of the problem
sand mass in which a thin horizontal layer, t0 thick, made of a weaker material than sand, is interposed at depth zt. Plane strain conditions are assumed.
The results of the tests are reported and discussed in the present paper.
Instrumentation
Centrifuge
The small-scale model tests were carried out using the Istituto Sperimentale Modelli Geotecnici (ISMGeo) seismic geotechnical centrifuge, which is a beam centrifuge made up of a symmetrical rotating arm with a diameter of 6 m, a height of 2 m and a width of 1 m, which gives it a nominal radius of 2 m. The arm holds two swinging platforms, one used to carry the model container and the other the counterweight. During the tests, the platforms lock horizontally onto the arm to prevent the transmission of the working loads to the basket suspensions. An outer fairing covers the arm; the arm and fairing concurrently rotate to reduce air resistance and perturbations during flight. The centrifuge has the potential of reaching an acceleration of 600g at a payload of 400 kg. Further details can be found in the papers of Baldi et al. (1988), Fioravante (1999) and Fioravante et al. (2012). The dimensions of the tested models are length = 0-62 m, height = 0-28 m and width = 0-16 m.
A picture of a model at the end of a test is shown in Figure 2. Load application
The axial load is applied by a mechanical actuator that pushes the footing into the sand at a constant rate of displacement of 0-5 mm/min. The load is measured by a 50 kN hydraulic cell.
Figure 2. The physical model fixed to the basket at rest after the conclusion of a test
Volume 4 Issue GR1 with a weak layer
Zlccarelll, Valore, Muscollno and Floravante
Data acquisition system and measurements
Data are recorded by an automatic six-channel system, with an acquisition frequency of one record every 2 s. The radial acceleration is measured by a piezoresistive accelerometer. The settlement of two opposite vertices of the top face of the footing and of a point of the soil surface located at a distance from the centroid of the footing base of 17-8cm (equal to 4-45B) are measured by means of linear displacement transducers (LDTs). Images of the frontal face of the model are taken through a poly(methyl methacrylate) (PMMA) window by a video-recording colour digital camera.
Model preparation
The model was prepared in the ISMGeo geotechnical laboratory adjoining the centrifuge room and then placed aboard the centrifuge.
The foundation soil, apart from the weaker layer, consists of silica sand and is hereafter called sand B. The sand grains are subrounded to angular. Each soil model was reconstituted at 1g to the target void ratio by pluviating, in air, the dry sand into a rigid container by using a travelling sand spreader. The target density was obtained by calibrating the height of fall and the size of the spreader hole. The height of fall of 1 m was kept constant during the deposition in order to achieve uniformity of the soil density within the sand bed. The main characteristics of sand B are summarised in Table 1, where the average void ratio, e0, the porosity, n, dry density, yd1 and relative density, Dr of the models before acceleration are also indicated.
Table 1. Characteristics and initial index properties of sand B
From the mineralogical point of view, sand B is essentially composed of silica (more than 95%). There are, however, traces of feldspars and calcite.
In order to compare the results of 1 g and centrifuge tests and to study the scale effects, sand B was also used in 1 g tests on 40 mm wide footings (Valore et al., 2017), according to a well-established practice (Altaee and Fellenius, 1994; Kimura et al., 1985; Schofield, 1980; Toyosawa et al, 2013; Yamaguchi et al., 1977). It is worth noting that the ratio of the footing width to the mean particle size, B/d50 = 40/0-45 = 88-9, is larger than 50, which is the minimum value beyond which the particle size effect can be considered negligible, as suggested by many researchers (e.g. Mikasa and Takasa (1973), Ovesen (1975), Gemperline and Ko (1984), Kutter et al. (1988), Tatsuoka et al. (1991), Kusakabe (1995), Herle and Tejchman (1997), Toyosawa et al. (2013)).
The crushing of the sand particles during the tests has been considered negligible due to the grading and the mineralogy of the sand and the mean stress level existing in the soil volume involved in failure (Valore and Ziccarelli, 2009).
The peak strength failure envelope of sand B, obtained from direct shear tests on dry sand (Figure 3), is strongly curvilinear in the range of low normal stresses. For very low normal stresses, the shear strength parameters vary as follows: c'1p = 0 and fp > 50° for d <20 kPa; c'1p = 0 and f p = 50° for d <50 kPa; and C1p = 0 and f 1 p = 45° for d >50 kPa, f 1 p being the secant angle of shearing resistance. The relation between f 1 p and the normal effective stress,
Gs1 jst: kN/m eo no emin emax Dr: % kN/m dmax: mm d60: mm d50: mm d1(): mm Cu = d6o/dio
2-65 26 0-647 0-393 0-634 0-897 95 15-8 0-85 0-47 0-45 0-33 1-42
2-65 26 0-647 0-393 0-634 0-897 95 15-8 0-85 0-47 0-45 0-33 1-42
Minimum and maximum void ratios were determined according to ASTM standards D 4253-00 (ASTM, 2004a) and D 4254-00 (ASTM, 2004b)
tfp: 43°
100 80 60
01p: 60'
o'v: kPa (a)
40 60 80 a'v: kPa (b)
100 200 300 400 a'v: kPa (c)
Figure 3. Peak strength failure envelope from direct shear results on dry silica sand B: (a) all tests and (b) tests at low stresses. (c) Secant angle of peak shearing resistance p as a function of the effective vertical normal stress oV
Volume 4 Issue GR1 with a weak layer
Ziccarelli, Valore, Muscollno and Floravante
sV, is shown in Figure 3(c). The peak dilation angle, y{p, was obtained from the results of direct shear tests and decreases with normal stress, o^. Within the range of sV relevant to the present research, y{p varies from 17 to 21° and is in good agreement with Bolton's relation: y'p = 1-25 (f1p - f1cv) (Bolton, 1986). The angle of shearing strength at constant volume f1cv is 32°.
The average values of peak strength parameters at the average stresses relevant to the models tested in the centrifuge, modified for the effects of the intermediate effective principal stress for plane strain conditions (Meyerhof, 1951; Roscoe, 1970; Rowe, 1969; Tatsuoka et al., 1986a, 1986b), are summarised in Table 2.
Materials of the weak layer
The material used for the weak layer is the 'CM3' dry talc powder that was also used in tests at single gravity (Valore et al., 2017). The angle of peak shear strength, f2p, of this material was determined by direct shear tests and is equal to 27°. The cohesion intercept and dilation angles of the CM3 talc powder are negligible. The initial (i.e. before centrifuge testing) weak layer thickness was 5 mm.
Footing
The footing is made of aluminium and can be considered rigid; its width, B, and length, L, are 40 and 160 mm, respectively. Sandpaper was glued onto the footing base. The settlements of the footing are uniform. The peak angle of shearing resistance, 8', of the interface between the silica sand and the sandpaper was 42° (Valore et al, 2017). As the average peak shearing resistance angle, f1p, of silica sand was about 48° at the stress level relevant to the model tests, the sand-footing contact can be considered perfectly rough (Hansen and Christensen, 1969; Kumar and Kouzer, 2007), since 8'/f1p > 0-7 (Fioravante, 2002; Garnier and Konig, 1998; Kishida and Useugi, 1987; Lings and Dietz, 2005).
The values of the base width of the models (Bm) and prototypes (Bp) are summarised in Table 3.
PMMA-sand friction
To minimise the friction, the sand-PMMA interface was lubricated with silicone oil. Although the PMMA-sand friction is not nil, it is believed that the actual deformation state can be approximately considered two dimensional.
Model formation
The sand was poured into the test box by using the dry pluviation procedure. Through a proper apparatus with a 2 mm slit, a
horizontal speed of 10 cm/s and a constant fall height of 100 cm, it was possible to obtain a dry soil with a uniform initial unit weight, g1d, of 15-8 kN/m3. The weak layer was formed with dry talc powder (type CM3). It was formed by pouring a known weight of talc powder into the test box, and then it was gently compressed by means of a flat wooden pestle in order to obtain a regular layer 5 mm thick.
Testing procedure
The following sequence was used.
■ The centrifuge was accelerated up to the selected acceleration, a (25g or 40g); in this phase, the footing was suspended over the soil.
■ The soil model densification due to self-weight at constant acceleration was monitored.
■ At the end of the in-flight densification (i.e. end of soil surface settlement, as measured by an LDT), a deceleration/ acceleration cycle of a/2 (half of testing acceleration) was carried out. At the end of this phase, the settlement of the top surface of the homogeneous sand bed was 0-5 cm for both tests at 25g and 40g. This settlement implies an increase in the dry unit weight from 15-8 to 16-1 kN/m3.
■ The footing was gently lowered until contact with the model surface was achieved. The loading test was then performed. The footing was pushed at a constant rate of displacement of about 0-5 mm/min, and the axial load Q was measured by a load cell until the failure load was attained (ultimate bearing capacity). Subsequently, the loading continued until a settlement of about 25 mm (=0-6B) was attained.
The duration (after the start of the loading proves) of each test ranged from 15 to 40 min.
The homogeneous sand bed test and the tests on sand containing a weak layer made of dry talc powder were, of course, drained since both the sand and the talc were dry.
Results
The results of tests are summarised in Tables 4 and 5 for a = 25g (N = a/g = 25) and a = 40g (N = 40), respectively.
Failure mechanisms
Observed failure mechanisms are sketched in Figure 4 for tests at 25g (N = 25) and in Figures 5 and 6 for tests at 40g (N = 40). Some photographs of the models are shown in Figures 7 and 8.
Table 3. Widths of the bases of footings for models (Bm) and prototypes (Bp)
N a Bm: m Bp: m
1 1g 0-04 0-04
25 25g 0-04 1-00
40 40g 0-04 1-60
Table 2. Shear strength parameters of sand B for plane strain conditions
c'lp: kPa eV ° elcv: ° Xp: °
0 46-49 32 17-21
yip, peak dilation angle
Table 4. Results of centrifuge tests performed at acceleration a = 25g on strip footing models resting on a sand bed containing a thin horizontal weak layer
Weak layer Test Zi/Bm Î0: mm /m: mm /m/Bm zm: mm Zm/Bm 6L:° qlim: kPa 9lim/9lim,0 qm,lim: mm qm,lim/Bm
Dry CM3 talc powder C02 1 7 106 265 40 1-00 39 36 2118-2 0-55 5-39 0-13
C03 2 5 161-4 4-04 80 2-00 36 — 2413-9 0-62 5-95 0-15
C08 3 5 116-9 2-92 57-3 1-43 — 39 3117-3 0-81 7-09 0-18
Homogeneous sand bed C01 — — 104 26 40 1 30 31 3869-7 1 5-39 0-15
Results of test C01 on homogeneous sand bed is also reported for comparison
Bm = 40 mm, footing model width; z„ depth of the top surface of the weak layer; t0, thickness of the weak layer; qlim, ultimate bearing capacity; qlim 0 = 3869-7 kPa, ultimate bearing capacity for the homogeneous sand case; pm lim, settlement of footing model corresponding to ultimate bearing capacity; /m, maximum lateral extent of failure mechanism; zm, maximum depth of failure mechanism; q and 0R, emersion angles of the failure surface on the left and right sides of the footing, respectively
Table 5. Results of centrifuge tests performed at acceleration a = 40g on strip footing models resting on a sand bed containing a horizontal thin weak layer
Weak layer Test Zi/Bm Î0: mm /m: mm /m/Bm zm: mm Zm/Bm 6L:° qlim: kPa 9lim/9lim,0 qm,lim: mm qm,lim/Bm
Dry CM3 talc powder C05 1 7 101 253 40 1-00 38 36 2575-1 0-52 538 0-13
C06 2 7 179 4-48 80 2-00 33 — 3163-9 0-62 5-95 0-15
C07 2-95 7 180 4-5 71 1-77 — 40 4285-1 0-87 8-59 0-21
C10 0-5 6 47-3 1-18 35 0-88 1 46 45 3144-7 0-64 1 7-35 0-18
Homogeneous sand bed C04 104 26 63 30 31 4937-2 662 0-17
C09 — — 155 387 57 1-43 — 39 4767 1 6-49 0-16
Results of tests C04 and C09 on homogeneous sand bed also reported for comparison
Bm = 40 mm, footing model width; z„ depth of the top surface of the weak layer; t0, thickness of the weak layer; qlim, ultimate bearing capacity; qlim0 = 4937 2 kPa, ultimate bearing capacity for the homogeneous sand case; pm,lim, settlement of footing model corresponding to ultimate bearing capacity; /m, maximum lateral extent of failure mechanism; zm, maximum depth of failure mechanism; q and 0R, emersion angles of the failure surface on the left and right sides of the footing, respectively
(a) (b)
Figure 4. Failure mechanisms observed in centrifuge tests performed at 25g (N = 25). Bm = 40 mm. (a) Homogeneous sand bed; (b) ratio z,-/6 = 1; (c) z,-/B = 2; (d) z,-/B = 3. The weak layer is shown as double-thickness lines; the thin lines are initially horizontally aligned coloured sand particles
Figure 5. Failure mechanisms observed in centrifuge tests performed at 40g (N = 40). Bm = 40 mm. Homogeneous sand bed. Thin lines are initially horizontally aligned coloured sand particles
The photograph of test B63 performed at 1g (Valore et al., 2017) is included in Figure 8 for comparison.
General shear failure mechanism (Vesic, 1973) was observed in all experiments. Some observed failure mechanisms were not
symmetric. Every effort was made to ensure initial geometric and mechanical symmetry of the tested physical models. But symmetry at the macroscopic scale does not imply an always-perfect symmetry at the microscopic level. Non-symmetric mechanisms could originate from initial non-symmetric texture and non-perfectly symmetric distribution of the pores. Results of many physical model experiments on strip footings reported in geotechnical literature yield non-symmetric failure mechanisms, particularly for foundation soil consisting of sands (e.g. Muhs (1965), Yamaguchi et al. (1976), Kimura et al. (1985), Tatsuoka et al. (1991), Kusakabe (1992), Aiban and Znidarcic (1995), Tatsuoka (2001), McMahon and Bolton (2011)). It is worth mentioning that Yamaguchi et al. (1976) discovered non-symmetry in centrifuge models at 40g (Bm = 20, 30 and 40 mm, Bp = 0-80, 1-2 and 1-60) using radiography, while Kimura et al. (1985) demonstrated, also by means of X-ray techniques, that at 20g (Bm = 30 mm, Bp = 0-6 m) that the slip lines are not symmetric when the base of the footing is rough, while for smooth bases the slip lines are very nearly symmetric.
In the case of the homogeneous sand bed, the failure surfaces resemble that of Prandtl (1920), but their lateral extent is smaller; see Figures 4(a), 5 and 7. Failure surfaces emerge at the ground level at an average angle 0 to the horizontal of about 45° - yíp/2. When the weak layer is present, 0 ranges from 33 to 39° and is
(a) (b)
Figure 6. Failure mechanisms observed in tests performed at 40g (N = 40). Bm = 40 mm: (a) ratio z,-/B = 0-5; (b) z,-/B = 1; (c) z,-/B = 2; (d) z,/B = 2-95. The weak layer is shown as double-thickness lines; the thin lines are initially horizontally aligned coloured sand particles
Volume 4 Issue GR1 with a weak layer
Ziccarelli, Valore, Muscolino and Fioravante
that the stress state, the stress paths and the strain paths in the sand bed with a weak layer differ from the corresponding ones in the homogeneous sand bed. If 0-5Bm < zt < 3Bm, the failure mechanisms run partly along the weak layer. It is to be noted that for a depth of the weak layer of zi/B = 2, the lateral extension of the failure mechanism relative to the 1 g test is greater than that at enhanced gravity (see Figure 8).
The results of centrifuge tests confirm the great influence of the weak layer on failure mechanisms and on ultimate bearing capacity.
Figure 7. Failure mechanism for footing resting on homogeneous sand bed. Test C09 at 40g. Bm = 40 mm
Ultimate bearing capacity
The bearing pressure-normalised settlement curves are shown in Figures 9 and 10 for accelerations of 25g (N = 25) and 40g (N = 40), respectively. Results of some 1g tests (cf. Valore et al. (2017)) are shown in Figure 11 for comparison.
very close to 45° - y'p/2 (y'p = 15-20°). The failure surface crosses the weak layer if it is located at depth zi < 0-5Bm, both at 1g and at enhanced gravity, since its shear strength is high enough to transfer shear stresses to the underlying soil. In contrast, if Zi > 3Bm, the mechanisms run entirely within the upper sand layer and are almost identical to those pertaining to the homogeneous sand bed but the ultimate bearing capacity is lower than that of the homogeneous case. This difference is probably due to the fact
In all cases, the curves are characterised by a distinct peak corresponding to the ultimate bearing capacity, qlim. Beyond the peak, the applied pressure, q, undergoes a conspicuous, but not abrupt, decrease. The normalised settlement, pm,lim/Bm, when the weak layer is lacking, is about 14 and 15-5% for the centrifuge tests at 25g and 40g, respectively; in the presence of the weak layer, pm,lim/Bm varies from 14 to 18%, respectively, for the centrifuge tests at 25g and from 15 to 21% for the centrifuge tests at 40g. The
Test C05 z, /B = 1 40g (N = 40)
Test C06 z,/B = 3 40g (N = 40)
Figure 8. Typical failure mechanisms observed in tests on sand bed containing a weak layer, at 1g and at enhanced gravity. Bm = 40 mm. Note that the scales of the photographs are not the same. White broken lines are failure mechanism
a = Ng Test C02 C03 C08
N = 25 z, /B 1 2 3
Bm = 0-04 m
Bp = 1 m C01 J C08
pm/Bm: %
Figure 9. Bearing pressure-normalised settlement curves for tests performed at a = 25g. Results of test C01 performed on homogeneous sand bed are reported for comparison. Bm = 40 mm; Bp = 1 m; z,-, depth of the weak layer. Weak layer made of CM3 talc powder
C04 and C09: C04
tests on C07
homogeneous
sand bed
a = Ng
N = 40 C09
Bm = 0-04 m //
Bp = 1-6 m /
p C06
Test C05 C06 C07
z¡ IB 1 2 3
0 5 10 15 20 25
Pm^m: %
Figure 10. Bearing pressure-normalised settlement curves for tests performed at a = 40g. Results of tests C04 and C09 performed on homogeneous sand bed are reported for comparison. Bm = 40 mm; Bp = 1-6 m. z,-, depth of the weak layer. Weak layer made of CM3 talc powder
settlements corresponding to the peak therefore range from about 0T5Bm to 0-2Bm; their influences on qlim are negligible according to Ovesen (1975), Pu and Ko (1988) and Dijkstra et al. (2013).
The normalised applied pressure, q/Bp, in the function of the normalised settlement, pp/Bp, is reported in Figure 12 for tests on homogeneous sand bed and in Figure 13 for tests on sand bed with a weak layer.
The results plotted in Figures 12 and 13 show that the normalised applied pressure decreases as the stresses increase with the centrifuge acceleration (a = Ng), both in homogeneous soil and in the presence of the weak layer.
Figure 14 shows the bearing pressure q-normalised settlement pm/ Bm curves for a = 25g (N = 25) and an equivalent width of the prototype footing 40 x 25 mm = 1 m and for a = 40g (N = 40) corresponding to an equivalent prototype width of 1-60 m. These curves show a rather small-scale effect in compliance with the modelling concepts. Figure 14 also shows that the curves, at the same zjB, are characterised by fairly comparable stiffnesses and that the stiffness in the presence of the weak layer is smaller than that for the homogeneous case (Figure 14(a)). Moreover, the stiffness increases with zJB (Figures 14(b)-14(d)).
The ultimate bearing capacity, qlim, has been normalised with respect to the ultimate bearing capacity, qlim>0, relative to the case of the homogeneous sand bed. The values of qiim>0 are 3870 kPa for N =25 and 4937 kPa for N =40. The ratio qlim/qlim>0 is plotted against zjB in Figure 15. In this figure, the results relative to those of 1g tests are also plotted for comparison. The minimum value of qlim/qlim 0 is attained at zJB = 1 with a reduction of the ultimate bearing capacity from 45 to 50% compared to the homogeneous sand case for tests at 25g and 40g, respectively. For the single-gravity tests, this reduction is 46% and occurs at zjB = 1. qlim tends to qljm>0 at zJB larger than 4 for both 1g and centrifuge tests.
The results presented earlier demonstrate the strong influence of the presence of a thin weak layer on the ultimate bearing capacity, qlim, which can undergo reductions as high as 48% relative to weak layers made of CM3 talc powder with f2p = 27°. Larger reductions are expected for f2p < 27°. These results confirm those relative to tests performed at a = 1g on the same sand B (d50 = 0-45 mm) and on a coarser sand A (d50 = 0-95 mm) (Valore et al, 2017).
Scale effects
It is well known that the ultimate bearing capacity reduces with increasing footing size and with increasing mean stress level (e.g. Bjerrum (1973), de Beer (1970), Shiraishi (1990), Briaud and
& 150 -
B60: test on homogeneous sand bed
a = Ng N = 1
Bm = 0-04 m
Bp = 0-04 m
B67 B64 —B63
Test B62 B63 B64 B65 B66 B67
z, /B 1-5 2 3 1 4 0-5
Pm/Bm: %
1g 25g
y" C01
Tests performed on homogeneous sand
P-Bp: %
Figure 11. Bearing pressure-normalised settlement curves for tests at a = 1 g on sand B. Test B60, on homogeneous sand bed, shown for comparison. Bm = Bp = 0-04 m. z,-, depth of the weak layer. Weak layer made of CM3 talc powder
Figure 12. Tests on homogeneous sand. Normalised bearing pressure-normalised settlement curves at 1g (N = 1), 25g (N = 25) and 40g (N = 40)
Jeanjean (1994), Cerato and Lutenegger (2007), Kumar and Katri (2008), White et al. (2008), Chakraborty and Kumar (2016)).
The stress effects originate, first of all, from the marked curvature of the dense sand failure envelope that is particularly relevant at low stress. In small-scale physical models, the self-weight stresses are very or extremely low under 'normal-gravity' conditions; as a consequence, the angle of shearing resistance is higher and variable along the failure surface, contrary to what happens for real footings. Stress effects may also depend on the heterogeneity of the foundation soil and progressive failure. According to Muhs (1965), Hettler and Gudheus (1988) and Lau and Bolton (2011a), progressive failure may be considered marginal in small-scale model tests such as the present ones.
N* for centrifuge tests was computed with reference to gd = 16-1 kN/m3 (instead of gd = 15-8 kN/m3 pertinent to single-gravity tests) in order to account for the densification undergone by the sand during the densification cycle phase.
The results shown in Figures 16 and 17 and in Figures 12 and 13 clearly confirm the well-known scale effects relative to homogeneous soils, observed both in single-gravity and in centrifuge tests (Kimura et al, 1985; Kutter et al., 1988; Tatsuoka et al, 1991; Ueno et al, 1998; Yamaguchi et al, 1977). They also demonstrate that the scale effects are present and are important for footings resting on sand bed in which a weak layer is present at depths smaller than a critical value, zcrit, depending on the ratio of its shear resistance angle (f2p) and that of the sand (f'1p). zcrit for the tested physical models ranges from 4Bm to 4-5Bm.
Figures 16 and 17 show the results of the model tests at different accelerations (1g, 25g and 40g) on the same sand B (cp. Valore et al. (2017) for 1g tests). In Figure 16(a), the trend of the 'equivalent ultimate bearing capacity factor' Ny = 2qlim/(gBp) is plotted against the prototype width, Bp, while in Figure 16(b), the experimental results relative to footing on homogeneous sand were compared with other experimental data (Kimura et al., 1985; Yamaguchi et al., 1976) and with some theoretical solution (Brinch Hansen, 1970; Kumar and Kouzer, 2007; Terzaghi, 1943; Vesic, 1973). In Figure 17, the ultimate bearing pressure, qlim, is plotted against Bp.
In all the analysed cases, for footings resting on either a homogeneous sand bed or one with a weak layer, at the same relative density of the sand, the general failure mechanisms are similar but the lateral extent of the failure surface in centrifuge tests is smaller than that in 1 g tests. This is due to the fact that higher mean angles of shearing resistance operate in 1 g tests.
Back-analysis
The main aims of the numerical analysis are to back-calculate the mobilised mean equivalent constant angle of shearing resistance, f'1p, of the sand (Lau and Bolton, 2011b) corresponding to the
CQ Ä 0-2
\ (ig) \
Tests performed on homogeneous sand
..............\ (25g)
/ B63 \
(40g) .. C03
pp/ßp: %
..-■-.....(ig)
Tests performed on homogeneous sand
pp/ßp: %
Figure 13. Normalised bearing pressure-normalised displacement curves in the function of the ratio z,-/B for tests at 1g (N =1), 25g (N = 25) and 40g (N = 40): (a) z,-/B = 0-5, (b) z,/B = 1, (c) z,/B = 2 and (d) z,-/B = 3. Dotted lines refer to tests carried out on homogeneous sand and are plotted for comparison
ultimate bearing pressure, qlim, and to compare the features of the computed failure mechanisms against the experimental ones.
The reference scheme for the finite-element (FE) analysis is shown in Figure 18 along with the boundary conditions. Plane strain state and drained conditions are assumed. To avoid mesh-related dissymmetries, only the half model is analysed. The unit weight of the materials and the angle of shearing resistance of the
weak layer are assumed to be known. The cohesion intercept is always considered nil. The numerical simulations are carried out using the FE code Plaxis 2D (Plaxis, 2008), considering the geometry of the reduced-scale physical model of the soil-footing system for single-gravity as well as for centrifuge tests. The unit weight of soils has been set equal to g = Ngp* (N = 1, 25 or 40, p* being the density). The footing is subjected to a vertical load, Q, corresponding to an average bearing pressure, q, on the
Tests performed on homogeneous sand
Pm/Bm: %
z, /B = 1
Pm/Bm: % (d)
Figure 14 Bearing pressure-normalised settlement curves as a function of the ratio z,-/B for tests at 25g (N = 25) and 40g (N = 40). Tests C01 and C04 performed on homogeneous sand. At 25g and 40g, the equivalent prototype footing widths Bp are 1 and 1-6 m, respectively
soil-footing interface. Actually, a uniform vertical settlement of the footing base is imposed rather than the vertical load, Q, so duplicating the true experimental procedure and accounting for the high stiffness of the footing and for the roughness of its base (Lee et al., 2013). The simple elastic-perfectly plastic Mohr-Coulomb constitutive model with non-associated flow rule is used for soils as many other authors have (e.g. Bolton and Lau (1993), Yin et al. (2001), Potts (2003), Mabrouki et al. (2010),
Kumar and Khatri (2011)). Geometric variations of the system and their effects on the stress state in the soil are not taken into account. This hypothesis and the assumption of perfect plasticity imply that pre-peak hardening, post-peak strain softening and the dependence of angle of shearing resistance, f'1p, on stress level variations within the relevant soil volume are not taken into account, although there are more sophisticated constitutive models available that allow modelling of the post-peak behaviour
0 1 2 3 4 5
Figure 15 Normalised ultimate bearing capacity qiim/qnm,0 against normalised depth z,/B for different values of N = a/g of the weak layer. qlim,0, ultimate bearing capacity of the footing on homogeneous sand. Bm = 40 mm, footing width; z,, depth of the weak layer. Data relative to 1g tests (sand B) reported for comparison. Weak layers made of CM3 talc powder with f2p = 27°
of the soil-footing system (e.g. Potts and Zdravkovic (1999), Yin et al. (2001), Potts (2003), Siddiquee et al. (1999, 2001), Cassidy et al. (2002), Salgado (2008), Loukidis and Salgado (2009, 2011)). An equivalent constant mean value of f1* has been sought (Lau and Bolton, 2011b; Lee et al., 2013). Of course, under 1 g conditions, the shear strength parameters in the low-stress range strongly depend on the effective stress level; consequently, they vary, within the relevant soil volume, from 'lower' values in the zone beneath the footing (where the effective normal stresses are relatively large) to higher values within the passive zone, where the stresses in tested 1 g physical models are low or extremely low (Lau and Bolton, 2011a, 2011b). In contrast, in the case of centrifuge tests at 25g or 40g, the average stress intensity within the relevant soil volume is high enough so that only modest variations in f1* occur along the failure mechanism.
The dilatancy angle was always related to the peak shear strength, fi* by Bolton's relation y* = 1-25 (f1* - f1cv) (Bolton, 1986), in which f1cv = 32°. Progressive failure is not taken into account as suggested by the results in 1 g small- and large-scale physical models and in centrifuge tests carried out by Muhs (1965), Hettler and Gudheus (1988) and Lau and Bolton (2011a). The preceding hypotheses do not permit the prediction of the behaviour of the soil-footing system beyond the peak bearing pressure (Potts and Zdravkovic, 1999, 2001).
Homogeneous sand bed (Ny)
zJB = 3
Weak layer: CM3 talc powder "T
. Kumar and Kouzer (2007) Vesic (1973)
. Terzaghi (1943)
Present
research Brinch Hansen (1970)
—#— This research (Dr = 95%) —□— Kimura et al. (1985)
Toyoura sand (Dr = 89-97%) _<<:>_ Yamaguchi et al. (1976)
Toyoura sand (Dr = 87%) -1-1-1-
Figure 16. Strip footing on sand B with or without a weak layer at depth z,. Weak layer made of CM3 talc powder. (a) Equivalent ultimate bearing capacity factor Ny* = 2qlim/(ydBp) against Bp equivalent prototype width. (b) Nyfor strip footing on homogeneous sand. Comparison of results of the present research with other experimental data and with theoretical solutions for f1p values according to Figure 3(c)
Homogeneous sand bed
/ y / '
" M- -
-El- z,/B = 1 -o- z /B = 2 _A- z, /b = 3 -1-
Figure 17. Ultimate bearing capacity qlim against prototype width Bp. z,, depth of the weak layer made of CM3 talc powder
Homogeneous sand bed
First, the homogeneous sand bed-footing systems were back-analysed. A good match of experimental and calculated results was reached as far as the ultimate bearing pressure, qlim,0, the bearing pressure-settlement curve (up to qlim,0) and the failure mechanisms are concerned.
The following parameters have been considered for the sand: Young's modulus: E' = 125 MPa, Poisson's ratio v' = 0-15 and coefficient of Earth pressure at rest K0 = 0-4. Results of the back-analysis for tests performed at acceleration a = 40g (N = 40) are shown in Figures 19 and 20. The failure mechanism (Figure 19) closely resembles Prandtl's (1920) except for the angle of emersion at ground surface that nearly equals 45° -y1p/2 (instead of 45° - f1p/2). The small instability of numerical results in the pre-peak phase is due to the non-associativity of the constitutive model (Frydman and Burd, 1997).
to be compared to the results of centrifuge tests at 40g: C10, C05 and C06. For the sake of simplicity of the numerical simulation, the values of Young's modulus, Poisson's ratio and the dry unit weight, gi2, of the weak layer were assumed equal to those of the sand. The calculations were performed assuming, always (i.e. irrespectively of the depth of the weak layer), for sand f1p = 47-6° obtained from the back-analysis of test C04 at 40g on homogeneous sand bed. The numerical analyses confirm that when the weak layer is located at depth zt = 0-5B, the failure mechanism crosses the weak layer, develops through a radial shear zone within the underlying sand and runs upwards along an inclined plane inclined 45° (Figure 6(a)), which crosses the weak layer again before emerging onto the ground
Figure 18. Reference scheme for FE analysis. Unit weight of soil g = Ngp" (N = 1, 25 or 40). Q, vertical load or resultant of the pressures transmitted to the soil by the footing
The values of the equivalent mean angle of mobilised shear strength are f* = 47-8° (y* = 19-5°) for tests at 25g (N = 25) and f1* =47-6° (y{p = 19-7°) for tests at 40g (N = 40).
The calculated values of f* are in good agreement with the experimental results of direct shear tests (Figure 3) pertaining to the range of normal stress from 250 to 350 kPa. This range has been selected according to Meyerhof (1951), who suggested that the value of the mean normal stress, o0, along the failure surface is about 1/10 of the ultimate bearing capacity, qlim, and according to de Beer (1965), who proposed the following relation: s0 = 0'25qlim(1 - sin f1p).
Sand bed with weak layer
Despite its simplifications, the numerical analysis allows accurate identification of the failure mechanisms that match the experimental ones quite well, as shown, for example, by Figures 21-23, which are
Figure 19. Footing on homogeneous sand bed, N = 40. Incremental shear strains at failure. Compare with test C04
5000 i
4000 -
o- 3000 -
.g 2000 -
1000 -
a = Ng
Test C04 (N = 40)
■N = 40
Test C01 (N = 25)
Experimental FE modelling
-1-1-1-1
0 2 4 6 8
Settlement, p: mm
Figure 20. Footing on homogeneous sand bed. Comparison between numerical and experimental load-settlement curves (test C04, N = 40; test C01, N = 25)
surface. In the cases shown in Figures 22 and 23, the failure mechanisms found by the numerical analysis develop in part along the weak layer similarly to experimental results.
The values of qiim/qiim,0 calculated by using the equivalent constant strength parameters of the sand back-calculated for the homogeneous case for N = 25 and N = 40 are plotted in Figures 24 and 25, respectively.
It can be observed that the results of numerical analysis match the experimental data very well. These figures prove that at enhanced
Figure 22. Incremental shear strains at failure for z,/B = 1. Weak layer made of CM3 talc powder. Compare with test C05 (Figure 6(b)). N = 40
gravity, the mean equivalent angle of shearing strength, f1*, is not appreciably affected by the stress-related variability of f1p and by the depth of the weak layer, in contrast with what occurs for 1g tests (Valore et al., 2017). The earlier-mentioned figures also confirm the remarkable effect of the presence of the weak layer on the ultimate bearing capacity, which may undergo reductions as high as 50% when the weak layer is made of CM3 talc powder with a shearing resistance angle, f2p, equal to 27°. For this latter value of f2p, the experimental results along with the experimental ones clearly suggest that the critical adimensionalised depth zi/B closely approaches 4.
Conclusions
The influence of a horizontal thin weak soil layer interposed in a dense sand bed on the behaviour of a shallow strip footing loaded to failure was investigated by means of centrifuge tests on small-scale physical models. From the test results, the following conclusions can be drawn.
Figure 21. Incremental shear strains at failure for z,-/B = 0-5. Weak layer made of CM3 talc powder. Compare with test C10 (Figure 6(a)). N = 40
Figure 23. Incremental shear strains at failure for z,/B = 2. Weak layer made of CM3 talc powder. Compare with test C06 (Figure 6(c)). N = 40
0-75 _
0-50 -
0-25 -
_______
O Experimental (N = 25)
< 1 • 1 A FE modelling = 47-8° i 1 i 1 i
Figure 24. Results of back-calculations. Comparison between numerical and experimental values of qlim/qlim,0 against z,/B for N = 25. Weak layer made of CM3 talc powder
varies as a function of the ratio f'1p/f2p between the angles of shearing resistance of the sand and the material making up the weak layer. The failure surface cuts through the weak layer when the latter is located at small depths, zt, beneath the footing (zJB < 0-5); at larger depths (0-5 < zjB < 3), the weak layer controls the maximum depth of the mechanism, forcing it to run partly horizontally along the weak layer before going up through the upper sand layer.
The ultimate bearing capacity, qlim, is always lower than qlim,0 pertinent to the homogeneous sand bed. The experiments show a reduction in qlim of up to 50% for weak layers made of talc powder with an angle of shearing resistance of 27°; larger reductions are expected for smaller values of f2p.
The presence of a weak layer reduces the stiffness of the
load-settlement curve before
is reached. Numerical
simulations of the reduced-scale centrifuge physical model tests by FE analysis are able to capture the failure mechanisms and the ultimate bearing capacity correctly, even if the very simple constitutive Mohr-Coulomb model is used. Moreover, they point out that the equivalent mean constant value of the sand angle of shearing resistance in tests at enhanced gravity is little influenced by the location and the properties of the weak layer, in contrast with what happens for single-gravity tests.
The weak layer strongly influences both the failure mechanism and the ultimate bearing capacity, qlim, if its depth, zt, does not exceed a critical value of about 4B for the tested materials (sand and talc powder making up the weak layer). In general, this critical value
0 1 2 3 4 5
Figure 25. Results of back-calculations. Comparison between numerical and experimental values of qlim/qlim,0 against z/B for N = 40. Weak layer made of CM3 talc powder
The test results confirm those relative to single-gravity tests, reported in a companion paper (Valore et al., 2017), also carried out on a coarser sand and using materials for the weak layer with a wide range of angles of shearing resistance.
Scale effects, well known for homogeneous sands, also operate in sand beds containing a thin horizontal weak layer.
Bearing capacity factors, Ng , derived from results of tests on footing on homogeneous sand bed decrease with the prototype width and are in good agreement with other published experimental results and with theoretical solutions by Brinch Hansen (1970), Terzaghi (1943) and Vesic (1973). An equivalent bearing capacity factor, N*, has been derived for a sand bed containing a thin weak layer; it is lower than Ng and depends on the location and shearing resistance of the weak layer.
The results of tests carried out at different accelerations and of the back-analysis encourage confident numerical predictions of the behaviour of actual soil-footing systems of the kind dealt with in the paper.
Acknowledgements
This research was partly supported by the Ministry of Education, Universities and Research (Progetti di Rilevante Interesse Nazionale (Prin): Gallerie in 'sezione mista', Project 9908328717_009). This support is gratefully acknowledged. The authors wish to thank the laboratory technicians Dr Eng. A. Casella and Dr Geol. G. Sapienza of the geotechnical laboratory of the University of Palermo and the technicians of the ISMGeo for their assistance during the tests.
Volume 4 Issue GR1 with a weak layer
Ziccarelli, Valore, Muscolino and Fioravante
REFERENCES
Aiban SA and Znidarcic D (1995) Centrifuge modeling of bearing capacity of shallow foundations on sand. Journal of Geotechnical Engineering 121(10): 704-712, http://dx.doi.org/10.1061/(ASCE)0733-9410(1995) 121:10(704).
Altaee A and Fellenius BH (1994) Physical modeling in sand. Canadian Geotechnical Journal 31(3): 420-431, http://dx.doi.org/10.1139/t94-049.
ASTM (2004a) D 4253-00: Standard test methods for maximum index density and unit weight of soils using a vibratory table. ASTM International, West Conshohocken, PA, USA.
ASTM (2004b) D 4254-00: Standard test methods for minimum index density and unit weight of soils and calculation of relative density. ASTM International, West Conshohocken, PA, USA.
Baldi G, Belloni G and Maggioni W (1988) The ISMES centrifuge. In Proceedings of International Conference on Geotechnical Centrifuge Modelling, Centrifuge '88, Paris (Corté JF (ed.)). Balkema, Rotterdam, the Netherlands, pp. 45-48.
Bjerrum L (1973) Geotechnical problems involved in foundations of structures in the North Sea. Géotechnique 23(3): 319-358, http://dx. doi.org/10.1680/geot.1973.23.3.319.
Bolton MD (1986) The strength and dilatancy of sands. Géotechnique 36(1): 65-78, http://dx.doi.org/10.1680/geot.1986.36.1.65.
Bolton MD and Lau CK (1989) Scale effects in the bearing capacity of granular soils. Proceedings of the 12th International Conference on Soil Mechanics and Foundation Engineering, Rio de Janeiro, Brazil, vol. 2, pp. 895-898.
Bolton MD and Lau CK (1993) Vertical bearing capacity factors for circular and strip footings on Mohr-Coulomb soil. Canadian Geotechnical Journal 30(6): 1024-1033, http://dx.doi.org/10.1139/t93-099.
Briaud JL and Jeanjean P (1994) Load-settlement curve method for spread footings on sand. In Vertical and Horizontal Deformations ofFoundations and Embankments 1994 (Yeung AT and Félio GY (eds)). American Society of Civil Engineers, Reston, VA, USA, Special Publication No. 40, vol. 2, pp. 1774-1804.
Brinch Hansen J (1970) A revised extended formula for bearing capacity. Danish Geotechnical Institute Bulletin 28: 5-11.
Cassidy MJ, Byrne BW and Houlsby GT (2002) Modelling the behaviour of circular footings under combined loading on loose carbonate sand. Géotechnique 52(10): 705-712, http://dx.doi.org/10.1680/geot.2002. 52.10.705.
Cerato AB and Lutenegger AJ (2007) Scale effects of shallow foundation bearing capacity on granular material. Journal ofGeotechnical and Geoenvironmental Engineering 133(10): 1192-1202, http://dx.doi.org/ 10.1061/(ASCE)1090-0241(2007)133:10(1192).
Chakraborty M and Kumar J (2016) The size effect of a conical footing on Ny. Computers and Geotechnics 76: 212-221, http://dx.doi.org/10. 1016/j.compgeo.2016.03.010.
de Beer EE (1965) The scale effect on the phenomenon of progressive rupture in cohesionless soils. Proceedings of the 6th ICSMFE, Montreal, QC, Canada, vol. 2, pp. 13-17.
de Beer EE (1970) Experimental determination of the shape factors and bearing capacity factors of sand. Géotechnique 20(4): 387-411, http:// dx.doi.org/10.1680/geot.1970.20.4.387.
Dijkstra J, White DJ and Gaudin C (2013) Comparison of failure modes below footings on carbonate and silica sands. International Journal of Physical Modelling in Geotechnics 13(1): 1-12, http://dx.doi.org/10. 1680/ijpmg.12.00004.
Fioravante V (1999) Sui principî della modellazione fisica con particolare riferimento alla centrifuga geotecnica. Università della Calabria, Rende, Italy (in Italian).
Fioravante V (2002) On the shaft friction modelling of non-displacement piles in sand. Soils and Foundations 42(2): 23-33, http://doi.org/10. 3208/sandf.42.2 23.
Fioravante V, Ghiretti D, Prearo C and Lai C (2012) Static and dynamic centrifuge modeling landslide stabilization with large-diameter shafts. Soils and Foundations 42(2): 23-33.
Frydman S and Burd HJ (1997) Numerical studies of the bearing capacity factor Ng. Journal of Geotechnical Engineering 123(1): 20-29, http:// dx.doi.org/10.1061/(ASCE)1090-0241(1997)123:1(20).
Garnier J and König D (1998) Scale effects in piles and nails loading tests in sand. In Proceedings of the International Conference Centrifuge '98, Tokyo (Kimura T, Kusakabe O and Takemura J (eds)). Balkema, Rotterdam, the Netherlands, vol. 1, pp. 205-210.
Gemperline MC and Ko HY (1984) Centrifugal model tests for ultimate bearing capacity of footings on steep slopes in cohesionless soils. In Proceedings ofApplication ofCentrifuge Modelling to Geotechnical Design, Manchester, UK (Craig W (ed.)). Balkema, Rotterdam, the Netherlands, pp. 206-225.
Hansen B and Christensen NH (1969) Discussion on 'Theoretical bearing capacity of very shallow footings' by Larkins AL. Journal of Soil Mechanics and Foundation Division 95(SM6): 1568-1573.
Herle I and Tejchman J (1997) Effects of grain size and pressure level on bearing capacity of footings on sand. In Deformation and Progressive Failure in Geomechanics, IS-Nagoya '97 (Asaoka A, Adachi T and Oka F (eds)). Pergamon, Elsevier Science Ltd, Oxford, UK, pp. 781-786.
Hettler A and Gudheus G (1988) Influence of the foundation width on the bearing capacity factor. Soils and Foundations 28(4): 81-92, http:// doi.org/10.3208/sandf1972.28.4_81.
Kimura T, Kusakabe O and Saitoh K (1985) Geotechnical model tests of bearing capacity problems in a centrifuge. Géotechnique 35(1): 33-45, http://dx.doi.org/10.1680/geot.1985.35.L33.
Kishida H and Useugi M (1987) Tests of the interface between sand and steel in the simple shear apparatus. Géotechnique 37(1): 45-52, http:// dx.doi.org/10.1680/geot.1987.37.1.45.
Kumar J and Kouzer KM (2007) Effect of footing roughness on bearing capacity factor Ng. Journal of Geotechnical and Geoenvironmental Engineering 133(5): 502-511, http://dx.doi.org/10.1061/(ASCE)1090-0241(2007)133:5(502).
Kumar J and Khatri VN (2008) Effect of footing width on Ng. Canadian Geotechnical Journal 45(12): 1673-1684, http://dx.doi.org/10.1139/ T08-113.
Kumar J and Khatri VN (2011) Bearing capacity factors of circular foundations for a general c-f soil using lower bound finite elements limit analysis. International Journal for Numerical and Analytical Method in Geomechanics 35(3): 393-405, http://dx.doi.org/10.1002/ nag.900.
Kusakabe O (1992) Large-scale loading tests of shallow footings in pneumatic caissons. Journal ofGeotechnical Engineering ASCE 118(11): 1681-1695, http://dx.doi.org/10.1061/(ASCE)0733-9410 (1992)118:11(1681).
Kusakabe O (1995) Foundations. In Geotechnical Centrifuge Technology (Taylor RN (ed.)). Blackie Academic & Professional, London, UK, pp. 118-167.
Kusakabe O, Yamaguchi H and Morikage A (1991) Experiment and
analysis on the scale effect of Ng for circular and rectangular footings. In Centrifuge 91 (Ko HY and Mclean FG (eds)). Balkema, Rotterdam, the Netherlands, pp. 179-186.
Kutter BL, Moquette O'Leary L, Thompson PY and Lather R (1988) Gravity-scaled tests on blast-induced soil-structure interaction. Journal of Geotechnical Engineering 114(4): 431-447, http://dx.doi. org/10.1061/(ASCE)0733-9410(1988)114:4(431).
Lau CK and Bolton MD (2011a) The bearing capacity of footings on granular soils. II: Experimental evidence. Géotechnique 61(8): 639-650, http://dx.doi.org/10.1680/geot7.00207.
Lau CK and Bolton MD (2011b) The bearing capacity of footings on granular soils: I: Numerical analysis. Géotechnique 61(8): 627-638, http://dx.doi.org/10.1680/geot7.00206.
Volume 4 Issue GR1 with a weak layer
Ziccarelli, Valore, Muscolino and Fioravante
Lee KK, Cassidy MJ and Randolph MF (2013) Bearing capacity on sand overlying clay soils: experimental and finite element investigation of potential punch-through failure. Géotechnique 63(15): 1271-1284, http://dx.doi.org/10.1680/geot.12.P. 175.
Leonards GA (1982) Investigation of failures. Journal of Geotechnical Engineering Division 108(GT2): 222-283.
Lings ML and Dietz MS (2005) The peak strength of sand-steel interfaces and the role of dilation. Soils and Foundations 45(6): 1-14, http://doi. org/10.3208/sandf.45.1.
Loukidis D and Salgado R (2009) Bearing capacity of strip and circular footings in sand using finite elements. Computers and Geotechnics 36(5): 871-879, http://dx.doi.org/10.1016/jxompgeo. 2009.01.012.
Loukidis D and Salgado R (2011) Effect of relative density and stress level on the bearing capacity of footings on sand. Géotechnique 61(2): 107-119, http://dx.doi.org/10.1680/geot8.R150.3771.
Mabrouki A, Benmeddour D, Frank R and Mellas M (2010) Numerical study of the bearing capacity for two interfering strip footings on sands. Computers and Geotechnics 37(4): 431-439, http://dx.doi.org/ 10.1016/j.compgeo.2009.12.007.
McMahon BT and Bolton MD (2011) Experimentally observed settlements beneath shallow foundations on sand. In Proceedings of the 15th European Conference on Soil Mechanics and Geotechnical Engineering (Anagnostopoulos A, Pachakis M and Tsatsanifos C (eds)). IOS Press, Amsterdam, the Netherlands, vol. 1, pp. 749-754.
Meyerhof GG (1951) The ultimate bearing capacity of foundations. Géotechnique 2(4): 301-332, http://dx.doi.org/10.1680/geot.1951.2A 301.
Mikasa M and Takasa N (1973) Significance of centrifugal model test in soil mechanics. Proceedings of the 8th International Conference on Soil Mechanics and Foundation Engineering, Moscow, Russia, vol. 1.2, pp. 273-278.
Muhs H (1965) Discussion on: 'The scale effect on the phenomenon of progressive rupture in cohesionless soils by de Beer EE, Proceedings of the 6th ICSMFE, Montreal, vol. 2, pp. 13-17'. Proceedings of the 6th ICSMFE, Montreal, QC, Canada, vol. 3, pp. 419-421.
Ng CWW (2014) The state-of-the-art centrifuge modelling of geotechnical problems at HKUST. Journal of Zhejiang University - Science A 15(1): 1-21.
Ovesen NK (1975) Centrifugal testing to bearing capacity problems of footings on sand. Géotechnique 25(2): 394-401, http://dx.doi.org/10. 1680/geot.1975.25.2.394.
Plaxis (2008) Plaxis 2D, Version 8.6. Plaxis, Delft, the Netherlands. See http://www.plaxis.nl/ (accessed 15/02/2017).
Potts DM (2003) Numerical analysis: a virtual dream or practical reality? Géotechnique 53(6): 535-573, http://dx.doi.org/10.1680/geot.2003.53. 6.535.
Potts DM and Zdravkovic L (1999) Finite Element Analysis in
Geotechnical Engineering: Theory. Thomas Telford, London, UK.
Potts DM and Zdravkovic L (2001) Finite Element Analysis in
Geotechnical Engineering: Application. Thomas Telford, London, UK.
Prandtl L (1920) Über die härte plastischer körper. Nachrichten von der Königlichen Gesellschaft der Wissenschaften zu Göttingen. In Mathematisch-Physikalische Klasse. Weidmannsche Buchhandlung, Berlin, Germany, pp. 74-85 (in German).
Pu JL and Ko HY (1988) Experimental determination of bearing capacity in sand by centrifuge model tests. In Proceedings ofInternational Conference on Geotechnical Centrifuge Modelling, Centrifuge '88, Paris (Corté JF (ed.)). Balkema, Rotterdam, the Netherlands, pp. 307-311.
Roscoe KH (1970) The influence of strains in soil mechanics. Géotechnique 20(2): 129-170, http://dx.doi.org/10.1680/geot.1970.20.2.129.
Rowe PW (1969) The relation between the shear strength of sands in triaxial compression, plane strain and direct shear. Géotechnique 19(1): 75-86, http://dx.doi.org/10.1680/geot.1969.19.1.75.
Rowe PW (1972) The relevance of soil fabric to site investigation
practice. Géotechnique 22(2): 195-300, http://dx.doi.org/10.1680/geot. 1972.22.2.195.
Salgado R (2008) The Engineering of Foundations. McGraw-Hill, New York, NY, USA.
Schofield AN (1980) Cambridge Geotechnical Centrifuge operations. Géotechnique 30(3): 227-268, http://dx.doi.org/10.1680/geot.1980.30. 3.227.
Shiraishi S (1990) Variation in bearing capacity factors of dense sand assessed by model loading tests. Soils and Foundations 30(1): 17-26, http://doi.org/10.3208/sandf1972.30.17.
Siddiquee MSA, Tanaka T, Tatsuoka F, Tani K and Morimoto T (1999) Numerical simulation of bearing capacity characteristics of strip footing on sand. Soils and Foundations 39(4): 93-109, http://doi.org/ 10.3208/sandf.39.4_93.
Siddiquee MSA, Tatsuoka F, Tanaka T et al. (2001) Model tests and FEM simulation of some factors affecting the bearing capacity of a footing on sand. Soils and Foundations 41(2): 53-76, http://doi.org/10.3208/ sandf.41.2_53.
Tatsuoka F (2001) Impacts on geotechnical engineering of several recent findings from laboratory stress-strain tests on geomaterials. In The 2000 Burmister Lecture at Columbia University, Geotechnics for Roads, Rail Tracks and Earth Structures (Gomez Correia A and Brandle H (eds)). Balkema, Rotterdam, the Netherlands, pp. 69-140.
Tatsuoka F, Goto S and Sakamoto M (1986a) Effect of some factors on strength and deformation characteristics of sand at low pressures. Soils and Foundations 26(1): 105-114, http://doi.org/10.3208/sandf1972.26. 105.
Tatsuoka F, Sakamoto M, Kawamura T and Fukushima S (1986b) Strength and deformation characteristics of sand in plane strain compression at extremely low pressures. Soils and Foundations 26(1): 65-84, http:// doi.org/10.3208/sandf1972.26.65.
Tatsuoka F, Okahara M, Tanaka T et al. (1991) Progressive failure and particle size effect in bearing capacity of a footing on sand. In Proceedings of ASCE Geotechnical Engineering Congress, 1991. American Society of Civil Engineers, Reston, VA, USA, ASCE Geotechnical Special Publication no. 27, pp. 788-802.
Terzaghi K (1929) Effects of minor geologic details on the safety of dams. In Geology and Engineering for Dams and Reservoirs. American Institute of Mining and Metallurgical EngineersTechnical Publication 215, pp. 31-44. Reprinted in Terzaghi K (1960) From Theory to Practice in Soil Mechanics: Selections from the Writings of Karl Terzaghi. Wiley, New York, NY, USA, pp. 119-132.
Terzaghi K (1943) Theoretical Soil Mechanics. Wiley, New York, NY, USA.
Toyosawa Y, Itoh K, Kikkawa N, Yang JJ and Liu F (2013) Influence of model footing diameter embedded depth on particle size effect in centrifugal bearing capacity test. Soils and Foundations 53(2): 349-356, http://dx.doi.org/10.1016/j.sandf.2012.11.027.
Ueno K, Miura K and Maeda Y (1998) Prediction of ultimate bearing capacity of surface footings with regard to size effects. Soils and Foundations 38(3): 165-178, http://doi.org/10.3208/sandf.38.3_165.
Ueno K, Miura K, Kusakabe O and Nishimura M (2001) Reappraisal of size effect of bearing capacity from plastic solution. Journal of Geotechnical and Geoenvironmental Engineering 127(3): 275-281, http://dx.doi.org/10.1061/(ASCE)1090-0241(2001)127:3(275).
Valore C and Ziccarelli M (2009) The evolution of grain-size distribution of sands under 1-D compression. In Proceedings of the 17th International Conference on Soil Mechanics and Geotechnical Engineering (Hamza M, Shahien M and El-Mossallamy Y (eds)). IOS Press, Amsterdam, the Netherlands, vol. 1, pp. 84-88.
Valore C, Ziccarelli M and Muscolino SR (2017) The bearing capacity of footings on sand with a weak layer. Geotechnical Research, http://dx. doi.org/10.1680/jgere.16.00020.
Vesic AS (1973) Analysis of ultimate loads of shallow foundations. Journal of the Geotechnical Engineering Division 99(SM1): 45-73.
Vesic AS (1975) Bearing capacity of shallow foundations. In Foundations Yamaguchi H, Kimura T and Fuji-i N (1977) On the scale effect of footings
Engineering Handbook (Winterkorn HF and Fang HY (eds)). Van in dense sand. Proceedings of the 9th ICSMFE, Tokyo, Japan, vol. 2,
Nostrand Reinhold, New York, NY, USA, pp. 121-147. pp. 795-798.
White DJ, Teh KL, Leung CF and Chow YK (2008) A comparison of the Yin JH, Wang YJ and Selvadurai APS (2001) Influence of nonassociativity
bearing capacity of flat and conical circular foundations on sand. on the bearing capacity of a strip footing. Journal of Geotechnical and
Géotechnique 58(10): 781-792, http://dx.doi.org.10.1680/geot.8.D. Geoenvironmental Engineering 127(11): 985-989, http://dx.doi.org/10.
024. 1061/(ASCE)1090-0241(2001)127:11(985).
Yamaguchi H, Kimura T and Fuji-i N (1976) On the influence of Zhu F, Clark JI and Phillips R (2001) Scale effect of strip and circular
progressive failure on the bearing capacity of shallow foundations in footings resting on dense sand. Journal ofGeotechnical Engineering
dense sand. Soils and Foundations 16(4): 11-22, http://doi.org/10. 127(7): 613-621, http://dx.doi.org/10.1061/(ASCE)1090-0241(2001)
3208/sandf1972.16.4_11. 127:7(613).
How can you contribute?
To discuss this paper, please submit up to 500 words to the editor atjournals@ice.org.uk. Your contribution will be forwarded to the author(s) for a reply and, if considered appropriate by the editorial board, it will be published as discussion in a future issue of the journal.
