• East Asian mathematical journal
• /
• v.14 no.1
• /
• pp.63-76
• /
• 1998

we consider the hp-version to solve non-constant coefficients elliptic equations $-div(a{\nabla}u)=f$ with Dirichlet boundary conditions on a bounded polygonal domain $\Omega$ in $R^2$. In [6], M. Suri obtained an optimal error-estimate for the hp-version: ${\parallel}u-u^h_p{\parallel}_{1,\Omega}{\leq}Cp^{(\sigma-1)}h^{min(p,\sigma-1)}{\parallel}u{\parallel}_{\sigma,\Omega}$. This optimal result follows under the assumption that all integrations are performed exactly. In practice, the integrals are seldom computed exactly. The numerical quadrature rule scheme is needed to compute the integrals in the variational formulation of the discrete problem. In this paper we consider a family $G_p=\{I_m\}$ of numerical quadrature rules satisfying certain properties, which can be used for calculating the integrals. Under the numerical quadrature rules we will give the variational form of our non-constant coefficients elliptic problem and derive an error estimate of ${\parallel}u-\tilde{u}^h_p{\parallel}_{1,\Omega}$.

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

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

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

• Advances in aircraft and spacecraft science
• /
• v.5 no.1
• /
• pp.51-71
• /
• 2018

This work proposes a theoretical and numerical study on the behavior of a tapered shaft rotor made of composite materials by the classical version h and the version p of the finite element method. Hierarchical form functions are used to define the model. The purpose of this paper is to determine the expressions of the kinetic and potential energies of the tree necessary for the results of the equations of motion. A comparison between the version h and the p version of the finite element method of the functions of polynomial and trigonometric hierarchical forms with six degrees of freedom per node, of a composite tapered and cylindrical shaft which rotates at a constant speed about its axis. It is found that when the number of functions of form (the version p) is increased, the solution converges. It is also observed that the conicity of the shaft increases the rigidity with respect to a uniform shaft having the same mechanical properties. The numerical simulation allowed us to determine the natural frequencies and the critical speeds of the composite shaft systems are compared with those available in the literature and the effectiveness of the methods used are discussed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.4
• /
• pp.387-407
• /
• 2015

This paper analyzes the $h{\times}p$ version of the finite element method for optimal control problems constrained by elliptic partial differential equations with random inputs. The main result is that the $h{\times}p$ error bound for the control problems subject to stochastic partial differential equations leads to an exponential rate of convergence with respect to p as for the corresponding direct problems. Numerical examples are used to confirm the theoretical results.

• 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
• /
• 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
• /
• 1995.10a
• /
• pp.145-152
• /
• 1995

A general stiffener element which includes transverse shear deformation(TSD) is formulated using the p-version of finite element method. Hierarchic C"-shape functions, derived from Integrals of Legendre polynomials, are used to define the assembled stiffness matrix of the stiffener and plate on the basis of 5 D.0.F displacement fields. The stiffness matrix for 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 of the finite element method are compared with the results in literatures, especially those by the h-version software, MICROFEAP-II.P-II.

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