# Influence of Intermediate Principal Stress on Undrained Behavior of Intact Clay under Pure Principal Stress Rotation

Authors: Jian Zhou; Jiajia Yan; Changjie Xu; Xiaonan Gong
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

This study presents the accumulations of the excess pore water pressure and the deformation as well as the noncoaxial behavior of intact soft clay subjected to pure principal stress rotation. Series of tests were carried out by using a dynamic hollow cylinder apparatus to highlight the influence of intermediate principal stress parameterb. It was found that the rate of PWP evolution was greatly influenced by b, but the influence was not monotonous. Specimens under the condition b = 0.75 had the highest accumulation of pore water pressure while under the condition b = 0 had the strongest resistance to the pore pressure generation. PWP accumulated mainly in the first cycle. The failure of specimens under principal stress rotation was controlled by the strain other than the pore pressure. The shear stiffness decreased more quickly with higher b value. The direction of the principal strain increment was strongly dependent on the principal stress increment orientation and less influenced by the b value and the number of cycles.


---

## Body

## 1. Introduction

Loading conditions in geotechnical engineering are very complex; therefore in laboratory test and numerical analysis the stresses on the soil elements need to well represent the real situation. In many scenarios, such as earthquakes, vehicular traffic, and ocean waves [1–3], not only the magnitude of principal stresses but also the direction will change. Studies on this kind of stress-induced anisotropy were stimulated and many were performed through laboratory test, where the volumetric strain and pore water pressure accumulation induced purely by principal stress rotation were highlighted [1–9].In the literature, the principal stress (major and minor stress) rotated in the plane perpendicular to the intermediate principal stress direction. So the influence of intermediate principal stress, normally expressed by the intermediate principal stress parameterb=(σ2-σ3)/(σ1-σ3), needs to be considered. The effect of b on the soil response has been widely studied through many fixed principal stress direction tests. And the results showed that the shear strength, the pore pressure parameter, and the frictional angle of soils are highly dependent on the magnitude of b [10–16].However, the influence ofb has not been well studied in principal stress rotation tests. Ishihara and Towhata [1] were among the first to conduct undrained test of continuous principal stress rotation on sand by using a hollow cylinder apparatus (HCA). In all their tests, b equaled to sin2(α), where α is the angle between the directions of major principal stress and the vertical axis. This was due to the limitation of their apparatus [17]. So their study on principal stress rotation was not independent of b. Since then, pure principal stress rotation tests with constant mean total stress, p, the shear stress, q, and intermediate principal stress parameter, b, were carried out by others on sand and clay [2, 4–8], but most of the tests were conducted under the condition of b=0.5.In recent years, the influence ofb on the sand behavior under the continuous principal stress rotation was studied experimentally by Yang et al. [2]. Undrained tests with pure principal stress rotation on the sand were conducted, b values of 0, 0.5, and 1.0 were studied. The test results showed that b has great influence on the pore pressure and the deformation response during the pure principal stress rotation. The shear stiffness of sand decreased more quickly with cyclic principal stress rotation under b=1.0. Tests with identical stress paths under drained condition were carried out by Tong et al. [7]. The volumetric strain induced by the cyclic principal stress rotation accumulated more quickly under b=1.0. However, the direction of principal strain increment was less influenced by b.In contrast with the relatively large number of studies on sand under the principal stress rotation, the response of intact clay to the pure principal stress rotation has received little attention. Tests on the intact Hangzhou soft clay under the pure principal stress rotation were presented in this study. A wide range ofb values (0, 0.25, 0.5, 0.75, and 1.0) were considered. The influence of b as well as the shear stress was studied through the pore pressure evolution and the strain accumulation. The noncoaxial behavior and the stiffness degradation were also examined with the purpose of formulating more realistic constitutive models for clayey soil.


## 2. Equipment and Test Scheme

### 2.1. Hollow Cylinder Apparatus

Hollow cylinder apparatuses have been successfully used in the soil anisotropy studies. Here, different types of principal stress rotation tests on the intact Hangzhou soft clay were conducted using an advanced Zhejiang University hollow cylinder apparatus (ZJU-HCA). In this HCA the external loadings of axial loadW, inner pressure pi, outer pressure po, and torque MT can be applied and controlled individually to simulate the complex stress paths. The big advantage of simultaneously applying dynamic axial load and torque with a frequency up to 5 Hz provides many opportunities for the complex test. The maximum values of axial load and torque of ZJU-HCA are 12 kN and 200 Nm.The ZJU-HCA can be used for tests on both intact and reconstituted samples with the size of 160 mm × 200 mm × 400 mm and 60 mm × 100 mm × 200 mm (inner diameter × outer diameter × height), respectively. The details of this apparatus and the sampling kit can refer to Shen [18].Corresponding to the axial load, inner cell pressure, outer cell pressure and torque, the stress components ofσz, σr, σθ, and τzθ on the specimen indicated in Figure 1(a) can be calculated according to Hight et al. [19]. These four stress components can be converted into an equivalent set of four stress-related parameters, p, q, α, and b as the following:
(1)p=(σ1+σ2+σ3)3=(σz+σr+σθ)3q=(σ1-σ3)2=(σz-σθ)2/4+τθz2α=12·tan-1(2τθzσz-σθ)b=σ2-σ3σ1-σ3=((2σr-σz-σθ)/2)+((σz-σθ)/2)2+τθz2(σz-σθ)2+4τθz2,
where p is the mean total stress, q is the shear stress, b is the intermediate principal stress parameter (0≤b≤1), and α is the angle between the directions of the major principal stress and the vertical axis (Figure 1(a)). In this study the principal stresses (σ1 and σ3) rotate in the plane perpendicular to σ2.(a) Stress state of hollow cylinder specimens; (b) a completed hollow cylindrical specimen.
(a)
(b)As the four stress components ofσz, σr, σθ, and τzθ can be controlled independently, stress paths with constants  p, q, and b while α increases continuously (pure principal stress rotation) can be achieved [18]. The main objective of this study is to investigate the influence of b on the responses of intact soft clay under the pure principal stress rotation. This study can be well carried out with the support of ZJU-HCA.


### 2.2. Test Material and Sample Preparation

The hand-cut clay blocks were taken from an excavation pit 4 m below the ground level in Hangzhou; the basic physical properties of the clay are summarized in Table1. After the hollow cylindrical specimen (Figure 1(b)) with the size of 100 mm in outer diameter, 60 mm in inner diameter and 200 mm in height was prepared by a special sampling kit [18]; filter paper strips were pasted around the outer surface of the specimen to equalize pore water pressure and reduce the consolidation time. After that, the specimen was assembled into the cell to be saturated with a backpressure of 100 kPa, ensuring the Skempton’s pore pressure parameter B higher than 0.98. After saturation the specimens were isotropically consolidated under an effective pressure of 150 kPa.Table 1
Physical properties of the material.
Specific gravity:Gs
Liquid limit:wl (%)
Plastic limit:wp (%)
Plasticity index:Ip
Natural density:ρ (g/cm3)
Natural water content:w (%)
2.68
49.8
22.8
27
1.78
42.2


### 2.3. Test Plan

In this studyα was increased up to 360° to highlight the first two cycles (principal stress rotated 180° in each cycle), since the soil response to the pure principal stress rotation was mainly embodied in the first several cycles, as manifested in the documents [2, 6, 7]. Nine tests were conducted and divided into three series according to the magnitudes of q and b. After the accomplishment of consolidation, b and shear stress q gradually increased to the values listed in Table 2 under undrained condition and constant p=150 kPa to point A in Figure 2. The torque was applied with an angular displacement rate of 0.2°/min (corresponding to a shear stress rate of 2.8 × 10−3kPa/sec) when pore water pressure equalized. All the tests followed the same procedure except test RS5010-2, which was sheared from α=90° (point B in Figure 2) and stopped earlier due to the failure. The stress path of pure principal stress rotation in the stress plane of (σz-σθ)/2~τzθ is illustrated in Figure 2, the setting stress path can be carried out by the HCA with an acceptable accuracy.Table 2
Test program of pure principal stress rotation.
Test no.
q/kPa
b
α/°
Series 1
RS3700
37.5
0
0~360
RS3725
37.5
0.25
0~360
RS3705
37.5
0.5
0~360
RS3775
37.5
0.75
0~360
RS3710
37.5
1
0~360
Series 2
RS2510
25
1
0~360
RS3710
37.5
1
0~360
RS5010
50
1
0~360
RS5010-2
50
1
90~270
Series 3
RS2505
25
0.5
0~360
RS3705
37.5
0.5
0~360Figure 2
Stress path of pure principal stress rotation.Stress components need to be well designed to perform pure principal stress rotation test. The corresponding stress paths were realized by simultaneous change ofσz, σθ, and τzθ with α, and keeping constant radial stress σr (Figure 3(a)). according to the equations of (1). The variation of stress components of σz, σr, σθ, and τzθ during principal stress rotation with constant p, q, and b can be determined
(2)σz=p-2bq-q3+qcos2α,(3)σr=p+2(2bq-q)3,(4)σθ=p-2bq-q3-qcos2α,(5)τzθ=qsin2α.Variation in stress components (a) and principal stresses (b) of RS3705.
(a)
(b)Figure3 shows the variation of the stress components and the principal stresses with rotation in test RS3705. The variations of stress components of σz, σθ, and τzθ all follow the equations of (2)–(5). Stress responses in Figure 3 also demonstrate the reliability and accuracy of the apparatus.


### 2.4. Repeatability of Test

Two parallel tests on the reconstituted clay were carried out despite the fact that the transducers in the ZJU-HCA have very good sensitivity and high resolution. The reconstituted clay was made from slurry being  isotropically consolidated with a vacuum pressure of 65 kPa in a patented reconstitution device [20]. Two specimens R1 and R2 were prepared and consolidated at an effective pressure of 150 kPa. Subsequently undrained shear test was performed under b=0.5. First the shear stress was increased to q=25 kPa (AB in Figure 4(a)), then the major principal stress axes rotated from 0° to 90° (BC in Figure 4(a)), and finally the specimens were sheared to failure with α fixed at 90° (CD in Figure 4(a)). Two tests had identical stress path, as illustrated in Figure 4(a). During the whole test, the magnitude of p was kept constant at 150 kPa.(a) Stress path ofR1 and R2; (b) Δσ1-Δε1 curves of R1 and R2.
(a)
(b)The stress-strain curves of the two tests had good consistency, especially at small strain stage, as shown in Figure4(b). Some divergence was observed as specimens approaching failure at large strain. This result provided strong evidence of the acceptable repeatability of this apparatus, and the negligible effect of the specimen preparation and assembling.


## 2.1. Hollow Cylinder Apparatus

Hollow cylinder apparatuses have been successfully used in the soil anisotropy studies. Here, different types of principal stress rotation tests on the intact Hangzhou soft clay were conducted using an advanced Zhejiang University hollow cylinder apparatus (ZJU-HCA). In this HCA the external loadings of axial loadW, inner pressure pi, outer pressure po, and torque MT can be applied and controlled individually to simulate the complex stress paths. The big advantage of simultaneously applying dynamic axial load and torque with a frequency up to 5 Hz provides many opportunities for the complex test. The maximum values of axial load and torque of ZJU-HCA are 12 kN and 200 Nm.The ZJU-HCA can be used for tests on both intact and reconstituted samples with the size of 160 mm × 200 mm × 400 mm and 60 mm × 100 mm × 200 mm (inner diameter × outer diameter × height), respectively. The details of this apparatus and the sampling kit can refer to Shen [18].Corresponding to the axial load, inner cell pressure, outer cell pressure and torque, the stress components ofσz, σr, σθ, and τzθ on the specimen indicated in Figure 1(a) can be calculated according to Hight et al. [19]. These four stress components can be converted into an equivalent set of four stress-related parameters, p, q, α, and b as the following:
(1)p=(σ1+σ2+σ3)3=(σz+σr+σθ)3q=(σ1-σ3)2=(σz-σθ)2/4+τθz2α=12·tan-1(2τθzσz-σθ)b=σ2-σ3σ1-σ3=((2σr-σz-σθ)/2)+((σz-σθ)/2)2+τθz2(σz-σθ)2+4τθz2,
where p is the mean total stress, q is the shear stress, b is the intermediate principal stress parameter (0≤b≤1), and α is the angle between the directions of the major principal stress and the vertical axis (Figure 1(a)). In this study the principal stresses (σ1 and σ3) rotate in the plane perpendicular to σ2.(a) Stress state of hollow cylinder specimens; (b) a completed hollow cylindrical specimen.
(a)
(b)As the four stress components ofσz, σr, σθ, and τzθ can be controlled independently, stress paths with constants  p, q, and b while α increases continuously (pure principal stress rotation) can be achieved [18]. The main objective of this study is to investigate the influence of b on the responses of intact soft clay under the pure principal stress rotation. This study can be well carried out with the support of ZJU-HCA.


## 2.2. Test Material and Sample Preparation

The hand-cut clay blocks were taken from an excavation pit 4 m below the ground level in Hangzhou; the basic physical properties of the clay are summarized in Table1. After the hollow cylindrical specimen (Figure 1(b)) with the size of 100 mm in outer diameter, 60 mm in inner diameter and 200 mm in height was prepared by a special sampling kit [18]; filter paper strips were pasted around the outer surface of the specimen to equalize pore water pressure and reduce the consolidation time. After that, the specimen was assembled into the cell to be saturated with a backpressure of 100 kPa, ensuring the Skempton’s pore pressure parameter B higher than 0.98. After saturation the specimens were isotropically consolidated under an effective pressure of 150 kPa.Table 1
Physical properties of the material.
Specific gravity:Gs
Liquid limit:wl (%)
Plastic limit:wp (%)
Plasticity index:Ip
Natural density:ρ (g/cm3)
Natural water content:w (%)
2.68
49.8
22.8
27
1.78
42.2


## 2.3. Test Plan

In this studyα was increased up to 360° to highlight the first two cycles (principal stress rotated 180° in each cycle), since the soil response to the pure principal stress rotation was mainly embodied in the first several cycles, as manifested in the documents [2, 6, 7]. Nine tests were conducted and divided into three series according to the magnitudes of q and b. After the accomplishment of consolidation, b and shear stress q gradually increased to the values listed in Table 2 under undrained condition and constant p=150 kPa to point A in Figure 2. The torque was applied with an angular displacement rate of 0.2°/min (corresponding to a shear stress rate of 2.8 × 10−3kPa/sec) when pore water pressure equalized. All the tests followed the same procedure except test RS5010-2, which was sheared from α=90° (point B in Figure 2) and stopped earlier due to the failure. The stress path of pure principal stress rotation in the stress plane of (σz-σθ)/2~τzθ is illustrated in Figure 2, the setting stress path can be carried out by the HCA with an acceptable accuracy.Table 2
Test program of pure principal stress rotation.
Test no.
q/kPa
b
α/°
Series 1
RS3700
37.5
0
0~360
RS3725
37.5
0.25
0~360
RS3705
37.5
0.5
0~360
RS3775
37.5
0.75
0~360
RS3710
37.5
1
0~360
Series 2
RS2510
25
1
0~360
RS3710
37.5
1
0~360
RS5010
50
1
0~360
RS5010-2
50
1
90~270
Series 3
RS2505
25
0.5
0~360
RS3705
37.5
0.5
0~360Figure 2
Stress path of pure principal stress rotation.Stress components need to be well designed to perform pure principal stress rotation test. The corresponding stress paths were realized by simultaneous change ofσz, σθ, and τzθ with α, and keeping constant radial stress σr (Figure 3(a)). according to the equations of (1). The variation of stress components of σz, σr, σθ, and τzθ during principal stress rotation with constant p, q, and b can be determined
(2)σz=p-2bq-q3+qcos2α,(3)σr=p+2(2bq-q)3,(4)σθ=p-2bq-q3-qcos2α,(5)τzθ=qsin2α.Variation in stress components (a) and principal stresses (b) of RS3705.
(a)
(b)Figure3 shows the variation of the stress components and the principal stresses with rotation in test RS3705. The variations of stress components of σz, σθ, and τzθ all follow the equations of (2)–(5). Stress responses in Figure 3 also demonstrate the reliability and accuracy of the apparatus.


## 2.4. Repeatability of Test

Two parallel tests on the reconstituted clay were carried out despite the fact that the transducers in the ZJU-HCA have very good sensitivity and high resolution. The reconstituted clay was made from slurry being  isotropically consolidated with a vacuum pressure of 65 kPa in a patented reconstitution device [20]. Two specimens R1 and R2 were prepared and consolidated at an effective pressure of 150 kPa. Subsequently undrained shear test was performed under b=0.5. First the shear stress was increased to q=25 kPa (AB in Figure 4(a)), then the major principal stress axes rotated from 0° to 90° (BC in Figure 4(a)), and finally the specimens were sheared to failure with α fixed at 90° (CD in Figure 4(a)). Two tests had identical stress path, as illustrated in Figure 4(a). During the whole test, the magnitude of p was kept constant at 150 kPa.(a) Stress path ofR1 and R2; (b) Δσ1-Δε1 curves of R1 and R2.
(a)
(b)The stress-strain curves of the two tests had good consistency, especially at small strain stage, as shown in Figure4(b). Some divergence was observed as specimens approaching failure at large strain. This result provided strong evidence of the acceptable repeatability of this apparatus, and the negligible effect of the specimen preparation and assembling.


## 3. Test Results and Discussion

### 3.1. Evolution of Pore Water Pressure

Although the shear stressq was constant in pure rotation tests, significant pore water pressure (PWP for short) accumulated with the principal stress rotation (Figure 5), which surely demonstrated the effect of the intermediate principal stress. Clear fluctuation was observed in PWP, which was in agreement with the test results from sand [2, 6]. The fluctuation was not significant in the first cycle as compared with the second. In the first cycle PWP accumulation was approximately double that of the second cycle. Table 3 shows PWP increments during α=0°~90° and α=90°~180° in each cycle of series 1. Most PWP accumulation generated during the stage when α increased from 90° to 180° and PWP reached a peak value at the end of each cycle (α=180°).Table 3
PWP increment of series 1.
PWP incrementΔu/kPa
Test no.
Cycles
α
=
0
°
–
90
°
α
=
90
°
–
180
°
∑
RS3700
Δ
u
cycle
1
3.03
28.96
31.99
Δ
u
cycle
2
−10.77
25.61
14.84
RS3725
Δ
u
cycle
1
12.7
32.78
45.48
Δ
u
cycle
2
−5.49
27.49
22
RS3750
Δ
u
cycle
1
6.33
33.04
39.37
Δ
u
cycle
2
−13.28
33.09
19.81
RS3775
Δ
u
cycle
1
14.28
32.73
47.01
Δ
u
cycle
2
−1.15
29.49
28.34
RS3710
Δ
u
cycle
1
9.46
34.82
44.28
Δ
u
cycle
2
−13.75
35.13
21.38Pore water pressure accumulation withα: (a) series 1; (b) series 2; (c) series 3; (d) R1, R2, and RS2505.
(a)
(b)
(c)
(d)The effect ofb values on PWP accumulation with α in the series 1 is presented in Figure 5(a). As it shows, the mode of PWP accumulation with α under different b was similar. However, the rate of accumulation was affected significantly by b. PWP accumulation under the conditions  b=0.75 and b=0.25 were larger than that under the conditions b=0,0.5 and 1.0. And specimen in the test under b=0 had the lowest accumulation of pore water pressure; this observation is also consistent with the test results from sand [2]. The fact that the PWP not monotonously increasing with b suggests that the influence of b is more complex; it could not be fully revealed if only a limited range of b values were considered.The effect of shear stress on PWP is much pronounced, as illustrated in Figures5(b) and 5(c). PWP increased rapidly with large shear stress. The difference in PWP accumulation between the tests at q=50 kPa and q=37.5 kPa was smaller than that between the tests at q=25 kPa and q=37.5 kPa in series 2 under the condition b=1.0. Undrained rotational shear tests on sand [2] had same observation that the influence of shear stress on PWP reduces when the value of shear stress increased.Test RS5010-2 had the same stress path as test RS5010 except that principal stress rotated fromα=90° instead of α=0°. It is interesting to note that the peak value of PWP in test RS5010-2 had nearly the same magnitude and occurred at the same direction (at α=180°) as the test RS5010. Moreover, the PWP generation in test RS5010 and test RS5010-2 were still very close when α=180°~250°, which clearly shows that PWP build-up during rotational shear is less dependent on the initial direction of principal stress. The PWP accumulation in RS5010-2 also supports the observations that PWP generated mainly during the period of α=90°~180°.Whenα=250° the PWP of RS5010-2 increased abruptly and the specimen failed with large deformation at α=270°. The PWP value at failure point was far smaller than the initial confining pressure (150 kPa), which implies that the failure in undrained pure principal stress rotation test is mainly controlled by strain, not PWP accumulation.PWP accumulation in reconstituted sample was quite different from the intact sample. Figure6(d) shows the comparison of PWP accumulation in tests R1,R2 (reconstituted), and RS2505 (intact). These tests started from the same stress state (b=0.5,q=25 kPa) and experienced the same rotation stress path (from 0° to 90°), but PWP in intact clay had higher developing rate and larger accumulation than the reconstituted ones. This is attributed to the inherent anisotropy of the intact sample. Here, it has to note that even in reconstituted specimens, which is normally treated as isotropic material [9], PWP still accumulated and was approximately 40% of the intact one. This phenomenon is very important and it indicates that even the changes in the direction of principal stress could lead to pore water pressure accumulation, which further causes soil degradation. This potential risk cannot be underestimated in practical design.Variation in strain components withα in series 1: (a) axial strain; (b) circumferential strain; (c) radial strain; (d) torsional shear strain.
(a)
(b)
(c)
(d)


### 3.2. Variation of Strain Components

Pure principal stress rotation test was performed under constantsp, q, and b; however strain development was still significant. Figure 6 presents strain development in series 1 test (take compression as positive). Strain components of εz,εθ, and γzθ had cyclic variation with α however the variation of stress and strain components was not synchronous. Development of εz, εθ, and γzθ lagged 20° behind the corresponding stresses. For example, the torsional shear stress τzθ reached the peak value at α=45° (Figure 3(a)), but γzθ approached peak value until α=65°. That was the evidence of the noncoaxiality between the directions of the principal stress and the principal strain increment also found in sand [4, 5, 7] and will be discussed later. It is also clear that strain accumulation in the second cycle was larger than that in the first cycle, which was different from the PWP evolution.The variation of strain components (εz, εθ, γzθ, and εr) with αwere very similar under smaller b, such as the cases b=0 and b=0.25. The effect of b increased with the magnitude, as clearly shown in εz in Figure 6(a). As b increased, that is the influence of radial stress σr increased (see, (3)), radial strain increased. Due to larger radial stress the influence of other stress on the corresponding strain decreased, and soil samples exhibited less compression, but more dilation. Similar situation also happened in εθ and γzθ (Figures 6(b) and 6(d)). Compared with εz and εθ, the impact of b on γzθ was relatively small (shown in Figure 6(d)).Other than the oscillation inεz, εθ and γzθ, approximately a linear relationship between εr and α was observed, which supported the fact that the variation of α was independent of σr. The increasing εr was induced by plastic deformation. In Figures 6(a), 6(b), and 6(d) irrecoverable strain at the end of each cycle was also noticed, which meant plastic deformation could be generated by the changes of principal stress direction alone. Results in this paper also indicate that the rate of strain development in clayey soil was much higher than sandy material [2, 5, 6] when soil samples experienced identical stress path and similar stress level.It needs to mention that whenb=0.5 very small radial strain developed during the whole rotation (seen in Figure 6(c)). This indicates that the assumption of plane strain condition is reasonable when b=0.5.Shear stress-strain relationship of series 1 is given in Figure7. The unclosed hysteretic loops manifest the accumulation of plastic strain. Shear stiffness G, defined as the secant modulus of shear stress-strain loops [2], also affected by b, as shown in Figure 7. Though specimens only experienced two cycles of principal stress rotation, significant stiffness degradation can be observed, and degradation was more significant as b increased. This observation agrees well with the tests on sand [2].Figure 7
Shear stress-strain relationship of series 1.The influence of shear stress on strain components is much larger thanb value. In series 2 the shear stress in test RS5010 was two times of test RS2510, but the magnitude of strain components was ten times larger, as seen in Figures 8(a) and 8(b). Horizontal symmetry was found in εz and εθ in the condition of low shear stress and small rotation range. Same phenomenon was observed in drained tests on Toyoura dense sand (Miura et al. [5]) and undrained tests on intact clay (Zhou and Xu [21]). However the behavior of symmetry was affected as α increased. With the rapid development of strain in high shear stress test, this phenomenon was not obvious, such as test RS5010-2. Specimen of RS5010-2 experienced identical stress path as RS5010 except the initial direction of principal stress. The different responses in strain components illustrate that the initial principal stress direction had significant influence on the deformation behavior of intact clay under pure principal stress rotation. Compared with test RS5010, RS5010-2 approached earlier failure and the radial and axial strain increased dramatically at the onset of failure (Figure 8(c)). This strain development is consistent with the development of pore pressure as given in Figure 5(b). However, the effect of initial rotation direction on the following behavior has not caused attention.Variation of strain components withα in series 2: (a) RS2510; (b) RS5010; (c) RS5010-2.
(a)
(b)
(c)


### 3.3. Noncoaxial Behavior of Intact Clay

As was observed in lots of experimental researches, the direction of the principal strain increment has no coincidence with the principal stress direction during the principal stress rotation [5, 7, 22–25]. The deviation between them is defined as the noncoaxial behavior of soils [23, 24], which can be described by the noncoaxial angle of β=(αdε-α). Where αdε is the direction of the principal strain increment from the vertical axis of the specimen, it can be expressed as
(6)αdε=12arctan(dγzθdεz-dεθ),
where dεz,dεθ, and dγzθ are the increment of axial strain, circumferential strain and torsional shear strain, respectively.The variation of noncoaxial angle in series 1 is presented in Figure9. It can be found that the noncoaxial behavior of intact clay was obvious in the pure principal rotation test. The variation of β with α in both cycles was similar (Figures 9(a) and 9(b)), and the noncoaxiality was less affected by the number of rotation cycles. Nevertheless there was some fluctuation in β. β decreased with α when α rotated from 0° to 40° and from 90° to 120°. Tong et al. [7] had the same observation from their drained test on sand. When α rotated from 45° to 80° and from 120° to 180°, results were very different; that is, β increased with α in this research, while β kept constant (β = 45°); in their tests. In general, it was found that noncoaxiality existed along the whole test procedure with clayey soil. Despite of the fluctuation between 22° and 30°, the average β value can be taken as 27°.Variation ofβ=(αdε-α) with α of series 1: (a) First cycle; (b) Second cycle.
(a)
(b)Experimental study in this research also shows that the noncoaxiality of intact clay was influenced by the shear stress rather than the intermediate principal stress (illustrated in Figure10). In series 2, β also fluctuated with α, but the amplitude was mainly controlled by the shear stress. The noncoaxial angles of specimens in the cases of q=37.5 kPa and q=50 kPa were very close in both cycles. However, the noncoaxial angles in the case of q=25 kPa were dramatic, this is consistent with the test result from sand by Cai [24]. As pointed out by Miura et al. [5], Pradel et al. [26], and Ohkawa et al. [27] that the direction of the principal strain increment was strongly dependent on the stress path (direction of stress increment), this study also concluded that noncoaxiality is very significant in complex stress path test and the associated flow rule is no longer applicable.Variation ofβ=(αdε-α) with α of series 2: (a) first cycle; (b) second cycle.
(a)
(b)


## 3.1. Evolution of Pore Water Pressure

Although the shear stressq was constant in pure rotation tests, significant pore water pressure (PWP for short) accumulated with the principal stress rotation (Figure 5), which surely demonstrated the effect of the intermediate principal stress. Clear fluctuation was observed in PWP, which was in agreement with the test results from sand [2, 6]. The fluctuation was not significant in the first cycle as compared with the second. In the first cycle PWP accumulation was approximately double that of the second cycle. Table 3 shows PWP increments during α=0°~90° and α=90°~180° in each cycle of series 1. Most PWP accumulation generated during the stage when α increased from 90° to 180° and PWP reached a peak value at the end of each cycle (α=180°).Table 3
PWP increment of series 1.
PWP incrementΔu/kPa
Test no.
Cycles
α
=
0
°
–
90
°
α
=
90
°
–
180
°
∑
RS3700
Δ
u
cycle
1
3.03
28.96
31.99
Δ
u
cycle
2
−10.77
25.61
14.84
RS3725
Δ
u
cycle
1
12.7
32.78
45.48
Δ
u
cycle
2
−5.49
27.49
22
RS3750
Δ
u
cycle
1
6.33
33.04
39.37
Δ
u
cycle
2
−13.28
33.09
19.81
RS3775
Δ
u
cycle
1
14.28
32.73
47.01
Δ
u
cycle
2
−1.15
29.49
28.34
RS3710
Δ
u
cycle
1
9.46
34.82
44.28
Δ
u
cycle
2
−13.75
35.13
21.38Pore water pressure accumulation withα: (a) series 1; (b) series 2; (c) series 3; (d) R1, R2, and RS2505.
(a)
(b)
(c)
(d)The effect ofb values on PWP accumulation with α in the series 1 is presented in Figure 5(a). As it shows, the mode of PWP accumulation with α under different b was similar. However, the rate of accumulation was affected significantly by b. PWP accumulation under the conditions  b=0.75 and b=0.25 were larger than that under the conditions b=0,0.5 and 1.0. And specimen in the test under b=0 had the lowest accumulation of pore water pressure; this observation is also consistent with the test results from sand [2]. The fact that the PWP not monotonously increasing with b suggests that the influence of b is more complex; it could not be fully revealed if only a limited range of b values were considered.The effect of shear stress on PWP is much pronounced, as illustrated in Figures5(b) and 5(c). PWP increased rapidly with large shear stress. The difference in PWP accumulation between the tests at q=50 kPa and q=37.5 kPa was smaller than that between the tests at q=25 kPa and q=37.5 kPa in series 2 under the condition b=1.0. Undrained rotational shear tests on sand [2] had same observation that the influence of shear stress on PWP reduces when the value of shear stress increased.Test RS5010-2 had the same stress path as test RS5010 except that principal stress rotated fromα=90° instead of α=0°. It is interesting to note that the peak value of PWP in test RS5010-2 had nearly the same magnitude and occurred at the same direction (at α=180°) as the test RS5010. Moreover, the PWP generation in test RS5010 and test RS5010-2 were still very close when α=180°~250°, which clearly shows that PWP build-up during rotational shear is less dependent on the initial direction of principal stress. The PWP accumulation in RS5010-2 also supports the observations that PWP generated mainly during the period of α=90°~180°.Whenα=250° the PWP of RS5010-2 increased abruptly and the specimen failed with large deformation at α=270°. The PWP value at failure point was far smaller than the initial confining pressure (150 kPa), which implies that the failure in undrained pure principal stress rotation test is mainly controlled by strain, not PWP accumulation.PWP accumulation in reconstituted sample was quite different from the intact sample. Figure6(d) shows the comparison of PWP accumulation in tests R1,R2 (reconstituted), and RS2505 (intact). These tests started from the same stress state (b=0.5,q=25 kPa) and experienced the same rotation stress path (from 0° to 90°), but PWP in intact clay had higher developing rate and larger accumulation than the reconstituted ones. This is attributed to the inherent anisotropy of the intact sample. Here, it has to note that even in reconstituted specimens, which is normally treated as isotropic material [9], PWP still accumulated and was approximately 40% of the intact one. This phenomenon is very important and it indicates that even the changes in the direction of principal stress could lead to pore water pressure accumulation, which further causes soil degradation. This potential risk cannot be underestimated in practical design.Variation in strain components withα in series 1: (a) axial strain; (b) circumferential strain; (c) radial strain; (d) torsional shear strain.
(a)
(b)
(c)
(d)


## 3.2. Variation of Strain Components

Pure principal stress rotation test was performed under constantsp, q, and b; however strain development was still significant. Figure 6 presents strain development in series 1 test (take compression as positive). Strain components of εz,εθ, and γzθ had cyclic variation with α however the variation of stress and strain components was not synchronous. Development of εz, εθ, and γzθ lagged 20° behind the corresponding stresses. For example, the torsional shear stress τzθ reached the peak value at α=45° (Figure 3(a)), but γzθ approached peak value until α=65°. That was the evidence of the noncoaxiality between the directions of the principal stress and the principal strain increment also found in sand [4, 5, 7] and will be discussed later. It is also clear that strain accumulation in the second cycle was larger than that in the first cycle, which was different from the PWP evolution.The variation of strain components (εz, εθ, γzθ, and εr) with αwere very similar under smaller b, such as the cases b=0 and b=0.25. The effect of b increased with the magnitude, as clearly shown in εz in Figure 6(a). As b increased, that is the influence of radial stress σr increased (see, (3)), radial strain increased. Due to larger radial stress the influence of other stress on the corresponding strain decreased, and soil samples exhibited less compression, but more dilation. Similar situation also happened in εθ and γzθ (Figures 6(b) and 6(d)). Compared with εz and εθ, the impact of b on γzθ was relatively small (shown in Figure 6(d)).Other than the oscillation inεz, εθ and γzθ, approximately a linear relationship between εr and α was observed, which supported the fact that the variation of α was independent of σr. The increasing εr was induced by plastic deformation. In Figures 6(a), 6(b), and 6(d) irrecoverable strain at the end of each cycle was also noticed, which meant plastic deformation could be generated by the changes of principal stress direction alone. Results in this paper also indicate that the rate of strain development in clayey soil was much higher than sandy material [2, 5, 6] when soil samples experienced identical stress path and similar stress level.It needs to mention that whenb=0.5 very small radial strain developed during the whole rotation (seen in Figure 6(c)). This indicates that the assumption of plane strain condition is reasonable when b=0.5.Shear stress-strain relationship of series 1 is given in Figure7. The unclosed hysteretic loops manifest the accumulation of plastic strain. Shear stiffness G, defined as the secant modulus of shear stress-strain loops [2], also affected by b, as shown in Figure 7. Though specimens only experienced two cycles of principal stress rotation, significant stiffness degradation can be observed, and degradation was more significant as b increased. This observation agrees well with the tests on sand [2].Figure 7
Shear stress-strain relationship of series 1.The influence of shear stress on strain components is much larger thanb value. In series 2 the shear stress in test RS5010 was two times of test RS2510, but the magnitude of strain components was ten times larger, as seen in Figures 8(a) and 8(b). Horizontal symmetry was found in εz and εθ in the condition of low shear stress and small rotation range. Same phenomenon was observed in drained tests on Toyoura dense sand (Miura et al. [5]) and undrained tests on intact clay (Zhou and Xu [21]). However the behavior of symmetry was affected as α increased. With the rapid development of strain in high shear stress test, this phenomenon was not obvious, such as test RS5010-2. Specimen of RS5010-2 experienced identical stress path as RS5010 except the initial direction of principal stress. The different responses in strain components illustrate that the initial principal stress direction had significant influence on the deformation behavior of intact clay under pure principal stress rotation. Compared with test RS5010, RS5010-2 approached earlier failure and the radial and axial strain increased dramatically at the onset of failure (Figure 8(c)). This strain development is consistent with the development of pore pressure as given in Figure 5(b). However, the effect of initial rotation direction on the following behavior has not caused attention.Variation of strain components withα in series 2: (a) RS2510; (b) RS5010; (c) RS5010-2.
(a)
(b)
(c)


## 3.3. Noncoaxial Behavior of Intact Clay

As was observed in lots of experimental researches, the direction of the principal strain increment has no coincidence with the principal stress direction during the principal stress rotation [5, 7, 22–25]. The deviation between them is defined as the noncoaxial behavior of soils [23, 24], which can be described by the noncoaxial angle of β=(αdε-α). Where αdε is the direction of the principal strain increment from the vertical axis of the specimen, it can be expressed as
(6)αdε=12arctan(dγzθdεz-dεθ),
where dεz,dεθ, and dγzθ are the increment of axial strain, circumferential strain and torsional shear strain, respectively.The variation of noncoaxial angle in series 1 is presented in Figure9. It can be found that the noncoaxial behavior of intact clay was obvious in the pure principal rotation test. The variation of β with α in both cycles was similar (Figures 9(a) and 9(b)), and the noncoaxiality was less affected by the number of rotation cycles. Nevertheless there was some fluctuation in β. β decreased with α when α rotated from 0° to 40° and from 90° to 120°. Tong et al. [7] had the same observation from their drained test on sand. When α rotated from 45° to 80° and from 120° to 180°, results were very different; that is, β increased with α in this research, while β kept constant (β = 45°); in their tests. In general, it was found that noncoaxiality existed along the whole test procedure with clayey soil. Despite of the fluctuation between 22° and 30°, the average β value can be taken as 27°.Variation ofβ=(αdε-α) with α of series 1: (a) First cycle; (b) Second cycle.
(a)
(b)Experimental study in this research also shows that the noncoaxiality of intact clay was influenced by the shear stress rather than the intermediate principal stress (illustrated in Figure10). In series 2, β also fluctuated with α, but the amplitude was mainly controlled by the shear stress. The noncoaxial angles of specimens in the cases of q=37.5 kPa and q=50 kPa were very close in both cycles. However, the noncoaxial angles in the case of q=25 kPa were dramatic, this is consistent with the test result from sand by Cai [24]. As pointed out by Miura et al. [5], Pradel et al. [26], and Ohkawa et al. [27] that the direction of the principal strain increment was strongly dependent on the stress path (direction of stress increment), this study also concluded that noncoaxiality is very significant in complex stress path test and the associated flow rule is no longer applicable.Variation ofβ=(αdε-α) with α of series 2: (a) first cycle; (b) second cycle.
(a)
(b)


## 4. Conclusions

A series of pure principal stress rotation tests on intact soft clay were conducted by using the ZJU-5 Hz hollow cylinder apparatus in this study. The magnitude of mean total stressp, intermediate principal stress parameter b and shear stress q were all maintained constant during rotation. The influence of b on the response of pore water pressure and deformation behavior was focused and the main conclusions can be drawn as follows.(1)
Significant pore water pressure accumulation was observed in pure principal stress rotation test. PWP accumulated mainly whenα increased from 90° to 180° and most of the accumulation generated in the first cycle.(2)
The rate of PWP evolution was greatly influenced byb, but the influence was not monotonous. Specimens under the condition b=0 had the lowest PWP accumulation, while the highest accumulation occurred under the condition b=0.75. PWP accumulation was less influenced by the previous stress path.(3)
Development of strain components induced by pure principal stress rotation was affected by the magnitudes ofb and q. Specimen can also fail in pure principal stress rotation test due to large deformation. Shear stiffness decreased with principal stress rotation, and it decreased quickly with large b.(4)
The deviation between the principal stress and the principal strain increment directions existed in the course of the test. The strain increment direction was strongly dependent on the direction of stress increment. Noncoaxial angles had some fluctuation and the average angle was 27° in this research. Noncoaxial character of intact soft clay under pure principal stress rotation was less influenced byb and the number of cycles.



---

*Source: 950143-2013-11-25.xml*