基于有限元极限平衡法的尾矿坝坝体稳定分析

在平面应变条件下,采用基于弹塑性应力应变分析的有限元极限平衡法对尾矿坝坝体的稳定性进行研究。该方法可较好地结合极限平衡分析方法和有限元应力分析方法的特点。通过弹塑性有限元分析,合理准确地考虑结构整体应力场的影响,进而结合Hooke-Jeeves法优化搜索最危险滑动面的位置及其对应的最小安全系数。以某尾矿坝实际工程为例,比较分析了考虑正常、洪水以及特殊工况运行下有限元极限平衡法与极限平衡条分法和有限元强度折减法在安全系数大小及滑动面形状和位置的差异,重点运用有限元极限平衡法对尾矿坝坝体的稳定性进行验算并作出安全评价,为尾矿坝的安全设计施工提供可靠依据和技术支撑。     

Numerical modeling of dredging effect on berthing structure

Piles and diaphragm wall-supported berthing structure on marine soils are loaded laterally from horizontal soil movements generated by dredging. The literature on the adequacy of the finite element method modeling of berthing structure to analyze their behavior during dredging is limited. This paper describes a finite element approach for analyzing the lateral response of pile and diaphragm wall during dredging. Piles are represented by equivalent sheet-pile walls and a plane strain analysis using the finite element method is performed. Results from the finite element method are compared with full-scale field test data. Full-scale field test was conducted on a bearing structure to measure the lateral deflection on pile and diaphragm wall for their full length using inclinometer during dredging in sequence. The finite element method results are in good agreement with full-scale field results. Conclusions are drawn regarding application of the analytical method to study the effect of dredging on piles and diaphragm wall-supported berthing structures.     

有限元法在滑坡稳定性分析中的应用

对比了滑坡稳定性分析的极限平衡法、数值分析法、概率法。而后,介绍了有限元法滑坡稳定性分析基本原理,并应用有限元法对中里滑坡进行了稳定性分析,得出了中里滑坡的主应力图和稳定系数,并将稳定系数与其他方法的结果进行了比较,得出有限元法是一种可靠的分析方法,其结果可应用于后期治理设计。      

Analysis of subsidence caused by underground mining

Summary Some aspects of subsidence caused by longwall coal mining are analysed using the finite element method. Results of the analysis are compared with a true mine panel, where measurements on subsidence were available. Rock deformations in the overburden were modelled by using an elasto-plastic constitutive model. The study indicates that the shape of the subsidence profile can be predicted reasonably well by using nonlinear finite element analysis.     

岩溶区高速公路路堤及溶洞顶板稳定性数值分析

采用大型通用有限元软件ABAQUS对某室内模型试验进行了数值模拟,验证了三维有限元法进行洞室结构稳定性分析的可靠性。结合永武（永安-武平）高速公路A11合同段岩溶地质资料,将溶洞近似为与实际相近的椭球状,建立岩溶区路堤的三维有限元模型,分析溶洞顶板厚度、填方高度、溶洞位置以及溶洞充填状况对路堤及溶洞顶板稳定性的影响。为与现行顶板安全厚度确定方法进行比较,进行二维平面应变有限元分析。研究表明,对于较完整的顶板岩层,按顶板受拉破坏稳定安全系数为4确定的安全厚度过于保守,按梁板理论弯曲受拉破坏确定的安全厚度也偏于保守,是由于实际顶板岩层并不完全符合梁板理论的简化假设,而以有限元分析所得最大弯矩按梁板弯曲受拉破坏确定的安全厚度最为经济;实际工程中三维效应明显的溶洞,现行理论简化和二维数值分析法都较为粗糙,宜采用三维数值分析法。通过对上述问题的分析,为该高速公路岩溶区填方路段的设计和溶洞的处治与否提出指导性意见,已创造较好的经济效益和社会效益,可供类似工程应用时参考。     

离岸深水全直桩码头动力简化计算方法研究

离岸深水全直桩码头结构长期承受波浪力、撞击力等动荷载作用,仅从静力角度分析不能满足工程设计要求。为研究离岸深水全直桩码头结构动力简化计算方法,利用ABAQUS软件建立了全直桩码头的嵌固点模型和结构-地基相互作用的三维弹塑性有限元模型。基于所建立的有限元模型,采用模态分析法研究了上述码头结构的自振特性。分析表明,结构基本周期较长,与波浪力、撞击力荷载周期接近,结构动力响应突出。进一步计算发现,全直桩码头结构在波浪力、撞击力荷载作用下的动力响应有限元计算结果与采用单自由度系统位移动力放大系数理论公式计算结果吻合较好,由此提出全直桩码头结构动力响应可按单自由度系统位移动力放大系数理论公式进行简化计算。在此基础上,结合规范中的m法、p-y曲线法、NL法提出了上述码头结构动力简化计算方法,通过与动力有限元计算结果对比,验证了所建立简化计算方法的正确性。     

基于三维多重网格法的抗滑稳定计算精度分析

在以有限元法来分析评价坝基、边坡的稳定安全度时,多重网格法将有限元分析网格和滑移块体的滑面网格作为独立的两套网格,滑面网格可随意切割有限元单元。这样一套有限元网格和一套应力成果就可以进行任意的滑块稳定分析。着重分析了多重网格法的精度。针对三维柱体和锦屏高边坡,研究两套网格如何匹配以提高计算精度进行处理,初步得出了一些有关多重网格法计算精度的知识。最后,应用基于多重网格法的第2类方法分析了锦屏左坝肩高边坡的稳定性问题。     

单桩承载力的一种直接动测法

应用高应变动测法确定基桩的承载力非常依赖测试人员的经验且存在多解性问题。为克服这种现象,介绍了一种新的直接动测法,并以一灌注桩的现场静载试验为对象,建立了详细的基桩有限元静力、动力分析模型。静力模拟比较分析表明,该有限元模型能准确地模拟文献的结果。在此基础上,通过有限元动力分析,得到锤击下桩的应力场、位移场、速度场和加速度场。根据有限元计算结果,分析了直接动测法的有效性。直接动测法确定的承载力和静载法一致。该直接动测法消除了现有高应变法存在的多解性问题,故可以得到单桩承载力的客观性结果。     

Computer modeling of folding in rocks

The results of numerical mathematical modeling of rock deformations under compression are given. The numerical solutions are obtained using discretization of the equations of the solid mechanics with the finite element method. The model of an ideal elastic-plastic material with a Huber-Mises yield surface was used in the calculations. The layered medium structure is taken into account in modeling of the compression of layers of the lower/middle crust on a local scale. The natural folds in strongly deformed metamorphic sequences were reproduced by the mathematical deformation models. It is shown that folding in the lower part of the Earth’s crust is possible, when the yield stresses of the host rocks are approximately two orders of magnitude less than those of the hard layers. The effect of the boundary conditions and thickness of the compressed rocks on the folding is shown.     

地震波场正演模拟中的有限单元方法

本文介绍了地震波场正演模拟中的有限单元方法。首先从粘弹性波动理论出发,推导了满足粘弹性介质的有限元方程,其次用不含阻尼矩阵的有限元方程模拟了地震记录。地表实例是自由界面的瑞雷面波。VSP记录有水平层、倾斜层;有零井源距和非零井源距;有线状和对称震源。 在程序中采用了超压缩存贮和分块计算措施,对降低有限元法要求计算机内存和运算速度的苛刻条件上做了有益的尝试。     

基于地应力边界条件的复杂边坡抗滑桩支护设计研究

复杂边坡稳定性问题一直是岩土工程界重点研究的问题之一,提高复杂边坡安全系数的方法有很多,其中采用抗滑桩支护是常用手段之一。以安徽省某高速公路复杂边坡为研究对象,采用套孔应力解除法对该边坡岩土体进行了地应力测试,以测试结果作为有限元分析的边界条件,分别对该边坡加固前和采用抗滑桩加固后两种工况进行数值分析,得到有代表性的有限元分析结果;并对有限元的计算出位移结果和实测位移进行了对比。分析表明,未采用抗滑桩支护的边坡,会产生局部应力过分集中、整体变形和水平向位移较大、塑性区贯通等不利于边坡稳定的危险状况。采用了抗滑桩支护后,贯通的塑性区消失,边坡变形和位移减少,应力分布趋于均匀,安全系数大幅度提高。有限元计算出的位移结果和实测位移十分接近,匹合度较好,从而说明了有限元分析手段的可靠性。     

Elasto-plastic analysis of an underground opening by FEM and coupled FEBEM

The elasto-plastic analysis of a circular underground opening using the finite element method (FEM), the coupled finite element and the boundary element method (FEBEM) is presented. The coupling procedure and elasto-plastic formulation are discussed. The effect of in situ stress ratio on yielded zone, displaced shape and principal stresses is presented. The results of FEBEM analysis are compared with those obtained from FEM. The computation time and the number of iterations as required by FEBEM and FEM are compared.     

非稳定温度场的三维p型有限元方法研究

非稳定温度场的边界往往是温度梯度较大的区域,需要划分较密的网格。但在三维温度应力的计算中太密的网格往往使得计算成本很高。将p型阶谱有限元方法引入温度场的计算,则可以开拓一条解决这类问题的有效途径。研究了p型有限元中的温度边界条件和初值,编制了有关程序,对同一个算例,采用不同网格密度的常规有限元法和p型有限元法进行了温度场计算分析和对比。研究结果表明,用p型有限元进行非稳定温度场计算,在网格稀疏的条件下能够达到较高的精度。     

滑坡稳定性分析与安全系数取值研究

在地质灾害中,滑坡灾害分布最广、发生频率最高、危害最大,是我国地质灾害防治的主要对象。滑坡研究的主要任务之一是稳定性分析。论文针对三板溪水电站进水口滑坡和水布垭水电站大岩淌滑坡,分别采用刚体极限平衡法、块体单元法、有限单元法对稳定安全系数和条间推力分布曲线进行分析和比较,研究各种方法的差异和计算精度,以及各种方法稳定安全系数的取值标准,希望为水利水电工程边坡设计规范有关有限单元法、块体单元法条目的编写和相应的允许安全系数取值标准的制定提供一定的参考。研究结果表明:①块体单元法计算的安全系数与刚体极限平衡法计算的安全系数可以采用相同的取值标准,而有限单元法计算的安全系数取值标准可适当降低;②块体单元法计算的条间水平推力、条间剪切力与刚体极限平衡法(尤其是Sarma法)计算的条间水平推力、条间剪切力非常相近,有限单元法计算的条间水平推力、条间剪切力较小;③分析方法、力学模型与参数取值应与安全系数取值标准相匹配。     

虚拟单元法求解渗流流量

提出了计算渗流场流量的虚拟单元法,推导了虚拟单元法有限元方程,并探讨了虚拟单元法的有限元实现方法,建立了虚拟单元法求解渗流流量的理论和方法体系。二维径向渗流和平底坝基算例比较分析表明,虚拟单元法具有计算精度高及有限元实现简单并容易编写通用程序等优点。虚拟单元法与中线法计算流量具有相同的精度,但在应用上虚拟单元法有一定优势。     

Finite element modeling of wave propagation problems in multilayered soils resting on a rigid base

A new finite element scheme is proposed, in this paper, for solving two-dimensional wave propagation problems in multilayered soils resting on a rigid base. The multilayered soils are treated as multiple horizontal layers of lateral infinite extension in geometry. Since these horizontal layers can be truncated by two artificially truncated vertical boundaries, two high-order artificial boundary conditions are applied for propagating the incoming waves from the interior domain into the far field of the system. Both the semi-analytical method and the truncated boundary migration procedure are used to derive the high-order artificial boundary conditions, which are comprised of a physically meaningful dashpot and a generalized energy absorber. The main advantage of using the proposed finite element scheme is that the derived artificial boundary condition can be straightforwardly implemented in the finite element analysis, without violating the band/sparse structure of the conventional finite element equation. The related numerical examples have demonstrated that the proposed finite element scheme is of high accuracy in dealing with wave propagation problems in multiple horizontal layers.     

Finite element analysis of underground openings

A highly efficient technique is presented for the finite element analysis of stresses around underground openings surrounded by an infinite extent of linearly elastic rock mass. The effect of unbounded rock is incorporated in the analysis by using elastic supports around the finite element model. Expressions for stiffnesses of the elastic supports are derived and these are found to depend on the location of elastic supports, the shear modulus and Poisson's ratio of the rock, and the ratio of horizontal to vertical initial stresses. With the use of the proposed technique, the extent of the finite domain to be considered in the analysis is highly reduced. This results in a great computational advantage. The other principal merit of the proposed technique is that a standard finite element code for stress analysis may be used without making any modification. Results of some numerical tests are reported to demonstrate the effectiveness and efficiency of the proposed method. The technique has the potential of being applied to more complex problems of unbounded domains in geomechanics.     

软土地基上路堤填筑的破坏性状分析

在软土地区建设高速公路、铁路等的填筑路堤越来越多,由于软土的不排水抗剪强度较低,软土地区路堤的填筑十分困难,而软土地基的渗透性很低,通常认为路堤的填筑是在不排水条件下进行。经典的设计方法是基于极限平衡法,无法考虑变形等因素,因此,能综合考虑稳定及变形等因素的有限元方法成为分析此类问题的一种有效手段。依据连云港铁路路堤在软土地基上的填筑试验资料,采用大变形有限元方法及Mohr-Coulomb模型,采用现场试验得到的土体强度值,对平面应变条件下路堤填筑的破坏性状及极限填筑高度等进行了分析,并对比了有限元分析结果及现场试验结果,结果表明,采用大变形有限元方法能够较准确的得到填筑路堤的极限高度,并有效地分析填筑路堤的破坏性状。     

增量位移反分析在水电地下洞室工程中的应用

反分析是确定计算模型参数的有效方法.通常多采用量测所得全量位移进行反演计算.但地下工程中许多量测数据为增量位移,且实际施工过程可以通过建立分步开挖的有限元模型来模拟.据此,结合某水电站地下洞室工程中地下厂房的开挖,建立了模拟动态施工的有限元模型,利用某一开挖步施工前后量测值之差,采用增量位移优化反分析方法对洞室附近初始地应力场及围岩弹性模量进行了反演.计算所得增量位移与实测值符合较好,表明了这种方法的可行性.同时,根据分析结果对该方法进行了评价.     

An improved search strategy for the critical slip surface using finite element stress fields

The finite element method can be used to advantage in slope stability problems. This paper proposes a technique to search for the critical slip surface as well as to define and calculate the factor of safety for the slope, when the finite element method is used to model its formation. First, stresses are estimated at each Gaussian point from the finite element analysis. Then, the global stress smoothing method is applied to get a continuous stress field. Based on this stress field, the factor of safety is calculated for a specified slip surface by a stress integration scheme. An improved search strategy is proposed for a noncircular critical surface which starts with a search method for a circular critical surface. During the search process, points defining a trial slip surface can freely move in the finite element mesh subject to some kinematical constraints. This method can be applied to both the limit equilibrium method and the finite element method. Effects of the slope stress history and soil parameters on the resulting critical surface are investigated.