• Journal of the Earthquake Engineering Society of Korea
• /
• v.4 no.2
• /
• pp.83-98
• /
• 2000

자오방향 및 주변방향으로 피르스트레스트 하중이 작용된 축대칭 쉘 구조물을 기하학적으로 축대칭인 구조물의 특성을 최대한으로 이용할 수 있도록 회전 링요소로 모델화하였다 보강링 요소의 모델은 축대칭 쉘요소를 이용하였으며 본체 구조물과 절점에서 부착되있는 것으로 가정하여 이의 편심을 고려하였다 유체-구조물의 상호관계는 접촉면에서 구조물의 가속도에 비례한 부가질량으로 표현하였으며 부가질량은 유체를 비점성 비압축 및 비회전을 가정하여 유한요소법에 의해 구하였다 이에 대한 수치해석을 통하여 고유진동해석 및 지진하중을 주하중으로 한 동적해석을 실시하였다 프로그램을 통하여 해석한 결과를 프리스트레스 하중 하에서 고유진동수에 대한 정해와 비교한 결과 20개의 요소로 모델링한 경우에서도 정해와 근접한 해를 얻을 수 있었다 또한 내부유체가 있는 경우와 링보강을 한 경우에 대한 고유진동수를 문헌과 비교한 결과 근접한 해를 얻을 수 있었다.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.12
• /
• pp.27-34
• /
• 2020

Use of axisymmetric shell element for the structure increases the efficiency and accuracy in finite element analysis of the interaction between the ground and the structure. This paper derived the force balance equation and the moment balance equation for an axisymmetric shell element based on Kirchhoff's theory. The governing equation for the axial deformation used the isoparametric shape function in the Galerkin formulation, and the governing equation for the shell bending used the higher-order shape function. The developed axisymmetric shell element was combined with Geo-COUS, a geotechnical finite element program for the coupled analysis with the ground. The accuracy of the developed element was confirmed through the example analyses of the circular plate and the liquid storage tank. And the energy balance equation for the axisymmetric shell element is presented.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.165-172
• /
• 2003

In this study, a new and efficient harmonic axisymmetric shell element for static and dynamic analysis Is proposed. The present element considering shear strain is based on a modified mixed variational principle in which the independent unknowns are only the Quantities prescribable at the shell edges. Unlike existing hybrid-mixed axisymmetric shell elements, the present element introduces additional nodeless degrees for displacement field Interpolation In order to enhance the numerical performance. The stress parameters are eliminated by the stationary condition and the nodeless degrees are condensed out by the Guyan reduction. Through several numerical examples, the hybrid-miked shell element with the additional nodeless degrees and the consistent stress parameters is shown to be efficient and yield very accurate results for static and vibration analysis.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1994.10a
• /
• pp.295-299
• /
• 1994

두꺼운 원통형 쉘은 공학적인 문제에서 많이 사용된다. 쉘 내부에 임피던스가 큰 유체와 구조물이 있을 때 쉘을 포함한 진동해석은 이론적인 해석이 매우 어렵다. 쉘 내부에 있는 유체의 임피던스가 공기에 비하여 매우 클 경우 쉘과 유체, 내부의 구조물과 유체사이의 구조물-유체 상호작용(structure-fluid interaction)이 고려되어야 한다. 얇은 원통형 쉘에 대해서는 상용 유한요소 코드를 이용하여 구조물-유체 상호작용을 고려한 진동해석이 많이 수행되었으나 축대칭 두꺼운 원통형 쉘에 대해서는 연구가 수행되지 않고 있다. 본 연구에서는 NASTRAN, ANSYS 같은 상용 유한요소 코드에서 지원되지 않는 축대칭 두꺼운 원통형 쉘 내부에 유체와 강체요소가 있을 경우 이에 대한 유한요소 코드를 개발하고, 구조물-유체 상호작용을 고려하여 진동해석을 하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.179-187
• /
• 2000

When one considers the property of the axisymmetry, an analysis of an axisymmetric shell subjected to unaxisymmetric loading can be employed to save time and computer memory space. If one considers the Fourier series of the circumference direction of loads and displacements, an axisymmetric tank subjected to a nonaxisymmetric load can be treated as a frame element. Using the Fourier series, the authors derived the stiffness matrix of the cylindrical shell subjected to unaxisymmetric loading by the usual finite element method, and converted the stiffness matrix of a frame element into a transfer matrix by rearranging the stiffness matrix to apply the transfer matrix method. Here the most significant purpose of this paper is to achieve the fewest number of simultaneous equations for analysing an axisymmetric shell subjected to a nonaxisymmetric load. The results of the proposed method of the analysis of the cylindrical shell subjected to a wind load and a water load show no differences when compared to the other methods.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.2
• /
• pp.193-201
• /
• 2004

The aim of this study is focused on getting an almost exact solution which is the simplicity and exactness of an axisymmetrically loaded cylindrical shell. This method replaces the finite element method which is a very powerful tool for analysis of any kind of structure which has an arbitrary shape, but is still a numerical analysis. Instead, this study uses the method of distribution of end actions which is a kind of iteration technique to implement the leading matrix method. The distribution and carry-over factors of a cylinder are calculated by the theory of a differential equation of a beam on an elastic foundation. The results are satisfactory when this method is applied to a cylinder that is subjected to a concentrated load and hydrostatic pressure when compared with the BEF analogy separately.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.9 no.1
• /
• pp.271-281
• /
• 2005

The effects of radiation damping is used to compensate the truncated boundary which is relatively close to the structure-fluid interface in the fluid element surrounding the submerged structures. An efficient ring element is presented to model the shell and fluid element which fully utilizes the characteristics of the axisymmetry. The computational model uses the technique which separate the meridional shape and circumferential wave mode and gets similar result with the exact solution in the eigenvalues and the earthquake analysis. The fluid-structure interaction techniques is developed in the finite element analysis of two dimensional problems using the relations between pressure, nodal unknown acceleration and added mass assuming the fluid to be invicid, incompressible and irrotational. The effectiveness and efficiency of the technique is demonstrated by analyzing the free vibration and seismic analysis using the added mass matrix considering the structural deformation effect.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.4
• /
• pp.71-82
• /
• 1999

An axisymmetric shell element which includes the effects of the meridional and circumferential cable prestresses is developed. The fluid-structure interaction is expressed as added mass effect which is in proportion to the acceleration of the structure in interface surface. The added mass is obtained by using finite element method under the assumption that the fluid is invicid, incompressible and irrotational. It is coded for personal computer by the maximum use of axisymmetic properties and the dynamic analysis are performed under seismic exitations. A ring element makes the characteristics of the axisymmetric shell to be fully utilized. The elgenvalue solutons under the initial prestresses and the internal fluid are well agreed with the exact solutions and references by using under 20 elements. The eigenvalues are decreased along the increasing the height of internal fluid and these effects are dominant under the lower wave numbers. The results of the seismic analysis show that the radial deflection under the meridional prestress is a little larger than that under the circumferential prestress.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.197-206
• /
• 2003

A three-dimensional (3-D) method of analysis is presented for determining the free vibration frequencies and mode shapes of solid and hollow hemispherical shells of revolution of arbitrary wall thickness having arbitrary constraints on their boundaries. Unlike conventional shell theories, which are mathematically two-dimensional (2-D), the present method is based upon the 3-D dynamic equations of elasticity. Displacement components μ/sub Φ/, μ/sub z/, and μ/sub θ/ in the meridional, normal, and circumferential directions, respectively, are taken to be sinusoidal in time, periodic in θ, and algebraic polynomials in the Φ and z directions. Potential (strain) and kinetic energies of the hemispherical shells are formulated, and the Ritz method is used to solve the eigenvalue problem, thus yielding upper bound values of the frequencies obtained by minimizing the frequencies. As the degree of the polynomials is increased, frequencies converge to the exact values. Novel numerical results are presented for solid and hollow hemispheres with linear thickness variation. The effect on frequencies of a small axial conical hole is also discussed. Comparisons are made for the frequencies of completely free, thick hemispherical shells with uniform thickness from the present 3-D Ritz solutions and other 3-D finite element ones.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.8 no.1
• /
• pp.195-202
• /
• 2004

The shell due to the effect of initial normal pressures on the shell surface was based on the assumption that the directions of the pressures are always normal to the undeformed shell surface, and that the change in the surface area of the shell is negligible. But the fact that the pressure are always normal to the deforming surface leads "follower force". The follower effect in the analysis can significantly alter the solution for natural frequency and buckling load as compared to the case when the direction of the pressures are assumed to be normal to the uniform shell surface. The expression for the part of strain energy contribution from normal pressure due to the effect of follower force was derived and added to the element stiffness matrix of axisymmetric shell. In the case of increasing external pressure, the natural frequencies decrease until one of them reaches zero. Theoretically the smallest applied load that reduces the frequency of any mode to zero, will have same magnitude as that of the buckling load. In order to determine the bucking load of the shell a few sets of frequencies are computed and the results considering the follower effects are well with the exact solution while the case without that are quite different. But in case of hemispherical dome, there are little difference in buckling pressure between with and without the effect of follower force.