• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• pp.481-486
• /
• 2007

The Zienkiewicz-Zhu(Z/Z) error estimate is slightly modified for the hierarchical p-refinement, and is then applied to L-shaped plates subjected to bending to demonstrate its effectiveness. An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the superconvergent patch recovery(SPR) technique. The modified Z/Z error estimate p-refinement is different from the conventional approach because the high order shape functions based on integrals of Legendre polynomials are used to interpolate displacements within an element, on the other hand, the same order of basis function based on Pascal's triangle tree is also used to interpolate recovered stresses. The least-square method is used to fit a polynomial to the stresses computed at the sampling points. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly or selectively. It is noted that the error decreases rapidly with an increase in the number of degrees of freedom and the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

• 한국전산구조공학회논문집
• /
• 제23권3호
• /
• pp.331-340
• /
• 2010

적응적 hp-세분화 기법과 그 기법의 효과적인 구성방법을 포함한 새로운 적응적 유한요소 알고리즘의 기초이론 및 적용이 이 연구를 통해 제시되었다. 적응적 hp-세분화 기초의 유한요소기법은 적분형 르장드르 형상함수와 요소별로 불균등한차수의 분배 및 비정형적인 절점연결과 관련된 연속조건을 만족시킬 수 있는 제약조건을 필요로 한다. 따라서 요소간의 접합부분에서 적응적 hp-유한요소망의 연속성이 중요한 문제로 대두된다. 이러한 문제를 요소경계에 연속성 제약조건을 절점연결 사상행렬을 적용하여 해결하였다. 또한, 적분형 르장드르 형상함수의 계층성질을 이용하여 제시된 알고리즘의 효율적 정식화 방안을 제시하였다. 간단한 캔틸레버문제가 h-세분화, p-세분화 그리고 hp-세분화 방법에 의해 계산되었다. hp-세분화의 결과는 다른 방식의 세분화에 비해 보다 빠른 수렴성을 보여 주는 것이 확인되었다. 그러므로 제시된 hp-세분화 알고리즘은 실제문제에 효율적으로 적용될 수 있을 것으로 생각된다.

• 한국공간구조학회논문집
• /
• 제8권4호
• /
• pp.81-90
• /
• 2008

이 논문에서는 적분형 르장드르 다항식을 사용한 3차원 계층적 고체요소의 유한요소 정식화를 보여준다. 제안하는 육면체 고체요소는 절점, 변, 면, 그리고 내부모우드를 포함한은 4개의 서로 다른 모우드로 구성되어 있다. 영에너지 모우드와 일정변형률 조건을 확인하기 위해 고유치 시험과 조각시험이 수행되었다. 여기에 추가되어, 적응적 p-유한요소해석을 위해 유한요소해석으로부터 구한 후처리 응력값의 평활화에 기초를 둔 사후오차평가 기법이 연구된다. 자유도가 증가함에 따라 수렴속도측면에서 균등 p-분배와 불균등 p-분배에 의한 유한요소해의 차이점이 비교된다. 제안된 요소의 성능을 보이기 위해 간단한 캔틸레버보가 테스트되었다.

• 한국전산구조공학회논문집
• /
• 제20권5호
• /
• pp.533-541
• /
• 2007

계층적 p-세분화를 위해 Zienkiewicz-Zhu 오차평가법이 약간 수정되었으며, 이 방법의 유효성을 보이기 위해 휨을 받는 개구부를 갖는 Reinssner-Mindlin $C^{\circ}$-평판에 적용하였다. 유한요소해석상의 적응적 체눈을 결정하는 단계는 초수렴 팻취 복구기법에 기초를 둔 사후오차평가자와 연계된 p-세분화에 의해 제안되었다. 요소내의 변위장을 정의하기 위해 적분형 르장드르 고차 형상함수가 사용되는 반면 복구응력을 보간하기 위해 파스칼의 삼각수에 기초를 둔 같은 차수의 고차다항식이 사용되는 이유로 수정 Z/Z 오차평가는 종래의 방법과 다소 차이를 보여준다. 가우스 적분점에서의 응력을 최적화하기 위해 필요한 다항식으로 표현되는 응력함수를 얻기 위해 최소제곱법이 사용되었다. 고정된 요소망에 거의 최적의 형상함수 차수의 분배를 찾기 위한 전략이 논의되었는데, 허용되는 정확도를 얻을 수 있을 때까지 각 요소마다 형상함수의 차수를 불균등하게 증가시키는 방법으로, 소위 최적의 선택적 p-분배를 자동으로 결정하도록 되어있다. 위의 사항들을 L-형 평판 해석에 적용한 결과, 적응적 p-체눈설계 단계가 진행됨에 따라 자유도의 증가에 따라 오차량은 급격히 감소되는 것을 알 수 있었고, 제안된 오차 지시자에 의한 적응적 p-체눈 세분화는 최적 p-분배 진행방향에 근접하는 것을 볼 수 있었다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제17권3호
• /
• pp.139-170
• /
• 2013

We extend a p-hierarchical decomposition of the second degree finite element space of N$\acute{e}$d$\acute{e}$lec for tetrahedral meshes in three dimensions given in [1] to meshes with hexahedral elements, and derive p-hierarchical decompositions of the second degree finite element space of Raviart-Thomas in two dimensions for triangular and quadrilateral meshes. After having proved stability of these subspace decompositions and requiring certain saturation assumptions to hold, we construct a local a posteriori error estimator for fem and bem coupling of a time-harmonic electromagnetic eddy current problem in $\mathbb{R}^3$. We perform some numerical tests to underline reliability and efficiency of the estimator and test its usefulness in an adaptive refinement scheme.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1992년도 가을 학술발표회 논문집
• /
• pp.39-44
• /
• 1992

This paper is concerned with formulations of the hierarchical $C^{o}$-plate element on the basis of Reissner-Mindlin plate theory. On reason for the development of the aforementioned element is that it is still difficult to construct elements based on h-version concepts which are accurate and stable against the shear locking effects. An adaptive mesh refinement and selective p-distribution of the polynomial degree using hp-version of the finite element method we proposed to verify the superior convergence and algorithmic efficiency with the help of the clamped L-shaped plate problems.s.

• 대한토목학회논문집
• /
• 제12권4_1호
• /
• pp.67-76
• /
• 1992

축대칭(軸對稱) 선형강성(線形彈性) 응력해석을 위해 p-version 유한요소법에 기초한 계층적(階層的) 정식화 과정이 제안되었다. 이 방식은 적분형 르장드르 다항식을 사용하여 절점좌표값을 갖지 않는 절점을 추가하여 형상함수의 조합형태로 변위함수(變位)를 근사시키는 방법이다. 형상함수(形狀函數)가 계층적 성질을 갖기 때문에 강성도(剛性度)행렬과 하중벡터도 계층적이 된다. 본 연구에서 제안된 요소(要所)의 장점(長點)은 다음과 같다. 첫째, 개선된 수치연산의 효율성이며 둘째, 요소간에 서로 다른 차수(次數)의 형상함수를 사용할 수 있고 셋째, p-세분화를 할 때 저차(低次)일 때 계산된 값을 그대로 사용할 수 있다. 수치예제를 통해 제안된 요소의 정확도(正確度), 효율성(效率性), 모델링의 간편성(簡便性), 적용성(適用性) 및 변위와 응력 그리고 에너지 Norm등을 사용하여 그 우월성을 입증하고 있다. 몇 가지 예제의 해석결과는 이미 발표된 논문과 아울러 해석적 방법에 의한 결과와 비교되었다.

• 한국전산구조공학회논문집
• /
• 제19권4호
• /
• pp.429-440
• /
• 2006

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

• 대한토목학회논문집
• /
• 제13권2호
• /
• pp.151-160
• /
• 1993

본 논문에서는 Reissner-Mindlin 평판이론에 근거한 계층적 $C^{\circ}$-평판요소가 제안되었다. 적분형 르장드르 형상함수에 근거한 계층요소를 제안하는 이유는 종래의 h-version 유한요소법의 개념 을 사용하여 전단구속 효과등에 대한 해의 정확도 및 수치안정성을 확보할 수 있는 요소를 만드는데 여전히 어려움이 수반되기 때문이다. 적응적 체눈 p-세분화와 선택적 형상함수 차수 p의 분포를 사용하는 hp-version 유한요소법을 사용하여 내부주변은 자유단의 개구부를 갖고, 외부주변이 단순지지된 L-형 평판해석을 수행하였는데 종래의 h-version 유한요소법에 비해 우월한 수렴성과 전단구속을 피할 수 있는 등의 알고리즘 효율성을 보여 주고 있다.