Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.     

三种电离层延迟多频修正算法的比较

电离层延迟多频修正算法在修正电离层延迟的同时会放大观测噪声等伪距误差的影响。分析了电离层延迟双频修正、三频一阶修正和三频二阶修正三种修正算法及观测噪声的影响。以Galileo系统的四个载波频率为例,对不同频率组合下三种修正算法进行了比较。结果表明电离层延迟多频修正中,并不是观测量的频点数越多或修正掉的高阶项越多修正精度越高,还与观测伪距噪声的大小以及采用的修正方法有关。为多频测距系统的电离层延迟修正算法和最佳频率组合设计提供了一种可行的分析方法。     

A smoothed particle hydrodynamics modelling of soil–water mixing and resulting changes in average strength

Soil–water interaction is a pivotal process in many underwater geohazards such as underwater landslides where soil sediments gradually evolve into turbidity currents after interactions with ambient water. Due to the large deformations, multiphase interactions and phase changes this involves, investigations from numerical modelling of the transition process have been limited so far. This study explores a simple numerical replication of such soil–water mixing with respect to changes in average strength using smoothed particle hydrodynamics (SPH). A uniform viscoplastic model is used for both the solid-like and fluid-like SPH particles. The proposed numerical solution scheme is verified by single-phase dam break tests and multiphase simple shear tests. SPH combinations of solid-like and fluid-like particles can replicate the clay–water mixture as long as the liquidity index of the solid-like particles is larger than unity. The proposed numerical scheme is shown to capture key features of an underwater landslide such as hydroplaning, water entrainment and wave generation and thus shows promise as a tool to simulate the whole process of subaquatic geohazards involving solid–fluid transition during mass transport.     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

载波/多普勒平滑伪距在Android手机中的定位实现

针对Android手机GNSS伪距定位精度较低的问题,利用手机端观测信息,通过载波相位/多普勒平滑伪距改善手机端伪距观测值的质量,从而达到提高定位精度的效果。首先给出了Android手机GNSS原始观测量的获取方法,然后推导了载波相位平滑伪距和多普勒平滑伪距算法模型,并设计合理有效的试验对算法的精度进行评定。试验结果表明：在手机端静态定位中载波相位和多普勒平滑算法均可提高原始伪距的定位精度,且多普勒平滑算法表现更优;在手机端动态定位中多普勒平滑算法可获得比原始伪距更优的定位精度,但是载波相位平滑算法较原始伪距更差;由于硬件的制约手机端周跳和信号失锁严重,占比超过50%,载波相位在手机端的可用性较低;多普勒平滑算法的最优平滑时间常数小于等于10 s,具有实时动态定位的巨大潜力。     

����ԳƷ��ڵ��㶨λ�е�Ӧ��

???GPSα????λ????????????????????μ?????????????λ??????????????ù????????????????????????????????棨???棩?????????????????????????????????α?????????????????????α???????????????????????λ??????????·??????????????????λ???????????????????????????????ν????????????????????????Ч????????λ?????????????????     

An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: Case of hydrostatic pore‐water pressure

The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd.     

多频数据组合在周跳探测和修复上的应用

利用模拟的L5观测数据，结合传统的伪距／载波组合法，对周跳的探测与修复进行了实例计算，并与双频组合观测值的周跳探测方法进行了比较。结果证明，不管对大周跳还是小周跳，该方法均能在单历元间准确地探测出各个频率上发生的周跳，较之双频组合观测值的探测方法具有更大的优越性。     

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model

Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.