Journal of the Computational Structural Engineering Institute of Korea (한국전산구조공학회논문집)

Computational Structural Engineering Institute of Korea (한국전산구조공학회)
권호 표지

Intrinsically Extended Moving Least Squares Finite Difference Method for Potential Problems with Interfacial Boundary

계면경계를 갖는 포텐셜 문제 해석을 위한 내적확장된 이동최소제곱 유한차분법

  • 윤영철 (명지전문대학 토목과) ;
  • 이상호 (연세대학교 사회환경시스템공학부)
  • Received : 2009.05.07
  • Accepted : 2009.08.18
  • Published : 2009.10.30
Download PDF
⟨ Previous Next ⟩

Abstract

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.

본 연구는 계면경계를 갖는 포텐셜 문제의 해석를 위한 이동최소제곱 기반의 확장된 유한차분법을 제시한다. 이동최소제곱법을 이용한 Taylor 전개로부터 얻어진 근사함수에 쐐기함수를 도입하여 계면경계의 특이성을 모사한다. 지배방정식은 요소나 그리드없이 절점만을 이용해 이산화한다. 계면경계의 특이성은 절점에서 구성되는 근사식에 매입되기 때문에 계면경계의 기하학적 모델링으로 발생하는 수치적인 어려움을 피할 수 있다. 계면경계 조건으로 인해 전체 계방정식에 추가되는 미지수는 없지만, 계방정식을 과결정 시스템으로 만드므로 강성도 행렬을 대칭화하여 미지수와 방정식의 개수를 일치시켰다. 이로 인한 계산량 증가는 계면경계 모델링의 간소화로 인한 수치적인 이득과 맞바꿀 수 있다. 다양한 수치적 검증을 통해 개발된 해석기법이 쐐기거동과 점프를 성공적으로 묘사할 뿐만 아니라 계면경계를 갖는 포텐셜 문제 효율적이고 정확하게 해석할 수 있음을 보였다.

Keywords

References

  1. 윤영철, 김도완 (2007) 이동최소제곱 유한차분법을 이용한 계면경계를 갖는 이종재료의 열전달문제 해석, 한국전산구조공학회 논문집, 20(6), pp. 779-787
  2. 윤영철, 김동조, 이상호 (2007a) 탄성균열해석을 위한 그리드 없는 유한차분법, 한국전산구조공학회 논문집, 20(3), pp.321-327
  3. 윤영철, 김효진, 김동조, 윙 캠 리우, 테드 벨리체코, 이상호 (2007b) 이동최소제곱 유한차분법을 이용한 응력집중문제해석(I) : 고체문제의 정식화, 한국전산구조공학회 논문집, 20(4), pp.493-499
  4. 윤영철, 김효진, 김동조, 윙 캠 리우, 테드 벨리체코, 이상호 (2007c) 이동최소제곱 유한차분법을 이용한 응력집중문제해석(II) : 균열과 국소화 밴드 문제로의 적용, 한국전산구조공학회 논문집, 20(4), pp.501-507
  5. Aluru, N.R. (2000) A point collocation method based on reproducing kernel approximations, International Journal for Numerical Methods in Engineering, 47, pp.1083-1121 https://doi.org/10.1002/(SICI)1097-0207(20000228)47:6<1083::AID-NME816>3.0.CO;2-N
  6. Belytschko, T., Fleming, M. (1999) Smoothing, enrichment and contact in the Element-Free Galerkin method, Computer and Structures, 71, pp.173-195 https://doi.org/10.1016/S0045-7949(98)00205-3
  7. Belytschko, T., Lu, Y.Y., Gu, L. (1994) Elementfree Galerkin methods. International Journal for Numerical Methods in Engineering, 37, pp.229-256 https://doi.org/10.1002/nme.1620370205
  8. Benzley, S.E. (1974) Representation of singularities with isoparametric finite elements, International Journal for Numerical Methods in Engineering, 8, pp.537-545 https://doi.org/10.1002/nme.1620080310
  9. Fleming, M., Chu, Y.A., Moran, B., Belytschko, T. (1997) Enrichment Element-Free Galerkin methods for crack tip fields, International Journal for Numerical Methods in Engineering, 40, pp.1483-1504 https://doi.org/10.1002/(SICI)1097-0207(19970430)40:8<1483::AID-NME123>3.0.CO;2-6
  10. Fries, T.P., Belytschko, T. (2006) The intrinsic XFEM：a method for arbitrary discontinuities without additional unknowns, International Journal for Numerical Methods in Engineering, 68, pp.1358 -1385 https://doi.org/10.1002/nme.1761
  11. Gifford, Jr., L.N., Hilton, P.D. (1978) Stress intensity factors by enriched finite elements, Engineering Fracture Mechanics, 10, pp.485-496 https://doi.org/10.1016/0013-7944(78)90059-0
  12. Kim, D.W., Yoon, Y.C., Liu, W.K., Belytschko, T. (2007a) Extrinsic meshfree approximation using asymptotic expansion for interfacial discontinuity of derivative, Journal of Computational Physics, 221, pp.370-394 https://doi.org/10.1016/j.jcp.2006.06.023
  13. Kim, D.W., Yoon, Y.C., Liu, W.K., Belytschko, T., Lee, S.H. (2007b) Meshfree collocation method with intrinsic enrichment for interface problems, Computational Mechanics, 40(6), pp.1037-1052 https://doi.org/10.1007/s00466-007-0162-1
  14. Lee, S. H, Yoon, Y.C. (2004) Meshfree point collocation method for elasticity and crack problem, International Journal for Numerical Methods in Engineering, 61, pp.22-48 https://doi.org/10.1002/nme.1053
  15. LeVeque, R.J., Li, Z. (1994) The immersed interface method for elliptic equations with discontinuous coefficients and singular sources, SIAM J. Numer. Anal., 31, pp.1019-1044 https://doi.org/10.1137/0731054
  16. Luo, Y., Häussler-Combe, U. (2002) A generalized finite-difference method based on minimizing global residual, Computer Methods in Applied Mechanics and Engineering, 191, pp.1421-1438 https://doi.org/10.1016/S0045-7825(01)00331-0
  17. Onate, E., Perazzo, F., Miquel, J. (2001) A finite point method for elasticity problems, Computer Methods in Applied Mechanics and Engineering, 79, pp.2151-2163 https://doi.org/10.1016/S0045-7949(01)00067-0
  18. Tu, C., Peskin, C.S. (1992) Stability and instability in the computation of flows with moving immersed boundaries：a comparison of three methods, SIAM J. SCi. Statist. Comput., 13, pp.1361-1376 https://doi.org/10.1137/0913077