A hybrid contact approach for modeling soil-structure interaction using the material point method

The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early contact and contact penetration can occur when the contact conditions are unsuitable,and(2)the method is not available for contact problems involving rigid-nonrigid materials,which can cause numerical instability.This study presents a new hybrid contact approach for the MPM to address these limitations to simulate the soil and structure interactions.The approach combines the advantages of point-point and point-segment contacts to implement contact detection,satisfying the impenetrability condition and smoothing the corner contact problem.The proposed approach is first validated through a disk test on an inclined slope.Then,several typical cases,such as granular collapse,bearing capacity,and deformation of a flexible retaining wall,are simulated to demonstrate the robustness of the proposed approach compared with FEM or analytical solutions.Finally,the proposed method is used to simulate the impact of sand flow on a deformable structure.The results show that the proposed contact approach can well describe the phenomenon of soil-structure interaction problems.

Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.

An Efficient Reliability-Based Optimization Method Utilizing High-Dimensional Model Representation and Weight-Point Estimation Method

The objective of reliability-based design optimization(RBDO)is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the efficiency of RBDO algorithm,which hinders their application to high-dimensional engineering problems.To address these issues,this paper proposes an efficient decoupled RBDO method combining high dimensional model representation(HDMR)and the weight-point estimation method(WPEM).First,we decouple the RBDO model using HDMR and WPEM.Second,Lagrange interpolation is used to approximate a univariate function.Finally,based on the results of the first two steps,the original nested loop reliability optimization model is completely transformed into a deterministic design optimization model that can be solved by a series of mature constrained optimization methods without any additional calculations.Two numerical examples of a planar 10-bar structure and an aviation hydraulic piping system with 28 design variables are analyzed to illustrate the performance and practicability of the proposed method.

Modeling footing resting on anisotropic sand using material point method

Sand typically exhibits anisotropic internal structure which may significantly influence its mechanical behavior. The material point method (MPM) can eliminate mesh distortion and thus is suitable for investigating geotechnical problems with large deformation. In this study, an advanced anisotropic critical state theory (ACST)-based soil model is implemented in MPM to study the response of strip footing resting on anisotropic sand. The capability of the model is verified by simulating several element tests and strip footing tests with different soil densities and fabric bedding plane orientations. For the footing problem with a vertical load, as the fabric bedding plane orientation increases, the bearing capacity decreases and its corresponding settlement increases. The failure pattern becomes asymmetrical when the bedding plane orientation or the loading direction is inclined. A comparison between the simulation results predicted by the anisotropic and isotropic models is made, which demonstrates that neglecting the fabric anisotropy may lead to the overestimation of the bearing capacity.

Modeling pipe-soil interaction under vertical downward relative offset using B-spline material point method

To analyze the pipeline response under permanent ground deformation,the evolution of resistance acting on the pipe during the vertical downward offset is an essential ingredient.However,the efficient simulation of pipe penetration into soil is challenging for the conventional finite element(FE)method due to the large deformation of the surrounding soils.In this study,the B-spline material point method(MPM)is employed to investigate the pipe-soil interaction during the downward movement of rigid pipes buried in medium and dense sand.To describe the density-and stress-dependent behaviors of sand,the J2-deformation type model with state-dependent dilatancy is adopted.The effectiveness of the model is demonstrated by element tests and biaxial compression tests.Afterwards,the pipe penetration process is simulated,and the numerical outcomes are compared with the physical model tests.The effects of pipe size and burial depth are investigated with an emphasis on the mobilization of the soil resistance and the failure mechanisms.The simulation results indicate that the bearing capacity formulas given in the guidelines can provide essentially reasonable estimates for the ultimate force acting on buried pipes,and the recommended value of yield displacement may be underestimated to a certain extent.

Coseismic site response and slope instability using periodic boundary conditions in the material point method

This paper proposed the explicit generalized-a time scheme and periodic boundary conditions in the material point method(MPM)for the simulation of coseismic site response.The proposed boundary condition uses an intuitive particle-relocation algorithm ensuring material points always remain within the computational mesh.The explicit generalized-a time scheme was implemented in MPM to enable the damping of spurious high frequency oscillations.Firstly,the MPM was verified against finite element method(FEM).Secondly,ability of the MPM in capturing the analytical transfer function was investigated.Thirdly,a symmetric embankment was adopted to investigate the effects of ground motion arias intensity(I_(a)),geometry dimensions,and constitutive models.The results show that the larger the model size,the higher the crest runout and settlement for the same ground motion.When using a Mohr-Coulomb model,the crest runout increases with increasing I_(a).However,if the strain-softening law is activated,the results are less influenced by the ground motion.Finally,the MPM results were compared with the Newmark sliding block solution.The simplified analysis herein highlights the capabilities of MPM to capture the full deformation process for earthquake engineering applications,the importance of geometry characterization,and the selection of appropriate constitutive models when simulating coseismic site response and subsequent large deformations.

A reasoning diagram based method for fault diagnosis of railway point system

Railway Point System(RPS)is an important infrastructure in railway industry and its faults may have significant impacts on the safety and efficiency of train operations.For the fault diagnosis of RPS,most existing methods assume that sufficient samples of each failure mode are available,which may be unrealistic,especially for those modes of low occurrence frequency but with high risk.To address this issue,this work proposes a novel fault diagnosis method that only requires the power signals generated under normal RPS operations in the training stage.Specifically,the failure modes of RPS are distinguished through constructing a reasoning diagram,whose nodes are either binary logic problems or those that can be decomposed into the problems of the binary logic.Then,an unsupervised method for the signal segmentation and a fault detection method are combined to make decisions for each binary logic problem.Based on the results of decisions,the diagnostic rules are established to identify the failure modes.Finally,the data collected from multiple real-world RPSs are used for validation and the results demonstrate that the proposed method outperforms the benchmark in identifying the faults of RPSs.

A Disturbance Localization Method for Power System Based on Group Sparse Representation and Entropy Weight Method

This paper addresses the problem of complex and challenging disturbance localization in the current power system operation environment by proposing a disturbance localization method for power systems based on group sparse representation and entropy weight method.Three different electrical quantities are selected as observations in the compressed sensing algorithm.The entropy weighting method is employed to calculate the weights of different observations based on their relative disturbance levels.Subsequently,by leveraging the topological information of the power system and pre-designing an overcomplete dictionary of disturbances based on the corresponding system parameter variations caused by disturbances,an improved Joint Generalized Orthogonal Matching Pursuit(J-GOMP)algorithm is utilized for reconstruction.The reconstructed sparse vectors are divided into three parts.If at least two parts have consistent node identifiers,the node is identified as the disturbance node.If the node identifiers in all three parts are inconsistent,further analysis is conducted considering the weights to determine the disturbance node.Simulation results based on the IEEE 39-bus system model demonstrate that the proposed method,utilizing electrical quantity information from only 8 measurement points,effectively locates disturbance positions and is applicable to various disturbance types with strong noise resistance.

Multi-Branch Fault Line Location Method Based on Time Difference Matrix Fitting

The distribution network exhibits complex structural characteristics,which makes fault localization a challenging task.Especially when a branch of the multi-branch distribution network fails,the traditional multi-branch fault location algorithm makes it difficult to meet the demands of high-precision fault localization in the multi-branch distribution network system.In this paper,the multi-branch mainline is decomposed into single branch lines,transforming the complex multi-branch fault location problem into a double-ended fault location problem.Based on the different transmission characteristics of the fault-traveling wave in fault lines and non-fault lines,the endpoint reference time difference matrix S and the fault time difference matrix G were established.The time variation rule of the fault-traveling wave arriving at each endpoint before and after a fault was comprehensively utilized.To realize the fault segment location,the least square method was introduced.It was used to find the first-order fitting relation that satisfies the matching relationship between the corresponding row vector and the first-order function in the two matrices,to realize the fault segment location.Then,the time difference matrix is used to determine the traveling wave velocity,which,combined with the double-ended traveling wave location,enables accurate fault location.

A Visual Indoor Localization Method Based on Efficient Image Retrieval

The task of indoor visual localization, utilizing camera visual information for user pose calculation, was a core component of Augmented Reality (AR) and Simultaneous Localization and Mapping (SLAM). Existing indoor localization technologies generally used scene-specific 3D representations or were trained on specific datasets, making it challenging to balance accuracy and cost when applied to new scenes. Addressing this issue, this paper proposed a universal indoor visual localization method based on efficient image retrieval. Initially, a Multi-Layer Perceptron (MLP) was employed to aggregate features from intermediate layers of a convolutional neural network, obtaining a global representation of the image. This approach ensured accurate and rapid retrieval of reference images. Subsequently, a new mechanism using Random Sample Consensus (RANSAC) was designed to resolve relative pose ambiguity caused by the essential matrix decomposition based on the five-point method. Finally, the absolute pose of the queried user image was computed, thereby achieving indoor user pose estimation. The proposed indoor localization method was characterized by its simplicity, flexibility, and excellent cross-scene generalization. Experimental results demonstrated a positioning error of 0.09 m and 2.14° on the 7Scenes dataset, and 0.15 m and 6.37° on the 12Scenes dataset. These results convincingly illustrated the outstanding performance of the proposed indoor localization method.

Comparison between Two Kinds of Methods of Locating Acupoints

ＣｏｍｐａｒｉｓｏｎｂｅｔｗｅｅｎＴｗｏＫｉｎｄｓｏｆＭｅｔｈｏｄｓｏｆＬｏｃａｔｉｎｇＡｃｕｐｏｉｎｔｓ￥ＤｉｎｇＹｉｄａｎ，ＬｉＪｉｎｇｍｉｎｇ（ＭｅｄｉｃａｌＣｏｌｌｅｇｅｏｆＪｉｎａｎＵｎｉｖｅｒｓｉｔｙ）Ｔｈｅｐｒｏｐｏｒｔｉｏｎａｌｍｅａ...

Precise integration method for a class of singular two-point boundary value problems

In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.

A Combined Homotopy Infeasible Interior-Point Method for Convex Nonlinear Programming

In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method （CHIIP） for convex nonlinear programming problems. For any convex nonlinear programming, without strict convexity for the logarithmic barrier function, we get different solutions of the convex programming in different cases by CHIIP method.

Simulation of Seawater Intrusion in Coastal Confined Aquifer Using a Point Collocation Method Based Meshfree Model

Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.

Three Dimensional Laser Point Cloud Slicing Method for Calculating Irregular Volume

Volume parameter is the basic content of a spatial body object morphology analysis.However,the challenge lies in the volume calculation of irregular objects.The point cloud slicing method proposed in this study effectively works in calculating the volume of the point cloud of the spatial object obtained through three-dimensional laser scanning(3DLS).In this method,a uniformly spaced sequent slicing process is first conducted in a specific direction on the point cloud of the spatial object obtained through 3DLS.A series of discrete point cloud slices corresponding to the point cloud bodies are then obtained.Subsequently,the outline boundary polygon of the point cloud slicing is searched one by one in accordance with the slicing sequence and areas of the polygon.The point cloud slice is also calculated.Finally,the individual point cloud section volume is calculated through the slicing areas and the adjacent slicing gap.Thus,the total volume of the scanned spatial object can be calculated by summing up the individual volumes.According to the results and analysis of the calculated examples,the slice-based volume-calculating method for the point cloud of irregular objects obtained through 3DLS is correct,concise in process,reliable in results,efficient in calculation methods,and controllable on accuracy.This method comes as a good solution to the volume calculation of irregular objects.

Coupling between finite element method and material point method for problems with extreme deformation

As a Lagrangian meshless method, the material point method （MPM） is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method （FEM） for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with that of both FEM and MPM

MESHFREE METHOD BASED ON POINT COLLOCATION FOR METAL FORMING SIMULATION

A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM （finite element method） and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.

An improved local radial point interpolation method for transient heat conduction analysis

The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.

Modified Subgradient Extragradient Method for Variational Inequality Problems and Fixed Point Problems

Many approaches inquiring into variational inequality problems have been put forward,among which subgradient extragradient method is of great significance.A novel algorithm is presented in this article for resolving quasi-nonexpansive fixed point problem and pseudomonotone variational inequality problem in a real Hilbert interspace.In order to decrease the execution time and quicken the velocity of convergence,the proposed algorithm adopts an inertial technology.Moreover,the algorithm is by virtue of a non-monotonic step size rule to acquire strong convergence theorem without estimating the value of Lipschitz constant.Finally,numerical results on some problems authenticate that the algorithm has preferable efficiency than other algorithms.