• Proceedings of the Korea Society for Simulation Conference
• /
• 1997.04a
• /
• pp.114-114
• /
• 1997

항만개발 사업은 매립이나 준설과 같은 해안선과 해저지형을 변화시키는 공사가 이루어지며 이는 조석조류 현상의 변화를 유발하게 된다. 특히 해안지형이 복잡한 항만이나 하구에서는 변화정도가 크게 나타나며 방파제나 선박접안 시설과 같은 항내 구조물을 축조 할 경우에도 조석조류 변화가 야기된다. 조석조류 현상의 변화는 해양 수질과 생물상에 영 향을 줄뿐만 아니라 인접한 육지에 이르기까지 크다란 피해를 유발할 수가 있다. 창조시 해 수를 육지로 범람하게 하거나 퇴적이나 침식조건을 변화시킬 수도 있으며 오염물질이 외해 로 확산되는 과정을 방해하고 해양생물의 서식조건을 악화시켜 부근 해역의 어업에 피해를 유발할 수도 있다. 또한 해안 공사 과정에 해수로 유입되는 토사나 항만 이용시에 배출되는 폐수 등이 항내 수질과 생태계를 변화시킬 수 있다. 따라서 조석조류와 수질변화를 사전에 예측하고 이를 최소화하는 방안이 강구되어야 한다. 조석조류와 수질 변화를 사전에 예측하 는 방법으로 축소모형 실험과 컴퓨터 시뮬레이션이 사용된다. 축소모형 실험은 현장을 재현 하는데 많은 경비와 시간이 소모되며 예측에 한계가 있다. 반면에 컴퓨터 시뮬레이션은 대 상해역의 지형과 조류 그리고 수질특성에 적합하게 모델을 개발할 경우 예측력이 매우 높기 때문에 현재 가장 널리 이용되고 있다. 본 연구는 유한 요소법에 기초하여 미국 공병단 (United States Army Corps of Engineers)에서 개발한 TABS-2 모델을 이용하여 온산항에 서 시행되는 부두건설에 따른 조석조류 및 수질변화를 시뮬레이션하였다. 측정된 자료로 모 델 계수를 보정하고 검증하였으며 항내 매립과 준설이 이루어질 때 나타나는 조위와 조속변 화를 예측하였다. 매립과 준설공사시 유출되는 토사가 조석조류현상에 따라 이동 확산되는 범위와 농도분포를 예측하였다.ophobic components)은 75~80%, 친수성 물질(hydrophilic components)은 30~33%정도의 분포를 보였고, 전염소 및 전오존 공정을 거친 처리수에서는 각각 62.2-62.8%, 43.9~49.0% 및 50~ 55%, 40~57% 정도의 분포를 보였다. 그리고 웅집-침전을 거친 처리수에서는 그 분포가 77~82%, 24-48%였다. 전주리 공정을 통하여 소수성 물질(byoghobic components)의 분포가 감 소하는 것을 볼 때 전염소 및 전오존 처리가 용존유기물의 응집에는 오히려 역 효과를 나타내는 것으로 판단된다. 것으로 판단된다.여 황토 2g에 대하여 Hieltijes and Lijklema 방법에 의해 Adsorbed-p, Nonapatite inorganic-P(NAI-P), Apatite-p, Organic-P로 구분하여 분석하고, 총인(Total Phosphorus)을 Standard Methods에 따라 Persulfate digestion후 0.45 m membrane 여지 여과하여 여액에 대해 PO3-4-P의 농도를 Ascorbic Acid 법으로 측정한 결과, NAI-P가 가장 큰 비율을 차지하였고, 부원료로 첨가된 금속 양이온 중 Fe3-이온이 흡착에 기여하는 정도가 가장 큰 것으로 평가되었다.다당류 T-AS의 보체 활성화 기작은 classical과 alternative complement pathway의 양 경로를 통해 활성화 되었다. T-AS 분획은 mouse내의 특정 혈청단백을 증가시켰으며, 항체 생성능의 증가가 관찰되어 effect T 세포의 활성화가 나타나고 있음을 알 수 있었다. T-AS는 생체내 투여시에 대식세포의 탐식능이 증진되었으며, 대식세포 기능 저해제에 의한 대식세포의 기능 저해 현상이 회복되었다. 이와 같은 결과들로부터, T-AS의 항암 활성은 활성

• Journal of the Korea Concrete Institute
• /
• v.21 no.1
• /
• pp.13-20
• /
• 2009

The influence of steel fiber volume on the tension softening behavior in steel fiber-reinforced ultra high strength concrete was investigated. Three-point bending test (TPBT) with notched beams was performed and inverse analysis method by Uchida et al. was adopted to obtain the tension softening behaviors from the results of TPBT. It could be found that the intial stiffness was constant regardless of steel fiber volume, the increase of steel fiber volume fraction made the tensile strength higher, but all of the curves converged on an asymptote with a crack width. It was proposed the equation of softening curve expressed by combination of plastic behavior part and exponential descending behavior part considering the steel fiber volume fraction and $\omega_0$, which is corresponding to the maximum crack width of plastic area. Thereafter, the crack propagation analysis using finite element method with smeared crack model was also carried out and it was confirmed that the proposed equation had a good agreement with the experimental results.

• Journal of the Korean Geotechnical Society
• /
• v.18 no.6
• /
• pp.5-16
• /
• 2002

A hydrogeomechanical numerical model is presented to evaluate rainfall impacts on groundwater flow in slopes and slope stability. This numerical model is developed based on the fully coupled poroelastic governing equations for groundwater flow in deforming variably saturated geologic media and the Galerkin finite element method. A series of numerical experiments using the model developed are then applied to an unsaturated slope under various rainfall rates. The numerical simulation results show that the overall hydromechanical slope stability deteriorates, and the potential failure nay initiate from the slope toe and propagate toward the slope crest as the rainfall rate increases. From the viewpoint of hydrogeology, the pressure head and hence the total hydraulic head increase as the rainfall rate increases. As a result, the groundwater table rises, the unsaturated zone reduces, the seepage face expands from the slope toe toward the slope crest, and the groundwater flow velocity increases along the seepage face. From the viewpoint of geomechanics, the horizontal displacement increases, and the vertical displacement decreases toward the slope toe as the rainfall rate increases. This may result from the buoyancy effect associated with the groundwater table rise as the rainfall rate increases. As a result, the overall deformation intensifies toward the slope toe, and the unstable zone, in which the factor of safety against shear failure is less than 1, becomes thicker near the slope toe and propagates from the slope toe toward the slope crest. The numerical simulation results also suggest that the potential tension failure is likely to occur within the slope between the potential shear failure surface and the ground surface.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.36D no.3
• /
• pp.85-91
• /
• 1999

$H:LiNbO_3$ optical modulator with Cu/parylene electrode layer, which has a merits in the bandwidth, power consumption and fabrication conditions as compared with conventional Au/Cr/$SiO_2$, is proposed and fabricated. Analysis and design of optical modulator is performed by finite element calculation. Various unit processes for fabricating the proposed modulator, 1550nm $H:LiNbO_3$ optical waveguide, parylene buffer layer, and CPW Cu electrode, were developed, After dicing and end-face polishing of fabricated modulator chip, optical modulation responses as sawtooth electrical driving voltage has been measured at low frequencies. Properties of optical waveguide had not been changed before and after Cu/parylene electrode processes, which make confirm the reproducible fabrication of optical modulator.

• Transactions of Materials Processing
• /
• v.15 no.6 s.87
• /
• pp.445-452
• /
• 2006

In most of automobile body panel manufacturing, trimming process is generally performed before flanging. To find feasible trimming line is crucial in obtaining accurate edge profile after flanging. Section-based method develops blank along manually chosen section planes and find trimming line by generating loop of end points. This method suffers from inaccurate results of edge profile. On the other hand, simulation-based method can produce more accurate trimming line by iterative strategy. In this study, new fast simulation-based method to find feasible trimming line is proposed. Finite element inverse method is used to analyze the flanging process because final shape after flanging can be explicitly defined and most of strain paths are simple in flanging. In utilizing finite element inverse method, the main obstacle is the initial guess generation for general mesh. Robust initial guess generation method is developed to handle genera] mesh with very different size and undercut. The new method develops final triangular mesh incrementally onto the drawing tool surface. Also in order to remedy mesh distortion during development, energy minimization technique is utilized. Trimming line is extracted from the outer boundary after finite element inverse method simulation. This method has many advantages since trimming line can be obtained in the early design stage. The developed method is verified by shrink/stretch flange forming and successfully applied to the complex industrial applications such as door outer flanging process.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.12
• /
• pp.1345-1350
• /
• 2014

This paper presented analysis methods for adapting E-glass fiber/epoxy composite (GFRP) materials to an automotive leaf spring. It focused on the static behaviors of the leaf spring due to the material composition and its fiber orientation. The material properties of the GFRP composite were directly measured based on the ASTM standard test. A reverse implementation was performed to obtain the complete set of in-situ fiber and matrix properties from the ply test results. Next, the spring rates of the composite leaf spring were examined according to the variation of material parameters such as the fiber angles and resin contents of the composite material. Finally, progressive failure analysis was conducted to identify the initial failure load by means of an elastic stress analysis and specific damage criteria. As a result, it was found that damage first occurred along the edge of the leaf spring owing to the shear stresses.

• Geotechnical Engineering
• /
• v.5 no.1
• /
• pp.47-60
• /
• 1989

Lateral earth Pressure distributions due to the ,randy soil backfill behind the rigid vertical walls for three different wall movement modes are obtained by the elasto-plastic finite element analys of soil deformation, and these earth pressures are compared with both Rankine's and Dubrova's active earth pressures. Thereby, the effects of the magnitude and the mode of wall displacement on the earth pressure distribution are investigated. Three different modes of wall movement considered in this study are the rotation about bottom, the rotation about top and the translation. For the case of the wall rotation about top, the earth pressure distribution is shown as a reverse S-curve-shaped distribution due to the arching effect. Consequently, the point of application of the lateral thrust is much higher than one-third of the wall height from the base. And, comparing the other modes of wall movement, the magnitude and the point of appliestion of the lateral thrust for the wall rotation about top are larger and higher, respectively. The wedge-shaped plastic zone in the backfill at active failure is developed only for the mode of wall rotation about bottom. The lateral earth pressure distributions on the walls with inclined backfill of several different slopes are shown for the mode of wall rotation about bottom.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.2 no.3
• /
• pp.1-13
• /
• 1982

Transition process of plastic region. displacements, stresses and strains ahead the flaw tips were analysed by the finite element method on the steel plate with the circular hole and the one with the elliptical hole under completely alternating load (repetition of tensile loading, unloading and compressive loading). As the results, the followings were obtained. Transition process of elastic failure (yielding) region was estimated. From this the tendency was confirmed that the fracture would be initiated from ahead the flaw tip, and propagated along the $45^{\circ}$ direction. The fundamental data available in estimating the stress intensity factor that was considered as the core in analysing the fracture mechanism of steel plates were obtained. It was indicated that when unloading after tension the effect of compressive loading, and even the compressive reyield, was occured ahead the flaw tip. Similarly it was indicated that when unloading after compression the effect of tensile loading, and even the tensile reyield, was occured ahead the flaw tip. It was considered that these phenomena were occured because the unloading effect was constrained by the residual strains when unloading. It was considered that the fatigue phenomenon was occured ahead, the flaw tip by repetition of tensile yield, the above compressive reyield, compressive yeild and the above tensile reyield. In addition, the tendency was confirmed that the fracture ahead the flaw tip was occured as the flaw was changed from the circular hole to the elliptical hole and became to be the crack lastly.

• Geophysics and Geophysical Exploration
• /
• v.1 no.2
• /
• pp.127-135
• /
• 1998

The Born approximation is widely used for solving the complex scattering problems in electromagnetics. Approximating total internal electric field by the background field is reasonable for small material contrasts as long as scatterer is not too large and the frequency is not too high. However in many geophysical applications, moderate and high conductivity contrasts cause both real and imaginary part of internal electric field to differ greatly from background. In the extended Born approximation, which can improve the accuracy of Born approximation dramatically, the total electric field in the integral over the scattering volume is approximated by the background electric field projected to a depolarization tensor. The finite difference and elements methods are usually used in EM scattering problems with a 2D model and a 3D source, due to their capability for simulating complex subsurface conductivity distributions. The price paid for a 3D source is that many wavenumber domain solutions and their inverse Fourier transform must be computed. In these differential equation methods, all the area including homogeneous region should be discretized, which increases the number of nodes and matrix size. Therefore, the differential equation methods need a lot of computing time and large memory. In this study, EM modeling program for a 2D model and a 3D source is developed, which is based on the extended Born approximation. The solution is very fast and stable. Using the program, crosshole EM responses with a vertical magnetic dipole source are obtained and the results are compared with those of 3D integral equation solutions. The agreement between the integral equation solution and extended Born approximation is remarkable within the entire frequency range, but degrades with the increase of conductivity contrast between anomalous body and background medium. The extended Born approximation is accurate in the case conductivity contrast is lower than 1:10. Therefore, the location and conductivity of the anomalous body can be estimated effectively by the extended Born approximation although the quantitative estimate of conductivity is difficult for the case conductivity contrast is too high.

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