Journal of Rock Mechanics and Geotechnical Engineering最新文献

Sand typically exhibits anisotropic internal structure which may significantly influence its mechanical behavior. The material point method (MPM) can eliminate mesh distortion and thus is suitable for investigating geotechnical problems with large deformation. In this study, an advanced anisotropic critical state theory (ACST)-based soil model is implemented in MPM to study the response of strip footing resting on anisotropic sand. The capability of the model is verified by simulating several element tests and strip footing tests with different soil densities and fabric bedding plane orientations. For the footing problem with a vertical load, as the fabric bedding plane orientation increases, the bearing capacity decreases and its corresponding settlement increases. The failure pattern becomes asymmetrical when the bedding plane orientation or the loading direction is inclined. A comparison between the simulation results predicted by the anisotropic and isotropic models is made, which demonstrates that neglecting the fabric anisotropy may lead to the overestimation of the bearing capacity.

The subgrade soil is generally in saturated or unsaturated condition. To analyze complex thermo-hydro-mechanical-chemical (THMC) behaviors of subgrade, it is essential to determine the soil–water characteristic curve (SWCC) that represents the relationship between matric suction and moisture content. In this study, a full-automatic rapid stress-dependent SWCC pressure-plate extractor was developed. Then, the influences of overburden stress and degree of compaction on the SWCC of subgrade soil such as high liquid limit silt (MH) and low liquid limit clay (CL) were analyzed. Accordingly, a new model taking into account the influences of overburden stress and degree of compaction based on the well-known Van Genuchten (VG) SWCC fitting model was presented and validated. The results show that with the increase of the degree of compaction and overburden stress, the saturated moisture content of subgrade soil decreases, while the air-entry value increases and the transition section curve becomes flat. The influences of the degree of compaction and overburden stress on the SWCC of MH is greater than that of CL. Meanwhile, there was a satisfactory agreement between the prediction and measurement, indicating a good performance of the new model for predicting the SWCC.

The threshold values of CO₂ gas stripped off membranous residual oil from the pore walls are not clear under different temperatures, pressures and wettability conditions. The extent to which temperature, pressure and wettability influence CO₂ flooding for enhancing the recovery of residual oil in membranous formations also remains uncertain. Therefore, further quantitative characterization is entailed. In this study, the molecular dynamics method was employed to explore CO₂ flooding under different temperatures, pressures and wettability conditions, aiming to enhance the production of membranous residual oil. The results reveal that the interaction energy between CO₂, decane molecules and pore walls exhibits a decrease with increasing temperature and an increase with increasing pressure, respectively, in distinct wettability scenarios. When the temperature was at or below 363 K and the pressure was not lower than 40 MPa, CO₂ gas could detach the membranous residual oil from the pore walls in the water-wet systems. When the temperature was equal to 363 K and the pressure remained under 40 MPa, or the temperature surpassed 363 K, CO₂ gas failed to detach the membranous residual oil from the pore walls in the water-wet systems. For the mixed-wet and oil-wet systems, CO₂ molecules could not detach the membranous residual oil from the pore walls. The hierarchy of influence regarding temperature, pressure and wettability on the competitive adsorption capacity of CO₂ and decane molecules on the pore walls emerged as follows: wettability > temperature > pressure. The findings of this study offer valuable insights into the application of CO₂ gas flooding for the exploitation of membranous residual oil on pore walls.

Defects in cast-in-situ piles have an adverse impact on load transfer at the pile‒soil interface and pile bearing capacity. In recent years, thermal integrity profiling (TIP) has been developed to measure temperature profiles of cast-in-situ piles, enabling the detection of structural defects or anomalies at the early stage of construction. However, using this integrity testing method to evaluate potential defects in cast-in-situ piles requires a comprehensive understanding of the mechanism of hydration heat transfer from piles to surrounding soils. In this study, small-scale model tests were conducted in laboratory to investigate the performance of TIP in detecting pile integrity. Fiber-optic distributed temperature sensing (DTS) technology was used to monitor detailed temperature variations along model piles in sand. Additionally, sensors were installed in sand to measure water content and matric suction. An interpretation method against available DTS-based thermal profiles was proposed to reveal the potential defective regions. It shows that the temperature difference between normal and defective piles is more obvious in wet sand. In addition, there is a critical zone of water migration in sand due to the water absorption behavior of cement and temperature transfer-induced water migration in the early-age concrete setting. These findings could provide important insight into the improvement of the TIP testing method for field applications.

In this study, a three-dimensional (3D) finite element modelling (FEM) analysis is carried out to investigate the effects of soil spatial variability on the response of retaining walls and an adjacent box culvert due to a braced excavation. The spatial variability of soil stiffness is modelled using a variogram and calibrated by high-quality experimental data. Multiple random field samples (RFSs) of soil stiffness are generated using geostatistical analysis and mapped onto a finite element mesh for stochastic analysis of excavation-induced structural responses by Monte Carlo simulation. It is found that the spatial variability of soil stiffness can be described by an exponential variogram, and the associated vertical correlation length is varied from 1.3 m to 1.6 m. It also reveals that the spatial variability of soil stiffness has a significant effect on the variations of retaining wall deflections and box culvert settlements. The ignorance of spatial variability in 3D FEM can result in an underestimation of lateral wall deflections and culvert settlements. Thus, the stochastic structural responses obtained from the 3D analysis could serve as an effective aid for probabilistic design and analysis of excavations.

To address the problems of strain localization, the exact Mohr-Coulomb (MC) model is used based on second-order cone programming (mpcFEM-SOCP) in the framework of micropolar continuum finite element method. Using the uniaxial compression test, we focused on the earth pressure problem of rigid wall segment involving non-associated plasticity. The numerical results reveal that when mpcFEM-SOCP is applied, the problems of mesh dependency can be effectively addressed. For geotechnical strain localization analysis involving non-associated MC plasticity, mpcFEM-SOCP in conjunction with the pseudo-time discrete scheme can improve the numerical stability and avoid the unreasonable softening issue in the pressure-displacement curves, which may be encountered in the conventional FEM. It also shows that the pressure-displacement responses calculated by mpcFEM-SOCP with the pseudo-time discrete scheme are higher than those calculated by mpcFEM-SOCP with the Davis scheme. The inclination angle of shear band predicted by mpcFEM-SOCP with the pseudo-time discrete scheme agrees well with the theoretical solution of non-associated MC plasticity.

To improve the accuracy of indirect tensile strength for a transversely isotropic rock in the Brazilian test, this study considered the three-dimensional (3D) deformation and the nonlinear stress–strain relationship. A parametric study of a numerical Brazilian test was performed for a general range of elastic constants, revealing that the 3D modeling evaluated the indirect tensile strength up to 40% higher than the plane stress modeling. For the actual Asan gneiss, the 3D model evaluated the indirect tensile strength up to 10% higher and slightly enhanced the accuracy of deformation estimation compared with the plane stress model. The nonlinearity in stress–strain curve of Asan gneiss under uniaxial compression was then considered, such that the evaluated indirect tensile strength was affected by up to 10% and its anisotropy agreed well with the physical intuition. The estimation of deformation was significantly enhanced. The further validation on the nonlinear model is expected as future research.

Reef limestone is a biogenic sedimentary rock widely distributed in coral reef areas, acting as an important foundation for coast construction. Due to its special biogenic origin, reef limestone is different from conventional rocks both in terms of rock structure and mechanical properties. In this study, mesoscale uniaxial compression experiments with five different loading directions were conducted on two kinds of reef limestones from the Maldives Islands and the South China Sea, respectively. The real-time high-resolution videos and images of failure processes were recorded simultaneously to investigate the fracture evolution and fracture surface roughness of reef limestones. It demonstrated that the reef limestones belonged to extremely soft to soft rocks, and their uniaxial compressive strength (UCS) values fluctuated with high discreteness. The mesoscale mechanical properties of reef limestones were highly anisotropic and mainly controlled by pore structure. The occurrence of dissolution pores in reef limestone tended to intensify mechanical anisotropy. With the integration of the fracture initiation and propagation features of reef limestones, it is supposed that the intrinsic mechanism of anisotropy was probably attributed to the differences in coral growth direction and dissolution. Furthermore, the quantified fracture surface roughness was revealed to have a good consistency with density and UCS for the reef limestones from the South China Sea. The findings are helpful for providing theoretical and experimental references for engineering construction in coral reef areas.