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.

Mobility and dynamic erosion process of granular flow:insights from numerical investigation using material point method

In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.

Low-order Wavefront Error Compensation for Multi-field of Lithography Projection Objective Based on Interior Point Method

Low-order wavefront error account for a large proportion of wave aberrations.A compensation method for low order aberration of projection lithography objective based on Interior Point Method is presented.Compensation model between wavefront error and degree of movable lens freedom is established.Converting over-determined system to underdetermined system,the compensation is solved by Interior Point Method(IPM).The presented method is compared with direct solve the over-determined system.Then,other algorithm GA,EA and PS is compared with IPM.Simulation and experimental results show that the presented compensation method can obtained compensation with less residuals compared with direct solve the over-determined system.Also,the presented compensation method can reduce computation time and obtain results with less residuals compare with AGA,EA and PS.Moreover,after compensation,RMS of wavefront error of the experimental lithography projection objective decrease from 56.05 nm to 17.88 nm.

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.

AHermitian C^(2) Differential Reproducing Kernel Interpolation Meshless Method for the 3D Microstructure-Dependent Static Flexural Analysis of Simply Supported and Functionally Graded Microplates

This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.

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.

Prediction of Flash Point Temperature of Organic Compounds Using a Hybrid Method of Group Contribution ＋ Neural Network ＋ Particle Swarm Optimization

The flash points of organic compounds were estimated using a hybrid method that includes a simple group contribution method (GCM) implemented in an artificial neural network (ANN) with particle swarm optimization (PSO). Different topologies of a multilayer neural network were studied and the optimum architecture was determined. Property data of 350 compounds were used for training the network. To discriminate different substances the molecular structures defined by the concept of the classical group contribution method were given as input variables. The capabilities of the network were tested with 155 substances not considered in the training step. The study shows that the proposed GCM+ANN+PSO method represent an excellent alternative for the estimation of flash points of organic compounds with acceptable accuracy (AARD = 1.8%; AAE = 6.2 K).

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.

VISCOSITY ITERATIVE METHODS FOR COMMON FIXED POINTS OF TWO NONEXPANSIVE MAPPINGS WITHOUT COMMUTATIVITY ASSUMPTION IN HILBERT SPACES

In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of the mappings which solves some variational inequality. The results presented extend the corresponding results of Shimizu and Takahashi IT. Shimizu, W. Takahashi, Strong convergence to common fixed point of families of nonexpansive mappings, J. Math. Anal. Appl. 211 （1997）, 71-83], and Yao and Chen [Y. Yao, R. Chert, Convergence to common fixed points of average mappings without commutativity assumption in Hilbert spaces, Nonlinear Analysis 67（2007）, 1758-1763].

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.

Machine Learning Enhanced Boundary Element Method:Prediction of Gaussian Quadrature Points

This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.

New Modification of Fixed Point Iterative Method for Solving Nonlinear Equations

In this paper, we have modified fixed point method and have established two new iterative methods of order two and three. We have discussed their convergence analysis and comparison with some other existing iterative methods for solving nonlinear equations.

Element-free Galerkin method for free vibration of rectangular plates with interior elastic point supports and elastically restrained edges

The element-free Galerkin method is proposed to solve free vibration of rectangular plates with finite interior elastic point supports and elastically restrained edges.Based on the extended Hamilton＇s principle for the elastic dynamics system,the dimensionless equations of motion of rectangular plates with finite interior elastic point supports and the edge elastically restrained are established using the element-free Galerkin method.Through numerical calculation,curves of the natural frequency of thin plates with three edges simply supported and one edge elastically restrained,and three edges clamped and the other edge elastically restrained versus the spring constant,locations of elastic point support and the elastic stiffness of edge elastically restrained are obtained.Effects of elastic point supports and edge elastically restrained on the free vibration characteristics of the thin plates are analyzed.

An ideal point method adopting grey relational entropy to estimate life cycle cost with effectiveness of torpedo

To making the decision of the developing blue prints,ideal point method was selected to estimate the life cycle cost with effectiveness of torpedo.At the same time,the concept of grey relational entropy of the grey system theory was adopted to compute the distance between each blue print and the ideal point(or negative ideal point).The blue print,nearest to the ideal point and farthest to the negative ideal point,is the best one.As an example,four blue prints of torpedo were estimated.The result indicates the practical value of this method.

“Three Methods and Three Points” regulates p38 mitogen-activated protein kinase in the dorsal horn of the spinal cord in a rat model of sciatic nerve injury

Tuina is a traditional Chinese treatment for sensory disturbances caused by peripheral nerve injury and related diseases. Our previous studies showed that tuina regulates relevant regions and indices of the spinal dorsal horn using the Dian, Bo, and Rou method in Yinmen（BL37）, Yanglingquan（GB34）, and Weizhong（BL40）. Treatment prevents muscle atrophy, protects spinal cord neurons, and promotes sciatic nerve repair. The mechanisms of action of tuina for treating peripheral nerve injury remain poorly understood. This study established rat models of sciatic nerve injury using the crushing method. Rats received Chinese tuina in accordance with the principle of ＂Three Methods and Three Points,＂ once daily for 20 days. Tuina intervention reduced paw withdrawal latency and improved wet weight of the gastrocnemius muscle, as well as promoting morphological recovery of sciatic nerve fibers, Schwann cells, and axons. The protein expression levels of phospho-p38 mitogen-activated protein kinase, tumor necrosis factor-α, and interleukin-1β also decreased. These findings indicate that ＂Three Methods and Three Points＂ promoted morphological recovery and improved behavior of rats with peripheral nerve injury.

Determination of work-hardening exponent for sheet metal using three-point method

Three characteristic points in the deformation history of a fractured tensile specimen are selected tocalculate two values of n( n1 and n2 ) to represent the hardening ability of material during the homogeneous plas-tic deformation and the following large plastic deformation. Experimental results obtained with mild streel andred copper show that n determined using the three-point method proposed is better to describe the hardening a-bility of material. It is therefore concluded that three-point method can be used to describe the hardening prop-erty of material during both homogeneous deformation and large plastic deformation.

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.

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.

Exploring the characteristics of acupoints in the treatment of stroke with complex network and point mutual information method

Objective:Explore the characteristics of acupoints in the treatment of stroke with complex network and point mutual information method.Methods:The complex network and point wise mutual information system-developed by Chinese academy of Chinese medical sciences wereused to analyze the specific acupoints,compatibility,frequency etc.Results:174 acumoxibustion prescriptions were collected,including 163 acupoints.among them eighteen acupoints were used more than 30 times such as Hegu(LI4),Zusanli(ST36),Quchi(LI11)and Fengshi(GB31).The combinations of 31 acupoints were used more than 15 times,such as the combination of Quchi(LI11)and Zusanli(ST36),the combination of acupoint Quchi(LI11)and Jianyu(LI15),Hegu and Quchi(LI11).The most commonly used treatment method for stroke treatment is to dredge the Yangming meridian and Shaoyang meridian through acupuncture the multiple acupoints located on these two meridians..The commonly used acupoints are mainly distributed in the limbs,head and face.The most commonly used specific acupoint is intersection acupoint.The usage frequency of specific acupoints are higher than that of non-specific acupoints.Conclusion:Dredging the collaterals,dispelling wind-evil and restoring consciousness are the main principle for the treatment of stroke.Specific acupoints in head,face and climbs maybe the main targeted acupoints.Combination of Yang meridians with other meridians is needed to improve the effects.The Yangming meridian and Shaoyang meridian are most used meridians and Hegu(LI4),Quchi(LI11)and Zusanli(ST36)are the most used acupionts.

Estimation of non-point source pollution loads with flux method in Danjiangkou Reservoir area,China

The estimation of non-point source pollution loads into the Danjiangkou Reservoir is highly significant to environmental protection in the watershed. In order to overcome the drawbacks of traditional watershed numerical models, a base flow separation method was established coupled with a digital filtering method and a flux method. The digital filtering method has been used to separate the base flows of the Hanjiang,Tianhe, Duhe, Danjiang, Laoguan, and Qihe rivers. Based on daily discharge, base flow, and pollutant concentration data, the flux method was used to calculate the point source pollution load and non-point source pollution load. The results show that:(1) In the year 2013, the total inflow of the six rivers mentioned above accounted for 95.9% of the total inflow to the Danjiangkou Reservoir. The total pollution loads of chemical oxygen demand(CODMn) and total phosphorus(TP) from the six rivers were 58.20 103 t and 1.863 10~3 t, respectively, and the non-point source pollution loads were 39.82 10~3 t and 1.544 10~3 t, respectively, indicating that the non-point source pollution is a major factor(with a contribution rate of 68.4% for CODMnand 82.9% for TP).(2) The Hanjiang River is the most significant contributor of pollution loads to the Danjiangkou Reservoir, and its CODMnand TP contribution rates reached 79.3% and 83.2%, respectively. The Duhe River took the second place.(3) Non-point source pollution mainly occurred in the wet season in 2013, accounting for 80.8% and 90.9% of the total pollution loads of CODMnand TP, respectively. It is concluded that the emphasis of pollution control should be placed on non-point source pollution.