Abstract
The stress-dependent K0, Vs, and Vs anisotropy and their correlations with sand for 1D consolidation stress were tested with a custom-designed floating-wall consolidometer-type Bender Element (BE) testing apparatus. K0 of a soil sample was calculated using stress measurements through soil pressure transducers installed at the midsection of the consolidometer. The Vs and Vs anisotropy were measured by the bender elements installed in three orthogonal directions in the consolidometer, i.e., vh, hv, and hh. Granular soils with different sizes and shapes were tested. The effects of the stress level, overconsolidation ratio (OCR), particle size and shape on the Vs anisotropy, and K0 of the granular soils during one-dimensional consolidation were investigated. The laboratory investigations suggested (1) the K0 showed a constant value during loading, while it increased as the OCR increased during unloading, (2) soils with smaller particle sizes, rough surfaces, and angular geometry tended to have a lower value of K0, and vice versa, (3) both the anisotropic stress state and the anisotropic fabric (geometry) could lead to the Vs anisotropy, but the Vs anisotropy was manifested due to the horizontal stress-lock during unloading stage, and (4) the published correlation between Vs and K0 was modified by introducing the influence of the OCR, which could effectively reduce the variation and improve the prediction accuracy. Therefore, the modified correlation could be used as a robust approach to estimate K0 for both normally consolidated and highly overly consolidated granular soils.
1. Introduction
K0 (equation (1)) is an intrinsic soil property that is related to soil characteristics, stress state, and stress history [1–5]. Therefore, it is used frequently in engineering practices such as retaining walls, embankment construction, mining and tunneling, and deep excavation:
Based on the stress distribution analysis within wedge-shaped sand dunes, Jaky [6] reported that the value of K0 was related to the internal friction angle of the sand. For normally consolidated soil, i.e., clay or loose sand, Jaky [6] derived a theoretical formula to calculate K0:where is the effective critical-state friction angle of the granular soil. There are theoretical flaws with predicting the lateral Earth pressure at rest in this equation because the equation indicates that K0 is related to the frictional force of the soil, while the stress state at rest is stress path dependent. Michalowski [4] revisited the original derivation of the equation of Jaky [6] and concluded that Jaky’s derivation was stress path dependent. It was a 1D strain condition that was associated with the K0 state. This finding indicates that K0 at the symmetry axis of a sand pile covers the entire range from the active stress state to the passive stress state. However, the stress state (the distribution of the active and passive stresses) in the 1D strain process may not be the same as that for a wedge-shaped sand pile [4]. In addition, K0 was also found to be related to the stress history (i.e., OCR) and physical properties of geomaterials [7]. For instance, preloading can induce the residual horizontal stresses to be locked, which will eventually lead to an increased K0. Therefore, in order to predict K0 during unloading and reloading, the overconsolidation ratio (OCR) should be taken into consideration. For this reason, Mayne and Kulhawy [8] made modifications to equation (2) that accounted for the effect of overconsolidation. However, this equation is still empirical in nature:where is the exponential parameter that defines the effect of the OCR on K0 and OCR is the overconsolidation ratio of soil, which is usually approximated by using . This modified correlation is good for overconsolidated soils during expansion but not for recompression. Therefore, for soils during unloading and reloading cycles, Mayne and Kulhawy [8] proposed another empirical equation:where can be interpreted with Jaky’s equation and is the maximum OCR of the soil in its stress history. Equations (1)–(3) indicated that K0 of a soil would stay quite constant during loading [8]. During the unloading process, K0 increases with the OCR [1]. Although the proposed estimations are widely used to estimate K0, there are still limitations associated with these previous published classical theoretical equations. For Jaky’s equation, the special stress state in the wedge-shaped sand pile is different from that in the one-dimensional consolidation, thus inhibiting its application and producing the inaccurate estimation of K0. Equations (2) and (3) could offer reasonable estimation of K0; however, they are empirical in nature. Therefore, the prediction of K0 still needs further investigation.
Except for the theoretical equations, K0 can be measured in a laboratory, i.e., using specially equipped oedometers [1, 3, 6, 9–17], a triaxial apparatus [18–22], a “soft” thin wall oedometer [2], and a plane-strain apparatus [23]. However, all these K0 test apparatuses have their own limitations, i.e., high soil-wall resistance and the so-called “zero” radial strain for K0 consolidation. K0 can also be obtained from an in situ test. However, in situ tests will inevitably cause soil displacement issues during penetration, and they may significantly influence the actual K0 state of stress. Kang et al. [24] manufactured a floating-wall oedometer bender element testing apparatus that overcame the limitations of the previous K0 devices. The testing system [24], which was integrated with stress sensors and the measurements of the shear wave velocity, claimed a truly nonintrusive approach. In addition, the use of Vs to predict the in situ horizontal stress, and K0 has also been reported by previous studies [5, 25]. Ku and Mayne [5] used different types of shear wave velocities to estimate K0. Once the Vs profiles (either Vs−hh and Vs−hv or Vs−hh and Vs−vh) are obtained, K0 may be calculated with either of the following equations:
Through a literature review and field measurements, Ku and Mayne [5] suggested that the “hh-vh” mode yielded good prediction of K0. The proposed fitting constants and were given as 0.78 and 2.55, respectively:
Soils often have strong anisotropy properties [26–33]. Therefore, Vs depends on the structure of the soil fabric and the measuring directions, such as longitudinal or transverse [24, 34–37]. By comparing Vs in three orthogonal directions, i.e., Vs−vh, Vs−hv, and Vs−hh, the Vs anisotropy of soils can be evaluated, where the first subscript letter represents the direction of the shear wave propagation direction and the second subscript letter represents the direction of the soil particle vibration ( = vertical and h = horizontal). The Vs−vh, Vs−hv, and Vs−hh are compared during the loading and unloading stages. The stiffness anisotropy is solely related to Vs because . Since there were three pairs of BE sets installed in three different directions inside a sample that could measure Vs in three directions, the term “Vs anisotropy” is considered the same as the term “small-strain stiffness anisotropy” in this paper. The measurement of Vs also offers a different way to investigate the anisotropy of geomaterials, which is related to the applied stress, the inherent soil fabric, and the physicochemical conditions [24]. Therefore, K0 of granular materials was experimentally investigated in this research, which was focused on the correlations with the shear wave velocity, stress levels, and different particle sizes and shapes. The main aims of this study were to (1) study the K0 and Vs characteristics of granular materials during one-dimensional consolidation, (2) evaluate the current correlations between K0 and the shear wave velocity of granular materials, and (3) characterize the Vs anisotropy of granular materials.
1.1. Test Setup
A custom-designed floating consolidometer bender element testing system (Figure 1(a)) was used in this study. The testing system consisted of a water chamber, a floating-wall oedometer, two soil pressure sensors (one sensor installed at the base, and the other sensor located at the center of the cell), three BEs, and the wave generation and processing apparatus. A piece of bridge-shaped steel was used to cover the benders and the soil pressure sensors as well as to apply the vertical load. The floating wall worked as a constrained ring so that the K0 condition was guaranteed. The horizontal stress was measured through a horizontally installed pressure sensor. A Teflon coating was put in the inner wall to reduce the soil-wall interface friction [38, 39]. A fluid chamber was used in which the floating wall was submerged in water during the test operation. This chamber maintained the saturation of the sample and controlled the physicochemical conditions of the bulk fluid, including the pH, temperature, and salt concentration. An Agilent signal generator 33210A was used to generate the source waves. A Krohn-Hite filter/signal conditioner 3364 was used to filter the received signals. All the received signals were read from the oscilloscope (54622A, Agilent). A single-cycle sine wave with a frequency ranging from 2 kHz to 10 kHz and an amplitude of 10 V was adopted. A BED-A soil pressure sensor was purchased from Civil Engineering and Construction Instruments, Kyowa Electronic Instruments Co., Ltd. The sensor had an outer diameter of 20 mm, and the sensing surface had a diameter of 27 mm. Its maximum capacity was 1 MPa. The BED-A sensors were able to work at a temperature range from −10 to 60°C and had a light weight, a good waterproof shield, and very good test sensitivity. All the soil pressure sensors were calibrated against a known load cell, and the test data was automatically read from the data logger. The floating-wall consolidometer design had the advantage of reducing the wall friction and protecting the side bender elements from detrimental bending. The weight of the oedometer could be counterbalanced through the built-in frictionless pulley system [24]. The bender element (Piezo Systems, Inc.) used in this study was cut into 12.7 mm × 8.0 mm × 0.6 mm (length × width × thickness) plates to fit the floating consolidometer bender element testing system. After trimming, polyurethane and silver powder were painted on the bender element surface to produce the physical and electrical shield. The first arrival time (t) was manually picked at the first zero-crossing of the first major peak by the cursor function in MathCAD (PTC, Needham, MA), and the tip-to-tip distance (l) was considered to be the shear wave propagation distance, so the shear wave velocity (Vs) could be calculated with Vs = l/t.
(a)
(b)
The schematic view of the floating-wall cell, soil pressure sensors, and BE sets is shown in Figure 1(a), and the detailed schematic setup of the floating-wall consolidometer cell is shown in Figure 1(b). The cell had a diameter of 11.4 cm, and its thickness was 1.3 cm. The height of the wall was 12.7 cm. Shear waves in different vibration directions could be generated in a soil specimen by the paired BE sets: Vs−ij, where “i” is the propagation direction and “j” is the polarization direction. One BE was installed in the center of the top platen vertically (represents vh direction), and two BEs were installed at the base platen with an alignment perpendicular to each other. Four horizontal BEs were installed in the side wall (Figure 1). One soil pressure transducer was installed in the center of the consolidometer cell with a screw plate fixed to the outside of the wall. Another soil sensor was located at the center of the base plate, right between the two bottom bender elements.
2. Materials and Methods
Glass beads (GB) with diameters of 0.5 mm, 3 mm, and 6 mm, Ottawa sands (20∼30, 50∼70, and 200∼325), and silica sand 30∼50 were used in this study to investigate the K0 and Vs anisotropy. The grain size distribution curves of all these materials (ASTM D 422-63) are shown in Figure 2. emax and emin are shown in Table 1 based on ASTM D 4253 and ASTM D 4254. The sand and the GB are classified as poorly graded sands, SP (ASTM D 2487). Sample preparation was performed using an air purification method [40, 41]. Dry specimens with a diameter of 114.3 mm (4.5 in) and a height of approximately 50.8 mm (2 in) were prepared in the consolidometer. A fixed falling distance of 80 cm from the funnel to the cell was predetermined to achieve approximately the same initial void ratio (relative density) for all the samples. A target mass of the soil was weighed and uniformly poured into the funnel. During each lift, a rubber hammer was used to tap the side of the cell gently in order to avoid any segregation and to internally build voids during falling. After setting the sample in the consolidometer, a steel tamper was used to carefully compact the top of the sand specimen and flatten the soil surface. Before putting on the top cap, measurements of the sample height were taken, and the average height was used to calculate the initial void ratio. A seating load of 16 kPa was then applied to push the top bender element into the soil and to maintain good contact between the soil and the sensors [42].
3. Results
One-dimensional consolidation tests were performed in this research. At the end of each loading stage, the base platen soil pressure and lateral Earth pressure were measured. At the end of consolidation, BE tests were carried out, during which the shear wave velocities at three orthogonal directions were measured and used to correlate with the K0. Each test material was subjected to triaxial compression tests to obtain the critical-state friction angles. For the specimens in the oedometer tests and CD triaxial compression tests, identical void ratios were maintained in order to avoid any differences that might have been induced from the initial density. The soil-wall interface friction coefficient, emax, emin, and e and the critical-state friction angle for each type of soil are shown in Table 1.
3.1. K0
In this study, the method proposed by Kang et al. [24] was applied to compute the effective normal stress, , of the specimens. By using the measured horizontal stress (), we calculated K0. K0 values of different materials are shown in Figure 3. All of the specimens were consolidated to a maximum vertical stress of 800 kPa and then unloaded to different stress levels. Loading (Figure 3(a)): in general, K0 was found to be constant during the loading stage. K0 of the glass beads was higher than K0 of the sand during loading (Figure 3(a)). For the glass beads, K0 decreased as the diameter decreased. As the consolidation loading increased, K0 increased slightly, as also observed by Mesri and Hayat [14] and Vardhanabhuti and Mesri [43]. However, the trend of the increase was not very large, which was similar to the trends reported by Wang and Gao [7]. The smallest sand (Ottawa sand 200–325) showed the highest value of K0. However, the Ottawa sand 20–30 (larger grains) exhibited the lowest K0. Unloading (Figure 3(b)): for all the test materials, K0 was found to increase as the vertical load was removed. In other words, K0 increased nonlinearly as the OCR increased. When the vertical stress decreased below 200 kPa (OCR greater than 3.0), K0 increased dramatically. Compared with the glass beads, the sand particles had a smaller K0 given the same vertical load. With the decreasing particle size, K0 of the glass beads decreased slightly. Unlike Ottawa sand 200–325, silica sand 30–50, Ottawa sand 50–70, and Ottawa sand 20–30 showed almost the same values of K0 during the entire loading range (Figure 3(b)).
(a)
(b)
3.2. Vs and Vs Anisotropy
Santamarina et al. [44] developed the following equation to uncover the stress-dependent characteristics of Vs:where is the mean normal stress, factor (m/s) is the velocity of a medium subject to 1 kPa confinement, and the exponent represents the amount of stress-dependent effect. For the vh and hv directions, . In the hh direction, . Figure 4 shows the shear wave velocities of the GB 3 mm, GB 0.5 mm, Ottawa sand 20–30, and Ottawa sand 50–70 samples during the loading and unloading stages in the vh, hv and hh directions. The curves with solid symbols represent the unloading behavior and the curves with open symbols represent the loading behavior. In general, three typical features were identified from figures. (1) The Vs values during unloading were slightly higher than the Vs values during the loading stage. (2) For the glass beads samples, Vs−hh was the smallest during loading. However, Vs−hh became the largest among the three orthogonal directions during the unloading stage (Figures 4(a) and 4(b)). (3) In contrast to the glass beads, the Ottawa sand displayed distinct behaviors during both loading and unloading (Figures 4(c) and 4(d)). Vs−hh was always the largest among the three directions in both the loading and unloading stages. The hierarchy of the Vs magnitude followed Vs−hh > Vs−hv > Vs−vh.
(a)
(b)
(c)
(d)
The Vs anisotropy of a soil is generally quantified by the ratio of Vs−hh/Vs−vh or Vs−hv/Vs−vh. The higher the ratio is, the higher the degree of anisotropy is. In order to compare the effect of the anisotropic stress induced Vs anisotropy and the inherent fabric induced Vs anisotropy, the ratios of Vs−hh/Vs−vh and Vs−hv/Vs−vh of the glass beads samples and the Ottawa sand samples were plotted, as shown in Figures 5(a) and 5(b). As observed in the figures, the stress-induced Vs anisotropy gradually increased when the vertical stress was gradually removed. The stress-induced Vs anisotropies of the glass beads and the Ottawa sand could reach up to 1.28 and 1.32, respectively. For the glass beads samples, the ratios of Vs−hh/Vs−vh and Vs−hv/Vs−vh were almost the same during unloading, which indicated that the anisotropy was largely due to the horizontal stress lock. For the Ottawa sand, however, the hh direction displaced the highest Vs anisotropy during unloading, indicating that fabric anisotropy governed the Vs anisotropy.
(a)
(b)
4. Discussion
4.1. Stress Dependency of Vs
The magnitude of Vs depends on the particle stiffness, effective stress conditions, and soil structure, as well as aging and other factors. The stress dependence of the shear wave velocity can be represented by equation (7), where the α term relates to the structural soil anisotropy. The value of the exponent β is affected by soil particle contact effects [44]. The fitted α and β were plotted as shown in Figure 6. Based on the observations, we found that most of the α and β values were plotted underneath the β = 1/4 area that represented the cone-to-plane contacts, and some of the data were plotted within the β = 1/6 area (Hertzian contact). As shown by the dotted line, a similar trend was observed for the fitted α and β values () and the hyperbolic relationship (equation (9)): ([44, 45])
The coefficient of determination was R2, and the values for the linear and hyperbola fittings were 0.75504 and 0.872071, respectively. In both the loading and unloading stages, equations (5) and (6) showed good fitting for the α values that were larger than 25 m/s. Equations (8) and (9) were no longer valid when the α value was smaller than 25 m/s because the governed contacts were not Hertzian elastic contacts.
4.2. Correlations between Vs and K0
The shear wave velocities from both the hh−vh and hh−hv modes and the K0 data during the unloading stages were fitted with equation (6). The results are presented in Figures 7 and 8. The 0.5 mm GB, silica sand 30–50, and Ottawa sand 50–70 were taken as examples to compare the fitted K0 and the direct measurements of the K0 from original test data. As indicated in Figure 8, the fitted K0 (lines with solid symbols) were very far off from the test data (lines with open symbols), which indicated that the fitting equation (6) might not have been applicable for the unloading behaviors. If the soil was highly overconsolidated, the estimated K0 might have been very far off from the real K0. The fitting parameters of “a” and “b” for all the test materials are presented in Figure 8. As indicated in Figure 8, the fitting parameter “a” was generally smaller than 1 and very constant (with a very small variation). However, the fitting parameter “b” showed very large scatter, which ranged from 0.5 to 6. The large variability of the fitting parameters and the “un-match” between the fitted K0 and the tested K0 indicated that the proposed fitting equation (6) might not be applicable for estimating K0 of granular soils during unloading conditions, especially for highly overly consolidated soils, because many other influencing factors were not accounted for in equation (6), e.g., the OCR.
Since K0 was found to be dependent on the OCR, the unsuccessful prediction of K0 with equation (6) might have been due to the omission of the OCR effects. The authors speculated that a good fitting between K0 and Vs might be achieved if the OCR was incorporated into the fitting equation. Therefore, the correlation (Equation (6)) that was developed by Ku and Mayne [5] was slightly modified so that the effect of the OCR was included, and the equation was rearranged:
By using the modified correlation (equation (10)), the entirety of the test data was replotted and compared with the tested K0 (see Figure 9). Similarly, the curves with solid symbols in the figure are the fitted K0, and the lines with open symbols are the tested K0. In contrast to the previous fitting, after incorporating the OCR into equation (6), all the fitted curves matched fairly well with the test data because all the estimated data points overlapped with the test data points (see Figure 9). The fitting parameters of “a” and “b” that were obtained from the modified fitting correlation (equation (9)) were plotted as shown in Figure 10. For the same vertical scale as that in Figure 8, the values of “a” and “b” from the modified correlation were very constant compared with the previous correlation (equation (6)). The average values of “a” and “b” for the modified correlation were 0.7 and 0.2, respectively. The good match between the fitted K0 and the tested K0 from the modified correlation and the smooth trends of the fitting parameters of “a” and “b” (see Figure 10) largely supported the concept that incorporating of the OCR could effectively reduce the variation and improve the correlation between Vs and K0 of highly overconsolidated granular soils.
4.3. Particle Shape and Size Effects on the K0 and Vs Anisotropy
In order to have a view of the microscopic surface texture of the test materials, scanning electron microscope (SEM) tests were conducted. The results are shown in Figure 11. The GBs were identical for the SEM, and all the GB particles were round in shape with smooth surfaces (Figure 11(a)). The relatively rough surfaces could be seen in the SEM photos of all the sands. The Ottawa sand particles were almost rounded to subrounded (Figures 11(b) and 11(d)). The silica sand 30–50 and Ottawa sand 200–325 particles were mostly subangular to angular (Figures 11(c) and 11(e)).
(a)
(b)
(c)
(d)
(e)
The glass beads samples tested in this study had almost the same initial density and packing. K0 was decreased slightly during unloading as the particle size decreased (Figure 3). The decrease of K0 implied that the degree of horizontal stress lock was decreased with the decreasing mean particle size. The force chain development among the granular materials has been evaluated in the literature [46, 47]. The smaller the particle was, the more uniform the stress distribution was inside the sample, and there was less chance of a “stress arch” occurring. Therefore, K0 might decrease slightly. In addition to the particle size, the angularity and the particle shape were found to influence K0 of the sand. The Ottawa sand 200–325 and the silica sand 30–50 both had high degrees of angularity, of which K0 were slightly higher than those of the other two sand samples under the same condition. The Ottawa sand 200–325 showed the highest K0 among all the sand specimens, which could be attributed to the relatively loose initial density and small critical-state friction angle. Mesri and Vardhanabhuti [3] reported that the looser samples normally have higher values of K0, which implied that K0 varied and was dependent on the initial density condition of the soils. By comparing the density index in Table 1, Ottawa sand 200–325 was found to have a relatively loose condition compared with the other sand samples. Thus, it exhibited a higher K0.
Stress anisotropy could cause the Vs anisotropy of the geomaterials composed of nonangular particles, such as round sand/spherical fly ash/glass beads [36], because the Vs anisotropy of coarse grained soils is majorly dependent on the stress components in the direction of shear wave propagation and polarization but relatively independent of the out-of-plane stress component. Under high overburden pressure, platy-clay particles often possessed preferential alignment which is the originality of the fabric anisotropy [48]. Therefore, is always higher than and [24]. The anisotropic initial state of stress, interparticle contact orientation, stress history of the medium, microcracks inside the materials, inclusions, and laminated media could also contribute to the fabric anisotropy of the clays, thus leading to the Vs anisotropy, i.e., the K0 consolidation induces the preferential alignment of clay fabric in a direction that is normal to the applied loading [24, 49].
During the one-dimensional consolidation, the anisotropic stress state could induce Vs anisotropy; e.g., Vs−hh was smaller than Vs−vh and Vs−vh (see Figures 4(a) and 4(b)) during loading. This has been proven by many researchers for the anisotropic loading condition [24, 34–36, 50]. In general, two types of Vs anisotropy could be observed, namely, the anisotropic stress state-induced Vs anisotropy and the inherent fabric anisotropy-induced Vs anisotropy. Since the glass beads were perfectly round, as indicated by the SEM images (Figure 11), the anisotropic stress state-induced Vs anisotropy was as shown in Figures 4(a) and 4(b) during the loading stage. However, during the unloading stage, the Vs anisotropic behavior was dramatically changed. Vs−hh gradually became the largest among the three orthogonal directions. This was related to the increased K0 during unloading, for which relatively high horizontal stresses were locked inside the consolidometer. The locked high horizontal stress could cause principal stress rotation inside the specimen when the horizontal direction became the major principal stress, and the vertical direction became the minor principal stress. The rotated principal stress states thus caused Vs−hh to evolve gradually to the highest Vs among the three orthogonal directions during unloading (the curves with solid symbols in Figures 4(a) and 4(b)).
The Ottawa sands had geometric anisotropy (indicated by the SEM images, where the particles had a subrounded shape with a longitudinal direction longer than the vertical direction). Due to the vertical stress, the sand grains might rotate and align more in the horizontal direction during loading, which was called the “preferred orientation.” Thus, the fabric might be anisotropic when the Vs anisotropy was manifested due to the fabric anisotropy that was induced by the anisotropic stress state. The anisotropic fabric of the Ottawa sand-induced Vs anisotropy during the loading stage when Vs−hh in general was higher than Vs−vh and Vs−hv (see Figures 4(c) and 5(b)). During unloading, the Vs anisotropy was manifested by the “stress rotation” when the horizontal stress became the major principal stress. Vs−hh remained the largest among the three orthogonal directions during the entire unloading stage, as indicated by the curves with solid symbols in Figures 4(c) and 4(d). The ratios of Vs−hh/Vs−vh and Vs−hv/Vs−vh could reach up to 1.32 and 1.10 during unloading.
5. Conclusions
The K0, Vs, and Vs anisotropies and their correlations with granular soils during one-dimensional consolidation loading were evaluated with a custom-designed floating-wall oedometer BE testing system. K0 stayed almost constant during loading. However, it increased dramatically when the consolidation loading was removed. The higher the OCR was, the greater the K0 was. K0 decreased slightly as the particle size decreased. The decrease of K0 could be attributed to the decrease of the “horizontal stress lock” due to the mean particle size decrease. The angularity and the particle shape were also found to influence K0 of the sand. Specimens with smaller particle sizes, rough surfaces, and angular geometry tended to have a lower value of K0. In addition, K0 also varied and it depended on the initial density condition of the soils. The Ottawa 200–325 sand was found to have a relatively loose condition compared with the other sand samples, thus exhibiting a higher value of K0. The stress dependence of Vs was evaluated, and Vs was found to increase with the increase of the norm stress. Both the anisotropic stress state and the fabric/geometry anisotropy of the sand had induced Vs anisotropy. During unloading, the Vs anisotropy was manifested by the “stress rotation and stress lock.” Furthermore, the directional Vs was validated to be a good tool for estimating the field stress state and predicting K0. The incorporation of OCR into Vs and K0 correlation effectively improved the prediction accuracy, which offers a robust way to predict K0 of both normally consolidated and highly overly consolidated granular soils in geotechnical practice.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
The work was supported by the National Natural Science Foundation of China (41977231) and Fundamental Research Funds for the Central Universities (CHD 300102269503).