• Structural Engineering and Mechanics
• /
• v.2 no.3
• /
• pp.245-256
• /
• 1994

A new axisymmetric crack model is proposed on the basis of p-version of the finite element method limited to theory of small scale yielding. To this end, axisymmetric stress element is formulated by integrals of Legendre polynomial which has hierarchical nature and orthogonality relationship. The virtual crack extension method has been adopted to calculate the stress intensity factors for 3-D axisymmetric cracked bodies where the potential energy change as a function of position along the crack front is calculated. The sensitivity with respect to the aspect ratio and Poisson locking has been tested to ascertain the robustness of p-version axisymmetric element. Also, the limit value that is an exact solution obtained by FEM when degree of freedom is infinite can be estimated using the extrapolation equation based on error prediction in energy norm. Numerical examples of thick-walled cylinder, axisymmetric crack in a round bar and internal part-thorough cracked pipes are tested with high precision.

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

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

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