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

• The Journal of the Acoustical Society of Korea
• /
• v.34 no.5
• /
• pp.335-342
• /
• 2015

Phononic crystals are composite materials consisting of a periodic arrangement of scattering inclusions in a host material. One of the most important properties of phononic crystals is the existence of band gaps, i.e., ranges of frequencies at which acoustic waves cannot propagate through the structure. The present study aims to investigate theoretically and experimentally the acoustic band structures in two-dimensional (2D) phononic crystals consisting of periodic square arrays of stainless steel solid cylinders with a diameter of 1 mm and a lattice constant of 1.5 mm in water. The theoretical dispersion relation that depicts the relationship between the frequency and the wave vector was calculated along the ${\Gamma}X$ direction of the first Brillouin zone using the finite element method to predict the band structures in the 2D phononic crystals. The transmission and the reflection coefficients were measured in the 2D phononic crystals with 1, 3, 5, 7, and 9 layers of stainless steel cylinders stacked in the perpendicular direction to propagation at normal incidence. The theoretical dispersion relation exhibited five band gaps at frequencies below 2 MHz, the first gap appearing around a frequency of 0.5 MHz. The location and the width of the band gaps experimentally observed in the transmission and the reflection coefficients appeared to coincide well with those determined from the theoretical dispersion relation.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.1
• /
• pp.99-112
• /
• 1999

By using 92 values of lateral earth pressure coefficient(K) measured in Korea, the tendency of K with varying depth is analyzed and compared with the range of K defined by Hoek and Brown. The horizontal stress is generally larger than the vertical stress in Korea : About 84 % of K values are above 1. In this study, the theory of elasto-plasticity is applied to analyze the variation of K values, and the results are compared with those of numerical analysis. This reveals that the erosion, sedimentation and weathering of earth crust are important factors in the determination of K values. Surface erosion, large lateral pressure and good rock mass increase the K values, but sedimentation decreases the K values. This study enable us to analyze the effects of geological processes on the K values, especially at shallow depth where underground excavation takes place. A neural network expert system using multi-layer back-propagation algorithm is developed to predict the K values. The neural network model has a correlation coefficient above 0.996 when it is compared with measured data. The comparison with 9 measured data which are not included in the back-propagation learning has shown an average inference error of 20% and the correlation coefficient above 0.95. The expert system developed in this study can be used for reliable determination of K values.

• Composites Research
• /
• v.27 no.4
• /
• pp.158-167
• /
• 2014

Final goal of this research is to establish the database for correlation factors which connects the test and analysis results of shear buckling allowables for composite plate. To accomplish the goal, extensive test and analysis works are required. In this paper, as the first step, a frame-type fixture for shear buckling test was designed and validated through the test and analysis. Final configuration of the fixture were determined via parametric study on the effect of specimen size, cross-sectional dimensions, and number of fastening bolts on the shear buckling load. Results of the study showed the designed frame-type fixture successfully induces the shear buckling of composite plate. However, there were deviations between the test results and analysis results for ideal case under pure shear load, which were mainly caused by the difference in plate sizes for both cases. The difference were larger in the plates with larger hole and simply supported boundary condition. It is concluded from the results that while the designed fixture can be used for the clamped plates with acceptable accuracy, it shows larger difference in the simply supported plates.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.9
• /
• pp.3176-3183
• /
• 2010

A multi-scale fatigue life prediction methodology of composite pressure vessels subjected to multi-axial loading has been proposed in this paper. The multi-scale approach starts from the constituents, fiber, matrix and interface, leading to predict behavior of ply, laminates and eventually the composite structures. The multi-scale fatigue life prediction methodology is composed of two steps: macro stress analysis and micro mechanics of failure based on fatigue analysis. In the macro stress analysis, multi-axial fatigue loading acting at laminate is determined from finite element analysis of composite pressure vessel, and ply stresses are computed using a classical laminate theory. The micro stresses are calculated in each constituent from ply stresses using a micromechanical model. Three methods are employed in predicting fatigue life of each constituent, i.e. a maximum stress method for fiber, an equivalent stress method for multi-axially loaded matrix, and a critical plane method for the interface. A modified Goodman diagram is used to take into account the generic mean stresses. Damages from each loading cycle are accumulated using Miner's rule. Monte Carlo simulation has been performed to predict the overall fatigue life of a composite pressure vessel considering statistical distribution of material properties of each constituent, fiber volume fraction and manufacturing winding angle.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.7-17
• /
• 2011

We have developed an inversion code for three-dimensional (3D) resistivity tomography including the anisotropy effect. The algorithm is based on the finite element approximations for the forward modelling and Active Constraint Balancing method is adopted to enhance the resolving power of the smoothness constraint least-squares inversion. Using numerical experiments, we have shown that anisotropic inversion is viable to get an accurate image of the subsurface when the subsurface shows strong electrical anisotropy. Moreover, anisotropy can be used as additional information in the interpretation of subsurface. This algorithm was also applied to the field dataset acquired in the abandoned old mine area, where a high-rise apartment block has been built up over a mining tunnel. The main purpose of the investigation was to evaluate the safety analysis of the building due to old mining activities. Strong electrical anisotropy has been observed and it was proven to be caused by geological setting of the site. To handle the anisotropy problem, field data were inverted by a 3D anisotropic tomography algorithm and we could obtain 3D subsurface images, which matches well with geology mapping observations. The inversion results have been used to provide the subsurface model for the safety analysis in rock engineering and we could assure the residents that the apartment has no problem in its safety after the completion of investigation works.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.15 no.3
• /
• pp.271-287
• /
• 2013

In this study, more convenient model(S-model) for predicting ground surface settlement is developed through comparing field monitoring data of the domestic subway applied shield TBM method with conventional equation & numerical analysis. Sample stations are chosen from whole of excavation section and lateral & vertical ground surface settlement characteristic with excavation are analysed. Based on analysis result, through the comparison with actual monitoring data, the model that is possible to compute maximum surface settlement and settlement influence area is suggested with assumption that lateral surface settlement forms are composed relaxed zone and elastic zone. In addition, vertical ground surface settlement patterns with excavation are similar to cubic-function and S-model with assumption that coefficients are function of tunnel diameter and depth is suggested. Consequently, the ground surface settlement patterns are significantly similar to actual monitoring data and numerical method result. Thus, as a result, when tunnels are excavated using sheild TBM through rather soft weathered soil & rock layer, prediction of ground surface settlement with excavation using convenient S-model is practicable.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.41 no.6
• /
• pp.25-31
• /
• 2004

In this paper, we introduce kinetic Monte Carlo (kMC) methods for simulating diffusion process in nano-scale device fabrication. At first, we review kMC theory and backgrounds and give a simple point defect diffusion process modeling in thermal annealing after ion (electron) implantation into Si crystalline substrate to help understand kinetic Monte Carlo methods. kMC is a kind of Monte Carlo but can simulate time evolution of diffusion process through Poisson probabilistic process. In kMC diffusion process, instead of. solving differential reaction-diffusion equations via conventional finite difference or element methods, it is based on a series of chemical reaction (between atoms and/or defects) or diffusion events according to event rates of all possible events. Every event has its own event rate and time evolution of semiconductor diffusion process is directly simulated. Those event rates can be derived either directly from molecular dynamics (MD) or first-principles (ab-initio) calculations, or from experimental data.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.5
• /
• pp.483-490
• /
• 2016

In plant industry, sulfur storage tank is made of steel and annular plate is connected with concrete foundation of ring wall type by anchor bolt. Due to keep sulfur at high temperature in tank by coil, sulfur storage tank is expanded larger than another tank stores fluid at room temperature. Generally, structural design of tank foundation is performed analysis with loading of temperature gradient between inner and outer surface, this method can't consider the phenomenon that load is intensively transferred to concrete foundation at anchor bolt. This means that temperature load is underestimated and causes crack of concrete near anchor bolt. In this study, evaluation formula considering temperature load transfer mechanism through anchor bolt is proposed and load acting on concrete foundation is rationally decided. For this purpose, it is analyzed variation of thermal load per various anchor bolt number using finite element model including tank annular plate and anchor bolt. Solution is proposed as specified term combining result of analysis and theoretical solution for evaluating load transferred by anchor bolt. For confirmation of validation of proposed formula, it is applied in design of sulfur storage tank at plant site, it shows that the formula can be practically applied.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.59-64
• /
• 1992

터빈 블레이드와 같이 회전하는 구조물의 파단은 공진 근처에서 진동이 발 생할 때에 이에 기인하는 피로에 의하여 발생한다. 그러므로 이와 같은 파단 을 피하기 위해서는 설계 단계에서 이론적인 계산에 의하여 구조물의 고유 진동수를 결정하는 것이 상당히 중요하다. 판이 회전을 받게 되면 원심력에 의하여 판의 강성이 증가하므로 고유진동수가 회전하지 않는 판의 고유진동 수보다는 상당히 증가하게 된다. 이에 대한 연구가 국내외에서 상당수 행하 여졌지만, 연구의 대부분이 회전의 영향을 고려하지 않은 정지판(stationary plate)에 대한 것이며 뢰전을 고려한 연구는 극히 제한되어 있다. 또한 회전 의 영향을 고려한 연구의 대부분이 해석 대상을 보로서 단순화 시켰고 해법 으로는 유한요소법과 Ritz법 등을 사용하였다. 이는 블레이드가 지니고 있는 기하학적인 형상과 진동 특성이 해석적인 방법으로 해결하는 데에는 상당한 어려움이 있기 때문이다. 실제적으로는 터빈 블레이드와 같은 회전체의 진동 특성이 설치각이나 비틀림각, 판의 형상비, 회전속도 등의 변화에 의하여 영 향을 받기 때문에 보와 같은 진동 거동을 보이기보다는 판이나 셀과 같은 진동 거동을 보이므로 보다 정확한 해석을 수행하기 위해서는 해석 대상을 판이나 셀로서 취급하는 것이 타당하다. 따라서 본 연구에서는 위와 같은 이 유 때문에 해석 대상을 등방성 사각판과 직교이방성 복합재료 사각판으로 선택하였으며, 구조물의 고유진동수에 영향을 미치는 다음과 같은 인자들을 해석에 고려하였다. 1. 회전속도 (rotational speed) 2. 설치각 (setting angle) 3. 허브의 반경 (hub radius) 4. 판의 형상비 (aspect ratio) 5. 적층순서 (stacking sequence)구조물에 대한 동적실험(dynamic test)을 통하여 단기간에 동적특성을 결정하고 SDM(structure dynamic modification)이나 FRS(force response simulation)를 수행하여 임의의 좌표 공간에 대한 진동수준을 해석적으로 예측할 뿐만 아니라 구조물의 진동제어 를 위한 동적인자를 변경시킬 수 있는 정보를 제공하며 장비를 방진할 경우 신뢰성 있는 전달률을 결정할 수 있다. 실험적으로 철교, 교량이나 건물의 철골구조 및 2층 바닥 등 대,중형의 복잡한 구조물에 대항 동특성을 나타내 는 모빌리티를 결정할 경우 충격 가진 실험이 사용되는 실험장비 측면에서 나 실험을 수행하는 과정이 대체적으로 간편하다. 그러나 이 경우 대상 구조 물을 충분히 가진시킬수 있는 용량의 대형 충격기(large impact hammer)가 필요하게 된다. 이러한 동적실험은 약 길이 61m, 폭 16m의 4경간 교량에 대 하여 동적실험을 수행하여 가능성을 확인하였다. 여기서는 실험실 수준의 평 판모델을 제작하고 실제 현장에서 이루어질 수 있는 진동제어 구조물에 대 한 동적실험 및 FRS를 수행하는 과정과 동일하게 따름으로써 실제 발생할 수 있는 오차나 error를 실험실내의 차원에서 파악하여 진동원을 있는 구조 물에 대한 진동제어기술을 보유하고자 한다. 이용한 해마의 부피측정은 해마경화증 환자의 진단에 있어 육안적인 MR 진단이 어려운 제한된 경우에만 실제적 도움을 줄 수 있는 보조적인 방법으로 생각된다.ofile whereas relaxivity at high field is not affected by τS. On the other hand, the change in τV does not affect low field profile but stron