西安建筑科技大学学报(自然科学版)

弱形式求积元法是一种新型的高阶方法,已经在岩土及结构分析中取得了成功的应用,本文对该方法进行进一步的发展和完善,将其应用于饱和土动力固结分析.首先基于Biot饱和土波动理论框架,以孔隙水压力和土体骨架位移为基本控制变量,建立饱和土波动问题控制方程弱形式描述,然后应用Lobatto积分和微分求积法进行数值积分和数值微分并采用Newmark方法进行时域逐步积分,建立饱和土波动问题求积元法求解列式.通过数值算例验证了本文方法的正确性,显示了方法的计算效率.

The weak form quadrature element method is a novel high order algorithm applied successfully in structural and geotechnical engineering. In the present study, it is reformulated for dynamic consolidation analysis of saturated soils. Based on Biot's theory of dynamic consolidation, the pore pressure and soil displacement are chosen to be the control variables to establish the weak form governing equations. Then the Lobatto integration rule and the differential quadrature method are employed respectively to numerically integrate and differentiate the weak form governing equations, and the Newmark scheme is used for time integration. The established formulation is verified by numerical examples and its efficiency is highlighted.

引言
1 动力固结问题控制方程
2 弱形式求积元法求解列式
3 数值算例
4 总结

图1 单相材料波动问题示意图<br/>Fig.1 Wave in the single phase material

图1 单相材料波动问题示意图
Fig.1 Wave in the single phase material

图2 一维波动问题计算结果对比<br/>Fig.2 Results of one-dimensional consolidation

图2 一维波动问题计算结果对比
Fig.2 Results of one-dimensional consolidation

图3 y=9 m处孔压时程曲线<br/>Fig.3 Pore pressure at y=9 m

图3 y=9 m处孔压时程曲线
Fig.3 Pore pressure at y=9 m

图4 t=0.07 s时孔压分布<br/>Fig.4 Pore pressure at t=0.07 s

图4 t=0.07 s时孔压分布
Fig.4 Pore pressure at t=0.07 s

[1] 徐洋, 谢康和, 卢廷浩. 二灰土桩复合地基三维固结有限元分析[J]. 岩土工程学报, 2002, 24(2):254-256.
[2] 王刚, 张建民. 波浪作用下某防沙堤的动力固结有限元分析[J]. 岩土力学, 2006, 27(4):555-560.
[3] 林丰, 陈环. 真空和堆载作用了砂井地基固结的边界元分析[J]. 岩土工程学报, 1987, 9(4):13-22.
[4] 欧贵宝, 饶意忠, 周博. 土体固结问题的一种边界元解法[J]. 哈尔滨工程大学学报, 2002, 23(6):39-42.
[5] WANG W, WANG J, WANG Z, et al. An unequal-order radial interpolation meshless method for Biot's consolidation theory[J]. Computers and Geotechnics, 2007, 34(2):61-70.
[6] 李冰河, 王奎华. 软粘土非线性一维固结有限差分法分析[J]. 浙江大学学报(工学版), 2000, 34(4):376-381.
[7] CHEN R, ZHOU W, WANG H, et al. One-dimensional nonlinear consolidation of multi-layered soil by differential quadrature method[J]. Computers and Geotechnics, 2005, 32(5):358-369.
[8] ZHENG G Y, LI P, ZHAO C Y. Analysis of non-linear consolidation of soft clay by differential quadrature method[J]. Applied Clay Science, 2013, 79:2-7.
[9] ZHONG H Z, YUE Z G. Analysis of thin plates by the weak form quadrature element method[J]. Science China Physics Mechanics & Astronomy, 2012, 55(5):861-871.
[10]HE R, ZHONG H. Large deflection elasto-plastic analysis of frames using the weak-form quadrature element method[J]. Finite Elements in Analysis & Design, 2012, 50(1):125-133.
[11]SHUAI Y, ZHONG H. Consolidation analysis of non-homogeneous soil by the weak form quadrature element method[J]. Computers & Geotechnics, 2014, 62:1-10.
[12]SHUAI Y, ZHONG H. Three dimensional analysis of unconfined seepage in earth dams by the weak form quadrature element method[J]. Journal of Hydrology, 2016, 533:403-411.
[13]BIOT M A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range[J]. The Journal of the acoustical Society of america, 1956, 28(2): 179-191.

备注

引言

1 动力固结问题控制方程

2 弱形式求积元法求解列式

3 数值算例

4 总结

期刊信息