2차원 균열판의 p-적응적 유한요소해석을 위한 정규크리깅 보간법의 적용

Application of Ordinary Kriging Interpolation Method for p-Adaptive Finite Element Analysis of 2-D Cracked Plates

  • 우광성 (영남대학교 건설환경공학부) ;
  • 조준형 (한국전력공사 서울전력구 건설처 토목1부) ;
  • 박미영 (제일엔지니어링 구조부)
  • 발행 : 2006.12.30

초록

이 연구의 목적은 두 가지로 대별할 수 있다. 첫째는, 베리오그램 모델링에 기초를 둔 정규크리깅 보간법의 p-적응적 유한요소법으로의 적용성이다. 둘째는, 수정된 초수렴 팻취복구 기법을 사용한 사후오차평가기와 연계된 계층적 p-체눈 세분화의 적응적 유한요소 과정을 제시하는 것이다. 가중치를 부여한 보간기법중의 하나인 정규크리깅 방법은 가우스 적분점에서의 응력데이타로 부터 소위 준정해를 얻는데 적용된다. 가중치를 동일하게 가정하는 종래의 보간기법과는 달리 실험적 및 이론적 베리오그램을 작성한 후 보간을 위한 가중치를 결정하게 된다. 한편, 적응적 p-체눈 세분화는 해석영역의 각 체눈에서 p-차수를 만족할만한 정확도를 얻을 수 있도록 프로그램내에서 자동으로 사후오차평가를 통해 불균등 또는 선택적으로 증가시킨다. 수정된 초수렴 팻취복구기법을 검증하기 위해 극한치를 사용한 새로운 오차평가기가 제안된다. 제안된 알고리즘의 정당성은 선형탄성파괴역학의 대표적 문제들인 중앙균열판, 일변균열 및 양변균열 해석을 통해 테스트되었다.

This paper comprises two specific objectives. The first is to examine the applicability of ordinary kriging interpolation(OK) to the p-adaptivity of the finite element method that is based on variogram modeling. The second objective Is to present the adaptive procedure by the hierarchical p-refinement in conjunction with a posteriori error estimator using the modified S.P.R. (superconvergent patch recovery) method. The ordinary kriging method that is one of weighted interpolation techniques is applied to obtain the estimated exact solution from the stress data at the Gauss points. The weight factor is determined by experimental and theoretical variograms for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. To verify the performance of the modified S.P.R. method, the new error estimator based on limit value has been proposed. The validity of the proposed approach has been tested with the help of some benchmark problems of linear elastic fracture mechanics such as a centrally cracked panel, a single edged crack, and a double edged crack.

키워드

참고문헌

  1. 김영인, 우광성, 이채규(1995) 선형탄성 파괴역학에서 J-적분법에 의한 p-Version 균열모델의 적용, 한국전산구조공학회, 8(4), pp.137-148
  2. 박병기, 신영식, 윤영필, 우광성(1993) 인장력을 받는 균열판의 응력확대계수 산정을 위한 P-version 균열모델, 한국전산구조공학 논문집, 6(4), pp.57-66
  3. 안재석, 우광성, 조준형(2003) 사후오차평가에 의한 적응적 p-체눈 세분화, 대한토목학회 논문집, 23(2A), pp.177-185
  4. 우광성, 이동진, 조준형 (2006) 적응적 p-Version 유한요소법에서 정규크리깅에 의한 응력복구기법, 대한토목학회 논문집. 26(4A), pp.677-687
  5. 최종근(2002) 공간정보 모델링 – 크리깅과 최적화기법, pp.83-163, 구미서관
  6. Abani Patra, Oden, J. T. (1995) A parallel adaptive strategy for hp finite element computations, Comput. Meth. Appl. Mech. Engrg., pp.449-470
  7. Atkinson, P.M., Lloyd. C. D. (2001) Assessing uncertainty in estimates with ordinary and indicator kriging, Computer & Geoscience, 27. pp.929-937 https://doi.org/10.1016/S0098-3004(00)00132-1
  8. Babuska, I., Szabo, B. (1991) Finite Element Analysis, Wiley, New York
  9. Chen, Y., Jiao. X. (2001) Semivariogram fitting with linear programming: Short Note. Computer & Geoscience, 27, pp.71-76 https://doi.org/10.1016/S0098-3004(00)00022-4
  10. Dai, K. Y., Gu, Y. T., Lim, K. M., Liu, G. R.(2003) Comparison between the radial point interpolation and the kriging interpolation used in meshfree method, Computational Mechanics, 32. pp.60-70 https://doi.org/10.1007/s00466-003-0462-z
  11. Goldberger, A. S. (1962) Best linear unbiased prediction in the generalized linear regression, J. Amer. Statist. Assoc., 57, pp.369-375 https://doi.org/10.2307/2281645
  12. Jian, X., Olea, R. A., Yu, Y. S. (1996) Semivariogram modeling by weighted least squares, Computer & Geosciences, 22(4), pp.387-397 https://doi.org/10.1016/0098-3004(95)00095-X
  13. Krige, D.G. (1951) A statiscal approach to some basic mine valuation problems on the Witwatersrand, Journal of the Chemical, Metallurgical and Mining Society of South Africa, 52, pp.119-139
  14. Matheron, G. (1973) The intrinsic random functions and their applications, J. Appl. Probab., 5, pp.439 -468 https://doi.org/10.2307/1425829
  15. Stein, M. L. (1999) Interpolation of Spatial Data- Some Theory for Kriging, Springer
  16. Woo, K. S., Hong. C. H., Shin. Y. S.(1998) An extended equivalent domain integral method for mixed mode fracture problems by the p-version of FEM, Int. J. Numer. Meth. Engrg., 42, pp.857- 884 https://doi.org/10.1002/(SICI)1097-0207(19980715)42:5<857::AID-NME390>3.0.CO;2-#
  17. Yamamoto, J. K. (2000) An Alternative Measure of the Reliability of Ordinary Kriging Estimates, Mathematical Geology, 32(4), pp.489-509 https://doi.org/10.1023/A:1007577916868
  18. Zhu, J. Z .. Zienkiewicz, O. C. (1992) The super-convergent patch recovery and a posteriori error estimate Part I : The recovery Technique, Int. J. Numer. Meth. Engrg., 33, pp.1331-1364 https://doi.org/10.1002/nme.1620330702
  19. Zhu, J. Z., Zienkiewicz, O. C.(1992) The super-convergent patch recovery and a posteriori error estimate Part II : Error Estimate and Adaptivity, Int. J. Numer. Meth. Engrg., 33, pp.1365-1382 https://doi.org/10.1002/nme.1620330703