• Journal of the Earthquake Engineering Society of Korea
• /
• v.15 no.5
• /
• pp.35-44
• /
• 2011

PGA (Peak Ground Acceleration) is the parameter which indicates the peak value for strong ground motion and is mainly due to the intensity of the seismic wave. Usually, seismic waves can consist of different characteristics and can have different effects on structures. Therefore, it may be undesirable that the effects of a seismic wave are evaluated only based on the PGA. In this study, time history analysis was executed with a single degree of freedom model for inelastic seismic analysis. The numerical model was assumed to be a perfect elasto-plastic model. Input accelerations were made with El Centro NS (1940), other earthquake records and artificial earthquakes. The displacement ductility demand and cumulative dissipated energy, which were calculated from other artificial earthquakes, were compared. As a result, different responses from other seismic waves which have the same PGA were identified. Therefore, an index which could reflect both seismic and structural characteristics is needed. The SI (Spectrum Intensity) scale which could be obtained from integration by parts of the velocity response spectrum could be an index reflecting the inelastic seismic response of structures. It can be possible to identify from correlation analysis among the SI scale, displacement ductility demand and cumulative dissipated energy that the SI scale is sufficient to be an index for the inelastic response of structures under seismic conditions.

• Journal of the Korean Geophysical Society
• /
• v.3 no.3
• /
• pp.161-174
• /
• 2000

Bayesian inversion is a stable approach to infer the subsurface structure with the limited data from geophysical explorations. In geophysical inverse process, due to the finite and discrete characteristics of field data and modeling process, some uncertainties are inherent and therefore probabilistic approach to the geophysical inversion is required. Bayesian framework provides theoretical base for the confidency and uncertainty analysis for the inference. However, most of the Bayesian inversion require the integration process of high dimension, so massive calculations like a Monte Carlo integration is demanded to solve it. This method, though, seemed suitable to apply to the geophysical problems which have the characteristics of highly non-linearity, we are faced to meet the promptness and convenience in field process. In this study, by the Gaussian approximation for the observed data and a priori information, fast Bayesian inversion scheme is developed and applied to the model problem with electric well logging and dipole-dipole resistivity data. Each covariance matrices are induced by geostatistical method and optimization technique resulted in maximum a posteriori information. Especially a priori information is evaluated by the cross-validation technique. And the uncertainty analysis was performed to interpret the resistivity structure by simulation of a posteriori covariance matrix.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.351-351
• /
• 2011

최근 지구온난화와 더불어 이상기후가 대두됨에 따라 기상 예측이 더욱더 중요시되고 있다. 또한 이전부터 가뭄 및 홍수와 같은 기상현상으로 인한 피해 사례가 빈번하였으며, 이로 인하여 물 관리의 어려움을 겪고 있다. 한 예로 이상기후가 유난히 잦았던 2010년 여름철 경우 평년보다 발달한 북태평양고기압의 영향으로 여름철 92일 가운데 81일의 전국 평균기온이 평년보다 높게 나타났다. 또한 강우 일수가 평년에 비해 7.4일 많은 44.2일을 기록하였으며, 국지성 집중호우 사례가 빈번하였다. 또한 8월 9일 발생한 태풍 `뎬무'를 포함해서 한 달 동안 3개의 태풍이 한반도에 영향을 끼치는 이례적인 사례가 발생하였다. 따라서 본 연구는 이러한 기상재해에 따른 물 관리를 장기적으로 대비하고자 고해상도 전지구 모델 GME를 이용하여 2010년 여름철 강수 예측을 실시하였다. 강수 예측에 사용된 전지구 모델 GME는 기존의 카테시안 격자체계를 가진 모델과 달리 전구를 삼각형으로 구성된 20면체로 격자화 한 Icosahedral-hexagonal grid 격자체계로 구성되어 있어, 해상도 증가에 용이할 뿐만 아니라, HPC(High Performance Computing)환경에서 효율성이 높은 장점을 가지고 있다. 본 계절 예측을 수행함에 있어 발생하는 잡음을 최소화하고자, Time-lag 기법을 이용하여 5개의 앙상블 멤버로 구성되어있으며, 이를 비교 분석하기위해 Climatology를 이용하여 총 10개의 앙상블 멤버로 규준실험을 수행하였다. 선행 연구에 따르면 1개월 이상의 장기 적분의 경우 초기조건보다 외부 강제력이 더 중요한 역할을 한다고 연구된 바 있다. (Yang et al., 1998) 특히 계절 변동성의 경우 대기-해양간의 상호작용에 의해 지배되며, 이를 고려하여 본 연구는 해수면 온도를 경계 자료로 사용하여 계절 예측을 수행하였다. 앞서 말한 실험 계획을 바탕으로 하여 나온 결과를 통해 동아시아지역 및 한반도 도별 강수 및 온도 변수에 대해 순별 및 월별 카테고리맵 분석을 실시하여 한눈에 보기 쉽게 나타냈다. 또한 주요 도시별 강수량 및 온도의 시계열 분석을 실시하여 시간이 지남에 따라 나타나는 변동성을 확인하였다. 계절 예측 결과에서 온도의 경우 평년보다 높게 나타났으며, 이는 실제 온도 예측과도 유사한 패턴을 가졌다, 강수의 경우 7월부터 8월 중순까지 평년보다 다소 적게 모의되었으며, 8월 하순경 회복하는 것으로 예측하였다. 따라서 본 계절 강수 예측은 다소 역학 모델이 가지는 한계를 가지고 있으나, 실제와 비교하여 어느 정도의 경향성이나 패턴에 있어 유사성을 보임을 확인하였으며, 이를 장기적 차원의 물관리를 함에 있어 참고 및 활용 가능할 것으로 예상한다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.1
• /
• pp.9-18
• /
• 2013

Two-phase vapor-liquid flows exist in many shell and tube heat exchangers such as condensers, evaporators, and nuclear steam generators. To understand the fluid dynamic forces acting on a structure subjected to a two-phase flow, it is essential to obtain detailed information about the characteristics of a two-phase flow. The characteristics of a two-phase flow and the flow parameters were introduced, and then, an experiment was performed to evaluate the pressure loss in the tube bundles and the fluid-dynamic force acting on the cylinder owing to the pressure distribution. A two-phase flow was pre-mixed at the entrance of the test section, and the experiments were undertaken using a normal triangular array of cylinders subjected to a two-phase cross-flow. The pressure loss along the flow direction in the tube bundles was measured to calculate the two-phase friction multiplier, and the multiplier was compared with the analytical value. Furthermore, the circular distributions of the pressure on the cylinders were measured. Based on the distribution and the fundamental theory of two-phase flow, the effects of the void fraction and mass flux per unit area on the pressure coefficient and the drag coefficient were evaluated. The drag coefficient was calculated by integrating the measured pressure on the tube by a numerical method. It was found that for low mass fluxes, the measured two-phase friction multipliers agree well with the analytical results, and good agreement for the effect of the void fraction on the drag coefficients, as calculated by the measured pressure distributions, is shown qualitatively, as compared to the existing experimental results.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.10a
• /
• pp.111-121
• /
• 1995

In Extracting the dynamic parameters for estimating the load carrying capacity and integrity of bridges, both the instrumentation and the processing the data plays important role . When the fixed point can not be secured, it is difficult and expensive to measure dynamic displacements. Even if the displacement is obtained through the integration of the acceleration data, the results can be quite different from the real behavior, because the main frequency contents can be leaked during discretized data processing. The instrumentation is used for measurements, and every measurement involves error and uncertainty, such as systematic, conformance, environmental, observational, sampling, and ranmom error. Systematic and conformance error can be remedied through the proper sellection and installation of the instruments, but sampling and random errors could not have been corrected properly and it becomes the limitation for using acceleration data. In this paper, the errors which can be occurred in numerical processing of dynamic data are referred, and the method to sellect proper sampling rate for the structural frequency range are proposed. Using the proposed method, the displacement response of the structures can be economically obtained from the measured acceleration record, and this procedure can be used properly to estimate the integrity of the bridges and infrastructures subjected to dynamic loads.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.199-206
• /
• 2014

In this paper, a scheme to evaluate the response variability for functionally graded material (FGM) beam with varying sectional area is presented. The randomness is assumed to appear in a spatial domain along the beam axis in the elastic modulus. The functionally graded material categorized as composite materials, however without the drawbacks of delamination and occurrence of cracks due to abrupt change in material properties between layers in the conventional composite materials. The functionally graded material is produced by the gradual solidification through thickness direction, which endows continuous variation of material properties, which makes this material performs in a smooth way. However, due to difficulties in tailoring the gradients, to have uncertainty in material properties is unavoidable. The elastic modulus at the center section is assumed to be random in the spatial domain along the beam axis. Introducing random variables, defined in terms of stochastic integration, the first and second moments of responses are evaluated. The proposed scheme is verified by using the Monte Carlo simulation based on the random samples generated employing the spectral representation scheme. The response variability as a function of correlation distance, the effects of material and geometrical parameters on the response variability are investigated in detail. The efficiency of the proposed scheme is also addressed by comparing the analysis time of the proposed scheme and MCS.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Journal of Navigation and Port Research
• /
• v.34 no.3
• /
• pp.195-203
• /
• 2010

This paper concerns about total dynamic analysis of integrated mining system. This system consists of vertical steel pipe, intermediate buffer station, flexible pipe and self-propelled miner. The self-propelled miner and buffer are assumed as rigid-body of 6-dof. Discrete models of vertical steel pipe and flexible pipe are adopted, which are obtained by means of lumped-parameter method. The motion of mining vessel is not considered. Instead, the motion of mining vessel is taken into account in form of various boundary conditions (e.g. forced excitation in slow motion and/or fast oscillation and so on). A terramechanics model of extremely cohesive soft soil is applied to the self-propelled miner. Hinged and ball constraints are used to define the connections between sub-systems (vertical steel pipe, buffer, flexible pipe, self-propelled miner). Equations of motion of the coupled model are derived with respect to the each local coordinates system. Four Euler parameters are used to express the orientations of the sub-systems. To solve the equations of motion of the total dynamic model, an incremental-iterative formulation is employed. Newmark-${\beta}$ method is used for time-domain integration. The total dynamic responses of integrated mining system are investigated.

• Elastomers and Composites
• /
• v.34 no.2
• /
• pp.135-146
• /
• 1999

The thermal degradation kinetics of SBR and tire were studied using a conventional thermogravimetric analysis in the stream nitrogen at a heating rate of 5, 10, 15, $20^{\circ}C/min$, respectively. Thermogravimetric curves and their derivatives were analyzed using various analytical methods to determine the kinetic parameters. The degradation of the SBR and tire was found to be a complex process which has multi-stages. The Friedman method gave average activation energies for the SBR and tire of 247.53kJ/mol and 230.00kJ/mol, respectively. Mean-while, the Ozawa method Eave 254.80kJ/mol and 215.76kJ/mol. It would appear that either. Friedman's differential method or Ozawa's integral method provided satisfactory mathematical approaches to determine the kinetic parameters for the degradation of the SBR and tire. Approximately 86% and 55% of oil products were obtained at a final temperature of $700^{\circ}C$ and a heating rate of $20^{\circ}C/min$ for the SBR and tire respectively.