共查询到20条相似文献，搜索用时 0 毫秒
1.
S. Suleau Ph. Bouillard 《International journal for numerical methods in engineering》2000,47(6):1169-1188
The standard finite element method (FEM) is unreliable to compute approximate solutions of the Helmholtz equation for high wave numbers due to the dispersion, unless highly refined meshes are used, leading to unacceptable resolution times. The paper presents an application of the element‐free Galerkin method (EFG) and focuses on the dispersion analysis in one dimension. It shows that, if the basis contains the solution of the homogenized Helmholtz equation, it is possible to eliminate the dispersion in a very natural way while it is not the case for the finite element methods. For the general case, it also shows that it is possible to choose the parameters of the method in order to minimize the dispersion. Finally, theoretical developments are validated by numerical experiments showing that, for the same distribution of nodes, the element‐free Galerkin method solution is much more accurate than the finite element one. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
2.
L. Gavete J. L. Cuesta A. Ruiz 《International journal for numerical methods in engineering》2002,53(3):677-690
In this paper, we present a procedure to estimate the error in elliptic equations using the element‐free Galerkin (EFG) method, whose evaluation is computationally simple and can be readily implemented in existing EFG codes. The estimation of the error works very well in all numerical examples for 2‐D potential problems that are presented here, for regular and irregular clouds of points. Moreover, it was demonstrated that this method is very simple in terms of economy and gives a good performance. The results show that the error in EFG approximation may be estimated via the error estimator described in this paper. The quality of the estimation of the error is demonstrated by numerical examples. The implemented procedure of error approximation allows the global energy norm error to be estimated and also gives a good evaluation of local errors. It can, thus, be combined with a full adaptive process of refinement or, more simply, provide guidance for redesign of cloud of points. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
3.
Gang Wang Guodong Zeng Xiangyang Cui Shizhe Feng 《International journal for numerical methods in engineering》2019,120(4):473-497
This paper reports a detailed analysis on the numerical dispersion error in solving one-, two-, and three-dimensional acoustic problems governed by the Helmholtz equation using the gradient weighted finite element method (GW-FEM) in comparison with the standard FEM and the modified methods presented in the literatures. The discretized system equations derived based on the gradient weighted operation corresponding to the considered method are first briefed. The discrete dispersion relationships relating the exact and numerical wave numbers defined in different dimensions are then formulated, which will be further used to investigate the dispersion effect mainly caused by the approximation of field variables. The influence of nondimensional wave number and wave propagation angle on the dispersion error is detailedly studied. Comparisons are made with the classical FEM and high-performance algorithms. Results of both theoretical and numerical experiments show that the present method can effectively reduce the pollution effect in computational acoustics owning to its crucial effectiveness in handing the dispersion error in the discrete numerical model. 相似文献
4.
L. Gavete M. L. Gavete B. Alonso A. J. Martín 《International journal for numerical methods in engineering》2003,58(15):2239-2263
Recently, considerable effort has been devoted to the development of the so‐called meshless methods. Meshless methods still require considerable improvement before they equal the prominence of finite elements in computer science and engineering. One of the paths in the evolution of meshless methods has been the development of the element free Galerkin (EFG) method. In the EFG method, it is obviously important that the ‘a posteriori error’ should be approximated. An ‘a posteriori error’ approximation based on the moving least‐squares method is proposed, using the solution, computed from the EFG method. The error approximation procedure proposed in this paper is simple to construct and requires, at most, nearest neighbour information from the EFG solution. The formulation is based on employing different moving least‐squares approximations. Different selection strategies of the moving least‐squares approximations have been used and compared, to obtain optimum values of the parameters involved in the approximation of the error. The performance of the developed approximation of the error is illustrated by analysing different examples for two‐dimensional (2D) potential and elasticity problems, using regular and irregular clouds of points. The implemented procedure of error approximation allows the global energy norm error to be estimated and also provides a good evaluation of local errors. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
5.
Xiaoying Zhuang Charles Augarde 《International journal for numerical methods in engineering》2010,81(3):366-380
The element‐free Galerkin (EFG) method is probably the most widely used meshless method at present. In the EFG method, shape functions are derived from a moving least‐squares approximation using a polynomial basis, a calculation involving the inversion of a small matrix. A new implementation of the EFG method was published soon after the original where an alternative approach using an orthogonal basis was proposed to avoid matrix inversion in the formulation of the shape functions. In this paper we revisit this topic and show that the difficulties associated with the use of a polynomial basis remain present in the orthogonal case. We also show that certain terms in the derivative expressions are omitted in the new implementation of the EFG, which can lead to errors. Finally, we propose a new approach that avoids inversion while maintaining accuracy. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
6.
Xiong Zhang Mingwan Lu J. L. Wegner 《International journal for numerical methods in engineering》2000,47(10):1649-1661
According to the characteristic structural features of jointed rock structures, a meshless model is proposed for the mechanics analysis of jointed rock structures based on the moving least‐squares interpolants. In this model, a jointed rock structure is regarded as a system of relatively intact rock blocks connected by joints or planes of discontinuity; these rock blocks are modelled by general shaped anisotropic blocks while these joints and planes of discontinuity are modelled by interfaces. The displacement field of each block is constructed by the moving least‐squares interpolants with an array of points distributed in the block. To deal with the discontinuities of rock structures, the displacement fields are constructed to be discontinuous between blocks. The displacement fields and their gradients are continuous in each block, hence no post processing is required for the output of strains and stresses. The finite element mesh is totally unnecessary, so the time‐consuming mesh generation is avoided. The rate of convergence can exceed that of finite elements significantly, and a high resolution of localized steep gradients can be achieved. Furthermore, the discontinuities of rock structures are also fully taken into consideration. The present method is developed for two‐dimensional linear elastic analysis of jointed rock structures, and can be extended to three‐dimensional and non‐linear analysis. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
7.
Arnaud Deraemaeker Ivo Babu&#x;ka Philippe Bouillard 《International journal for numerical methods in engineering》1999,46(4):471-499
For high wave numbers, the Helmholtz equation suffers the so‐called ‘pollution effect’. This effect is directly related to the dispersion. A method to measure the dispersion on any numerical method related to the classical Galerkin FEM is presented. This method does not require to compute the numerical solution of the problem and is extremely fast. Numerical results on the classical Galerkin FEM (p‐method) is compared to modified methods presented in the literature. A study of the influence of the topology triangles is also carried out. The efficiency of the different methods is compared. The numerical results in two of the mesh and for square elements show that the high order elements control the dispersion well. The most effective modified method is the QSFEM [1,2] but it is also very complicated in the general setting. The residual‐free bubble [3,4] is effective in one dimension but not in higher dimensions. The least‐square method [1,5] approach lowers the dispersion but relatively little. The results for triangular meshes show that the best topology is the ‘criss‐cross’ pattern. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
8.
Leopoldo P. Franca Antonini P. Macedo 《International journal for numerical methods in engineering》1998,43(1):23-32
A two-level finite element method is introduced and its application to the Helmholtz equation is considered. The method retains the desirable features of the Galerkin method enriched with residual-free bubbles, while it is not limited to discretizations using elements with simple geometry. The method can be applied to other equations and to irregular-shaped domains. © 1998 John Wiley & Sons, Ltd. 相似文献
9.
10.
Assad A. Oberai Peter M. Pinsky 《International journal for numerical methods in engineering》2000,49(3):399-419
A new residual‐based finite element method for the scalar Helmholtz equation is developed. This method is obtained from the Galerkin approximation by appending terms that are proportional to residuals on element interiors and inter‐element boundaries. The inclusion of residuals on inter‐element boundaries distinguishes this method from the well‐known Galerkin least‐squares method and is crucial to the resulting accuracy of this method. In two dimensions and for regular bilinear quadrilateral finite elements, it is shown via a dispersion analysis that this method has minimal phase error. Numerical experiments are conducted to verify this claim as well as test and compare the performance of this method on unstructured meshes with other methods. It is found that even for unstructured meshes this method retains a high level of accuracy. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
11.
Antonio Huerta Sonia Fernndez‐Mndez 《International journal for numerical methods in engineering》2001,51(11):1361-1383
Volumetric locking (locking in the incompressible limit) for linear elastic isotropic materials is studied in the context of the element‐free Galerkin method. The modal analysis developed here shows that the number of non‐physical locking modes is independent of the dilation parameter (support of the interpolation functions). Thus increasing the dilation parameter does not suppress locking. Nevertheless, an increase in the dilation parameter does reduce the energy associated with the non‐physical locking modes; thus, in part, it alleviates the locking phenomena. This is shown for linear and quadratic orders of consistency. Moreover, the biquadratic order of consistency, as in finite elements, improves the locking behaviour. Although more locking modes are present in the element‐free Galerkin method with quadratic consistency than with standard biquadratic finite elements. Finally, numerical examples are shown to validate the modal analysis. In particular, the conclusions of the modal analysis are also confirmed in an elastoplastic example. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
12.
Frank Ihlenburg Ivo Babu&#x;ka 《International journal for numerical methods in engineering》1995,38(22):3745-3774
When applying numerical methods for the computation of stationary waves from the Helmholtz equation, one obtains ‘numerical waves’ that are dispersive also in non-dispersive media. The numerical wave displays a phase velocity that depends on the parameter k of the Helmholtz equation. In dispersion analysis, the phase difference between the exact and the numerical solutions is investigated. In this paper, the authors' recent result on the phase difference for one-dimensional problems is numerically evaluated and discussed in the context of other work directed to this topic. It is then shown that previous error estimates in H1-norm are of nondispersive character but hold for medium or high wavenumber on extremely refined mesh only. On the other hand, recently proven error estimates for constant resolution contain a pollution term. With certain assumptions on the exact solution, this term is of the order of the phase difference. Thus a link is established between the results of dispersion analysis and the results of numerical analysis. Throughout the paper, the presentation and discussion of theoretical results is accompanied by numerical evaluation of several model problems. Special attention is given to the performance of the Galerkin method with a higher order of polynomial approximation p(h-p-version). 相似文献
13.
Gye‐Hee Lee Heung‐Jin Chung Chang‐Koon Choi 《International journal for numerical methods in engineering》2003,56(3):331-350
In this paper, an adaptive analysis of crack propagation based on the error estimation by the element‐free Galerkin (EFG) method is presented. The adaptivity analysis in quasi‐static crack propagation is achieved by adding and/or removing the nodes along the background integration cells, those are refined or recovered according to the estimated errors. These errors are obtained basically by calculating the difference between the values of the projected stresses and original EFG stresses. To evaluate the performance of the proposed adaptive procedure, the crack propagation behaviour is investigated for several examples. The results of these examples show the efficiency and accuracy of the proposed scheme in crack propagation analysis. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
14.
W. Kanok‐Nukulchai W. Barry K. Saran‐Yasoontorn P. H. Bouillard 《International journal for numerical methods in engineering》2001,52(7):705-725
In this study, a method for completely eliminating the presence of transverse shear locking in the application of the element‐free Galerkin method (EFGM) to shear‐deformable beams and plates is presented. The matching approximation fields concept of Donning and Liu has shown that shear locking effects may be prevented if the approximate rotation fields are constructed with the innate ability to match the approximate slope (first derivative of displacement) fields and is adopted. Implementation of the matching fields concept requires the computation of the second derivative of the shape functions. Thus, the shape functions for displacement fields, and therefore the moving least‐squares (MLS) weight function, must be at least C1 continuous. Additionally, the MLS weight functions must be chosen such that successive derivatives of the MLS shape function have the ability to exactly reproduce the functions from which they were derived. To satisfy these requirements, the quartic spline weight function possessing C2 continuity is used in this study. To our knowledge, this work is the first attempt to address the root cause of shear locking phenomenon within the framework of the element‐free Galerkin method. Several numerical examples confirm that bending analyses of thick and thin beams and plates, based on the matching approximation fields concept, do not exhibit shear locking and provide a high degree of accuracy for both displacement and stress fields. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
15.
Christina Wenterodt Otto von Estorff 《International journal for numerical methods in engineering》2009,77(12):1670-1689
When numerical methods such as the finite element method (FEM) are used to solve the Helmholtz equation, the solutions suffer from the so‐called pollution effect which leads to inaccurate results, especially for high wave numbers. The main reason for this is that the wave number of the numerical solution disagrees with the wave number of the exact solution, which is known as dispersion. In order to obtain admissible results a very high element resolution is necessary and increased computational time and memory capacity are the consequences. In this paper a meshfree method, namely the radial point interpolation method (RPIM), is investigated with respect to the pollution effect in the 2D‐case. It is shown that this methodology is able to reduce the dispersion significantly. Two modifications of the RPIM, namely one with polynomial reproduction and another one with a problem‐dependent sine/cosine basis, are also described and tested. Numerical experiments are carried out to demonstrate the advantages of the method compared with the FEM. For identical discretizations, the RPIM yields considerably better results than the FEM. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
16.
17.
Tinh Quoc Bui Tan Nhat Nguyen Hung Nguyen‐Dang 《International journal for numerical methods in engineering》2009,77(10):1371-1395
This paper mainly proposes an alternative way for numerical implementation of thin plates bending based on a new improvement of meshless method, which is combined between the standard element‐free Galerkin method and one different shape functions building technique. The moving Kriging (MK) interpolation is applied instead of the traditional moving least‐square approximation in order to overcome Kronecker's delta property where the standard method does not satisfy. Obviously, the deflection of the thin plates is approximated via the MK interpolation. To illustrate this approach, numerical analysis is examined in both regular and irregular systems. Three examples with different geometric shapes of thin plates undergoing a simply supported boundary are performed. In addition, two important parameters of the present method are also analyzed. A good agreement can be found among the proposed, analytical and finite element methods. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
18.
Eran Grosu Isaac Harari 《International journal for numerical methods in engineering》2009,78(11):1261-1291
19.
Weilong Ai Charles E. Augarde 《International journal for numerical methods in engineering》2016,108(13):1626-1648
Computational modelling of fracture has been attempted in the past with a range of numerical approaches including finite element, extended finite element and meshless methods. The cracking particle method (CPM) of Rabczuk is a pragmatic alternative to explicit modelling of crack surfaces in which a crack is represented by a set of cracking particles that can be easily updated when the crack propagates. The change of cracking angle is recorded in discrete segments of broken lines, which makes this methodology suitable to model discontinuous cracks. In this paper, a new CPM is presented that improves on two counts: firstly, crack path curvature modelling is improved by the use of bilinear segments centred at each particle and secondly, efficiency for larger problems is improved via an adaptive process of both refinement and recovery. The system stiffness is calculated and stored in local matrices, so only a small influenced domain should be recalculated for each step while the remainder can be read directly from storage, which greatly reduces the computational expense. The methodology is applied to several 2D crack problems, and good agreement to analytical solutions and previous work is obtained. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
20.
Harm Askes Jerzy Pamin Ren de Borst 《International journal for numerical methods in engineering》2000,49(6):811-832
Gradient‐dependent damage formulations incorporate higher‐order derivatives of state variables in the constitutive equations. Different formulations have been derived for this gradient enhancement, comparison of which is difficult in a finite element context due to higher‐order continuity requirements for certain formulations. On the other hand, the higher‐order continuity requirements are met naturally by element‐free Galerkin (EFG) shape functions. Thus, the EFG method provides a suitable tool for the assessment of gradient enhanced continuum models. Dispersion analyses have been carried out to compare different gradient enhanced models with the non‐local damage model. The formulation of the additional boundary conditions is addressed. Numerical examples show the objectivity with respect to the discretization and the differences between various gradient formulations with second‐ and fourth‐order derivatives. It is shown that with the same underlying internal length scale, very different results can be obtained. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献