• 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.12 no.4
• /
• pp.537-549
• /
• 1999

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

• 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.4 no.2
• /
• pp.77-85
• /
• 1991

A new finite element technology based on the p-version of E.F.M. is discussed with reference to its potential for application to stress intensity factor computations in linear elastic fracture mechanics, especially cracked cylindrical shells. It is shown that the p-version model is far better suited for computing the stress intensity factors than the conventional h-version models with the help of three test problems. The main advantage of this technology is that the accuracy of approximation can be established without mesh refinement or the use of special procedures such as crack-tip element and mixed variational approach.

• Computational Structural Engineering
• /
• v.5 no.3
• /
• pp.97-103
• /
• 1992

Since many design problems in the railroad, aerospace and machine structures involve considerations of the effect of cyclic loading, manufacturing and quality control processes much fully account for fatigue of critical components. Due to the sensitivity of the Paris law, it is very important to calculate .DELTA.K numerically to minimize the error of predicted fatigue life in cycles. However, it is shown that the p-version of FEM based on LEFM analysis is far better suited for computing the stress intensity factors than the conventional h-version. To demonstrate the proficiency of the proposed scheme, the welded T-joint with crack problems of box car body bolster assembly and a crack problem emanating from a circular hole in finite strip have been solved.

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

• 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) 그 효율성을 극대화할 수 있겠다.

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

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