Earthquake-induced soil liquefaction causes ground and foundation failures, and it challenges the scientific community to explore the liquefaction problem in deep deposit under strong shaking. Due to the capacity limitation of physical modelling facility, it is difficult to reproduce soil liquefaction response of deep sand ground by centrifuge shaking table test. To address this problem, a suite of centrifuge model tests were conducted with the aid of Iai's Type III generalized scaling law (i.e., GSL) to observe the liquefaction response of deep sand ground, where Models 1 and 2 were used to validate the GSL and Model 3 with the validated GSL stands for the deep sand ground with prototype depth of 80 m. The test results of Models 1 and 2 indicate that GSL generally performs well for small-strain shear modulus, nonlinear dynamic response of acceleration and the generation of excess pore water pressure, but leaves considerable errors for post-shaking dissipation process and ground settlement with large plastic strain. The test results of Model 3 indicate that liquefaction is also possible in depth of 30-40 m under shaking event of PBA = 0.4 g and Mw = 7.5. For deeper depth without triggering of liquefaction, a depthdependent power function relationship between the peak excess pore water pressure and Arias intensity has been established. The test results also revealed that consolidation and earthquake shaking history contribute to the development of soil anisotropy in a deep ground, leading to a continuous increase of anisotropy degree, which could be evaluated using the small-strain shear moduli in different stress planes under orthogonal stress conditions.

Seepage problems in half-space domains are crucial in hydrology, environmental, and civil engineering, involving groundwater flow, pollutant transport, and structural stability. Typical examples include seepage through dam foundations, coastal aquifers, and levees under seepage forces, requiring accurate numerical modeling. However, existing methods face challenges in handling complex geometries, heterogeneous media, and anisotropic properties, particularly in multi-domain half-spaces. This study addresses these challenges by extending the modified scaled boundary finite element method (SBFEM) and using this method to explore steady seepage problems in complex half-space domain. In the modified SBFEM framework, segmented straight lines or curves, parallel to the far-field infinite boundary, are introduced as scaling lines, with a one-dimensional discretization applied to them, thereby reducing computational costs.Then the weighted residual method is applied to obtain the modified SBFEM governing equations and boundary conditions of steady-state seepage problem according to the Laplace diffusion equation and Darcy's law. Furthermore, the steady seepage matrix at infinity is obtained by solving the eigenvalue problem of Schur decomposition and then the 4th-order Runge-Kutta algorithm is used to iteratively solve until the seepage matrix at the boundary lines is reached. Comparisons between the present numerical results and solutions available in the published work have been conducted to demonstrate the efficiency and accuracy of this method. At the same time, the influences of the geometric parameters and complex half-space domain on the seepage flow characteristics in complex half-space domain are investigated in detail.

The dynamic response of piles is a fundamental issue that significantly affects the performance of pile foundations under vertical cyclic loading, yet it has been insufficiently explored due to the limitations of experimental methods. To address this gap, a hydraulic loading device was developed for centrifuge tests, capable of applying loads up to 2.5 kN and 360 Hz. This device could simulate the primary loading conditions encountered in engineering applications, such as those in transportation and power machinery, even after the amplification of the dynamic frequency for centrifuge tests. Furthermore, a design approach for model piles that considers stress wave propagation in pile body and pile-soil dynamic interaction was proposed. Based on the device and approach, centrifuge comparison tests were conducted at 20 g and 30 g, which correspond to the same prototype. The preliminary results confirmed static similarity with only a 1.25% deviation in ultimate bearing capacities at the prototype scale. Cyclic loading tests, conducted under various loading conditions that were identical at the prototype scale, indicated that dynamic displacement, cumulative settlement, and axial forces at different burial depths adhered the dynamic similarity of centrifuge tests. These visible phenomena effectively indicate the rationality of centrifuge tests in studying pile-soil interaction and provide a benchmark for using centrifuge tests to investigate soil-structure dynamic interactions.

In subsurface projects where the host rock is of low permeability, fractures play an important role in fluid circulation. Both the geometrical and mechanical properties of the fracture are relevant to the permeability of the fracture. To evaluate this relationship, we numerically generated self-affine fractures reproducing the scaling relationship of the power spectral density (PSD) of the measured fracture surfaces. The fractures were then subjected to a uniform and stepwise increase in normal stress. A fast Fourier transform (FFT)-based elastic contact model was used to simulate the fracture closure. The evolution of fracture contact area, fracture closure, and fracture normal stiffness were determined throughout the whole process. In addition, the fracture permeability at each step was calculated by the local cubic law (LCL). The influences of roughness exponent and correlation length on the fracture hydraulic and mechanical behaviors were investigated. Based on the power law of normal stiffness versus normal stress, the corrected cubic law and the linear relationship between fracture closure and mechanical aperture were obtained from numerical modeling of a set of fractures. Then, we derived a fracture normal stiffness-permeability equation which incorporates fracture geometric parameters such as the root-mean-square (RMS), roughness exponent, and correlation length, which can describe the fracture flow under an effective medium regime and a percolation regime. Finally, we interpreted the flow transition behavior from the effective medium regime to the percolation regime during fracture closure with the established stiffness-permeability function. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

The underground concrete silo, designed as a hollow cylinder with a large aspect ratio and thin walls, is highly susceptible to failure caused by intentional or accidental soil explosions. To enhance its protection, this study investigates the dynamic tensile responses and failure mechanisms of underground concrete silos subjected to high-yield soil explosions. The concept of nominal crack width is proposed to quantitatively describe the degree of overall bending-induced tensile responses and failure of the concrete silo. The influences of explosive weights, standoff distances, and the aspect ratios and thicknesses of the underground concrete silo are quantitatively explored first. On this basis, a dimensionless number combining these major influencing factors is derived using dimensional analysis. The derived dimensionless number has a clear physical meaning, reflecting three aspects: the inertia of the blast loading, the resistance ability of concrete material to bending responses and failure, and the resistance ability of silo structure to bending responses and failure. The results demonstrate that the proposed dimensionless number effectively correlates with the overall bending-induced tensile responses and failure of silo structures across various geometries and explosion scenarios, exhibiting a good linear relation with the dimensionless nominal crack width of the concrete silo. With its solid physical foundation, the dimensionless number offers practical applications in scaling analysis and fast damage assessment. Specific examples of these applications are presented and discussed in this study.

This paper employs the discrete element method (DEM) to simulate the nanoindentation creep of calcium- silicate-hydrate (C-S-H), focusing on indentation deformation, particle interactions, and stress transmission paths. The Rate Process Theory (RPT), previously utilized in the creep modeling of cohesive soils and other granular materials, is proposed to simulate C-S-H creep. Due to the nanometer size of C-S-H particles, the critical time step in DEM simulations is very small. Therefore, a time-scaling algorithm is used to match the DEM simulation time with the physical time in laboratory tests, accelerating the simulation time by a factor of 1 x 108. C-S-H particle assemblies with specific packing densities are generated using Particle Flow Code (PFC3D, version 5.0), with coordination numbers and cohesion forces controlled by the stress-servo of PFC walls. Virtual nano- indentations using a Berkovich indenter are conducted on C-S-H particle assemblies with three different packing densities (0.74, 0.64, and 0.58), followed by parameters calibration. Results show that the DEM + RPT method can capture the scaling relations between the indentation modulus, hardness, and contact creep modulus of C-S- H particle assemblies and the packing density. Furthermore, DEM simulations reveal particle rearrangement under Berkovich and flat-tip indenters, highlighting that different indenter types lead to distinct creep kinetics in C-S-H, with the Berkovich indenters experimentally capturing long-term creep and flat-tip indenters measuring short-term creep.

Remotely sensed top-of-the-canopy (TOC) SIF is highly impacted by non-physiological structural and environmental factors that are confounding the photosystems' emitted SIF signal. Our proposed method for scaling TOC SIF down to photosystems' (PSI and PSII) level uses a three-dimensional (3D) modeling approach, capable of accounting physically for the main confounding factors, i.e., SIF scattering and reabsorption within a leaf, by canopy structures, and by the soil beneath. Here, we propose a novel SIF downscaling method that separates the structural component from the functional physiological component of TOC SIF signal by using the 3D Discrete Anisotropic Radiative Transfer (DART) model coupled with the leaf-level fluorescence model Fluspect-CX, and estimates the Fluorescence Quantum Efficiency (FQE) at photosystem level. The method was first applied on in- situ diurnal measurements acquired at the top of the canopy of an alfalfa crop with a near-distance point- measuring FloX system. The retrieved photosystem-level FQE diurnal courses correlated significantly with photosynthetic yield of PSII measured by an active leaf florescence instrument MiniPAM (R = 0.87, R2 = 0.76 before and R =-0.82, R2 = 0.67 after 2.00 pm local time). Diurnal FQE trends of both photosystems jointly were descending from late morning 9.00 am till afternoon 4.00 pm. A slight late-afternoon increase, observed for three days between 4.00 and 7.00 pm, could be attributed to an increase in FQE of PSI that was retrieved separately from PSII. The method was subsequently extended and applied to airborne SIF images acquired with the HyPlant imaging spectrometer over the same alfalfa field. While the input canopy SIF radiance computed by two different methods, i) a spectral fitting method (SFM) and ii) a spectral fitting method neural network (SFMNN), produce broad and irregularly shaped (skewed) histograms (spatial coefficients of variation: CV = 29-35 % and 14-20 %, respectively), the retrieved HyPlant per-pixel FQE estimates formed significantly narrower and regularly bell- shaped near-Gaussian histograms (CV = 27-34 % and 14-17 %, respectively). The achieved spatial homogeneity of resulting FQE maps confirms successful removal of the TOC SIF radiance confounding impacts. Since our method is based on direct matching of measured and physically modelled canopy SIF radiance, simulated by 3D radiative transfer, it is versatile and transferable to other canopy architectures, including structurally complex canopies such as forest stands.

Two new structure-specific scalar intensity measures for plane steel frames under far-field earthquakes are proposed. These intensity measures of the spectral acceleration and spectral displacement type are multi-modal as they take into account the effect of the first four natural periods and multi-level as they are defined for four performance levels and consider inelasticity and period elongation up to the collapse prevention level. This is accomplished with the aid of the equivalent modal damping ratios of a structure previously developed by the authors for performance-based seismic design purposes. These modal damping ratios are period, soil type and deformation dependent and associate the equivalent linear structure to the original nonlinear one. The proposed intensity measures are conceptually simple, elegant and include all the aforementioned features in a rational way without artificially combining terms, defining period ranges and adding coefficients to be determined by optimization procedures as it is the case for all the existing measures, which try to take into account more than one mode and inelasticity. Comparison of the proposed intensity measures against some of the most popular ones existing in the literature, with respect to efficiency (beta), practicality (b), proficiency (zeta), sufficiency in terms of seismic magnitude (M) and source-to-site distance (R), scaling robustness and the range of their values at any damage or performance level demonstrates their very good performance as indicators of the destructive power of an earthquake.

The conventional similarity theory derived from dimensional analysis struggles with the well-known issue of non-scalability of material strain-rate effects between scaled models and prototypes. This limitation has significantly hindered the application of scaled model tests, particularly small-scale centrifugal model tests, in the study of structures against blast loading. To overcome this challenge, this study proposes a rate-dependent similarity theory for scaling the dynamic tensile responses and failure of large-scale underground concrete silos (46 m in height) subjected to large-yield soil explosions. The proposed theory includes a correction method derived from a verified dimensionless number, Dcs, which accurately reflects the overall bending-induced tensile response and failure mechanism of concrete silos. The correction strategy involves maintaining an equal Dcs between the scaled model and the prototype by adjusting the explosive weight and the concrete's static tensile strength in the scaled model to account for differences in strain-rate effects. To verify the theory, a series of geometrically similar silo models with scaling factors beta = 1, 1/2, 1/5, 1/10, 1/20, 1/50, and 1/100 were designed. High-fidelity numerical simulations were performed using a fully coupled numerical model encompassing the explosive-soil-silo system. The results demonstrate that, with the conventional dimensional analysisbased similarity theory, the tensile damage and failure of the scaled silo models differ significantly from those of the prototype. However, with the proposed rate-dependent similarity theory, the failure patterns of the silo models with beta = 1 similar to 1/100 are almost identical, indicating that the proposed theory can effectively address the troublesome issue of dissimilar material strain-rate effects between scaled models and prototypes. This similarity theory offers a solid theoretical foundation for designing scaled models that accurately reflect prototype behavior, thereby advancing the application of scaled model tests in the study of structures against blast loading.

Low-saturation liquid-containing granular materials are commonly encountered in both natural and industrial settings, where interstitial liquids significantly affect the motion of particles, while particle size polydispersity plays a crucial role in determining the level of system cohesion. In this study, the collapse of wet polydisperse granular columns is numerically investigated based on the developed discrete element model, with corresponding dam-break experiments performed to validate our numerical model and methodology. The dependence of the dynamics and flow mobility on particle size distribution is primarily examined, and the underlying mechanisms are also explored by analyzing particle path lengths and average fidelity. Building upon the effective Bond number proposed using the mixing theory, a macroscopic cohesion parameter at the material scale is defined by considering the dependence of the collapse on the system size effect. The relevance of this cohesion parameter in describing different wet polydisperse granular collapses is further validated based on our designed experimental tests and DEM simulations. The approach of constructing the cohesion parameters at different scales can be extended to characterize cohesion effects in more complex wet polydisperse granular flows and describe their associated rheological behaviors.