Three Dimensional Numerical Analysis of Underground Bifurcated Tunnel

As far as the bifurcated tunnel of underground engineering is concerned, it is usually used in the water conveyance system. Due to the complexity of underground rock masses and concrete lining, researches on mechanical characteristic and stability of the bifurcation tunnel have attracted more and more attention in the geotechnical field. In order to understand bifurcated tunnel in detail, three-dimensional (3D) numerical method is applied to solve the above key subjects by simulating a practical project. Furthermore, sub-model technology is applied to analyze the intersection position, corresponding deformation and stress results in the practical condition. Meanwhile, 3D excavation and support calculation under four conditions have been simulated based on 3D self-compiled code and Ansys software. In addition, the paper plays emphasis on the stress, displacement analysis considering different stress releasing ratio instead of rheological analysis, and the results corresponding with the fact indicate the feasibility of 3D elasto–visco–plastic code.     

Bifurcation modeling in geomaterials: From the second‐order work criterion to spectral analyses

The present paper investigates bifurcation analysis based on the second‐order work criterion, in the framework of rate‐independent constitutive models and rate‐independent boundary‐value problems. The approach applies mainly to nonassociated materials such as soils, rocks, and concretes. The bifurcation analysis usually performed at the material point level is extended to quasi‐static boundary‐value problems, by considering the stiffness matrix arising from finite element discretization. Lyapunov's definition of stability (Annales de la faculté des sciences de Toulouse 1907; 9 :203–274), as well as definitions of bifurcation criteria (Rice's localization criterion (Theoretical and Applied Mechanics. Fourteenth IUTAM Congress, Amsterdam, 1976; 207–220) and the plasticity limit criterion are revived in order to clarify the application field of the second‐order work criterion and to contrast these criteria. The first part of this paper analyses the second‐order work criterion at the material point level. The bifurcation domain is presented in the 3D stress space as well as 3D cones of unstable loading directions for an incrementally nonlinear constitutive model. The relevance of this criterion, when the nonlinear constitutive model is expressed in the classical form (dσ = Mdε) or in the dual form (dε = Ndσ), is discussed. In the second part, the analysis is extended to the boundary‐value problems in quasi‐static conditions. Nonlinear finite element computations are performed and the global tangent stiffness matrix is analyzed. For several examples, the eigenvector associated with the first vanishing eigenvalue of the symmetrical part of the stiffness matrix gives an accurate estimation of the failure mode in the homogeneous and nonhomogeneous boundary‐value problem. Copyright © 2008 John Wiley & Sons, Ltd.     

基于分叉理论的岩体局部化现象研究

为了研究岩体局部化现象,基于轴向压杆稳定实例,运用数理和力学知识,根据Drucker-Prager（简称D-P）准则,推导了基于分叉理论的通用表达式。在此基础上,结合有限差分理论得出FLAC3D计算中岩体内发生分叉的判定条件,之后在自带的FISH编译语言的辅助下,在FLAC3D中建立数值分析模型,采用D-P本构关系,研究了岩体局部化现象及其发展过程。研究结果表明：在FLAC3D中以20 000步为间隔所得的不平衡力比率模型系统逐渐趋于稳定,不平衡力比率逐渐降低,局部化现象的规模逐渐扩大至稳定,清晰度逐渐增加,证实了分叉理论能够有效地分析研究岩体局部化现象。给出的岩体局部分叉判定条件和数值验证,较为科学合理地阐述了常用分叉判定条件的由来,为解释岩体局部化现象的机制和过程提供了理论支撑,从而对岩体局部化现象研究起到了一定的推进作用。     

A two-mechanism soil–structure interface model for three-dimensional cyclic loading

A three-dimensional (3D) soil–structure interface model is proposed within the two-mechanism constitutive theory and bounding surface theory originally established for soils. The proposed model has two main characteristics: first, the model is formulated based on two different and superposed deformation mechanisms. The first mechanism is due to the stress ratio increment, and the second is due to the normal stress increment. Each mechanism induces a shear strain component and a normal strain component. The proposed model can be reduced to the conventional single-mechanism interface model. Second, the plastic modulus and stress dilatancy are defined using the bounding surface theory. The plastic flow rule under cyclic loading is modified and assumed to be dependent on both the stress state of the mapping point and the stress reversal loading direction. The proposed model was validated against the available 3D interface tests and was found to satisfactorily reflect the salient features of the interfaces under monotonic and cyclic loading paths with different normal boundaries. The problem in which the elastic normal stiffness in conventional single-mechanism interface models is often underestimated to enhance the simulation performance under varying normal stress conditions is solved by incorporating the second mechanism. And the effect of the second mechanism on the modeling behavior is discussed. The modified plastic flow direction accurately simulates the 3D cyclic shear response, and the difference between the model simulation and test result increases with the number of cycles by use of the plastic flow direction defined in conventional bounding surface theory.     

Material stability analysis of rock joints

For prediction of rockfalls, the failure of rock joints is studied. Considering these failures as constitutive instabilities, a second‐order work criterion is used because it explains all divergence instabilities (flutter instabilities are excluded). The bifurcation domain and the loading directions of instabilities, which fulfill the criterion, are determined for any piecewise linear constitutive relation. The instability of rock joints appears to be ruled by coupling features of the behavior (e.g., dilatancy). Depending on the loading parameters, instabilities can lead to failure, even before the plastic limit criterion. Results for two given constitutive relations illustrate the approach. Some given loading paths are especially considered. Constant volume (undrained) shear and τ‐constant paths are stable or not depending on the link between the deviatoric stress and strain along undrained paths, as found for soils. Some unstable loading paths are illustrated. Along these paths, failure before the plastic limit criterion is possible. The corresponding failure rules are determined. Copyright © 2012 John Wiley & Sons, Ltd.     

Diffuse failure in geomaterials: Experiments,theory and modelling

This paper presents a synthesis of the works performed by various teams from France, Italy and Canada around the question of second‐order work criterion. Because of the non‐associative character of geomaterials plastic strains, it is now recognized that a whole bifurcation domain exists in the stress space with various possible modes of failure. In a first part these failure modes are observed in lab experimental tests and in discrete element modelling. Then a theoretical study of second‐order work allows to establish a link with the kinetic energy, giving a basis to explain the transition from a prefailure (quasi)static regime to a postfailure dynamic regime. Eventually the main features of geomaterials failure are obtained by applying second‐order work criterion to five different constitutive rate‐independent models—three being phenomenological and two micromechanical. As a whole this paper tries to gather together all the elements for a proper understanding and use of second‐order work criterion in geomechanics. Copyright © 2010 John Wiley & Sons, Ltd.     

Non‐coaxial elasto‐plasticity model and bifurcation prediction of shear banding in sands

Numerous constitutive models built on coaxial theory and validated under axi‐symmetric condition often describe the stress–stain relationships and predict the inceptions of shear banding in sands inaccurately under true triaxial condition. By adopting an elaborated Mohr–Coulomb yield function and using non‐coaxial non‐associated flow rule, a 3D non‐coaxial elasto‐plasticity model is proposed and validated by a series of true triaxial tests on loose sands. The bifurcation analysis of true triaxial tests on dense sands predicts the influence of the intermediate principal stress ratio on the onset of shear band accurately. The failure of soils is shown to be related to the formation of shear band under most intermediate principal stress ratio conditions except for those which are close to the axi‐symmetric compression condition. Copyright © 2009 John Wiley & Sons, Ltd.     

Transient bifurcation condition of partially saturated porous media at finite strain

Bifurcation of unsaturated soils into a localized shear band is a ubiquitous failure mode of partially saturated soils. The density and degree of saturation have major impacts on the inception of localized deformations in unsaturated soils. Unsaturated fluid flow may dramatically change the density and degree of fluid saturation of unsaturated soils. Therefore, the unsaturated fluid flow is a potential trigger for shear banding in such materials. In this paper, we derive a simplified bifurcation condition of localized deformation in unsaturated soils under the local transient condition at finite strain. This transient bifurcation condition is implemented into a nonlinear finite element code to study the inception of localized deformation in unsaturated soil specimens. Numerical simulations are conducted to study the impact of soil fabrics of density, a ‘bonding’ variable, and intrinsic permeability on the inception of localized failures via the transient bifurcation criterion. Mesh sensitivity analysis is performed to demonstrate the viscosity effect of unsaturated fluid flow on the localized deformation. Numerical simulations demonstrate that the transient bifurcation condition can detect the localized deformation triggered by the internal unsaturated fluid flow process in unsaturated soils. Copyright © 2016 John Wiley & Sons, Ltd.     

Uniqueness and existence in plasticity models for unsaturated soils

The paper investigates the mathematical structure of plasticity models for unsaturated soils and provides a strategy to capture the loss of uniqueness of the incremental solution upon loading and/or wetting paths. To derive bifurcation conditions in simple analytical form, the analysis is restricted to isotropic stress states. This choice has allowed the inspection of the most common classes of constitutive models through a unified notation, as well as the study of different forms of coupling between plasticity and state of saturation. It is shown that, similar to saturated soil plasticity, the loss of admissibility of the plastic solution is governed by critical values of the hardening modulus. At variance with the classical case, however, these moduli can be positive even if the plastic flow rule is associated (bifurcation in the hardening regime). The paper shows that such non-trivial features derive from hydro-mechanical coupling, i.e. they depend on the approach used to reproduce suction effects and evolving retention properties. In other words, although the problem of loss of uniqueness affects all classes of plasticity models for unsaturated soils, different constitutive assumptions may not have the same outcome in terms of bifurcation potential. As a result, new concepts are introduced to compare the mathematical robustness of the different constitutive approaches, as well as to interpret their predictions in the light of precise bifurcation criteria.     

节理岩体能量损伤本构模型与工程应用

本文根据Betti能量互易定理,节理岩体能量损伤演化方程,广义正交法则和塑性损伤一致性条件建立了节理岩体的能量损伤本构模型,并将该模型编制成三维非线性有限元计算程序应用于溪洛渡地下厂房洞室群施工开挖稳定性分析,获得了较为满意的计算结果。     

On the transition between 2D and 3D dunes

Sediment transport in sand-bedded alluvial channels is strongly conditioned by bedforms, the planimetric morphology of which can be either two- or three-dimensional. Experiments were undertaken to examine the processes that transform the bed configuration from two-dimensional (2D) dunes to three-dimensional (3D) dunes. A narrowly graded, 500 μm size sand was subjected to a 0·15 m deep, non-varying mean flow ranging from 0·30 to 0·55 m sec⁻¹ in a 1 m wide flume. Changes in the planimetric configuration of the bed were monitored using a high-resolution video camera that produced a series of 10 sec time-lapsed digital images. Image analysis was used to define a critical value of the non-dimensional span (sinuosity) of the bedform crestlines that divides 2D forms from 3D forms. Significant variation in the non-dimensional span is observed that cannot be linked to properties of the flow or bedforms and thus appears random. Images also reveal that, once 2D bedforms are established, minor, transient excesses or deficiencies of sand are passed from one bedform to another. The bedform field appears capable of absorbing a small number of such defects but, as the number grows with time, the resulting morphological perturbations produce a transition in bed state to 3D forms that continue to evolve, but are pattern-stable. The 3D pattern is maintained by the constant rearrangement of crestlines through lobe extension and starving downstream bedforms of sediment, which leads to bifurcation. The experiments demonstrate that 2D bedforms are not stable in this calibre sand and call into question the reliability of bedform phase diagrams that use crestline shape as a discriminator.     

形成湿地的“淹埋深-历时-频率”阈值研究进展

湿地水文是形成湿地的发生学因素,“淹埋深-历时-频率”(S, D, F)是表征湿地水文的特征指标。为了推进(S, D, F)阈值研究理论和方法的发展,总结了(S, D, F)阈值研究的一致结论:S阈值确定的原则是水饱和至地表;D阈值确定的原则是从水饱和开始到厌氧环境形成的时滞;F阈值应取≥50%。通过对这些一致结论的分析后认为:D阈值确定的原则存在着科学性问题;目前(S, D, F)阈值研究理论针对正常情况下湿地类型,不适合非正常情况下的湿地类型。为了科学地解决(S, D, F)阈值问题,必须构建新的(S, D, F)阈值确定的理论和方法,包括正常情况下和非正常情况(S, D, F)阈值确定的理论和方法以及(S, D, F)阈值科学性检验的理论和方法。     

饱和悬浮塑料砂流动变形可视环剪试验研究

为了分析剪切条件下零有效应力状态砂土的流动变形规律,对传统环剪仪进行了试样可视化改造,研制了透明环形剪切盒,通过对剪切盒膨胀性能分析及与标准砂的剪切试验对比,验证了环剪装置改造的合理性。通过分析不排水条件下饱和悬浮塑料砂的剪应力-应变曲线,发现其剪切强度具有应变软化的特性。通过分析环剪仪中饱和悬浮塑料砂试样的有效应力,证实了可视环剪试验中的土体基本处于零有效应力状态。开展了饱和悬浮塑料砂的流动变形可视环剪试验研究,结果表明：饱和悬浮塑料砂在固结不排水条件下剪切变形不连续,直接在剪切面发生断裂;在不固结不排水的条件下饱和悬浮塑料砂的剪切变形表现出流动形态,且与剪切速率有关：在低剪切速率下,剪切变形仅在剪切面处形成具有曲线轨迹的流动变形而在其他区域不发生变形;而在高剪切速率下剪切变形为整体的倾斜变形,符合黏性流体的流动变形特征。     

Second‐order work analysis for granular materials using a multiscale approach

It has been established that the second‐order work criterion is a general necessary condition for all instabilities by divergence in rate‐independent granular materials. The relation between the values of discrete second‐order work at the intergranular contact point level and its global macroscopic value is recalled at the beginning of this paper. Then, the basic purpose of the paper is tackled by an analysis of the main features of second‐order work criterion in relation with the granular microstructure. For that, it is considered a novel micromechanical model (the so‐called ‘H‐microdirectional model’), which has the property to involve three scales: grain scale, mesoscale with a specific granular configuration and continuum mechanics macroscale. Eventually, these exhibited features (a bifurcation stress domain and some instability cones) are qualitatively compared with the ones provided by direct numerical simulations issued from a discrete element model. The ultimate goal is to analyse what happens at the granular scale, when the macrosecond‐order work is vanishing at the macrolevel. Copyright © 2013 John Wiley & Sons, Ltd.     

模拟土体本构特性的热力学方法

简要介绍了建立岩土材料弹塑性本构模型的热力学方法。它不仅具有紧凑的数学结构,而且自动满足热力学定律,仅从两个热力学势函数,即自由能函数与耗散增量函数出发,就足以导出弹塑性理论必须的屈服条件,流动法则,硬化定律和弹性定律。通过理论证明指出,只要耗散增量函数依赖于当前应力,流动法则必然是非关联的,岩土材料的摩擦特性与非关联流动法则密不可分。介绍该方法在三维模型,岩土材料的微细观力学特性,应力应变的均匀化以及剪胀和各向异性方面应用的主要研究进展,并对一些重要的概念,诸如“储存的塑性功”,“Reynolds-Taylor状态”等,进行分析与解释。最后给出近期需进一步深入研究的几点建议。     

海南石碌铁矿床构造变形特征及其与铁多金属成矿富集的关系

海南石碌铁矿曾被誉为“亚洲最大的富铁矿”,其形成、定位与褶皱变形及伴随的剪切、塑性流动等构造密切相关。石碌矿区构造变形大致分为早期(D1)复式向斜的形成时期、晚期(D2)褶皱叠加和剪切变形时期。后者又进一步分为韧-脆性变形(D2a),脆-韧性变形、层间滑脱断层形成(D2b)和脆性变形、矿体破坏(D2c)3个阶段。石碌铁矿床受NW—SE向复式向斜及其与NE—NNE向褶皱叠加所形成的横跨或斜跨褶皱的严格控制;褶皱过程所伴随的剪切变形和高温塑性流动是富铁矿形成的重要机制。复式向斜轴部,尤其是2期向斜褶皱轴的叠加部位往往可发现厚大的富铁矿体。     

Bifurcation and second‐order work in geomaterials

In this paper, the ability of a material rate‐independent system to evolve toward another mechanical state from an equilibrium configuration, with no change in the control parameters, is investigated. From a mechanical point of view, this means that the system can spontaneously develop kinetic energy with no external disturbance from an equilibrium state, which corresponds to a particular case of bifurcation. The existence of both conjugate incremental strain and stress such that the second‐order work vanishes is established as a necessary and sufficient condition for the appearance of this bifurcation phenomenon. It is proved that this fundamental result is independent of the constitutive relation of the rate‐independent material considered. Then the case of homogeneous loading paths is investigated, and, as an illustration, the subsequent results are applied to interpret the well‐known liquefaction observed under isochoric triaxial loading conditions with loose granular materials. Copyright © 2006 John Wiley & Sons, Ltd.     

三峡深水高土石围堰工程的渗流研究

为了对三峡工程二期高土石围堰防渗设施的布置方案及其阻渗效果进行比较,采用有限元法对二期高土石围堰在不利的运行工况下的渗流场进行数值分析.比较了双排混凝土防渗墙方案、单排塑性混凝土低防渗墙加土工膜斜墙和单排厚塑性混凝土防渗墙等3种方案,分别采用恒定与非恒定模型计算,对立面二维和三维绕渗及防渗墙局部开裂等不利工况分别进行了数值模拟.结果表明:三方案均可有效抑制渗流场,双排混凝土防渗墙的防渗效果最好;墙体的局部开裂仅对局部区域的流场有影响.非恒定数值分析表明,堰体、基础不均匀沙石料及基坑抽水速度对渗流场影响极大,为保证堰体稳定,应限制基坑水位降落速度小于2 m/d.     

考虑脆性损伤和渗流的圆形水工隧洞解析解

针对脆性岩石的损伤特性引入脆性损伤模型,导出三维各向同性损伤情况下的损伤变量表达式。在考虑损伤、不同工况主应力对应关系和渗流作用影响下,推导出圆形水工隧洞洞周围岩的弹塑性应力和塑性半径的解析解。对岩石的脆性特性的研究表明,考虑脆性损伤时的洞周围岩塑性半径比不考虑时有一定程度提升;围岩脆性特征越明显,塑性半径越大。无论洞周围岩应力状态如何,考虑渗流影响后,洞周围岩塑性区范围均会有所扩展。在内水压从0逐渐增加过程中,塑性半径随内水压增加而逐渐减小,直至塑性区消失,在超出临界内压后,塑性区重新出现,塑性半径逐渐增大。     

The effect of flood intermittency on bifurcations in fluviodeltaic systems: Experiment and theory

The dimensions and organization of deltaic islands and channels dictate delta morphology. This study presents experimental results modelling deposition at a river mouth and flow bifurcation around delta islands. Mouth bar formation and channel bifurcation is achieved in a laboratory setting by alternating input of suspended load transport and bedload transport. These two modes of transport produce two characteristic deposits with different advection lengths. Suspended load transport creates a steep deposit far from the inlet, while bedload creates a low angle, levéed deposit near the inlet. This study found that flow bifurcations occur where the proximal and distal deposits encroach on one another; and determined that there is a relationship between the frequency of suspended load transport and the length to channel bifurcation. Frequent flooding causes shorter length to bifurcations, whereas infrequent flooding causes greater length to bifurcations. This work overturns the hitherto understood mechanism of bifurcation location as a function of only high-transport conditions. Instead, the interactions between the sediment transport and deposition from normal flow and large-scale flooding events dictate delta island morphology.