• Computational Structural Engineering
• /
• v.2 no.4
• /
• pp.15-20
• /
• 1989

이글에서는 2가지 예제를 통해 h-version과 p-version의 비교를 살펴보면서 p-version 해석이 h-version에 비해 상대적으로 많은 장점들을 가지고 있으며, 신뢰도, 정확도, 효율성, 경제성, 용장성 등 측면에서 우월함을 증명해 보였다. 특히 응력집중(stress concentration)이 일어나는 crack-tips, cut-outs, reentrant corners, presence of stiffners, mixed boundary conditions 등 많은 특이성(singularity) 문제에 더욱 적합함을 본 예제 외의 발표된 많은 논문들을 통해 알 수 있으며, 모델링의 단순성에 기인하여 사용이 매우 쉽다는 것도 무엇보다 큰 이점이라 하겠다. p-version은 h-version의 비효율성을 차수 p를 1, 2 또는 3으로 줄인 후 이 값을 고정시키고 다시 요소분할을 통해 진해(true solution)에 접근시키는 방식을 위하면 다시 종래의 h-version으로 환원되는 호환성을 갖고 있다는 것이다. 고로 구조해석에서 h-p version이 가장 이상적인 유한요소해석 방법이라 할 수 있겠는데, 다시 말하면 균열문제의 경우 균열선단(crack-tip)에서는 p-level을 높이고 (p=8, 9 or 10) 비교적 응력집중이 낮은 영역에서는 p-level을 낮춤으로써 (p=3, 4 or 5) 그 효율성을 극대화할 수 있겠다.

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

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.57-66
• /
• 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 high accuracy and the rapid 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 that the approximations enter 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 to tension.

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

• Computational Structural Engineering
• /
• v.5 no.1
• /
• pp.75-82
• /
• 1992

In the h-version of F.E.M., all piecewisely smooth curved boundaries can be approximated by a sufficient number of straight-sided elements. However, in the p-version the size of the element is usually large and hence the probability of distortions is more. An attempt has been made to generate a curved boundary by using a transfinite interpolation technique to avoid the discretization errors. In the following sections, it will be shown how to construct transfinite interpolants both in h-version and in p-version over polygonal and nonpolygonal regions. Three numerical tests are shown to validate the applicability and superior capability of transfinite interpolation technique.

• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.101-111
• /
• 1990

The elastic analysis of floor slabs using the p-version of finite element method encounters stress singularities at certain types of reentrant corners, openings and cut-outs. Results obtained using the computer code based on C.deg. - hierarchic plate element formulated by Reissner-Mindlin theory are compared with theoretical predictions and with computational results reported in the literature. The convergence rate of h-, p- and hp-version can be estimated on the basis of the energy norm in global sense. If accuracy in terms of the number of degree-of-freedom is used as a criterion, the solutions presented here are the most efficient that have been published up to date. Examples are the rhombic plate with the obtuse angle of 150.deg. and the square plate with cut-outs.

• Korean Journal of Mathematics
• /
• v.6 no.1
• /
• pp.17-26
• /
• 1998

In this paper, we briefly study the condition number of stiffness matrix with $h$-version and analyze it with $p$-version of the finite element method.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.217-228
• /
• 1997

The high precision analysis by the p-version of the finite element method are fairly well established as highly efficient method for linear elastic problems, especially in the presence of stress singularity. It has been noted that the merits of the p-version are accuracy, modeling simplicity, robustness, and savings in user's and CPU time. However, little has been done to exploit their benefits in elasto-plastic analysis. In this paper, the p-version finite element model is proposed for the materially nonlinear analysis that is based on the incremental theory of plasticity using the constitutive equation for work-hardening materials, and the associated flow rule. To obtain the solution of nonlinear equation, the Newton-Raphson method and initial stiffness method, etc are used. Several numerical examples are tested with the help of the square plates with cutout, the thick-walled cylinder under internal pressure, and the circular plate with uniformly distributed load. Those results are compared with the theoretical solutions and the numerical solutions of ADINA

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.4
• /
• pp.537-549
• /
• 1999

선형탄성 파괴해석은 균열을 갖는 변형률 경화재료의 파괴거동을 예측하는데 불충분하기 때문에 최근에는 균열 선단 부에서 대규모 소성 역을 갖는 균열 체에 적용할 수 있는 많은 파괴역학개념이 제안되고 있다. 따라서, 본 연구에서는 대규모항복 조건하의 연성파괴를 보이는 평판을 정확하게 해석할 수 있는 새로운 유한요소모델을 제시하고자 한다. 균열 선단 부의 응력 장을 정의하는데 가장 지배적인 파괴매개변수인 J-적분 값과 소성 역의 크기 및 형상을 J-적분법과 등가영역적분법을 통해 파괴거동을 설명할 수 있도록 증분소성이론에 기초를 둔 p-version 유한요소해석이 채택되었다. 제안된 유한요소모델에 의한 수치해석결과는 이론 해와 h-version 유한요소해석과 비교되었다.