Similar articles

Abstract: The weak form quadrature element method is a high-order algorithm which has been applied successfully to consolidation analysis of saturated and unsaturated soils. Its superiority over the conventional finite element method has been verified. However, in consolidation analysis, pore pressure oscillations will appear when small time increments are used due to the very small permeability and near-incompressibility of the pore water. This can produce almost zero values of main diagonal elements in the coefficient matrix. To overcome the pressure oscillations, we propose a coupled composite quadrature element in which different orders of integration are employed for the pore pressure and the displacement. Its performance is evaluated and compared with that of the standard element through 1D and 2D numerical tests. Our results show that pressure oscillations can be effectively alleviated and stability and accuracy can be significantly enhanced by using the proposed element.

一种非等阶积分的求积元孔压单元

[1]Aguilar G, Gaspar F, Lisbona F, et al., 2008. Numerical stabilization of Biot’s consolidation model by a perturbation on the flow equation. International Journal for Numerical Methods in Engineering, 75(11):1282-1300.

[2]Azari B, Fatahi B, Khabbaz H, 2015. Numerical analysis of vertical drains accelerated consolidation considering combined soil disturbance and visco-plastic behaviour. Geomechanics & Engineering, 8(2):187-220.

[3]Bellman R, Casti J, 1971. Differential quadrature and long-term integration. Journal of Mathematical Analysis and Applications, 34(2):235-238.

[4]Biot MA, 1941. General theory of three-dimensional consolidation. Journal of Applied Physics, 12(2):155-164.

[5]Boal N, Gaspar FJ, Lisbona FJ, et al., 2011. Finite-difference analysis of fully dynamic problems for saturated porous media. Journal of Computational and Applied Mathematics, 236(6):1090-1102.

[6]Chen R, Zhou W, Wang H, et al., 2005. One-dimensional nonlinear consolidation of multi-layered soil by differential quadrature method. Computers and Geotechnics, 32(5):358-369.

[7]Darkshanamurthy V, Fredlund DG, Rahardjo H, 1984. Coupled three-dimensional consolidation theory of unsaturated porous media. Fifth International Conference on Expansive Soils, p.99.

[8]Fei K, Liu HL, 2009. Implementation and application of bounding surface model in ABAQUS. Journal of PLA University of Science and Technology (Natural Science Edition), 10(5):447-451 (in Chinese).

[9]Gaspar FJ, Lisbona FJ, Vabishchevich PN, 2003. A finite difference analysis of Biot’s consolidation model. Applied Numerical Mathematics, 44(4):487-506.

[10]Gaspar FJ, Lisbona FJ, Vabishchevich PN, 2006. Staggered grid discretizations for the quasi-static Biot’s consolidation problem. Applied Numerical Mathematics, 56(6):888-898.

[11]Gaspar FJ, Lisbona FJ, Matus P, et al., 2016. Numerical methods for a one-dimensional non-linear Biot’s model. Journal of Computational and Applied Mathematics, 293: 62-72.

[12]Ho L, Fatahi B, 2015. Analytical solution for the two-dimensional plane strain consolidation of an unsaturated soil stratum subjected to time-dependent loading. Computers & Geotechnics, 67:1-16.

[13]Ho L, Fatahi B, Khabbaz H, 2016. Analytical solution to axisymmetric consolidation in unsaturated soils with linearly depth-dependent initial conditions. Computers & Geotechnics, 74:102-121.

[14]Hua L, 2011. Stable element—free Galerkin solution procedures for the coupled soil—pore fluid problem. International Journal for Numerical Methods in Engineering, 86(8):1000-1026.

[15]Le TM, Fatahi B, 2016. Trust-region reflective optimisation to obtain soil visco-plastic properties. Engineering Computations, 33(2):410-442.

[16]Le TM, Fatahi B, Khabbaz H, et al., 2016. Numerical optimisation applying trust-region reflective least squares algorithm with constraints to optimise the non-linear creep parameters of soft soil. Applied Mathematical Modelling, 16:1-21.

[17]Manzari MT, Dafalias YF, 1997. A critical state two-surface plasticity model for sands. Geotechnique, 47(2):255-272.

[18]Mira P, Pastor M, Li T, et al., 2003. A new stabilized enhanced strain element with equal order of interpolation for soil consolidation problems. Computer Methods in Applied Mechanics and Engineering, 192(37-38):4257-4277.

[19]Preisig M, Prévost JH, 2011. Stabilization procedures in coupled poromechanics problems: a critical assessment. International Journal for Numerical and Analytical Methods in Geomechanics, 35(11):1207-1225.

[20]Reed M, 1984. An investigation of numerical errors in the analysis of consolidation by finite elements. International Journal for Numerical and Analytical Methods in Geomechanics, 8(3):243-257.

[21]Rodrigo C, Gaspar F, Hu X, et al., 2016. Stability and monotonicity for some discretizations of the Biot’s consolidation model. Computer Methods in Applied Mechanics and Engineering, 298:183-204.

[22]Samimi S, Pak A, 2012. Three-dimensional simulation of fully coupled hydro-mechanical behavior of saturated porous media using element free Galerkin (EFG) method. Computers & Geotechnics, 46(4):75-83.

[23]Sandhu RS, Liu H, Singh KJ, 1977. Numerical performance of some finite element schemes for analysis of seepage in porous elastic media. International Journal for Numerical and Analytical Methods in Geomechanics, 1(2):177-194.

[24]Simo JC, Rifai M, 1990. A class of mixed assumed strain methods and the method of incompatible modes. International Journal for Numerical Methods in Engineering, 29(8):1595-1638.

[25]Wang W, Wang J, Wang Z, et al., 2007. An unequal-order radial interpolation meshless method for Biot’s consolidation theory. Computers and Geotechnics, 34(2):61-70.

[26]Wang X, Gu H, 1997. Static analysis of frame structures by the differential quadrature element method. International Journal for Numerical Methods in Engineering, 40(4):759-772.

[27]Yuan S, Zhong HZ, 2014. Consolidation analysis of non-homogeneous soil by the weak form quadrature element method. Computers and Geotechnics, 62:1-10.

[28]Yuan S, Zhong HZ, 2016a. Three dimensional analysis of unconfined seepage in earth dams by the weak form quadrature element method. Journal of Hydrology, 533:403-411.

[29]Yuan S, Zhong HZ, 2016b. A weak form quadrature element formulation for coupled analysis of unsaturated soils. Computers and Geotechnics, 76:1-11.

[30]Yuan S, Zhong HZ, 2017. Finite deformation elasto-plastic consolidation analysis of soft clay by the weak form quadrature element method. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering), 18(12):942-957.

[31]Zhong HZ, Yu T, 2009. A weak form quadrature element method for plane elasticity problems. Applied Mathematical Modelling, 33(10):3801-3814.