In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
An efficient Galerkin meshfree formulation for three dimensional simulation of large deformation failure evolution in soils is presented. This formulation utilizes the stabilized conforming nodal integration, where fo...An efficient Galerkin meshfree formulation for three dimensional simulation of large deformation failure evolution in soils is presented. This formulation utilizes the stabilized conforming nodal integration, where for the purpose of stability and efficiency a Lagrangian smoothing strain at nodal point is constructed and thereafter the internal energy is evaluated nodally. This formulation ensures the linear exactness, efficiency and spatial stability in a unified manner and it makes the conventional Galerkin meshfree method affordable for three dimensional simulation. The three dimensional implementation of stabilized conforming nodal integration is discussed in details. To model the failure evolution in soil medium a coupled elasto-plastic damage model is used and an objective stress integration algorithm in combination of elasto-damage predictor and plastic corrector method is employed for stress update. Two typical numerical examples are shown to demonstrate the effectiveness of the present method for modeling large deformation soil failure.展开更多
A three-dimensional two-level gradient smoothing meshfree method is presented for rainfall induced landslide simulations. The two-level gradient smoothing for meshfree shape function is elaborated in the threedimensio...A three-dimensional two-level gradient smoothing meshfree method is presented for rainfall induced landslide simulations. The two-level gradient smoothing for meshfree shape function is elaborated in the threedimensional Lagrangian setting with detailed implementation procedure. It is shown that due to the successive gradient smoothing operation without the requirement of derivative computation in the present formulation, the two-level smoothed gradient of meshfree shape function is capable of achieving a given influence domain more efficiently than the standard gradient of meshfree shape function. Subsequently, the two-level smoothed gradient of meshfree shape function is employed to discretize the weak form of coupled rainfall seepage and soil motion equations in a nodal integration format, as provides an efficient three-dimensional regularized meshfree formulation for large deformation rainfall induced landslide simulations. The exponential damage and pressure dependent plasticity relationships are utilized to describe the failure evolution in landslides. The plastic response of soil is characterized by the true effective stress measure, which is updated according to the rotationally neutralized objective integration algorithm. The effectiveness of the present threedimensional two?level gradient smoothing meshfree method is demonstrated through numerical examples.展开更多
In order to overcome the possible singularity associated with the Point Interpolation Method(PIM),the Radial Point Interpolation Method(RPIM)was proposed by G.R.Liu.Radial basis functions(RBF)was used in RPIM as basis...In order to overcome the possible singularity associated with the Point Interpolation Method(PIM),the Radial Point Interpolation Method(RPIM)was proposed by G.R.Liu.Radial basis functions(RBF)was used in RPIM as basis functions for interpolation.All these radial basis functions include shape parameters.The choice of these shape parameters has been and stays a problematic theme in RBF approximation and interpolation theory.The object of this study is to contribute to the analysis of how these shape parameters affect the accuracy of the radial PIM.The RPIM is studied based on the global Galerkin weak form performed using two integration technics:classical Gaussian integration and the strain smoothing integration scheme.The numerical performance of this method is tested on their behavior on curve fitting,and on three elastic mechanical problems with regular or irregular nodes distributions.A range of recommended shape parameters is obtained from the analysis of different error indexes and also the condition number of the matrix system.All resulting RPIM methods perform very well in term of numerical computation.The Smoothed Radial Point Interpolation Method(SRPIM)shows a higher accuracy,especially in a situation of distorted node scheme.展开更多
利用添加多项式项的(Radial Point Interpolation Method,RPIM)形函数,形成了动力学问题的无网格全局弱形式.生成了伴随于域节点的Voronoi图,利用基于应变光滑稳定方案的稳定相容节点积分得到了改进后的总体刚度矩阵离散化形式,并利用...利用添加多项式项的(Radial Point Interpolation Method,RPIM)形函数,形成了动力学问题的无网格全局弱形式.生成了伴随于域节点的Voronoi图,利用基于应变光滑稳定方案的稳定相容节点积分得到了改进后的总体刚度矩阵离散化形式,并利用直接法施加位移边界条件.自由振动分析得到了与有限元参考解吻合良好的数值解,受迫振动分析采用了无条件稳定的Newmark法,从而验证了本方法在求解动力学问题所展现的稳定性、精确性及收敛性.展开更多
This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear a...This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear analysis encompasses the fully geometric nonlinearities due to large deflection,large deformation as well as finite rotation.The incremental equilibrium equation is obtained by the consistent linearization of the nonlinear variational equation.The Lagrangian meshfree shape function is utilized to discretize the variational equation.Subsequently to resolve the shear and membrane locking issues and accelerate the computation,the method of stabilized conforming nodal integration is systematically implemented through the Lagrangian gradient smoothing operation.Numerical results reveal that the present formulation is very effective.展开更多
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10972188, 10602049)the Program for New Century Excellent Talents in University from China Education Ministry (Grant No. NCET-09-0678)the Fundamental Research Funds for the Central Universities of China (Grant No. 2010121073)
文摘An efficient Galerkin meshfree formulation for three dimensional simulation of large deformation failure evolution in soils is presented. This formulation utilizes the stabilized conforming nodal integration, where for the purpose of stability and efficiency a Lagrangian smoothing strain at nodal point is constructed and thereafter the internal energy is evaluated nodally. This formulation ensures the linear exactness, efficiency and spatial stability in a unified manner and it makes the conventional Galerkin meshfree method affordable for three dimensional simulation. The three dimensional implementation of stabilized conforming nodal integration is discussed in details. To model the failure evolution in soil medium a coupled elasto-plastic damage model is used and an objective stress integration algorithm in combination of elasto-damage predictor and plastic corrector method is employed for stress update. Two typical numerical examples are shown to demonstrate the effectiveness of the present method for modeling large deformation soil failure.
基金the National Natural Science Foundation of China (Grant Nos.11772280 and 11472233)he Program for Scientific and Technological Innovation Leading Talents of Fujian Province is gratefully acknowledged.
文摘A three-dimensional two-level gradient smoothing meshfree method is presented for rainfall induced landslide simulations. The two-level gradient smoothing for meshfree shape function is elaborated in the threedimensional Lagrangian setting with detailed implementation procedure. It is shown that due to the successive gradient smoothing operation without the requirement of derivative computation in the present formulation, the two-level smoothed gradient of meshfree shape function is capable of achieving a given influence domain more efficiently than the standard gradient of meshfree shape function. Subsequently, the two-level smoothed gradient of meshfree shape function is employed to discretize the weak form of coupled rainfall seepage and soil motion equations in a nodal integration format, as provides an efficient three-dimensional regularized meshfree formulation for large deformation rainfall induced landslide simulations. The exponential damage and pressure dependent plasticity relationships are utilized to describe the failure evolution in landslides. The plastic response of soil is characterized by the true effective stress measure, which is updated according to the rotationally neutralized objective integration algorithm. The effectiveness of the present threedimensional two?level gradient smoothing meshfree method is demonstrated through numerical examples.
文摘In order to overcome the possible singularity associated with the Point Interpolation Method(PIM),the Radial Point Interpolation Method(RPIM)was proposed by G.R.Liu.Radial basis functions(RBF)was used in RPIM as basis functions for interpolation.All these radial basis functions include shape parameters.The choice of these shape parameters has been and stays a problematic theme in RBF approximation and interpolation theory.The object of this study is to contribute to the analysis of how these shape parameters affect the accuracy of the radial PIM.The RPIM is studied based on the global Galerkin weak form performed using two integration technics:classical Gaussian integration and the strain smoothing integration scheme.The numerical performance of this method is tested on their behavior on curve fitting,and on three elastic mechanical problems with regular or irregular nodes distributions.A range of recommended shape parameters is obtained from the analysis of different error indexes and also the condition number of the matrix system.All resulting RPIM methods perform very well in term of numerical computation.The Smoothed Radial Point Interpolation Method(SRPIM)shows a higher accuracy,especially in a situation of distorted node scheme.
文摘利用添加多项式项的(Radial Point Interpolation Method,RPIM)形函数,形成了动力学问题的无网格全局弱形式.生成了伴随于域节点的Voronoi图,利用基于应变光滑稳定方案的稳定相容节点积分得到了改进后的总体刚度矩阵离散化形式,并利用直接法施加位移边界条件.自由振动分析得到了与有限元参考解吻合良好的数值解,受迫振动分析采用了无条件稳定的Newmark法,从而验证了本方法在求解动力学问题所展现的稳定性、精确性及收敛性.
基金supported by the National Natural Science Foundation of China (10972188)the Program for New Century Excellent Talents in University from China Education Ministry (NCET-09-0678)
文摘This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear analysis encompasses the fully geometric nonlinearities due to large deflection,large deformation as well as finite rotation.The incremental equilibrium equation is obtained by the consistent linearization of the nonlinear variational equation.The Lagrangian meshfree shape function is utilized to discretize the variational equation.Subsequently to resolve the shear and membrane locking issues and accelerate the computation,the method of stabilized conforming nodal integration is systematically implemented through the Lagrangian gradient smoothing operation.Numerical results reveal that the present formulation is very effective.
