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

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.445-454
• /
• 1991

In this thesis, Mindlin plate element with nine nodes and three degrees-of-freedom at each node is formulated and is employed in eigen-analysis of a rectangular plates in order to alleviate locking phenomenon of eigenvalues. Eigenvalues and their modes may be locked if conventional $C_{0}$-isoparametric element is used. In order to reduce stiffness locking phenomenon, two methods (1, the general reduced and selective integration, 2, the new element that use of modified shape function) are studied. Additionally in order to reduce the error due to mass matrix, two mass matrixes (1, Gauss-Legendre mass matrix, 2, Gauss-Lobatto mass matrix) are considered. The results of eigen-analysis for two models (the square plate with all edges simply-supported and all edges built-in), computed by two methods for stiffness matrix and by two mass matrixes are compared with theoretical solutions and conventional numerical solutions. These comparisons show that the performance of the two methods with Gauss-Lobatto mass matrix is better than that of the conventional plate element. But, by considering the spurious rigid body motions, the element which employs modified shape function with full integration and Gauss-Lobatto mass matrix can elevate the accuracy and convergence of numerical solutions.

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

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.78.2-78.2
• /
• 2019

We develop an Alcock-Paczynski (AP) test method that uses the evolution of redshift-space two-point correlation function (2pCF) of galaxies. The method improves the AP test proposed by Li et al. (2015) in that it uses the full two-dimensional shape of the correlation function. Similarly to the original method, the new one uses the 2pCF in redshift space with its amplitude normalized. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF. This is because the 2pCF should not change apart from the expected small non-linear evolution if galaxy clustering is not distorted by incorrect choice of cosmology used to convert redshift to comoving distance. Our new method decomposes the redshift difference of the 2-dimensional correlation function into the Legendre polynomials whose amplitudes are modelled by radial fitting functions. The shape of the normalized 2pCF suffers from small intrinsic time evolution due to non-linear gravitational evolution and change of type of galaxies between different redshifts. It can be accurately measured by using state of the art cosmological simulations. We use a set of our Multiverse simulations to find that the systematic effects on the shape of the normalized 2pCF are quite insensitive to change of cosmology over \Omega_m=0.21 - 0.31 and w=-0.5 - -1.5. Thanks to this finding, we can now apply our method for the AP test using the non-linear systematics measured from a single simulation of the fiducial cosmological model.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.10a
• /
• pp.38-45
• /
• 1994

P-version finite element model for the computation of stress intensity factors in two dimensional cracked panels by J-integral method is presented. The proposed model is based on high order theory and hierarchical shape function. The displacements fields are defined by integrals of Legendre polynomials which can be classified into three part such as basic mode, side mode, integral mode. The stress intensity factors are computed by J-integral method. The example models for validating the proposed p-version model are centrally cracked panel, single and double edged crack in a rectangular panel under pure Mode I. And the analysis results are compared with those by the h-version of FEM and empirical solutions in literatures. Very good agreement with the existing solution are shown.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.04a
• /
• pp.62-69
• /
• 1996

In this study, an equivalent domain integral (EDI) method is presented to estimate the track-till integral parameter, J-value, for two dimensional cracked elastic bodies which may quantify the severity of the crack-tit) stress fields. The conventional J-integral method based on line integral has been converted to equivalent area or domain integrals by using the divergence theorem. It is noted that the EDI method is very attractive because all the quantities necessary for computation of the domain integrals are readily available in a finite element analysis. The details and its implementation are extened to both h-version finite element model with 8-node isoparametric element and p-version finite element model with high order hierarchic element using Legendre type shape fuctions. The variations with respect to the different path of domain integrals from the crack-tip front and the choice of 5-function have been tested by several examples.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.319-326
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed for the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted for in the sense of von Karman hypothesis. The material model Is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized for anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The Integrals of Legendre Polynomials we used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version finite element model is demonstrated through several comparative points of view in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic zone

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• 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.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.48 no.1
• /
• pp.61-68
• /
• 2006

The p-version nonlinear finite element model has been developed to analyze the nonlinear behavior of simply supported RC slabs strengthened with carbon fiber reinforced plastic sheets. The shape function is adopted with integral of Legendre polynomials. The compression model of concrete is based on the Kupfer's yield criterion, hardening rule, and crushing condition. The cracking behavior is modeled by a smeared crack model. In this study, the fixed crack approach is adopted as being geometrically fixed in direction once generated. Each steel layer has a uniaxial behavior resisting only the axial force in the bar direction. Identical behavior is assumed fur tension and compression of steel according to the elastic modulus. The carbon fiber reinforced plastic sheets are considered as reinforced layers of equivalent thickness with uniaxial strength and rigidity properties in the present model. It is shown that the proposed model is able to adequately predicte the displacement and ultimate load of nonlinear simply supported RC slabs by a patch with respect to reinforcement ratio, thickness and angles of CFRP sheets.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.