• Structural Engineering and Mechanics
• /
• v.23 no.5
• /
• pp.487-504
• /
• 2006

Several new versatile two-dimensional p-version finite elements are developed. The element matrices are integrated analytically to guarantee the accuracy and monotonic convergence of the predicted solutions of the proposed p-version elements. The analysis results show that the convergence rate of the present elements is very fast with respect to the number of additional Fourier or polynomial terms in shape functions, and their solutions are much more accurate than those of the linear finite elements for the same number of degrees of freedom. Additionally, the new p-version plate elements without the reduced integration can overcome the shear locking problem over the conventional h-version elements. Using the proposed p-version elements with fast convergent characteristic, the elasto-plastic impact of the structure attached with the absorber is simulated. Good agreement between the simulated and experimental results verifies the present p-version finite elements for the analyses of structural dynamic responses and the structural elasto-plastic impact. Further, using the elasto-plastic impact model and the p-version finite element method, the absorber of the T structure on the Qiantang River is designed for its collision protection.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1992.10a
• /
• pp.39-44
• /
• 1992

This paper is concerned with formulations of the hierarchical $C^{o}$-plate element on the basis of Reissner-Mindlin plate theory. On reason for the development of the aforementioned element 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 we proposed to verify the superior convergence and algorithmic efficiency with the help of the clamped L-shaped plate problems.s.

• Journal of the Korean GEO-environmental Society
• /
• v.18 no.2
• /
• pp.53-58
• /
• 2017

In this paper, we propose a new method to debond between mixed finite elements for porous media in ABAQUS (2014). ABAQUS just provides debonding algorithm for the u-p model using cohesive elements in standard version. However, this approach has a drawback that it is hard to simulate complex debonding problems like element separation, rigid body motion, and contact between separated elements in standard version. ABAQUS-explicit can resolve these complex problems, but cohesive elements for the u-p model cannot be applied. We introduce a new algorithm for debonding for porous media instead of using cohesive elements. In this method, subroutines VUMAT to apply constitutive models and VDISP to separate elements in ABAQUS are used to simulate debonding problems. In addition, a simple 2-D example is demonstrated in the ABAQUS-explicit solver.

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

The p-version crack model based on integrals of Legendre polynomial and virtual crack extension method is proposed with its potential for application to stress intensity factor computations in linear elastic fracture mechanics. The main advantage of this model is that the data preparation effort is minimal because only a small number of elements are used and the high accuracy and the rapid rate of convergence can be achieved in the vicinity of crack tip. There are two important findings from this study. Firstly, the limit value, the strain energy of the exact solution can be estimated with successive three p-version approximations by ascertaining the approximations is entered the asymptotic range. Secondly, the rate of convergence of p-version model is almost twice that of h-version model on the basis of uniform or quasiuniform mesh refinement for the cracked panel problem subjected tension.

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

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

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.345-350
• /
• 2000

Adaptive finite element analysis, which its solution error meets with the user defined allowable error, is recently used far improving reliability of finite element analysis results. This adaptive analysis is composed of two procedures; one is the error estimation of an analysis result and another is the reconstruction of finite elements. In the rp-method, an element size is controlled by relocating of nodal positions(r-method) and the order of an element shape function is determined by the hierarchical polynomial(p-method) corresponding to the element solution error. In order to show the effectiveness and accuracy of the suggested rp-method, various numerical examples were analyzed and these analysis results were examined by comparing with those obtained by the existed methods. As a result of this study, following conclusions are obtained. (1) rp-method is more accurate and effective than the r- and p-method. (2) The solution convergency of the rp-method is controlled by means of the iterative calculation numbers of the r- and p- method each other.

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