Abstract: To address the issue of spurious pore pressure oscillations in equal-order finite element analyses of the dynamic responses of low-permeability saturated porous media, we compare two primary formulations: (1) a novel unequal-order interpolation element enriched with cover functions; (2) a stabilized equal-order interpolation element. To circumvent inherent linear dependence issues in the cover function-enriched formulation, specialized quadrilateral and triangular elements are developed for 1D and 2D problems, respectively. The performance of each element type is evaluated for both compressible and incompressible materials. Our results demonstrate that the cover element consistently and accurately captures pore pressure evolution for both forms of material. In contrast, the stabilized equal-order element achieves acceptable accuracy only for incompressible materials with high skeleton stiffness and exhibits degraded performance for compressible materials. Despite efforts to mitigate the influence of the stabilization term, the numerical results remain prone to either oscillatory or excessively damped responses.

低渗透性饱和多孔介质动力分析有限元格式对比评估

[1]AguilarG, GasparF, LisbonaF, 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]AnXM, LiLX, MaGW, et al., 2011. Prediction of rank deficiency in partition of unity-based methods with plane triangular or quadrilateral meshes. Computer Methods in Applied Mechanics and Engineering, 200(5-8):665-674.

[3]BabuškaI, 1971. Error-bounds for finite element method. Numerische Mathematik, 16(4):322-333.

[4]BatheKJ, 2001. The inf–sup condition and its evaluation for mixed finite element methods. Computers & Structures, 79(2):243-252.

[5]BiotMA, 1956a. Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range. The Journal of the Acoustical Society of America, 28(2):168-178.

[6]BiotMA, 1956b. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range. The Journal of the Acoustical Society of America, 28(2):179-191.

[7]ChaiYB, BatheKJ, 2020. Transient wave propagation in inhomogeneous media with enriched overlapping triangular elements. Computers & Structures, 237:106273.

[8]ChaiYB, LiW, LiuZY, 2022. Analysis of transient wave propagation dynamics using the enriched finite element method with interpolation cover functions. Applied Mathematics and Computation, 412:126564.

[9]CorciuloS, ZanoliO, PisanòF, 2017. Transient response of offshore wind turbines on monopiles in sand: role of cyclic hydro–mechanical soil behaviour. Computers and Geotechnics, 83:221-238.

[10]de BoerR, EhlersW, LiuZF, 1993. One-dimensional transient wave propagation in fluid-saturated incompressible porous media. Archive of Applied Mechanics, 63(1):59-72.

[11]de-PouplanaI, OñateE, 2017. A FIC-based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems. International Journal for Numerical and Analytical Methods in Geomechanics, 41(1):110-134.

[12]FeulvarchE, BrosseA, VincentY, 2023. An LBB-stable P1/RNP0 finite element based on a pseudo-random integration method for incompressible and nearly incompressible material flows. International Journal for Numerical Methods in Engineering, 124(24):5558-5573.

[13]GavagninC, SanaviaL, de LorenzisL, 2020. Stabilized mixed formulation for phase-field computation of deviatoric fracture in elastic and poroelastic materials. Computational Mechanics, 65(6):1447-1465.

[14]GuiQ, LiW, ChaiYB, 2023a. The enriched quadrilateral overlapping finite elements for time-harmonic acoustics. Applied Mathematics and Computation, 451:128018.

[15]GuiQ, LiW, ChaiYB, 2024. Improved modal analyses using the novel quadrilateral overlapping elements. Computers & Mathematics with Applications, 154:138-152.

[16]GuiQ, WangYW, YouXY, et al., 2023b. Performance of the quadrilateral overlapping elements combined with implicit time integration method in solving the transient wave propagations. Engineering Analysis with Boundary Elements, 155:1013-1027.

[17]GuoGL, FallM, 2019. Modelling of preferential gas flow in heterogeneous and saturated bentonite based on phase field method. Computers and Geotechnics, 116:103206.

[18]HafezM, SolimanM, 1993. Numerical solution of the incompressible Navier–Stokes equations in primitive variables on unstaggered grids. In: Gunzburger MD, Nicolaides RA (Eds.), Incompressible Computational Fluid Dynamics. Cambridge University Press, Cambridge, UK, p.183-201.

[19]HamS, BatheKJ, 2012. A finite element method enriched for wave propagation problems. Computers & Structures, 94-95:1-12.

[20]HaoQS, YuJL, WangXN, et al., 2023. Stabilized low-order finite-element formulation for static and dynamic simulation of saturated soils based on a hybrid integration scheme. Computers and Geotechnics, 161:105596.

[21]HeLS, JingLR, FengMF, 2025. New stabilized mixed finite element methods for two-field poroelasticity with low permeability. Applied Mathematics and Computation, 494:129285.

[22]HidanoS, YamaguchiY, TakaseS, et al., 2024. Semi-implicit material point method for simulating infiltration-induced failure of unsaturated soil structures. International Journal for Numerical and Analytical Methods in Geomechanics, 48(10):2661-2699.

[23]HuangMS, ZienkiewiczOC, 1998. New unconditionally stable staggered solution procedures for coupled soil-pore fluid dynamic problems. International Journal for Numerical Methods in Engineering, 43(6):1029-1052.

[24]HuangMS, WuSM, ZienkiewiczOC, 2001. Incompressible or nearly incompressible soil dynamic behaviour—a new staggered algorithm to circumvent restrictions of mixed formulation. Soil Dynamics and Earthquake Engineering, 21(2):169-179.

[25]HuangMS, YueZQ, ThamLG, et al., 2004. On the stable finite element procedures for dynamic problems of saturated porous media. International Journal for Numerical Methods in Engineering, 61(9):1421-1450.

[26]KementzetzidisE, CorciuloS, VersteijlenWG, et al., 2019. Geotechnical aspects of offshore wind turbine dynamics from 3D non-linear soil-structure simulations. Soil Dynamics and Earthquake Engineering, 120:181-199.

[27]KimKT, ZhangLB, BatheKJ, 2018. Transient implicit wave propagation dynamics with overlapping finite elements. Computers & Structures, 199:18-33.

[28]KontoeS, ZdravkovicL, PottsDM, 2008. An assessment of time integration schemes for dynamic geotechnical problems. Computers and Geotechnics, 35(2):253-264.

[29]KularathnaS, LiangWJ, ZhaoTC, et al., 2021. A semi-implicit material point method based on fractional-step method for saturated soil. International Journal for Numerical and Analytical Methods in Geomechanics, 45(10):1405-1436.

[30]LiWT, WeiCF, 2018. Stabilized low-order finite elements for strongly coupled poromechanical problems. International Journal for Numerical Methods in Engineering, 115(5):531-548.

[31]LiuFS, 2020. Modeling hydraulic fracture propagation in permeable media with an embedded strong discontinuity approach. International Journal for Numerical and Analytical Methods in Geomechanics, 44(12):1634-1655.

[32]LiuFS, BorjaRI, 2010. Stabilized low-order finite elements for frictional contact with the extended finite element method. Computer Methods in Applied Mechanics and Engineering, 199(37-40):2456-2471.

[33]MaJW, DuanQL, TianR, 2022. A generalized finite element method without extra degrees of freedom for large deformation analysis of three-dimensional elastic and elastoplastic solids. Computer Methods in Applied Mechanics and Engineering, 392:114639.

[34]MarkertB, HeiderY, EhlersW, 2010. Comparison of monolithic and splitting solution schemes for dynamic porous media problems. International Journal for Numerical Methods in Engineering, 82(11):1341-1383.

[35]McGannCR, ArduinoP, Mackenzie-HelnweinP, 2015. A stabilized single-point finite element formulation for three-dimensional dynamic analysis of saturated soils. Computers and Geotechnics, 66:126-141.

[36]MelenkJM, BabuškaI, 1996. The partition of unity finite element method: basic theory and applications. Computer Methods in Applied Mechanics and Engineering, 139(1-4):289-314.

[37]MonforteL, NavasP, CarbonellJM, et al., 2019. Low-order stabilized finite element for the full Biot formulation in soil mechanics at finite strain. International Journal for Numerical and Analytical Methods in Geomechanics, 43(7):1488-1515.

[38]OñateE, 2000. A stabilized finite element method for incompressible viscous flows using a finite increment calculus formulation. Computer Methods in Applied Mechanics and Engineering, 182(3-4):355-370.

[39]PastorM, LiTC, MerodoJAF, 1997. Stabilized finite elements for harmonic soil dynamics problems near the undrained-incompressible limit. Soil Dynamics and Earthquake Engineering, 16(3):161-171.

[40]RajendranS, ZhangBR, 2007. A “FE-meshfree” QUAD4 element based on partition of unity. Computer Methods in Applied Mechanics and Engineering, 197(1-4):128-147.

[41]ShanZD, LingDS, DingHJ, 2011. Exact solutions for one-dimensional transient response of fluid-saturated porous media. International Journal for Numerical and Analytical Methods in Geomechanics, 35(4):461-479.

[42]ShiGH, 1991. Manifold method of material analysis. Transactions of the 9th Army Conference on Applied Mathematics and Computing, Minneapolis, USA, p.57-76.

[43]SonntagA, WagnerA, EhlersW, 2023. Dynamic hydraulic fracturing in partially saturated porous media. Computer Methods in Applied Mechanics and Engineering, 414:116121.

[44]StaubachP, MachačekJ, 2023. Separating fluid and solid contact constraints for hydro-mechanically coupled finite elements discretising fluid displacement. Computer Methods in Applied Mechanics and Engineering, 416:116451.

[45]SunW, OstienJT, SalingerAG, 2013. A stabilized assumed deformation gradient finite element formulation for strongly coupled poromechanical simulations at finite strain. International Journal for Numerical and Analytical Methods in Geomechanics, 37(16):2755-2788.

[46]TangTT, ZhouWX, LuoK, et al., 2024. Dynamic modeling and analysis of discontinuous wave propagation in a rod. Journal of Sound and Vibration, 569:117991.

[47]TaylorC, HoodP, 1973. A numerical solution of the Navier–Stokes equations using the finite element technique. Computers & Fluids, 1(1):73-100.

[48]TianR, YagawaG, TerasakaH, 2006. Linear dependence problems of partition of unity-based generalized FEMs. Computer Methods in Applied Mechanics and Engineering, 195(37-40):4768-4782.

[49]TrutyA, ZimmermannT, 2006. Stabilized mixed finite element formulations for materially nonlinear partially saturated two-phase media. Computer Methods in Applied Mechanics and Engineering, 195(13-16):1517-1546.

[50]WangZY, JinYF, YinZY, et al., 2023. A dynamic SNS-PFEM with generalized-α method for hydro-mechanical coupled geotechnical problems. Computers and Geotechnics, 159:105466.

[51]WhiteJA, BorjaRI, 2008. Stabilized low-order finite elements for coupled solid-deformation/fluid-diffusion and their application to fault zone transients. Computer Methods in Applied Mechanics and Engineering, 197(49-50):4353-4366.

[52]WuW, ZhengH, YangYT, 2019a. Enriched three-field numerical manifold formulation for dynamics of fractured saturated porous media. Computer Methods in Applied Mechanics and Engineering, 353:217-252.

[53]WuW, ZhengH, YangYT, 2019b. Numerical manifold method for dynamic consolidation of saturated porous media with three-field formulation. International Journal for Numerical Methods in Engineering, 120(6):768-802.

[54]XingJQ, ZhaoC, 2023. A hydro-mechanical phase field model for hydraulically induced fractures in poroelastic media. Computers and Geotechnics, 159:105418.

[55]XueHJ, ZhangXY, ChengJA, 2024. Transient analyses of underwater acoustic propagation with the modified meshfree radial point interpolation method and newmark time integration techniques. Ocean Engineering, 304:117842.

[56]YuJD, LiangWJ, ZhaoJD, 2025. Enhancing dynamic modeling of porous media with compressible fluid: a THM material point method with improved fractional step formulation. Computer Methods in Applied Mechanics and Engineering, 444:118100.

[57]YuanWH, ZhengHG, ZhengCX, et al., 2023. An improved semi-implicit material point method for simulating large deformation problems in saturated geomaterials. Computers and Geotechnics, 161:105614.

[58]YuanWH, ZhuJX, LiuK, et al., 2022. Dynamic analysis of large deformation problems in saturated porous media by smoothed particle finite element method. Computer Methods in Applied Mechanics and Engineering, 392:114724.

[59]YunG, LeeJ, KimDN, 2023. Stability of mixed overlapping elements in incompressible analysis. Computer Methods in Applied Mechanics and Engineering, 412:116104.

[60]ZhangHW, ZhouL, 2006. Numerical manifold method for dynamic nonlinear analysis of saturated porous media. International Journal for Numerical and Analytical Methods in Geomechanics, 30(9):927-951.

[61]ZhaoYD, ChooJ, 2020. Stabilized material point methods for coupled large deformation and fluid flow in porous materials. Computer Methods in Applied Mechanics and Engineering, 362:112742.

[62]ZienkiewiczOC, ChangCT, BettessP, 1980. Drained, undrained, consolidating and dynamic behaviour assumptions in soils. Géotechnique, 30(4):385-395.

[63]ZienkiewiczOC, HuangMS, WuJ, et al., 1993. A new algorithm for the coupled soil–pore fluid problem. Shock and Vibration, 1(1):3-14.