• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.729-740
• /
• 1994

In the shape optimal design based on h-version of FEM, the ideal mesh for the initial geometry most probably will not be suitable for the final analysis. Thus, it is necessary to remesh the geometry of the model at each stage of optimization. However, the p-version of FEM appears to be a very attractive alternative for use in shape optimization. The main advantages are as follows; firstly, the elements are not sensitive to distortion for interpolation polynomials of order $p{\geq}3$; secondly, even singular problems can be solved more efficiently with p-version than with the h-version by proper mesh design; thirdly, the initial mesh design are identical. The 2-D p-version model for shape optimization is presented on the basis of Bezier's curve fitting, gradient projection method, and integrals of Legendre polynomials. The numerical results are performed by p-version software RASNA.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.3
• /
• pp.329-336
• /
• 2000

A method for automatic extraction of 2-dimensional finite element connectivities by searching the shortest closed path from a certain node to the starting node was developed. Only the best path among the possible paths was probed. The uniqueness and validity of the extracted path were examined. The proposed method was proved to be complete. Examples show that the proposed method can extract elements exactly from the irregular mesh which can not be handled easily by the conventional automatic mesh generation.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.165-172
• /
• 1996

For thin elastic structures such as beams, arches, plates and shells, there may exist the boundary layer in the narrow thin region neighborhood of boundaries, where the solution displays the singular behavior exponentially decaying in the normal direction to the boundary. In the finite element analysis of these structures, finite element mesh patterns have a significant role to capture this singularity. This paper introduces the analytic study of this problem and provides a guideline to construct optimal mesh patterns together with numerical experiments.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.4 s.74
• /
• pp.429-440
• /
• 2006

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.

• KOREAN JOURNAL OF CROP SCIENCE
• /
• v.54 no.4
• /
• pp.366-374
• /
• 2009

To determine the optimum operating conditions of the indented cylinder length grader based on varietal characteristics, broken rice removal capabilities were measured on the milled rice of 41 varieties at different conditions consisted of 3 types of indented cylinders with different indent-opening diameters and 3 different collecting-angles of broken rice. The broken rice removal capabilities were swayed by the indent's opening diameter and depth as well as the angle of the collecting trough of broken rice on the point of instrument, and by the kernel length and width as well as 1000 grain weight of milled rice on the point of rice variety. When the angle of broken rice collecting-trough reached to near the horizontal center line of the indented cylinder on the direction of upward turning side, which was referred $0^{\circ}$ in this paper, the amount of rice collected in trough increased, whereas the loss of head rice also increased. Considering the removal rate of broken rice as well as loss of head rice, it was thought that the suitable angle of trough for broken rice collecting was located $5^{\circ}$ to $15^{\circ}$ depending on varietal characteristic and indent opening diameter. It was thought that 4.2 mm or more of indent opening diameter was recommendable for the rice varieties having heavier 1000 grain weight than 22.3g, as well as larger sizes than 2.9 and 5.2 mm in width and length of rice kernel, respectively; 3.8 mm for the small-sized thin kernels, and a proper diameter between 3.8 and 4.2 mm for short to middle kernels. The varieties with relatively shorter length compared to width of kernel were more difficult to separate the broken rice than the opposite ones. For effective separation of that, it seems that some specific indent shapes such as wider opening and shallow depth etc. are required. When the broken rice content were excessively high, wider diameters of indent openings than specified sizes were thought to be the better.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.533-541
• /
• 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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1993.04a
• /
• pp.3-10
• /
• 1993

A new three-dimensional transition solid elements was presented for the automated three-dimensional adaptive h-refinement where the steep stress gradient exists. To be consistent with 8-node solid element with nonconforming modes in accuracy, these transition elements were improved through the addition of the associated nonconforming modes. Numerical examples show that the performance of the element and the applicability to 3D adaptations are satisfactory.

• Computational Structural Engineering
• /
• v.11 no.1
• /
• pp.90-95
• /
• 1998

본 고에서는 비압축성 점성유체의 유한요소해석 기법을 소개하였다. 대류항의 상류화 기법으로 안정된 해를 도출할 수 있으며 Penalty 방법에 기반하여 압력항을 지배방정식으로부터 소거함으로써 해석시간과 요구저장공간을 감소시켰다. 실린더 주변의 유동장을 해석하여 와의 방출을 성공적으로 묘사하였으며 항력계수를 17%정도의 오차로 계산하였다. 적응적 요소세분화 기법에 대한 연구를 통해 적절한 오차평가 기법 및 최적의 체눈을 형성하는 기법을 제시하였다. 또한 동적 해석에 적합한 요소재결합 알고리즘에 대한 연구가 진행중이다. 본 고의 결과는 직접적으로 풍공학분야에 사용하기에는 아직 계산 시간의 효율성이나 해의 정확도 및 안정성면에서 무리가 있으나 추가적인 연구를 통하여 해석기법의 개선을 도모하고 컴퓨터 등 계산장비의 급속한 발전으로 장래에 경쟁력을 획득할 수 있을 것으로 기대된다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.3
• /
• pp.331-340
• /
• 2010

The basic theory and application of new adaptive finite element algorithm have been proposed in this study including the adaptive hp-refinement strategy, and the effective method for constructing hp-approximation. The hp-adaptive finite element concept needs the integrals of Legendre shape function, nonuniform p-distribution, and suitable constraint of continuity in conjunction with irregular node connection. The continuity of hp-adaptive mesh is an important problem at the common boundary of element interface. To solve this problem, the constraint of continuity has been enforced at the common boundary using the connectivity mapping matrix. The effective method for constructing of the proposed algorithm has been developed by using hierarchical nature of the integrals of Legendre shape function. To verify the proposed algorithm, the problem of simple cantilever beam has been solved by the conventional h-refinement and p-refinement as well as the proposed hp-refinement. The result obtained by hp-refinement approach shows more rapid convergence rate than those by h-refinement and p-refinement schemes. It it noted that the proposed algorithm may be implemented efficiently in practice.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.4A
• /
• pp.677-687
• /
• 2006

Kriging interpolation is one of the generally used interpolation techniques in Geostatistics field. This technique includes the experimental and theoretical variograms and the formulation of kriging interpolation. In contrast to the conventional least square method for stress recovery, kriging interpolation is based on the weighted least square method to obtain the estimated exact solution from the stress data at the Gauss points. The weight factor is determined by variogram modeling for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. In addition to this, the p-level is increased non-uniformly or selectively through a posteriori error estimation based on SPR (superconvergent patch recovery) technique, proposed by Zienkiewicz and Zhu, by auto mesh p-refinement. The cut-out plate problem under tension has been tested to validate this approach. It also provides validity of kriging interpolation through comparing to existing least square method.