白鹤滩水电站拱坝推力对坝肩的影响分析

对于地质条件复杂的高拱坝,坝肩的变形及稳定性对拱坝的受力状态有很大的影响,在研究过程中必须考虑拱坝与坝肩之间的相互作用。为降低三维有限元建模的难度和更好地模拟坝肩复杂的地质构造,分别建立坝肩和拱坝的有限元模型,采用试算迭代的方法考虑二者之间的相互作用;并运用卸荷岩体力学理论的方法分析坝肩岩体开挖卸荷后,在拱坝推力作用下的应力应变情况。结果表明,考虑拱坝推力后坝肩槽上、下游边坡位移变化不大,变幅在1 cm左右;拱坝推力对拱端边坡有比较明显的影响,在推力作用下,边坡向坡内有一定的变形。     

一种基于安全因子判据的坝肩裂隙岩体动力抗滑稳定分析方法

基于某混凝土拱坝,结合现行拱坝设计规范,研究提出了一种基于安全因子判据的坝肩岩体动力抗滑稳定安全度分析方法。该方法综合利用有限元和（刚体）极限平衡两种方法的优点,考虑了凝聚力和摩擦力是不同的抗力,其不确定性具有较大差异性,定义了抗滑稳定安全因子,据此来判别坝肩岩体的抗滑稳定安全度,安全因子随着凝聚力与摩擦力所占权重不同而变化,这种方法既考虑了坝肩岩体与坝体的联动作用,又能够反映拱坝体系真实的工作性态。通过某混凝土高拱坝工程实例,分析了坝肩裂隙岩体的抗震稳定性,结果表明,该方法具有较高的实用价值。     

基于超载法的拱坝坝肩地震动力稳定性方法研究

拱坝坝肩地震动力稳定性一直是工程界关注的一大热点,然而对坝肩动力稳定性的评判方法和标准的研究还不成熟。为完善坝肩动力稳定分析方法,将有限元法与刚体极限平衡法的优势相结合,提出了地震动力超载方法以及极限累积位移评判标准。首先,基于有限元法开展实时动力分析,对坝肩可能的滑移块体采用多重网格应力插值求任意时刻的阻滑力和下滑力,确定其主滑方向,而后采用刚体极限平衡法计算每个滑块的动安全系数;同时,对瞬时动安全系数小于1.0的时段,进行加速度积分,求出滑块累积失稳位移。结合滑块尺寸及滑动面强度参数提出极限累积位移计算公式作为失稳判据,通过比较累积位移与极限累积位移判定滑块的稳定性。基于以上思路开发了一套坝肩动力稳定分析系统SAFEDAM。以沙牌拱坝经受反映5.12地震的人工波为例,开展了沙牌拱坝坝肩动力超载稳定分析,验证了极限累积位移判据的合理性。研究表明,沙牌左岸坝肩超载稳定安全倍数Kp0 =3.2,右岸坝肩超载稳定安全倍数Kp0 =4.2。工程实例计算表明,该方法计算坝肩动力稳定性使用较为方便可靠,具有较大的工程实用价值,为坝肩动力稳定性分析和评价提供了一条新路径。     

Application of Numerical Method to Soil Mechanics

This is a summing-up report based on 48 papers on the applications of numerical method to soil mechanics. Recently, there is a new development of coupling of finite element with other numerical ones in soil mechanics. The dynamic numerical analyses have found further application to nonlinear stress-strain relationship, effective stress, and three dimensional problems. The suitable contact element or conjunct element could be used for contact conditions in finite element calculation. It is an important direction to combine usual numerical method with probabilistic method, which has been used in evaluations of consolidation effect, slope stability, dam, seepage, deformation calculation, interaction, and the performance analyses of pole and foundation.     

数值计算方法在土力学中的应用

本文是对48篇有关数值方法在土力学方面应用的论文进行综述的报告。近年来,在土力学方面有限元和其它数值方法的耦合有了新的发展,动力分析在非线性应力应变关系、三维问题和有效应力分析中得到了进一步的应用,接触单元或连接单元可用来处理有限元计算时的接触条件。数值方法和概率方法的结合是一个值得重视的发展方向。这些方法可应用于估计固结效应、土坡稳定、土石坝、渗流、变形计算、共同作用、桩和地基性状分析等方面。     

拱坝联合坝肩岩体三维计算模型构建

坝址区复杂的地质条件与岩体力学环境变化,直接影响坝肩岩体抗滑稳定及拱坝的安全。我国西南乌东德超高拱坝坝址区工程地质条件极为复杂。坝体载荷后,区域内地层岩性、断层构造、岩溶系统、拱座内不连续裂隙等自然地质条件,成为影响坝肩岩体稳定性的主要因素。基于可研阶段地质勘探数据、二维地质剖面和工程设计数据,采用面向对象技术、多种软件协同作业,建立自然地质条件和工程对象的几何模型。构建适用于地应力反演分析,坝体载荷、工程开挖等因素影响下坝肩稳定性分析的三维精细有限元计算模型。该三维模型最大程度地模拟了工程岩体复杂结构,为高拱坝坝肩岩体稳定性综合分析提供良好的基础。计算结果表明:河谷两侧近地表地应力带及河床深部地应力区应力水平与实际量测值接近。坝肩槽开挖完成后,右岸拱肩槽800 m高程附近位移值最大达到12.9 mm。两岸9个可能滑动楔块在不同荷载组合条件下有一定差异,整体上左岸坝肩稳定性强于右岸。组合地震荷载作用下,所有楔块安全系数均大幅下降,但仍具有一定的安全裕度。     

高拱坝施工初-中期导流风险模型及应用

面向工程设计阶段,采用高拱坝施工动态仿真技术获取施工初-中期挡水度汛面貌数据,综合考虑水文、水力随机性因素,构建高拱坝施工初-中期导流风险模型,提出采用Monte Carlo方法耦合挡水度汛面貌数据和主要随机因素进行风险模型求解的方法。基于风险分析原理提出了导流洞设计的风险判别方法,给出导流洞尺寸设计优化的数学模型和具体步骤。通过金沙江上游某高拱坝工程案例分析的结果表明:所提风险模型及求解方法是适用的、有效的,该模型能够得到整个施工初-中期导流风险率,较为客观地反映高拱坝施工中期度汛可能存在的两种挡水情况,克服了初期导流风险模型的局限性;施工中期导流风险率随导流洞尺寸增大而减小,导流洞尺寸设计的可行方案集存在界限,即优化方案。研究成果可为高拱坝施工导流的风险决策和设计优化提供理论支撑。     

A dual mortar contact method for porous media and its application to clay-core rockfill dams

The clay-core rockfill dam is a multibody contact system in which the hydromechanical response of the clay core plays a crucial role. This complex problem is highly challenging to model numerically. We present a numerical approach that considers the multibody contact, consolidation, and strong geometric and material nonlinearities for the modeling of clay-core rockfill dams. Within the framework of the dual mortar finite element method, the presented approach considers the contact bodies as independent porous media continuums. The nonlinear contact conditions are derived based on the effective contact traction on contact interfaces and pore pressure continuity. The weak forms are obtained by introducing Lagrange multipliers as additional unknowns, which are then condensed through an extended general transformation. The presented method is first validated with a patch test considering the contact between two porous media. Then, a three-dimensional analysis of the Rumei clay-core rockfill dam is performed. The main numerical analysis concerns are the two observation galleries planned for construction inside the clay core. The galleries consist of dozens of tunnel-like concrete blocks, giving rise to complex concrete-concrete and concrete-clay contacts. The discontinuous separation and sliding between concrete blocks are investigated. For the concrete-concrete contact, both hard and soft joint approaches are evaluated and compared. The pore pressure results of the concrete structures are also analyzed.     

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

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

Dynamics of pile driving by the finite element method

It is proposed that the dynamics of pile driving based on the one-dimensional theory of wave propagation be analyzed by the finite element method. This approach enables one to continuously interpolate the displacement, velocity and acceleration profiles throughout the pile length. In contrast to the finite difference technique presented by E. A. L. Smith [1] and the finite element procedure presented by I. M. Smith [2] in which initial conditions are defined on the basis of a prescribed hammer velocity, the method presented herein defines initial conditions on the basis of a prescribed impact force versus time curve at the pile/hammer point of contact. Applications of the proposed technique to typical pile driving problems on an elastoplastic soil and using an implicit time-integration scheme are discussed using a numerical example.     

碾压混凝土高拱坝坝肩稳定及坝体开裂静动力分析

强震区高拱坝抗震问题的研究十分重要,高拱坝静力稳定研究和极限抗震能力复核,是确保地震后不发生库水失控下泄的关键课题之一。采用物理模型试验与数值计算的方法,根据汶川震后新核定参数,并考虑降强因素的影响,对沙牌碾压混凝土拱坝进行静动力分析和复核。通过三维地质力学模型综合法静力破坏试验,研究拱坝坝肩在正常荷载作用下的稳定安全性,得到沙牌拱坝坝肩综合稳定安全系数为3.76,建议对两坝肩开裂较严重区域进行加固处理。在此基础上,基于反应谱理论进行三维有限元动力分析,将应力成果导入开发程序由Drucker-Prager（D-P）和Mohr-Coulomb（M-C）准则对加固后的沙牌拱坝进行动力复核计算,并研究大坝的开裂情况和极限抗震能力,计算结果表明,沙牌拱坝的整体抗震稳定性能良好,仅在万年一遇地震工况下发生坝肩浅表层失稳。研究成果表明,沙牌拱坝坝肩稳定性较高、坝肩加固效果良好,从而在“5.12”汶川地震中表现出超强的整体稳定性和抗震能力,其全面开展的研究工作和采取的工程措施可供国内强震区同类型拱坝工程建设及运行借鉴和参考。     

模拟三维裂纹问题的扩展有限元法

扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其计算网格与不连续面相互独立,因此模拟移动不连续面时无需对网格进行重新剖分。给出了模拟三维裂纹问题的扩展有限元法。在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性。用两个水平集函数表示裂纹。采用线性互补法求解裂纹面非线性接触条件,不需要迭代,提高了计算效率。采用两点位移外推法计算裂纹前缘应力强度因子。给出了3个三维弹性静力问题算例,其结果显示了所提方法能获得高精度的应力强度因子,并能有效地处理裂纹面间的接触问题,同时表明扩展有限元结合线性互补法求解不连续问题具有较好的前景。     

基于随机有限元法的溶蚀坝基系统建坝适应性评价

依托溶蚀地基上建坝的典型工程问题,在研究坝基溶蚀空间展布特征的基础上,建立溶蚀带力学参数随机场模型,采用随机有限元方法研究坝体工作性态的统计特性,评价了溶蚀坝基建坝适应性。采用溶蚀带力学参数随机场模型,等效模拟了空间内溶蚀体的无序分布系统,排除地质统计中线溶蚀率的人为假定,避免了以往确定性分析模型中溶蚀体缺失现象。研究方法和成果亦可以为其他相似复杂地基工程的数值分析提供新的思路和参考。     

Stress and Deformation Analysis of Concrete-Facing Sand–Gravel Dam Based on Inversion Parameters

With the long-term operation of the project, the material parameters of concrete-facing sand–gravel dam will change, which brings great difficulty to the scientific and effective stress and deformation analysis. Combining with the measured displacement data, the finite element analysis model of the concrete-facing sand–gravel dam of Heiquan reservoir was established, and the modulus of elasticity and internal friction angle of the dam body were inverted by the measured displacement of the dam, then the simulation analysis of the filling construction process and the reservoir storage process of dam was carried out, and the stress and deformation values of the dam during the construction period and the impoundment period were calculated. The results showed that the parameters obtained from the inversion are smaller than the original parameters, but there is little difference between them. The displacement calculated by finite element inversion was close to the measured displacement value, the overall displacement and stress distribution of the dam body and panel were in line with the general law, and the calculated displacement and stress values were at the normal level. This study provides a reference for parameter inversion and stress and deformation analysis of concrete slab dam through monitoring data analysis.

     

考虑桩-土非线性接触的自平衡桩基测试有限元分析

自平衡桩基测试由于不受场地限制,在桥梁桩基检测和水上试桩领域得到较为广泛应用,但由于自身承受荷载机理、桩土相互作用与传统桩基静载测试不同,目前其实际应用仍然存在一些问题,因此需要在自平衡加载方式下建立合理数值模型来为进一步理论分析奠定模型基础。在考虑桩体非线性接触面基础上,研究了决定非线性接触性质的摩擦系数计算原则及相关模型参数取值,基于ANSYS接触面单元建立了自平衡桩基有限元模型,以武汉工程为依托,与实测结果进行对比,结果表明基于接触面单元建立的自平衡桩基模型具有更高的精确度,能够很好反映真实荷载作用下上下段桩位移曲线图。     

Finite element modelling of frictional instability between deformable rocks

Earthquakes are recognized as resulting from a stick–slip frictional instability along faults. Based on the node‐to‐point contact element strategy (an arbitrarily shaped contact element strategy applied with the static‐explicit algorithm for modelling non‐linear frictional contact problems proposed by authors), a finite element code for modelling the 3‐D non‐linear friction contact between deformable bodies has been developed and extended here to analyse the non‐linear stick–slip frictional instability between deformable rocks with a rate‐ and state‐dependent friction law. A typical fault bend model is taken as an application example to be analysed here. The variations of the normal contact force, the frictional force, the transition of stick–slip instable state and the related relative slip velocity along the fault between the deformable rocks and the stress evolution in the total bodies during the different stages are investigated, respectively. The calculated results demonstrate the usefulness of this code for simulating the non‐linear frictional instability between deformable rocks. Copyright © 2003 John Wiley & Sons, Ltd.     

隧道围岩弹塑性随机有限元分析及可靠度计算

在弹塑性有限元分析的初应力法基础上,推导了弹塑性随机有限元增量初应力法的计算迭代格式,结合可靠度理论,给出了围岩稳定的可靠度分析方法,在计算出可靠度的同时还给出了围岩和锚喷支护结构的应力特征值。并对一工程算例进行了详细的分析和讨论,得出了一些有益的结论。     

Effect of particle friction and polydispersity on the macroscopic stress–strain relations of granular materials

The macroscopic mechanical behavior of granular materials inherently depends on the properties of particles that compose them. Using the discrete element method, the effect of particle contact friction and polydispersity on the macroscopic stress response of 3D sphere packings is studied. The analytical expressions for the pressure, coordination number and fraction of rattlers proposed for isotropically deformed frictionless systems also hold when the interparticle coefficient of friction is finite; however, the numerical values of the parameters such as the jamming volume fraction change with varying microscopic contact and particle properties. The macroscopic response under deviatoric loading is studied with triaxial test simulations. Concerning the shear strength, our results agree with previous studies showing that the deviatoric stress ratio increases with particle coefficient of friction μ starting from a nonzero value for μ = 0 and saturating for large μ. On the other hand, the volumetric strain does not have a monotonic dependence on the particle contact friction. Most notably, maximum compaction is reached at an intermediate value of the coefficient of friction μ ≈ 0.3. The effect of polydispersity on the macroscopic stress–strain relationship cannot be studied independent of initial packing conditions. The shear strength increases with polydispersity when the initial volume fraction is fixed, but the effect of polydispersity is much less pronounced when the initial pressure of the packings is fixed. Finally, a simple hypoplastic constitutive model is calibrated with numerical test results following an established procedure to ascertain the relation between particle properties and material coefficients of the macroscopic model. The calibrated model is in good qualitative agreement with simulation results.     

基于抗转模型的颗粒材料宏−细观关系研究

接触模型的宏?细观参数标定是成功使用离散元方法的关键。在离散元的接触模型中线性接触模型与抗转线性接触模型均可用于模拟砂性土的力学行为,其中抗转线性接触模型在模拟密砂的剪胀性方面具备优势。采用抗转线性接触模型对室内密实砂土三轴试验进行了离散元模拟,验证了抗转线性接触模型的可靠性;进而系统分析了颗粒间摩擦系数、刚度比和抗转动系数等细观参数与砂土峰值内摩擦角、残余内摩擦角、峰值剪胀角等宏观参数的相关关系并进行了验证;揭示了偏应力作用下,细观参数对密实砂土试样内部剪切带宽度与倾角变化的影响规律,提出了考虑剪胀角的剪切带倾角经验公式。通过研究建立了抗转线性接触模型宏?细观参数的量化关系并给出了标定参数的具体流程图,提出了快速标定宏观参数的方法并应用实例进行了验证,为采用抗转线性接触模型精准模拟密实砂土的力学特性提供依据。     

龙井水库大坝的渗流和稳定分析

本文以龙井水库大坝为例,运用有限单元法对粘土心墙坝的渗流进行分析计算,同时运用瑞典圆弧滑动法和简化Bishop法对粘土心墙坝的坝体稳定进行分析计算,在此基础上对水库大坝进行了渗流安全评价和稳定性安全评价。