边坡稳定上限极限分析通常采用三节点三角形单元,并对屈服准则线性化,因而其计算精度较低。针对此问题,基于六节点三角形单元和二阶锥规划,假设材料严格满足Mohr-Coulomb屈服准则并考虑孔隙水压力与地震荷载的作用,发展了一种应用于边坡稳定性分析的上限有限元法,扩大了上限极限分析的应用范围。根据上限定理,在满足屈服条件、流动法则、功能平衡以及相应的边界条件的基础上,将边坡稳定的上限分析形成二阶锥规划数学模型,并用先进的内点法进行求解。通过2个算例的计算分析,并与已有计算方法进行对比,验证了本文方法的正确性与可行性,并表明该方法的计算结果并不明显依赖于有限元网格的密度,且材料内摩擦角较大时,仍具有较高的计算精度。
Abstract
Three-node triangular elements together with the linearized yield criterion are commonly used in the finite element upper bound limit analysis. Therefore, the method is of low calculation precision. Aiming at this problem, a method of upper bound limit analysis using six-node triangular elements and second-order cone programming is developed to investigate the slope stability subjected to pore water pressure and earthquake loads. The proposed method formulates the slope stability problem as a second-order cone programming with constraints based on the yield criterion, flow rule, boundary conditions, and the energy-work balance equation. The optimization problem is solved by a state-of-the-art interior-point method, and the strict upper bound solutions can be obtained. Finally, the results of two numerical examples are compared with published solutions, which demonstrate the validity of the proposed method. The results also indicate that the mesh-dependence phenomenon is overcome and the calculation precision is improved even for large internal friction angel of materials.
关键词
边坡稳定性 /
极限分析 /
上限定理 /
二阶锥规划 /
有限元
Key words
slope stability /
limit analysis /
upper bound theorem /
second-order cone programming /
finite element
{{custom_sec.title}}
{{custom_sec.title}}
{{custom_sec.content}}
参考文献
[1] SLOAN S W, KLEEMAN P W. Upper Bound Limit Analysis Using Discontinuous Velocity Fields[J]. Computer Methods in Applied Mechanics and Engineering, 1995, 127(1/4): 293-314.
[2] KIM J M, SALGADO R, YU H S. Limit Analysis of Soil Slopes Subjected to Pore-water Pressures[J]. Journal of Geotechnical and Geoenvironmental Engineering, 1999, 125(1): 49-58.
[3] LOUKIDIS D, BANDINI P, SALGADO R. Stability of Seismically Loaded Slopes Using Limit Analysis[J]. Geotechnique, 2003, 53(5): 463-479.
[4] YU H S, SLOAN S W, KLEEMAN P W. A Quadratic Element for Upper Bound Limit Analysis[J]. Engineering Computations, 1994, 11(3): 195-212.
[5] 杨 峰,阳军生,李昌友,等. 基于六节点三角形单元和线性规划模型的上限有限元研究[J]. 岩石力学与工程学报, 2012, 31(12): 2556-2563.
[6] 孙 聪,李春光,郑 宏,等. 基于四边形网格的边坡上限有限元法[J]. 岩石力学与工程学报, 2015, 34(1): 114-120.
[7] 李 泽,王均星. 基于非线性规划的岩质边坡有限元塑性极限分析下限法研究[J]. 岩石力学与工程学报, 2007, 26(4): 747-753.
[8] 杨昕光,迟世春. 土石坝坡极限抗震能力的下限有限元法[J]. 岩土工程学报, 2013, 35(7): 1202-1209.
[9] LYAMIN A V, SLOAN S W. Upper Bound Limit Analysis Using Linear Finite Elements and Non-linear Programming[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2002, 26(2): 181-216.
[10]LI H X, YU H S. Kinematic Limit Analysis of Frictional Materials Using Nonlinear Programming[J]. International Journal of Solids and Structures, 2005, 42(14): 4058-4076.
[11]周建烽,王均星,陈 炜,等. 线性与非线性强度的土石坝坝坡稳定分析下限法[J]. 岩土力学, 2015, 36(1): 233-239.
[12]MAKRODIMOPOULOS A, MARTIN C M. Upper Bound Limit Analysis Using Simplex Strain Elements and Second-order Cone Programming[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2007, 31(6): 835-865.
[13]MAKRODIMOPOULOS A, MARTIN C M. Upper Bound Limit Analysis Using Discontinuous Quadratic Displacement Fields[J]. Communications in Numerical Methods in Engineering, 2007, 24(11): 911-927.
[14]MICHALOWSKI R L. Slope Stability Analysis: A Kinematic Approach[J]. Geotechnique, 1995, 45(2): 283-293.
[15]李国英,沈珠江. 下限原理有限单元法及其在土工问题中的应用[J]. 岩土工程学报, 1997, 19(5): 86-91.
基金
国家自然科学基金项目(51509019);raybet体育在线
创新团队项目(CKSF2015051/YT);水利部土石坝破坏机理与防控技术重点实验室开放研究基金项目(YK914017)
PDF(1016 KB)
