Abstract
In order to probe into one simplified method to predict the shear strength of Shaoxing unsaturated silty clay, the test method combining unsaturated soil consolidation instrument and conventional direct shear instrument is used to study the shear strength, and the method is compared and verified with the results of equal suction direct shear test. The research results show that the soil water characteristic curve fitted by the measured data points and VG model has obvious stage characteristics in the range of 0~38 kPa, 38~910 kPa, and 910~10000 kPa. The shear strength of unsaturated soil measured by consolidation meter combined with conventional direct shear test is in good agreement with that measured by equal suction direct shear test in the range of 0~500 kPa. The results show that the shear strength, total cohesion, and effective internal friction angle of soil increase slightly with the increase of matric suction in the range of 0~38 kPa. When the matric suction increases from 38 kPa to 500 kPa, the shear strength and total cohesion force of the soil have similar stage characteristics with the SWCC, which first increases and then tends to be stable, while the effective internal friction angle changes slightly. Finally, taking the air-entry value as the demarcation point, an improved model of unsaturated shear strength is proposed by analyzing the error value. Compared with the measured value, the absolute value of relative error is basically kept in the range of 5%~10%, which is close to the measured value.
1. Introduction
Unsaturated soils widely exist in nature, and their engineering properties are quite different from saturated soils. The soils above the groundwater level are generally considered as unsaturated soils. With the change of seasons and rainfall, the unsaturated state is constantly changing. The failure of foundation and slope occurs easily, and it brings great threat to the safety of buildings and structures [1, 2]. Ciabatta et al. [3] have shown that a large amount of rainfall caused by changing global climate may cause more frequent slope failure. It has serious impact on the safety of human life and property. Rahardjo et al. [4] studied the residual soil slope, and they found that the high and steep residual soil slope could be safe in the dry period and only collapsed in the rainy season. The phenomena indicate that the stability of the slope is related to the additional shear strength of the unsaturated zone above the groundwater level. The infiltration of rainwater during rainfall affects the negative pressure of the pore water which decreases the additional shear strength of the unsaturated soils and leads to slope failure easily [1, 5–7].
In order to improve the properties of saturated and unsaturated soils, some scholars try to add inorganic cementitious materials such as cement into the soils [8, 9] or add cement, nano, fiber, and other mixed materials [10–12]. Obviously, in theory and in practice, it is very important to understand and investigate the characteristics especially the shear strength of unsaturated soils directly, which exist widely in nature. Mohr-Coulomb strength theory has been confirmed by a large number of tests and engineering practice. It can accurately determine the shear strength of saturated soil. However, at the beginning of the 20th century, the strength of unsaturated soil was similarly predicted by the Mohr-Coulomb shear strength formula for saturated soil. Until later, with the development of the mechanics of unsaturated soil, many scholars have done a lot of research on the relationship between the shear strength of unsaturated soil and matric suction [13–18]. Idinger and Wu [19] showed that the matric suction reduced the shrinkage effect of granular soil and increased the total cohesion, while the effective angle of internal friction did not change too much; Rasool and Kuwano [20] studied the influence of matric suction on the shear strength of unsaturated silty soil, and they found the relationship between matric suction and shear strength was nonlinear; Pujiastuti et al. [21] pointed out that the total cohesion and effective angle of internal friction of unsaturated sand will increase with the increase of matric suction, but the shear strength increase a little. Khalili et al. [22] studied the mechanical behavior of unsaturated soils and found that the total cohesion increased linearly within the range of high matric suction. Çokça and Tilgen [23] studied the relationship between the shear strength of Ankara compacted clay and matric suction and found that the matric suction and shear strength were correlated positively. Inevitably, Mohr-Coulomb shear strength formula will lead to errors when it is used to predict the shear strength of unsaturated soils.
In geotechnical engineering practice, permeability coefficient, shear strength, and other relevant indexes describing the behavior of unsaturated soils can be obtained indirectly by using the soil-water characteristic curves [24–30]. According to the influence of matric suction on the cohesion and internal friction angle of unsaturated soil, many scholars used the Mohr-Coulomb strength failure criterion of saturated soil and the measured shear strength parameters of saturated soil and soil-water characteristic curve to predict the shear strength of unsaturated soil. For example, Vanapalli et al. [31] proposed a nonlinear shear strength calculation model related to soil-water characteristic curve of soil mass. Oberg and Sallfors [32] presented a formula for calculating the shear strength of unsaturated soils in combination with soil saturation and matric suction. Tang et al. [33] discussed the feasibility of using conventional direct shear instrument combined with the initial soil-water characteristic curve of 0 pressure to estimate the strength of unsaturated soil.
This paper conducted three experiments, respectively, for the unsaturated soil consolidation test, conventional direct shear test, and equal suction direct shear test. The unsaturated soil consolidation test can obtain soil-water characteristic curve, which describes the relationship between matric suction and water content of soils and reflects the water-holding capacity of soils under suction [33]. The conventional direct shear test is to obtain the shear strength of the soil with different water contents, which reflects the relationship between the different water content and the shear strength; the equal suction direct shear test can directly obtain the shear strength of the soil under different matric suction states. It is difficult to measure the required test analysis data in a short time because of the time-consuming of the constant suction direct shear test. However, the unsaturated soil consolidation test combined with conventional direct shear test can achieve the same test purpose as the direct shear test with equal suction, and the time consumption is relatively short. In this paper, the influence of matric suction on the shear strength of soil is analyzed by combining the unsaturated soil consolidation test with the conventional direct shear test and by comparing and checking the partial equal suction direct shear test. At the same time, based on the soil-water characteristic curve fitted by VG model, the shear strength prediction values of three typical models are calculated, and the model is improved through error analysis.
2. Test Materials and Preparation
2.1. Test Materials
The test soil sample is taken from the bottom of a cutting foundation pit in Keqiao District, Shaoxing, and the basic physical property indexes of the test materials are shown in Table 1. It can be seen from the table that the particle size of the soil sample is mainly between 0.005 mm and 0.075 mm, among which the average sand content is 21.43%, the average silt content is 64.52%, the clay content is 14.05%, and plastic index is 12.5, which is between 10 and 17. According to the code of Engineering Classification Standard of Soil (GB/T50145-2007) [34], it is judged to be low liquid limit silty clay.
2.2. Sample Preparation
In this experiment, remolded soil samples are dried, crushed, and sieved for 2 mm, and then, the natural stacking quartering method is used for diagonal sampling. Through theoretical calculation, the required distilled water is evenly sprayed on the surface of remolded soil sample and fully mixed, and then, the soil sample is prepared into a loose shape, sealed in a moisture retaining plastic bag, and placed in a moisture container for more than 48 h, so that the moisture in the soil can be fully transported and mixed evenly. According to the requirements for soil sample preparation in Standards for Geotechnical Test Methods [35], a disc-shaped sample with a diameter of 61.8 mm and a height of 20 mm is prepared by the layered compaction method with different dry density.
2.3. Test Scheme
The unsaturated soil consolidation instrument (SDSWCC) used in the test is made in Nanjing TKA company, and it is shown in Figure 1. According to the code of Technical specifications for highway subgrade construction [36], the compactness of urban secondary trunk roads is not less than 93%, and branch roads are not less than 90%. The water content is close to the optimal water content %. Therefore, this test selects the lowest compaction degree () of the silty clay of the cutting base. In various natural unsaturated soils, the matric suction range that affects the shear strength of the soil is usually 0~500 k Pa [17, 37–39]. Therefore, in this study, the matric suction range of the unsaturated soil consolidation test is taken as 0~500 kPa. The saturated sample is put into the unsaturated soil consolidation instrument, and apply the matric suction path of 1 kPa, 25 kPa, 50 kPa, 75 kPa, 100 kPa, 150 kPa, 200 kPa, 250 kPa, 300 kPa, 400 kPa, and 500 kPa to the sample in sequence to dehydration test. The stability standards of deformation and drainage during the test are as follows: the vertical deformation and displacement do not exceed 0.01 mm/2 h and 0.01 cm3/2 h, respectively, and the suction time of each stage is not less than 48 h. In order to reduce the error of displacement calculation, the pipeline at the bottom of the ceramic plate is bubble washed every 1 d during the test, so as to eliminate the influence of the air at the bottom of the ceramic plate.

This test is unsaturated soil direct shear test with the same degree of compaction (, g·cm-3) under different water content of the test instrument is made from Nanjing TKA production of fully automatic saturated soil direct shear instrument, as shown in Figure 2. Based on the soil-water characteristic curve obtained in the above test (Figure 2), samples are taken according to the volume water content corresponding to the matric suction of 25 kPa, 38 kPa, 50 kPa, 75 kPa, 100 kPa, 150 kPa, 200 kPa, 300 kPa, 400 kPa, and 500 kPa, and the volume water content corresponding to the suction path is 43.3%, 41.8%, 40.4%, 37.4%, 34.8%, 30.8%, 28.8%, 24.4%, 22.3%, and 20.9%, respectively. Three samples of each water content are prepared, and vertical stress of 50, 100, and 200 kPa is imposed for the direct shear test. The shear rate in the test is 0.0032 mm/min. The water content is measured by the drying method. To ensure the water content retaining constant, an impermeable plastic film is placed between the permeable stone and the soil sample during the test to reduce water loss.

The equal suction direct shear test equipment is from unsaturated soil direct shear instrument produced by Nanjing TKA LTD., as shown in Figure 3, which uses the axis translation technology to realize the direct shear test of matric suction control for samples with different saturations. In order to carry out comparative study, the samples are prepared with the optimal water content of 18.4% and compactness of 0.90 ( g·cm-3). The test consists of three stages, respectively, suction balance stage, equal suction consolidation stage, and equal suction shear stage. Combined with the regional characteristics of Shaoxing unsaturated silty clay, the matric suction is set at 0 kPa, 25 kPa, 50 kPa, 75 kPa,100 kPa, 150 kPa, 200 kPa, 300 kPa, 400 kPa, and 500 kPa, and the net vertical stress is 50 kPa, 100 kPa, and 200 kPa, respectively. In order to adapt to the characteristics of homogeneous and slow dissipation of pore pressure and pore water pressure in unsaturated soil during direct shear test, the shear rate of the sample is 0.0032 mm/min and the horizontal displacement is 6 mm. After each sample is sheared, part of the soil is taken from the shear plane, and the water content after shearing is measured by drying method.

3. Test Results and Analysis
3.1. Unsaturated Soil Consolidation Test Results and Analysis
Under the action of different matric suction, the mass water content corresponding to the equilibrium and stable of the sample by the consolidation test of unsaturated soils is shown in Table 2. The measured mass water content is converted into the corresponding volume water content to establish the relationship between the volume water content and the matric suction, as shown in Figure 4. It can be seen that the volume water content shows a significant downward trend with the increase of matric suction. At the beginning, the volume water content drops slowly, then the curve steepens, and the rate of decline increases obviously. However, with the further increase of matric suction, the decreasing rate of moisture content in the sample volume gradually slows down.

Since most of the model equations describing the soil-water characteristic curve are based on experience and the shape of the curve, many scholars have proposed different model equations. The representative soil-water characteristic curve equations include Gardner equation [40], Brooks-Corey model [41], Van Genuchten model [29], and Fredlund and Xing model [30]. Among them, the parameters of the Van Genuchten model have a dear physical meaning, and its fitting effect is perfect, so it is widely used in practice. Therefore, the curve between matric suction and volume water content of unsaturated silty clay in this paper is fitted with the Van Genutchen model, and its expression is as follows: where is the volumetric water content (%); is the saturated volumetric water content; is the residual volumetric water content; is the matric suction (kPa); is the parameter related to the inlet value (kPa-1); is parameters related to the degree of drainage; is a parameter related to residual water content, where .
The logarithm of matric suction is made as abscissa, and volume water content is taken as ordinate, and VG model is used to fit the test data. The soil-water characteristic curve fitted by measured data points and VG model is shown in Figure 4. Among them, matric suction 0~500 kPa is the measured value fitting range, and 500~10000 kPa is the model prediction range. Basic parameters of the fitting are shown in Table 3. Among them, the parameters , , , , and are obtained by inputting the VG model formula and the corresponding values of matric suction and water content measured in the experiment into Origin Software for fitting. As can be seen from Figure 4, the soil-water characteristic curve fitted by VG model shows an inverse “S” shape with obvious three stages. This result is similar to the law of dividing the soil-water characteristic curve dehydration process into three stages: boundary effect zone, transition zone, and unsaturated residual zone. The air-entry value of soil sample is 38 kPa, the residual volume water content is 12.4%, and the matric suction corresponding to the turning point of soil-water characteristic curve is 150 kPa. Within the value of 500 kPa matric suction, the soil-water characteristic curve fitted by the measured value shows obvious two-stage characteristics, which corresponding to the boundary area and transition area of the typical soil-water characteristic curve, respectively. In the boundary area ( kPa), the volume water content of soil did not decrease obviously with the increase of matric suction. This is because the initial water content of the sample is close to saturation, and the pore pressure is small which the gas in the soil can only be suspended in water in the form of closed bubbles and flow together with the water. At this time, the soil loses less water and is close to the saturated soil. When the matric suction increases to the air-entry value and reaches the transition zone (), the soil sample begins to lose water rapidly at a certain rate. This is due to the increase of pore pressure, and the water in the soil will be gradually discharged by gas and occupied a larger pore channel. With the further increase of pore pressure, more water will be discharged from the soil, and the soil water content will decrease rapidly. When the matric suction is in the unsaturated residual zone ( kPa), and the pore pressure continues to increase and the water in the soil is little or no longer discharged, it is considered that the soil has reached the residual water content.
Compared with the typical soil water characteristic curve, due to the limitation of the test instrument, the maximum air intake value of the clay plate is 500 kPa, and the data corresponding to the residual zone cannot be measured by the instrument, so the residual area of the sample can only be predicted by fitting. At present, many definitions of residual matric suction are based on experience. They do not have clear physical meaning and theoretical support. For example, Sillers and Fredlund and Sillers et al. [42, 43] think that the matric suction corresponding to the residual water content is 3000 kPa. They also think that the residual state is the state that pore water in soil changes from capillary action to adsorption force. This definition is quite popular at present. Based on the recommendations of Sillers and Fredlund and Sillers et al. [42, 43], this paper assumes that the residual matric suction is 3000 kPa for fitting, as shown in the curve of the end of measured value to point C in Figure 4. Obviously, the volume water content corresponding to the matric suction of 3000 kPa is not the residual volume water content. Through repeated fitting calculation, the residual matric suction is finally assumed to be 10000 kPa, as shown in the curve to point D in Figure 4, the fitting effect of this region is better, and the corresponding residual volume water content is 12.4%.
3.2. Direct Shear Test and Result Analysis
3.2.1. Relationship between Shear Strength and Matric Suction
Table 4 shows the corresponding shear strength values of the conventional direct shear test under different volume water contents and the equal suction direct shear test under different matric suction. Among them, the corresponding matric suction in the data column of the conventional direct shear test is obtained from the matric suction corresponding to the volume water content in the soil-water characteristic curve.
The corresponding shear strengths of the conventional direct shear test and equal suction direct shear test under different matric suction are shown in the same figure, as shown in Figure 5. It can be seen from the figure that the shear strength of unsaturated silty clay samples changes with matric suction at different vertical normal stresses, and the change trend of boundary effect zone and transition zone is basically consistent with the soil water characteristic curve. In the boundary effect zone (0~38 kPa), a small matric suction can still slightly increase the shear strength of the soil, when the matric suction is less than the air-entry value comparing to the fully saturated state. The soil is in an unsaturated state, in the transition zone (38~500 kPa). The shear strength changes in two stages with the change of the matric suction. When the matric suction increases from 38 kPa to 200 kPa, the shear strength increases rapidly and the corresponding volume water content at this stage is 41.8%~28%, and the lower limit volume water content is close to that of plastic limit. The increase rate of shear strength gradually decreases after the matric suction exceeds 200 kPa. Combining with the soil-water characteristic curve, it can be found that the transition zone can also be divided into two stages, namely, the main transition zone and the secondary transition zone. Taking the matric suction as 200 kPa as the demarcation point, the shear strength in the main transition zone increases rapidly with the matric suction comparing with the secondary transition zone. This is similar to the results of Çokça and Tilgen [44], Bai and Liu [45], Pandya et al. [13], and Khalili et al. [22].

It can be found that the shear strength change trends measured by the two methods are relatively consistent comparing the two test results. When the matric suction is less than the air-entry value, the shear strength of the equal suction direct shear test is basically equal to that of the conventional direct shear test. When the matric suction is in the transition zone, the shear strength of equal suction direct shear test is generally higher than that of conventional direct shear test. Due to the limited range of the matric suction force applied by this test instrument, the results of the matric suction force exceeding 500 kPa cannot be compared.
3.2.2. Relationship between Shear Strength and Matric Suction
In order to avoid the limitation of Bishop equation, Sillers et al. [43] used two independent stress state variables to study the mechanical properties of unsaturated soils. Through a lot of laboratory tests and theoretical analysis, the theoretical formula of the shear strength of unsaturated soils with net normal stress and matric suction as two variables: where represents the rate at which the shear strength increases with the matric suction , which is simply referred to as the internal friction angle of adsorption. If is replaced by in the formula (1), the shear strength formula of unsaturated soil can be rewritten as formula (2): where , in which is the intercept of total cohesion, which is referred to as total cohesion in this article. The total cohesive force of unsaturated soil is composed of the cohesive force in the saturated state and the shear strength caused by the matric suction in the unsaturated state.
Based on the above formula, the total cohesion and effective internal friction angle of soil samples at different volume water contents are obtained according to the results of conventional direct shear tests. By transforming the soil-water characteristic curve into the corresponding matric suction, the relationship between total cohesion and effective internal friction angle with matric suction can be obtained. The strength index of the equal suction direct shear test is obtained by applying the corresponding shear strength of different matric suction. Through calculation, the corresponding shear strength parameters of the conventional direct shear test and equal suction direct shear test under different volume water content and matrix suction are shown in Table 5. Among them, the corresponding matric suction in the data column of the conventional direct shear test is obtained from the matric suction corresponding to the volume water content in the soil-water characteristic curve.
Figures 6 and 7 show the relationship between the total cohesive and effective internal friction angle with the matric suction of equal suction direct shear test and conventional direct shear test. It can be seen from the figures that the total cohesion and effective internal friction angle have different rules of change with the matric suction. With the increase of matric suction, the total cohesion increases nonlinearly and tends to be stable gradually, and its characteristics are synchronous with the shear strength. The effective internal friction angle increases slightly with the increase of the matric suction, and the size is close to the internal friction angle when the soil is saturated, and its change characteristics are not synchronized with the shear strength. According to the geotechnical engineering investigation report of Shaoxing area, the typical silty clay in Shaoxing area is marine sedimentary soil, and the particles are in thin sheet shape. The particles are connected by edge to edge and edge to surface, and the arrangement is directional, with the characteristics of flocculent structure. It contains plant debris, a small amount of flaky particles filled in the coarse particles, clay minerals, silt particles, and humus which adsorb by each other. But microscopically, there is still a very small amount of gas present in the tiny corners of individual pores, the resulting very low suction force plays a role in connecting soil particles, and the lubrication of water between the soil particles is reduced. Therefore, compared with the fully saturated state, the total cohesion and effective internal friction angle of the soil are slightly increased. In the transition zone, the soil is in unsaturated state, the air entering the pores increases gradually, the moisture in the pores decreases correspondingly, the matric suction generated by the difference of pore pressure and water pressure increases, and with the change of air-water interface shrinkage membrane, the surface tension generated increases gradually with the increase of air in the pore, so the total cohesion of soil increases. At the same time, as the pore gas entering the soil increases, the moisture gradually decreases, and the lubrication between soil particles is reduced due to moisture, so the effective internal friction angle of the soil increases slightly. However, when the matric suction increased to a certain extent, the water content in the soil would not be discharged any more, and the water content in the soil would remain unchanged, eventually making the total cohesion in the soil gradually stabilized.


By comparing the data points obtained from the two tests, it can be found that in the boundary region ( kPa), the total cohesion and effective internal friction angle under the equal suction direct shear test and the conventional direct shear test have little change, which is basically the same as that under the saturated state. In the transition zone (), when the matric suction force is lower than 100 kPa, the total cohesion measured by the equal suction direct shear test is smaller than the total cohesion measured by the conventional direct shear test, and when the suction direct shear test exceeds 100 kPa, the total cohesion is greater than that measured by the conventional direct shear test. However, the effective internal friction angle measured by the equal suction direct shear test is always greater than the effective internal friction angle measured by the conventional direct shear test.
From the above analysis, it can be concluded that with the increase of matric suction, the shear strength and total cohesion of the soil increase significantly, and the effective internal friction angle changes less. It can be considered that the contribution of matric suction to the shear strength is more reflected in the increase in total cohesion.
4. Prediction Model and Improved Model of Shear Strength Based on Soil-Water Characteristic Curve
According to the above analysis, as the matric suction increases, the effective internal friction angle of the soil changes only slightly, while the cohesive force changes significantly, so the change of the effective internal friction angle is ignored. In order to consider only the contribution of matric suction to shear strength, the effect of net vertical normal stress on shear strength should be excluded. Therefore, if the net vertical normal stress is 0, the increment of cohesion obtained from saturated state to unsaturated state is the contribution value of matric suction to shear strength. Based on the research of unsaturated soils, Khalili and Khabbaz [46], Vanapalli et al. [31], Lamborn [47], and others earlier proposed different prediction models for the shear strength of unsaturated soils, which have wide applicability. Without considering the influence of the net vertical normal stress and only considering the contribution of the matric suction to the shear strength, three models of the contribution of the shear strength can be obtained, namely, model 1, model 2, and model 3, as (4)~(6): where is the contribution value of the matric suction to the shear strength.
According to the results of the direct shear test, the measured shear strength contribution value can be obtained. The internal friction angle at saturated state of soil and the soil-water characteristic curve parameters kPa, %, and % fitted by VG model are substituted into equations (4)~(6) to calculate the contribution value of shear strength under model 1~model 3, as shown in Table 6. According to the results in Table 6, the calculated shear strength contribution values of the three models are compared with the measured values, as shown in Figure 8. It can be clearly seen from the bar graph that the shear strength contribution value calculated by model 1, model 2, and model 3 differs greatly from the measured value and is significantly lower than the measured value. The error analysis of the three models is performed, and the calculation results are shown in Table 6. The relationship between the absolute relative error of the three models and the change of the matric suction is represented by a bar graph, as shown in Figure 9. It can be seen from Figure 9 that the absolute value of the relative error calculated by the three models has a large dispersion. The absolute relative error calculated by model 1 and model 2 reached the maximum in the boundary area, 67.2% and 54.8%, respectively, while the absolute relative error calculated by model 3 is close to 0 in the boundary area, and the fitting effect is good. However, in the transition region, the absolute values of relative errors calculated by model 1 and model 2 are mostly in the range of 20%~35%, while those of model 3 are mostly in the range of 35%~60%. The average errors of the three models are 33.2%, 26.5%, and 40.5%. Due to the large error, the three calculation models cannot be directly used in the transition zone to calculate the contribution of the matric suction to the shear strength, and a new improved model needs to be established.


According to the literature [48], taking the matric suction value corresponding to the air-entry value and the residual volume water content as the dividing line, the soil can be divided into the boundary zone, the transition zone, and the unsaturated residual zone. According to the comparison of the absolute relative error calculated by the above model 1 to model 3, in the boundary zone, that is, the matric suction is less than the air-entry value, only the absolute relative error of model 3 is about 5%, which is in line with the actual allowable error range. Therefore, taking the air-entry value of the soil 38 kPa as the dividing line and based on the formula proposed by Lamborn [47], when the matric suction is less than the air intake value, the relationship between the matric suction and the shear strength contribution value is formula (6); when the matric suction is greater than 38 kPa, as the matric suction continues to increase, the contribution rate of the matric suction to the shear strength gradually decreases. Therefore, 38 kPa is used as the initial zero point, and the subsequent matric suction minus 38 kPa is used as the matric suction increase. According to the conversion relationship between matric suction and volume water content in reference [31], the expression of the contribution of matric suction to shear strength can be given as follows:
Combining equations (6) and (7), the calculation model of the contribution value of matric suction to shear strength can be obtained as follows:
According to the foregoing, substituting °, kPa, , , and into improved model (8), the shear strength under different matric suction can be obtained contribution value and compare with the measured value and predicted values of three models, as shown in Figure 10. The relationship between the absolute value of the relative error of model 1, model 2, model 3, and the improved model and matric suction is shown in Figure 11. It can be seen from the figure that, except for the absolute relative error between the actual measured value and the calculated value of the improved model, which exceeds 10%, the absolute relative error of the remaining measured value and the calculated value of the model is within 10%. The average absolute value of relative error of the improved model is 5.4% through calculation within the range of matric suction measured. Compared with the absolute value of relative error of model 1 to model 3, the mean value of absolute value of relative error obtained by the improved model decreases by 27.8%, 21.1%, and 35.1%, respectively, which shows the effectiveness of the improved model.


Therefore, the improved model formula for the shear strength of Shaoxing unsaturated silty clay is as follows:
5. Conclusions
(1)The soil water characteristic curve fitted by VG model combined with the strength measured by conventional direct shear test is in good agreement with that measured by equal suction direct shear test of unsaturated soil. Within a certain range, the shear strength and total cohesive force of the soil increase nonlinearly with the increase of matric suction, which has similar phase characteristics of the soil-water characteristic curve, and the effective internal friction angle changes little(2)By analyzing the physical meaning of the air-entry value and taking the air-entry value as the demarcation point, an improved model of the shear strength of unsaturated silty clay in this paper is proposed. The absolute value of the error between the shear strength contribution value calculated by the improved model and the measured value is basically within 10%. Compared with the average absolute relative error of model 1 to model 3, the absolute relative error of the improved model decreases by 27.8%, 21.1%, and 35.1%, respectively, which shows the effectiveness of the improve model
The research in this paper can simplify the testing method of unsaturated soil shear strength, save time, and provide new ideas for the popularization and application of unsaturated soil shear strength theory in engineering practice.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This research was funded by the Zhejiang Provincial Natural Science Foundation of China (Grant number Q20E080042) and the Scientific research project of Zhejiang education department (Grant number Y201839206).