Published: February 15, 2026

Penetration mechanism of CPT probe in the dense sand assisted by ultrasonic vibration

Faquan Wu1
Zhehong Shen2
Fang Zhang3
Bolong Liu4
Yuanlong Li5
Jie Wu6
Zhong Qing Chen7
Yuqi Lang8
Zihao Zhu9
Enjie Xu10
Weihong Chen11
Zhongli Xu12
1, 2, 3, 4, 5, 7Zhejiang-Greece International Joint Laboratory for Digital Geotechniques and Intelligent Exploration, Shaoxing, 312000, China
1, 2, 3, 4, 5, 7Key Laboratory of Rock Mechanics and Geohazards of Zhejiang Province, School of Civil Engineering, Shaoxing, 312000, China
1, 2, 3, 4, 5, 7, 8, 9, 10State Key Laboratory of Intelligent Deep Mental Mining and Equipment, Shaoxing University, Shaoxing, 312000, China
11, 12Tongchuang Engineering Design Co. Ltd., Shaoxing, 312000, China
6Zhejiang Yanchuang Technology Co, Ltd., Shaoxing, 312000, China
Corresponding Author:
Fang Zhang
Article in Press
Views 2
Reads 1
Downloads 11

Abstract

The cone penetration test (CPT) is widely used because of its minimal soil disturbance, simplicity, low-cost and dual functionality for exploration and testing. However, conventional CPT technology encounters significant difficulty in dense sands, which limits its applicability. To address this challenges, a novel ultrasonic-assisted CPT probe was developed based on an ultrasonic-vibration drilling technique, a series of penetration tests were performed in sand with relative densities (Dr) ranging from 30 % to 90 % to quantify the effects of Dr and ultrasonic vibration on penetration resistance and particle breakage. Besides, the particle breakage mechanism of sand and interaction between probe and sand under ultrasonic vibration were revealed. The results showed that the penetration load increases with Dr. Ultrasonic vibration effectively reduced the penetration load, however, the reduction decreased from 60 % to 18 % as Dr increased from 30 % to 90 %. As Dr increases, particle breakage around the cone tip becomes more pronounced, with a higher breakage rate closer to the cone tip. Under the ultrasonic-static coupling, particle breakage was mainly caused by ultra-high-frequency impacts between particles and the cone tip, particles collided with each other, vertical and shear stress around the probe, The stress state of the sand and soil around the ultrasonic CPT probe was divided into five regions. This study introduces a novel ultrasonic vibration CPT probe, quantifies its penetration benefits in the dense sand, and reveals the particle breakage mechanisms under the ultrasonic-static coupling. The investigation will provide insights for applications of CPT in dense sands.

1. Introduction

The cone penetration test (CPT) is an efficient in-situ testing technology in geotechnical engineering and is widely used to assess foundation bearing capacity, classify soil strata, and evaluate liquefaction potential using real-time measurements of cone tip resistance (qc) and sleeve friction resistance (fs) [1]. Its standardized operating procedures and wealth of empirical correlations make CPT a core tool for site investigation under complex geological conditions [2]. However, conventional CPT faces significant limitations in dense sand or hard strata [3]. On the one hand, high penetration resistance causes probe wear and may even interrupt testing, compromising data continuity. On the other hand, the effect of particle breakage on the soil’s physical and mechanical properties is not yet well understood, resulting in systematic deviations in parameter interpretation based on classical elastoplastic theory.

In recent years, the vibratory pile driving technology have been developed rapidly [4-8]. Many studies have shown that vibratory pile driving can effectively reduce the penetration resistance of piles and increase the driving rate. Researchers have introduced vibration factors into soil-structure interaction model experiments to explore the impact of vibration on the mechanical properties of granular materials such as sand. For example, Denies et al. [9] conducted sphere penetration tests in vibrated sand. It is found that vibration significantly increases the sinking rate of the penetrator in the sand samples. Baptiste et al. [10] performed small-scale vibratory pile penetration tests with glass beads, it was found that pile base vibration caused localized fluidization of the surrounding glass material, reducing penetration resistance by a factor of about ten. Denies [11], [12] pointed out that when the vibration amplitude exceeds 1 g, vibration causes fluidization in the sand and leads to a complete loss of shear strength as the acceleration increases. The experimental results of Johnson [13], [14] showed that under vibration, local fluidization occurs in the shear band of granular materials, resulting in strength degration. Jiang [15] developed a soil-structure interface shear testing device capable of coupling high-frequency vibration with static shear, revealing the strength degradation behavior of granular material-structure interfaces under vibration. Ye [16] improved a triaxial apparatus and applied axial high-frequency vibration to sand specimens, observing a significant reduction in deviatoric stress accompanied by additional axial strain. In summary, the interaction between granular materials and structures under vibration conditions will change. However, in CPT, vibration is rarely taken into consideration to improve the feasibility in dense sands.

Until now, the ultrasonic rock-breaking technology in geotechnical engineering has demonstrated its feasibility for assisting penetration [17]. Researchers have developed a range of devices, including direct-drive ultrasonic drills, rotary ultrasonic-vibration drills, and high-frequency acoustic rotary-impact rock-breaking drills [18-20]. These devices operate on the principle that an ultrasonic transducer converts electrical energy into high-frequency mechanical vibrations, enabling efficient drilling. For example, Li [21] designed a rotary-impact ultrasonic sampler for lunar soil and permafrost resource development and conducted drilling tests in simulated lunar permafrost. The results showed that the average drilling force was reduced by 52 % under the application of ultrasonic vibration. These achievements provide a clear pathway for integrating ultrasonic technology into static penetration testing, it is expected to break through the technical limitations of CPT equipment in dense strata. Given the difference in loading mode, the mechanisms of resistance reduction in ultrasonic vibration CPT remain unclear.

In terms of the research on the penetration mechanism in sandy soil layer, the influence of particle breakage on the mechanical behavior of soil during the piling has been widely investigated. Many researchers have systematically investigated the distribution characteristics of the particle breakage zone at the pile tip through pile penetration tests. Based on visual observation technology and particle image velocimetry (PIV) detection, it was found that there is a conical “nose cone” particle breakage zone below the pile tip, the breakage in the zone is uneven, and the breakage mainly concentrated in the shear zone. Especially, the maximum breakage degree directly below the pile tip [22], [23]. The breakage zones in the sand can be classified into a severely fragmented shear zone, significantly breakage non-shear zone, and a moderately fragmented non-shear zone [24]. According to the Digital image correlation (DIC) and X-ray tomography, the fine particle shear bands and radial flow within the soil-pile interface was observed [25], [26]. The relationship between particle fragmentation and penetration resistance is a key research focus in pile foundations and static penetration testing (CPT). Areshed [27] conducted static penetration tests on three different types of crushable silica sands, whereas Peng [28] performed pile penetration tests on coral reef sand foundations. Both studies found that the penetration load was lower was lower for crushable sands than for non-crushable sands. In addition, Mao et al. [29-32] demonstrated experimentally that acoustic emission (AE) can effectively monitor the bearing behavior of piles. The time histories of multiple AE parameters were found to closely resemble the load-settlement curves of model piles, and the location and evolution characteristics of fragmentation during penetration could also be pinpointed. In recent years, the Discrete Element Method (DEM) has been widely used to simulate CPT of granular materials, enabling the reproduction of soil discontinuities and large deformations [33]. O. Falagush and colleagues [34] used DEM to study the effect of particle fragmentation on the cone resistance in CPT and found that the cone resistance in the fragmentable model was significantly lower than in non-fragmented model. Abouzar [35], Ge [36], Huang [37], and Song [38] conducted numerical simulations of static penetration tests in a calibration chamber and analyzed the factors influencing cone resistance from both macro-and micro-scales. These findings provide valuable insights for further understanding the penetration mechanism in CPT. However, most existing studies have focused on pile penetration. The combined effects of static loading and vibration- especially high-frequency vibration, have not been explored in the context of ultrasonic and static penetration (CPT). However, existing studies mainly focus on the static penetration, the particle-breakage mechanism under coupled static loading and high-frequency vibration remains unclear.

To enhance CPT penetration capacity, a novel ultrasonic CPT probe was developed in the study. Subsequently, ultrasonic CPT were conducted on the sand with different Dr, the effectiveness of the device is verified. Subsequently, the effects of Dr and ultrasonic vibration on the penetration load and the particle breakage were quantified. The particle-breakage mechanism under ultrasonic-static coupling and the probe-sand interaction under ultrasonic vibration were revealed. The research process is shown in Fig. 1. The main contributions of this study are as follows: (1) A new ultrasonic CPT probe was developed, achieving dual functions of dynamic and static cone penetration. (2) The influence of relative density of the sand and ultrasonic vibration on the penetration load of the ultrasonic probe was analyzed. (3) The particle breakage mechanism and stress state of the sand under the ultrasonic-static coupling were revealed.

Fig. 1Research process. The photos were taken by Yuanlong Li at Shaoxing University on March 10, 2024

Research process. The photos were taken by Yuanlong Li  at Shaoxing University on March 10, 2024

2. Experimental methodology

This section describes the operating principle and procedures of the ultrasonic CPT probe. To verify the device’s feasibility, we conducted, a series of ultrasonic CPTs in sand at different relative densities.

2.1. The design of the ultrasonic CPT probe

The structure of the ultrasonic CPT probe is shown in Fig. 2. The external dimensions of this probe are identical to those of the conventional dual-bridge CPT probe, with a total length of 475 mm and a diameter of 43.7 mm. The cross-sectional area and the cone angle of the tip are 15 cm2 and 60°, The outer diameter of the thread is 39 mm. The ultrasonic vibration unit is integrated within the CPT probe. The ultrasonic CPT probe consists of the cone tip, extension sleeve, ultrasonic transducer, ultrasonic generator, friction sleeve, metal rod and a sensor for measuring the total penetration resistance. The extension sleeve is a 110-mm-long stainless steel tube that connects the cone tip to the main body of the probe, with the ultrasonic transducer installed inside the extension sleeve. The ultrasonic transducer consists of front and rear cover plates, four PZT-8-type piezoelectric ceramics with an outer diameter of 25 mm, electrode slices, wires, and bolts. When powered, the ultrasonic generator rectifies and filters the 220 V (50/60 Hz) AC supply to 310 V DC, which is then converted into a specific high-frequency AC drive signal. The high-frequency signal is amplified to several kilovolts by the high-frequency transformer to drive the ultrasonic transducer. Under high-voltage, high-frequency excitation, the piezoelectric ceramic discs deform via the inverse piezoelectric effect, electrical energy is converted into mechanical vibrational energy. The transducer transmits the ultra-high frequency axial vibrations to the cone tip, thereby intensifying soil agitation around the tip during penetration, reducing the penetration resistance. The ultrasonic generator is equipped with a control module that tunes the operating frequency and power for penetration. To prevent the ultrasonic wave transmit to the probe body, a nylon threaded was installed at the outer layer of the cone tip, which concentrates the energy at the cone tip. Thereby minimizing measurement interference and improving the accuracy of the force sensor [40, 41].

Fig. 2Structure of the ultrasonic CPT probe. The photos were taken by Zhehong Shen at Shaoxing University on October 23, 2024

Structure of the ultrasonic CPT probe. The photos were taken  by Zhehong Shen at Shaoxing University on October 23, 2024

To measure the total penetration resistance, a sensor is mounted at the top of the probe. The sensor has a 50 kN capacity with an accuracy of ±0.5 %. The sensor operates based on resistive strain-gauge force measurement. When an axial load is applied to the load-bearing member in the sensor, the electric resistance is changed. According to the relationship between axial strain (Eq. 1), the strain is obtained. Subsequently, according to Eqs. (2-3), the penetration resistance is measured:

1
ΔRR=GF·ε,
2
σ=Eε,
3
F=Aσ,

where R is the nominal resistance, R is the change in electric resistance, GF represents the gauge factor. A is the effective area of load-bearing member, E is Young’s modulus of the load-bearing member. ε and σ are strain and stress, respectively. F is the total penetration resistance. To convert the small resistance variation into a measurable electrical signal, the gauges are wired in a Wheatstone bridge. With an excitation voltage Vex, the bridge output voltage Vo can be approximated for small strains as:

4
V0=VexKbGFε,

where Kb is a configuration-dependent coefficient (for a full-bridge configuration, Kb 1). Substituting Eq. (3) yields the force-voltage relation:

5
F=EAVexKbGFVO=CF,

where, CF is the conversion coefficient provided by the manufacturer through static-load calibration, which is used to convert the output voltage of the acquisition system (after signal conditioning/amplification and data acquisition) into the total penetration load F.

Due to the presence of the cone tip, the vertical penetration load is resolved on the cone surface into a normal component and a frictional component; the cone surface is subjected to the normal stress σn and friction characterized by the cone-soil friction coefficient μ, as shown in Fig. 3. The probe shaft is subjected to the horizontal normal stress σh and the corresponding shaft friction fs. Under simplified assumptions (Coulomb friction and an equivalent mean stress acting on the contact surface), the force equilibrium together with the probe geometry yields the following relationship between the cone resistance qc and the normal stress on the cone surface:

6
qc=σn1+3μ.

Considering an effective contact area at the cone tip (taken as the circular area with probe diameter D), the total cone-tip resistance is:

7
Qc=σn(1+3μ)4πD2.

The unit shaft friction can be expressed in a Coulomb form as:

8
fs=μσh,
9
QS=0LfszπDdz.

Fig. 3Forces acting on the probe

Forces acting on the probe

Accordingly, the total penetration load F acting on the probe can be expressed as the sum of the cone-tip and shaft contributions:

10
F=Qc+Qs.

It should be noted that the force sensor in this study directly measures the total penetration load F (in N), which is used to evaluate the resistance-reduction effect of ultrasonic vibration.

The operating procedure of ultrasonic CPT is similar to that of the conventional CPT. During the test, the ultrasonic CPT probe is connected to the ultrasonic generator, and the CPT probe is steadily advanced into the soil with a rate of (20±5) mm/s. Total penetration load are monitored by sensors. When hard soil layers are encountered, the ultrasonic generator is activated to induce high-frequency vibrations at the cone tip. Depending on the geological conditions, parameters such as ultrasonic frequency, amplitude, and power can be adjusted accordingly to optimize penetration performance.

2.2. Ultrasonic CPT in dense sandy soils

2.2.1. Preparation of the sand specimen

The sand sample was collected from a river in Shaoxing, Zhejiang Province. Based on sieve test, the grain size distribution curve was obtained, as shown in Fig. 4. The median particle diameter was D50= 0.447 mm, the coefficients of uniformity was Cu= 4.67, and the coefficient of curvature was Cc= 1.26. The maximum and minimum dry densities were measured using the vibration-hammer method and funnel-cylinder method, respectively. The maximum dry density was 1.96 g/cm3, and the minimum was 1.45 g/cm3. In the experiment, the sand was compacted in layers within a model box to simulate a sand foundation. To ensure uniform compaction, the sand was compacted in 10 layers, each 4 cm thick. Five different relative densities (Dr) were prepared for sand: 30 %, 45 %, 60 %, 75 %, and 90 %. The dry densities corresponding to the five relative densities were 1.57, 1.64, 1.72, 1.8, and 1.89 g/cm3, respectively. According to the density-volume relationship, the required masses for each target relative density were 78.92, 82.44, 86.46, 90., and 95 kg, respectively.

2.2.2. Experimental equipment

In this study, a series of ultrasonic CPTs were conducted to examine the effects of sand relative density and ultrasonic vibration on the penetration load of the ultrasonic probe and particle-breakage mechanism in sand. The experimental setup comprises an ultrasonic CPT probe, a model box, a loading system and load-monitoring system, as illustrated in Fig. 5. The model box is 40 cm in diameter and 50 cm in height. The loading system includes a reaction frame and a hydraulic jack with a maximum stroke of 600 mm. The hydraulic jack applies normal force to the ultrasonic CPT probe. The load-monitoring system consists of pressure sensors, a computer, and data-acquisition software.

Fig. 4The particle distribution curve of the river sand

The particle distribution curve of the river sand

Fig. 5Schematic diagram of the experimental setup

Schematic diagram of the experimental setup

2.2.3. Experimental process

The ultrasonic CPT procedure for dense sandy was as follows:

(1) The sand was compacted in 10 layers to the target relative density until the total height reached 400 mm. After that, the sample was then allowed to stand for 12 hours to achieve a uniform stress distribution in the sand (Fig. 6(a)).

(2) The ultrasonic CPT probe, ultrasonic generator, and computer were connected. The sensor sampling frequency was set to 10 Hz, and the ultrasonic frequency was set to 30.71 kHz before loading.

(3) After the preparation, the ultrasonic CPT probe was advanced into the sand by the hydraulic jack at a rate of 2 mm/s. At the same time, the ultrasonic generator and load-acquisition software were activated simultaneously. The sum of the cone resistance and sleeve friction was recorded. When the penetration depth reaches 160 mm, the test terminated, and the ultrasonic generator and load-monitoring system were switched off (Fig. 6(b)).

(4) After penetration, sand samples surrounding the ultrasonic CPT probe were collected (Fig. 6(c)). As shown in Fig. 7, regions around the cone tip were defined at distances of 85-125 mm (A), 45‐85 mm (B), 0‐45 mm (C), –20‐0 mm (D), and –20 to –40 mm (E) from the tip. The sand in A, B, C, D, and E were collected using ring cutters with a diameter of 61.8 mm. To facilitate sampling, the sand was moistened before the penetration of the ring cutters. The collected samples from each region were then placed in an oven at 70 °C and dried for 24 hours until completely dry.

(5) Sieving tests were carried out on the sand samples. The microscopic of adhering to the probe sidewall in regions A, B, and C was observed under a stereomicroscope (Figs. 6(d-e)).

Fig. 6Schematic of ultrasonic CPT in sand

Schematic of ultrasonic CPT in sand

a) Layered compaction

Schematic of ultrasonic CPT in sand

b) Penetration

Schematic of ultrasonic CPT in sand

c) Sampling the sand around the ultrasonic CPT probe

Schematic of ultrasonic CPT in sand

d) Sieving tests

Schematic of ultrasonic CPT in sand

e) observation of the microstructure of sand grains

Fig. 7Schematic diagram of excavation regions

Schematic diagram of excavation regions

3. Effect of relative densities of the sand on the penetration load and particle breakage under the penetration of the ultrasonic CPT probe

This section examines the penetration load behavior from the perspectives of Dr and ultrasonic vibration. First, the variation of penetration load-depth curves and penetration mechanisms under ultrasonic penetration at different Dr is analyzed. Subsequently, conventional penetration and ultrasonic-assisted penetration are compared at the same relative density to quantify the resistance-reduction effect induced by ultrasonic vibration and its dependence on Dr.

3.1. Penetration load under different relative densities

The variations of the penetration load with penetration depth for different relative densities are shown in Fig. 8(a), All curves exhibit a similar pattern trend across different relative densities, and can be divided into a linear stage and nonlinear stage. As Dr increases, the penetration load increases gradually, with a greater rate of increase at higher Dr. Taking a depth of 140 mm as an example, Fig. 8(b) plots penetration load versus Dr. When Dr= 30 %, 45 %, 60 %, 75 %, and 90 %, the corresponding penetration loads were 215.6 N, 519.4 N, 1372 N, 2508.88 N, and 3635.8 N, respectively.

The results indicate that the penetration load increases with the increasing Dr. This is because both the internal friction angle and the shear strength of sand increase with relative density [42, 43], and higher Dr strengthens the soil’s force-chain network, thereby increasing the penetration load. At the initial stage of loading, the penetration load increases linearly with depth, and the sand particles around the probe are compressed [31, 44]. Sandy with low relative density has higher compressibility, providing more space for particle rearrangement and reorganization during penetration, resulting in a longer linear stage. In contrast, sandy with high relative density has lower compressibility, and more tightly packed particles, making particle motion difficult and shortening the duration of the linear stage. As the penetration depth increases, the penetration load rises nonlinearly, with a gradually increasing slope. During this stage, the shear load applied by the probe causes plastic deformation and particle fragmentation in the surrounding sand [30]. According to Coulomb’s law, the shear strength of sand is proportional to the normal stress on the shear plane. Sliding friction between particles and the interlocking of rough surfaces generate frictional resistance, which increases with Dr. As Dr increases, soil porosity decreases, the contact area between the soil and the probe increases, and the particles become more tightly packed. This increases the normal stress between particles and the internal friction angle, thereby enhancing the shear strength of the soil. Therefore, a greater penetration load is required to penetrate sand at high Dr. Additionally, fluctuations are observed in the penetration curves for Dr= 60 %, 75 %, and 90 %, which can be attributed to the layered artificial compaction of the sand. At higher Dr, compaction is more challenging, make it difficult to ensure uniformity.

Fig. 8Penetration load acting on the CPT probe at different relative densities

Penetration load acting on the CPT probe at different relative densities

a) Penetration load- depth curve

Penetration load acting on the CPT probe at different relative densities

b) Relative density-penetration load curve

3.2. Penetration load under static load and ultrasonic vibration

The variations of the total penetration load with penetration depth for both the CPT and ultrasonic CPT probes in the sand at different relative densities are shown in Fig. 9. From the figure, it can be observed that the trend of the total penetration load under ultrasonic vibration is similar to that of the static CPT, the total penetration load increases with increasing the penetration depth. Furthermore, as shown in Fig. 10(a), at the same penetration depth, a higher relative density corresponds to a greater total penetration load. Additionally, under the ultrasonic vibration, the total penetration load is significantly reduced. However, as shown in Fig. 10(b), the reduction in penetration load decreases as the relative density of the sand increases. At a penetration depth of 140 mm, the total penetration load for ultrasonic CPT is reduced by 60.0 % for Dr= 30 %, 41.0 % for Dr= 45 %, 47.0 % for Dr= 60 %, 24.0 % for Dr= 75%, and 18.0 % for Dr= 90 %.

Fig. 9Curve of penetration load with depth for ultrasonic vibration probe and conventional probe

Curve of penetration load with depth for ultrasonic vibration probe and conventional probe

a)Dr= 30 %

Curve of penetration load with depth for ultrasonic vibration probe and conventional probe

b)Dr= 45 %

Curve of penetration load with depth for ultrasonic vibration probe and conventional probe

c)Dr= 60 %

Curve of penetration load with depth for ultrasonic vibration probe and conventional probe

d)Dr= 75 %

Curve of penetration load with depth for ultrasonic vibration probe and conventional probe

e)Dr= 90 %

The experimental results indicate that ultrasonic vibration effectively reduces penetration load, thereby enhancing the penetration capability of the CPT. However, compared to loose sand, the reduction in penetration resistance is smaller in high relative density sands. In both penetration modes, the probe is subjected to normal and shear stress, causing compression and shear deformation of the surrounding sand. Gong [39] utilized the DEM to clarify that under vibration, the contact normal anisotropy, normal contact force anisotropy, and tangential contact force anisotropy parameters of granular soil all decrease, leading to a reduction in the shear strength of the sand. During the vibratory penetration process, the high-frequency vibration of the cone tip causes particles in the contact area to rearrange, reducing the normal contact force between particles and the frictional force at the probe interface. The sand particles undergo fluidization [40-42], and the force chains between sand particles break [11], resulting in a local decrease in stiffness [43]. Therefore, the penetration load under vibratory conditions was observed to be lower than that under conventional penetration conditions.

Fig. 10Comparison of the penetration load of acting on the CPT probe and ultrasonic CPT probe at penetration depth of 140 mm

Comparison of the penetration load of acting on the CPT probe  and ultrasonic CPT probe at penetration depth of 140 mm

a) Penetration load-relative density curve

Comparison of the penetration load of acting on the CPT probe  and ultrasonic CPT probe at penetration depth of 140 mm

b) Reduction in the penetration load-relative density curve

Research has found that the degree of shear strength attenuation is significantly affected by the sample density. As the sample density increases, the weakening effect of vibration on the strength of sandy soil diminishes [44]. Denies [12] suggested that vibration causes particle collisions, generating a “shaking pressure” effect, which macroscopically resembles excess pore-water pressure in saturated soils. This effect reduces the internal friction angle of the sand and weakens the effective normal stress. In loose sand, the friction between the soil and the probe is lower, and the larger pores allow particles around the cone tip to move and rotate. During the penetration process, soil particles flow around the probe and collide, resulting in fragmentation. This combined effect significantly reduces the penetration load. In the penetration of high relative density sand, the friction between the soil and the probe is greater, and the influence of the cone tip vibration is limited to a smaller area. The smaller pore space and tighter particle contact restrict particle movement and rotation, increasing the strength of the force chains between particles [45-47]. This suppresses the strength attenuation caused by vibration, leading to limited particle fragmentation. As a result, the reduction in penetration load is less compared to that in low relative density sand.

4. Effect of ultrasonic vibration and relative densities on the broken mechanism of sand particles

Under ultrasonic vibration-static coupling, sand particles are broken, thereby reducing the penetration load. To reveal the particle breakage mechanism at different relative densities and ultrasonic vibration, we first investigated the microscopic structure of the sands grains adhering to the surface of the ultrasonic CPT cone tip, owing to the high stresses and intense vibration in this region. Subsequently, the sand particle breakage along the ultrasonic CPT probe was studied by the breakage index Bg. Finally, a quantitative comparison of particle breakage between ultrasonic-assisted CPT and conventional CPT is conducted.

Fig. 11Microstructural characteristics and average particle size of the sands grains adhering to the cone tip at different relative densities. The photos were taken by Zhehong Shen at Shaoxing University on November 8, 2024

Microstructural characteristics and average particle size of the sands grains adhering  to the cone tip at different relative densities. The photos were taken  by Zhehong Shen at Shaoxing University on November 8, 2024

a)Dr= 30 %

Microstructural characteristics and average particle size of the sands grains adhering  to the cone tip at different relative densities. The photos were taken  by Zhehong Shen at Shaoxing University on November 8, 2024

b)Dr= 45 %

Microstructural characteristics and average particle size of the sands grains adhering  to the cone tip at different relative densities. The photos were taken  by Zhehong Shen at Shaoxing University on November 8, 2024

c)Dr= 60 %

Microstructural characteristics and average particle size of the sands grains adhering  to the cone tip at different relative densities. The photos were taken  by Zhehong Shen at Shaoxing University on November 8, 2024

d)Dr= 75 %

Microstructural characteristics and average particle size of the sands grains adhering  to the cone tip at different relative densities. The photos were taken  by Zhehong Shen at Shaoxing University on November 8, 2024

e)Dr= 90 %

Microstructural characteristics and average particle size of the sands grains adhering  to the cone tip at different relative densities. The photos were taken  by Zhehong Shen at Shaoxing University on November 8, 2024

f) Average particle size (mean ± standard deviation) of sand grains adhered to the cone tip

4.1. Particle breakage of the sand under different relative densities

4.1.1. Particle breakage of the sand adhered to the ultrasonic CPT probe tip

Figs. 11(a-e) show the microscopic structure of the sand grains adhered to the surface of the ultrasonic CPT cone tip at different relative densities. Based on the experimental results, approximately 100 particles were measured for each relative density using ImageJ, from which the mean particle size and the standard deviation were obtained (Fig. 11(f)). At Dr= 30 %, both broken and unbroken grains are present at the tip. Beneath the larger unbroken grains, numerous small broken grains (0.1-0.25 mm) are closely attached to the cone surface. The broken grains exhibit sharp edges (Fig. 11(a)), with an average particle size of 331.8±41.6 μm. At Dr= 45 %, the overall characteristics of the sand particles are similar to those observed at Dr= 30 %, with broken and unbroken particles coexisting, but the average particle size decreases to 200.3±20.3 μm (Figs. 11(b, f)). For Dr= 60 %-90 %, particle size of the sand decreases significantly, and the grains appear more extensively crushed (Figs. 11(c-e)), with mean sizes of 71.8±11.1 μm (Dr= 60 %), 64.6±14.8 μm (Dr= 75 %), and 39.7±4.4 μm (Dr= 90 %). Apart from a few larger grains, most particles have been crushed into fine powder, indicating a more uniform degree of breakage.

The results indicate that dense sand facilitates particle breakage around the cone tip. As sand relative density increases, the contact area between the sand and the CPT probe increases, Ultrasonic vibration effectively crushes the sand surrounding the probe, thereby enhancing the penetration capability of the CPT probe in sand. XRD analysis revealed that the sand consists primarily of quartz, feldspar, and mica. Contacts between the mineral particles tight packing (frictional contacts) and cemented bonds. In dense sand, particles are packed more tightly than in loose sand. Under ultrasonic vibration loading, stress is transmitted through the cementing materials and particles, leading to local stress concentrations. When the stress exceeds the tensile or shear strength of the cementing bonds, particles are broken. Additionally, repeated ultrasonic loading induces fatigue damage in the cementing bonds, ultimately increasing particle breakage (Fig. 12).

Fig. 12Breakage of sand grains under ultrasonic vibration

Breakage of sand grains under ultrasonic vibration

a) Loose sand

Breakage of sand grains under ultrasonic vibration

b) Dense sand

4.1.2. Sand particle breakage along the ultrasonic CPT probe

Under the combined effects of the penetration of CPT probe and ultrasonic wave, the sands surrounding the probe are broken. Currently, the breakage index Bg is commonly used to define the degree of particle breakage. This index was first introduced by Marsal [48]:

11
Bg=i=1n|ai-bi|,

where ai and bi are the mass percentage of particles in the size range i before and after the breakage, respectively, n is the total number of size ranges for a particle sample.  stands for adding up the mass percentage from the largest particle to the smallest one.

After sieving the sand samples from different regions, the breakage index (Bg) was computed for regions along the ultrasonic CPT probe. Fig. 13(a) illustrates the variation of Bg with the mean depth of regions A-E at different realtive densities. The results indicate that Bg generally increases with depth, with the highest value occurring in region D, adjacent to the cone tip. Additionally, the average breakage index for regions A-E was calculated at different relative densities. As shown in Fig. 13(b), the average breakage index generally decreases as relative density increases.

Fig. 13Particle breakage rate of sandy soil

Particle breakage rate of sandy soil

a) Breakage indexes of sandy soil at different depths

Particle breakage rate of sandy soil

b) Average breakage indexes of sandy soil adjacent to the cone tip under different relative densities

Fig. 14The microstructure of sand grains along the ultrasonic CPT probe (Dr= 75 %). There photos were taken by Zhehong Shen at Shaoxing University on November 8, 2024

The microstructure of sand grains along the ultrasonic CPT probe (Dr= 75 %).  There photos were taken by Zhehong Shen at Shaoxing University on November 8, 2024

a) Region A

The microstructure of sand grains along the ultrasonic CPT probe (Dr= 75 %).  There photos were taken by Zhehong Shen at Shaoxing University on November 8, 2024

b) Region B

The microstructure of sand grains along the ultrasonic CPT probe (Dr= 75 %).  There photos were taken by Zhehong Shen at Shaoxing University on November 8, 2024

c) Region C

The microstructure of sand grains along the ultrasonic CPT probe (Dr= 75 %).  There photos were taken by Zhehong Shen at Shaoxing University on November 8, 2024

d) Region D

Since particle fragmentation primarily occurs near the CPT probe, we further analyzed its characteristics by examining the microstructure of sand grains adhering to regions A-C and the cone tip with a stereomicroscope after the penetration test. As shown in Fig. 14, using sand at Dr= 75 % as an example, in region A, unbroken sand grains with diameters of 0.2-0.5 mm were adhered to the probe surface (Fig. 14(a)). In region B, the adhering particles are generally smaller than those in region A, and broken sand grains with diameters of 0.074-0.25 mm were observed (Fig. 14(b)). In region C, numerous broken sand grains with diameters of 0.074-0.25 mm were observed, indicating significant breakage (Fig. 14(c)). At the cone tip, abundant sand powder was observed, exhibiting the smallest average particle size (Fig. 14(d)). Therefore, breakage is concentrated at the cone tip owing to the intense ultrasonic vibration.

4.2. Particle breakage characteristics under ultrasonic-static coupling

Similar to section 4.1.2, the breakage indexes Bg for zones A, B, C, D, and E under CPT and ultrasonic CPT were calculated by Eq. (11), with the results shown in Figs. 15(a-e). Additionally, the average breakage indexes for zones A, B, C, D, and E for sand with different relative densities are presented in Fig. 14(f). Overall, under ultrasonic vibration, Bg in zones A, B, C, D, and E gradually increases, which is consistent with the microstructures of the sand adhering to the CPT probe. However, the Bg for ultrasonic CPT and CPT differ. Specifically, when Dr= 30 % and 45 %, the trends of Bg under ultrasonic CPT are similar to those under CPT, Bg increases with depth and reaching its peak in region D. Additionally, Bg under ultrasonic CPT is generally higher than that under CPT (Figs. 15(a-b)). However, when Dr= 45 %, the Bg in zone D (beneath the cone tip) is slightly lower under ultrasonic CPT than under traditional CPT. When Dr= 60 %, Bg in zones D and E (beneath the cone tip) under ultrasonic CPT is slightly lower than that under CPT (Fig. 15(c)). When Dr= 75 % and 90 %, Bg under ultrasonic CPT is significantly lower than that under CPT (Figs. 15(d-e)). As shown in Fig. 15 (f), when Dr= 30 % and 45 %, the average Bg in regions A, B, C, D, and E under ultrasonic CPT are higher than those under CPT. In contrast, when Dr= 60 %, 75 % and 90 %, the average Bg under ultrasonic CPT are lower than those under CPT.

5. Discussions

Base on ultrasonic CPT conducted in sand with different relative densities, we investigated the effects of relative density and ultrasonic vibration on the penetration load and particle breakage. Essentially, under the ultrasonic-static coupling, particles are prone to break, the interaction between the ultrasonic CPT probe and the sand is altered. Thus, cone-tip resistance and the friction between the probe and the sand decrease, thereby reducing load acting on the probe. In this section, we elucidate the particle breakage mechanism under the ultrasonic-static coupling. Subsequently, based on the experimental results, we investigated the probe-sand interaction under ultrasonic vibration.

5.1. Particle breakage mechanism under the ultrasonic-static coupling

According to the experimental results, particle breakage mechanism under the ultrasonic-static coupling are illustrated in Fig. 16.

(1) Ultrahigh-frequency impacts between particles and the cone tip. The large voids between particles allow high-frequency vibrations of the sand under the ultrasonic excitation, causing particle movement and rotation (Fig. 16(a)). During motion, sand particles collide with the cone tip at high frequency, promoting particle breakage. The maximum acceleration of the vibrating cone tip is given by: a=A2πf2. Substituting A= 1 μm and f= 30.71 kHz, the calculated maximum acceleration a= 3.7×104 m/s2, which is far greater than gravitational acceleration g. Under the influence of high acceleration, sand particles around the cone tip experience ultrahigh-frequency impacts, leading to localized breakage. The breakage increases sand compressibility, resulting in a reduction in penetration load.

Fig. 15Comparison of breakage index between CPT and ultrasonic CPT

Comparison of breakage index between CPT and ultrasonic CPT

a)Dr= 30 %

Comparison of breakage index between CPT and ultrasonic CPT

a)Dr= 45 %

Comparison of breakage index between CPT and ultrasonic CPT

c)Dr= 60 %

Comparison of breakage index between CPT and ultrasonic CPT

d)Dr= 75 %

Comparison of breakage index between CPT and ultrasonic CPT

e)Dr= 90 %

Comparison of breakage index between CPT and ultrasonic CPT

f) Average breakage index

Particle breakage caused by ultrasonic vibration can be attributed to: (a) Stress-wave propagation: High-frequency stress waves generated by ultrasonic vibration propagate through the sand, creating localized high-pressure and low-stress zones between particles. The rapid variation in local stress weakens the frictional and cohesive forces between sand particles, thereby further promoting particle breakage. (b) Resonance effect: When the ultrasonic vibration frequency aligns with the sand’s resonance frequency of the sand, the internal vibration amplitude increases significantly, further reducing frictional forces between particles, making them more susceptible to breakage. (c) Fatigue-damage accumulation: Prolonged ultrasonic vibration induces fatigue damage in sandy soil. Under repeated vibrational loading, microcracks and pores within the sand particles gradually expand, eventually leading to structural failure of the soil. This cumulative fatigue damage process makes the sand particles more prone to breakage.

(2) Reduction in local sandy density near the cone tip due to ultrahigh-frequency vibration. Assuming the penetration velocity v0 of the CPT probe. Under the vibration, the velocity of the probe varies between (v0-2πfA, v0+2πfA). Given the ultrasonic vibration frequency f= 30.71 kHz and the amplitude A= 1 μm, the maximum vibration velocity of the cone tip can be calculated, v=2πfA= 0.19 m/s. When the penetration velocity v02πfA), the relative velocity range between the cone tip and the surrounding sand falls within (v0-2πfA, 0), indicating periodic separation between the cone tip and the sand. This separation, along with particle movement and rotation, results in a decrease in local sand density. Since the experimental penetration velocity satisfies the separation condition (v0 0.19 m/s), the penetration load is reduced.

(3) Vertical stress and shear stresses from the CPT probe cause sand-particle breakage. In sand with medium to high relative density, the voids between particles are smaller, and particles are difficult to move and rotate under ultrasonic vibration. Consequently, high-frequency collisions between sand particles and the cone tip are suppressed. In high-density sand, the propagation range of ultrasonic waves is reduced, meaning that the outer soil layers are less affected by ultrasonic vibrations, leading to a lower breakage rate. At this stage, the primary effect of ultrasonic vibration is to loosen the compacted sand surrounding the probe (Fig. 16(b)).

Fig. 16Schematic diagram of the interaction mechanism between ultrasonic CPT probe

Schematic diagram of the interaction mechanism between ultrasonic CPT probe

a) Loose sand

Schematic diagram of the interaction mechanism between ultrasonic CPT probe

b) Dense sand

5.2. Interaction between probe and sand under ultrasonic vibration

During penetration, the stress state of the sand surrounding the probe varies. According to observations by Yang et al. [24], [54], on particles in the cone-tip region, the area can be divided into three zones: Zone I (high-normal-stress-shear coexisting zone), Zone II (high-normal-stress zone) and Zone III (low-normal-stress zone) (Fig. 17(a)). Tovar-Valencia [54] found that fragmentation mainly occurs in Zones I and II during penetration. In Zone I, the sand is simultaneously subjected to both normal and shear stresses. Due to the relative movement of the surrounding sand, a shear band forms around the cone tip, within which significant fragmentation occurs [24. 54]. In Zone II, the normal stress is lower than in Zone I, and because sidewall friction is reduced, the sand in this zone is scarcely affected by shear stress. Previous studies have shown that particles are more prone to fragmentation under shear than under compression [55-57], consequently, the fragmentation rate in Zone II is lower than. In Region III, the normal stress further decreases, and it is primarily composed of unbroken sand and debris, which are migrated from Regions I and II through particle migration [54].

However, the stress state of the sand surrounding the ultrasonic CPT probe differs, and the can be further subdivided into: Zone I-1 (high normal stress-shear-vibration zone), Zone I-2 (high-normal-stress-shear zone), Zone II-1 (high-normal-stress-vibration zone), Zone II-2 (high-normal-stress zone) and Zone III (low-normal-stress zone) (Fig. 17(b)). In Zone I-1, sand beneath the cone tip experiences high normal and shear stress but also subjected to high-frequency ultrasonic vibration induced, leading to further breakage. Therefore, this zone is defined as the high-normal-stress-shear-vibration zone and corresponds to experimental regions C and D, where the overall sand breakage rate is relatively high. In Zone I-2, the stresses are similar to those in Zone I around the conventional CPT probe, sand is subjected to high normal and shear stress, but ultrasonic vibration is negligible. In Zone II-1, the sand is affected by both high normal stress and ultrasonic vibration, while the shear stress is negligible. The sand in Zone II-2 is subjected to high normal stress. Due to the limited propagation of ultrasonic waves, the stress conditions in Zone III are similar to those around a conventional CPT probe, normal stress is low and ultrasonic vibration is negligible.

Fig. 17The stress around the CPT probe during the penetration

The stress around the CPT probe during the penetration

a) Conventional CPT probe

The stress around the CPT probe during the penetration

b) Ultrasonic CPT probe

6. Conclusions

To overcome the limitation of the conventional CPT in dense sand layer, a novel ultrasonic CPT probe was developed based on ultrasonic drilling technology. To verify its effectiveness, a series of laboratory CPTs were conducted in sand at different relative densities, and the effects of relative density and ultrasonic vibration on penetration load and particle breakage were investigated. The key conclusions are as follows:

1) The ultrasonic CPT probe integrates an ultrasonic vibration unit into a dual-bridge CPT probe. The ultrasonic CPT probe consists of the cone tip, extension sleeve, ultrasonic transducer, ultrasonic generator, friction sleeve, metal rod and sensors. During penetration, ultrasonic vibration is transmitted to the cone tip, reducing the penetration load acting the probe.

2) The penetration load for both conventional and ultrasonic vibration penetration increases with the relative density of the sand. For conventional penetration, the load increases from 539 N at Dr= 30 % to 4410 N at Dr= 90 %; for ultrasonic-assisted penetration, it increases from 215.6 N at Dr= 30 % to 3635.8 N at Dr= 90 %. Ultrasonic vibration effectively reduces penetration load, thereby enhancing the CPT penetration capability. However, the load reduction is lower in dense sand compared to loose sand, the load reduction decreases from 60 % to 18 % when Dr increases from 30 % to 90 %.

3) Increasing the relative density of the sand facilitates the particle breakage around the cone tip, the closer to cone tip, the higher the breakage rate. However, the average breakage along the probe is limited. Image-J statistics show that the mean size of grains adhered to the cone tip decreases from 331.9±41.6 μm (Dr= 30 %) to 39.7±4.4 μm (Dr= 90 %), indicating progressively stronger and more uniform breaking at higher relative densities. Under ultrasonic vibration, the breakage index exceeds that of conventional CPT probe at Dr= 30 % and 45 %, but becomes lower when Dr 60 %.

4) The particle breakage mechanism under the ultrasonic-static coupling involves: (i) ultrahigh-frequency impacts between particles and the cone tip, (ii) a reduction in the local sand density of near the cone tip due to ultrahigh-frequency vibration, (iii) the effects of normal stress and shear stress around the probe.

5) The stress of the sand around the ultrasonic CPT probe can be further divided into: Zone I-1 (high-normal-stress-shear-vibration zone), Zone I-2 (high-normal-stress-shear zone), Zone II-1 (high-normal-stress-vibration zone), Zone II-2 (high normal stress zone) and Zone III (low normal stress zone).

Of course, there are still shortcomings in this study, and further works can be done:

1) Optimizing ultrasonic excitation parameters to further enhance resistance reduction in dense sands, such as frequency, amplitude, and power.

2) Particle Image Velocimetry (PIV), X-ray tomography and discrete element method (DEM) can be employed to quantitatively characterize the spatial evolution of the breakage zone and stress near the cone tip under vibration. Therefore, micro penetration mechanisms and macro penetration resistance can be combined together.

3) In-situ tests can be conducted in different soil types (silty sand and calcareous/coral sand) and varying saturation conditions to validate the engineering applicability of ultrasonic CPT probe.

References

  • Z. Wang, “The Chinese CPT and the prospect of penetration tests,” (in Chinese), Chinese Journal of Geotechnical Engineering, Vol. 22, No. 5, pp. 517–522, 2000.
  • Z. Cao et al., “A efficient optimization identification method for soil stratification based on cone penetration test and joint posterior distribution of variable dimensionality,” (in Chinese), Chinese Journal of Geotechnical Engineering, Vol. 47, No. 2, pp. 346–354, 2025, https://doi.org/10.11779/cjge20230715
  • L. Song, F. Y. Liu, and N. Li, “On mechanism of rotary cone penetration test,” (in Chinese), Rock and Soil Mechanics, Vol. 32, 2011, https://doi.org/10.16285/j.rsm.2011.s1.071
  • M. Martinelli, A. Tsetas, A. B. Fărăgău, A. Metrikine, and A. Tsouvalas, “Soil resistance during vibratory pile installation: Experimental findings from lab-scale tests,” Soil Dynamics and Earthquake Engineering, Vol. 199, p. 109692, Dec. 2025, https://doi.org/10.1016/j.soildyn.2025.109692
  • K. Rainer Massarsch, C. Wersäll, and B. H. Fellenius, “Vibratory driving of piles and sheet piles – state of practice,” Proceedings of the Institution of Civil Engineers – Geotechnical Engineering, Vol. 175, No. 1, pp. 31–48, Feb. 2022, https://doi.org/10.1680/jgeen.20.00127
  • S. Moriyasu, S.-I. Kobayashi, and T. Matsumoto, “Experimental study on friction fatigue of vibratory driven piles by in situ model tests,” Soils and Foundations, Vol. 58, No. 4, pp. 853–865, Aug. 2018, https://doi.org/10.1016/j.sandf.2018.03.010
  • Z. Qin, L. Chen, C. Song, and L. Sun, “Field tests to investigate the penetration rate of piles driven by vibratory installation,” Shock and Vibration, Vol. 2017, No. 1, pp. 1–10, Jan. 2017, https://doi.org/10.1155/2017/7236956
  • Z. Qin, L. Chen, C. Song, and J. Zhang, “Field tests for investigating the extraction rate of piles using a vibratory technique,” Journal of Shanghai Jiaotong University (Science), Vol. 23, No. 4, pp. 482–489, Aug. 2018, https://doi.org/10.1007/s12204-018-1971-4
  • N. Denies, J. Canou, J. Roux, and A. Holeyman, “Sphere penetration experiments in vertically vibrated sand,” in International Conferences on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics, Vol. 16, May 2010.
  • B. Darbois Texier, A. Ibarra, and F. Melo, “Low-resistive vibratory penetration in granular media,” PLOS ONE, Vol. 12, No. 4, p. e0175412, Apr. 2017, https://doi.org/10.1371/journal.pone.0175412
  • N. Denies, J. Canou, J.-N. Roux, and A. Holeyman, “Vibrocompaction properties of dry sand,” Canadian Geotechnical Journal, Vol. 51, No. 4, pp. 409–419, Apr. 2014, https://doi.org/10.1139/cgj-2012-0436
  • N. Denies and A. Holeyman, “Shear strength degradation of vibrated dry sand,” Soil Dynamics and Earthquake Engineering, Vol. 95, pp. 106–117, Apr. 2017, https://doi.org/10.1016/j.soildyn.2017.01.039
  • P. A. Johnson, H. Savage, M. Knuth, J. Gomberg, and C. Marone, “Effects of acoustic waves on stick-slip in granular media and implications for earthquakes,” Nature, Vol. 451, No. 7174, pp. 57–60, Jan. 2008, https://doi.org/10.1038/nature06440
  • P. A. Johnson et al., “Dynamically triggered slip leading to sustained fault gouge weakening under laboratory shear conditions,” Geophysical Research Letters, Vol. 43, No. 4, pp. 1559–1565, Feb. 2016, https://doi.org/10.1002/2015gl067056
  • H. Jiang, L. Fu. S. Zhou, P. Guo, and W. Ye, “Experimental study on critical state strength characteristics of granular material-structure interface under high-frequency vibration,” (in Chinese), Rock and Soil Mechanics, No. 3, pp. 810–820, Mar. 2023, https://doi.org/10.16285/j.rsm.2022.0496
  • W. Ye et al., “Experimental study on dynamic characteristics of granular materials under axial high-frequency vibration,” Acta Geotechnica, Vol. 17, No. 8, pp. 3211–3227, Feb. 2022, https://doi.org/10.1007/s11440-022-01474-z
  • X. Wang, X. Wang, J. Wang, and Z. Tian, “Feasibility study and prospects of rock fragmentation using ultrasonic vibration excitation,” Applied Sciences, Vol. 10, No. 17, p. 5868, Aug. 2020, https://doi.org/10.3390/app10175868
  • M. Wiercigroch, J. Wojewoda, and A. M. Krivtsov, “Dynamics of ultrasonic percussive drilling of hard rocks,” Journal of Sound and Vibration, Vol. 280, No. 3-5, pp. 739–757, Feb. 2005, https://doi.org/10.1016/j.jsv.2003.12.045
  • J. Huang, H. Wang, G. Ji, F. Zhao, Q. Ming, and Y. Hao, “The rock breaking mechanism of ultrasonic high frequency rotary-percussive drilling technology,” (in Chinese), Petroleum Drilling Techniques, No. 4, pp. 23–29, Mar. 2018, https://doi.org/10.11911/syztjs.2018097
  • Z. Lu, J. Zheng, Z. Jiang, and F. Zhao, “An experimen tal study on rock breaking efficiency with ultrasonic high-frequency rotary-percussive drilling technology,” (in Chinese), Petroleum Drilling Techniques, No. 1, pp. 20–15, Mar. 2021, https://doi.org/10.11911/syztjs.2020126
  • Y. Li, L. Yang, X. Zhang, Y. Pang, and W. Bai, “Research on rotational-impact ultrasonic sampling technology for simulated lunar regolith permafrost,” (in Chinese), Journal of Space Science and Experiment, No. 3, pp. 35–45, Mar. 2025, https://doi.org/10.19963/j.cnki.2097-4302.2025.03.004
  • W. Mao, H. Hamaguchi, and J. Koseki, “Discrimination of particle breakage below pile tip after model pile penetration in sand using image analysis,” International Journal of Geomechanics, Vol. 20, No. 1, p. 04019, Jan. 2020, https://doi.org/10.1061/(asce)gm.1943-5622.0001535
  • D. J. White and M. D. Bolton, “Displacement and strain paths during plane-strain model pile installation in sand,” Géotechnique, Vol. 54, No. 6, pp. 375–397, Aug. 2004, https://doi.org/10.1680/geot.2004.54.6.375
  • Z. X. Yang, R. J. Jardine, B. T. Zhu, P. Foray, and C. H. C. Tsuha, “Sand grain crushing and interface shearing during displacement pile installation in sand,” Géotechnique, Vol. 60, No. 6, pp. 469–482, Jun. 2010, https://doi.org/10.1680/geot.2010.60.6.469
  • I. M. Silva, G. Combe, P. Foray, F. Flin, and B. Lesaffre, “Postmortem analysis of sand grain crushing from pile interface using X-ray tomography,” in Powders and Grains 2013: Proceedings of the 7th International Conference on Micromechanics of Granular Media, Vol. 1542, No. 1, pp. 297–300, Jan. 2013, https://doi.org/10.1063/1.4811926
  • J. Doreau-Malioche et al., “Characterising processes at sand-pile interface using digital image analysis and X-ray CT,” Géotechnique Letters, Vol. 9, No. 4, pp. 254–262, Dec. 2019, https://doi.org/10.1680/jgele.18.00232
  • M. I. Arshad, F. S. Tehrani, M. Prezzi, and R. Salgado, “Experimental study of cone penetration in silica sand using digital image correlation,” Géotechnique, Vol. 64, No. 7, pp. 551–569, Jun. 2014, https://doi.org/10.1680/geot.13.p.179
  • Y. Peng, Z.-Y. Yin, and X. Ding, “Micromechanical analysis of the particle corner breakage effect on pile load performance in coral sand,” Acta Geotechnica, Vol. 18, No. 12, pp. 6353–6370, Jun. 2023, https://doi.org/10.1007/s11440-023-01975-5
  • W. Mao, S. Aoyama, and I. Towhata, “A study on particle breakage behavior during pile penetration process using acoustic emission source location,” Geoscience Frontiers, Vol. 11, No. 2, pp. 413–427, Mar. 2020, https://doi.org/10.1016/j.gsf.2019.04.006
  • W. Mao, J. Chen, and L. Pan, “Acoustic emission characteristics during pile penetration process in coral sand,” (in Chinese), Journal of Engineering Geology, Vol. 32, No. 5, pp. 1872–1879, Oct. 2024, https://doi.org/10.13544/j.cnki.jeg.2024-0370
  • W. Mao, S. Aoyama, and I. Towhata, “Feasibility study of using acoustic emission signals for investigation of pile spacing effect on group pile behavior,” Applied Acoustics, Vol. 139, pp. 189–202, Oct. 2018, https://doi.org/10.1016/j.apacoust.2018.05.001
  • W. Mao, Y. Yang, W. Lin, S. Aoyama, and I. Towhata, “High frequency acoustic emissions observed during model pile penetration in sand and implications for particle breakage behavior,” International Journal of Geomechanics, Vol. 18, No. 11, p. 04018, Nov. 2018, https://doi.org/10.1061/(asce)gm.1943-5622.0001287
  • M. J. Jiang, D. Harris, and H. H. Zhu, “Future continuum models for granular materials in penetration analyses,” Granular Matter, Vol. 9, No. 1-2, pp. 97–108, Nov. 2006, https://doi.org/10.1007/s10035-006-0026-y
  • O. Falagush, G. R. Mcdowell, H. S. Yu, and J. P. de Bono, “Discrete element modelling and cavity expansion analysis of cone penetration testing,” Granular Matter, Vol. 17, No. 4, pp. 483–495, Jun. 2015, https://doi.org/10.1007/s10035-015-0570-4
  • A. Sadrekarimi and S. Hashemi, “Effect of boundary conditions on CPT calibration chamber tests,” International Journal of Geomechanics, Vol. 25, No. 3, p. 04025, Mar. 2025, https://doi.org/10.1061/ijgnai.gmeng-10376
  • Y. Ge, A. Zhou, M. Nazem, and Y. Deng, “Numerical simulation of cone penetration test by using CFD-DEM coupled analysis,” Acta Geotechnica, Vol. 19, No. 11, pp. 7635–7653, Jul. 2024, https://doi.org/10.1007/s11440-024-02369-x
  • X. Huang et al., “Effects of strain rate and softening on the cone-tip resistance of CPT in soft clay,” Marine Georesources and Geotechnology, Vol. 43, No. 10, pp. 1876–1887, Oct. 2025, https://doi.org/10.1080/1064119x.2024.2439047
  • Y. Song, X. Gu, J. Hu, and X. Z., “DEM-FDM coupled simulation of cone penetration tests in a virtual calibration chamber with sand,” Chinese Journal of Geotechnical Engineering, Vol. 47, No. 6, pp. 1249–1258, 2025, https://doi.org/10.11779/cjge20240172
  • J. Gong et al., “Exploring the microscale mechanics of shear strength weakening in granular soils induced by vibrations: A DEM study,” Powder Technology, Vol. 456, p. 120819, Apr. 2025, https://doi.org/10.1016/j.powtec.2025.120819
  • H. Jiang, S. Zhou, L. Fu, Y. Shan, N. Dai, and P. Guo, “Vibro-induced weakening of interface friction between granular materials and textured surfaces: An experimental study with a modified direct shear apparatus,” Soil Dynamics and Earthquake Engineering, Vol. 158, p. 107289, Jul. 2022, https://doi.org/10.1016/j.soildyn.2022.107289
  • H. Lastakowski, J.-C. Géminard, and V. Vidal, “Granular friction: triggering large events with small vibrations,” Scientific Reports, Vol. 5, No. 1, p. 13455, Sep. 2015, https://doi.org/10.1038/srep13455
  • K. A. Taslagyan, D. H. Chan, and N. R. Morgenstern, “Vibrational fluidization of granular media,” International Journal of Geomechanics, Vol. 16, No. 3, p. 04015, Jun. 2016, https://doi.org/10.1061/(asce)gm.1943-5622.0000568
  • E. Bauer, “Modelling of the pressure and density sensitive behaviour of sand within the framework of hypoplasticity,” Task Quarterly, Vol. 4, No. 3, pp. 367–387, Sep. 2000.
  • C. M. Bingham, D. A. Stone, N. Schofield, D. Howe, and D. Peel, “Amplitude and frequency control of a vibratory pile driver,” IEEE Transactions on Industrial Electronics, Vol. 47, No. 3, pp. 623–631, Jun. 2000, https://doi.org/10.1109/41.847903
  • B. Yuan et al., “Effect of sand relative density on response of a laterally loaded pile and sand deformation,” Journal of Chemistry, Vol. 2015, No. 1, p. 89121, Jul. 2015, https://doi.org/10.1155/2015/891212
  • L. Fu, S. Zhou, P. Guo, S. Wang, and Z. Luo, “Induced force chain anisotropy of cohesionless granular materials during biaxial compression,” Granular Matter, Vol. 21, No. 3, p. 52, Jun. 2019, https://doi.org/10.1007/s10035-019-0899-1
  • Y. Guo, C. Zhao, V. Markine, C. Shi, G. Jing, and W. Zhai, “Discrete element modelling of railway ballast performance considering particle shape and rolling resistance,” Railway Engineering Science, Vol. 28, No. 4, pp. 382–407, Aug. 2020, https://doi.org/10.1007/s40534-020-00216-9
  • R. J. Marsal, “Large scale testing of rockfill materials,” Journal of the Soil Mechanics and Foundations Division, Vol. 93, No. 2, pp. 27–43, Mar. 1967, https://doi.org/10.1061/jsfeaq.0000958
  • H. Z. Khaleel et al., “Measurement enhancement of ultrasonic sensor using pelican optimization algorithm for robotic application,” Indonesian Journal of Science and Technology, Vol. 9, No. 1, pp. 145–162, Nov. 2023, https://doi.org/10.17509/ijost.v9i1.64843
  • A. Taher Azar et al., “Optimized model predictive control for unmanned aerial vehicles with sensor uncertainties,” in 2022 International Conference on Control, Automation and Diagnosis (ICCAD), pp. 1–6, Jul. 2022, https://doi.org/10.1109/iccad55197.2022.9869421
  • F. Liu, L. Zhong, and Y. Zhong, “Effects of different particle arrangements on mechanical properties of sand,” (in Chinese), Chinese Journal of Geotechnical Engineering, Vol. 43, No. S1, pp. 155–160, Jul. 2021.
  • D. Hao, C. Yue, R. Chen, J. R., and F. C., “Shear characteristics and stress-dilation relation of medium sand under normal to high pressures.,” (in Chinese), Chinese Journal of Geotechnical Engineering, Vol. 42, No. 4, pp. 765–772, Apr. 2020.
  • D. Su, Z. Wu, G. Lei, and M. Zhu, “Numerical study on the installation effect of a jacked pile in sands on the pile vertical bearing capacities,” Computers and Geotechnics, Vol. 145, p. 104690, May 2022, https://doi.org/10.1016/j.compgeo.2022.104690
  • R. D. Tovar-Valencia, A. Galvis-Castro, R. Salgado, M. Prezzi, and D. Fridman, “Experimental measurement of particle crushing around model piles jacked in a calibration chamber,” Acta Geotechnica, Vol. 18, No. 3, pp. 1331–1351, Sep. 2022, https://doi.org/10.1007/s11440-022-01681-8
  • G. Ma, W. Zhou, X.-L. Chang, T.-T. Ng, and L.-F. Yang, “Formation of shear bands in crushable and irregularly shaped granular materials and the associated microstructural evolution,” Powder Technology, Vol. 301, pp. 118–130, Nov. 2016, https://doi.org/10.1016/j.powtec.2016.05.068
  • G. Ma, R. A. Regueiro, W. Zhou, Q. Wang, and J. Liu, “Role of particle crushing on particle kinematics and shear banding in granular materials,” Acta Geotechnica, Vol. 13, No. 3, pp. 601–618, Jan. 2018, https://doi.org/10.1007/s11440-017-0621-6
  • F. Yu, “Particle breakage and the drained shear behavior of sands,” International Journal of Geomechanics, Vol. 17, No. 8, p. 04017, Aug. 2017, https://doi.org/10.1061/(asce)gm.1943-5622.0000919

About this article

Received
November 10, 2025
Accepted
January 28, 2026
Published
February 15, 2026
SUBJECTS
Mechanical vibrations and applications
Keywords
CPT
ultrasonic vibration
particle breakage
penetration load
Acknowledgements

The study was supported by the National Natural Science Foundation of China (U2244228), Major Project of the National Natural Science Foundation of China (41831290), the Central Government's Guidance Fund Projects for Local Science and Technology Development (2024ZY01041, 2025ZY01028), the Zhejiang Provincial Natural Science Foundation of China (LTGS23E040001, LY22D020001), Science and Technology Plan Project of Shaoxing (2024A11014) and Construction Research Project of Zhejiang Provincial Department of Housing and Urban Rural Development (2023K263), which is gratefully acknowledged.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author Contributions

Faquan Wu, Fang Zhang, Bolong liu, Jie Wu: conceptualization. Weihong Chen, Zhongli Xu: methodology. Zhehong Shen, Yuanlong Li, Yuqi Lang: writing-original draft. Bolong Liu, Zhongqing Chen, Zihao Zhu, Enjie Xu: writing-review and editing. All authors read and approved the manuscript.

Conflict of interest

The authors declare that they have no conflict of interest.