• Journal of Korean Society of Steel Construction
• /
• v.12 no.3 s.46
• /
• pp.303-315
• /
• 2000

The main object of this study is to analyze mechanical behaviors as anisotropic three-dimensional body under various static loads. This paper presents the applicability of the finite difference method to three dimensional problem of anisotropic body. The finite difference method as applied here is generalized to anisotropic three-dimensional problem of elastic body where the governing differential equations of equilibrium of such bodies are expressed in terms of the displacement u, v, and w in the coordinates axes x, y and z, care being taken to modify the finite difference expressions to satisfy the appropriate boundary conditions. By adopting a new three dimensional finite difference modelling including elimination of pivotal difference points in the case of free boundary condition, the three dimensional problem of anisotropic body was successfully completed. Several numerical results show quick convergence and numerical validity of finite difference technique in three dimensional problem.

• Tunnel and Underground Space
• /
• v.18 no.3
• /
• pp.202-213
• /
• 2008

In this study, theoretical approach of 3D seismic traveltime tomography was investigated. To guarantee the successful field application of 3D tomography, appropriate control of problem associated with blind zone is pre-requisite. To overcome the velocity distortion of the reconstructed tomogram due to insufficient source-receiver array coverage, the algorithm of 3D seismic traveltime tomography based on the Fresnel volume was developed as a technique of ray-path broadening. For the successful reconstruction of velocity cube, 3D traveltime algorithm was explored and employed on the basis of 2nd order Fast Marching Method(FMM), resulting in improvement of precision and accuracy. To prove the validity and field application of this algorithm, two numerical experiments were performed for globular and layered models. The algorithm was also found to be successfully applicable to field data.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.36 no.6
• /
• pp.474-482
• /
• 2016

A Kelvin foam plate - widely used in the energy and transport industries as a lightweight structural material - was examined to estimate its Young's modulus using ultrasound. An isotropic tetrakaidecahedron foam structure was designed in SolidWorks and printed using 3D printer with an ABS plastic material. The 3D printed foam structure was used to build a foam plate with a 14 mm thickness ($50mm{\times}100mm$ in size) for the ultrasonic test. The Kelvin foam plate, a significantly porous medium, was completely filled with paraffin wax to enable the ultrasound to penetrate through the porous medium. The acoustic wave velocity of the wax-filled Kelvin foam was measured using the time of flight (TOF) method. Furthermore, the elastic modulus of the Kelvin foam was estimated based on an elastic structural model developed in this study. The Young's modulus of the produced Kelvin foam was observed to be approximately 3.4% of the bulk value of the constituent material (ABS plastic). This finding is consistent with experimental and theoretical results reported by previous studies.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.192-195
• /
• 2010

본 논문에서는 해상 크레인과 중량물의 동적 거동을 시뮬레이션하기 위해, 유한 요소 정식화(finite element formulation)를 이용하여 해상 크레인의 붐(boom)을 탄성체로 모델링 하였다. 붐은 3차원 탄성 빔(beam) 요소로 가정하고, 각 요소의 변형에 의한 변위는 형상 함수(shape function)과 절점 좌표(nodal coordinate)를 이용하여 정의하였다. 변형 변위를 이용하여 탄성 붐의 강성 행렬(stiffnes matrix)을 유도하고, 탄성 변위를 포함하는 위치 벡터를 이용하여 질량 행렬을 유도한다. 해상 크레인과 중량물로 이루어진 운동 방정식에 탄성 붐을 포함하여 유연 다물체계(flexible multibody system) 운동 방정식을 구성한다. 외력으로는 선박 유체정역학적 힘, 유체동역학적 힘, wire rope의 장력, 중력 그리고 계류력(mooring force)이 고려되었다. 먼저 요소의 개수를 변경하며 탄성 붐의 동적 거동을 시뮬레이션 하여, 유한 요소 정식화를 이용한 모델링의 타당성을 검증하였다. 그리고 해상 크레인과 중량물의 동적 거동 시뮬레이션에 탄성 붐 모델을 적용하였다.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.6 no.3
• /
• pp.247-258
• /
• 2004

Two 3-dimensional data processing techniques to predict the fractured zone ahead of a tunnel face by the tunnel seismic survey were proposed so that the geometric formation of the fractured zone could be estimated. The first 3-dimensional data processing technique was developed based on the principle of ellipsoid, The input data needed for the 3D migration can be obtained from the 2-dimensional tunnel seismic prediction (TSP) test where the TSP test should be performed in each sidewall of a tunnel. The second 3-dimensional migration technique that was developed based on the concept of wave travel plane was proposed. This technique can be applied when the TSP is operated with sources in one sidewall of a tunnel while the receivers are installed in both sidewalls. New migration technique was applied to an in-situ tunnelling site. The 3-dimensional migration was performed using measured TSP data and its results were compared with the geological investigation results that were monitored during tunnel construction. This comparison revealed that the proposed migration technique could reconstruct the discontinuity planes reasonably well.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.1
• /
• pp.39-48
• /
• 1994

An equivalent homogeneous 3-D linear composite analysis and accomponying finite element program is presented for elastomeric bearings. This study is limited to the 3-D layered system with linear, elastic, isoparametric small deformation. And we used method of multiscale to model the 3-dimensional configurations and overall response of the layered elastomeric bearings with global and local coordinates. The primary dependent variables for the theory have been selected that require only $C_o$ continuity of the finite element analysis. As a result, it is very simple and computationally economical. The presented theory can also be applied easily to the analysis of nonlinear behavior of layered systems. And those of past are not applicable to nonlinear analysis, because it uses superposition theory. Numerical examples are presented to verify the theory and to illustrate potential applications of the analysis.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.4 no.4
• /
• pp.301-317
• /
• 2002

Recently, a geophysical exploration technology is frequently utilized in the civil engineering field as well as in the resource exploration. It might be important for civil engineers to understand the fundamental theory of seismic survey and limitation of the technique when utilizing these techniques in the civil engineering field. A 3-dimensional migration technique based on the principle of ellipsoid to predict the fractured zone ahead of tunnel face utilizing the tunnel seismic survey was proposed so that the geometry of the fractured zone can be estimated, i.e. the angle between tunnel axis and discontinuity zone, and the dip. Moreover, a numerical analysis technique to simulate the TSP (Tunnel Seismic Prediction) test was proposed in this paper. Based on parametric studies, the best element size, the analysis time step, and the dynamic characteristics of pressure source were suggested to guarantee the stability and accuracy of numerical solution. Example problems on a hypothetical site showed the possibility that the 3-dimensional migration technique proposed in this paper appropriately estimate the 3D-geometry of fractures ahead of tunnel face.

• Computational Structural Engineering
• /
• v.4 no.1
• /
• pp.20-27
• /
• 1991

본 논문에서는 유체-구조 상호작용해석의 일종의 3차원 접수구조물의 진동해석을 효과적으로 수행하기 위한 해석방법을 제시하기 위하여 동적재해석기법을 검토하였다. 접수구조물의 유한구조 상호작용해석 결과는 구조진동의 관심 주파수역에서는 3차원 연성 부가수질량으로 표현되는 관성력으로 나타난다. 따라서 구조질량행렬에 부가수질량 행렬이 더해져서 전체 관성력으로 표현된다. 이 부가수질량을 추가질량으로 보고 재해석기법을 응용하는 방법을 수치실험을 통해 검증하였다. 이 때 재해석기법이 갖추어야 할 조건은 원구조의 질량과 거의 같은 정도의 질량이 추가되고 또한 완전 연성질량이 추가된 경우에도 정확한 해를 주어야 한다는 것이다. 이를 검증하기 위해 직접재해석기법과 섭동법을 이용한 재해석기법으로 4질량 스프링지지구조에 대한 해석을 수행한 결과 직접재해석기법의 응용이 적합함을 쉽게 입증할 수 있었다. 접수구조물의 예로는 3차원 잠수주상체에 대해 접수진동해석을 수행하였으며 그 결과 선체진동해석에 전통적으로 이용되고 있는 2차원 부가수질량과 3차원 수정계수를 사용한 기준차수법에서는 수지모드와 수평-비틔 연성모드와 같이 서로 독립적인 모드에 대해서는 따로 진동해석을 수행해 주어야 하는 단점이 발견되었다. 이 단점을 보완한 각 모드의 3차원 수정계수행렬을 이용한 재해석기법을 도입하여 모드에 상관없이 동시에 해를 구할 수 있었다. 그러나, 이 방법은 3차원 수정계수가 구해져 있는 경우에 한해서만 적용가능하며 실제 선체진동의 경우에는 10Hz 미만의 저차 주선체 진동에 한해서만 적용가능한 방법이다. 고차의 진도옴드에는 3차원 수정계수를 구할 수 없기 때문에 유체-구조 상호작용 해석결과로부터 얻은 3차원 연성 부가수질량을 이용하게 되며 이 때 이 행렬이 접수구조 표면의 전 자유도와 연성되어 있기 때문에 방대한 방정식을 푸어야 하지만 직접재해석기법을 적용함으로써 정확한 해를 구할 수 있었다. 또한 3차원 부가수질량을 이용한 직접재해석기법은 종래의 2차원 부가수질량과 3차원 수정계수를 이용한 방법에 비해 해석시간 면에서도 전혀 불리한 점이 없는 경제적 방법임이 밝혀졌다. 앞으로 Slamming 혹은 수중폭파 등의 충격하중에 의한 천이 구조응답 해석을 위한 효과적인 방법에 대해서도 연구결과를 발표할 계획이다.

• International Journal of Highway Engineering
• /
• v.10 no.4
• /
• pp.213-224
• /
• 2008

Flexible pavement responses to vehicular loading, such as critical stresses and strains, in each pavement layer, could be predicted by the multilayered elastic analysis. However, multilayered elastic theory suffers from major drawbacks including spatial dimension of a numerical model, material properties considered in the analysis, boundary conditions, and ill-presentation of tire-pavement contact shape and stresses. To overcome these shortcomings, three-dimensional finite element (3D FE) models are developed and numerical analyses are conducted to calculate pavement responses to moving load in this study. This paper introduces a methodology for an effective 3D FE to simulate flexible pavement structure. It also discusses the mesh development and boundary condition analysis. Sensitivity analyses of flexible pavement response to loading are conducted. The infinite boundary conditions and time-dependent history of calculated pavement responses are considered in the analysis. This study found that the outcome of 3D FE implicit dynamic analysis of flexible pavement that utilizes appropriate boundary conditions, continuous moving load, viscoelastic hot-mix asphalt model is comparable to field measurements.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.