循环荷载作用下不均质地层中盾构隧道纵向响应研究

修建在纵向不均质地层中的地铁隧道,由于列车循环荷载的作用,会导致隧道下部的土体产生不均匀沉降,对既有隧道产生不利的影响。针对这一问题,提出考虑隧道剪切效应的地基不均匀沉降对既有隧道竖向变形影响的解析解。既有隧道简化为搁置在Winkler地基上的Timoshenko梁,通过两阶段分析法,分析下卧地层不均匀沉降引起的隧道响应。首先确定列车荷载引起的动偏应力,并运用土层的力学指标计算出静偏应力和破坏偏应力。然后运用累积应变的经验公式计算出隧道下部土体的累计沉降,将土体的沉降转化为力施加在隧道上。基于Timoshenko梁理论,建立考虑隧道剪切效应的隧道竖向变形微分方程,求解得到隧道变形的解析解,进一步可以得到隧道的弯矩、剪力、转角、错台。     

地震荷载作用下兰州地铁隧道结构动力响应

根据《兰州轨道交通1号线一期工程地震安全性评价报告》所给出的100年超越概率63%、10%和2%的场地基岩地震加速度时程,利用有限差分软件进行地下隧道硐室的地震反应分析。在模型底部施加基岩地震动,设置监测点监测衬砌结构的弯矩、轴力及剪力随时间的变化过程,得到100年超越概率63%、10%及2%工况下的隧道结构地震响应。结果表明:隧道衬砌结构最大弯矩位于拱顶处,最大轴力位于拱顶和拱底处,最大剪力位于上侧壁或下侧壁处;隧道结构内力随着超越概率的降低而增大;以超越概率63%的结构最大内力为基准值,在超越概率10%和2%时,弯矩分别增大1.2和1.7倍,轴力分别增大1.3和1.5倍,剪力分别增大1.5和2.9倍,增幅最大。这可能预示着隧道结构在强地震动作用下会发生剪切破坏。     

地震作用下预应力混凝土管桩运动响应三维数值分析

基于黏弹性人工边界,建立上部结构-桩-土的共同作用三维有限元模型,分析地震作用下预应力混凝土管桩的运动响应特性。分别针对预应力混凝土管桩的桩径、双层软硬土剪切波速比值、上覆土层厚度、上部结构荷载等影响因素进行数值计算。参数分析表明:在地震作用下,桩径的增大会导致桩身整体弯矩相应增加,特别是桩身土层分界面处增大明显;软硬土层剪切波速比及上覆土层厚度的增加,引起土层分界面处桩身峰值弯矩增加;固定桩头条件下,桩头与桩身软硬土层分界面处均会产生较大的运动弯矩;上部结构的惯性荷载对固定桩头的内力有着较大影响,对桩身深处段弯矩影响较小。本文研究结论可为预应力混凝土管桩抗震设计提供有益的理论参考。     

基于动力BNWF法的冻土-桩相互作用模型研究

基于水平循环荷载作用下不同负温冻土环境中单桩动力特性模型试验结果,在已有分析桩-土-结构相互作用的动力BNWF模型的基础上,提出改进的冻土-桩基动力相互作用非线性反应分析模型。在该模型中,利用改进的双向无拉力多段屈服弹簧考虑桩侧冻土的水平非线性力学特性,同时兼顾桩侧与冻土间的竖向非线性摩擦效应、桩尖土的挤压与分离作用以及远场土体阻尼对桩基动力特性的影响。其中桩侧水平多段屈服弹簧参数根据冻土非线性p-y关系获得,该关系曲线以三次函数曲线段及常值函数段共同模拟,并由室内冻土压缩试验结果确定。最后基于改进的动力BNWF模型,提取动位移荷载作用下该桩顶力-位移滞回曲线及桩身不同埋深处的弯矩动响应数值分析结果,并与相应的模型试验结果对比,二者具有较好的拟合度,表明本文所提出的改进模型在分析冻土-桩动力相互作用时有较好的适用性。     

Kinematic internal forces in deep foundations with inclined piles

This paper presents a parametric study that looks into the influence of pile rake angle on the kinematic internal forces of deep foundations with inclined piles. Envelopes of maximum kinematic bending moments, shear forces and axial loads are presented along single inclined piles and 2 × 2 symmetrical square pile groups with inclined elements subjected to an earthquake generated by vertically incident shear waves. Inclination angles from 0° to 30° are considered, and three different pile–soil stiffness ratios are studied. These results are obtained through a frequency–domain analysis using a boundary element–finite element code in which the soil is modelled by the boundary element method as a homogeneous, viscoelastic, unbounded region, and the piles are modelled by finite elements as Euler–Bernoulli beams. The rotational kinematic response of the pile foundations is shown to be a key factor on the evolution of the kinematic internal forces along the foundations. Copyright © 2015 John Wiley & Sons, Ltd.     

Soil-structure interaction effects on the seismic response of tall chimneys

The effects of soil-structure interaction on the seismic response of tall (>100m) steel and reinforced concrete chimneys are described. Detailed models of a 130m high steel chimney and a 150m high reinforced concrete chimney are used as structural models. The foundations are represented as rigid blocks resting on a uniform viscoelastic soil model. Perfect bonding between the foundation and the soil is assumed. Parametric studies of the interaction effects on the magnitude and distribution of bending moments and shear forces include four soil rigidities and two seismic excitations characterized by very different frequency contents. The results obtained indicate strong interaction effects for intermediate and soft soils (V_s500 m/sec). The extent of the interaction effects are highly dependent on the characteristics of the seismic excitation.     

Static equivalent method for the kinematic interaction analysis of single piles

This paper presents a static equivalent approach to estimate the maximum kinematic interaction effects on piles subjected to lateral seismic excitation. Closed-form expressions are reported for the evaluation of the maximum free-field soil movements and for the computation of maximum pile shear force and bending moments. Firstly, modal analysis, combined with a suitable damped response spectrum, is used to evaluate the maximum free-field response. Secondly, the pile is schematised as a Winkler's beam subjected to equivalent static forces defined according to soil vibration modal shapes and amplitude. The method may be applied by using response spectra suggested by National Standards or those obtained with accelerograms. The procedure proposed may be conveniently implemented in simple spreadsheets or in commercial finite element programs and easily used by practicing engineers. Method accuracy is demonstrated by comparing the results with those obtained with a more rigorous model. Good results may be achieved by considering only the first soil vibration mode making the procedure straightforward for practical design purposes.     

Experimental p-y loops for estimating seismic soil-pile interaction

Seismic soil-pile interaction is evaluated in this study based on back-calculated p-y loops constructed from sampled data of pile bending moments. Fundamental properties of p-y loops are implemented to derive distributed springs and dashpots, thereby quantifying soil-pile interaction in the realm of a Beam on Dynamic Winkler Foundation modeling. The procedure is validated by means of well-documented centrifuge tests of a single pile supported structure founded on a two-layer soil profile that comprises of soft clay overlying dense sand. Two shaking levels of a real earthquake motion applied at the base of the soil profile were examined and the generated seismic p-y loops were compared to cyclic p-y curves commonly used in pile design practice. The results demonstrate the strong influence of intensity of the input motion on seismic p-y loops while cyclic p-y curves established for soft clays tend to overestimate soil stiffness under strong excitation. Typical sets of recorded and computed structural response are presented, denoting the ability of the BDWF model related to p-y loops in reproducing adequately fundamental aspects of seismic soil-pile interaction.     

Centrifuge modeling of interaction between reverse faulting and tunnel

In this study, a series of centrifuge tests, modeling reverse fault rupture with 60° dip angle, were conducted in a dry sandy soil with a tunnel embedded in the soil layer. The test results showed that the tunnel and soil responses depended on the tunnel position, soil relative density and tunnel rigidity. Tunnels appeared be able to deviate the fault rupture path, while this deviation may be associated with significant rotation and displacement of the tunnel. However, a deeper tunnel was able to diffuse the shear deformation to a wider zone with an unsmooth surface displacement which may cause severe damage to the surface structures. Finally, the tunnel rotation, the location of the fault outcropping, the vertical displacement of the ground surface, the effect of tunnel rigidity on fault rupture path and surface displacement and the effect of soil relative density on fault–tunnel interaction were reported and discussed in this study.     

超长混凝土框架梁弹性内力行波效应研究

鉴于超长混凝土框架结构行波效应规律尚缺乏系统研究,借助建筑结构通用有限元软件MI-DAS/GEN,采用一致与非一致两种激励方式对165～1100m/s地震波速作用下长度为33～264m的3层钢筋混凝土框架结构进行弹性时程分析。通过对比一致与非一致激励下的梁弯矩、剪力和轴力,详细讨论了其随波速、长度以及不同部位的变化规律。研究表明:波速越慢,模型越长,行波效应越明显。框架梁弯矩的增大作用主要集中在264m模型第1层两端,减小作用表现在各模型第1层中部和第2、3层的梁上;剪力与弯矩有相似的变化规律;梁中存在拉力与压力,中部各梁的拉力和压力较两端的梁要大。     

Vulnerability of floating tunnel shafts for increasing earthquake loading

Fragility curves constitute the cornerstone in seismic risk evaluations and performance-based earthquake engineering. They describe the probability of a structure to experience a certain damage level for a given earthquake intensity measure, providing a relationship between seismic hazard and vulnerability. In this paper a numerical approach is applied to derive fragility curves for tunnel shafts built in clays, a component that is found in several critical infrastructure such as urban metro networks, airport facilities or water and waste water projects. The seismic response of a representative tunnel shaft is assessed using tridimensional finite difference non-linear analyses carried out with the program FLAC^3D, under increasing levels of seismic intensity. A hysteretic model is used to simulate the soil non-linear behavior during the seismic event. The effect of soil conditions and ground motion characteristics on the soil-structure system response is accounted for in the analyses. The damage is defined based on the exceedance of the concrete wall shaft capacity due to the developed seismic forces. The fragility curves are estimated in terms of peak ground acceleration at a rock or stiff soil outcrop, based on the evolution of damage with increasing earthquake intensity. The proposed fragility models allows the characterization of the seismic risk of a representative tunnel shaft typology and soil conditions considering the associated uncertainties, and partially fill the gap of data required in performing a risk analysis assessment of tunnels shafts.     

Dynamic response of a tunnel buried in a saturated poroelastic soil layer to a moving point load

The dynamic response of a tunnel buried in a two-dimensional poroelastic soil layer subjected to a moving point load was investigated theoretically. The tunnel was simplified as an infinite long Euler–Bernoulli beam, which was placed parallel to the traction-free ground surface. The saturated layer was governed by Biot’s theory. Combined with the specified boundary conditions along the beam and saturated poroelastic layer, the coupled equations of the system were solved analytically in the frequency–wavenumber domain based on Fourier transform. The time domain responses were obtained by the fast inverse Fourier transform. The critical velocity of the considered structure was determined from the dispersion curves. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium subjected to the underground moving load were investigated. It is concluded that, for coarse materials or fine materials subjected to the high-velocity loading, models ignoring the coupling effects between the pore fluid and the soil skeleton may cause errors. The shear modulus and the permeability coefficients of the saturated soil as well as the load moving velocity had significant influence on the displacement and pore pressure responses.     

Effect of a forced harmonic vibration pile to its adjacent pile in layered elastic soil with double-shear model

A new model named double-shear model based on Pasternak foundation and Timoshenko beam theory is developed to evaluate the effect of a forced harmonic vibration pile to its adjacent pile in multilayered soil medium. The double-shear model takes into account the shear deformation and the rotational inertia of piles as well as the shear deformation of soil. The piles are simulated as Timoshenko beams, which are embedded in a layered Pasternak foundation. The differential equation of transverse vibration for a pile is solved by the initial parameter method. The dynamic interaction factors for the layered soil medium are obtained by the transfer matrix method. The formulation and the implementation have been verified by means of several examples. The individual shear effects of soil and piles on the interaction factors are evaluated through a parametric study. Compared to Winkler model with Euler beam, the present model gives much better results for the dynamic interaction of piles embedded in stiff soil with small slenderness ratios. Finally, the effect of a forced long pile to a short pile embedded in multilayered soil medium is studied in detail.     

Kinematic response of single piles for different boundary conditions: Analytical solutions and normalization schemes

Kinematic pile–soil interaction is investigated analytically through a Beam-on-Dynamic-Winkler-Foundation model. A cylindrical vertical pile in a homogeneous stratum, excited by vertically-propagating harmonic shear waves, is examined in the realm of linear viscoelastic material behaviour. New closed-form solutions for bending, displacements and rotations atop the pile, are derived for different boundary conditions at the head (free, fixed) and tip (free, hinged, fixed). Contrary to classical elastodynamic theory where pile response is governed by six dimensionless ratios, in the realm of the proposed Winkler analysis three dimensionless parameters suffice for describing pile–soil interaction: (1) a mechanical slenderness accounting for geometry and pile–soil stiffness contrast, (2) a dimensionless frequency (which is different from the classical elastodynamic parameter a₀=ω d/V_s), and (3) soil material damping. With reference to kinematic pile bending, insight into the physics of the problem is gained through a rigorous superposition scheme involving an infinitely-long pile excited kinematically, and a pile of finite length excited by a concentrated force and a moment at the tip. It is shown that for long piles kinematic response is governed by a single dimensionless frequency parameter, leading to a unique master curve pertaining to all pile lengths and pile–soil stiffness ratios.     

Evaluation of seismic performance of pile‐supported models in liquefiable soils

The seismic performance of four pile‐supported models is studied for two conditions: (i) transient to full liquefaction condition, i.e. the phase when excess pore pressure gradually increases during the shaking; (ii) full liquefaction condition, i.e. defined as the state where the seismically induced excess pore pressure equalises to the overburden stress. The paper describes two complementary analyses consisting of an experimental investigation, carried out at normal gravity on a shaking table, and a simplified numerical analysis, whereby the soil–structure interaction (SSI) is modelled through non‐linear Winkler springs (commonly known as p–y curves). The effects of liquefaction on the SSI are taken into account by reducing strength and stiffness of the non‐liquefied p–y curves by a factor widely known as p‐multiplier and by using a new set of p–y curves. The seismic performance of each of the four models is evaluated by considering two different criteria: (i) strength criterion expressed in terms of bending moment envelopes along the piles; (ii) damage criterion expressed in terms of maximum global displacement. Comparison between experimental results and numerical predictions shows that the proposed p–y curves have the advantage of better predicting the redistribution of bending moments at deeper elevations as the soil liquefies. Furthermore, the proposed method predicts with reasonable accuracy the displacement demand exhibited by the models at the full liquefaction condition. However, disparities between computed and experimental maximum bending moments (in both transient and full liquefaction conditions) and displacement demands (during transient to liquefaction condition) highlight the need for further studies. Copyright © 2016 The Authors Earthquake Engineering & Structural Dynamics Published by John Wiley & Sons Ltd.     

Winkler model for lateral response of rigid caisson foundations in linear soil

A generalized spring multi-Winkler model is developed for the static and dynamic response of rigid caisson foundations of circular, square, or rectangular plan, embedded in a homogeneous elastic. The model, referred to as a four-spring Winkler model, uses four types of springs to model the interaction between soil and caisson: lateral translational springs distributed along the length of the caisson relating horizontal displacement at a particular depth to lateral soil resistance (resultant of normal and shear tractions on the caisson periphery); similarly distributed rotational springs relating rotation of the caisson to the moment increment developed by the vertical shear tractions on the caisson periphery; and concentrated translational and rotational springs relating, respectively, resultant horizontal shear force with displacement, and overturning moment with rotation, at the base of the caisson. For the dynamic problem each spring is accompanied by an associated dashpot in parallel. Utilising elastodynamic theoretical available in the literature results for rigid embedded foundations, closed-form expressions are derived for the various springs and dashpots of caissons with rectangular and circular plan shape. The response of a caisson to lateral static and dynamic loading at its top, and to kinematically-induced loading arising from vertical seismic shear wave propagation, is then studied parametrically. Comparisons with results from 3D finite element analysis and other available theoretical methods demonstrate the reliability of the model, the need for which arises from its easy extension to multi-layered and nonlinear inelastic soil. Such an extension is presented in the companion papers by the authors [Gerolymos N, Gazetas G. Development of Winkler model for lateral static and dynamic response of caisson foundations with soil and interface nonlinearities. Soil Dyn Earthq Eng. Submitted companion paper; Gerolymos N, Gazetas G. Static and dynamic response of massive caisson foundations with soil and interface nonlinearities—validation and results. Soil Dyn Earthq Eng. Submitted companion paper.].     

开挖过程及围护结构对地铁车站抗震性能的影响

为了研究开挖工序以及围护结构作为永久结构在明挖车站的抗震设计校核中的作用,基于以灌注桩作为支挡结构的郑州市某明挖结构地铁站,采用Plaxis2D建立二层单柱双跨矩形框架式有限元模型,采取硬化土本构,建立考虑分步开挖、和不考虑开挖的模型,进行有无围护墙及不同墙厚抗震验算.研究结果表明:侧墙作为主要抗侧力构件,在地震荷载作...     

Evaluation of substructuring method for seismic soil-structure interaction analysis of bridges

This paper evaluates the commonly used substructuring method for analysis of bridge systems where the bridge is divided into two sub-systems: the bridge superstructure and the substructure including the pile foundations, abutments, and soil. Modeling of the soil-structure interaction (SSI) in the system is simplified by replacing the pile foundations, abutments, and soil with sets of independent equivalent linear springs and dashpots at the base of the superstructure. The main objective of the paper is to examine how well the substructuring method simulates the seismic response of a bridge system. The baseline data required for the evaluation process is derived from analyzing a fully-coupled continuum bridge model, already validated for the instrumented two-span Meloland Road Overpass. The same bridge system is also simulated using the substructuring method. The results from both approaches are compared, and it is shown that the differences between them can be significant. The substructuring method consistently overestimates the pier base shear forces and bending moments and the pier top deflections. Moreover, the spectral response of the bridge structure is mispredicted. The analyses are repeated for a three-span bridge system subjected to several ground motions, leading to a similar observation as before. Hence, the current state of practice for simulating seismic SSI in bridges using the substructure model is shown to be too simplified to capture the major mechanisms involved in SSI.     

Nonlinear dynamic analysis of Meloland Road Overpass using three-dimensional continuum modeling approach

This paper presents a three-dimensional (3D) continuum nonlinear analysis of the Meloland Road Overpass (MRO) near El Centro, California. The modeling methodology and the computational tools are discussed in detail. The performance of the computational model is evaluated by comparing the computed responses with the responses recorded at the bridge site during the 1979 Imperial Valley and 2010 El Mayor-Cucapah earthquakes. Amongst the recorded earthquake events at the bridge site, these two events caused the strongest shaking. The comparison shows that the 3D model is potentially an effective tool for detailed analysis of a full bridge system including foundation soils, pile foundations, embankments, supporting columns, and the bridge structure itself in a unified system without relying on any ancillary models such as Winkler springs. Additional response parameters such as displacements, rockings, and bending moments are also evaluated although none of these was measured during the seismic events.     

Seismic pressures on rigid cantilever walls retaining elastic continuously non-homogeneous soil: An exact solution

The dynamic response of an elastic continuously nonhomogeneous soil layer over bedrock retained by a pair of rigid cantilever walls to a horizontal seismic motion and the associated seismic pressure acting on these walls are determined analytically–numerically. The soil non-homogeneity is described by a shear modulus increasing nonlinearly with depth. The problem is solved in the frequency domain under conditions of plane strain and its exact solution is obtained analytically. This is accomplished with the aid of Fourier series along the horizontal direction and solution of the resulting system of two ordinary differential equations with variable coefficients by the method of Frobenius in power series. Due to the complexity of the various analytical expressions, the final results are determined numerically. These results include seismic pressures, resultant horizontal forces and bending moments acting on the walls. The solution of the problem involving a single retaining wall can be obtained as a special case by assuming the distance between the two walls to be very large. Results are presented in terms of numerical values and graphs using suitable dimensionless quantities. The effect of soil non-homogeneity on the system response is assessed through comparisons for typical sets of the parameters involved.