• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.3
• /
• pp.237-244
• /
• 2016

This paper presents a nonlinear Moving Least Squares(MLS) difference method for material nonlinearity problem. The MLS difference method, which employs strong formulation involving the fast derivative approximation, discretizes governing partial differential equation based on a node model. However, the conventional MLS difference method cannot explicitly handle constitutive equation since it solves solid mechanics problems by using the Navier's equation that unifies unknowns into one variable, displacement. In this study, a double derivative approximation is devised to treat the constitutive equation of inelastic material in the framework of strong formulation; in fact, it manipulates the first order derivative approximation two times. The equilibrium equation described by the divergence of stress tensor is directly discretized and is linearized by the Newton method; as a result, an iterative procedure is developed to find convergent solution. Stresses and internal variables are calculated and updated by the return mapping algorithm. Effectiveness and stability of the iterative procedure is improved by using algorithmic tangent modulus. The consistency of the double derivative approximation was shown by the reproducing property test. Also, accuracy and stability of the procedure were verified by analyzing inelastic beam under incremental tensile loading.

• 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.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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.5
• /
• pp.577-588
• /
• 2011

This paper provides a novel extended Moving Least Squares(MLS) difference method for the potential problem with weak and strong discontinuities. The conventional MLS difference method is enhanced with jump functions such as step function, wedge function and scissors function to model discontinuities in the solution and the derivative fields. When discretizing the governing equations, additional unknowns are not yielded because the jump functions are decided from the known interface condition. The Poisson type PDE's are discretized by the difference equations constructed on nodes. The system of equations built up by assembling the difference equations are directly solved, which is very efficient. Numerical examples show the excellence of the proposed numerical method. The method is expected to be applied to various discontinuity related problems such as crack problem, moving boundary problem and interaction problems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.5
• /
• pp.279-286
• /
• 2020

In this study, we investigate the accuracy of higher order derivatives in the moving least square (MLS) difference method. An interpolation function is constructed by employing a Taylor series expansion via MLS approximation. The function is then applied to the mixed variational theorem in which the displacement and stress resultants are treated as independent variables. The higher order derivatives are evaluated by solving simply supported beams and cantilevers. The results are compared with the analytical solutions in terms of the order of polynomials, support size of the weighting function, and number of nodes. The accuracy of the higher order derivatives improves with the employment of the mean value theorem, especially for very high-order derivatives (e.g., above fourth-order derivatives), which are important in a classical asymptotic analysis.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.179-182
• /
• 2010

본 연구에서는 MLS 차분법을 이용하여 동역학 문제를 해석하기 위한 explicit 동적해석 알고리즘을 제시한다. 격자망이 없는 장점을 부각시키기 위해 이동최소제곱법에 근거한 Taylor 전개로부터 미분근사를 얻고 차분식을 구성했다. 지배 미분방정식의 시간항을 CDM(Central difference Method) 차분하여 빠른 속도로 동적해석을 수행하였다. 수치결과를 통해 본 연구에서 제시한 알고리즘의 정확성과 안정성을 확인할 수 있었다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.583-590
• /
• 2016

This paper presents a dynamic crack propagation algorithm with Rayleigh damping effect based on the MLS(Moving Least Squares) Difference Method. Dynamic equilibrium equation and constitutive equation are derived by considering Rayliegh damping and governing equations are discretized by the MLS derivative approximation; the proportional damping, which has not been properly treated in the conventional strong formulations, was implemented in both the equilibrium equation and constitutive equation. Dynamic equilibrium equation including time relevant terms is integrated by the Central Difference Method and the discrete equations are simplified by lagging the velocity one step behind. A geometrical feature of crack is modeled by imposing the traction-free condition onto the nodes placed at crack surfaces and the effect of movement and addition of the nodes at every time step due to crack growth is appropriately reflected on the construction of total system. The robustness of the proposed numerical algorithm was proved by simulating single and multiple crack growth problems and the effect of proportional damping on the dynamic crack propagation analysis was effectively demonstrated.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.2-7
• /
• 2008

A highly efficient moving least squares finite difference method (MLS FDM) for heat transfer analysis of composite material with interface. In the MLS FDM, governing differential equations are directly discretized at each node. No grid structure is required in the solution procedure. The discretization of governing equations are done by Taylor expansion based on moving least squares method. A wedge function is designed for the modeling of the derivative jump across the interface. Numerical examples showed that the numerical scheme shows very good computational efficiency together with high aocuracy so that the scheme for heat transfer problem with different heat conductivities was successfully verified.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.12
• /
• pp.1260-1268
• /
• 2009

Three methods(the stepped sine method, the statistical method(random excitation method) and the maximum-length sequence(MLS) method) for head-related transfer functions(HRTFs) are experimentally compared in view point of accuracy and efficiency. First, the stepped sine method has high signal-to-noise ratio, but low efficiency. Second, the statistical method is fast measurement speed, but weak to noise than the other methods. Finally, the MLS method shows both good efficiency and high signal-to-noise ratio, but it needs additional software or equipment such as MLS signal generator. For comparison of measurement accuracy, HRTFs of KEMAR dummy are measured for various azimuths and elevations. Error norms for magnitude and phase of HRTFs are defined and calculated for the measured HRTFs. The calculated error norms show that the methods give similar results in magnitude and phase except a little phase difference in the MLS method.