NUMERICAL ANALYSIS OF MINDLIN SHELL BY MESHLESS LOCAL PETROV-GALERKIN METHOD

The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.

Quasi Ellipsoid Gear Surface Reconstruction Based on Meshless Local Petrov-Galerkin Method and Transmission Characteristic

Special transmission 3D model simulation must be based on surface discretization and reconstruction, but special transmission usually has complicated tooth shape and movement, so present software can＇t provide technical support for special transmission 3D model simulation. Currently, theoretical calculation and experimental method are difficult to exactly solve special transmission contact analysis problem. How to reduce calculation and computer memories consume and meet calculation precision is key to resolve special transmission contact analysis problem. According to 3D model simulation and surface reconstruction of quasi ellipsoid gear is difficulty, this paper employes meshless local Petrov-Galerkin （MLPG） method. In order to reduce calculation and computer memories consume, we disperse tooth mesh into finite points--sparseness points cloud or grid mesh, and then we do interpolation reconstruction in some necessary place of the 3D surface model during analysis. Moving least square method （MLSM） is employed for tooth mesh interpolation reconstruction, there are some advantages to do interpolation by means of MLSM, such as high precision, good flexibility and no require of tooth mesh discretization into units. We input the quasi ellipsoid gear reconstruction model into simulation software, we complete tooth meshing simulation. Simulation transmission ratio during meshing period was obtained, compared with theoretical transmission ratio, the result inosculate preferably. The method using curve reconstruction realizes surface reconstruction, reduce simulation calculation enormously, so special gears simulation can be realized by minitype computer. The method provides a novel solution for special transmission 3D model simulation analysis and contact analysis.

The complex variable meshless local Petrov-Galerkin method of solving two-dimensional potential problems

Based on the complex variable moving least-square（CVMLS） approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin（CVMLPG） method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square（MLS） approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin（MLPG） method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.

A MESHLESS LOCAL PETROV-GALERKIN METHOD FOR GEOMETRICALLY NONLINEAR PROBLEMS

Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.

Structural Reliability Assessment by a Modified Spectral Stochastic Meshless Local Petrov-Galerkin Method

This study presents a new tool for solving stochastic boundary-value problems. This tool is created by modify the previous spectral stochastic meshless local Petrov-Galerkin method using the MLPG5 scheme. This modified spectral stochastic meshless local Petrov-Galerkin method is selectively applied to predict the structural failure probability with the uncertainty in the spatial variability of mechanical properties. Except for the MLPG5 scheme, deriving the proposed spectral stochastic meshless local Petrov-Galerkin formulation adopts generalized polynomial chaos expansions of random mechanical properties. Predicting the structural failure probability is based on the first-order reliability method. Further comparing the spectral stochastic finite element-based and meshless local Petrov-Galerkin-based predicted structural failure probabilities indicates that the proposed spectral stochastic meshless local Petrov-Galerkin method predicts the more accurate structural failure probability than the spectral stochastic finite element method does. In addition, generating spectral stochastic meshless local Petrov-Galerkin results are considerably time-saving than generating Monte-Carlo simulation results does. In conclusion, the spectral stochastic meshless local Petrov-Galerkin method serves as a time-saving tool for solving stochastic boundary-value problems sufficiently accurately.

Weakly-Singular Traction and Displacement Boundary Integral Equations and Their Meshless Local Petrov-Galerkin Approaches

The general meshless local Petrov-Galerkin (MLPG) weak forms of the displacement and trac- tion boundary integral equations (BIEs) are presented for solids undergoing small deformations. Using the directly derived non-hyper-singular integral equations for displacement gradients, simple and straight- forward derivations of weakly singular traction BIEs for solids undergoing small deformations are also pre- sented. As a framework for meshless approaches, the MLPG weak forms provide the most general basis for the numerical solution of the non-hyper-singular displacement and traction BIEs. By employing the various types of test functions, several types of MLPG/BIEs are formulated. Numerical examples show that the pre- sent methods are very promising, especially for solving the elastic problems in which the singularities in dis- placements, strains, and stresses are of primary concern.

A Novel Localized Meshless Method for Solving Transient Heat Conduction Problems in Complicated Domains

This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.

Meshless Local Discontinuous Petrov-Galerkin Method with Application to Blasting Problems

A meshless local discontinuous Petrov-Galerkin (MLDPG) method based on the local symmetric weak form (LSWF) is presented with the application to blasting problems. The derivation is similar to that of mesh-based Runge-Kutta Discontinuous Galerkin (RKDG) method. The solutions are reproduced in a set of overlapped spherical sub-domains, and the test functions are employed from a partition of unity of the local basis functions. There is no need of any traditional non-overlapping mesh either for local approximation purpose or for Galerkin integration purpose in the presented method. The resulting MLDPG method is a meshless, stable, high-order accurate and highly parallelizable scheme which inherits both the advantages of RKDG and meshless method (MM), and it can handle the problems with extremely complicated physics and geometries easily. Three numerical examples of the one-dimensional Sod shock-tube problem, the blast-wave problem and the Woodward-Colella interacting shock wave problem are given. All the numerical results are in good agreement with the closed solutions. The higher-order MLDPG schemes can reproduce more accurate solution than the lower-order schemes.

LOCAL PETROV-GALERKIN METHOD FOR A THIN PLATE

The meshless local Petrov_Galerkin (MLPG) method for solving the bending problem of the thin plate were presented and discussed. The method used the moving least_squares approximation to interpolate the solution variables, and employed a local symmetric weak form. The present method was a truly meshless one as it did not need a finite element or boundary element mesh, either for purpose of interpolation of the solution, or for the integration of the energy. All integrals could be easily evaluated over regularly shaped domains (in general, spheres in three_dimensional problems) and their boundaries. The essential boundary conditions were enforced by the penalty method. Several numerical examples were presented to illustrate the implementation and performance of the present method. The numerical examples presented show that high accuracy can be achieved for arbitrary grid geometries for clamped and simply_supported edge conditions. No post processing procedure is required to computer the strain and stress, since the original solution from the present method, using the moving least squares approximation, is already smooth enough.

h-ADAPTIVE ANALYSIS BASED ON MESHLESS LOCAL PETROV-G ALERKIN METHOD WITH B SPLINE WAVELET FOR PLATES AND SHELLS

Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as the basis of the moving least square method to construct the meshless interpolation function. Multi-resolution analysis is used to decompose the field variables into high and low scales and the high scale component can commonly represent the gradient of the solution according to inherent characteristics of wavelets. The high scale component in the present method can directly detect high gradient regions of the field variables. The developed adaptive refinement scheme has been applied to simulate actual examples, and the effectiveness of the present adaptive refinement scheme has been verified.

ADAPTIVE MESHLESS METHOD BASED ON LOCAL FIT TECHNOLOGY

An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement. New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.

RESEARCH ON THE COMPANION SOLUTION FOR A THIN PLATE IN THE MESHLESS LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method(GFEM), boundary element method(BEM) and element free Galerkin method(EFGM), and is a truly meshless method possessing wide prospects in engineering applications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.

A meshless local Petrov–Galerkin method for solving the neutron diffusion equation

The goal of this study is to solve the neutron diffusion equation by using a meshless method and evaluate its performance compared to traditional methods. This paper proposes a novel method based on coupling the meshless local Petrov–Galerkin approach and the moving least squares approximation. This computational procedure consists of two main steps. The first involved applying the moving least squares approximation to construct the shape function based on the problem domain. Then, the obtained shape function was used in the meshless local Petrov–Galerkin method to solve the neutron diffusion equation.Because the meshless method is based on eliminating the mesh-based topologies, the problem domain was represented by a set of arbitrarily distributed nodes. There is no need to use meshes or elements for field variable interpolation. The process of node generation is simply and fully automated, which can save time. As this method is a local weak form, it does not require any background integration cells and all integrations are performed locally over small quadrature domains. To evaluate the proposed method,several problems were considered. The results were compared with those obtained from the analytical solution and a Galerkin finite element method. In addition, the proposed method was used to solve neutronic calculations in thesmall modular reactor. The results were compared with those of the citation code and reference values. The accuracy and precision of the proposed method were acceptable. Additionally, adding the number of nodes and selecting an appropriate weight function improved the performance of the meshless local Petrov–Galerkin method. Therefore, the proposed method represents an accurate and alternative method for calculating core neutronic parameters.

A Local Meshless Method for Two Classes of Parabolic Inverse Problems

A local meshless method is applied to find the numerical solutions of two classes of inverse problems in parabolic equations. The problem is reconstructing the source term using a solution specified at some internal points;one class is that the source term is time dependent, and the other class is that the source term is time and space dependent. Some numerical experiments are presented and discussed.

基于Voronoi结构的无网格局部Petrov-Galerkin方法

基于自然邻结点近似位移函数提出了一种用于求解弹性力学平面问题的无网格局部Petrov-Galerkin方法.这种方法在结构的求解域Ω内任意布置离散的结点,并且利用需求结点的自然邻结点和Voronoi结构来构造整体求解的近似位移函数.对于构造好的近似位移函数,在局部的Delaunay三角形子域上采用局部Petrov-Galerkin方法建立整体求解的平衡控制方程,这样平衡方程的积分可在背景三角形积分网格的形心上解析计算得到,而采用标准Galerkin方法的自然单元法需要三个数值积分点.该方法能够准确地施加边界条件,得到的系统矩阵是带状稀疏矩阵,对软件用户来说,它还是一种完全的、真正的无网格方法.所得计算结果表明,该方法的计算精度与有限元法四边形单元相当,但计算和形成系统平衡方程的时间比有限元法四边形单元提高了将近一倍,是一种理想的数值求解方法.

薄板的局部Petrov-Galerkin方法

利用薄板控制微分方程的等效积分对称弱形式和对变量(挠度)采用移动最小二乘近似函数进行插值,研究了薄板弯曲问题的无网格局部Petrov＿Galerkin方法· 这是一种真正的无网格方法,它不需要任何有限元或边界元网格,不管这种网格是用于能量积分还是进行插值的目的· 所有的积分都在规则形状的子域及其边界上进行,并用罚因子法施加本质边界条件· 数值例子表明,无网格局部Petrov＿Galerkin法不但能够求解二阶微分方程的边值问题,而且求解四阶微分方程的边值问题也很有效,也具有收敛快、稳定性好。

用无网格局部Petrov-Galerkin方法分析Winkler弹性地基板

利用Winkler弹性地基板控制微分方程的等效积分对称弱形式,同时对变量(挠度)采用移动最小二乘近似函数进行插值,研究了无网格局部Petrov Galerkin方法在弹性地基板弯曲问题中的应用.它不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,并用罚因子法施加本质边界条件.数值算例说明,无网格局部Petrov Galerkin法不但能够求解弹性静力学问题,而且在求解弹性地基板问题时仍具有收敛快、稳定性好和精度高的特点.

大变形问题分析的局部Petrov-Galerkin法

在微机电系统(MEMS)的建模和模拟研究中,大变形或大移动要充分予以考虑。用有限元法分析这类问题,由于难以避免的网格畸变,使模拟效率精度降低甚至失效,无网格方法(Meshless Method)则能在分析这类问题时显示出明显的优势,无网格局部Petrov-Galerkin(MLPG)法被誉为是一种有发展前景的真正无网格法。本文进一步发展了MLPG法,通过对任意的离散分布节点采用局部径向基函数构造插值形函数和Heaviside权函数,分析方程采用局部加权弱形式离散,建立了变量仅依赖于初始构型的完全Lagrange分析格式,最后用Newton-Raphson法迭代求解。文中分析了悬臂梁典型算例和微机电开关非线性大变形问题,通过与有限元结果的比较,表明本文提出的大变形问题无网格局部Petrov-Galerkin法具有稳定性好及收敛性快等优点。

无网格局部Petrov-Galerkin方法在弹塑性断裂力学问题中的应用

采用无网格局部Petrov-Galerkin方法来分析弹塑性断裂力学问题.这种无网格方法采用移动最小二乘法(MLS)来构造近似试函数和采用Heaviside函数作为加权残值法中的权函数,由于近似函数不满足Kronecker Delta条件,因此采用直接插值法来施加本质边界条件.如果不考虑体力,所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分.采用增量Newton-Raphson迭代法来求解弹塑性增量形式的局部Petrov-Galerkin方程.数值算例结果表明,该文方法对于弹塑性断裂力学问题的求解是可行的和有效的,并且所得到的结果具有较好的精度.

层合板分析的无网格局部Petrov-Galerkin方法

基于Kirchhoff均匀各向异性板控制方程的等效积分弱形式和对挠度函数采用移动最小二乘近似函数进行插值,进一步研究无网格局部Petrov-Galerkin方法在纤维增强对称层合板弯曲问题中的应用。该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,其问题的本质边界条件采用罚因子法来施加。通过数值算例和与其他方法的结果比较,表明无网格局部Petrov-Galerkin法求解层合薄板弯曲问题具有解的精度高、收敛性好等一系列优点。