ABSTRACT

Materiales de Construcción

0465-2746

Consejo Superior de Investigaciones Cientificas

MC201915_e185

10.3989/mc.2019.03718

Articles

Pile Side Resistance in Sands for the Unloading Effect and Modulus Degradation

Resistencia de la superficie lateral del pilote a la descarga y a la degradación del módulo cortante en suelos arenosos

Pile Side Resistance in Sands for the Unloading Effect and Modulus Degradation

Zhao

C. F.

a Wu

a Zhao

a b * Zhang

Q. Z.

a Liu

F. M.

a Liu

aDepartment of Geotechnical Engineering, Tongji University, (Shanghai, China) bSchool of Engineering, Tibet University, Lhasa, (Tibet, China)

*zhaocheng@tongji.edu.cn

ORCID ID: C. F. Zhao (https://orcid.org/0000-0002-0010-1819); Y. Wu (https://orcid.org/0000-0002-8366-9871); C. Zhao (https://orcid.org/0000-0002-8973-9164); Q. Z. Zhang (https://orcid.org/0000-0002-7007-5871); F. M. Liu (https://orcid.org/0000-0002-5517-162X); F. Liu (https://orcid.org/0000-0001-9349-2400)

30 06 2019

2019

69 334

10.3989/mc.2019.03718

08 04 2018 01 10 2018 09 04 2019

2019

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0) License.

ABSTRACT

A total of 36 groups of sand-concrete interface loading and unloading direct shear tests were used to analyze the mechanical properties of the pile side-soil interface. The test results show that the interface residual shear stress for the same applied normal stress tends to be constant for the rough sand-concrete interface. The initial shear modulus and peak shear stress of the interface both decrease with the degree of unloading and increase with the interface roughness. The maximum amount of interface shear dilatancy increases with the degree of unloading, and the maximum amount of interface shear shrinkage decreases with unloading for the same interface roughness. A pile side resistance-displacement model is established using the shear displacement method. The proposed function considers both the radial unloading effect and modulus degradation of soil around the pile. The effect of radial unloading and interface roughness on the degradation of the equivalent shear modulus is analyzed using a single fitting parameter b. Good agreement of the proposed model is confirmed by applying the direct shear tests of the 36 groups.

RESUMEN

Resistencia de la superficie lateral del pilote a la descarga y a la degradación del módulo cortante en suelos arenosos. Se han realizado 36 series de ensayos de corte directo con carga y descarga para analizar el comportamiento mecánico del interfaz entre la superficie del lateral del pilote y el suelo. Los resultados demuestran que al aplicar un determinado esfuerzo normal, el esfuerzo cortante residual en el interfaz rugoso entre el suelo arenoso y el hormigón tiende a ser constante. Tanto el módulo cortante inicial como el esfuerzo de cizalla máximo disminuyen al incrementar el grado de descarga y aumentan con la rugosidad del interfaz. Dada una misma rugosidad interfacial, el esfuerzo cortante interfacial máximo crece más rápidamente y su reducción máxima disminuye más lentamente al aumentar el grado de descarga. Se ha empleado el «método del desplazamiento del esfuerzo cortante» para desarrollar un modelo que describe la relación entre la resistencia de la superficie lateral del pilote y el desplazamiento. La función propuesta tiene en cuenta tanto la influencia de la descarga radial como la degradación del módulo cortante del suelo que rodea el pilote. Se analiza la influencia de la descarga radial y la de la rugosidad del interfaz sobre la degradación del módulo cortante equivalente mediante un solo parámetro de ajuste, b. Se ha confirmado que el acuerdo entre el modelo y los datos empíricos de los ensayos de corte directo de las 36 series analizadas es bueno.

KEYWORDS Concrete Composite Stability Mechanical properties Characterization

PALABRAS CLAVE Hormigón Composite Estabilidad Propiedades mecánicas Caracterización

1. INTRODUCTION

Piles are generally used to transfer loads from the superstructure to a competent soil or rock layer (1). The use of pile side resistance analyses is becoming more important compared to solely using bearing capacity analyses. Extensive effort has been expended to develop theoretical methods to analyze pile side resistance behavior. These methods fall into four categories: (i) the load transfer method (1–3), which uses the pile side load-transfer function to describe the relationship between pile side resistance and pile-soil interface displacement; (ii) the shear displacement method (4–9), which assumes that the vertical displacement of the soil at any point around the pile is only related to the shear stress at this point and consider the vertical displacement of the soil induced by the shaft shear stress as a logarithmic relationship of the radial distance away from the pile side; (iii) the elastic theory method (10–12), which employs the Mindlin solution under the concentrated load in the elastic half-space to calculate the displacement of the soil. The equilibrium equation is established by the coordination conditions between the displacement of the pile body and soil to obtain the displacement and side resistance of the pile; (iv) the numerical analysis method (13–18), which includes the finite element method (17), boundary element method (13), discrete element method (16) and infinite layer method (14), and generally needs high computational requirements that caused them are not commonly used in practice.

The shear displacement method assumes that the vertical displacement of the soil at any point around the pile is only related to the shear stress at this point, the shear stress transfer causes the settlement of the surrounding soil, and thus the stress and deformation characteristics of the pile-soil system can be obtained (19). Although the conventional shear displacement method has some defects. For example, this method cannot reflect the pile side resistance softening phenomenon (20), but it is still widely used in the analysis of pile side load transfer mechanism due to its easy calculation process and theoretically reflecting the shear deformation properties of the soil around the pile. (7–8) derived the pile side resistance-displacement functions based on the shear displacement method. Guo et al. (21) modified the pile side resistance-displacement functions that derived from the shear displacement, considering the variations of soil stiffness and limiting pile side resistance with depth. Zhu at al. (22) considered the nonlinear elastic properties and modulus degradation characteristics of the soil into the pile side resistance-displacement functions. The stress-strain behavior of natural soils during shear is highly nonlinear, and the elastic modulus generally decreases with increase in shear strain (22).

On one hand, although above scholars considered many factors that affect the pile side resistance developed, the soil radial unloading for the non-displacement pile during boring stage almost has never been addressed. On other hand, the mechanical properties of soil around a non-displacement pile are very different compared with those of loading soil (23, 24). In particular, the radial unloading effect and modulus degradation characteristics of the soil are not considered comprehensively to a great degree in the pile side resistance-displacement model. As a result, it is difficult to quantify the radial unloading effect of the soil around the pile side on the pile side resistance developed through the field full scale pile tests. Furthermore, the static field full scale pile load tests are extensive and time-consuming although it is the most reliable means of assessing a single pile load transfer mechanism (25). However, the indoor model test has the advantages of parameter control, test repeatability, and lower cost (26). Hence, it can be achieved through the concrete-soil interface large direct shear tests to simulate the radial unloading during the boring stage and pile side-soil shear process (27).

In this paper, a series of sand-concrete interface loading and unloading direct shear tests are used to analyze the mechanical properties of the pile side-soil interface. First, the effect of the unloading process and interface roughness on the mechanical properties of the interface are discussed. Then, a pile side resistance-displacement model is established using the shear displacement method. The proposed τ–s function considers the effect of radial unloading during the non-displacement pile boring stage and modulus degradation of soil around the pile. The direct shear test data is used to validate the proposed model. At the end of this paper, the effect of radial unloading and sand-concrete interface roughness on the degradation of the equivalent shear modulus G_eq is analyzed using a single fitting parameter b.

2. LARGE DIRECT SHEAR TEST 2.1. Test apparatus

Although the large-scale direct shear apparatus has some inherent defects, it is frequently used in interface research due to its simplicity in principle and operation. The test apparatus used in this experiment is the large-scale multi-function direct shear test apparatus SJW-200, which is independently researched and developed by Tongji University, Shanghai, China. The test apparatus has a large size shear box with net size 600 mm × 400 mm × 200 mm (length × width × height) and wall thickness of 40 mm. The test apparatus in both the normal and tangential directions is equipped with an advanced server and control system; the range of the test displacement and applied load are enough to satisfy the requirements of this paper. Compared with a conventional direct shear device, this test apparatus has more accurate test results because it can effectively reduce the boundary effect (7). Figure 1 is the schematic diagram of the large-scale direct shear test apparatus.

Figure 1

Schematic of the large direct shear apparatus.

2.2. Soil specimen

The soil specimen was gray silt, which was obtained from a construction project in Shanghai City. The main properties of the soil are listed in Table 1; the grain-size distribution of the soil is shown in Figure 2.

Table 1

Parameters of the silt in tests.

Unit weight γ/(kN·m^–3)	Cohesion c/kPa	Internal friction angle φ/(°)	Water content ω/%	Void ratio e	Compression modulus E_S1-2/MPa
19.5	4	31.5	20	0.754	11.23

Figure 2

Grain-size distribution of gray silt soil (accumulative granulometrical grapher).

2.3. Concrete plate specimen

Because the bored pile side surface is not smooth in actual engineering, it is necessary to consider the influence of the physical form of the sand-concrete interface. The key is how to simulate and define the roughness of the concrete plate. The quantitative criteria for the interface roughness are as follows.

Dove et al. (28) proposed using the surface regular sawtooth to quantify the soil-concrete interface roughness, as shown in Figure 3. The roughness of the interface can be changed by changing the sawtooth angle, height and space position relationship.

Figure 3

Illustration of the interaction between soil particles and a concrete plate (after (28)).

Zhang et al. (29), (30) proposed using the peak-to-valley distance R to quantify the steel plate interface roughness, as shown in Figure 4. The interface roughness can be changed by changing the peak-to-valley distance.

Figure 4

Illustration of the rough steel plate shape and the definition of roughness (after (29), (30)).

Considering the advantages and disadvantages of the above interface roughness definition and the actual situation of the pile wall, the interface roughness model chosen in this paper is shown in Figure 5. Figure 5 shows that the shape of the rough concrete plate surface is a standard tooth, and its profile is trapezoidal. The tooth angle α = 45 ° remains constant, and the length of the bottom edge S₂ of the trapezoidal convex portion is equal to the bottom edge S₃ of the trapezoidal concave portion. The volume of the concave portion is always equal to the convex portion. By changing the tooth height h and keeping the other conditions unchanged, including tooth angle, space position relationship and so on, the roughness of the concrete plate interface can be adjusted. Hence, the interface roughness R can be expressed by the size of h.

Figure 5

Illustration of the interface roughness model in this paper.

Three variations of concrete plates are used in this experiment, namely, R = 0 mm, 10 mm and 20 mm according to the above definition of concrete interface roughness model. The concrete plate for testing was 600 mm long by 400 mm wide by 50 mm thick and is made of C25 concrete, and bidirectional bars with a diameter of φ8 @ 100 are added into it.

2.4. Test procedures

A total of 36 groups of sand-concrete interface direct shear tests with three different degrees of interface roughness and initial normal stress were used. The test soil is first consolidated under an initial normal stress σ_c for 1 hour and then unloaded to a specific applied normal stress σ_s to shear. Both the loading and unloading rates are 0.5 kPa/s. The three programs of loading and unloading for three different roughness of concrete-soil interface are: (i) σ_c = 300 kPa, and σ_s = 300, 250, 200, 150, 100, 50 kPa, respectively; (ii) σ_c = 200 kPa, and σ_s = 200, 150, 100, 50 kPa, respectively; (iii) σ_c = 100 kPa, and σ_s = 100, 50 kPa, respectively. The shear stage begins when the curve of the normal displacement to time is stable, which means the normal displacement is less than 0.001 mm within 1 minute under the applied normal stress. The shear rate is controlled at a constant rate of 2 mm/min. All data regarding the test are collected by a computerized data logging system; the results are monitored and saved using the computer software TEST.

3. TEST RESULTS ANALYSIS 3.1. Effect of loading and unloading on the sand-concrete interface mechanical properties

Figure 6 shows the interface shear stress-shear displacement curves for different initial and applied normal stresses to study the effect of loading and unloading on the interface mechanical properties. For the sake of simplicity, the legend is expressed as initial normal stress σ_c-applied normal stress σ_s. For example, the legend showing 300–100 means the interface is consolidated at the initial normal stress of 300 kPa, and sheared at the applied normal stress of 100 kPa. Figure 6 also shows that the interface peak shear stress τ_f increases with the initial normal stress under the unloading conditions when the interface roughness R = 0 mm and 10 mm, but both τ_f are less than the loading conditions. For R = 20 mm, the interface peak shear stress increases with the initial normal stress. For the rough soil-concrete interface, the interface residual shear stress tends to be a constant regardless of the loading and unloading conditions and increases with the interface roughness. However, the interface residual shear stress tends to be a constant only for the unloading condition in the smooth soil-concrete interface.

Figure 6

Curves of shear stress-shear displacement for loading and unloading conditions: (a), (b) and (c) show the interface roughness at 0 mm, 10 mm and 20 mm, respectively.

The above can be explained by these points. Higher initial normal stress contributes to a higher density of the interface, which results in soil particles around the interface needing a higher shear stress to make them move along the roughness interface and each other. Therefore, interface peak shear stress increases with the initial normal stress at high roughness. If the interface roughness is low, the interface peak shear stress is mainly influenced by the water content of the interface soils rather than the density. Higher initial normal stress lowers the water content of the interface soils, and if the water content is below the optimal water content, the interface peak shear stress under unloading may be less than the loading condition.

The shapes of the curves for loading and unloading are both logarithmic, and the strain softening phenomena of the curves are also consistent between loading and unloading. These phenomena show that the effect of loading and unloading on the sand-concrete interface mechanical properties change the value of interface peak shear stress, but there is no fundamental change in the failure mode of the interface.

3.2. Effect of unloading degree on the sand-concrete interface mechanical properties

Figure 7 shows the shear stress-shear displacement curves for the initial normal stress at 300 kPa unloaded with different applied normal stress to shear. The unloading degree is defined as the ratio of the increment between the initial and applied normal stress to the initial normal stress. Figure 7 also shows that the unloading degree has a significant effect on the sand-concrete interface mechanical properties. The peak shear stress, the shear displacement corresponding to the peak shear stress and the interface initial shear modulus G₀ all decrease with the unloading degree for the same interface roughness R.

Figure 7

Curves of shear stress-shear displacement for different unloading degree conditions: (a), (b) and (c) show the interface roughness at 0 mm, 10 mm and 20 mm, respectively.

3.3. Effect of interface roughness on the sand-concrete interface mechanical properties

Figure 8 shows the fitting curves between peak shear stress and applied normal stress for different interface roughness and initial normal stress. Table 2 shows the interface equivalent frictional angle φ calculated from the fitting curves. Figure 8 and Table 2 also show that the interface equivalent friction angle φ increases positively with roughness for the same initial normal stress; the rate of increase is positive with the initial normal stress. The interface equivalent friction angle φ is mainly affected by the interface roughness and the initial normal stress.

Table 2

The interface equivalent frictional angle φ for different roughness and initial normal stress.

Roughness R /mm	Initial normal stress /kPa
Roughness R /mm	100	200	300
0	45.7°	46.5°	47.1°
10	38.4°	44.5°	47.7°
20	36.5°	43.3°	51.7°

Figure 8

Fitting curves of peak shear stress and applied normal stress for different roughness conditions: (a), (b) and (c) show the initial normal stress at 100 kPa, 200 kPa and 300 kPa, respectively.

3.4. Effect of loading and unloading on the sand-concrete interface shear dilatancy and shrinkage

Figure 9 shows the variation curves of maximum interface, shear dilatancy and shrinkage amount for different conditions. In this paper, the sign of normal displacement corresponding to the interface shear dilatancy is negative, and the sign of normal displacement corresponding to the interface shrinkage is positive. Figure 9 also shows that the maximum amount of interface shear dilatancy decreases with the applied normal stress σ_s, while the maximum amount of interface shear shrinkage increases with the applied normal stress σ_s as the σ_s increases to a value for the same interface roughness and initial normal stress. The maximum amount of interface shear dilatancy increases with the initial normal stress, and the maximum amount of interface shear shrinkage decreases with the initial normal stress for the same interface roughness and applied normal stress.

Figure 9

Variation curves of maximum interface shear dilatancy and shrinkage amount for different conditions: (a), (b) and (c) show the interface roughness at 0 mm, 10 mm and 20 mm, respectively.

4. MODEL DESCRIPTION

The proposed pile side resistance-displacement model curve, as shown in Figure 10, consists of two parts: (i) a nonlinear pre-failure portion and (ii) a perfectly plastic after-failure portion. Figure 10 shows that the pile side resistance increases nonlinearly with the pile-soil interface displacement increased at first. When the pile-soil interface displacement reaches the value s_u, the pile side resistance achieves its peak value τ_f. The pile side resistance will keep the fixed value τ_f with the pile-soil interface displacement increased. The softening or hardening behavior of the soil, as illustrated in Figure 10, is not considered in this paper.

Figure 10

Schematic of the proposed model.

4.1. Theory function of pile side resistance-displacement

The vertical displacement of the soil at any point around the pile is only related to the shear stress at this point, according to the shear displacement method assumptions, and can be expressed as [1] (19):

s ( z ) = f ( τ s ) 1

where s(z) is the vertical displacement of the soil at a point where the depth is z; τ_s is the shear stress of the soil at that point.

Randolph et al. (7) proposed an approximate analytical solution for the pile side resistance-displacement based on the shear displacement method and expressed the solution as [2]:

s ( z ) = τ 0 r 0 ∫ r 0 r m d r r ⋅ G ( τ s , r ) 2

where r₀ and r_m are the radius of the pile and the limiting radius where the shear stress becomes negligible, respectively; τ₀ and G(τ_s, r) are the shear stress of the pile side and the shear modulus of the soil mass, respectively.

Soil Modulus Degradation

The stress-strain behavior of most geomaterials is highly nonlinear at all phases of loading (22). There is a significant reduction in stiffness with increasing strain level. The larger value of soil shear strain, the lower value of soil shear modulus during the shearing process. According to the assumptions of the shear displacement method, the shear strain of the soil around the pile is decreased with the radial distance from the pile increased. Hence, the degradation model of the soil shear modulus G(τ_s, r) along the pile radial direction can be assumed in Figure 11. The equation of the soil shear modulus G(τ_s, r) is expressed as [3]:

Figure 11

Degradation model of the G (τ_s, r).

{ G ( τ s , r ) = G s 0 ⋅ [ M + ( 1 − M ) ⋅ r − r 0 r 1 − r 0 ] ( r 0 ≤ r < r 1 ) G ( τ s , r ) = G s 0 ( r ≥ r 1 ) 3

where r₁ represents the turning point where the soil shear modulus G(τ_s, r) begins to degenerate and the value of r₁ ≤ r_m; M represents the soil shear modulus degradation coefficient.

By introducing Eq. [3] into Eq. [2] and integrating them, Eq. [2] becomes [4]:

s ( z ) = τ 0 r 0 ( ∫ r 0 η d r r ⋅ G s 0 ⋅ [ M + ( 1 − M ) ⋅ r − r 0 r 1 − r 0 ] + ∫ η r m d r r ⋅ G s 0 ) 4

The concept of equivalent shear modulus G_eq is adopted to facilitate the solution of Eq. [4] where the shear modulus varying in Eq. [4] is equivalent to a constant shear modulus G_eq expressed as [5] and [6]:

s ( z ) = τ 0 r 0 ( ∫ r 0 η d r r ⋅ G s 0 ⋅ [ M + ( 1 − M ) ⋅ r − r 0 r 1 − r 0 ] + ∫ η r m d r r ⋅ G s 0 ) = τ 0 r 0 G e q ∫ r 0 r m d r r 5

G e q = ∫ r 0 r m 1 r ⋅ d r ∫ r 0 η 1 r ⋅ G s 0 ⋅ [ M + ( 1 − M ) ⋅ r − r 0 r 1 − r 0 ] ⋅ d r + ∫ η r m 1 r ⋅ G s 0 ⋅ d r 6

The shear modulus G (τ_s, r) (or G_eq) is related to the soil shear stress τ_s around the pile (or pile side resistance τ₀) and soil initial shear modulus G_s0. Fahey et al. (31) proposed a modified hyperbolic model to express the degradation of soil sear modulus G (τ_s, r) expressed as [7]

G ( τ s , r ) G s 0 = 1 − a ⋅ ( τ s τ max ) b 7

where the coefficient a and b are the degree and rate, respectively, of the shear modulus degradation. τ_s and τ_max are the soil shear stress and maximum soil shear stress around the pile. For the equivalent shear modulus G_eq, it is necessary to replace the τ_s and τ_max of the pile side resistance τ₀ and peak pile side resistance τ_f in Eq. [7], expressed as [8]

G e q G s 0 = 1 − a ⋅ ( τ 0 τ f ) b 8

By introducing Eq. (8) into Eq. (5) and integrating them, Eq. [5] becomes [9]

s ( z ) = τ 0 ⋅ r 0 G s 0 [ 1 − a ⋅ ( τ 0 τ f ) b ] ⋅ ln ( r m / r 0 ) 9

Evaluation of Peak Pile Side Resistance <italic>τ</italic> <italic>f</italic>

There are three methods for determining the peak pile side resistance τ_f : (i) the α method belonging to the total stress method (32); (ii) the β method belonging to the effective stress method (33); and (iii) the λ method belonging to the mixing method (34). For the object of this paper, the β method is modified to determine the peak pile side resistance τ_f considering the unloading effect.

The β method can be used to calculate the pile side resistance of both cohesive soil and non-cohesive soil. It was proposed by Chandler et al. (33), and the calculated formula is expressed as [10]

τ f = σ V ' ⋅ K 0 ⋅ tan δ = β ⋅ σ V ' 10

where σ’_v is the average vertical effective stress of the pile side-soil layer; K₀ is the static lateral pressure coefficient; δ is the pile-soil interface friction angle.

For normal consolidated soil, K₀=1-sinφ′, where φ’ is the effective internal friction angle of the soil δ=φ′, and β = (1- sinφ')∙tanφ'. For over-consolidated soil, the effect of over-consolidation OCR should be considered. Mayne et al. (35) proposed a modified formula for calculating the value of β for over-consolidated soil [11]

β = ( 1 − sin φ ′ ) ⋅ O C R sin φ ′ ⋅ tan δ 11

The over-consolidation ratio OCR is calculated as [12]

O C R = σ c / σ s 12

where σ_c and σ_s are the initial consolidation normal stress and applied normal stress in the process of shearing, respectively. The unloading ratio ξ is introduced when considering the unloading effect and is expressed as [13]

ξ = ( σ c − σ s ) / σ c 13

Introducing Eq. [11], Eq. [12] and Eq. [13] into Eq. [10] becomes [14]

τ f = ( 1 − sin φ ′ ) . ( 1 − ξ ) − sin φ ′ ⋅ tan δ . σ v ' 14

Eq. [22] is the peak pile side resistance calculation function that considering the unloading effect.

Evaluation of the influence area of radial limiting radius r <italic>m</italic>

Randolph et al. (23) proposed the following equation as an estimation method to calculate the r_m [15]

r m = C ⋅ L ⋅ ρ ⋅ ( 1 − v ) 15

where C is the empirical coefficient, and its value is 2.5 for a pile in the semi-infinite space and 2.0 for the rigid layer below the pile end 2.5 times the pile length; L is the depth of the pile in the soil; ρ is the ratio of soil shear modulus at mid-depth to that at the pile tip, and ρ = 1.0 for uniform soil; v is the Poisson ratio, and its value generally takes v = 0.2 ~ 0.4 for Shanghai sand.

From Eq. [15], it is relatively simple to determine the r_m, but this equation does not consider the impact of the load level of the pile and the change of r_m with the pile buried depth. Hence, this paper uses the following method to determine the value of r_m.

The downward displacement of the soil at r = r_m is 0, according to the shear displacement method assumptions. Thus, differentiating Eq. [2] on both sides is expressed as [16]

d s ( z ) | r = r m = τ 0 r 0 r m ⋅ G ( τ s , r m ) = τ 0 r 0 r m ⋅ G s 0 = 0 16

Eq. [16] shows that the hold condition is r_m→∞. In the actual engineering, it is impossible to meet the conditions of r_m→∞. Therefore, there exists a treatment measure as follows [17]

d s ( z ) = τ 0 r 0 r m ⋅ G s 0 = η 17

where η is a finite small value. Solving r_m from Eq. (17) can be expressed as [18]

r m = τ 0 r 0 η ⋅ G s 0 18

Eq. [18] shows that the r_m is related to the pile side resistance τ₀ (reflecting the pile load level), the soil initial shear modulus G_s0 (reflecting the soil properties) and the pile radius r₀ (reflecting the pile shape).

By introducing Eq. [14] and Eq. [18] into Eq. [9], the pile side resistance-displacement model considering the unloading effect and soil modulus degradation is obtained as follows [19]:

s ( z ) = τ 0 ⋅ r 0 G s 0 ⋅ ln [ τ 0 / η ⋅ G s 0 ] τ 0 [ 1 − a ⋅ ( τ 0 σ v ' ⋅ ( 1 − sin φ ′ ) ⋅ ( 1 − ξ ) − sin φ ′ ⋅ tan δ ) b ] 19

4.2. Determination of the values of parameters <italic>a</italic>, <italic>b</italic>, <italic>η</italic>, <italic>δ</italic> and <italic>G</italic>s0 Evaluation of Parameters a, b and η

The range of the coefficient a is from 0 to 1.0, and (36) suggests that the results of most cases can be well fitted when a = 0.98. Mayne et al. (37) suggests that the range of the coefficient b is from 0.2 to 0.4 for most types of soil. In this paper, the range of empirical coefficient b is extended to b = 0.02 ~ 0.4.

The value of parameter η can be obtained by experiment. Randolph et al. (7) suggests that the range of value for η is 1 × 10^–7 ~ 1 × 10^–5, its value can take 1 × 10^–6.

Evaluation of Parameters δ

The parameter δ is the pile-soil interface friction angle, and its value varies with the type of pile and soil properties and generally takes as Table 3 recommended.

Table 3

Recommended interface friction angle (after (38)).

Pile material	δ
Rough concrete	φ’
Smooth concrete	0.8φ’ to φ’
Steel	0.5φ’ to 0.9φ’
Timber	0.8φ’ to 0.9φ’

Evaluation of Initial Soil Shear Modulus G <italic>s0</italic>

The initial shear modulus of the soil G_s0 is generally calculated by the E-v model in actual engineering [20]

G s 0 = E 2 ( 1 + v ) 20

where E and v are the elastic modulus and Poisson’s ratio of the soil, respectively. For Shanghai soil, E = 3.5E_s1-2 according to Wang et al. (39). v is the Poisson ratio, and its value generally is v = 0.2 ~ 0.4 for Shanghai sand.

4.3. Algorithm for analysis of the pile side resistance-displacement model

The algorithm for the analysis of the pile side resistance-displacement model can be summarized as follows:

The integral relation of τ-s function is derived by using the shear displacement method, as shown in Eq. [5].

The degradation model of the shear modulus is obtained by using the modified hyperbolic function, as shown in Eq. [8].

Solve the peak pile side resistance τ_f using Eq. [14].

Solve the limiting radius r_m using Eq. [18].

Obtain the formula of the pile side resistance-displacement model by introducing Eq. [8], Eq. [14] and Eq. [18] into Eq. [5], thus becoming Eq. [19].

5. MODEL VALIDATION AND ANALYSIS

The test results obtained by the large-scale direct shear test are used to validate the pile side resistance-displacement model. The parameter r₀ in Eq. [19] cannot be reflected in the direct shear test. The pile-soil interface shear strength τ_f can be derived directly from the test results and do not need to be calculated by Eq. [14]. Hence, Eq. [19] can be transformed as follows [21]

s r 0 = τ 0 G s 0 ⋅ ln [ τ 0 / ( η ⋅ G s 0 ) ] [ 1 − a ⋅ ( τ 0 / τ f ) b ] 21

where τ₀ and τ_f are the current shear stress and peak shear stress in the direct shear test, respectively. Using the parameter τ to replace τ₀ in Eq. [21], s is the shear displacement.

5.1. Comparison curves between the model calculation and test values with different σc

The curves of s/r₀-τ/τ_f are plotted and compared with the τ-s curves obtained by the large-scale direct shear test. The τ/τ_f -s/r₀ curves are plotted on the main coordinate axis, and the τ-s curves are plotted on the sub-coordinate axis. The main coordinate axis and sub-coordinate axis in the figure have strict correspondence. The curves of both model calculation and test values are the pre-failure portion; after-failure portion is not considered for it is assumed perfectly plastic. The range of empirical coefficient b is extended to b = 0.02 ~ 0.4.

Figures 12 a, b and c are the comparison curves for the interface roughness R = 0 mm and σ_c = 100 kPa, 200 kPa and 300 kPa, unloaded to σ_s = 50 kPa. Figure 12 shows that the test values fit well with the model calculation, and the value of empirical coefficient b gradually changed from 0.05 to 0.3 as the σ_c increased. This phenomenon shows that the value of empirical coefficient b increases positively with the initial normal stress σ_c in the same applied normal stress and interface roughness.

Figure 12

Comparison curves of the model calculation and test values for R = 0 mm and σ_s = 50 kPa: (a), (b) and (c) are the initial normal stress σ_c = 100 kPa, 200 kPa and 300 kPa, respectively.

5.2. Comparison curves between the model calculation and test values under different σs

Figures 13 a, b and c are the comparison curves for the interface roughness R = 0 mm and σ_c = 300 kPa unloaded to σ_s = 300 kPa, 200 kPa and 100 kPa, respectively. Figure 13 show that the test values in the loading conditions fit the best with the model calculation values when the empirical coefficient a = 0.98 and b = 0.05. The test curves that fit with the model calculations gradually shift up with unloading degree increase; the value of empirical coefficient b is gradually changed from 0.05 to 0.2 with the unloading degree increased. This phenomenon shows that the value of empirical coefficient b increases with the unloading degree increased.

Figure 13

Comparison curves of the model calculation and test values under R = 0 mm and σ_c = 300 kPa: (a), (b) and (c) are the applied normal stress σ_s = 300 kPa, 200 kPa and 100 kPa, respectively.

5.3. Comparison curves between the model calculation and test values under different <italic>R</italic>

Figures 14 a and b, and Figures 15 a and b are both the comparison curves of the interface roughness R = 10 mm and R = 20 mm, respectively. Figure 14 (a) ~ (b) show that the test values fit best with the calculated values at b = 0.05~0.1 for R = 10 mm and b = 0.02 for R = 20 mm, respectively. Figures 15 a and b show that b = 0.1~0.2 and b = 0.05~0.1 for R = 10 mm and R = 20 mm, respectively. Combined with Figures 13 b and c, Figure 14 and Figure 15 show that the value of empirical coefficient b decreases with interface roughness.

Figure 14

Comparison curves of the model calculation and test values for σ_c = 300 kPa and σ_s = 200 kPa: (a) and (b) show the interface roughness at 10 mm and 20 mm, respectively.

Figure 15

Comparison curves of the model calculation and test values for σ_c = 300 kPa and σ_s = 100 kPa: (a) and (b) show the interface roughness at 10 mm and 20 mm, respectively.

5.4. Effect of unloading and roughness on the equivalent shear modulus degradation

Figure 16 shows the degradation curves of the equivalent shear modulus G_eq, which is calculated from the modified hyperbolic function given by Eq. (8) where a = 0.98 and b is 0.02 to 1. Figure 16 shows that the degradation rate of G_eq decreases with the value of b during the early shearing stage (about τ < 0.05τ_f) and increases with the value of b after the early shearing stage.

Figure 16

Degradation curves of the equivalent shear modulus.

The above analysis (section 5.1 to 5.3) shows that the value of the empirical coefficient b increases with the unloading degree but decreases with interface roughness. Unloading will slow down the degradation of G_eq during the early shearing stage; the degradation rate of G_eq decreases with the unloading degree. After the early shearing stage, unloading accelerates the degradation of G_eq; however, the effect of interface roughness on the degradation of G_eq is opposite to the unloading degree. The increase of interface roughness will accelerate the degradation of G_eq during the early shearing stage and decelerate the degradation of G_eq after the early shearing stage.

6. CONCLUSIONS

In this study, large scale direct shear tests are conducted to analyze the mechanical properties of the pile side-soil interface. The effects of the unloading process and interface roughness on the mechanical properties of the interface are discussed. A pile side resistance-displacement model is developed using the shear displacement method. The proposed function considers both the radial unloading effect and modulus degradation of soil around the pile. The main conclusions are drawn as follows:

The shapes of the interface shear stress-displacement curves for loading and unloading are both logarithmic and the effect of loading and unloading on the sand-concrete interface mechanical properties mainly change the value of interface peak shear stress, but there is no fundamental change in the failure mode of the interface. The interface peak shear stress is mainly influenced by the soil density around the high roughness interface (R > 10 mm) and the peak shear stress increases with the initial normal stress increased. If the interface roughness is low (R < 10 mm), the interface peak shear stress is mainly influenced by the water content of the interface soils and the peak shear stress may decrease with the increased of initial normal stress.

The interface peak shear stress, shear displacement corresponding to the peak shear stress and the interface initial shear modulus G₀ all decrease with the unloading degree for the same interface roughness R. The interface equivalent friction angle φ increases with roughness increased for the same initial normal stress and the rate of increase is positive with the initial normal stress. The maximum amount of interface shear dilatancy increases with the unloading degree increased, while the maximum amount of interface shear shrinkage decreases with the unloading degree increased for the same interface roughness.

The proposed pile side resistance- displacement model considers both the radial unloading effect and modulus degradation of soil around the pile has been validated by the direct shear tests of 36 groups; The coefficient a takes fixed value of 0.98 and b takes value varies from 0.02 to 0.3 can best fit existing τ–s curves from direct shear tests for different unloading degree and interface roughness.

The effect of radial unloading and interface roughness on the interface mechanical properties can be attributed to the degradation of the equivalent shear modulus G_eq, which can be analyzed by a single fitting parameter b. The value of the empirical coefficient b increases with the unloading degree but decreases with interface roughness. Unloading decelerates the degradation of G_eq during the early shearing stage (about τ < 0.05τ_f) and the degradation rate of G_eq decreases with the unloading degree increased. After the early shearing stage (about τ > 0.05τ_f), unloading accelerates the degradation of G_eq. However, the increase of interface roughness will accelerate the degradation of G_eq during the early shearing stage (about τ < 0.05τ_f) but decelerate the degradation of G_eq after the early shearing stage (about τ > 0.05τ_f).

ACKNOWLEDGEMENTS

This work was supported by the National Key Research and Development Plan [grant number 2017YFC0806000], the National Natural Science Foundation of China under Grant [No. 41672265, No. 41572262, No. 41502275]; and Shanghai Rising-Star Program under Grant No. 17QC1400600. The authors are deeply grateful for this support. The anonymous reviewers’ comments have improved the quality of this paper and are also greatly acknowledged.

REFERENCES 1

P. P.

Song

L. H.

Mei

G. X.

Zhao

Y. L.

Generalized nonlinear softening load-transfer model for axially loaded piles

Int. J. Geomech. 2017 17 [8] 04017019.

https://doi.org/10.1061/(ASCE)GM.1943-5622.0000899

Guo

W. D.

Pile capacity in nonhomogeneous softening soil

Soils Found. 2001 41 2 111 120

https://doi.org/10.3208/sandf.41.2_111

Seed

H. B.

Reese

L. C.

The action of soft clay along friction piles

Transactions of the ASCE 1957 122 1 731 754

Chow

Y. K.

Analysis of vertically loaded pile groups

Int. J. Numer. Anal. Meth. Geomech. 1986 10 1 59 72

https://doi.org/10.1002/nag.1610100105

Chow

Y. K.

Axially loaded piles and pile groups embedded in a cross-anisotropic soil

Geotechnique 1989 39 2 203 212

https://doi.org/10.1680/geot.1989.39.2.203

L. L.

Chen

Q. S.

Huang

M. S.

Basack

Hybrid approach for rigid piled-raft foundations subjected to coupled loads in layered soils

Int. J. Geomech. 2017 17 5 04016122

https://doi.org/10.1061/(ASCE)GM.1943-5622.0000825

Randolph

M. F.

Wroth

C. P.

Analysis of deformation of vertically loaded piles

J. Geotech. Eng. Div. 1978 104 12 1465 1488

https://doi.org/10.1016/0148-9062(79)90663-6

Randolph

M. F.

Wroth

C. P.

An analysis of the vertical deformation of pile groups

Geotechnique 1979 29 4 423 439

https://doi.org/10.1680/geot.1979.29.4.423

Zhang

Q. Q.

Zhang

Z. M.

Study on interaction between dissimilar piles in layered soils

Int. J. Numer. Anal. Meth. Geomech. 2011 35 1 67 81

https://doi.org/10.1002/nag.893

Poulos

H. G.

Davis

E. H.

The settlement behaviour of single axially loaded incompressible piles and piers

Geotechnique 1968 18 3 351 371

https://doi.org/10.1680/geot.1968.18.3.351

Wang

Yang

An improved elastic analysis method of pile foundation under vertical loading

Roc. Soil. Mech. 2006 27 8 1403 1406 in Chinese

https://doi.org/10.3969/j.issn.1000-7598.2006.08.036

K. J.

Poulos

H. G.

General elastic analysis of piles and pile groups

Int. J. Numer. Anal. Meth. Geomech. 2000 24 15 1109 1138

https://doi.org/10.1002/1096-9853(20001225)24:15%3C1109::AID-NAG72%3E3.0.CO;2-N

Z. Y.

Han

Boundary element analysis of axially loaded piles embedded in a multi-layered soil

Comput. Geotech. 2009 36 3 427 434

https://doi.org/10.1016/j.compgeo.2008.06.001

Cheung

Y. K.

Tham

L. G.

Guo

D. J.

Analysis of pile group by infinite layer method

Geotechnique 1988 38 3 415 431

https://doi.org/10.1680/geot.1988.38.3.415

Desai

C. S.

Numerical design analysis for piles in sands

J. Geotech. Eng. Div. 1974 100 6 613 635

https://doi.org/10.1016/0148-9062(74)91242-X

Lee

C. Y.

Discrete layer analysis of axially loaded piles and pile groups

Comput. Geotech. 1991 11 4 295 313

https://doi.org/10.1016/0266-352X(91)90014-7

Sheng

D. C.

Eigenbrod

K. D.

Wriggers

Finite element analysis of pile installation using large-slip frictional contact

Comput. Geotech. 2005 32 1 17 26

https://doi.org/10.1016/j.compgeo.2004.10.004

Trochanis

A. M.

Bielak

Christiano

Three - dimensional nonlinear study of piles

J. Geotech. Eng. 1991 117 3 429 447

https://doi.org/10.1061/(ASCE)0733-9410(1991)117:3(429)

Everett

J. P.

Load transfer functions and pile performance monitoring

Int. J. Roc. Mech. Min. Sciences. 1991 30 1 229 234

https://doi.org/10.1016/0148-9062(93)90572-U

Y. H.

Wang

W. D.

J. B.

Bearing deformation of large-diameter and super-long bored piles based on pile shaft generalized shear model

Chin. J. Geotech. Eng. 2015 37 12 2157 2166 in Chinese

https://doi.org/10.11779/CJGE201512004

Guo

W. D.

Randolph

M. F.

Vertically loaded piles in non-homogeneous media

Int. J. Numer. Anal. Meth. Geomech. 1997 21 8 507 532

https://doi.org/10.1002/(SICI)1096-9853(199708)21:8%3C507::AID-NAG888%3E3.0.CO;2-V

Zhu

Chang

M. F.

Load transfer curves along bored piles considering modulus degradation

J. Geotech. Geoenviron. Eng. 2002 128 9 764 774

https://doi.org/10.1061/(ASCE)1090-0241(2002)128:9(764)

Lade

P. V.

Duncan

J. M.

Stress-path independent behavior of cohesionless soil

J. Geotech. Eng. Div. 1976 102 1 51 68

Lambe

T. W.

Stress path method

J. Soil. Mech. Found. Div. 1967 93 6 309 331

Wang

Z. J.

Xie

X. Y.

Wang

J. C.

A new nonlinear method for vertical settlement prediction of a single pile and pile groups in layered soils

Comput. Geotech. 2012 45 118 126

https://doi.org/10.1016/j.compgeo.2012.05.011

Yang

J. J.

2005 Similarity theory and structural model test Wuhan University of Technology Press

Wuhan

in Chinese

Zhao

C. F.

Zhao

Tao

G. X.

Experimental research on the clay-concrete interface shear behaviors considering the roughness and unloading Effect 2018

Proc. of GeoShanghai 2018 International Conf. Fundamentals of soil behaviours

522 530

https://doi.org/10.1007/978-981-13-0125-4_58

Dove

J. E.

Jarrett

J. B.

Behavior of dilative sand interfaces in a geotribology framework

J. Geotech. Geoenviron. Eng. 2002 128 1 25 37

https://doi.org/10.1061/(ASCE)1090-0241(2002)128:1(25)

Zhang

J. M.

Development and application of cyclic shear apparatus for soil-structure interface

Chin. J. Geotech. Eng. 2003 25 2 149 153 in Chinese

https://doi.org/10.3321/j.issn:1000-4548.2003.02.005

Zhang

J. M.

Experimental study on monotonic behavior of interface between soil and structure

Chin. J. Geotech. Eng. 2004 26 1 21 25 in Chinese

https://doi.org/10.3321/j.issn:1000-4548.2004.01.003

Fahey

Carter

J. P.

A finite element study of the pressuremeter test in sand using a nonlinear elastic plastic model

Can. Geotech. J. 1993 30 2 348 362

https://doi.org/10.1139/t93-029

Tomlinson

M. J.

Woodward

1977 Pile design and construction practice 4

A View Point Publications

London

Chandler

R. J.

The shaft friction of piles in cohesive soils in terms of effective stress

Civil Eng and Public Works Review 1968 63 738 48 51

Focht

Vijayvergiya

A new way to predict capacity of piles in clay

Proc. 4th Offshore Technology Conf. 1972 2 865 874

https://doi.org/10.4043/1718-MS

Mayne

P. W.

Kulhawy

F. H.

K-₀OCR relationships in soil

J. Geotech. Eng. Div. 1982 108 6 851 872

Ruiz

M. E.

Study of axially loaded post grouted drilled shafts using CPT based load transfer curves 2005

University of Puerto Rico

Master Thesis, Civil Engineering

Mayne

P. W.

Schneider

J. A.

Evaluating axial drilled shaft response by seismic cone

Foundations and ground improvement, Geotechnical Special Publication 2001 113 655 669

Kulhawy

F. H.

Limiting Tip and Side Resistance: fact or fallacy

Proc. Conf. on Analysis and Design of Pile Foundations 1984

ASCE Convention

San Francisco, California

80 98

https://doi.org/10.1016/0148-9062(86)91224-6

Wang

W. B.

Yang

Elasto-plastic analysis for vertical pile based on extended compatibility of deformation

Chin. J. Geotech. Eng. 2005 27 12 1442 1446 in Chinese

https://doi.org/10.3321/j.issn:1000-4548.2005.12.014