DOI QR코드

DOI QR Code

무요소법을 이용한 균열진전 문제의 형상 최적설계

Shape Design Optimization of Crack Propagation Problems Using Meshfree Methods

  • 김재현 (서울대학교 아이소-지오메트릭 최적설계 창의연구단) ;
  • 하승현 (한국해양대학교 해양공학과) ;
  • 조선호 (서울대학교 아이소-지오메트릭 최적설계 창의연구단)
  • Kim, Jae-Hyun (National Creative Research Initiatives(NCRI) Center for Isogeometric Optimal Design, Department of Naval Architecture and Ocean Engineering, Seoul National University) ;
  • Ha, Seung-Hyun (Department of Ocean Engineering, Korea Maritime and Ocean University) ;
  • Cho, Seonho (National Creative Research Initiatives(NCRI) Center for Isogeometric Optimal Design, Department of Naval Architecture and Ocean Engineering, Seoul National University)
  • 투고 : 2014.03.04
  • 심사 : 2014.09.24
  • 발행 : 2014.10.31

초록

본 논문에서는 재생 커널 기법을 사용하여 혼합모드 균열진전 문제에 대한 연속체 기반의 형상 설계민감도 해석을 수행하였다. 재생 커널 기법은 기존의 유한요소법과 달리 요소망을 재구성할 필요가 없어, 커널 함수의 연속성을 증가시켰을 때 높은 정밀도의 형상함수를 얻을 수 있다는 장점을 가지고 있다. 균열선단 주변에서 J-적분을 수행하기 위해 선형탄성 조건이 고려되었다. 변위장과 응력 확대 계수의 설계변수에 대한 감도해석을 위하여 물질도함수를 도입하였으며 직접 미분법보다 효율적인 애조인 방법을 사용하여 설계민감도를 유도하였다. 수치 예제들을 통해서 재생 커널 기법을 이용한 균열진전 해석결과의 타당성을 확인하였으며 애조인 방법을 이용한 형상 설계민감도 해석 결과를 유한차분법과 비교하여 매우 정확하고 효율적인 결과를 얻을 수 있음을 알 수 있었다. 이를 바탕으로 간단한 모델에 대하여 형상 최적설계를 수행하여 균열이 발생될 수 있는 구조물에 대해서 균열에 의한 피해를 최소화할 수 있도록 균열을 제어할 수 있는 최적의 형상을 도출하였다.

This paper presents a continuum-based shape design sensitivity analysis(DSA) method for crack propagation problems using a reproducing kernel method(RKM), which facilitates the remeshing problem required for finite element analysis(FEA) and provides the higher order shape functions by increasing the continuity of the kernel functions. A linear elasticity is considered to obtain the required stress field around the crack tip for the evaluation of J-integral. The sensitivity of displacement field and stress intensity factor(SIF) with respect to shape design variables are derived using a material derivative approach. For efficient computation of design sensitivity, an adjoint variable method is employed tather than the direct differentiation method. Through numerical examples, The mesh-free and the DSA methods show excellent agreement with finite difference results. The DSA results are further extended to a shape optimization of crack propagation problems to control the propagation path.

키워드

참고문헌

  1. Chen, J.S., Pan, C., Wu, C.T. (1997) Large Deformation Analysis of Rubber Based on a Reproducing Kernel Particle Method, Comput. Mech.., 19, pp.211-217. https://doi.org/10.1007/s004660050170
  2. Ha, S.H., Cho, S. (2008) Level Set Based Topological Shape Optimization of Geometrically Nonlienar Structures using Unstructured Mesh, Comput. & Struct., 86, pp.1447-1455. https://doi.org/10.1016/j.compstruc.2007.05.025
  3. Ha, Y., Cho, S. (2006) Design Sensitivity Analysis and Topology Optimization of Eigenvalue Problems for Piezoelectric Resonators, Smart Mater. & Struct., 15, pp.1513-1524. https://doi.org/10.1088/0964-1726/15/6/002
  4. Haber, R.B., Koh, H.M. (1985) Explicit Expressions for Energy Release Rates Using Virtual Crack Extensions, International, J. Numer. Method in Eng., 21, pp.301-315. https://doi.org/10.1002/nme.1620210208
  5. Haug, E.J., Choi, K.K., Komkov, V. (1986) Design Sensitivity Analysis of Structural System, Academic Press, New York, NY.
  6. Kim, M.-G., Ha, S.H., Cho, S. (2012) Level Set Based Topological Shape Optimization of Hyper-Elastic Nonlinear Structures using Topological Derivatives, J. Comput. Struct. Eng. Inst. Korea, 25(6), pp.559-567. https://doi.org/10.7734/COSEIK.2012.25.6.559
  7. Park, S.M. (2004) Crack Propagation Analysis using Reproducing Kernel Method, Master Thesis, Seoul National University.