钢套管灌注群桩施工顺序对多隧道的影响

通过理论研究和数值分析方法,分析了钢套管桩不同的施工顺序对邻近既有隧道的影响,并结合现场实测数据对分析结果进行了验证。首先,采用理论方法分析钢套管旋转下压过程中土塞柱的受力情况,建立微分方程计算土塞柱底端的应力、土塞柱受到的剪切扭矩和自身抗剪扭矩,进而分析得到了土塞柱发生剪切破坏和土塞效应的临界深度。在此基础上按是否产生闭塞效应,分两种情况研究了钢套管旋压过程中产生的土体膨胀率。然后,利用既得的土体膨胀率,通过理论分析和数值模拟方法分析了群桩不同的施工顺序对既有多条盾构隧道的影响,得出了相对较优的群桩施工顺序方案。工程现场实测数据与数值模拟结果的对比分析表明,上述计算结果与实际情况能较好吻合。该研究对类似工程具有理论指导意义和实践参考价值。     

液化流动条件下隧道结构动力分布场研究

将液化流动的土体视为流体,运用流体力学原理,基于矢量符号运算法进行液化场地的动力场分析,求解出动力场解析解。采用ABAQUS/CFD进行液化场地流体动力学有限元分析,根据流动基本控制方程,计算出动力场的数值解。分析结果表明,液化土体横向流动时,隧道结构周围的应力场既包括由表面压力引起的压力阻力,也包括由剪应力引起的摩擦阻力;隧道结构周围的液化土体流动速度很小,但在隧道结构下方存在流动速度加强区;理论计算的动力场解析解大于有限元计算的数值解,但总体的分布变化趋势基本一致;隧道结构附近的应力场变化较为密集,所受到的应力主要分布在液化土体流动的迎面位置。     

砂土中开口管桩沉桩过程的颗粒流模拟研究

基于颗粒流理论,采用PFC2D程序,模拟再现不同型号开口管桩在沉桩过程中土塞的形成演化规律、土颗粒细观结构变化以及桩周土应力场分布情况,并通过分析土体细观变化模式揭示沉桩过程中宏观力学响应的内在机制。计算结果表明,管桩直径对土塞效应影响很大,外径为30 mm的开口管桩,沉桩过程中土塞增量填充率（IFR）值较小,土塞效应明显,土塞高度小,类似闭口桩;随着管桩直径的增大,土塞效应迅速减小,大直径管桩在砂土中沉桩全部闭塞的可能性很小。细观因素（孔隙率和滑动比例）与土体宏观位移表现之间存在着明显的相互对应关系,并依此将桩周土划分3个区域。桩周土体水平应力、竖向应力和剪应力都在桩底附近形成“应力核”,不同型号管桩桩周土应力场分布相近。     

Two-dimensional Dynamics Simulation of Two-phase Debris Flow

To investigate the movement mechanism of debris flow, a two‐dimensional, two‐phase, depth‐integrated model is introduced. The model uses Mohr‐Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second‐order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.     

Numerical modelling of the driving response of thin‐walled open‐ended piles

The driving response of thin‐walled open‐ended piles is studied using numerical simulation of the wave propagation inside the soil plug and the pile. An elastic finite element analysis is carried out to identify the stress wave propagation in the vicinity of the pile toe. It is found that the shear stress wave has the highest magnitude above the bottom of the soil plug. Below the bottom of the soil plug, the vertical stress wave has the highest magnitude. Although the shear stress wave propagating in the radial direction is similar in magnitude to the vertical stress wave at the bottom of the soil plug, it decays rapidly while travelling downwards. The highest vertical stress at the bottom of the soil plug appears after the vertical stress wave interacts with the shear stress wave travelling in the radial direction. Initially, the vertical stress wave propagates with the dilation wave velocity in both the radial and vertical directions. After it interacts with the shear stress wave, the vertical stress wave starts to propagate with the shear wave velocity in the radial direction and with the axial wave velocity downwards. It is concluded that at the bottom of the soil plug, the interaction between the waves travelling in radial and vertical directions is important. The capabilities of several one‐dimensional pile‐in‐pile models to reproduce the driving response given by a two‐dimensional axisymmetric finite element model is studied. It is seen that when the base of the soil plug fails, a one‐dimensional pile‐in‐pile model can be used to achieve results in agreement with the finite element model. However, when the pile is unplugged, where the base of the soil plug does not fail, a reduced finite element mesh that permits the radial wave propagation inside the soil plug must be used. Copyright © 2001 John Wiley & Sons, Ltd.     

井底压差对垂直井井底应力场的影响研究

井底压差是影响钻井机械钻速的重要因素之一。实践证明欠平衡钻井与常规钻井相比机械钻速有显著提高,不同钻井条件下的井底应力场需要深入地研究。在分析井底岩石所受载荷的基础上,建立了轴对称条件下井底岩石流固耦合模型,并利用有限元方法进行求解,通过对不同井底压差下的井底岩石应力场进行数值模拟计算,针对井壁岩石的理论解和数值模拟结果进行对比验证模型的合理性,并对井底待破碎岩石进行了力学分析。结果表明,数值模拟结果与理论结果相符,验证了模型的正确性;井底待破碎岩石按照应力状态分为3个区域,分别为3向拉伸区、2向压缩区和3向压缩区。井底压差对井底应力场的影响研究为快速高效破岩提供理论基础。     

泥石流危险性分区及其在泥石流减灾中的应用

泥石流是一种突发性山地灾害,至今尚缺乏准确有效的预报方法。减灾工程也只能对一定规模的泥石流起到防御作用。泥石流危险性分区在泥石流减灾中具有重要作用。以泥石流运动数值模拟为基础,以数值模拟获得的流速和流深等参量为分区指标的危险性分区是泥石流危险性分区研究的重点和发展方向。其中,泥石流危险性动量和动能分区充分反映了泥石流的破坏能力,可以提供更为精确的定量化的分区结果。该类分区方法在城镇等有重要危害对象的泥石流减灾中具有广泛的应用。不仅可以应用于泥石流危险区和安全区的划定、泥石流灾害预估、泥石流临灾预案制定、泥石流抢险救灾方案制定和泥石流灾情评估等,还可以应用到山区土地利用规划、山区城镇建设规划和财产保险评估等领域,并起到防灾和减灾的作用。     

A numerical simulation of debris flow and its application

Debris flow is the flow of solid-fluid mixture and was treated as the flow of a continuum in routing in this study. A mathematical model was proposed to describe debris flow including deposition process and then solved numerically with suitable boundary conditions. Laboratory experiments were also conducted for comparison and calibration of the numerical results as well as for investigation of debris flow phenomena. The numerical model was also applied to simulate the debris flow caused by heavy rainfall in Tungman village of Hualien County located in the east of Taiwan on 23 June 1990. The simulated bed topographies in alluvial fan were in good agreement with those obtained from laboratory experiments and field observations.     

A numerical study of particle motion and two‐phase interaction in aeolian sand transport using a coupled large eddy simulation – discrete element method

Wind‐blown sand movement, considered as a particle‐laden two‐phase flow, was simulated by a new numerical code developed in the present study. The discrete element method was employed to model the contact force between sand particles. Large eddy simulation was used to solve the turbulent atmospheric boundary layer. Motions of sand particles were traced in the Lagrangian frame. Within the near‐surface region of the atmospheric boundary layer, interparticle collisions will significantly alter the velocity of sand. The sand phase is quite dense in this region, and its feedback force on fluid motion cannot be ignored. By considering the interparticle collision and two‐phase interaction, four‐way coupling was achieved in the numerical code. Profiles of sand velocity from the simulations were in good agreement with experimental measurements. The mass flux shows an exponential decay and is comparable to reported experimental and field measurements. The turbulence intensities and shear stress of sand particles were estimated from particle root‐mean‐square velocities. Distributions of slip velocity and feedback force were analysed to reveal the interactions between sand particles and the continuous fluid phase.     

Discrete-Element Computation of Averaged Tensorial Fields in Sand Piles Consisting of Polygonal Particles

This work is a contribution to the understanding of the mechanical properties of non-cohesive granular materials in the presence of friction and a continuation of our previous work (Roul et al. 2010) on numerical investigation of the macroscopic mechanical properties of sand piles. Besides previous numerical results obtained for sand piles that were poured from a localized source (“point source”), we here consider sand piles that were built by adopting a “line source” or “raining procedure”. Simulations were carried out in two-dimensional systems with soft convex polygonal particles, using the discrete element method (DEM). First, we focus on computing the macroscopic continuum quantities of the resulting symmetric sand piles. We then show how the construction history of the sand piles affects their mechanical properties including strain, fabric, volume fraction, and stress distributions; we also show how the latter are affected by the shape of the particles. Finally, stress tensors are studied for asymmetric sand piles, where the particles are dropped from either a point source or a line source. We find that the behaviour of stress distribution at the bottom of an asymmetric sand pile is qualitatively the same as that obtained from an analytical solution by Didwania and co-workers (Proc R Soc Lond A 456:2569–2588, 2000).     

Numerical modelling of fluid flow in microscopic images of granular materials

A program for the simulation of two‐dimensional (2‐D) fluid flow at the microstructural level of a saturated anisotropic granular medium is presented. The program provides a numerical solution to the complete set of Navier–Stokes equations without a priori assumptions on the viscous or convection components. This is especially suited for the simulation of the flow of fluids with different density and viscosity values and for a wide range of granular material porosity. The analytical solution for fluid flow in a simple microstructure of porous medium is used to verify the computer program. Subsequently, the flow field is computed within microscopic images of granular material that differ in porosity, particle size and particle shape. The computed flow fields are shown to follow certain paths depending on air void size and connectivity. The permeability tensor coefficients are derived from the flow fields, and their values are shown to compare well with laboratory experimental data on glass beads, Ottawa sand and silica sands. The directional distribution of permeability is expressed in a functional form and its anisotropy is quantified. Permeability anisotropy is found to be more pronounced in the silica sand medium that consists of elongated particles. Copyright © 2001 John Wiley & Sons, Ltd.     

基于FLAC^3D的三眼峪小流域危险性分析

通过对甘肃省舟曲县城三眼峪特大泥石流灾害的现场调查,采用有限差分FLAC^3D方法对该流域重力条件下侵蚀机理进行了探索,对流域应力场、位移场以及塑性区进行了数值模拟,对流域天然状态下危险性进行分析。结果表明,流域位移主要以压缩变形为主,最大位移主要出现在沟道两侧梁峁顶处的松散碎石土层内,局部区域已经破坏,且表现为竖向位移最大;由于碎石土、灰岩分界面的存在,使其附近区域的最大主应力方向较其他区域变化大而且迅速;流域内以剪切、张拉并存的破坏的塑性屈服模式为主,沟道两侧的表层碎石土层梁峁顶和梁峁坡上部同时处于张拉塑性屈服区和剪切塑性屈服区,破坏程度较为严重。建模方法和模拟结果可为泥石流预测提供科学依据。     

Arching within hydraulic fill stopes

Arching is a well known phenomenon, which effects stress developments which were investigated and compared using analytical and numerical solutions. Marston’s (1930) solution was extended to a generalised 3-dimensional rectangular stope and later modified for square and circular stopes for comparison with FLAC results. Aubertin et al. (2003) & Li et al. (2003) models were improved significantly by placing the backfill within narrow stopes as lifts or layers in numerical modelling where the normal stress variation with depth were found to be more realistic. The FLAC results were compared with analytical solutions which were developed by previous researchers and modified by the authors to evaluate the arching effects in backfilled placed in narrow and circular stopes. It appeared from the investigation herein that δ = 0.67 ϕ and K = K _o condition gives a very close match with the numerical model solutions obtained from FLAC. Many laboratory tests were conducted to find out friction angles for four Australian mines, which were between 30 and 49 degrees.     

Thresholds of intrabed flow and other interactions of turbidity currents with soft muddy substrates

Controlled laboratory experiments reveal that the lower part of turbidity currents has the ability to enter fluid mud substrates, if the bed shear stress is higher than the yield stress of the fluid mud and the density of the turbidity current is higher than the density of the substrate. Upon entering the substrate, the turbidity current either induces mixing between flow‐derived sediment and substrate sediment, or it forms a stable horizontal flow front inside the fluid mud. Such ‘intrabed’ flow is surrounded by plastically deformed mud; otherwise it resembles the front of a ‘bottom‐hugging’ turbidity current. The ‘suprabed’ portion of the turbidity current, i.e. the upper part of the flow that does not enter the substrate, is typically separated from the intrabed flow by a long horizontal layer of mud which originates from the mud that is swept over the top of the intrabed flow and then incorporated into the flow. The intrabed flow and the mixing mechanism are specific types of interaction between turbidity currents and muddy substrates that are part of a larger group of interactions, which also include bypass, deposition, erosion and soft sediment deformation. A classification scheme for these types of interactions is proposed, based on an excess bed shear stress parameter, which includes the difference in the bed shear stress imposed by the flow and the yield stress of the substrate and an excess density parameter, which relies on the density difference between the flow and the substrate. Based on this classification scheme, as well as on the sedimentological properties of the laboratory deposits, an existing facies model for intrabed turbidites is extended to the other types of interaction involving soft muddy substrates. The physical threshold of flow‐substrate mixing versus stable intrabed flow is defined using the gradient Richardson number, and this method is validated successfully with the laboratory data. The gradient Richardson number is also used to verify that stable intrabed flow is possible in natural turbidity currents, and to determine under which conditions intrabed flow is likely to be unstable. It appears that intrabed flow is likely only in natural turbidity currents with flow velocities well below ca 3·5 m s^?1, although a wider range of flows is capable of entering fluid muds. Below this threshold velocity, intrabed flow is stable only at high‐density gradients and low‐velocity gradients across the upper boundary of the turbidity current. Finally, the gradient Richardson number is used as a scaling parameter to set the flow velocity limits of a natural turbidity current that formed an inferred intrabed turbidite in the deep‐marine Aberystwyth Grits Group, West Wales, United Kingdom.     

Closed-Form Solutions for a Circular Tunnel in Elastic-Brittle-Plastic Ground with the Original and Generalized Hoek–Brown Failure Criteria

The paper presents a closed-form solution for the convergence curve of a circular tunnel in an elasto-brittle-plastic rock mass with both the Hoek–Brown and generalized Hoek–Brown failure criteria, and a linear flow rule, i.e., the ratio between the minor and major plastic strain increments is constant. The improvement over the original solution of Brown et al. (J Geotech Eng ASCE 109(1):15–39, 1983) consists of taking into account the elastic strain variation in the plastic annulus, which was assumed to be fixed in the original solution by Brown et al. The improvement over Carranza-Torres’ solution (Int J Rock Mech Min Sci 41(Suppl 1):629–639, 2004) consists of providing a closed-form solution, rather than resorting to numerical integration of an ordinary differential equation. The presented solution, by rigorously following the theory of plasticity, takes into account that the elastic strain components change with radial and circumferential stress changes within the plastic annulus. For the original Hoek–Brown failure criterion, disregarding the elastic strain change leads to underestimate the convergence by up to 55%. For a rock mass failing according to the generalized Hoek–Brown failure criterion, using the original failure criterion leads to a high probability (97%) of underestimating the convergence by up to 100%. As a consequence, the onset or degree of squeezing may be underestimated, and the loading on the support/reinforcement calculated with the convergence/confinement method may be largely underestimated.     

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.     

统一弹塑性本构模型在FLAC3D中的计算格式

统一强度理论是一个有着坚实理论基础、可表达各种岩土材料强度特性且被广泛应用的新的强度理论。FLAC3D是一个具有强大的计算分析功能且专门针对岩土工程问题开发并被广泛应用的数值分析软件。若能将二者结合起来,无疑会大大促进岩土工程领域相关问题的解决。针对这一问题,研究了统一弹塑性本构模型在有限差分方法中应用的基本理论格式。根据FLAC3D软件中UDM接口计算格式的要求,将统一强度理论引入其中,详细推导了统一弹塑性本构模型在FLAC3D软件中应用的计算公式。由于统一屈服面在应力空间内由12个面组成,在应力角为[0,π/3]的范围内,统一屈服面由两个相交面组成,为计算塑性因子及处理计算中应力超出屈服面的应力调整问题,利用应力角分析了应力空间的分区,推导了区域分界面的公式。最后,编制了相应的UDM接口程序对一圆形隧洞进行了弹塑性分析,对比了计算结果与解析解,结果很好地验证了此计算格式及相应接口程序的正确性。统一强度理论和FLAC3D软件的结合,将使二者的优点得以充分发挥,以解决更多的岩土工程问题。     

砂质碎屑流和底流改造——部分传统浊积岩成因新解

以鲍马层序为代表的古典浊积岩,在地质学界曾广泛流行。近年来的研究表明,真正的浊积岩是很少发生的。由于浊流和碎屑流具有不同的流变学特点,人们对部分传统浊积岩提出了砂质碎屑流和底流改造成因及其鉴别标志,并对传统海底浊积扇模式在进行重新研究的基础上提出了海底碎屑流模式。     

Hydraulic and numerical study on the generation of a subaqueous landslide-induced tsunami along the coast

By carrying out the hydraulic experiments in a one-dimensional open channel and two-dimensional basin, we clarified the process of how a landslide on a uniform slope causes the generation of a tsunami. The effect of the interactive force that occurs between the debris flow layer and the tsunami is significant in the generation of a tsunami. The continuous flow of the debris into the water makes the wave period of the tsunami short. The present experiments apply numerical simulation using the two-layer model with shear stress models on the bottom and interface, and the results are compared. The simulated debris flow shows good agreement with the measured results and ensures the rushing process into the water. We propose that the model use a Manning coefficient of 0.01 for the smooth slope and 0.015 for the rough slope, and a horizontal viscosity of 0.01 m²/s for the landslide; an interactive force of 0.2 for each layer is recommended. The dispersion effect should be included in the numerical model for the propagation from the shore.     

Quantifying a critical marl thickness for vertical fracture extension using field data and numerical experiments

In fractured reservoirs characterized by low matrix permeability,fracture networks control the main fluid flow paths.However,in layered reservoirs,the vertical extension of fractures is often restricted to single layers.In this study,we explored the effect of changing marl/shale thickness on fracture extension using comprehensive field data and numerical modeling.The field data were sampled from coastal exposures of Liassic limestone-marl/shale alternations in Wales and Somerset(Bristol Channel Basin,UK).The vertical fracture traces of more than 4000 fractures were mapped in detail.Six sections were selected to represent a variety of layer thicknesses.Besides the field data also thin sections were analyzed.Numerical models of fracture extension in a two-layer limestone-marl system were based on field data and laboratory measurements of Young's moduli.The modeled principal stress magnitude σ_3 along the lithological contact was used as an indication for fracture extension through marls.Field data exhibit good correlation(R~2=0.76) between fracture extension and marl thickness,the thicker the marl layer the fewer fractures propagate through.The model results show that almost no tensile stress reaches the top of the marl layer when the marls are thicker than 30 cm.For marls that are less than 20 cm,the propagation of stress is more dependent on the stiffness of the marls.The higher the contrast between limestone and marl stiffness the lower the stress that is transmitted into the marl layer.In both model experiments and field data the critical marl thickness for fracture extension is ca.15-20 cm.This quantification of critical marl thicknesses can be used to improve predictions of fracture networks and permeability in layered rocks.Up-or downsampling methods often ignore spatially continuous impermeable layers with thicknesses that are under the detection limit of seismic data.However,ignoring these layers can lead to overestimates of the overall permeability.Therefore,the understanding of how fractures propagate and terminate through impermeable layers will help to improve the characterization of conventional reservoirs.