A Novel Meshfree Analysis of Transient Heat Conduction Problems Using RRKPM

By introducing the radial basis functions(RBFs)into the reproducing kernel particle method(RKPM),the calculating accuracy and stability of the RKPM can be improved,and a novel meshfree method of the radial basis RKPM(meshfree RRKPM)is proposed.Meanwhile,the meshfree RRKPM is applied to transient heat conduction problems(THCP),and the corresponding equations of the meshfree RRKPM for the THCP are derived.The twopoint time difference scheme is selected to discretize the time of the THCP.Finally,the numerical results illustrate the effectiveness of the meshfree RRKPM for the THCP.

Precision of meshfree methods and application to forward modeling of two-dimensional electromagnetic sources

Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method （EFGM）, the point interpolation method （PIM）, and the radial point interpolation method （RPIM）. Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.

A Size-Dependent Functionally Graded Higher Order Plate Analysis Based on Modified Couple Stress Theory and Moving Kriging Meshfree Method

A size-dependent computational approach for bending,free vibration and buckling analyses of isotropic and sandwich functionally graded(FG)microplates is in this study presented.We consider both shear deformation and small scale effects through the generalized higher order shear deformation theory and modified couple stress theory(MCST).The present model only retains a single material length scale parameter for capturing properly size effects.A rule of mixture is used to model material properties varying through the thickness of plates.The principle of virtual work is used to derive the discrete system equations which are approximated by moving Kriging interpolation(MKI)meshfree method.Numerical examples consider the inclusions of geometrical parameters,volume fraction,boundary conditions and material length scale parameter.Reliability and effectiveness of the present method are confirmed through numerical results.

Computational Study of Collective Cell Migration By Meshfree Method

The collective cell migration behavior on a substrate was studied using RKPM meshfree method.The cells were modeled as nematic liquid crystal with hyperelastic cell nucleus.The cell-substrate and cell-cell interactions were modeled by coarse-grained potential forces.Through this study,the pulling and pushing phenomenon during collective cell migration process was observed and it was found that the individual cell mobility significantly influenced the collective cell migratory behavior.More self-propelled cells are in the system along the same direction,the faster the collective group migrates toward coordinated direction.The parametric study on cell-cell adhesion strength indicated that as the adhesion strength increases,the collective cell migration speed increases.It also showed that the mechanical stress in leader cell is higher than stress in follower cells.

A meshfree-based local Galerkin method with condensation of degree of freedom for elastic dynamic analysis

Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.

An efficient Galerkin meshfree formulation for shear deformable beam under finite deformation

This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear analysis encompasses the fully geometric nonlinearities due to large deflection,large deformation as well as finite rotation.The incremental equilibrium equation is obtained by the consistent linearization of the nonlinear variational equation.The Lagrangian meshfree shape function is utilized to discretize the variational equation.Subsequently to resolve the shear and membrane locking issues and accelerate the computation,the method of stabilized conforming nodal integration is systematically implemented through the Lagrangian gradient smoothing operation.Numerical results reveal that the present formulation is very effective.

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.

基于MIDAS MeshFree的门式起重机结构分析

介绍了隐式边界法基本原理和MIDAS MeshFree软件特点,以480 t-146 m门式起重机为建模对象,在不进行模型简化和网格划分的条件下,对其进行强度和刚度计算,得到在自重工况下门式起重机主梁的挠度以及工况501下主梁跨中最大Von-mises应力,以判断门式起重机结构的合理性。MIDAS MeshFree仿真软件计算过程简单、高效,有助于起重机的设计和结构改进,可缩短设计周期。

基于MeshFree的铁路货车转K6型摇枕仿真计算分析

介绍了隐式边界法基本原理和MIDAS MeshFree软件特点.摇枕是铁路货车关键部件之一,其仿真计算的准确性对于疲劳寿命分析具有很大的影响.文章以转K6摇枕为建模对象,在不进行模型简化和网格划分的条件下,建立基于摇枕静强度试验台架的仿真计算模型,在试验件与工装接触处设置多个不同的摩擦系数进行计算分析,基于摇枕静强度试验测试结果,分析了接触条件对测点计算应力的影响.结果表明,在转K6摇枕结构的仿真计算模型中,接触处摩擦系数的影响明显,其最佳摩擦系数为0.15.

A meshfree method and its applications to elasto-plastic problems

Standard finite element approaches are still ineffective in handling extreme material deformation, such as cases of large deformations and moving discontinuities due to severe mesh distortion. Among meshfree methods developed to overcome the ineffectiveness, Reproducing Kernel Particle Method (RKPM) has demonstrated its great suitability for structural analysis.This paper presents applications of RKPM in elasto-plastic problems after a review of meshfree methods and an introduction to RKPM. A slope stability problem in geotechnical engineering is analyzed as an illustrative case. The corresponding numerical simulations are carried out on an SGI Onyx3900 supercomputer. Comparison of the RKPM and the FEM under identical conditions showed that the RKPM is more suitable for problems where there exists extremely large strain such as in the case of slope sliding.

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.

Factors affecting accuracy of radial point interpolation meshfree method for 3-D solid mechanics

Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.

Three-Dimensional Meshfree Analysis of Interlocking Concrete Blocks for Step Seawall Structure

This study adapts the flexible characteristic of meshfree method in analyzing three-dimensional(3D)complex geometry structures,which are the interlocking concrete blocks of step seawall.The elastostatic behavior of the block is analysed by solving the Galerkin weak form formulation over local support domain.The 3D moving least square(MLS)approximation is applied to build the interpolation functions of unknowns.The pre-defined number of nodes in an integration domain ranging from 10 to 60 nodes is also investigated for their effect on the studied results.The accuracy and efficiency of the studied method on 3D elastostatic responses are validated through the comparison with the solutions of standard finite element method(FEM)using linear shape functions on tetrahedral elements and the well-known commercial software,ANSYS.The results show that elastostatic responses of studied concrete block obtained by meshfree method converge faster and are more accurate than those of standard FEM.The studied meshfree method is effective in the analysis of static responses of complex geometry structures.The amount of discretised nodes within the integration domain used in building MLS shape functions should be in the range from 30 to 60 nodes and should not be less than 20 nodes.

ICM Method Combined with Meshfree Approximation for Continuum Structure

The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of the essential boundary condition and derivative of various sensitivities,a singular weight function in element free Galerkin method is introduced.Material point variable is defined to illustrate the condition of material point and its vicinity instead of element or node.The topological variables field is constructed by moving least square approximation which inherits the continuity and smoothness of the weight function.Due to reciprocal relationships between the topological variables and design variables,various structural responses sensitivities are derived according to the method for calculating the partial derivatives of compound functions.Numerical examples indicate that checkerboard pattern and mesh-dependence phenomena are overcome without additional restriction methods.

Meshfree First-order System Least Squares

We prove convergence for a meshfree first-order system least squares(FOSLS) partition of unity finite element method(PUFEM).Essentially,by virtue of the partition of unity,local approximation gives rise to global approximation in H(div)∩H(curl). The FOSLS formulation yields local a posteriori error estimates to guide the judicious allotment of new degrees of freedom to enrich the initial point set in a meshfree dis- cretization.Preliminary numerical results are provided and remaining challenges are discussed.

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of this paper.Firstly,we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space,which is competitive for high-dimensional random inputs.Secondly,the a priori error estimates are derived for the state,the co-state and the control variables.Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.

An Updated Lagrangian Particle Hydrodynamics (ULPH)-NOSBPD Coupling Approach forModeling Fluid-Structure Interaction Problem

A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.

A RELAXED HSS PRECONDITIONER FOR SADDLE POINT PROBLEMS FROM MESHFREE DISCRETIZATION＊

In this paper, a relaxed Hermitian and skew-Hermitian splitting （RHSS） preconditioner is proposed for saddle point problems from the element-free Galerkin （EFG） discretization method. The EFG method is one of the most widely used meshfree methods for solving partial differential equations. The RHSS preconditioner is constructed much closer to the coefficient matrix than the well-known HSS preconditioner, resulting in a RHSS fixed-point iteration. Convergence of the RHSS iteration is analyzed and an optimal parameter, which minimizes the spectral radius of the iteration matrix is described. Using the RHSS pre- conditioner to accelerate the convergence of some Krylov subspace methods （like GMRES） is also studied. Theoretical analyses show that the eigenvalues of the RHSS precondi- tioned matrix are real and located in a positive interval. Eigenvector distribution and an upper bound of the degree of the minimal polynomial of the preconditioned matrix are obtained. A practical parameter is suggested in implementing the RHSS preconditioner. Finally, some numerical experiments are illustrated to show the effectiveness of the new preconditioner.

Three dimensional efficient meshfree simulation of large deformation failure evolution in soil medium

An efficient Galerkin meshfree formulation for three dimensional simulation of large deformation failure evolution in soils is presented. This formulation utilizes the stabilized conforming nodal integration, where for the purpose of stability and efficiency a Lagrangian smoothing strain at nodal point is constructed and thereafter the internal energy is evaluated nodally. This formulation ensures the linear exactness, efficiency and spatial stability in a unified manner and it makes the conventional Galerkin meshfree method affordable for three dimensional simulation. The three dimensional implementation of stabilized conforming nodal integration is discussed in details. To model the failure evolution in soil medium a coupled elasto-plastic damage model is used and an objective stress integration algorithm in combination of elasto-damage predictor and plastic corrector method is employed for stress update. Two typical numerical examples are shown to demonstrate the effectiveness of the present method for modeling large deformation soil failure.

A three-dimensional two-level gradient smoothing meshfree method for rainfall induced landslide simulations

A three-dimensional two-level gradient smoothing meshfree method is presented for rainfall induced landslide simulations. The two-level gradient smoothing for meshfree shape function is elaborated in the threedimensional Lagrangian setting with detailed implementation procedure. It is shown that due to the successive gradient smoothing operation without the requirement of derivative computation in the present formulation, the two-level smoothed gradient of meshfree shape function is capable of achieving a given influence domain more efficiently than the standard gradient of meshfree shape function. Subsequently, the two-level smoothed gradient of meshfree shape function is employed to discretize the weak form of coupled rainfall seepage and soil motion equations in a nodal integration format, as provides an efficient three-dimensional regularized meshfree formulation for large deformation rainfall induced landslide simulations. The exponential damage and pressure dependent plasticity relationships are utilized to describe the failure evolution in landslides. The plastic response of soil is characterized by the true effective stress measure, which is updated according to the rotationally neutralized objective integration algorithm. The effectiveness of the present threedimensional two?level gradient smoothing meshfree method is demonstrated through numerical examples.