Introduction

With the promotion of the new era of Western Development and the strategy of “Build China’s Strength in Transportation”, the construction of railways, highways, airports, and other types of projects is gradually expanding from the plains and hills to the highlands and mountainous areas, especially in the high mountainous and extremely high mountainous areas of the Qinghai-Tibet Plateau and the Yunnan-Guizhou Plateau [1-3]. This type of area in the Quaternary Pleistocene period, through several glacial and interglacial period alternations, the end of the interglacial period glacial meltwater amount decreased, carrying capacity is weakened, the formation of fine-grained clay soil layer, and the accumulation of gravelly soils, block gravel soils during the water-abundant period in the form of layers of alternating distribution, by the influence of the dry and warm river valley climate, this type of area is low in rainfall, the soil body is dry, and in a long-term state of non-saturation [4, 5]. When these unsaturated ice-water- deposited clay soils are used as foundation soils, often cause foundation settlement, slope instability, and other engineering problems due to an insufficient foundation-bearing capacity. The Lijiang to Shangri-La Railway, which passes through the snow-covered plateau region in the middle of the Heng-Duan Mountains at the south-east edge of the Qinghai-Tibet Plateau, is built on the upper Pleistocene hard-plasticized ice-water deposited clay soils (Q3fgl) along the line, the unsaturated bearing capacity was not clearly explained during the selection of the parameters at the design stage. In the line design stage, the selection of the roadbed section and the foundation treatment program selection were affected [6]. Leshan-Ya’an Expressway along the distribution of complex engineering properties of ice and water deposition accumulation of soil, accumulation of soil in the lower part of the clay, due to the majority of sections belong to the high-fill roadbed and fill material for local materials, as the filler of the ice and water deposition of clay soil in a long-term non-saturated state, by the distribution of pore space, non-saturated soil in the pore microscopic connectivity characteristics of the different (Fig. 1), resulting in the difference between the pore water, air transport characteristics, and thus causing its bearing capacity is seriously insufficient. This results in a serious lack of bearing capacity, which has an important impact on the settlement of the roadbed and the stability of the high-fill roadbed slopes [7]. KangDing Airport in GanZi, Sichuan Province, is a high plateau-airport. During the investigation and design stage, because the foundation soil of the airport track area contains part of ice-water cohesive clay soil and is unsaturated state, due to the lack of understanding of the drainage characteristics and the bearing capacity parameters of the ice-water deposited cohesive soils when it is unsaturated state at that time, it resulted in several remedial repairs to the foundation during the operation and maintenance stage [8]. Therefore, an accurate interpretation of the bearing capacity of unsaturated ice-water deposited clay soils has important theoretical value for the construction operation and maintenance of actual projects. At the same time, according to the engineering practice, these three methods are often selected to carry out the planning and design of airport engineering in the actual project, but there is still a lack of relevant comparative studies on the suitability of these three methods. Therefore, this paper chooses these three methods to compare, in order to find a reasonable research method for the specific soil research. This research is of great significance to the design and construction of airport engineering, especially the high-altitude airport engineering. The comparative study of the three methods is also very reasonable for practical engineering.

Material and Methods

Related Theories

In all kinds of projects, the bearing capacity and settlement deformation of foundation soil involves the bearing capacity of the foundation soil, which depends on the strength of the soil itself. The fragmentation and multiphase nature of the soil body and the variability formed in the long-term geological history make its strength characteristics show a special nature different from that of other materials [9]. The damage to the soil body is mainly shear damage, and what determines its shear strength is often the cohesion and friction between soil particles. When unsaturated soils are used as foundation soils for various types of projects, matrix suction is generated at the shrinkage film inside the unsaturated soil unit under external loading. The presence of matrix suction makes the bearing characteristics of unsaturated foundation soils quite different from those of conventional saturated soils. A large number of studies have shown that in unsaturated soils, due to the presence of pore gas, which in turn generates the corresponding matrix suction, additional suction friction is generated in the soil, resulting in an unusually complex analysis of the strength characteristics of unsaturated soils [10].

The research in the area of strength characterization of unsaturated soils taking into account matrix suction began in the 1950’s. In 1959, Bishop proposed an effective stress formula for unsaturated soils,

(1)
σ′=(σ−ua)+χs

Where, σʹ is the effective stress of the unsaturated soil, σ is the total stress, ua is the pore gas pressure, s is the matrix suction (s = ua–uw), uw is the pore water pressure, and x is the parameter related to the degree of saturation.

In 1961, the effective stress formula was proposed at the Colloquium on Pore Pressure and Suction in Soils held in England:

(2)
σ′=σ+φp*

Where, p* is the pore water pressure difference; φ is the parameter.

Jennings (1961) [11] proposed to utilize the following equation to express the effective stress formula for unsaturated soils:

(3)
σ′=σ+βp**

Where p** is the negative pore water pressure; β is a statistical coefficient reflecting the contact area between grains, which can be determined by test.

The above three formulas are most widely used with Bishop’s effective stress formula, and they have in common the use of a single strain-effective stress to analyze the strength and deformation of unsaturated soils. However, Jennings and Burland (1961) [12] pointed out that Equation (3) could not explain the phenomenon of wetting and deformation of unsaturated soil by water immersion. Fredlund et al. considered the effects of net stress (σ–ua ) and matric suction (s = ua–uw) separately in unsaturated mechanics and proposed the following bivariate strength equations for unsaturated soils:

(4)
τf=c′+(σ−ua)tanϕ+(ua−uw)tanϕb
Fig. 1.

Three states of connectivity within unsaturated soils. a) Completely closed state of the gas phase; b) Partially contact state of the gas phase; c) Fully connected state of the gas phase.

https://www.pjoes.com/f/fulltexts/175058/PolJEnvironStud-2024-33-3-2603-2611_Fig_1_min.jpg

Where c' is the effective cohesion, tan ϕ is the effective internal friction coefficient for net normal stress, and tan ϕb is the effective internal friction coefficient for matrix suction.

According to Fredlund’s double stress variable theory, scholars have proposed related unsaturated soil constitutive models, and Spanish scholars Alonso et al. (1990) first proposed an elastoplastic constitutive model for unsaturated soils to describe the basic mechanical properties of unsaturated soils, and the most important part of this model is that it can predict the wetting deformation of unsaturated soils. Sundan et al. (2000) [13] verified this model through laboratory tests and proposed a three-dimensional elastic-plastic model for unsaturated soils. In recent years, many scholars have obtained many theoretical models on the strength characteristics of unsaturated soils through different laboratory tests and theoretical analyses: Zhao Yuxin et al. [14] studied the laws of cohesion, internal friction angle, and shear expansion angle of unsaturated soils with the saturation degree in the existing literature, and obtained an evolution model considering the shear strength indexes of unsaturated soils in a wide range of saturation degrees. Hu Xiaorong et al. [15] combined the principle of generalized effective stress for unsaturated soils with the three-shear strength criterion and proposed a three-shear strength criterion for generalized effective stress for unsaturated soils.

Although there are many theoretical analysis models, in actual engineering, unsaturated soils are strongly influenced by the actual engineering environment, the same kind of soil in different engineering environments may show different bearing capacity characteristics, and the bearing capacity characteristics of the soil are a comprehensive reflection in the complex engineering environment. Therefore, the theoretical analysis should be combined with the relevant field test and laboratory test results for comparison, and then arrive at the final bearing capacity index of the corresponding unsaturated soil body. At present, a large number of scholars at home and abroad have given their suggestions for the determination method of bearing capacity of unsaturated clay soil. In summary, there are three main determination methods: standard penetration test, laboratory test and static probe test, but the applicability of each method varies for specific soils. For the Neoproterozoic Upper Neogene semimorphic clay soils in Baoji area, especially in the diffuse beach and terrace fabric on the north and south sides of Weihe River, it is suggested to be safe and reliable when the characteristic value of foundation bearing capacity is 500 kPa through the relevant experimental research [16]. Zihao Zhang analyzed the method of taking the characteristic value of bearing capacity of soft clay soil in Hangzhou by studying the relationship between the microstructural characteristics of the soil and the macroscopic mechanical properties [17]. Wang Xiaofeng Using the survey data of Hulushan Bay sea area, Changxing Island, Dalian, the main physical indexes of 742 groups of clay soils were correlated and analyzed with the comparative data of bearing capacity, and the method of evaluating the characteristic value of the bearing capacity in the actual survey work was proposed [18].

Fig. 2.

Field soil sample.

https://www.pjoes.com/f/fulltexts/175058/PolJEnvironStud-2024-33-3-2603-2611_Fig_2_min.jpg
Table 1.

Basic physical parameters.

Natural DensityGranular DensityNatural Moisture ContentNatural Pore RatioPorosityLiquid LimitPlastic LimitPlasticity IndexDry DensityFree Expansion Rate
ρ (g/cm3)ρS (g/cm3)ω (%)en (%)WL (%)Wp (%)Ipρd (g/cm3)Fs (%)
1.822.7633.31.12150.564.228.136.11.3720

In view of this, this paper takes an unsaturated clay soil of a highland airport in the mountainous area of southwest China as the research object, and for the deep clay soil appearing in this range, three methods of standard penetration test, laboratory test, and static probe test are adopted to determine its bearing capacity, and the bearing capacity determined by the three methods are compared and studied, to select the optimal method suitable for determining the bearing capacity of unsaturated clay soil of a highland airport.

Overview of Soil Samples

Soil samples were selected from the mileage range of K12+000~K12+840 of a highland airport in the mountainous area of southwest China, and the field soil samples are shown in Fig. 2. The average basic physical property parameters of the soil are shown in Table 1.

Results and Discussion

Standard Penetration Test

The representative standard penetration test data of this section were selected for analysis, and to compare with the other two methods, the data were selected from soil samples with a depth of about 6~7 m for analysis, and a total of 22 groups were selected. The test data are shown in Table 2.

Referring to the relevant industry standard: When the standard value of foundation bearing capacity is determined according to the standard penetration test hammer blow number N, the field test hammer blow number shall be corrected by the following formula to take the N value to the nearest whole digit,

(5)
N=u−1.645σ

Where u stands for the average number of hits in the number of holes after correction for rod length; the calculation is corrected for rod length according to the method of the Engineering Geology Manual, and the correction factor is shown in Table 3. The σ stands for the standard deviation, and the interpolation calculations were then performed according to Table 4.

From Table 3, it can be seen that there is a relationship between the rod length correction coefficient and the rod length as shown in Fig. 3.

According to the relationship curve in Fig. 3, the corrected number of strikes at each test point can be obtained. According to the above specification method, according to the requirements of Table 4 (relevant parameters stipulated by national industry standards), the characteristic value of the bearing capacity of the cohesive soil body in this section obtained according to the standard penetration test can be interpolated and calculated as 179.18 kPa, and the results of the calculation are shown in Table 5.

Table 2.

Standard penetration test data.

NumberMileageStarting Depth (m)End Depth (m)Measured Number of HitsTotal Rod Length (m)Rod Length Correction FactorCorrection Number of Hits
1K12+0006.156.45207.680.87182982917.44
2K12+04013.4513.751514.50.77522940511.63
3K12+0807.057.35149.20.8443810711.82
4K12+1205.756.05168.70.852874913.65
5K12+1606.857.15149.30.84273781111.80
6K12+2005.756.05148.30.86002916212.04
7K12+2407.758.051310.10.83019461610.79
8K12+2809.259.551311.70.80784249610.50
9K12+3206.656.95149.20.8443810711.82
10K12+3605.856.15168.40.85820878113.73
11K12+4006.456.75138.90.84942020611.04
12K12+4406.456.75129.10.84604228910.15
13K12+4807.157.45149.70.83633686511.71
14K12+5206.156.45238.150.86280127519.84
15K12+5606.556.85198.550.85551844516.25
16K12+6007.157.45179.150.8452094114.37
17K12+6407.657.95169.650.83712239713.39
18K12+6806.757.05208.750.85200383817.04
19K12+7207.357.65159.350.84192279612.63
20K12+7607.758.05179.750.83555537314.20
21K12+8008.158.451710.150.82944399714.10
22K12+8405.956.25167.950.86657786713.87
Table 3.

Standard penetration test rod length correction factor.

Rod Length (m)≤369121518212530405075
α1.000.920.860.810.770.730.700.700.680.640.600.50
Table 4.

The standard value of bearing capacity for cohesive soil.

N357911131517192123
fk(kPa)105145190235280325370430515600680

[i] fk(kPa) 105 145 190 235 280

Table 5.

The test data and the recommended value of bearing capacity.

Rock Soil NameNumber of Tests (n)Basic ValueStandard Deviation σN (Calculated by N = μ-1.645σ)Eigenvalue of Bearing Capacity (kPa)
maxminμ
Cohesive Soil2222.755.2512.143.216.87179.18
Fig. 3.

The correction coefficient varies with the length of the rod.

https://www.pjoes.com/f/fulltexts/175058/PolJEnvironStud-2024-33-3-2603-2611_Fig_3_min.jpg

Laboratory Test

When analyzing according to this method, due to the concentration of the laboratory test data for each parameter, four representative groups of data were selected for analysis, and the raw data are shown in Table 6.

Referring to the relevant industry codes, when determining the standard value of foundation bearing capacity based on the average value of physical and mechanical indexes of the soil obtained from laboratory tests, the basic value of bearing capacity in the codes shall be multiplied by the regression correction factor by the following provisions:

Regression correction factors, as specified in the relevant industry codes, i.e., Equation (6).

(6)
Ψf=1−(2.884n+7.918n2)δ

Where Ψf: - regression correction coefficient.

n-Number of participating statistics for the soil indicators on which the tables are based.

δ-Coefficient of variation.

The results of the data analysis calculations are shown in Table 7.

According to the calculation results in Table 7, it is possible to analyze the pore ratio and liquidity index as two indicators for the evaluation of the bearing capacity of clay soils (Table 8).

Referring to Table 8, taking the pore ratio e = 1.15 and IL = 0.08, it can be interpolated to obtain that the bearing capacity of the clay soil under study is 187.2 kPa.

Table 6.

Basic physical parameters of cohesive soil.

Test No.Test MileageTest DepthNatural DensityParticle DensityNatural Moisture ContentNatural Pore RatioPorositySaturationLiquid Limit (10 mm)Plastic LimitPlasticity IndexLiquidity IndexWater Content RatioDry Density
h(m)ρρSωenSrWLWpIpILawρd
(g/cm3)(g/cm3)(%)(%)(%)(%)(%)(g/cm3)
1K12+000 Center7.00~7.401.732.7643.61.29156.493.288.24345.20.010.491.20
2K12+220 Center7.00~7.401.752.7638.61.18654.389.889.644.345.3-0.130.431.26
3K12+440 Center7.50~7.901.822.7633.31.02150.590.064.228.136.10.140.521.37
4K12+600 Center7.60~8.001.802.7337.51.08552.094.455.630.425.20.280.671.31
Table 7.

Experimental data processing calculation results.

Name of IndexNumber of Tests (n)Basic ValueStandard Deviation σCoefficient of Variation δRegression Correction Coefficient ψf
MaxMinAverage Value μ
Natural Pore Ratio e41.291.021.150.120.100.80
Natural Pore Ratio IL40.28-0.130.080.182.34-3.53

[i] Ratio IL

Table 8.

Basic value of bearing capacity of cohesive soil (kPa).

Second Indicator Liquidity Index IL First Indicator Porosity Ratio e00.250.500.751.001.20
0.5430390(360)
0.6475360325295(265)
0.7400295265240210170
0.8325240220200170135
0.9275210190170135105
1.0230180160135115
1.1200160135115105
Table 9.

Static cone penetration tests data.

NumberMileageDepthDesignationBearing Capacity (kPa)Compression Modulus (MPa)
J-12-03K12+0007Unsaturated clay soil322.1416.6
J-12-04K12+0402.7Unsaturated clay soil327.7417.1
J-12-06K12+0803.2Unsaturated clay soil355.5719.7
J-12-08K12+1203.4Unsaturated clay soil326.0617
J-12-09K12+1604Unsaturated clay soil338.8118.1
J-12-10K12+2003.6Unsaturated clay soil360.520.1
J-12-11K12+2400.8Unsaturated clay soil168.526
J-12-12K12+2801.6Unsaturated clay soil221.59
J-12-13K12+3206.8Unsaturated clay soil355.6719.7
J-12-14K12+3605.5Unsaturated clay soil237.2310
J-12-16K12+4002.1Unsaturated clay soil276.0412.8
J-12-18K12+4400.4Unsaturated clay soil209.468.3
J-12-19K12+4801.4Unsaturated clay soil226.29.3
J-12-20K12+5202Unsaturated clay soil28613.6
J-12-21K12+5600.5Unsaturated clay soil379.4222
J-12-22K12+6001.3Unsaturated clay soil257.6611.5
J-12-23K12+6400.5Unsaturated clay soil314.6616
J-12-24K12+6800.4Unsaturated clay soil246.6910.7
J-12-25K12+7203.9Unsaturated clay soil350.9819.2
J-12-26K12+7604.6Unsaturated clay soil335.2917.8
J-12-27K12+8009.9Unsaturated clay soil251.7211
J-12-28K12+8405Unsaturated clay soil307.9215.4
Table 10.

The comparative of results for the three test methods.

Test Method Bearing Capacity(kPa)Standard Penetration TestLaboratory TestStatic Cone Penetration Test
179.18187.20214.19

Static Cone Penetration Test

By the field static cone penetration tests, the resulting data for the cohesive foundation soils in this section are shown in Table 9.

According to Table 9, it can be concluded that the average value of the bearing capacity of the cohesive foundation soil according to the static cone penetration test is 214.19 kPa.

The results of the three test methods are compared in Table 10.

As can be seen from Table 10, among the three test methods, the bearing capacity value determined by the static cone penetration test method is the largest, the bearing capacity determined by the laboratory test method is the second largest, and the bearing capacity value determined by the standardized penetration test method is the smallest. Among them, the bearing capacity values determined by the standardized penetration test method and the laboratory test method are close to each other, but the values determined by the laboratory test method are slightly larger than those determined by the standardized penetration test method, which is because that the laboratory test is often a process of taking samples from the field and then sending them to the test laboratory, where the moisture content of the soil, dry density, etc., changes, which in turn results in a larger value of the tested bearing capacity, while the standardized penetration test is an instantaneous on-site test. But even so, as can be seen from Table 10, the difference does not exceed 5%, indicating that the values determined by these two methods can be used as the recommended values for foundation bearing capacity of unsaturated cohesive soils for high plateau airports. From Table 10, it can be seen that the difference between the static cone penetration test method over the other two methods is greater than 5%, which indicates that for the soil samples under study, this method is not recommended as a bearing capacity analysis when the test data are analyzed mathematically and statistically.

Conclusions

For the unsaturated clay soils occurring in the high plateau airports in the mountainous areas of southwest China, three methods of standard penetration test, laboratory test and static cone penetration test were used to determine their bearing capacity, and the bearing capacity determined by the three methods were compared and studied, and the following conclusions were obtained:

(1) Among the three test methods, the maximum value of bearing capacity was determined by the static cone penetration test method, followed by the laboratory test method, and the minimum value of bearing capacity was determined by the standard penetration test method. The bearing capacity values determined by the standard penetration test method and the laboratory test method are close to each other, but the value determined by the laboratory test method is slightly larger than that of the standard penetration test method.

(2) For the studied soil, it is recommended that the bearing capacity be determined by the standard penetration test method and the laboratory test method rather than the static cone penetration test method.