Engineering,Environment
Evaluation ofbearing capacity and settlement of foundations
Bunyamin Anigilaje SALAHUDEEN 1* and Ja'afarAbubakar SADEEQ 2
1 Samaru College ofAgriculture, Ahmadu Bello University, Zaria, Nigeria
2 Department ofCivil Engineering, Ahmadu Bello University, Zaria, Nigeria
Emails: * basalahudeen@gmail.com;jaafaras@yahoo.com
*Corresponding Author phone:+2348058565650
Received: February 12, 2016/ Accepted: December 8, 2016 /Published: December 30, 2016
Abstract
The standard penetration test (SPT) results (SPT N-values)used in this study were corrected to the standard average energy of 60% (N60)before using them to correlate soil properties and evaluate foundation settlementcharacteristics and bearing capacity in the North Central zone of Nigeria. Basedon the corrected N-values, some geotechnical design parameters including the allowablebearing pressure and elastic settlement of foundations were predicted atvarying applied foundation pressures of 50, 100, 200, 300 and 500 kN/m2using conventional analytical models and numerical modelling. The numericalanalysis results using Plaxis 2D, a finite element code, shows that Meyerhof’sand Peck’s et al. analytical/empirical methods of estimating the allowablebearing pressure of shallow foundations provide acceptable results. Resultsobtained show that an average bearing capacity value of 150– 350 kN/m2 can be used for shallowfoundations at embedment depth of 0.6 to 3.6 m in the North Central zone. Basedon recommendation of Eurocode 7 which allows a maximum total settlement of 25mm for serviceability limit state, it is recommended that raft or deepfoundations to be considered for applied foundation pressures exceeding 300 kN/m2in the North-Central zone to avoid excessive settlement.
Keywords
Numerical modelling; Plaxis 2D; Finite element method; StandardPenetration Test; Bearing capacity; Elastic settlement
Introduction
One of the mostsignificant components of any structure is its foundation. Foundations areintegral to overall structural performance. They help in bearing andtransmitting the structural loads to the soil, reducing settlements (total anddifferential), preventing possible movement of structures due to periodicshrinkage and swelling of subsoils, allow building over water or water-loggedgrounds, resist uplifting or overturning forces due to wind, and resist lateralforces due to soil movement and control water penetration and dampness. Toperform satisfactorily, foundations must have two main characteristics: theyhave to be safe against overall shear failure in the soil that supports themand they must not undergo excessive settlement [1].
Probably the most difficult of the problemsthat a soil engineer is asked to solve is the accurate prediction of thesettlement of a loaded foundation [2]. The problem is in two distinct parts:the value of the total settlement that will occur and the rate at which thisvalue will be achieved. The design of shallow foundations is generallycontrolled by settlement rather than bearing capacity [3]. As a consequence,settlement prediction is a major concern and is an essential criterion in thedesign process of foundations. Consistent and accurate prediction of settlementis yet to be achieved by the use of a variety of analytical methods [4].
Comparative studies of the availablemethods by engineers/researchers [5-8] indicate inconsistent prediction of themagnitude of the calculated settlements. This may be attributed to the factthat most of these methods are model driven, in which the form of the model hasto be selected in advance and the unknown model parameters are determined byminimising an error function between model predictions and known historicalvalues. Consequently, prior knowledge regarding the relationship between modelinputs and corresponding outputs is needed. In case of settlement of shallowfoundations on granular soils, such knowledge is not yet entirely understood [3].
The finite elementmethod can be particularly useful for identifying the patterns of deformationsand stress distribution during deformation and at the ultimate state. Becauseof these capabilities of the finite element method, it is possible to model theconstruction method and investigate the behaviour of shallow footings and thesurrounding soil throughout the construction process, not just at the limitequilibrium conditions [9]. The finite element method (FEM) allows modellingcomplicated non linear soil behaviour through a constitutive model, variousgeometrics with different boundary conditions and interfaces. It can predictthe stresses, deformations and pore pressures of a specified soil profile [10].
The high level of demandsfor housing in Nigeria due tothe growing population and migration of people to urban areas requirealternative construction methods that provide fast, safe and affordable qualityhousing for the citizens [11-13]. The aim of this study is to explore numericalmodelling method that better represents soil constitutive behaviour to developan improved approximation of foundation soil bearing capacity and settlementand compare the results with those of empirical methods. Also, there is a needto investigate and determine the most appropriate methods to Nigerian soilpeculiarities and distinctions based on SPT results, being the most common andeconomical geotechnical field test used in Nigeria.The study focused on the prediction of foundation soil bearing capacity andsettlement based on Standard Penetration Test (SPT) N-values using empirical/analytical(deterministic) models and numerical modelling in the North-Central zone (i.e. Benue, Federal Capital Territory, Kogi, Kwara, Nassarawa, Nigerand Plateau States)of the Federal Republic of Nigeria.
The objectives of thisresearch were to estimate the bearing capacity and settlement of foundation soilsfrom measured penetration resistance in terms of the SPT corrected N-values atvarying depths, to evaluate design equations for foundation soil bearingcapacity and settlements using different constitutive models based on SPTresults, to model foundation settlement numerically using PLAXIS 2D software andto compare the results of empirical/analytical methods with those of numericalanalysis.
Research Methodology
The research made use ofstandard penetration test (SPT) data (using Donut hammer type) collected from 592test holes (5328 data set) distributed over the study area. Computations weredone based on the average that reliably represents each State and the averageof the States was used for the North-Central zone. Bearing capacity andfoundation settlement estimations were made at depths of 0.6, 2.1, 3.6, 5.1,6.6, 8.1, 9.6, 11.1 and 12.6 m and settlement was determined at varying appliedfoundation pressures of 50, 100, 200, 300 and 500 kN/m2.
Based on the analyticalmethods, bearing capacity and foundation settlement were evaluated using somecarefully selected models listed in Table A1 and A2, respectively (see Appendix).On the other hand, numerical analysis of foundation settlement and bearingcapacity were performed using a non-linear finite element analysis with a finiteelement code, Plaxis 2D, which is a software for deformation analysis andmodelling of geotechnical problems.
The input data in Plaxisare index, elastic and strength parameters obtained from the processed SPTN-values. The software portfolio includes simulation of soil and soil-structureinteraction. Plaxis 2D is an axisymmetric finite element package used fortwo-dimensional analysis of deformation and stability in geotechnicalengineering. It uses advanced soil constitutive models for the simulation ofthe non-linear, time dependent and anisotropic behaviour of soils and rocks.Plaxis 2D portfolio models the structure, the soil and the interaction betweenthe structure and the soil. Soil layers and foundation structure parameters areinputted into Plaxis and the construction stages, loads and boundary conditionsare defined in an already defined geometry cross-section containing the soilmodel then Plaxis automatically generates the unstructured 2D finite elementmeshes with options of global and local mesh refinements. Using its calculationfacilities, Plaxis 2D undergoes a calculation process and presents thecalculation and model outputs which can be accessed in animation and/ornumerical forms. The input data in numerical modelling are index, elastic andstrength parameters obtained from the processed SPT N-values [14]. The steps involved in developing the numerical model can bedepicted by the chart shown in Figure 1.
Figure 1. Chartdepicting the steps involved in developing the numerical models
Results and Discussion
The results of thisstudy were computed using the most common and conventional empirical methods inthe literature based on the average values of input data obtained from theavailable field information at the time of the study. It is pertinent to statethat result presented herein can approximately represent soil conditions in theNorth-Central zone of Nigeria considered. The elastic settlement and allowablebearing capacity results of the empirical/analytical methods were compared withthose of numerical modelling output using Plaxis.
Corrected N-values (N60)
According to Bezgin [15]a correction to average energy ratio of 60% (N60) is required to SPTN-values because of the greater confinement caused by the increasing overburdenpressure. The correction factors used in the study are those proposed by Das [1]to standardize the field penetration number as a function of the input drivingenergy and its dissipation around the sampler into the surrounding soil. Thevariation of N60 with depth of test is shown in Figure 2.
Figure 2. Variation of corrected N-values withboring depth
N60 increasedwith depth having the highest value of 89.25 on the average in the North-Centralzone. The increase in N60 value with depth is due to increasedoverburden. This confirms the conclusions of Salahudeen et al. [13] that thesoils in the northern part of Nigeriaare crystalline in nature (which has higher N60 values compared withsedimentary formations) from the basement complex. N60 values areneeded for more accurate design analyses and have less variability or scatterdue to the test method.
Bearing capacity
Based on field testresults, the bearing capacities of shallow foundations are determined in termsof the allowable bearing pressures while those of deep foundations (piles) aregiven in terms of the ultimate bearing capacity. This isbecause settlement (service limit) controls the allowable bearing capacity indesign of shallow foundations while the ultimate limit (shear failure) usuallycontrols the allowable bearing capacity in deep foundations design [13].For the allowable bearing pressures of shallow foundations, footing plandimensions of 2 m by 2 m by 0.4 m for length, breadth and depth, were respectivelyassumed with a safety factor of 3. Variations of allowable bearing capacity ofshallow foundations and bearing capacity of piles with boring depth are shownin Figures 3 and 4, respectively.
Figure 3. Variation of allowable bearingpressure with depth
Figure 4. Variation of bearing capacity of pilewith depth
Based on the methodproposed by Meyerhof [16] and Plaxis, foundation pressures in the range of 150– 350 kN/m2 were proposed for use in the North-Centralzone at shallow depths (depths in the range of 0.6 - 3.6 m).
In evaluating thebearing capacity of piles, assumed dimensions of 0.3 m by 0.3 m cross-sectionalarea with embedment length of 10 m were used. Although the results according toBriaud et al. [21] are high compared to those proposed by Meyerhof [22],based on 60 pile case histories and SPT borehole results comprising 43 fullscale pile load tests and 17 dynamic tests with Case Pile Wave Analysis Program(CAPWAP) collected from 18 sources reporting data from 26 sites and from 7countries, Shariatmadari et al. [23] reported up to 10,505 kN/m2pile bearing capacity at 25 m depth. Maximum pile bearing capacity values of3748.50 and 10006.36 kN/m2were obtained for the North-Central zone using Meyerhof [22] and Briaud etal. [21] methods respectively, at foundation embedment depth of 12.6 m.
Elastic settlement of foundations
For the elasticsettlement of shallow foundations, plan dimensions of 2 m by 2 m by 0.4 m forlength, breadth and depth were respectively assumed. Variations of elasticsettlement of foundations with embedment depth for various applied pressuresare shown in Figures 5 - 9.
Figure 5. Variation of elastic settlement withdepth for 50 kN/m2 foundation pressure
Figure 6. Variation of elastic settlement withdepth for 100 kN/m2 foundation pressure
Figure 7. Variation of elastic settlement withdepth for 200 kN/m2 foundation pressure
Figure 8. Variation of elastic settlement withdepth for 300 kN/m2 foundation pressure
Figure 9. Variation of elastic settlement withdepth for 500 kN/m2 foundation pressure
The figures show thedifferent empirical/analytical models commonly used in computing elasticsettlement of shallow foundations. The N60 values indicate thatsettlement values will be low due to high N60 values in the regionas confirmed in the elastic settlement results. The recorded trend isconsistent with observations reported by Rasin [24] and Salahudeen et al.[13].
The numerical analysisresults of soil body deformation, stress distribution and settlementrespectively at collapse of the soil body for the North-Central zone at 0.6 and12.6 m depths of embedment are shown in Figures 10 - 15.
Figure 10. Numerical analysis mesh showingdeformation of the soil body at collapse at 0.6 m embedment depth in the North-Centralzone
Figure 11. Numerical analysis result of stressdistribution up to the collapse of the soil body at 0.6 m embedment depth inthe North-Central zone
Figure 12. Numerical analysis result ofsettlement up to the collapse of the soil body at 0.6 m embedment depth in the North-Centralzone
Figure 13. Numerical analysis mesh showingdeformation of the soil body at collapse at 12.6 m embedment depth in the North-Centralzone
Figure 14. Numerical analysis result of stressdistribution up to the collapse of the soil body at 12.6 m embedment depth inthe North-Central zone
Figure 15. Numerical analysis result of settlementup to the collapse of the soil body at 12.6 m embedment depth in the North-Centralzone
A comparison among theanalytical methods used in this study with the results of numerical modellingshow that the methods proposed by Schmertmann et al. [29], Burland andBurbidge [31], Canadian Foundation Engineering Manual (CFEM) [34] as well asthat of Mayne and Poulos [36] gave good estimations of foundation settlement.
Total settlement of piles
Variation of settlementof pile (bored) under a vertical working load with depth based on methodsproposed by Vesic [38] and that of Das [1] is shown in Figure 16.
Figure 16. Variation of total settlement ofpiles with depth
The total settlement isthe sum of elastic settlement of pile, settlement of pile caused by the load atthe pile tip and settlement of pile caused by the load transmitted along thepile shaft. The wide margin between the methods used is an indication that thesettlements need to be modelled in order to come up with an appropriate methodthat will be more suitable for Nigerian soils. However, it is pertinent tostate that numerical modelling of piles was not included in the scope of thestudy carried out.
Conclusions
The study consideredN-values corrected to the standard average energy of 60% (N60) asinput data in analytical/empirical and numerical models used to predictfoundation settlement and bearing capacity in the North-Central zone of Nigeriafor footing of 2 m by 2 m by 0.4 m size and varying pressures of 50, 100, 200,300 and 500 kN/m2. Based on the results of this study, the followingconclusions can be taken:
¸Allowable bearing pressures of 150-350kN/m2 at depths between 0.6 and 3.6 m obtained using the Meyerhof methodare adequate for North-Central soils. The values are very close with those ofnumerical analysis using Plaxis 2D.
¸The maximum elastic settlement valuesrespectively for all the applied foundation pressures show that the soils inthe North-Central zone are on the average, less susceptibility to compression.
¸It was observed that settlement increased withincreased value of applied foundation pressure. Settlements of footingsembedded at depths in the range 0.6-3.6 m and pressures above 300 kN/m2exceeded the limiting value of 25 mm value of allowable total settlement stipulatedby Eurocode 7.
¸A comparison of the fifteen empirical/analyticalmethods considered in this study, showed that the Schmertmann et al.[29], Burland and Burbidge [31], Canadian Foundation Engineering Manual (CFEM) [34]as well as the Mayne and Poulos [36] methods gave good estimations offoundation settlement.
¸Plaxis 2D tend to overestimate the elasticsettlement of footings for embedment depths up to 2 m and at applied foundationpressure greater than 300 kN/m2.
Recommendations
Based on the results ofthis study, the following are hereby recommended for the North-Central zone ofNigeria:
¸Foundations should be placed at a minimum depthof 1.0 m.
¸Deep foundations should be considered forapplied foundation loads exceeding 300 kN/m2 to avoid excessivesettlement.
¸Results of the study can be used as firstapproximation of foundation bearing capacity and settlement but does notpreclude the use of site specific data.
Acknowledgements
The authors wish toacknowledge the assistance of the Management of In-depth Engineering Limited, Kaduna, Nigeria that provided all standardpenetration test data used for the study. For the assistance of Dr. M. Jaliliof Islamic Azad University, Semnan, Iran with respect to training on the use ofPlaxis software.
Appendix
Table A1. Empirical/analyticalmodels for soil bearing capacity analysis
Property | Model | Reference |
Corrected N-value | | Seed et al. (1985) and Skempton (1986) |
Allowable bearing pressure of shallow foundations | | Teng (1969) |
| Meyerhof (1974) | |
| Peck etal. (1974) | |
| Bowles (1996) | |
| Terzaghi et al. (1996) | |
Bearing capacity of piles | | Meyerhof (1976) |
| Briaud et al. (1985) | |
Total settlement of piles | | *Vesic (1977) **Das (2011) |
Table A2. Empirical/analyticalmodels for elastic settlement analysis
S/N | Model | Reference |
1 | | Janbu et al. (1956) |
2 | | Terzaghi and Peck (1967) |
3 | | Schmertmann (1970) |
4 | | Schultze and Sherif (1973) |
5 | | Meyerhof (1974) |
6 | | Schmertmann et al. (1978) |
7 | | Timoshenko and Goodier (1982) |
8 | | Burland and Burbidge (1985) |
9 | | Bowles (1987) |
10 | | Anagnostropolous et al. (1991) |
11 | | Can. Found. Eng. Man. (CFEM) (1992) |
12 | | Papadopoulos (1992) |
13 | | Terzaghi et al. (1996) |
14 | | Mayne and Poulos (1999) |
15 | | Anderson et al. (2007) |
Abbreviations
(N1)60= N60 correction for overburden pressure
B =Width of foundation (m)
BR= Reference width = 0.3 m
Df= Depth of embedment (m)
Es =Appropriate value of elastic modulus of soil (kN/m2)
Es= Elastic modulus of soil
H =Thickness of the compressible layer (m)
L =Length of foundation (m)
N =Measured penetration number (N-value)
N60= Corrected standard penetration number for field conditions
N60(a)= Adjusted N60 value
Pa= Atmospheric pressure = 100 kN/m2
q =Applied foundation pressure (kN/m2)
q = Neteffective pressure applied at the level of the foundation (kN/m2)
qn=Net pressure on the foundation (kN/m2)
Se= Elastic settlement (mm)
Se(1)= Elastic settlement of piles
Se(2)= Settlement of pile caused by the load at the pile tip
Se(3)= Settlement of pile caused by the load transmitted along the pile shaft
μS= Poisson’s ratio of soil
σ10= Effective overburden pressure in kN/m2
References
1.Das, B. M., Principles of FoundationEngineering, SI, Seventh Edition, Cengage Learning. USA,2011.
2.Klemencic R., McFarlane I. S., Hawkins N. M., NikolaouS., Seismic design of reinforced concrete mat foundations, NEHRP SeismicDesign Technical Brief, NIST GCR 12-917-22, 7, 2012.
3.Shahin M. A., Maier H. R., Jaksa M. B., Predictingsettlement of shallow foundations using Neural Networks, Journal ofGeotechnical and Geoenvironmental Engineering, ASCE, 2002, 128 (9), p. 785-793.
4.Poulos H. G., Common procedures forfoundation settlement analysis-Are they adequate?, Proc. 8th Australia NewZealand Conf. on Geomechanics, Hobart, 1999, p. 3-25.
5.Kalhor M., Azadi M., Assessment of changes inelasticity modulus and soil Poisson’s Ratio on applied forces to tunnel liningduring seismic loading, International Research Journal of Applied and BasicSciences, 2013, 4 (9), p. 2756-2761.
6.Hussein H. M. A., Effects of FlexuralRigidity and Soil Modulus on the Linear Static Analysis of Raft Foundations,Journal of Babylon University, Pure and Applied Science, 2011, 19(2), p.740-752.
7.Hooshmand A., Aminfar M. H., Asghari E., AhmadiH., Mechanical and Physical Characterization of TabrizMarls, IranSpringer Science+Business Media B.V. Geotech. Geol. Eng.,2011, 30, p. 219-232.
8.Banadaki A. D., Ahmad K., Ali N., Initialsettlement of mat foundation on group of cement columns in peat: numericalanalysis, Electronic Journal of Geotechnical Engineering (EJGE), 2012, 17,p. 2243-2253.
9.Laman M., Yildiz A., Numerical studies ofring foundations on geogrid-reinforced sand, Geosynthetics International,2007, 14 (2), p. 1–13.
10.Subramaniam P., Reliability based analysis ofslope, foundation and retaining wall using finite element method, PublishedMaster of Technology Thesis, Department of Civil Engineering, NationalInstitute of Technology, Rourkela – 769008, India, 2011.
11.Osinubi K. J., A method for estimating settlementof weak soils in reclaimed bases, The Nigerian Engineer, 1992, 27 (4), p.41-50.
12.Osinubi K. J., Foundation options forbuildings erected on hydraulic sand fills, Proc. 11th RegionalConference on Soil Mechanics and Foundation Engineering, Cairo,11 - 15 December, 1995, 2, p. 85-102.
13.Salahudeen A. B., Ijimdiya T. S., Eberemu A. O.,Osinubi K. J., Prediction of bearing capacity and settlement of foundationsusing standard penetration data in the South-South geo-political zone ofNigeria, Book of Proceedings, International Conference on ConstructionSummit, Nigerian Building and Road Research Institute (NBRRI) May 2016.
14.PLAXIS 2D manual, Plaxis 2D-Version 8 edited byBrinkgreve R. B. J. Delft University of Technology and PLAXIS b.v., TheNetherland, 2012.
15.Bezgin O., An insight into the theoreticalbackground of: Soil structure interaction analysis of deep foundations, Atechnical report, Istanbul, 2010.
16.Meyerhof G. G., Penetration Testing Outside Europe:General Report, Proceedings of the European Symposium onPenetration Testing. Available from National Swedish Institute for BuildingResearch, 785, S-801-29-GAVLEÄ, Sweden,1974, 2.1, p 40-48.
17.Teng W. C., Foundation Design,Prentice-Hall Inc., New Jersey, 1969.
18.Peck R. B., Hanson W. E., Thornburn T. H., FoundationEngineering, 2nd edition, Wiley, New York,1974.
19.Bowles, J. E., Foundation analysis and design,5th Edition, McGraw-Hill Book Co., New York, USA, 1996.
20.Terzaghi K., Peck R.B., Mesri G., SoilMechanics in Engineering Practice, (Third edition) John Wiley & Sons, New York, 549, 1996.
21.Briaud J. L., Tucker L., Lytton R. L., Coyle H.M., Behaviour of Piles and Pile Groups, Report No. Federal Highway Administration, Washington, DC, 1985.
22.Meyerhof G. G., Bearing Capacity andSettlement of Pile Foundations, The Eleventh Terzaghi Lecture, Journal ofGeotechnical Engineering Division, ASCE, 1976, 102 (GT3), p. 195-228.
23.Shariatmadari N., Eslami A., Karimpour-Fard M., Bearingcapacity of driven piles in sands from SPT applied to 60 case histories,Iranian Journal of Science & Technology, Transaction B, Engineering, 2008, 32(B2), p 125-140.
24.Rasin D., Observed and predicted settlementof shallow foundation, 2nd International Conference on New Developments in SoilMechanics and Geotechnical Engineering, Near East University,Nicosia, North Cyprus. 2009, p. 590-597.
25.Janbu N., Bjerrum L., Kjaernsli, B., VeiledningvedLosningavFundamentering-soppgaver,Publication 16, Norwegian Geotechnical Institute, Oslo, 1956, p. 30-32.
26.Terzaghi K., Peck R. B., Soil Mechanics in EngineeringPractice, 2nd Ed. John Wiley & Sons, New York, 1967.
27.Schmertmann J. H., Static cone to computestatic settlement over sand., Journal of Soil Mechanics and FoundationsDivision ASCE 96(SM3), 1970, p. 7302–1043.
28.Schultze E., Sherif G., Prediction ofsettlements from evaluated settlement observations for sand, In Proc., 8thInt. Conf. On Soil Mech. & Found. Engrg., 1973, 1(3), p.225-230.
29.Schmertmann J. H., Hartman J.P., Brown P. R., Improvedstrain influence factor diagrams, Journal of Geotechnical Engineering,Division, ASCE, 1978, 104 (8), 1131-1135.
30.Timoshenko, S. P., Goodier, J. N., Theory ofelasticity, Third edition, pp. 398-409, New York, McGraw-Hill, USA, 1982.
31.Burland J. B., Burbidge M. C., Settlement offoundations on sand and gravel, Proceedings of Institution of CivilEngineers, 1985, 78 (part. 1), p. 1325-1381.
32.Bowles J. E., Elastic foundation settlementon sand deposits, Journal of Geotechnical Engineering, ASCE, 1987, 113(8),p. 846-860.
33.Anagnostopoulos A.G., Papadopoulos, B.P., KavvadasM. J., SPT and compressibility of cohesionlesssoils, Proceedings of the2nd European Symposium on Penetration Testing, Amsterdam, 1991.
34.Canadian Foundation Engineering Manual, Thirdedition, BiTech, Publishers Ltd. Richmond, Canada, 1992.
35.Papadopoulos B. P., Settlements of shallowfoundations on cohesionless soils, J. Geotech. Eng.,ASCE, 1992, 118(3), p. 377-393.
36.Mayne P. W., Poulos H. G., Approximatedisplacement influence factors for elastic shallow foundations, Journal ofGeotechnical and Geo-environmental Engineering, ASCE, 1999, 125 (6), p. 453-460.
37.Anderson B.J., Townsend F. C., Rahelison L., Loadtesting and settlement prediction of shallow foundation, Journal ofGeotechnical and Geoenvironmental Engineering, ASCE. 2007, 133(12), p.1494-1502.
38.Vesic A. S., Design of pile foundations,National Cooperative Highway Research Program Synthesis of Practice,Transportation Research Board, Washington, DC, 1977, 42.