首页|非均质地层网格加密与稀疏并举的高阶单元自适应上限有限元研究

非均质地层网格加密与稀疏并举的高阶单元自适应上限有限元研究

扫码查看
岩土极限上限有限元网格自适应策略常以"后验"方式捕捉塑性破坏区,可能引起刚性与塑性缓冲区过度加密.同时,自适应上限有限元初始网格不宜过于稀疏,以避免初始解误差过大及网格加密路径偏离方向.为此,在应用六节点三角形高阶单元和二阶锥规划模型的基础上,以耗散能密度指标区分活跃与非活跃单元,提出一种应用于上限有限元加密与稀疏并举的网格自适应策略;即每次自适应过程针对性地加密与稀疏不同区域网格,达到降低模型单元总数,同时提高计算精度的目的.进一步地,将网格自适应策略推广到可考虑强度随深度线性增加的非均质地层,给出上限有限元的编程实现流程,最后以非均质地层隧道开挖面稳定性和主动开拓门失稳破坏算例,开展多参数的对比分析,验证网格自适应方法高效、可靠.
A simultaneous adaptive mesh refining and coarsening algorithm in high-order element-based upper bound finite element method in non-homogeneous soils
In geotechnical engineering,upper bound finite element method with adaptive mesh refinement often evaluates the plastic failure area in a posterior manner,which may cause an excessive refinement in the buffer zone between rigid and plastic zones and lead to a high computational cost.Meanwhile,the initial mesh for the adaptive upper bound finite element method should not be too sparse to avoid the large error caused by the over-rigidity of model,as it can induce subsequent adaptive refinement to deviate from the correct direction.In this study,with a combination of six-node triangular elements and a second-order cone programming model,an upper bound finite element method with simultaneous adaptive mesh refining and coarsening algorithm is proposed.This method can automatically determine both active elements and inactive elements using plastic power dissipation-based estimators.For each adaptive step,those active elements are automatically refined and those inactive ones are thus coarsened,which can achieve the purpose of reducing the total number of elements and improving computational accuracy.The proposed method is then extended to the stability analysis in non-homogeneous soils,and numerical implementation of the proposed method is discussed.Two examples including tunnel face stability and failure of active trapdoors in non-homogeneous soils are analyzed to evaluate the validity of the proposed method through a series of parameter analysis.

numerical analysisadaptive upper bound finite element methodmesh refinementmesh coarseningnon-homogeneous soilsix-node triangular elementsecond-order cone programmingfailure mode

郑响凑、秦傲韩、杨峰、阳军生

展开 >

中南大学土木工程学院,湖南长沙 410075

数值分析 自适应上限有限元 网格加密 网格稀疏 非均质地层 六节点三角形单元 二阶锥规划 破坏模式

2024

岩石力学与工程学报
中国岩石力学与工程学会

岩石力学与工程学报

CSTPCD北大核心
影响因子:2.589
ISSN:1000-6915
年,卷(期):2024.43(z2)