• Interaction and multiscale mechanics
• /
• v.1 no.2
• /
• pp.251-278
• /
• 2008

In this study, we investigate the discontinuous-derivative treatment at the gas-liquid interface in underwater explosion (UNDEX) problems by using the Moving Least Squares-Smoothed Particle Hydrodynamics (MLS-SPH) method, which is known as one of the particle methods suitable for problems where large deformation and inhomogeneity occur in the whole domain. Because the numerical oscillation of pressure arises from derivative discontinuity in the UNDEX analysis using the standard SPH method, the MLS shape function with Discontinuous-derivative Basis Function (DBF) that is able to represent the derivative discontinuity of field function is utilized in the MLS-SPH formulation in order to suppress the nonphysical pressure oscillation. The effectiveness of the MLS-SPH with DBF is demonstrated in comparison with the standard SPH and conventional MLS-SPH though a shock tube problem and benchmark standard problems of UNDEX of a trinitrotoluene (TNT) charge.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1089-1093
• /
• 2003

The moving least square (MLS) basis is implemented for the real-space band-structure calculation of 2D photonic crystals. The value-periodic MLS shape function is thus used in order to represent the periodicity of crystal lattice. Any periodic function can properly be reproduced using this shape function. Matrix eigenequations, derived from the macroscopic Maxwell equations, are then solved to obtain photonic band structures. Through numerical examples of several lattice structures, the MLS-based method is proved to be a promising scheme for predicting band gaps of photonic crystals.

• Structural Engineering and Mechanics
• /
• v.25 no.1
• /
• pp.91-106
• /
• 2007

A new class of finite elements is described for dealing with mesh gradation. The approach employs the moving least square (MLS) scheme to devise a class of elements with an arbitrary number of nodal points on the parental domain. This approach generally leads to elements with rational shape functions, which significantly extends the function space of the conventional finite element method. With a special choice of the nodal points and the base functions, the method results in useful elements with polynomial shape functions for which the $C^1$ continuity breaks down across the boundaries between the subdomains comprising one element. Among those, (4 + n)-noded MLS based finite elements possess the generality to be connected with an arbitrary number of linear elements at a side of a given element. It enables us to connect one finite element with a few finite elements without complex remeshing. The effectiveness of the new elements is demonstrated via appropriate numerical examples.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.335-341
• /
• 2009

In this paper, review of developed MLS-based finite elements and a proposal for their applications are described. The shape functions and their derivatives of MLS-based finite elements are constructed using Moving-Least Square approximation. In MLS-based finite element, using the adequate influence domain of weight function used in MLS approximation, kronecker delta condition could be satisfied at the element boundary. Moreover, because of the characteristics of MLS approximation, we could easily add extra nodes at an arbitrary position in MLS-based finite elements. For these reasons, until now, several variable-node elements(2D variable element for linear case and quadratic case and 3D variable-node elements) and finite crack elements are developed using MLS-based finite elements concept. MLS-based finite elements could be extended to 2D variable-node triangle element, 2D finite crack triangle element, variable-node shell element, finite crack shell element, and 3D polyhedron element. In this paper, we showed the feasibility of 3D polyhedron element at the case of femur meshing.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.12 s.255
• /
• pp.1626-1633
• /
• 2006

In the present paper, in the context of the meshless interpolation of a moving least squares (MLS) type, a novel method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The support domains for the shape functions in the MLS approximation are defined from the primary nodes, and the secondary nodes use the same support domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and moved without an additional mesh. Several numerical examples are presented to illustrate the effectiveness of the present method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.6
• /
• pp.331-337
• /
• 2018

This paper presents dynamic algorithm of the MLS(moving least squares) difference method using first order differential Approximation. The governing equations are only discretized by the first order MLS derivative approximation. The system equation consists of an assembly of the approximate function, so the shape of system equation is similar to FEM(finite element method). The CDM(central difference method) is used for time integration of dynamic equilibrium equation. The natural frequency analyses of the MLS difference method and FEM are performed, and two analysis results are compared. Also, the accuracy of the proposed numerical method is verified by displaying the dynamic analysis results together with the results by the existing second order differential approximation. In the process of assembling the first order MLS derivative approximation, the oscillation error was suppressed and the stress distribution was interpreted as relatively uniform.

• Steel and Composite Structures
• /
• v.12 no.1
• /
• pp.1-11
• /
• 2012

In this paper static analysis of FGM cylinders subjected to internal and external pressure was carried out by a mesh-free method. In this analysis MLS shape functions are used for approximation of displacement field in the weak form of equilibrium equation and essential boundary conditions are imposed by transformation method. Mechanical properties of cylinders were assumed to be variable in the radial direction. Two types of cylinders were analyzed in this work. At first cylinders with infinite length were considered and results obtained for these cylinders were compared with analytical solutions and a very good agreement was seen between them. Then the proposed mesh-free method was used for analysis of cylinders with finite length and two different types of boundary conditions. Results obtained from these analyses were compared with results of finite element analyses and a very good agreement was seen between them.

• Structural Engineering and Mechanics
• /
• v.11 no.2
• /
• pp.123-132
• /
• 2001

An improved version of the Element-free Galerkin method (EFGM) is presented here for addressing the problem of transverse shear locking in shear-deformable beams with a high length over thickness ratio. Based upon Timoshenko's theory of thick beams, it has been recognized that shear locking will be completely eliminated if the rotation field is constructed to match the field of slope, given by the first derivative of displacement. This criterion is applied directly to the most commonly implemented version of EFGM. However in the numerical process to integrate strain energy, the second derivative of the standard Moving Least Square (MLS) shape functions must be evaluated, thus requiring at least a $C^1$ continuity of MLS shape functions instead of $C^0$ continuity in the conventional EFGM. Yet this hindrance is overcome effortlessly by only using at least a $C^1$ weight function. One-dimensional quartic spline weight function with $C^2$ continuity is therefore adopted for this purpose. Various numerical results in this work indicate that the modified version of the EFGM does not exhibit transverse shear locking, reduces stress oscillations, produces fast convergence, and provides a surprisingly high degree of accuracy even with coarse domain discretizations.

• Wind and Structures
• /
• v.27 no.4
• /
• pp.235-245
• /
• 2018

In this study the stress analysis of orthotropic thin plate with arbitrary shapes for different boundary conditionsis investigated. Meshfreemethod is applied to static analysis of thin plates with various geometries based on the Kirchhoff classical plate theory. According to the meshfree method the domain of the plates are expressed through a set of nodes without using mesh. In this method, a set of nodes are defined in a standard rectangular domain, then via a third order map, these nodes are transferred to the main domain of the original geometry; therefore the analysis of the plates can be done. Herein, Meshless local Petrov-Galerkin (MLPG) as a meshfree numerical method is utilized. The MLS function in MLPG does not satisfy essential boundary conditions using Delta Kronecker. In the MLPG method, direct interpolation of the boundary conditions can be applied due to constructing node by node of the system equations. The detailed parametric study is conducted, focusing on the arbitrary geometries of the thin plates. Results show that the meshfree method provides better accuracy rather than finite element method. Also, it is found that trend of the figures have good agreement with relevant published papers.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.15 no.2
• /
• pp.101-108
• /
• 2007

This paper presents an efficient meshless method for analyzing cracked piezoelectric structures subjected to mechanical and electrical loading. The method employs an element free Galerkin (EFG) formulation and an enriched basic function as well as special shape functions that contain discontinuous derivatives. Based on the moving least squares (MLS) interpolation approach, The EFG method is one of the promising methods for dealing with problems involving progressive crack growth. Since the method is meshless and no element connectivity data are needed, the burdensome remeshing procedure required in the conventional finite element method (FEM) is avoided. The numerical results show that the proposed method yields an accurate near-tip stress field in an infinite piezoelectric plate containing an interior hole. Another example is to study a ceramic multilayer actuator. The proposed model was found to be accurate in the simulation of stress and electric field concentrations due to the abrupt end of an internal electrode.