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

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

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

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

This paper presents a p-version finite element approach for modeling the stress distribution around a circular hole in a finite strip subjected to membrane and flexural behaviors. Also, same problem with a crack emanating from a perforated tension strip was solved by virtual crack extension method. The p-version of the finite element method based on integrals of Legendre polynomials is shown to perform very well for modeling geometries with very steep stress gradients in the vicinity of a circular cutout. Here, the transfinite mapping technique for circular boundaries was used to avoid the discretization errors. The numerical results from the proposed scheme have a good comparison with those by Nisida, Howland, Newman etc. and the conventional finite element approach.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.33-45
• /
• 1996

A general stiffener element which includes transverse shear deformation is formulated using the p-version finite element method. Hierarchic C/sup o/-shape functions, derived from Integrals of Legendre polynomials, are used to define the assembled stiffness matrix of the stiffener with respect to the local reference frame is transformed to the plate reference system by applying the appropriate transformation matrices in order to insure compatibility of displacements at the junction of the stiffener and plate. The transformation matrices which account for the orientation and the eccentricity effects of the stiffener with respect to the plate reference axes are used to find local behavior at the junction of the stiffener and the relative contributions of the plate and stiffener to the strength of the composite system. The results obtained by the p-version finite element method are comared with the results in literatures, especially those by the h-version finite element analysis program, MICROFEAP-II.