TY - JOUR
T1 - An unequal-order radial interpolation meshless method for Biot's consolidation theory
AU - Wang, W.D.
AU - Wang, Jian-Guo
AU - Wang, Z.L.
AU - Nogami, T.
PY - 2007
Y1 - 2007
N2 - Due to its second-order accuracy and unconditional stability, the Crank-Nicholson scheme is widely used to obtain numerical solutions of Biot's consolidation theory. However, spurious oscillation may occur when interpolations have equal orders for displacement and pore water pressure. This spurious oscillation may substantially degrade numerical accuracy. In order to avoid or alleviate spurious oscillation and improve numerical accuracy, a radial point interpolation meshless technique is proposed, with a polynomial for displacement that is one-order higher than for pore pressure. Three examples are numerically investigated with this meshless technique: a typical Terzaghi consolidation problem, a 2D consolidation problem, and a multi-layer consolidation problem. The numerical results are compared with the closed-form solution, the equal-order radial point interpolation method or FEM. The effects of the parameters of the current meshless method on spurious oscillation are studied. It is found that such an interpolation has consistent derivatives for both displacement and pore water pressure in an equilibrium equation. Spurious oscillation can thus be alleviated or even avoided, leading to high accuracy. (C) 2006 Elsevier Ltd. All rights reserved.
AB - Due to its second-order accuracy and unconditional stability, the Crank-Nicholson scheme is widely used to obtain numerical solutions of Biot's consolidation theory. However, spurious oscillation may occur when interpolations have equal orders for displacement and pore water pressure. This spurious oscillation may substantially degrade numerical accuracy. In order to avoid or alleviate spurious oscillation and improve numerical accuracy, a radial point interpolation meshless technique is proposed, with a polynomial for displacement that is one-order higher than for pore pressure. Three examples are numerically investigated with this meshless technique: a typical Terzaghi consolidation problem, a 2D consolidation problem, and a multi-layer consolidation problem. The numerical results are compared with the closed-form solution, the equal-order radial point interpolation method or FEM. The effects of the parameters of the current meshless method on spurious oscillation are studied. It is found that such an interpolation has consistent derivatives for both displacement and pore water pressure in an equilibrium equation. Spurious oscillation can thus be alleviated or even avoided, leading to high accuracy. (C) 2006 Elsevier Ltd. All rights reserved.
U2 - 10.1016/j.compgeo.2006.10.010
DO - 10.1016/j.compgeo.2006.10.010
M3 - Article
SN - 0266-352X
VL - 34
SP - 61
EP - 70
JO - Computers and Geotechnics
JF - Computers and Geotechnics
IS - 2
ER -