Three-dimensional upper bound limit analysis of underground cavities using nonlinear Baker failure criterion

A generalized nonlinear Baker failure criterion is employed with the upper bound limit analysis to study the surrounding rock stability of underground cavities. A three-dimensional(3D) failure mode is established by extending the two-dimensional(2D) failure mode, which offers an upper bound expression of the surrounding rock pressure. This method is validated with a series of examples before the influence of four parameters of scale parameter, curvature parameter, shift parameter and lateral pressure coefficient, on the surrounding rock pressure is analyzed. According to these results, failure ranges of the underground cavities are determined. The following conclusions are reached:(1) the proposed approach is more accurate to predict surrounding rock pressure than the Mohr-Coulomb failure criterion;(2) the surrounding rock with large scale parameter, curvature parameter, shift parameter, and lateral pressure coefficient can lead to a more stable underground cavity;(3) the failure range in 3D mode can be predicted according to the upper bound solutions.

Upper-bound limit analysis based on the natural element method

The natural element method （NEM） is a newly- developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functions in the framework of Galerkin method. Owing to its distinctive advantages, the NEM is used widely in many problems of computational mechanics. Utilizing the NEM, this paper deals with numerical limit analysis of structures made up of perfectly rigid-plastic material. According to kinematic the- orem of plastic limit analysis, a mathematical programming natural element formulation is established for determining the upper bound multiplier of plane problems, and a direct iteration algorithm is proposed accordingly to solve it. In this algorithm, the plastic incompressibility condition is handled by two different treatments, and the nonlinearity and nons- moothness of the goal function are overcome by distinguishing the rigid zones from the plastic zones at each iteration. The procedure implementation of iterative process is quite simple and effective because each iteration is equivalent to solving an associated elastic problem. The obtained limit load multiplier is proved to monotonically converge to the upper bound of true solution. Several benchmark examples are investigated to validate the significant performance of the NEM in the application field of limit analysis.

LOWER BOUND LIMIT ANALYSIS OF THREE-DIMENSIONAL ELASTOPLASTIC STRUCTURES BY BOUNDARY ELEMENT METHOD

Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.

Upper bound limit analysis of roof collapse in shallow tunnels with arbitrary cross sections under condition of seepage force

The analytical solutions for predicting the exact shape of collapse mechanisms in shallow tunnels with arbitrary excavation profiles were obtained by virtue of the upper bound theorem of limit analysis and variation principle according to Hoek-Brown failure criterion. The seepage force was included in the upper bound limit analysis, and it was computed from the gradient of excess pore pressure distribution. The seepage was regarded as a work rate of external force. The numerical results of roof collapse in square and circular tunnels with different rock parameters were derived and discussed, which proves to be valid in comparison with the previous work. The influences of different parameters on the shape of collapsing blocks were also discussed.

Combined influence of nonlinearity and dilation on slope stability evaluated by upper-bound limit analysis

The combined influence of nonlinearity and dilation on slope stability was evaluated using the upper-bound limit analysis theorem.The mechanism of slope collapse was analyzed by dividing it into arbitrary discrete soil blocks with the nonlinear Mohr–Coulomb failure criterion and nonassociated flow rule.The multipoint tangent(multi-tangent) technique was used to analyze the slope stability by linearizing the nonlinear failure criterion.A general expression for the slope safety factor was derived based on the virtual work principle and the strength reduction technique,and the global slope safety factor can be obtained by the optimization method of nonlinear sequential quadratic programming.The results show better agreement with previous research result when the nonlinear failure criterion reduces to a linear failure criterion or the non-associated flow rule reduces to an associated flow rule,which demonstrates the rationality of the presented method.Slope safety factors calculated by the multi-tangent inclined-slices technique were smaller than those obtained by the traditional single-tangent inclined-slices technique.The results show that the multi-tangent inclined-slices technique is a safe and effective method of slope stability limit analysis.The combined effect of nonlinearity and dilation on slope stability was analyzed,and the parameter analysis indicates that nonlinearity and dilation have significant influence on the result of slope stability analysis.

Comparative study of material point method and upper bound limit analysis in slope stability analysis

The upper bound limit analysis(UBLA)is one of the key research directions in geotechnical engineering and is widely used in engineering practice.UBLA assumes that the slip surface with the minimum factor of safety(FSmin)is the critical slip surface,and then applies it to slope stability analysis.However,the hypothesis of UBLA has not been systematically verified,which may be due to the fact that the traditional numerical method is difficult to simulate the large deformation.In this study,in order to systematically verify the assumption of UBLA,material point method(MPM),which is suitable to simulate the large deformation of continuous media,is used to simulate the whole process of the slope failure,including the large-scale transportation and deposition of soil mass after slope failure.And a series of comparative studies are conducted on the stability of cohesive slopes using UBLA and MPM.The proposed study indicated that the slope angle,internal friction angle and cohesion have a remarkable effect on the slip surface of the cohesive slope.Also,for stable slopes,the calculation results of the two are relatively close.However,for unstable slopes,the slider volume determined by the UBLA is much smaller than the slider volume determined by the MPM.In other words,for unstable slopes,the critical slip surface of UBLA is very different from the slip surface when the slope failure occurs,and when the UBLA is applied to the stability analysis of unstable slope,it will lead to extremely unfavorable results.

SOME PROPERTIES OF BOUNDED SETS IN INDUCTIVE LIMITS

Let (E,ξ）=indlim (En,ξn） be an inductive limit of a sequence of locally convex spaces,For brevity,denote by (DS) each set Bbounded in (E,ξ） is contained in some En； and (DST) each set B bounded in (E,ξ） is contained and bounded in some (En,ξn）. Theovem 1.(DS) holds provided that (i) for each n∈N,there is a neighborhood Un of o in (En,ξn） and m(n)∈ such that -↑Un^E包含于Em(n),and (ii) for any neighborhood V n of o in (En,ξn）,∞↑Un=1 Vn absorbs every bounded set in (E,ξ）. theorem 2 Let all (En,ξn） be metrizable and (DS) hold,then for each bounded set B IN (E,ξ）and each n ∈N thcrc is a neighborhood U k of o in (Ek,ξk）, 1≤k≤n ,and m(n)∈N such that ——↑(B+U1+U2+…+Un)^E包含于 Em(n). theorem 3. Let all (En,ξn） be Frechet spaces.Then (DST) holds if and only if (i) for each n ∈N,there is u neighborhood U n of in (En,ξn） and m(n)∈N such that 0↑Un^E包含于Em(n),and (ii) for each each closed ,absosed,absolutely conuex,bounded set B in (E,ξ）,∞↑Un=1(（εnB)∩Un）absorbs B,where U n is any neighborhood of o in (En,ξn） and εn is any positive number for every n ∈N。

INCOMPRESSIBLE LIMIT OF THE NON-ISENTROPIC MAGNETOHYDRODYNAMIC EQUATIONS IN BOUNDED DOMAINS

The incompressible limit of the non-isentropic magnetohydrodynamic equations with zero thermal coefficient, in a two dimensional bounded domain with the Dirichlet condi- tion for velocity and perfectly conducting boundary condition for magnetic field, is rigorously justified.

A High Order Accurate Bound-Preserving Compact Finite Difference Scheme for Two-Dimensional Incompressible Flow

For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.

Upper bound analysis of ultimate pullout capacity of shallow 3-D circular plate anchors based on nonlinear Mohr-Coulomb failure criterion

Based on the nonlinear Mohr-Coulomb failure criterion and the associated flow rules,the three-dimensional(3-D)axisymmetric failure mechanism of shallow horizontal circular plate anchors that are subjected to the ultimate pullout capacity(UPC)is determined.A derivative function of the projection function for projecting the 3-D axisymmetric failure surface on plane is deduced using the variation theory.By using difference principle,the primitive function of failure surface satisfying boundary condition and numerical solution to its corresponding ultimate pullout capacity function are obtained.The influences of nonlinear Mohr-Coulomb parameters on UPC and failure mechanism are studied.The result shows that UPC decreases with dimensionless parameter m and uniaxial tensile strength increases but increases when depth and radius of plate anchor,surface overload,initial cohesion,geomaterial density and friction angle increase.The failure surface is similar to a symmetrical spatial funnel,and its shape is mainly determined by dimensionless parameter m;the surface damage range expands with the increase of radius and depth of the plate anchor as well as initial cohesion but decreases with the increase of dimensionless parameter m and uniaxial tensile strength as well as geomaterial density.As the dimensionless parameter m=2.0,the numerical solution of UPC based on the difference principle is proved to be feasible and effective through the comparison with the exact solution.In addition,the comparison between solutions of UPC computed by variation method and those computed by upper bound method indicate that variation method outperforms upper bound method.

Limit variation analysis of shallow rectangular tunnels collapsing with double-layer rock mass based on a three-dimensional failure mechanism

In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) velocity discontinuity surfaces. According to the virtual work principle, the difference theorem and the variation method, the collapse surface of double-layer rock mass is determined based on the Hoek-Brown failure criterion. The formula can be degenerated to a single-layer rock collapsing problem when the rock mass is homogeneous. To estimate the validity of the result, the numerical simulation software PLAXIS 3D is used to simulate the collapse of shallow tunnels with double-layer rock mass, and the comparative analysis shows that numerical results are in good agreement with upper-bound solutions. According to the results of parametric analysis, the potential range of collapse of a double-layer rock mass above a shallow cavity decreases with a decrease in A1/A2,σci1/σci2 and σtm1/σtm2 and an increase in B1/B2,γ1/γ2. The range will decrease with a decrease in support pressure q and increase with a decrease in surface overload σs. Therefore, reinforced supporting is beneficial to improve the stability of the cavity during actual construction.

极限分割估计矩形的Bounding-Box定位算法

针对无线传感器网络中Bounding-Box算法定位精度低的问题,在Bounding-Box算法的基础上提出一种极限分割估计矩形的方法来改进定位算法,分割后的估计矩形产生一个待选质心,当满足分割条件时,对估计矩形继续分割,并不断产生待选质心,当满足终止分割条件时,将上一次分割得到的待选质心坐标作为未知节点的最终位置,通过极限分割的方法可以修正未知节点的定位误差。仿真结果表明,在无需增加额外通信开销的情况下,改进的算法在一定程度上降低了算法的平均相对定位误差。

LOW MACH NUMBER LIMIT OF A COMPRESSIBLE NON-ISOTHERMAL NEMATIC LIQUID CRYSTALS MODEL

In this paper, we study the low Mach number limit of a compressible nonisothermal model for nematic liquid crystals in a bounded domain. We establish the uniform estimates with respect to the Mach number, and thus prove the convergence to the solution of the incompressible model for nematic liquid crystals.

LIMITING WEAK-TYPE BEHAVIORS FOR CERTAIN LITTLEWOOD-PALEY FUNCTIONS

In this paper, we establish the following limiting weak-type behaviors of Littlewood-Paley g-function g_φ: for nonnegative function f∈ L^1(R^n),■and ■where f_t(x) =t^(-n)f(t^(-1) x) for t > 0. Meanwhile, the corresponding results for Marcinkiewicz integral and its fractional version with kernels satisfying L_α~q-Dini condition are also given.

大截面矩形顶管施工对既有管线影响研究

顶管施工会扰动邻近建筑物,危及工程安全。该文以苏州市城北路综合管廊元和塘段顶管工程为背景,通过ABAQUS有限元分析软件建立三维模型,研究顶管下穿地下管线施工对管周土及既有管线影响,并基于极限分析方法探讨了管线失稳模式。结果表明:大直径矩形顶管顶进方向的影响范围约为2倍顶管宽度,两侧的影响范围约为1倍顶管宽度,地下管线到顶管轴线的水平距离为0.5倍~1倍顶管宽度时,应注意管线受到的水平方向附加荷载是否会导致管线破坏;管线与顶管相对位置对管线破坏影响较为显著,越靠近顶管的管线破坏越严重,滑移线由顶管两侧经由管线向地面发展,破坏部位由管线顶部两侧延伸至底部;除受与顶管相对位置影响外,管线破坏模式还受矩形顶管形状影响,管线顶部及靠近顶管顶角区域破坏风险最大。

考虑水力效应的裂缝边坡稳定性非线性能耗分析

现有考虑水力效应的裂缝边坡稳定性分析大多基于线性破坏准则,而岩土体破坏往往呈现非线性特征,因此开展水力效应影响下的裂缝边坡稳定性非线性极限上限分析具有重要意义。基于极限分析上限定理及强度折减技术,结合“外切线法”非线性破坏准则,构建坡顶含竖直裂缝的边坡对数螺旋线破坏模式,根据虚功原理推导出裂缝边坡安全系数解析式,通过MATLAB优化计算,结合边坡工程实例探讨了典型因素对裂缝边坡稳定性、临界裂缝及滑动面位置的影响规律。研究结果表明:随着地下水位h的持续上升,边坡安全系数不断降低,临界裂缝深度逐渐增大且临界裂缝位置逐渐向坡顶缘偏移;非线性系数m明显影响边坡的稳定性,边坡安全系数随着非线性系数的增大显著减小,采用线性破坏准则会高估地下水位变化对边坡稳定性的影响;随着非线性系数的增大,临界裂缝深度随之增大,临界裂缝位置距离坡顶缘越来越远,滑坡体体积逐渐增大;裂缝与坡顶缘距离l m随着非线性系数m的增大而增大(当m增大0.2时,l m增大约1 m),且随着裂缝与坡顶缘距离的增大,边坡安全系数无明显变化,临界裂缝深度先逐渐减小后趋于稳定。

开洞刚性墙后土体破坏模式运动单元上限分析

超载作用下开洞地下刚性连续墙后土体稳定性问题具有广泛的工程应用背景。考虑墙体矩形开洞宽度足够大,且刚性墙与土体完全光滑接触等条件,参考构造刚性多块体上限法思路,预设墙后土体极简破坏机构并转化为初始三角网格,提出交替循环式新增节点、单元缩减、对角互换及成层加密等网格更新措施,同步开展基于刚体平动运动单元上限方法的刚性墙后土体破坏分析,通过多层级非线性规划运算获得高精度临界超载系数σ_(s)/c上限解,界定墙后土体滑移线网破坏模式及其定量范围与扩展边界。进一步依前述流程计算获得多参数影响下的临界超载系数σ_(s)/c上限解图表,结果表明:1)刚性墙后土体典型破坏特征表现为土体上方和开洞处的整体下沉刚性区域,开洞前方核心区域以密集滑移线网区隔的刚性滑块群,以及外围主要滑动面;2)临界超载系数σ_(s)/c随土体内摩擦角ϕ增大而增大,且后续增幅显著增加;3)σ_(s)/c随开洞深度比C/D的增大呈现出近似线性增大的趋势;4)相较于ϕ和C/D,土体重度参数γD/c的变化对于σ_(s)/c的影响不甚明显,σ_(s)/c随γD/c增大呈现近似线性减小的趋势;5)σ_(s)/c受剪胀角参数ψ/ϕ的影响与ϕ的增大程度密切相关。本研究结果可为开洞地下刚性连续墙后土体失稳破坏模式定量分析、破坏机理研究及工程措施制定等工作提供有力的数据支持。

极限下限分析的区域光滑径向点插值法

极限分析的高效数值计算方法在结构设计和安全评定中具有非常重要的作用.为了更加有效地求解极限分析问题,将区域光滑径向点插值法与二阶锥规划相结合,提出了理想弹塑性结构极限下限分析的一种新方法.将问题域离散为简单的三角形背景单元,每个单元进一步划分成若干个光滑域.为了将复杂的域积分转化为简单的边界积分,并且避免计算形函数的导数,采用广义梯度光滑技术对每个光滑域进行应变光滑处理.由于径向点插值法构造的形函数满足Kronecker delta性质,本质边界条件可以直接施加.依据下限定理,在满足以等效积分弱形式表达的自平衡应力场平衡条件的基础上,用二阶锥规划成功构建了极限分析下限法的计算模型,从而可方便地通过基于原始-对偶内点法的数学规划求解器MOSEK直接求解该问题.数值算例结果表明,本文所提方法有效地克服了维数障碍问题,具有较高的计算精度,并且计算结果对网格畸变十分不敏感.

考虑基质吸力影响的非饱和路堤三维稳定性上限分析

采用极限上限分析定理构建了路堤失稳的三维旋转破坏机构,继而根据Bishop非饱和土抗剪强度理论,考虑非饱路堤内部基质吸力的空间分布及其随地下水位的变化,建立了路堤破坏土体外力功率与内部能量耗散功率的能量守恒方程,并利用遗传算法编写了非饱和路堤最小上限解的高效搜索算法。通过将非饱和土路堤基底破坏模式退化为边坡破坏模式,并与现有非饱和土边坡稳定性计算结果对比,验证了所提上限解的正确性和遗传搜寻算法的准确性。进一步地,对路堤三维稳定性的关键影响因素展开了系统分析,研究了路堤填土孔径分布、基质吸力、进气值倒数、路堤倾角以及有效内摩擦角等因素对路堤三维稳定性的影响规律。研究表明,非饱和土质路堤的稳定性不仅取决于路堤填土性质,而且依赖于影响土体吸力大小与分布的填土孔径参数和进气值等因素。地下水位升降引起的基质吸力变化对路堤稳定性存在显著影响。研究结果为路堤稳定性精细化分析提供了重要的理论依据。

海啸作用下滨水挡土墙抗震滑动稳定性上限分析

针对沿海地区滨水挡土墙在海啸和地震联合作用下易发生被动破坏的问题,结合极限分析上限法和拟静力法确定地震荷载,研究了海啸作用下滨水挡土墙的抗震滑动稳定性。通过建立墙体滑动与破坏区土楔体整体滑动的墙-土系统及速度容许场,推导出水平地震加速度系数(kh)的表达式。分析了海啸、海平面与地下水高度及墙-地和墙-土摩擦角对地震屈服加速度系数的影响。结果表明,控制地下水位在适当高度并增大墙-土摩擦角可以有效提高滨水挡土墙的抗震滑动稳定性。与经典极限平衡理论相比,该方法考虑了土体的塑性流动及挡土墙的位移模式,且无需计算地震被动土压力的值。