• Structural Engineering and Mechanics
• /
• 제17권3_4호
• /
• pp.379-391
• /
• 2004

This paper deals with the development of h-version adaptive mesh refinement and recovery strategy using variable-node elements and its application to various engineering field problems with 2D quadrilateral and 3D hexahedral models. The variable-node elements which have variable mid-side nodes on edges or faces are effectively used in overcoming some problems in connecting the different layer patterns of the transition zone between the refined and coarse mesh. A modified recovery technique of gradients adequate for variable-node elements and proper selection of error norms for each engineering field problems are proposed. In the region in which the error is greater than the permissible refinement error, the mesh is locally refined by subdivision. Reversely, in some parts of the domain having the error smaller than the permissible recovery error, the mesh is locally recovered (coarsened) by combination. Hierarchical structures (e.g. quadtrees and octrees) and element-based storage structures are composed to perform this adaptive process of refinement and recovery. Some numerical examples of a 3D heat conduction analysis of the concrete with hydration heat and a 2D flow analysis of vortex shedding show effectiveness and validity of the proposed scheme.

• Wind and Structures
• /
• 제1권1호
• /
• pp.111-125
• /
• 1998

This paper deals with the development of variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined element and efficiently used for the construction of a refined mesh without generating distorted elements. A modified Guassian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of the independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems in which the locations to be refined are changed in accordance with the dynamic distribution of velocity gradient, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

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

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2000년도 춘계학술발표대회 논문집
• /
• pp.89-92
• /
• 2000

In this paper, an effective procedure is presented for the local recovery of displacements and stresses in multilayered composite panels, which incorporate the local refinement using mesh superposition. The mesh superposition method is used to refine the global coarse mesh by superimposing refined mesh to the localized zone of interest without transition zones. The finite element model used is a solid element based on the Hellinger-Reissner variational principle. The a posteriori computation of the through-the-thickness distributions of displacements and stresses is achieved using a predictor-corrector procedure. The procedure utilizes the superconvergent stresses and nodal displacements of the finite element patch. The element patch is generated by locally superimposing a refined local mesh to the coarse global mesh.

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

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

• 대한기계학회논문집A
• /
• 제25권6호
• /
• pp.1004-1014
• /
• 2001

The fully automatic algorithm from initial finite element mesh generation to remeshing in two dimensional geometry is introduced using bubble packing method (BPM) for finite element analysis. BPM determines the node placement by force-balancing configuration of bubbles and the triangular meshes are made by Delaunay triangulation with advancing front concept. In BPM, we suggest two node-search algorithms and the adaptive/recursive bubble controls to search the optimal nodal position. To use the automatically generated mesh information in FEA, the new enhanced bandwidth minimization scheme with high efficiency in CPU time is developed. In the remeshing stage, the mesh refinement is incorporated by the control of bubble size using two parameters. And Superconvergent Patch Recovery (SPR) technique is used for error estimation. To verify the capability of this algorithm, we consider two elasticity problems, one is the bending problem of short cantilever beam and the tension problem of infinite plate with hole. The numerical results indicate that the algorithm by BPM is able to refine the mesh based on a posteriori error and control the mesh size easily by two parameters.

• Structural Engineering and Mechanics
• /
• 제3권4호
• /
• pp.325-340
• /
• 1995

The iterative procedure to use the nonconforming elements in elasto-plastic problems is established and applied to the variable node transition solid element developed for the automated three-dimensional local mesh refinement. Through numerical tests, the validity and performance of the element are examined. As the nonlinear iterative procedure presented in this paper is accomplished for the general three-dimensional case, it can also be easily applied to the two-dimensional elements such as membranes, plates and shells.

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

The finite element linear buckling analysis of folded plate structures using adaptive h-refinement methods is presented in this paper. The variable-node flat shell element used in this study possesses the drilling D.O.F. which, in addition to improvement of the element behavior, permits an easy connection to other elements with six degrees of freedom per node. The Box-typed structures can be analyzed using these developed flat shell elements. By introducing the variable node elements some difficulties associated with connecting the different layer patterns, which are common in the adaptive h-refinement on quadrilateral mesh, can be overcome. To obtain better stress field for the error estimation, the super-convergent patch recovery is used. The convergent buckling modes and the critical loads associated with these modes can be obtained.

• 대한토목학회논문집
• /
• 제26권4A호
• /
• pp.677-687
• /
• 2006

크리깅 보간법은 지구통계학 분야에 주로 사용되는 보간법의 하나이다. 이 방법은 실험적 베리오그램과 이론적 베리오그램의 작성과 크리깅 보간법의 정식화에 관한 연구를 포함하고 있다. 종래의 응력복구를 위한 최소제곱법과 대조적으로, 가우스적분점에서의 응력데이타로부터 준정해를 얻기 위해 가중 최소제곱법에 기초를 둔다. 즉, 동일한 가중치를 사용하는 종래의 방식들과는 달리 가우스적분점에서의 응력값의 보간을 위하여 베리오그램 모델링을 통한 가중치가 결정된다. 한편, 분할된 요소망에 Zienkiewicz와 Zhu에 의해 제안된 SPR기법에 기초를 둔 사후오차평가를 통해 p-차수를 균등 또는 선택적으로 증가시키는 자동체눈 방식이 도입되었다. 이 방법의 정당성을 보기위해 인장력을 받는 개구부를 갖는 평판문제를 해석하였다. 또한, 기존의 최소제곱법과의 비교를 통한 크리깅보간법의 정당성을 보여 주었다.