This paper establishes a novel full-process numerical simulation framework for analyzing the 3D seismic response of mountain tunnels induced by active faults. The framework employs a two-step approach to achieve wavefield transmission through equivalent seismic load: first, a highly efficient and accurate FMIBEM (Fast multipole indirect boundary element method) is used for large-scale 3D numerical simulations at the regional scale to generate broadband ground motions (1-5 Hz) for specific sites; subsequently, using the FEM (Finite element method), a refined simulation of the plastic deformation of surrounding rock and the elastoplastic behavior of the tunnel structure was conducted at the engineering scale. The accuracy of the framework has been validated. To further demonstrate its effectiveness, the framework is applied to analyze the impact of different fault movement mechanisms on the damage to mountain tunnels based on a scenario earthquake (Mw 6.7). By introducing tunnel structure damage classification and corresponding damage indicators, the structural damage levels of tunnels subjected to active fault movements are quantitatively evaluated. The findings demonstrate that the framework successfully simulates the entire process, from fault rupture and terrain amplification to the seismic response of tunnel structures. Furthermore, the severity of tunnel damage caused by different fault types is ranked as follows: reverse fault > normal fault > strike-slip fault.

This study presents a hierarchical multiscale approach that combines the finite-element method (FEM) and the discrete-element method (DEM) to investigate tunneling-induced ground responses in coarse-grained soils. The approach considers both particle-scale physical characteristics and engineering-scale boundary value problems (BVPs) simultaneously, accurately reproducing typical tunneling-induced mechanical responses in coarsegrained soils, including soil arching and ground movement characteristics observed in laboratory tests and engineering practice. The study also unveils particle-scale mechanisms responsible for the evolution of soil arching through the underlying DEM-based RVEs. The results show that the rearrangement of microstructures and the deflection of strong contact force chains drive the rotation of macroscopic principal stress and the formation of soil arch. The microscopic fabric anisotropy direction can serve as a quantitative indicator for characterizing soil arching zones. Moreover, the effects of particle size distributions (PSD) and soil densities on ground deformation patterns are interpreted based on the stress-strain responses and contact network characteristics of DEM RVEs. These multiscale insights enrich the knowledge of tunneling-induced ground responses and the same approach can be applied to other geotechnical engineering analyses in coarse-grained soils.

Accurately capturing the seismic response of underground structures subjected to obliquely incident seismic waves, particularly when the angle of incidence surpasses the critical value, remains a challenging task in earthquake engineering. To address this gap, this paper presents a three-dimensional (3D) nonlinear seismic analysis of subway stations embedded in a layered site, specifically in response to obliquely incident shear (SV) waves at arbitrary angles. An innovative procedure, termed the coupled dynamic stiffness matrix-finite element method (DSM-FEM), is introduced to enable seismic input by transforming responses induced by arbitrarily incoming SV waves into equivalent nodal loads. To accurately simulate wave propagation within the site, a viscous-spring artificial boundary is utilized, while a nonlinear generalized Masing model that incorporates modified damping is employed. Using the Daikai subway station as a benchmark, the research examines the effects of varying oblique incident angles on the structural response, taking into account dynamic soil-structure interaction. The results reveal that the maximum response, including peak deformation, internal forces, Mises stress, occurs when the incident angle approaches the critical value. Beyond this critical angle, the seismic response notably diminishes. Additionally, the influence of horizontal incident angles is found to be noticeable, leading to variations in deformation patterns and internal forces across different structural components. Specifically, it has been observed that the drift ratio, displacement, shear force, acceleration, and Mises stress exhibit a decreasing trend as the horizontal incident angles increase. These findings highlight the significance of considering non-vertical input ground motion in seismic analysis, and offer valuable insights for the structural design and safety evaluation of underground structures.

This study evaluates dykes stability of bauxite residue storage facility using limit equilibrium (LEM) and finite element methods (FEM), considering diverse construction phases. In LEM, steady state seepage is simulated using piezometric line while factor of safety (FOS) is determined by Morgenstern-Price method using SLOPE/W. In FEM, actual loading rates and time dependent seepage is modelled by coupled stress-pore water pressure analysis in SIGMA/W and dyke stability is assessed by stress analysis in SLOPE/W, referencing SIGMA/W analysis as a baseline model. Both the analysis incorporated suction and volumetric water content functions to determine FOS. FEM predicted pore pressures are validated against in-situ piezometer data. The results highlight that coupled hydro-mechanical analysis offers accurate stability assessment by integrating stress-strain behaviour, pore pressure changes, seepage paths, and dyke displacements with time. It is found that inclusion of unsaturated parameters in Mohr-Coulomb model improved the reliability in FOS predictions.

The current investigation examines the fluctuating behaviour of stiff pavement built on a two-parameter base and is influenced by aircraft loading impacts. This investigation is driven by the necessity for an accurate evaluation of pavement behaviour under elevated stress scenarios caused by aircraft, which can guide pavement design and upkeep. A stochastic numerical model, the vehicle-pavement interaction model (VPI), was created using a comprehensive 3D dynamic model of an aircraft vehicle and stationary runway roughness profiles. The rigid pavement is simulated using a computationally efficient 1D finite element mathematical model incorporating six DOF. The Pasternak model represents the soil medium, incorporating shear interaction between the spring elements. The pavement's irregularities are considered and replicated using a power spectral density (PSD) function. This assembled model was used to investigate the dynamical reaction of concrete pavement vibrations caused by the passing of an aircraft vehicle using MATLAB code. The dynamic governing differential equations of the aircraft's motion are developed and coupled with the pavement system equations. The coupled system is then solved in the time domain using the direct computational integration approach with the Newmark-Beta integration scheme, explicitly utilizing the linear average acceleration method. This approach is employed to resolve the equations that govern and assess the performance of the connected system. The current findings are being compared to existing analytical outcomes to verify the precision of the current coding. The research examined the impact of various pavement and aircraft vehicle behaviors and factors on the dynamic response of pavement, including the speed, main and auxiliary suspension components, mass and the load position of the aircraft, also the damping, random roughness, thickness, span length and elastic constant of the pavement, even, the modulus of subgrade of the foundation, the rigidity modulus of the shear layer. The findings demonstrate notable influences of aircraft speed and pavement surface roughness on various response parameters. Specifically, the results reveal that a higher subgrade modulus leads to decreased deflection, rotation, and bending moments. Conversely, longer span lengths tend to elevate response parameters while simultaneously reducing shear force. In conclusion, the results highlight the significance of critical factors, including velocity and subgrade modulus, in forecasting the performance of pavement subjected to aircraft loads. The present research is confined to the investigation of the dynamic's performance of the VPI simulation of airfield rigid pavement. The findings from this study can be expanded on by paving engineers to improve the structural effectiveness and reliability of the pavement, serving as a basis for subsequent fatigue analysis in response to diverse dynamic loads such as earthquake, temperature and vehicle load.

This paper presents a method for analyzing slope stability in anisotropic and heterogeneous clay using a strength reduction finite element method (SRFEM) integrated with the level set method (LSM). Anisotropy refers to the inherent anisotropy in the clay's strength, while heterogeneity describes the spatial variability in strength parameters. The static LSM uses a zero level set function to model heterogeneous clay slopes. The method is validated through undrained slope stability analyses on different types of anisotropic clay and heterogeneous fields, showing its effectiveness in modeling anisotropic shear strength and capturing the characteristics of heterogeneous regions. The results indicate that the proposed method accurately predicts factors of safety and slip surfaces across various soil conditions, accounting for both anisotropic and heterogeneous characteristics.

This study presents experimental results from scale model tests on laterally loaded bridge pile foundations in soils subjected to seasonal freezing. A refined finite-element model (FEM) was established and calibrated based on data obtained from the experiments. Furthermore, the model was utilized to investigate the impact of soil scouring depth on the lateral behavior of bridge pile foundations embedded in seasonally frozen soils. The findings indicate that soil freezing significantly enhances the lateral bearing capacity of the pile-soil interaction (PSI) system while reducing lateral deflection of the pile foundation. However, soil freezing results in increased damage to the pile foundation and upward movement of the plastic zone toward the ground surface. Under unfrozen conditions, significant plastic deformations occur on the ground surface and even inside the piles due to the extrusion effect. Additionally, increasing soil scouring depth significantly reduces the lateral bearing capacity of the PSI system while also increasing lateral deflection of the pile foundation for a given load level. Notably, when the scouring depth exceeds 2 m in unfrozen soils, the entire pile experiences obvious deformation and inclination, exhibiting a short-pile behavior that negatively affects the lateral stability of the pile under lateral loads.

In the cotton fields in Xinjiang, residual film is present in the soil for a long period of time, leading to a decrease in the tensile strength of the residual film and increasing the difficulty of recycling. Existing technologies for residual film recovery focus on mechanical properties and ignore the dragging and tearing of residual film by cotton stubble. The effect of cotton straw-root stubble on residual film recovery can only be better determined by appropriate machine operating parameters, which are essential to improving residual film recovery. Through analyses of the pickup device, key parameters were identified, and a model was built by combining the FEM and SPH algorithms to simulate the interaction of nail teeth, residual film, soil and root stubble. The simulation revealed the force change law of residual film in root stubble-containing soil and the influence of root stubble. By simulating the changes in the characteristics of the residual film during the process, the optimum operating parameters for the nail teeth were determined: a forward speed of 1849.57 mm/s, a rotational speed of 5.5 r/s and a soil penetration angle of 30 degrees. Under these optimized conditions, the maximum shear strain, pickup height (maximum deformation) and average peak stress of the residual film were 1293, 363.81 mm and 3.42 MPa, respectively. Subsequently, field trials were conducted to verify the change in the impact of the nail teeth at the optimized speed on the recovery of residual film in plots containing root stubble. The results demonstrated that when the root stubble height was 5-8 cm, the residual film averaged a recovery rate of 89.59%, with a dragging rate of only 4.10% at crossings. In contrast, 8-14 cm stubble plots showed an 82.86% average recovery and an 11.91% dragging rate. In plots with a root stubble height of 5-8 cm, compared with plots with a root stubble height of 8-14 cm, the recovery rate increased by 6.73%, and the dragging rate of residual film on root stubble decreased by 7.81%. The percentage of entangled residual film out of the total unrecovered film was 30.10% lower in the 5-8 cm stubble plots than in the 8-14 cm stubble plots. It was confirmed that the effect of cotton root stubble on residual film recovery could be reduced under appropriate machine operating parameters. This provides strong support and a theoretical and practical basis for future research on the correlation between root stubble and residual film and how to improve the residual film recovery rate.

Most natural soils exhibit a certain degree of soil structure which, in general, leads to increased strength and stiffness properties. However, the mechanical characterization of these soils based on conventional laboratory testing proves difficult in many cases due to sample disturbance. The present work aims to characterize the microstructure of a postglacial, normally consolidated, fine-grained deposit in Seekirchen, Austria, adopting in situ testing, laboratory testing on high-quality samples, and numerical analysis. The latter involves recalculating in situ piezocone penetration tests (CPTu) using an advanced constitutive model for structured soil. In contrast to existing in situ interpretation methods, the results of the numerical study, the mineralogical and hydrochemical testing, as well as the oedometer and bender element testing on undisturbed and reconstituted samples suggest that the soil is characterized by a significant amount of structure. It is demonstrated that the difference in shear wave velocity measured in situ and through bender element testing on reconstituted samples can be used as an indicator for soil structure. Ignoring the effects of structure may lead to inaccurate parameter determination for advanced constitutive models, which are subsequently employed to solve complex boundary value problems in geotechnical practice. As a consequence, the prediction of expected displacement may not be reliable.

The study presents a comprehensive study on the assessment of the bearing capacity of closely spaced strip footings on c-& oslash; soil, considering spatial variability in soil properties. A linear elastic model is employed for footings and elastic-perfect plastic soil behaviour via the Mohr-Coulomb yield criterion. Soil properties obtained from extensive field investigations of Taranto Blue Clay (TBC) in Italy are modelled as stationary random fields (RFs) generated using the Fourier series method. The cohesion and friction angle RFs are integrated with the Z-soil FE code. The final results are obtained according to the random finite element method (RFEM). The study investigates the influence of spacing distances between footings and spatial correlation lengths of soil parameters on the bearing capacity. Results show how spacing distance affects bearing capacity. Moreover, it indicates that neighbouring footing bearing capacity is strongly correlated with investigated parameters. In the case of small spatial correlation lengths, the patterns were obtained as non-symmetrical, transitioning to more symmetrical patterns at larger lengths. The manuscript concludes by presenting reliability-based design considerations for the ultimate bearing capacity, considering the horizontal spatial scale of fluctuation (SOF). The findings emphasize the importance of evaluating allowable design bearing capacity for proximity structures using RFEM and provide valuable insights into the interplay between spacing distances and spatial variability in soil properties. To this end, the study underscores the critical interplay between spacing distance, spatial correlation lengths, and random soil properties in assessing neighbouring footing-bearing capacities.