The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess v...The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.展开更多
A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The ...A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.展开更多
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.展开更多
This paper investigates the strength and critical state(CS)behaviours of granular materials via DEM simulations of true triaxial drained tests under three different loading modes including constant b(in-termediate str...This paper investigates the strength and critical state(CS)behaviours of granular materials via DEM simulations of true triaxial drained tests under three different loading modes including constant b(in-termediate stress ratio)tests with constantσ'1(major principal stress),constant p(mean pressure)and constantσ'3(minor principal stress)respectively.To this end,a series of samples are generated with the same particle size distribution,and with the confining stresses ranging from 100 kPa to 900 kPa.The CS is achieved for all samples.Both the macroscopic behaviours and the microscopic behaviours are examined and compared considering different loading modes,confining stresses and intermediate stress ratios(b).The critical state lines(CSLs)are found to be unique and independent of the loading modes,but dependent on the b values.The CSLs with b=0 and b=1 form the two boundaries of CSLs respectively beyond which CSLs under all other b tests cannot go beyond.Six different strength criteria are examined and compared in terms of both peak and CS failures.The Mohr-Coulomb strength criterion is found to be only suitable for axisymmetric loading conditions.The Lade-Duncan criterion is only suitable for describing peak strengths,which is dependent on the loading modes and confining stresses.The Satake criterion and Matsuoka-Nakai criterion are the more appropriate strength criteria for describing CS failures,indicating that the CS values of both the Satake parameter and the Matsuoka-Nakai parameter describe an inherent property that characterizes the CS failure for a given type of soil.The CS mechanical coordination number is fitted by a curved line for a given b value,which is unique regardless of the loading modes.The peak and CS values of both major and minor principal fabric tensors decrease with increasing b values,while they increase with increasing b values for the intermediate principal fabric tensor.展开更多
This article studies the influences of particle morphology on the behaviors of granular materials at both macroscopic and microscopic levels based on the discrete element method(DEM).A set of numerical tests under dra...This article studies the influences of particle morphology on the behaviors of granular materials at both macroscopic and microscopic levels based on the discrete element method(DEM).A set of numerical tests under drained triaxial compression was performed by controlling two morphological descriptors,i.e.ratio of the smallest to the largest pebble diameter,x,and the maximum pebbleepebble intersection angle,b.These descriptors are vital in generating particle geometry and surface textures.It was found that the stress responses of all assemblies exhibited similar behavior and showed post-peak strainsoftening.The normalized stress ratio and volumetric strains flatten off and tended to reach a steady value after an axial strain of 40%.While the friction angles at peak state varied with different morphological descriptors,the friction angles at critical state showed no significant variation.Moreover,evolution of the average coordination numbers showed a dramatic exponential decay until an axial strain of about 15%after which it stabilized and was unaffected by further increase of axial strain.In addition,stress ratio q/p and strong fabric parameter fs d=fs m were found to follow an approximately linear relationship for each assembly.These findings emphasized the significance of the influences of particle morphology on the macroscopic and microscopic responses of granular materials.展开更多
For geohazards and geotechnics,numerous problems involve large deformation,such as the installation of foundations(Jin YF et al.,2018a),landslides(Jin YF et al.,2020b),debris flow(Dai et al.,2017),collapse of undergro...For geohazards and geotechnics,numerous problems involve large deformation,such as the installation of foundations(Jin YF et al.,2018a),landslides(Jin YF et al.,2020b),debris flow(Dai et al.,2017),collapse of underground structures(Zhang et al.,2019),and the formation of sinkholes(Baran-diaran Villegas,2018).Benefitting from the sustained development of computing power,numerical simulations have become useful analytical methods in geomechanics and related fields.展开更多
基金part of the TPS projecta Vied-Newton PhD scholarship+1 种基金a Dixon scholarship from Imperial College London,UKthe Dean’s Fund from Imperial College London for financial support(2017-2020)。
文摘The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.
基金This is part of the TPS projecta Vied-Newton PhD scholarship and a Dixon scholarship from Imperial College London, UK, for supporting his studies at Imperial College Londonthe Dean’s Fund from Imperial College London for financial support (2017-2020).
文摘A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.
基金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.
基金The financial support from Xi'an Jiaotong-Liverpool University(grant No.RDF 18-01-23,PGRS 1906002 and REF-20-01-01)is gratefully acknowledgedThe first and second authors would also like to appreciate the support by the Key Program Special Fund at Xi'an Jiaotong-Liverpool University(grant No.KSF-E-19).
文摘This paper investigates the strength and critical state(CS)behaviours of granular materials via DEM simulations of true triaxial drained tests under three different loading modes including constant b(in-termediate stress ratio)tests with constantσ'1(major principal stress),constant p(mean pressure)and constantσ'3(minor principal stress)respectively.To this end,a series of samples are generated with the same particle size distribution,and with the confining stresses ranging from 100 kPa to 900 kPa.The CS is achieved for all samples.Both the macroscopic behaviours and the microscopic behaviours are examined and compared considering different loading modes,confining stresses and intermediate stress ratios(b).The critical state lines(CSLs)are found to be unique and independent of the loading modes,but dependent on the b values.The CSLs with b=0 and b=1 form the two boundaries of CSLs respectively beyond which CSLs under all other b tests cannot go beyond.Six different strength criteria are examined and compared in terms of both peak and CS failures.The Mohr-Coulomb strength criterion is found to be only suitable for axisymmetric loading conditions.The Lade-Duncan criterion is only suitable for describing peak strengths,which is dependent on the loading modes and confining stresses.The Satake criterion and Matsuoka-Nakai criterion are the more appropriate strength criteria for describing CS failures,indicating that the CS values of both the Satake parameter and the Matsuoka-Nakai parameter describe an inherent property that characterizes the CS failure for a given type of soil.The CS mechanical coordination number is fitted by a curved line for a given b value,which is unique regardless of the loading modes.The peak and CS values of both major and minor principal fabric tensors decrease with increasing b values,while they increase with increasing b values for the intermediate principal fabric tensor.
基金The authors would like to express their gratitude for the financial support from Xi’an Jiaotong-Liverpool University(XJTLU)(Grant Nos.RDF 15-01-38 and RDF 18-01-23)Also,the support by the Key Program Special Fund at XJTLU(Grant No.KSF-E-19)Natural Science Foundation of Jiangsu Province(Grant No.BK20160393)is greatly appreciated.
文摘This article studies the influences of particle morphology on the behaviors of granular materials at both macroscopic and microscopic levels based on the discrete element method(DEM).A set of numerical tests under drained triaxial compression was performed by controlling two morphological descriptors,i.e.ratio of the smallest to the largest pebble diameter,x,and the maximum pebbleepebble intersection angle,b.These descriptors are vital in generating particle geometry and surface textures.It was found that the stress responses of all assemblies exhibited similar behavior and showed post-peak strainsoftening.The normalized stress ratio and volumetric strains flatten off and tended to reach a steady value after an axial strain of 40%.While the friction angles at peak state varied with different morphological descriptors,the friction angles at critical state showed no significant variation.Moreover,evolution of the average coordination numbers showed a dramatic exponential decay until an axial strain of about 15%after which it stabilized and was unaffected by further increase of axial strain.In addition,stress ratio q/p and strong fabric parameter fs d=fs m were found to follow an approximately linear relationship for each assembly.These findings emphasized the significance of the influences of particle morphology on the macroscopic and microscopic responses of granular materials.
基金Projects (RDF18-01-23,PGRS1906002,REF-20-01-01) supported by Xi’an Jiaotong-Liverpool University,ChinaProject (KSF-E-19) supported by the Key Program Special Fund。
基金the Research Impact Fund(RIF)Project of Hong Kong Special Administrative Region Government of China(No.R5037-18)。
文摘For geohazards and geotechnics,numerous problems involve large deformation,such as the installation of foundations(Jin YF et al.,2018a),landslides(Jin YF et al.,2020b),debris flow(Dai et al.,2017),collapse of underground structures(Zhang et al.,2019),and the formation of sinkholes(Baran-diaran Villegas,2018).Benefitting from the sustained development of computing power,numerical simulations have become useful analytical methods in geomechanics and related fields.