• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.94-109
• /
• 2006

An improved meshfree method called the Petrov-Galerkin natural element (PG-NE) method is introduced in order to secure the numerical integration accuracy. As in the Bubnov-Galerkin natural element (BG-NE) method, we use Laplace interpolation function for the trial basis function and Delaunay triangles to define a regular integration background mesh. But, unlike the BG-NE method, the test basis function is differently chosen, based on the Petrov-Galerkin concept, such that its support coincides exactly with a regular integration region in background mesh. Illustrative numerical experiments verify that the present method successfully prevents the numerical accuracy deterioration stemming from the numerical integration error.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1274-1279
• /
• 2003

In this paper, a new meshfree technique which improves the numerical integration accuracy is introduced. This new method called the Petrov-Galerkin natural element(PG-NE) is based on the Voronoi diagram and the Delaunay triangulation which is based on the same concept used for conventional natural element method called the Bubnov-Galerkin natural element(BG-NE). But, unlike BG-NE method, the test shape function is differently chosen from the trial shape function. The proposed technique ensures that the numerical integration error is remarkably reduced.

• Structural Engineering and Mechanics
• /
• v.56 no.4
• /
• pp.589-603
• /
• 2015

The mixed-mode stress intensity factors of 2-D angled cracks are evaluated by Petrov-Galerkin natural element (PG-NE) method in which Voronoi polygon-based Laplace interpolation functions and CS-FE basis functions are used for the trial and test functions respectively. The interaction integral is implemented in a frame of PG-NE method in which the weighting function defined over a crack-tip integral domain is interpolated by Laplace interpolation functions. Two Cartesian coordinate systems are employed and the displacement, strains and stresses which are solved in the grid-oriented coordinate system are transformed to the other coordinate system aligned to the angled crack. The present method is validated through the numerical experiments with the angled edge and center cracks, and the numerical accuracy is examined with respect to the grid density, crack length and angle. Also, the stress intensity factors obtained by the present method are compared with other numerical methods and the exact solution. It is observed from the numerical results that the present method successfully and accurately evaluates the mixed-mode stress intensity factors of 2-D angled cracks for various crack lengths and crack angles.