Al-Salakh 2019 IOP Conf. Ser. Mater. Sci. Eng. 584 012056
Al-Salakh 2019 IOP Conf. Ser. Mater. Sci. Eng. 584 012056
1. Introduction
Problem of the seismic bearing capacity of shallow foundations is far from being considered solved.
Many researchers stated that the response of a footing to dynamic loads is affected by:
(1) Nature and magnitude of dynamic loads,
(2) Number of pulses and
(3) The strain rate response of soil.
Shallow foundations for seismic loads are usually designed by the equivalent static approach. The
foundations are considered as eccentrically loaded and the ultimate bearing capacity is accordingly
estimated. To account for the effect of dynamic nature of the load, the bearing capacity factors are
determined by using dynamic angle of internal friction which is taken as 2-degrees less than its static
value [1].
International Building Code permits an increase of 33 % in allowable bearing capacity when earthquake
loads in addition to static loads are used in design of the foundation. This recommendation may be
reasonable for dense granular soils, stiff to very stiff clays or hard bedrocks but is not applicable for
friable rock, loose soils susceptible to liquefaction or pore water pressure increase [2].
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
horizontal and vertical accelerations to the center of gravity of the structure, the problem is reduced to
a static case of bearing capacity with inclined and eccentric loads. Nevertheless, in these methods the
inertial effects on the soil mass beneath the foundation are not considered and this could have important
effects on the overall system response.
Isolated column footings, strip footings, mat footings, and even pile foundations all may fail during
seismic events. Such failures are generally attributed to liquefaction. However, a number of failures
have occurred where field conditions indicate there was only partial saturation or a dense soil and
therefore liquefaction alone is a very unlikely explanation. Rather, the reasons for the seismic
settlements of these foundations seem to be that the bearing capacity was reduced [3]. Richards et al [3]
proposed a simplified approach to estimate the dynamic bearing capacity que and seismic settlement
SEq of a strip footing for assumed failure surfaces as in equation (1). There is hardly any experimental
verification of these theoretical solutions. There is a need for such validation [1].
V2 𝑘 −4
𝑆𝑒𝑞 (m) = 0.174 A.g ( 𝐴ℎ ) . 𝑡𝑎𝑛𝛼𝐴𝐸 (1)
Where: SEq is the settlement due to earthquake, V = peak velocity for the design earthquake (m/sec),
A=acceleration coefficient for the design earthquake, g = acceleration due to gravity (9.81 m/sec2), the
values of kh∗ and tanαAE can be obtained from some well-known figures which related the required
factor with angle of internal friction and bearing capacity as mentioned in [4] as well as in table 1 [4].
Generally, the behavior of any structure during earthquake depend on many items such as earthquake
magnitude, acceleration, duration, height of the structure/width of the foundation (aspect ratio),
foundation boundaries, soil properties and condition, and probability of liquefaction triggering [1].
In present study the specified structure has an aspect ratio < 0.8 according to this ratio and as cited in
many researches the structure will suffer from vertical settlement only without any visible tilting [5],
while in dry soil the foundation will suffer from excessive settlement and sliding in the wards of
shaking.
The propagation of the acceleration amplitude through soil depth will cause additional stress and
rearrangement to the soil particles to produce more dense soil, which causes a settlement and sliding
failure to the shallow structures and infrastructures.
Table 1. Variation of tanαAE with kh and soil friction angle φ (Compiled from [3])
In the settlement equation, it is clear that it neglects the effect of saturation or PWP generation, changes
in soil density due to shaking, mitigation, and changes in effective stress or bearing capacity failure. To
understand the effect of these factors, it is very important to understand the effect of each item on soil
behavior individually and together.
As well as the settlement formula [3] suggested a simplified approach to evaluate the dynamic bearing
capacity que and seismic settlement SEq of a strip footing for assumed failure surfaces. The proposed
seismic bearing capacity (quE) is given by equation (2) as follow:
2
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
NcE, NqE, and NγE are the factors of seismic bearing capacity.
These factors are related to the magnitude of dynamic angle of internal friction of the soil (φdyn) and
tanψ = kh / (1-kv)
kh and kv are the horizontal and vertical coefficients of acceleration due to earthquake.
For static case, kh = kv, = 0 and equation (1) becomes:
qu = cNc + qNq + ½ γBNγ (3)
3
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
The comparison between acceleration values at different levels in the shaking table test shows the effect
of foundation level on its behavior during shaking.
The acceleration at the foundation level will become larger than its value at the source due to the
following:
- Effect of inertia mass motion.
- Effect of loose soil strata behavior during shaking.
- Effect of generation of PWP and variation of stress states.
Figure 2 presents the variation of acceleration with the depth of soil. It is clear that when the
acceleration at hard strata is 0.2g the resulting surface acceleration at foundation level become more
than 0.35g.
When saturated soil (the water table is close to ground level) exposed to an earthquake, the pore water
pressure will increase. If pore water pressure continues to increase, liquefaction will take place in a few
seconds as shown in figure 3. This effect will result in a decrease in bearing capacity of the soil.
Therefore, it is proposed to include the effect of pore water pressure in the determination of the dynamic
bearing capacity of the soil.
Figure 1. Shaking table schematic diagram explain soil layers and distribution of transducers.
4
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
(a)
(b)
(c)
Figure 2. (a) Acceleration (g) Vs time (sec.) at bottom saturated sand layer for 1Hz.; (b) Acceleration
(g) vs time (sec.) at 1st saturated sand layer for 1Hz.; (c) Acceleration (g) Vs time (sec.) at foundation
for 1Hz.
5
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
Figure 3. Generation of PWP (bar) Vs time (sec.) at top layer of saturated soil model.
The more suitable and meaningful expression which relate the effects of PWP effective stress changes
is the pore water pressure ratio (ru), this factor can be introduced in the bearing capacity equation as a
control parameter for all cases as follow:
1
qu = (1 − 𝑟𝑢) ∗ (𝐶. 𝑁𝑐𝑒𝑞. +𝑞. 𝑁𝑞𝑒. . + 2 ᵞ𝐵. 𝑁ᵞ𝑒) (4)
When the soil is dry ru value become equal to zero thus there is no effect on the original equation, at r u
equal to 0.5 which mean that the resulting bearing capacity will reduce to the half magnitude due to the
generation of PWP and reduction in effective stress, and when ru become equal or near to 1 which mean
that the liquefaction has occurred in this case the bearing capacity will has zero value which compatible
with actual or field situation. When liquefaction mitigation methods applied on the liquefiable soils the
following changes and improvement recognized:
- The soil strength increases to resist the shaking.
- There are specific reductions in the generation of pore water pressure during shaking.
- There is an increase in soil density due to filling voids as shown in table 2.
These effects will lead to maintain the propagation of acceleration along the soil column during shaking
which leads to fixing the acceleration value from bottom to top layers.
- The above changes in mitigated soil produce a noticeable increase in soil stability and the bearing
capacity of the mitigated soils.
Table 2. Effect of mitigation method (mixing percent and dilution fraction) on soil unit weight.
Sat. Sand + Bentonite slurry grout unit weight, Sand + CKD slurry grout unit weight,kN/m3 Granular
sand kN/m3 pile
unit 3% 5% 10% 3% 5% 10% equivalent
weight unit weight
kN/m3 kN/m3
3/1 6/1 3/1 6/1 3/1 6/1 3/1 6/1 3/1 6/1 3/1 6/1
18.5 19.3 18.92 19.5 19.5 19.24 19.2 19.3 19.4 19.25 19 18.65 19 18.65
The presence of water within soil strata will make the soil settle due to lubrication of the particles which
tend to become closer to each other. This phenomenon recognized in the experiment when the saturated
soil left for one day leading to settlement. When the shaking or earthquakes are applied there are
noticeable generations of PWP which cause soil liquefaction.
The change in soil unit weight due to shaking in the dry state is differing from a change in soil density
due to liquefaction mitigation process. In the first case, the densification causes a specific settlement
due to soil particles rearrangement while in soil mitigation the change in density is due to replacement
of the water inside the voids by grouting materials without and rearrangement or densification of soil
particles. This behavior exactly recognized during shaking tests at all applied acceleration. The main
difference between the non-mitigated soil and the mitigated one is the magnitude of the settlement.
6
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
It is important to note that the simplified approach proposed by [3] is derived to estimate settlement for
non-liquefied soils. This may explain the divergence between experimental and the simplified approach,
since the settlement measured experimentally for liquefied soil post liquefaction.
Also, at mitigated soil the estimation of settlement according to the same formula required a
modification so as to adjust the settlement magnitude.
Settlement formula didn’t take into account the soil density, degree of saturation, generation of PWP,
and reduction in effective stress during shaking. A modified formula for settlement estimation
suggested.
This modification makes Richards formula takes the effect of soil density and pore water pressure
generation during an earthquake in consideration, See equations (5) and (6).
For frequency amplitude ≤ 0.75 Hz.
𝛾𝑠𝑜𝑖𝑙 𝑎 𝑏
𝑆𝐸𝑞𝑚 = 𝑆𝐸𝑞 ∗ ( ) ∗( ) ( ru ≤0.90) (5)
𝛾𝑛𝑒𝑤 (𝑐−𝑟𝑢 𝑑 )
a, b, c, and d: are an empirical coefficient from statistical analysis by SPSS software and it equal to (
58.8, -1.85, 1, and -0.133 respectively) the regression coefficient (R2 = 0.95).
For frequency amplitude > 0.75 Hz.
The empirical coefficients become as follow: (-35.1, -321.6, -468.8 and -3.83 respectively) the
regression coefficient (R2 = 0.86).
Where ᵞnew is the new or improved soil unit weight after mitigation or improvement process. In case of
ru (0.9-1) there are additional formula can be introduce which cover the effect of pore water pressure
generation and the reduction in effective stress. This formula can be written as follow:
𝛾𝑠𝑜𝑖𝑙 𝑒 1 𝑔
𝑆𝐸𝑞𝑚 = 𝑆𝐸𝑞 ∗ ( ) ∗( ) (6)
𝛾𝑛𝑒𝑤 𝑟𝑢
Where e, is an empirical coefficient from statistical analysis which equal to 1.5 and g = (γsoil / (1- (γnew
- γsoil)) /100, the regression coefficient (R2 = 0.89).
The above-modified formula has been obtained using statistical analysis by using SPSS software. A
good convergence between modified formula and experimental measurement noticed in table 3 and
Figures 4, 5, 6 and 7.
The settlement of the foundation during shaking before and after soil mitigation were recorded for some
cases and calculated using original and modified formula.
It is observed that the original formula is suitable for high acceleration levels without probability of
liquefaction triggering.
Table 3. Settlements of saturated pure, mitigated soils at different acceleration levels depending [3],
modified, and experimental measurements.
Settlement formula Saturated sand Sand +5% Sand +CKD +6/1 Granular piles
Bentonite+6/1
0.05g 0.1g 0.2g 0.05g 0.1g 0.2g 0.05g 0.1g 0.2g 0.05g 0.1g 0.2g
Formula of [3] 0.05 0.83 11.6 0.05 0.83 11.6 0.05 0.83 11.6 0.05 0.83 11.6
Modified Ritchard 2.0 17.6 25.3 0.064 1.56 22.2 0.065 1.03 11.6 0.05 1.4 15.15
Experimental 2.5 19.1 26.0 1.2 2.8 31.0 0.1 0.75 10.2 0.3 3.5 15.35
results
7
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
Figure 4. Relationship between settlement of foundation and acceleration on saturated pure sand.
8
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
Figure 7. Relationship between settlement of foundation and acceleration on saturated sand+ granular
pile.
4. Conclusions
A shallow footing resting on liquefiable poorly graded saturated sand was subjected to horizontal
acceleration by a manufactured controlled and instrumented shaking table. The acceleration of the
foundation and of the sand deposit is quite different from the input table motion. A modification on
bearing capacity equation was introduced for seismic design taking into account the buildup in pore
water pressure. This modification makes the equation applicable on liquefiable soil. Another
modification was prepared for [3] equation of settlement to take into account the liquefaction triggering
due to the buildup in pore water pressure. Good convergence was observed between the predicted and
measured settlement.
References
[1] Vijay K. Puri and Shamsher Prakash 2007 ON Foundations under Seismic Loads 4th
International Conference on Earthquake Geotechnical Engineering June 25-28, 2007, Paper No.
1118.
[2] Robert W. Day 2006 Foundation Engineering Handbook I edition, McGraw-Hill Companies
Inc.
[3] R. Richards Jr., D. G. Elms, and M. Budhu 1993 Seismic Bearing Capacity and Settlements of
Foundations (Journal of geotechnical engineering 1993 Vol.119) pp. 627-662.
[4] Braja M. Das and G. V. Ramana 2011 Principles of Soil Dynamics Second Edition (Stamford,
CT 06902, Cengage Learning press.).
[5] Gazetas,G. M. Apostou and J.Anasta- sopoular 2004 Seismic Bearing Capacity Failure and
Overtiming of Terveler Building in Adapazari 1999 Proc.5th. International conference on case
histories in Geo. Engg. (New York) SOAP11(1-51).
[6] Gopal S. P. Madabhushi and Stuart K. Haigh 2010 Liquefaction Induced Settlement of
Structures Fifth International Conference on Recent Advances in Geotechnical Earthquake
Engineering and Soil Dynamics, 2010.
[7] H.H. Schwab 1981 Settlement and Pore Pressure Due to Cyclic Loading (Technical University
of Darmstadt, Federal Republic of Germany).
[8] Indrajit Chowdhury and Shambhu P. Dasgupta 2016 Dynamic Bearing Capacity of Shallow
Foundation Under Earthquake Force (Indian Geotechnical Jouranal, 2016).
[9] Jonathan D. Bray and Jorge Macedo 2017 Simplified Procedure for Estimating Liquefaction-
Induced Building Settlement Proceedings of the 19th International Conference on Soil
Mechanics and Geotechnical Engineering, Seoul 2017 pp 101-118.
[10] Juan Carlos Tiznado A. and Danilo Paillao V. 2014 Analysis of the Seismic Bearing Capacity
of Shallow Foundations (Journal of Construction, 2014, Vol. 13(2)), pp. 40 – 48.
9
International Conference on Civil and Environmental Engineering Technologies IOP Publishing
IOP Conf. Series: Materials Science and Engineering 584 (2019) 012056 doi:10.1088/1757-899X/584/1/012056
10