A difference scheme in curvilinear coordinates is put forward for calculation of salinity in estuaries and coastal waters, which is based on Eulerian-Lagrangian method. It combines first-order and second-order Lagrang...A difference scheme in curvilinear coordinates is put forward for calculation of salinity in estuaries and coastal waters, which is based on Eulerian-Lagrangian method. It combines first-order and second-order Lagrangian interpolation to reduce numerical dispersion and oscillation. And the length of the curvilinear grid is also considered in the interpolation. Then the scheme is used in estuary, coast and ocean model, and several numerical experiments for the Yangtze Estuary and the Hangzhou Bay are conducted to test it. These experiments show that it is suitable for simulations of salinity in estuaries and coastal waters with the models using curvilinear coordinates.展开更多
The simulation of large-strain geotechnical laboratory tests with conventional Lagrangian finite element method(FEM)techniques is often problematic due to excessive mesh distortion.The multiple reversal direct shear(M...The simulation of large-strain geotechnical laboratory tests with conventional Lagrangian finite element method(FEM)techniques is often problematic due to excessive mesh distortion.The multiple reversal direct shear(MRDS)test can be used to measure the residual shear strength of soils in a laboratory setting.However,modelling and simulation generally require advanced numerical methods to accommodate the large shear strains concentrated in the shear plane.In reality,when the standard direct shear(DS)apparatus is used,the MRDS method is prone to two major sources of measurement error:load cap tilting and specimen loss.These sources of error make it difficult or even impossible to correctly determine the residual shear strength.This paper presents a modified DS apparatus and multi-reversal multi-stage test method,simulated using the coupled Eulerian-Lagrangian(CEL)method in a finite element environment.The method was successful in evaluating equipment and preventing both load cap tilting and specimen loss,while modelling large-deformation behaviour that is not readily simulated with the conventional FEM or arbitrary Lagrangian-Eulerian(ALE)analysis.Thereafter,a modified DS apparatus was created for the purpose of analysing mixtures of organic materials found in an Australian clay.The results obtained from the modified DS CEL model in combination with laboratory tests show a great improvement in the measured residual shear strength profiles compared to those from the standard apparatus.The modified DS setup ensures that accurate material residual shear strengths are calculated,a factor that is vital to ensure appropriate soil behaviour is simulated for numerical analyses of large-scale geotechnical projects.展开更多
The expansion of a thick-walled hollow cylinder in soil is of non-self-similar nature that the stress/deformation paths are not the same for different soil material points.As a result,this problem cannot be solved by ...The expansion of a thick-walled hollow cylinder in soil is of non-self-similar nature that the stress/deformation paths are not the same for different soil material points.As a result,this problem cannot be solved by the common self-similar-based similarity techniques.This paper proposes a novel,exact solution for rigorous drained expansion analysis of a hollow cylinder of critical state soils.Considering stress-dependent elastic moduli of soils,new analytical stress and displacement solutions for the nonself-similar problem are developed taking the small strain assumption in the elastic zone.In the plastic zone,the cavity expansion response is formulated into a set of first-order partial differential equations(PDEs)with the combination use of Eulerian and Lagrangian descriptions,and a novel solution algorithm is developed to efficiently solve this complex boundary value problem.The solution is presented in a general form and thus can be useful for a wide range of soils.With the new solution,the non-self-similar nature induced by the finite outer boundary is clearly demonstrated and highlighted,which is found to be greatly different to the behaviour of cavity expansion in infinite soil mass.The present solution may serve as a benchmark for verifying the performance of advanced numerical techniques with critical state soil models and be used to capture the finite boundary effect for pressuremeter tests in small-sized calibration chambers.展开更多
针对海上设施的安全防护问题,提出了一种网-桁架式海上拦截装置,目的是在有效拦截来袭小艇撞击的同时尽量降低对拦截装置本身的损伤.为验证这一装置的有效性,基于显式动力学和欧拉-拉格朗日耦合方法对拦截装置拦截小艇过程进行数值模拟...针对海上设施的安全防护问题,提出了一种网-桁架式海上拦截装置,目的是在有效拦截来袭小艇撞击的同时尽量降低对拦截装置本身的损伤.为验证这一装置的有效性,基于显式动力学和欧拉-拉格朗日耦合方法对拦截装置拦截小艇过程进行数值模拟,设定小艇垂直撞击支撑柱、支撑柱间隙和45°撞击支撑柱、支撑柱间隙4种工况,每种工况包含10、20、30 m/s 3种速度.通过对这12组仿真的碰撞力、能量以及破损状态分析,全方位分析了该装置的防撞性能和对小艇的拦截效果.所有工况下该装置都能完成对小艇的拦截,说明这种装置具有优秀的拦截效果.展开更多
This paper provides an analysis on the effects of exact and inexact integrations on stability, convergence, numerical diffusion, and numerical oscillations for the Eulerian- Lagrangian method （ELM）. In the finite el...This paper provides an analysis on the effects of exact and inexact integrations on stability, convergence, numerical diffusion, and numerical oscillations for the Eulerian- Lagrangian method （ELM）. In the finite element ELM, when more accurate integrations are used for the right-hand-side, less numerical diffusion is introduced and better approximation is obtained. When linear interpolation is used for numerical integrations, the resulting ELM is shown to be unconditionally stable and of first-order accuracy. When Gauss quadrature is used, conditional stability and second-order accuracy are established under some mild constraints for the convection-diffusion problems. Finally, numerical experiments demonstrate that more accurate integrations lead to better approximation, and spatial adaptivity can substantially reduce numerical oscillations and smearing that often occur in the ELM when inexact numerical integrations are used.展开更多
This paper is devoted to the study of the Eulerian-Lagrangian method(ELM)for convection-diffusion equations on unstructured grids with or without accurate numerical integration.We first propose an efficient and accura...This paper is devoted to the study of the Eulerian-Lagrangian method(ELM)for convection-diffusion equations on unstructured grids with or without accurate numerical integration.We first propose an efficient and accurate algorithm to calculate the integrals in the Eulerian-Lagrangian method.Our approach is based on an algorithm for finding the intersection of two non-matching grids.It has optimal algorithmic complexity and runs fast enough to make time-dependent velocity fields feasible.The evaluation of the integrals leads to increased precision and the unconditional stability.We demonstrate by numerical examples that the ELM with our proposed algorithm for accurate numerical integration has the following two features:first it is much more accurate and more stable than the ones with traditional numerical integration techniques and secondly the overall cost of the proposed method is comparable with the traditional ones.展开更多
基金This project was supported by the Major State Basic Research Program under Contract Grant No. G1999043803the University Fund for Mainstay Teachers of State Ministry of Education and the Opening Fund of Open Laboratory of Marine Dynamic Process and Sa
文摘A difference scheme in curvilinear coordinates is put forward for calculation of salinity in estuaries and coastal waters, which is based on Eulerian-Lagrangian method. It combines first-order and second-order Lagrangian interpolation to reduce numerical dispersion and oscillation. And the length of the curvilinear grid is also considered in the interpolation. Then the scheme is used in estuary, coast and ocean model, and several numerical experiments for the Yangtze Estuary and the Hangzhou Bay are conducted to test it. These experiments show that it is suitable for simulations of salinity in estuaries and coastal waters with the models using curvilinear coordinates.
文摘The simulation of large-strain geotechnical laboratory tests with conventional Lagrangian finite element method(FEM)techniques is often problematic due to excessive mesh distortion.The multiple reversal direct shear(MRDS)test can be used to measure the residual shear strength of soils in a laboratory setting.However,modelling and simulation generally require advanced numerical methods to accommodate the large shear strains concentrated in the shear plane.In reality,when the standard direct shear(DS)apparatus is used,the MRDS method is prone to two major sources of measurement error:load cap tilting and specimen loss.These sources of error make it difficult or even impossible to correctly determine the residual shear strength.This paper presents a modified DS apparatus and multi-reversal multi-stage test method,simulated using the coupled Eulerian-Lagrangian(CEL)method in a finite element environment.The method was successful in evaluating equipment and preventing both load cap tilting and specimen loss,while modelling large-deformation behaviour that is not readily simulated with the conventional FEM or arbitrary Lagrangian-Eulerian(ALE)analysis.Thereafter,a modified DS apparatus was created for the purpose of analysing mixtures of organic materials found in an Australian clay.The results obtained from the modified DS CEL model in combination with laboratory tests show a great improvement in the measured residual shear strength profiles compared to those from the standard apparatus.The modified DS setup ensures that accurate material residual shear strengths are calculated,a factor that is vital to ensure appropriate soil behaviour is simulated for numerical analyses of large-scale geotechnical projects.
基金funding support from the National Key Research and Development Program of China(Grant No.2023YFB2604004)the National Natural Science Foundation of China(Grant No.52108374)the“Taishan”Scholar Program of Shandong Province,China(Grant No.tsqn201909016)。
文摘The expansion of a thick-walled hollow cylinder in soil is of non-self-similar nature that the stress/deformation paths are not the same for different soil material points.As a result,this problem cannot be solved by the common self-similar-based similarity techniques.This paper proposes a novel,exact solution for rigorous drained expansion analysis of a hollow cylinder of critical state soils.Considering stress-dependent elastic moduli of soils,new analytical stress and displacement solutions for the nonself-similar problem are developed taking the small strain assumption in the elastic zone.In the plastic zone,the cavity expansion response is formulated into a set of first-order partial differential equations(PDEs)with the combination use of Eulerian and Lagrangian descriptions,and a novel solution algorithm is developed to efficiently solve this complex boundary value problem.The solution is presented in a general form and thus can be useful for a wide range of soils.With the new solution,the non-self-similar nature induced by the finite outer boundary is clearly demonstrated and highlighted,which is found to be greatly different to the behaviour of cavity expansion in infinite soil mass.The present solution may serve as a benchmark for verifying the performance of advanced numerical techniques with critical state soil models and be used to capture the finite boundary effect for pressuremeter tests in small-sized calibration chambers.
文摘针对海上设施的安全防护问题,提出了一种网-桁架式海上拦截装置,目的是在有效拦截来袭小艇撞击的同时尽量降低对拦截装置本身的损伤.为验证这一装置的有效性,基于显式动力学和欧拉-拉格朗日耦合方法对拦截装置拦截小艇过程进行数值模拟,设定小艇垂直撞击支撑柱、支撑柱间隙和45°撞击支撑柱、支撑柱间隙4种工况,每种工况包含10、20、30 m/s 3种速度.通过对这12组仿真的碰撞力、能量以及破损状态分析,全方位分析了该装置的防撞性能和对小艇的拦截效果.所有工况下该装置都能完成对小艇的拦截,说明这种装置具有优秀的拦截效果.
文摘This paper provides an analysis on the effects of exact and inexact integrations on stability, convergence, numerical diffusion, and numerical oscillations for the Eulerian- Lagrangian method （ELM）. In the finite element ELM, when more accurate integrations are used for the right-hand-side, less numerical diffusion is introduced and better approximation is obtained. When linear interpolation is used for numerical integrations, the resulting ELM is shown to be unconditionally stable and of first-order accuracy. When Gauss quadrature is used, conditional stability and second-order accuracy are established under some mild constraints for the convection-diffusion problems. Finally, numerical experiments demonstrate that more accurate integrations lead to better approximation, and spatial adaptivity can substantially reduce numerical oscillations and smearing that often occur in the ELM when inexact numerical integrations are used.
文摘This paper is devoted to the study of the Eulerian-Lagrangian method(ELM)for convection-diffusion equations on unstructured grids with or without accurate numerical integration.We first propose an efficient and accurate algorithm to calculate the integrals in the Eulerian-Lagrangian method.Our approach is based on an algorithm for finding the intersection of two non-matching grids.It has optimal algorithmic complexity and runs fast enough to make time-dependent velocity fields feasible.The evaluation of the integrals leads to increased precision and the unconditional stability.We demonstrate by numerical examples that the ELM with our proposed algorithm for accurate numerical integration has the following two features:first it is much more accurate and more stable than the ones with traditional numerical integration techniques and secondly the overall cost of the proposed method is comparable with the traditional ones.
