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

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

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 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.

• Structural Engineering and Mechanics
• /
• 제11권5호
• /
• pp.535-546
• /
• 2001

A new p-adaptive analysis scheme for hp-clouds method is presented. In the scheme, refined global equations are resolved into two parts, one of them being related to the newly appended dof's. The solution obtained in previous analysis step is reflected in the force vector. The size of the p-adaptive equation consisting of the newly appended dof's is much smaller than the original equation. Consequently, the computational cost is drastically decreased. Through numerical examples, the efficiency and efficacy of the method in comparison with the existing p-refinement scheme of the hp-clouds have been demonstrated.

• Structural Engineering and Mechanics
• /
• 제39권6호
• /
• pp.767-793
• /
• 2011

An efficient edge-based smoothed finite element method (ES-FEM) has been recently developed for solving solid mechanics problems. The ES-FEM uses triangular elements that can be generated easily for complicated domains. In this paper, the complexity study of the ES-FEM based on triangular elements is conducted in detail, which confirms the ES-FEM produces higher computational efficiency compared to the FEM. Therefore, the ES-FEM offers an excellent platform for adaptive analysis, and this paper presents an efficient adaptive procedure based on the ES-FEM. A smoothing domain based energy (SDE) error estimate is first devised making use of the features of the ES-FEM. The present error estimate differs from the conventional approaches and evaluates error based on smoothing domains used in the ES-FEM. A local refinement technique based on the Delaunay algorithm is then implemented to achieve high efficiency in the mesh refinement. In this refinement technique, each node is assigned a scaling factor to control the local nodal density, and refinement of the neighborhood of a node is accomplished simply by adjusting its scaling factor. Intensive numerical studies, including an actual engineering problem of an automobile part, show that the proposed adaptive procedure is effective and efficient in producing solutions of desired accuracy.

• International Journal of Aeronautical and Space Sciences
• /
• 제18권4호
• /
• pp.601-613
• /
• 2017

An adjoint-based error estimation and mesh adaptation study is conducted for two-dimensional viscous flows on unstructured hybrid meshes. The error in an integral output functional of interest is estimated by a dot product of the residual vector and adjoint variable vector. Regions for the mesh to be adapted are selected based on the amount of local error at each nodal point. Triangular cells in the adaptive regions are refined by regular refinement, and quadrangular cells near viscous walls are bisected accordingly. The present procedure is applied to single-element airfoils such as the RAE2822 at a transonic regime and a diamond-shaped airfoil at a supersonic regime. Then the 30P30N multi-element airfoil at a low subsonic regime with a high incidence angle (${\alpha}=21deg.$) is analyzed. The same level of prediction accuracy for lift and drag is achieved with much less mesh points than the uniform mesh refinement approach. The detailed procedure of the adjoint-based mesh refinement for the multi-element airfoil case show that the basic flow features around the airfoil should be resolved so that the adjoint method can accurately estimate an output error.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2003년도 봄 학술발표회 논문집
• /
• pp.288-295
• /
• 2003

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 recovery technique. In case of the recovery technique, the SPR(superconvergent patch recovery) approach has been modified for p-adaptive mesh refinement. 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. To verify the proposed algorithm, the limit value approach is proposed which utilizes the exact strain energy computed from the extrapolation equation. A new pre-processor is developed for the p-version finite element program in which the vector graphic editor is used for the automatic generation of node connection and coordinate by halfedge solid data structure according to uniform or nonuniform p-distribution. The general 2-D algorithm is also developed to generate face modes and internal modes in accordance with different mesh types. The quality of the error estimator is investigated with the help of two mumerical examples. The results show that the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

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

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

• 한국항만학회지
• /
• 제10권2호
• /
• pp.61-68
• /
• 1996

The main problems in structural analysis by Finite Eelement Method are difficulty in making data file and error estimation. For decreasing these problems' pays. have been suggesting the adaptive mesh refinement and error estimation method. Posteriory error estimation methods suggested by Jang[1], Babuska[2,3], Ohtsubo[8,9], and this paper. Comparing these methods and examine their properties. According this paper, In the problem supposed having singularity, the method suggested by this paper is good, But the problem supposed having no singularity, the method suggested by Jang[1] is good. For decreasing the effect of initial mesh in p-refinement, make application h-refinement at first and apply p-refinement, and confine polynomial's degree to two, for making program simply by plural mesh models are not needed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.219-222
• /
• 2006

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

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