• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.501-507
• /
• 2007

In the first part of this study, the moving least squares finite difference method for solving solid mechanics problems was formulated. This second part verified the accuracy, robustness and effectiveness of the developed method through several numerical examples. It was shown that the method gives excellent convergence rate for elasticity problem. The solution process of elastic crack problems showed the easiness in discontinuity modeling and demonstrates the accuracy and efficiency in finding singular stress solution based on adaptive node distribution. The applicability to the engineering problem with abrupt change in displacement and stresses gradient fields is verified through a localization band problem. The developed method is expected to be extended to the various special engineering problems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.411-420
• /
• 2009

This study presents an extended finite difference method based on moving least squares(MLS) method for solving potential problems with interfacial boundary. The approximation constructed from the MLS Taylor polynomial is modified by inserting of wedge functions for the interface modeling. Governing equations are node-wisely discretized without involving element or grid; immersion of interfacial condition into the approximation circumvents numerical difficulties owing to geometrical modeling of interface. Interface modeling introduces no additional unknowns in the system of equations but makes the system overdetermined. So, the numbers of unknowns and equations are equalized by the symmetrization of the stiffness matrix. Increase in computational effort is the trade-off for ease of interface modeling. Numerical results clearly show that the developed numerical scheme sharply describes the wedge behavior as well as jumps and efficiently and accurately solves potential problems with interface.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.321-327
• /
• 2007

This study presents a new gridless finite difference method for solving elastic crack problems. The method constructs the Taylor expansion based on the MLS(Moving Least Squares) method and effectively calculates the approximation and its derivatives without differentiation process. Since no connectivity between nodes is required, the modeling of discontinuity embedded in the domain is very convenient and discontinuity effect due to crack is naturally implemented in the construction of difference equations. Direct discretization of the governing partial differential equations makes solution process faster than other numerical schemes using numerical integration. Numerical results for mode I and II crack problems demonstrates that the proposed method accurately and efficiently evaluates the stress intensity factors.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.493-499
• /
• 2007

The Taylor expansion expresses a differentiable function and its coefficients provide good approximations for the given function and its derivatives. In this study, m-th order Taylor Polynomial is constructed and the coefficients are computed by the Moving Least Squares method. The coefficients are applied to the governing partial differential equation for solid problems including crack problems. The discrete system of difference equations are set up based on the concept of point collocation. The developed method effectively overcomes the shortcomings of the finite difference method which is dependent of the grid structure and has no approximation function, and the Galerkin-based meshfree method which involves time-consuming integration of weak form and differentiation of the shape function and cumbersome treatment of essential boundary.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.2
• /
• pp.139-148
• /
• 2012

The MLS(Moving Least Squares) Difference Method is a numerical scheme that combines the MLS method of Meshfree method and Taylor expansion involving not numerical quadrature or mesh structure but only nodes. This paper presents an dynamic algorithm of MLS difference method for solving transient solid mechanics problems. The developed algorithm performs time integration by using Newmark method and directly discretizes strong forms. It is very convenient to increase the order of Taylor polynomial because derivative approximations are obtained by the Taylor series expanded by MLS method without real differentiation. The accuracy and efficiency of the dynamic algorithm are verified through numerical experiments. Numerical results converge very well to the closed-form solutions and show less oscillation and periodic error than FEM(Finite Element Method).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.215-222
• /
• 2019

The heat conduction problem with discontinuous material coefficients generally consists of the conservative equation, boundary condition, and interface condition, which should be additionally satisfied in the solution procedure. This feature often makes the development of new numerical schemes difficult as it induces a layered singularity in the solution fields; thus, a special approximation is required to capture the singular behavior. In addition to the approximation, the construction of a total system of equations is challenging. In this study, a wedge function is devised for enriching the approximation, and the interface condition itself is embedded in the moving least squares(MLS) derivative approximation to consistently satisfy the interface condition. The heat conduction problem is then discretized in a strong form using the developed derivative approximation, which is named as the interface immersed MLS difference method. This method is able to efficiently provide a numerical solution for such interface problems avoiding both numerical quadrature as well as extra difference equations related to the interface condition enforcement. Numerical experiments proved that the developed numerical method was highly accurate and computationally efficient at solving the heat conduction problem with interfacial jump as well as the problem with a geometrically induced interfacial singularity.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.17-26
• /
• 2014

This paper presents a dynamic crack propagation algorithm based on the Moving Least Squares(MLS) difference method. The derivative approximation for the MLS difference method is derived by Taylor expansion and moving least squares procedure. The method can analyze dynamic crack problems using only node model, which is completely free from the constraint of grid or mesh structure. The dynamic equilibrium equation is integrated by the Newmark method. When a crack propagates, the MLS difference method does not need the reconstruction of mode model at every time step, instead, partial revision of nodal arrangement near the new crack tip is carried out. A crack is modeled by the visibility criterion and dynamic energy release rate is evaluated to decide the onset of crack growth together with the corresponding growth angle. Mode I and mixed mode crack propagation problems are numerically simulated and the accuracy and stability of the proposed algorithm are successfully verified through the comparison with the analytical solutions and the Element-Free Galerkin method results.