Non‐local damage modelling of concrete: a procedure for the determination of model parameters

This paper presents a procedure for the determination of parameters of non‐local damage models. This is to assure a consistent response of a non‐local damage model, as choice of the internal length and other parameters of the model are varied. Correlations between the internal length and other parameters governing the local constitutive behaviour of the model are addressed and exploited. Focus is put on the relationship between the internal length of the non‐local model and the width of the fracture process zone. Numerical examples are used to demonstrate the rigour of the proposed method. Copyright © 2006 John Wiley & Sons, Ltd.     

Stability analysis of membrane‐reinforced curved earth retaining walls,using a multiphase approach

The stability analysis of curved earth retaining walls, stabilized by reinforcing membranes, is investigated by means of a multiphase model developed in the framework of the yield design approach. This model is an extension of that previously developed for soils reinforced by linear inclusions. It combines the advantage of a homogenization approach in terms of improved computational efficiency, with its capability to account for a specific soil–reinforcement failure condition, in a rational and systematic way. Application of this model is performed on the illustrative example of a cylindrical‐reinforced retaining wall by means of the kinematic approach of yield design, which provides upper bound estimates for the retaining wall stability factor. Nondimensional charts are finally presented assessing the influence of relevant parameters such as the curvature of the wall, the length of the reinforcing membranes or the reinforcement pull‐out resistance. Copyright © 2009 John Wiley & Sons, Ltd.     

劈裂注浆复合土体三维等效弹性模型理论研究

劈裂注浆可以有效改善土体的变形参数,大大降低土体在受力状态改变时的变形量,对劈裂注浆后复合土体的等效变形参数进行研究十分重要。在综合分析劈裂注浆扩散机制和工程应用实际的基础上,基于均质化理论提出了劈裂注浆后复合土体的三维单元体几何模型,按等效原则给出了浆-土体积及受力分配关系模型图;接着基于横向各向同性本构关系推导了模型的等效弹性模量和等效泊松比的解析解。然后采用有限元方法取得了模型特定条件下的等效弹性模量和等效泊松比,并与解析结果进行对比分析。最后把模型和相应的解析结果引入Flac3D岩土工程专业分析软件,结合一个热力隧道工程实例对隧道劈裂注浆后关键位置的沉降进行预测分析,并与实测值进行了对比。研究表明：对所提出的计算模型,解析计算与有限元方法计算结果吻合度较高,说明了解析结果的正确性;基于该模型及其解析结果得到的隧道开挖后的沉降预测值与实测值具有良好的一致性,说明所提出的模型和相应的解析计算方法具有一定的可靠性和实用性。     

锚杆参数对锚固边坡抗震性能影响研究

分析了锚固边坡锚杆与岩土体在动力荷载下的黏结机制,提出了一种简化的锚杆岩土体相互作用的动力学模型,推导了锚杆岩土体系统的动力响应解析解。结合该简化的锚杆岩土体动力学模型,借助ANSYS有限元程序建立锚固边坡数值模型,分析了在地震荷载作用下,锚杆参数（弹性模量、横截面积、长度、间距）变化对锚固边坡抗震性能（位移、加速度、应力）的影响。分析结果表明：锚杆参数对锚固边坡抗震性能影响强弱不一,增大横截面积可以有效提高边坡抗震稳定性,减小间距也可以改善边坡抗震性能;而弹性模量和长度对边坡抗震性能提高不大,甚至还会起反作用。这些结论对于锚固边坡的抗震设计有一定的借鉴作用。     

基于间隙参数模型的盾构隧道周围土体位移分析

基于一种盾构加载模型并引进间隙参数的概念,对盾构隧道周围土体位移的机理及预测方法进行了总结。盾构相对于设计轴线的偏移是影响土体变形的主要原因之一。确定间隙参数是应用解析方法的关键步骤,提出了一个半经验、半解析方法,并针对一个工程实例把现场测量值和理论计算值进行了对比分析。     

Study of bond stress–slip relationship and radial dilation in prestressed concrete

This paper presents two test procedures for evaluating the bond stress–slip and the slip–radial dilation relationships when the prestressing force is transmitted by releasing the steel (wire or strand) in precast prestressed elements. The bond stress–slip relationship is obtained with short length specimens, to guarantee uniform bond stress, for three depths of the wire indentation (shallow, medium and deep). An analytical model for bond stress–slip relationship is proposed and compared with the experimental results. The model is also compared with the experimental results of other researchers. Since numerical models for studying bond‐splitting problems in prestressed concrete require experimental data about dilatancy angle (radial dilation), a test procedure is proposed to evaluate these parameters. The obtained values of the radial dilation are compared with the prior estimated by numerical modelling and good agreement is reached. Copyright © 2012 John Wiley & Sons, Ltd.     

Homogenization framework for three‐dimensional elastoplastic finite element analysis of a grouted pipe‐roofing reinforcement method for tunnelling

This paper deals with the grouted pipe‐roofing reinforcement method that is used in the construction of tunnels through weak grounds. This system consists on installing, prior to the excavation of a length of tunnel, an array of pipes forming a kind of ‘umbrella’ above the area to be excavated. In some cases, these pipes are later used to inject grout to strengthen the ground and ‘connect’ the pipes. This system has proven to be very efficient in reducing tunnel convergence and water inflow when tunnelling through weak grounds. However, due to the geometrical and mechanical complexity of the problem, existing finite element frameworks are inappropriate to simulate tunnelling using this method. In this paper, a mathematical framework based on a homogenization technique to simulate ‘grouted pipe‐roofing reinforced ground’ and its implementation into a 3‐D finite element programme that can consider stage construction situations are presented. The constitutive model developed allows considering the main design parameters of the problem and only requires geometrical and mechanical properties of the constituents. Additionally, the use of a homogenization approach implies that the generation of the finite element mesh can be easily produced and that re‐meshing is not required as basic geometrical parameters such as the orientation of the pipes are changed. The model developed is used to simulate tunnelling with the grouted pipe‐roofing reinforcement method. From the analyses, the effects of the main design parameters on the elastic and the elastoplastic analyses are considered. Copyright © 2004 John Wiley & Sons, Ltd.     

Scratch hardness–strength solutions for cohesive‐frictional materials

In this paper we develop analytical solutions for scratch hardness–strength relations for cohesive‐frictional materials of the Mohr–Coulomb and Drucker–Prager type. Based on the lower bound yield design approach, closed‐form solutions are derived for frictionless scratch devices, and validated against computational upper bound and elastoplastic finite element solutions. The influence of friction at the blade–material interface is also investigated, for which a simple computational optimization is proposed. Illustrated for scratch tests on cement paste, we show that the proposed solutions provide a convenient way to determine estimates of cohesion and friction parameters from scratch data, and may serve as a benchmark to identify the relevance of strength models for scratch test analysis. Copyright © 2011 John Wiley & Sons, Ltd.     

An analytical solution for stresses induced by a post-tensioned anchor in rocks containing two perpendicular joint sets

A new analytical approach for the analysis of stress around a post-tensioned anchor in rock with two perpendicular joint sets is presented. The solution considers nonlinear shear stresses developed along the anchor bond length as well as debonding at the tendon–grout interface. The bearing plate effects are also considered in the analyses. The following assumptions are made: (1) homogeneous and orthotropic elastic rock; (2) half-space rock mass with plane strain conditions; (3) trilinear bond-slip model for the behavior of the tendon–grout interface; (4) an elastic anchor. The employed methodology is to decompose the anchor problem into two problems of simpler loadings: stresses produced under a rigid bearing plate and stresses induced by interfacial shear stresses mobilized along the bond length. Based on the proposed solution, an illustrative example is given and the results show that a large compressive region is formed around the anchor free length. Tensile stress concentrations are also observed around the bond length and near the rock surface outside the bearing plate. A parametric study is also conducted to investigate the effects of joint properties on the induced stresses. The results show that as the number of joints intersecting the anchor axis increases, the magnitude of compressive stresses in the free zone decreases, while the size of compression zone around the anchor increases. To verify the results of the analytical approach, a comparison is made with numerical results obtained by using the finite element method. The analytical solutions compare very well with the results obtained by numerical method.     

The determination of the epicentre by a vectorial modelling of macroseismic intensity distribution

A vectorial modelling of observed macroseismic intensity aimed at the analytical determination of the epicentre is proposed here. The methodology is based on the determination of a plane system of vectors which characterises the macroseismic intensity distribution. The epicentre of each seismic event considered is determined as the centre of this vector system by an analytical expression which is independent from all possible directions of seismic energy propagation. The analysis of the intensity distribution is carried out by a new model called a macroseismic plane, different from the one known as macroseismic field, formed by a set of small areas built around the observed intensity points; hence its name.With the proposed methodology, some earthquakes in southern Italy and eastern Sicily are analysed calculating their epicentres, also for distributions of observed intensity which are particularly complex.     

复合衬层中变系数有机污染物迁移规律分析

土工膜和下伏压实黏土组成的复合衬层已被广泛用作填埋场的防渗屏障系统。有机污染物在复合衬层中迁移时,其主要运移机制是扩散作用。假设有机污染物在土工膜中稳态扩散,并假设在下伏衬层中扩散系数为迁移距离的线性函数,且考虑降解作用的影响,建立了有机污染物在复合衬层中的一维扩散模型,针对零浓度下边界条件,获得了模型解析解。基于该解析解,分析讨论了相关参数的敏感性。结果表明,相关参数对计算结果影响很大,降解半衰期对污染物运移同样存在较大影响。该解析模型仅适用于有机污染物在两层复合衬垫中的迁移问题,可为填埋场的初步设计提供参考     

Finite element simulations of the behavior of piled‐raft foundations using a multiphase model

Finite element simulations of the behavior of a piled raft foundation have been carried out using a multiphase model conceived as an improved homogenization approach. According to this model, the ground reinforced by a group of piles is treated as a homogeneous continuous medium. In this approach, no specific interface elements are necessary to account for the mechanical interaction between the piles and the ground: this interaction is described by means of two scalar parameters, one stiffness parameter and one which can easily be derived from the maximum ground‐pile friction. The implementation of the model into a finite element code provides an efficient tool for the analysis of the influence of the pile number or length on the settlement and bearing capacity of a square piled raft foundation and of the way the total applied load is shared between the raft and the piles. Results are compared with a standard 3D finite element analysis. The comparison highlights the fact that the proposed approach remains to be improved to account for tip resistance. Copyright © 2012 John Wiley & Sons, Ltd.     

粘弹性地基中有限长桩横向受迫振动问题解析解

基于动力Winkler模型(BDWF)对粘弹性地基中有限长桩的动力响应进行求解,得到了其在桩顶受水平循环荷载作用下的简单形式的时域解析解。通过迭加原理得到了半正弦激振条件下的瞬态解析解,经过简单的数学运算得到了桩头动力阻抗系数。通过对稳态激振时幅频曲线和相频曲线的分析,讨论了桩侧土阻尼系数和刚度系数对共振频率和共振峰值的影响规律。并讨论了桩底边界条件对水平动力响应的影响。本文为单桩的动力分析提供了一种简单的解析分析方法,可为工程设计提供初步参考。     

基于Drucker-Prager准则的扩孔器单齿正交切削岩土三维力学模型

扩孔器是定向穿越扩孔钻进阶段进行岩土切削的主要工具,合理的结构设计是提高扩孔时效、降低工程复杂问题出现几率的主要手段,且由于扩孔阶段是决定定向穿越工程工期的最主要阶段,因此,对扩孔器切削齿的切削效率进行深入研究是非常必要的。从单齿正交切削岩土的机制入手,综合考虑岩土与切削齿面之间的接触摩擦及岩土剪切区的剪切作用,分析了单齿正交切削岩土的力学特征,基于能够考虑中间主应力对岩土破坏影响的Drucker-Prager强度准则,推导了单齿正交切削岩土的主切削力与主切削力矩计算公式;在利用有限元计算程序对所推导公式进行验证的基础上,优选了剪切角计算模型,获得了不同岩土材料时的主切削力解析解,并分析了不同齿前角对主切削力的影响规律,从而获得单齿正交切削岩土的最优齿前角,为进一步研究扩孔器整体结构设计并针对不同岩土介质设计出适应性更强的扩孔器结构奠定基础。     

Updating the models and uncertainty of mechanical parameters for rock tunnels using Bayesian inference

Rock mechanical parameters and their uncertainties are critical to rock stability analysis, engineering design, and safe construction in rock mechanics and engineering. The back analysis is widely adopted in rock engineering to determine the mechanical parameters of the surrounding rock mass, but this does not consider the uncertainty. This problem is addressed here by the proposed approach by developing a system of Bayesian inferences for updating mechanical parameters and their statistical properties using monitored field data, then integrating the monitored data, prior knowledge of geotechnical parameters,and a mechanical model of a rock tunnel using Markov chain Monte Carlo(MCMC) simulation. The proposed approach is illustrated by a circular tunnel with an analytical solution, which was then applied to an experimental tunnel in Goupitan Hydropower Station, China. The mechanical properties and strength parameters of the surrounding rock mass were modeled as random variables. The displacement was predicted with the aid of the parameters updated by Bayesian inferences and agreed closely with monitored displacements. It indicates that Bayesian inferences combined the monitored data into the tunnel model to update its parameters dynamically. Further study indicated that the performance of Bayesian inferences is improved greatly by regularly supplementing field monitoring data. Bayesian inference is a significant and new approach for determining the mechanical parameters of the surrounding rock mass in a tunnel model and contributes to safe construction in rock engineering.     

An integrative approach to robust design and probabilistic risk assessment for CO<Subscript>2</Subscript> storage in geological formations

CO₂ storage in geological formations is currently being discussed intensively as a technology with a high potential for mitigating CO₂ emissions. However, any large-scale application requires a thorough analysis of the potential risks. Current numerical simulation models are too expensive for probabilistic risk analysis or stochastic approaches based on a brute-force approach of repeated simulation. Even single deterministic simulations may require parallel high-performance computing. The multiphase flow processes involved are too non-linear for quasi-linear error propagation and other simplified stochastic tools. As an alternative approach, we propose a massive stochastic model reduction based on the probabilistic collocation method. The model response is projected onto a higher-order orthogonal basis of polynomials to approximate dependence on uncertain parameters (porosity, permeability, etc.) and design parameters (injection rate, depth, etc.). This allows for a non-linear propagation of model uncertainty affecting the predicted risk, ensures fast computation, and provides a powerful tool for combining design variables and uncertain variables into one approach based on an integrative response surface. Thus, the design task of finding optimal injection regimes explicitly includes uncertainty, which leads to robust designs with a minimum failure probability. We validate our proposed stochastic approach by Monte Carlo simulation using a common 3D benchmark problem (Class et al., Comput Geosci 13:451–467, 2009). A reasonable compromise between computational efforts and precision was reached already with second-order polynomials. In our case study, the proposed approach yields a significant computational speed-up by a factor of 100 compared with the Monte Carlo evaluation. We demonstrate that, due to the non-linearity of the flow and transport processes during CO₂ injection, including uncertainty in the analysis leads to a systematic and significant shift of the predicted leakage rates toward higher values compared with deterministic simulations, affecting both risk estimates and the design of injection scenarios.     

土钉支护边坡动力模型的建立及地震响应分析

参照上部结构抗震设计的方法,对土钉支护边坡的抗震设计提出了3个水准的抗震设防要求和两阶段设计法。综合考虑输入地震波的特性、边坡土体特性、坡高因素,以及土体边坡在水平地震激励下,其位移响应主要为剪切型变形,建立了土钉支护边坡地震动力简化模型。该模型将同一水平层不稳定土体、稳定土体在水平地震作用下运动趋势假想为同步,两者之间相互作用很小,因此,取不稳定土体和面板为研究对象。由于面板很柔,其运动随着土体的运动而变化,其刚度只考虑剪切刚度,将土钉处理成弹性支座,建立了土钉支护边坡振动方程,求解了在简谐地震作用下土钉支护边坡动力响应的解析解,通过该模型可以得到滑移面附近土钉轴力动力响应和土钉支护边坡的弹性动位移响应。最后,结合案例进行了分析,比较了按此分析方法与用ADINA按弹性有限元方法计算的结果,两者吻合较好。故这种方法可用于土钉支护边坡的第1阶段的抗震分析与设计。     

An analytical solution for wave propagation in a square pile due to transient point load

This paper presents an analytical solution for wave propagation in a square pile due to transient point load. The differential equation of dynamic equilibrium is established considering propagation of waves in both vertical and transverse directions. The soil resistance is simulated by Voigt model. The three-dimensional analytical solution is deduced by using Fourier transform and the separation of variable method. The arithmetical results of the proposed solution show that the velocity responses along the radial direction at the pile top are highly non-uniform. In addition, Young’s modulus and the pile side length exert undisputable influences on the velocity responses.     

Effects of fractures on the well production in a coalbed methane reservoir

A water-gas two-phase flow model of a fractured vertical well in a coalbed methane (CBM) reservoir is developed based on the Buckley-Leverett equation. A program is developed based on the proposed model to predict the production of a CBM well. Compared with the monitoring results of an actual vertical fractured well, the accuracy and reliability of the analytical model are approved. The impact of reservoir parameters and fracture parameters on the CBM production are quantitatively studied. The results indicate that the production is affected by the three major factors: gas supply capacity of the matrix (viz., CBM volume released by unit volume of matrix during unit time), supply capacity of the fracture network to the major fracture, and gas exhausting ability of the major fracture to the wellbore. The permeability, matrix radius, and coal seam thickness with larger values are good for the gas supply capacity of the matrix. Higher flow conductivity of major fracture shows a stronger impact on the gas exhausting ability of the major fracture to the wellbore. A smaller fracture interval is beneficial to the supply capacity of the fracture network to the major fracture to obtain higher production. To improve the efficiency of the actual well, the flow conductivity of the major fracture should be no more than 200 mD m and the crack interval should be between 20 and 30 m. The achievements of this study provide a case reference for the production prediction of vertical fractured wells and the optimization design of hydraulic fracture parameters of CBM reservoirs.