饱和砂土动静力特性的网格自适应数值模拟

›› 2014, Vol. 35 ›› Issue (7): 2079-2087.

饱和砂土动静力特性的网格自适应数值模拟

邵琪^{1, 2}，唐小微^{1, 2}

1. 大连理工大学海岸与近海工程国家重点实验室，辽宁大连 116024； 2. 大连理工大学土木工程学院岩土工程研究所，辽宁大连 116024

收稿日期:2013-07-16 出版日期:2014-07-10 发布日期:2014-07-11
通讯作者: 唐小微，男，1968年生，博士，副教授，博导，主要从事土动力学、海洋土动力学等方面的研究。E-mail: tangxw@dlut.edu.cn
作者简介:邵琪，男，1983年生，博士研究生，主要从事地震液化数值模拟、网格自适应方法等方面的研究。
基金资助:
国家高技术研究发展计划863资助项目(No. 2012AA112510)。

Numerical simulation of static and dynamic behaviors of saturated sand by an adaptive mesh scheme

SHAO Qi^{1, 2}，TANG Xiao-wei^{1, 2}

1. State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian, Liaoning 116024, China; 2. Institute of Geotechnical Engineering, School of Civil Engineering, Dalian University of Technology, Dalian, Liaoning 116024, China

Received:2013-07-16 Online:2014-07-10 Published:2014-07-11

摘要/Abstract

摘要： 在有限元数值计算平台上，建立了一套基于后验误差评估和Delaunay三角剖分的网格自适应方案，针对饱和砂土静力受压和地震液化的特性进行模拟。验证了超收敛单元片回归（SPR）误差评估中原用于四边形单元的双线性回归函数在用于三角形单元网格时的适用性和可靠性；在饱和砂土动、静力算例中，网格自适应计算获得的变形、应变、超孔压比等的变化规律与常规有限元结果趋势一致。随着网格的再生成，参考点的位移和全域的平均相对误差逼近精确值。对于初始网格，讨论了合理的自适应程度并应用于地震液化的自适应数值模拟中，也对Delaunay三角剖分实施了一些改进。最终证明该自适应方案在提高计算效率的同时，亦可以保证计算所需的精准度。

关键词: 后验误差评估, Delaunay三角剖分, 网格自适应, 饱和砂土, 静力受压, 地震液化

Abstract: Based on the finite element method with a posterior error estimation and the Delaunay triangulation, an adaptive mesh scheme is set up to simulate the saturated sand behaviors when being compressed statically and when being liquefied during earthquakes. The bilinear recovered function, usually being used for quadrilateral elements in superconvergent patch recovery (SPR) error estimation, is confirmed to be available and reliable for triangular element mesh. In both the static and dynamic examples of saturated sand, the changing laws of deformation, strain, excess pore water pressure ratio, etc. due to the adaptive mesh calculation are maintained in the same trend as the normal finite element calculation. The displacement of the reference point is approaching to the exact value as the mesh is regenerated, so is the average relative error of the whole area. The suitable adaptation degree for an initial mesh is discussed and applied to the adaptive numerical simulation of the seismic liquefaction. Also, some improvement for the initial mesh is implemented in the Delaunay triangulation. This adaptive mesh scheme is proved to improve calculating efficiency, and to ensure the desired accuracy as well.

Key words: posterior error estimation, Delaunay triangulation, adaptive remesh, saturated sand, static compression, seismic liquefaction

中图分类号:

TU 435

邵琪，唐小微 , . 饱和砂土动静力特性的网格自适应数值模拟[J]. , 2014, 35(7): 2079-2087.

SHAO Qi ，TANG Xiao-wei , . Numerical simulation of static and dynamic behaviors of saturated sand by an adaptive mesh scheme[J]. , 2014, 35(7): 2079-2087.

参考文献

相关文章 15

[1]	董建勋, 刘海笑, 李洲. 适用于砂土循环加载分析的边界面塑性模型[J]. 岩土力学, 2019, 40(2): 684-692.
[2]	高源, 刘海笑, 李洲, . 适用于饱和砂土循环动力分析边界面塑性模型的显式积分算法[J]. 岩土力学, 2019, 40(10): 3951-3958.
[3]	周燕国，谭晓明，梁甜，黄博，凌道盛，陈云敏，. 利用地震动强度指标评价场地液化的离心模型试验研究[J]. , 2017, 38(7): 1869-1877.
[4]	孙奇，董全杨，蔡袁强，王军, . 扭剪荷载作用下饱和南京砂不排水特性试验研究[J]. , 2017, 38(1): 41-48.
[5]	王成龙，刘汉龙，孔纲强，吴迪, . 工作荷载下温度循环对桩基变形与应力的影响分析[J]. , 2016, 37(S1): 317-322.
[6]	刘汉龙，王成龙，孔纲强，吴迪,. U型、W型和螺旋型埋管形式能量桩热力学特性对比模型试验[J]. , 2016, 37(S1): 441-447.
[7]	周葛，刘汉龙，井合进，陈育民，飞田哲男,. 基于应变空间多机构CG模型的地震液化大变形数值分析[J]. , 2016, 37(8): 2417-2425.
[8]	胡记磊，唐小微，裘江南，. 基于贝叶斯网络的地震液化概率预测分析[J]. , 2016, 37(6): 1745-1752.
[9]	武军，廖少明，霍晓波,. 地铁振动荷载对下穿越盾构开挖面孔隙水压力的影响[J]. , 2015, 36(S1): 496-500.
[10]	孙奇，董全杨，蔡袁强，王军，胡秀青，. 偏应力比和应力主轴偏转角同时变化下饱和砂土变形特性研究[J]. , 2015, 36(8): 2261-2269.
[11]	阳军生，张箭，杨峰. 浅埋隧道掌子面稳定性二维自适应上限有限元分析[J]. , 2015, 36(1): 257-264.
[12]	金书成，张永涛，杨炎华，黎冰 . 饱和砂土地基中吸力式桶形基础水平承载力研究[J]. , 2013, 34(S1): 221-227.
[13]	杨永香，周健，贾敏才，胡金虎. 饱和砂土液化特性的可视化试验研究[J]. , 2011, 32(6): 1643-1648.
[14]	夏志凡，叶冠林，王建华，叶斌，张锋. 防水措施对堤坝地震液化影响的数值分析[J]. , 2010, 31(8): 2682-2688.
[15]	张孟喜，陈高峰，朱引，李晋杰. H-V加筋饱和砂土性状的三轴试验研究[J]. , 2010, 31(5): 1345-1351.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed