• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.1
• /
• pp.43-50
• /
• 1992

The purpose of this study is to ascertain the robustness of p-version model with hierarchic intergrals of Legendre shape functions in various applications including plane stress/strain, axisymmetric and shell problems. The most important symptoms of accuracy failure in modern finite elements are spurious mechanisms and a phenomenon known as locking which are exhibited for incompressible materials and irregular shapes which contain aspect ratios(R/t, a/b), tapered ratio(d/b), and skewness. The condition numbers and energy norms are used to estimate numerical errors, convergence characteristics and algorithmic efficiencies for verifying the aforementioned symptoms of accuracy failure. Numerical results from p-version models are compared with those from NASTRAN, SAP90, and Cheung's hybrid elements.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1289-1298
• /
• 1994

The new p-version crack model is proposed to estimate the stress intensity factors of the thick cracked plate under flexure. The proposed model is based on high order theory and $C^{\circ}$-plate element including shear deformation. The displacements fields are defined by integrals of Legendre polynomials which can be classified into three groups such as basic mode, side mode and internal mode. The computer implementation allows arbitrary variations of p-level Up to a maximum value of 10. The stress intensity factors are computed by virtual crack extention approach. The effects of ratios of thickness to crack length(h/a), crack length to width(a/W) and boundary conditions are investigated. Very good agreement with the existing solution in the literature are shown for the uncracked plate as well as the cracked plate.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1991.10a
• /
• pp.75-80
• /
• 1991

The purpose of this study is to ascertain the robustness of p-version model with hierarchic intergrals of Legendre shape functions in various applications including plane stress/strain, axisymmetric and shell problems. The most important symptoms of accuracy failure in modern finite elements are spurious mechanisms and a phenomenon known as locking which are exhibited for incompressible materials and irregular shapes which contain aspect ratios(R/t, a/b), tapered ratio(d/b), and skewness. The condition numbers and energy norms are used to estimate numerical errors, convergence characteristics and algorithmic efficiencies for verifying the aforementioned symptoms of accuracy failure. Numerical results from p-version models are compared wi th those from NASTRAN, SAP90, and Cheung's hybrid elements.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.4 s.35
• /
• pp.82-97
• /
• 1999

A deep sea vehicle must be designed to ensure its safety under ultra-high pressure circumstances. If a pressure housing of a deepsea vehicle is collapsed by ultra-high pressure, the deepsea vehicle may be lost. The objective of this paper is to introduce a design collapse pressure for the deep sea pressure vessel which is composed of one cylinder and two hemispheres. Especially the collapse pressure of hemispherical shell with a hole at top is analyzed by a variational approach (weighted residual method). And for the purpose of design, the salty factor of collapse pressure is presented which is analyzed by interpolation method.

• Communications for Statistical Applications and Methods
• /
• v.17 no.6
• /
• pp.811-827
• /
• 2010

This paper deals with a density estimation method in binary choice models that can be regarded as a statistical inverse problem. We use an orthogonal basis to estimate density function and consider the choice of an appropriate truncation parameter to reflect the model complexity and the prediction accuracy. We propose a data-dependent rule to choose the truncation parameter in the context of binary choice models. A numerical simulation is provided to illustrate the performance of the proposed method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.151-160
• /
• 1993

This paper is concerned with formulations of the hierarchical $C^{\circ}$-plate element on the basis of Reissner-Mindlin plate theory. On reason for the development of the aforementioned element based on Integrals of Legendre shape functions 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 are proposed to verify the superior convergence and algorithmic efficiency with the help of the simply supported L-shaped plate problems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.24 no.2
• /
• pp.183-191
• /
• 2006

The accuracy of the upward continuation is assessed through the gravity modeling using an ultra-high degree spherical harmonic expansion. The difficulties in the numerical calculation of Legendre function with ultra-high degree, underflow and/or overflow, is successfully resolved in 128 bit calculation scheme. Using the generated Legendre function, the gravity anomaly with spatial resolution of $1'{\times}1'$ on the geoid is calculated. The generated gravity anomaly is degraded and extracted with various noise levels and data intervals, then upward continuation is applied to each data sets. The comparison between the upward continued gravity disturbances and the directly calculated from the spherical harmonics showed that the accuracy on the direct method was significantly better than that of Poisson method. In addition, it is verified that the denser and less noised gravity data on the geoid generates better gravity disturbance vectors at an altitude. Especially, it is found that the gravity noise level less than 5mGal, and the data interval less than 2arcmin is necessary for next generation precision INS navigation which requires the accuracy of 5mGal or better at an altitude.