Study on the tunnel shape and soil-lining interaction influencing the lining behavior under seismic loading

The response of tunnels subjected to seismic loading is a complex mechanism and depends not only on the seismic nature but also on tunnel structure and surrounding soil properties.The individual behavior of circular,rectangular,and sub-rectangular tunnels subjected to seismic loadings has already been studied in the literature.In the present research,two case scenarios of circular,rectangular tunnels and four sub-rectangular shaped tunnels,with similar cross-section areas,were adopted to perform a comprehensive numerical investigation.The purpose of the study was to determine the mechanical behavior of tunnels of different shapes,depending upon seismic conditions.Analyses were performed by considering the influence of soil-lining interaction,soil parameters,and lining thickness,as well as lining rigidity.Computations were performed for no-slip and full-slip conditions.The results indicate that the tunnel shape design is of great importance when regarding the mechanical behavior of the surrounding soil.This concerns no-slip as well as full-slip soil-lining interaction,especially when the lining is subjected to seismic loading.Moreover,it is shown that changes in incremental bending moments for circular,rectangular and sub-rectangular tunnels that depend upon the soil-lining interaction conditions differ significantly.

Seismic wave input method for three-dimensional soil-structure dynamic interaction analysis based on the substructure of artificial boundaries

The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.

Three-dimensional seismic isolation bearing and its application in long span hangars

Based on the seismic response characteristics of space frame structures,a new type of seismic isolation bearing defined as a three-dimensional seismic isolation bearing（3DSIB） is developed in this paper.The bearing offers excellent properties such as multi-dimensional seismic isolation,reasonable rotation capability,good ability to resist lifting load,uncoupled stiffness in horizontal and vertical directions,etc.In the 3DSIB,the horizontal dimension is designed by combining the Teflon sliding device and helical spring,while the vertical dimension is developed by introducing disk springs or helical springs.The mathematical model of the 3DSIB was established and its performance with the critical parameters was tested on a shaking table.Furthermore,the 3DSIB was applied in a 120 m span hangar structure and simulated using SAP2000 software to evaluate its performance in practical structures.The performance of the structures with and without 3DSIB was compared.It is shown that the hangar structure with 3D bearings achieves a better performance.The axial force and acceleration response of the structures with 3DSIB are effectively reduced,while the displacement response of the bearing is within the predetermined range.

Limit analysis for the seismic stability of three-dimensional rock slopes using the generalized Hoek-Brown criterion

The parameters that influence slope stability and their criteria of failure are fairly understood but over-conservative design approaches are often preferred,which can result in excessive overburden removal that may jeopardize profitability in the context of open pit mining.Numerical methods such as finite element and discrete element modelling are instrumental to identify specific zones of stability,but they remain approximate and do not pinpoint the critical factors that influence stability without extensive parametric studies.A large number of degrees of freedom and input parameters may make the outcome of numerical modelling insufficient compared to analytical solutions.Existing analytical approaches have not tackled the stability of slopes using non-linear plasticity criteria and threedimensional failure mechanisms.This paper bridges this gap by using the yield design theory and the Hoek-Brown criterion.Moreover,the proposed model includes the effect of seismic forces,which are not always taken into account in slope stability analyses.The results are presented in the form of rigorous mathematical expressions and stability charts involving the loading conditions and the rock mass properties emanating from the plasticity criterion.

Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

Probabilistic seismic stability of three-dimensional slopes by pseudo-dynamic approach

Probabilistic analysis is a rational approach for engineering design because it provides more insight than traditional deterministic analysis. Probabilistic evaluation on seismic stability of three dimensional (3D) slopes is studied in this paper. The slope safety factor is computed by combining the kinematic approach of limit analysis using a three-dimensional rotational failure mechanism with the pseudo-dynamic approach. The variability of input parameters, including six pseudo-dynamic parameters and two soil shear strength parameters, are taken into account by means of Monte-Carlo Simulations (MCS) method. The influences of pseudo-dynamic input variables on the computed failure probabilities are investigated and discussed. It is shown that the obtained failure probabilities increase with the pseudo-dynamic input variables and the pseudo-dynamic approach gives more conservative failure probability estimates compared with the pseudo-static approach.

Research and application on three-dimensional seismic and vibration isolation for building

This paper presents the study of a three-dimensional(3D) isolation system.Firstly,the authors investigated the effects of an innovative 3D isolator,which was composed of a connecting plate,a rubber pad for vibration isolation in the vertical direction and a horizontal rubber bearing for seismic isolation in both horizontal directions.Secondly,the authors designed such a vibration isolation system and installed it underneath two specific residential buildings which were built directly over an existing subway communication hub platform in Beijing.These buildings required good performance vibration and seismic isolation system to reduce the impact from the running of nearby subway trains.Finally,in situ tests were conducted for both the isolated and the non-isolated buildings for the purpose of comparison.The test results showed that the maximum acceleration response level of the isolated superstructure is reduced by 10% as compared to that of the platform.The maximum attenuation of vibration reaches up to 25 dB.The 3D system explored in this paper is very effective in control and suppression of building vibration induced by earthquakes or running of trains.

THREE-DIMENSIONAL ELLIPTIC CRACK UNDER IMPACT LOADING

The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.

Three-dimensional seismic velocity tomography of the upper crust in Tengchong volcanic area, Yunnan Province

Based on data collected by deep seismic sounding carried out in 1999, a three-dimensional P wave velocity structure is determined with tomographic inversion. The tomographic result shows that there is a P wave low velocity zone (LVZ) in the upper crust beneath the Tengchong volcanic area. The LVZ is in the depth of 7~8 km and may be a smgma chamber or a partial melting body. The result also shows that the LVZ is in the northeastern side of the Rehai hydrothermal field, which is located in another LVZ near the surface. The shallow LVZ may represent a well-developed fracture zone. The strong hydrothermal activity in Rehai area can attribute to the existence of fractures between two LVZs. These fractures are the channels for going upwards of the deep hot fluid.

Numerical study on wave loads and motions of two ships advancing in waves by using three-dimensional translating-pulsating source

A frequency domain analysis method based on the three-dimensional translating-pulsating （3DTP） source Green function is developed to investigate wave loads and free motions of two ships advancing on parallel course in waves. Two experiments are carried out respectively to mea- sure the wave loads and the free motions for a pair of side-by- side arranged ship models advancing with an identical speed in head regular waves. For comparison, each model is also tested alone. Predictions obtained by the present solution are found in favorable agreement with the model tests and are more accurate than the traditional method based on the three dimensional pulsating （3DP） source Green function. Numer- ical resonances and peak shift can be found in the 3DP pre- dictions, which result from the wave energy trapped in the gap between two ships and the extremely inhomogeneous wave load distribution on each hull. However, they can be eliminated by 3DTP, in which the speed affects the free sur- face and most of the wave energy can be escaped from the gap. Both the experiment and the present prediction show that hydrodynamic interaction effects on wave loads and free motions are significant. The present solver may serve as a validated tool to predict wave loads and motions of two ves- sels under replenishment at sea, and may help to evaluate the hydrodynamic interaction effects on the ships safety in replenishment operation.

Investigation of control effects of end-wall selfadaptive jet on three-dimensional corner separation of a highly loaded compressor cascade

To overcome the limitations posed by three-dimensional corner separation,this paper proposes a novel flow control technology known as passive End-Wall(EW)self-adaptive jet.Two single EW slotted schemes(EWS1 and EWS2),alongside a combined(COM)scheme featuring double EW slots,were investigated.The results reveal that the EW slot,driven by pressure differentials between the pressure and suction sides,can generate an adaptive jet with escalating velocity as the operational load increases.This high-speed jet effectively re-excites the local low-energy fluid,thereby mitigating the corner separation.Notably,the EWS1 slot,positioned near the blade leading edge,exhibits relatively low jet velocities at negative incidence angles,causing jet separation and exacerbating the corner separation.Besides,the EWS2 slot is close to the blade trailing edge,resulting in massive low-energy fluid accumulating and separating before the slot outlet at positive incidence angles.In contrast,the COM scheme emerges as the most effective solution for comprehensive corner separation control.It can significantly reduce the total pressure loss and improve the static pressure coefficient for the ORI blade at 0°-4° incidence angles,while causing minimal negative impact on the aerodynamic performance at negative incidence angles.Therefore,the corner stall is delayed,and the available incidence angle range is broadened from -10°--2°to -10°-4°.This holds substantial promise for advancing the aerodynamic performance,operational stability,and load capacity of future highly loaded compressors.

Slope stability analysis under seismic load by vector sum analysis method

The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method （DFEM）. Combining with the DFEM, the vector sum analysis method （VSAM） is employed in seismic stability analysis of a slope in this paper. Different from other conventional methods, the VSAM is proposed based on the vector characteristic of force and current stress state of the slope. The dynamic stress state of the slope at any moment under seismic load can he obtained by the DFEM, thus the factor of safety of the slope at any moment during earthquake can be easily obtained with the VSAM in consideration of the DFEM. Then, the global stability of the slope can be estimated on the basis of time-history curve of factor of safety and reliability theory. The VSAM is applied to a homogeneous slope under seismic load. The factor of safety of the slope is 1.30 under gravity only and the dynamic factor of safety under seismic load is 1.21. The calculating results show that the dynamic characteristics and stability state of the slope with input ground motion can be actually analyzed. It is believed that the VSAM is a feasible and practical approach to estimate the dynamic stability of slopes under seismic load.

Earth slope reliability analysis under seismic loadings using neural network

A new method was proposed to cope with the earth slope reliability problem under seismic loadings. The algorithm integrates the concepts of artificial neural network, the first order second moment reliability method and the deterministic stability analysis method of earth slope. The performance function and its derivatives in slope stability analysis under seismic loadings were approximated by a trained multi-layer feed-forward neural network with differentiable transfer functions. The statistical moments calculated from the performance function values and the corresponding gradients using neural network were then used in the first order second moment method for the calculation of the reliability index in slope safety analysis. Two earth slope examples were presented for illustrating the applicability of the proposed approach. The new method is effective in slope reliability analysis. And it has potential application to other reliability problems of complicated engineering structure with a considerably large number of random variables.

Dynamic behavior of single-layer latticed cylindrical shells subjected to seismic loading

The single-layer latticed cylindrical shell is one of the most widely adopted space-fl'amed structures.In this paper,free vibration properties and dynamic response to horizontal and vertical seismic waves of single-layer latticed cylindrical shells are analyzed by the finite element method using ANSYS software.In the numerical study,where hundreds of cases were analyzed,the parameters considered included rise-span ratio,length-span ratio,surface load and member section size.Moreover,to better define the actual behavior of single-layer latticed shells,the study is focused on the dynamic stress response to both axial forces and bending moments.Based on the numerical results,the effects of the parameters considered on the stresses are discussed and a modified seismic force coefficient method is suggested.In addition,some advice based on these research results is presented to help in the future design of such structures.

Topology optimization of shear wall structures under seismic loading

A topology optimization formulation is developed to find the stiffest structure with desirable material distribution subjected to seismic loads. Finite element models of the structures are generated and the optimality criteria method is modified using a simple penalty approach and introducing fictitious strain energy to simultaneously consider both material volume and displacement constraints. Different types of shear walls with/without opening are investigated. Additionally, the effects of shear wall-frame interaction for single and coupled shear walls are studied. Gravity and seismic loads are applied to the shear walls so that the definitions provide a practical approach for locating the critical parts of these structures. The results suggest new viewpoints for architectural and structural engineering for placement of openings.

A Numerical Study of Wave Propagation and Cracking Processes in Rock-Like Material under Seismic Loading Based on the Bonded-Particle Model Approach

An earthquake is usually followed by a considerable number of aftershocks that play a significant role in earthquake-induced landslides,During the aftershock,the cracking process in rocks becomes more complex because of the formation of faults.In order to investigate the effects of seismic loading on the cracking processes in a specimen containing a single flaw,a numerical approach based on the bonded-particle model(BPM)was adopted to study the seismic loading applied in two orthogonal directions.The results reveal that no transmission and reflection phenomena were observable in the small specimens(76 mm×152 mm)because they were considerably smaller than the wavelength of the P-wave.Furthermore,under seismic loading,the induced crack was solely tensile in nature.Repeated axial seismic loading did not induce crack propagation after the first axial seismic loading.Cracks began to propagate only when the seismic loading direction was changed from axial to lateral,and then back to axial,ultimately resulting in the failure of the specimen.

Three-dimensional P velocity structure in Beijing area

A detail three-dimensional P wave velocity structure of Beijing, Tianjin and Tangshan area (BTT area) was deter-mined by inverting local earthquake data. In total 16 048 P wave first arrival times from 16048 shallow and mid-depth crustal earthquakes, which occurred in and around the BTT area from 1992 to 1999 were used. The first arrival times are recorded by Northern China United Telemetry Seismic Network and Yanqing-Huailai Digital Seismic Network. Hypocentral parameters of 1 132 earthquakes with magnitude ML=1.7~6.2 and the three-dimensional P wave velocity structure were obtained simultaneously. The inversion result reveals the com-plicated lateral heterogeneity of P wave velocity structure around BTT area. The tomographic images obtained are also found to explain other seismological observations well.

Three-dimensional consolidation deformation analysis of porous layered soft soils considering asymmetric effects

Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on the long-term deformation for underground structures. A three-dimensional consolidation analysis method under the asymmetric loads is developed for porous layered soil based on Biot's classical theory. Time-displacement effects can be fully considered in this work and the analytical solutions are obtained by the state space approach in the Cartesian coordinate. The Laplace and double Fourier integral transform are applied to the state variables in order to reduce the partial differential equations into algebraic differential equations and easily obtain the state space solution. Starting from the governing equations of saturated porous soil, the basic relationship of state space variables is established between the ground surface and the arbitrary depth in the integral transform domain. Based on the continuity conditions and boundary conditions of the multi-layered pore soil model, the multi-layered pore half-space solutions are obtained by means of the transfer matrix method and the inverse integral transforms. The accuracy of proposed method is demonstrated with existing classical solutions. The results indicate that the porous homogenous soils as well as the porous non-homogenous layered soils can be considered in this proposed method. When the consolidation time factor is 0.01, the value of immediate consolidation settlement coefficient calculated by the weighted homogenous solution is 27.4% bigger than the one calculated by the non-homogeneity solution. When the consolidation time factor is 0.05, the value of excess pore water pressure for the weighted homogenous solution is 27.2% bigger than the one for the non-homogeneity solution. It is shown that the material non-homogeneity has a great influence on the long-term settlements and the dissipation process of excess pore water pressure.

Autologous nerve graft repair of different degrees of sciatic nerve defect:stress and displacement at the anastomosis in a three-dimensional finite element simulation model

In the repair of peripheral nerve injury using autologous or synthetic nerve grafting, the mag- nitude of tensile forces at the anastomosis affects its response to physiological stress and the ultimate success of the treatment. One-dimensional stretching is commonly used to measure changes in tensile stress and strain; however, the accuracy of this simple method is limited. There- fore, in the present study, we established three-dimensional finite element models of sciatic nerve defects repaired by autologous nerve grafts. Using PRO E 5.0 finite element simulation software, we calculated the maximum stress and displacement of an anastomosis under a 5 N load in 10-, 20-, 30-, 40-mm long autologous nerve grafts. We found that maximum displacement increased with graft length, consistent with specimen force. These findings indicate that three-dimensional finite element simulation is a feasible method for analyzing stress and displacement at the anas- tomosis after autologous nerve grafting.

A computational platform for considering the effects of aerodynamic and seismic load combination for utility scale horizontal axis wind turbines

The wide deployment of wind turbines in locations with high seismic hazard has led engineers to take into account a more comprehensive seismic design of such structures. Turbine specific guidelines usually use simplified methods and consider many assumptions to combine seismic demand with the other operational loads effecting the design of these structures. As the turbines increase in size and capacity, the interaction between seismic loads and aerodynamic loads becomes even more important. In response to the need for a computational tool that can perform coupled simulations of wind and seismic loads, a seismic module is developed for the FAST code and described in this research. This platform allows engineers working in this industry to directly consider interaction between seismic and other environmental loads for turbines. This paper details the practical application and theory of this platform and provides examples for the use of different capabilities. The platform is then used to show the suitable earthquake and operational load combination with the implicit consideration of aerodynamic damping by estimating appropriate load factors.