• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.4
• /
• pp.1349-1360
• /
• 2013

The post-tensioned anchorage zones of normal concrete have larger cross sections because of congested reinforcements to resist high bearing and bursting stresses. The high compressive and tensile strength of newly developed UHPC (Ultra High Performance Concrete) may reduce the cross sectional dimensions and simplify the reinforcement details, if used for post-tensioned members. The Finite Element Analysis was performed to evaluate the mechanical behavior of post-tensioned anchorage zones using UHPC without anchorage plates and confining reinforcements. The results show that the maximum bursting stresses are less than the values given in current design code without failure due to vertical cracks. The location of maximum bursting stresses were at 0.2 times of width of the models. The bursting force from FEA is less than that is obtained using simplified formular in Korean Bridge Design Code.

• Economic and Environmental Geology
• /
• v.27 no.2
• /
• pp.191-199
• /
• 1994

We have developed an algorithm based on the method of integral equations to simulate the magnetotelluric (MT) responses of three-dimensional (3-D) bodies in a layered half-space. The inhomogeneities are divided into a number of cells and are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric Green's tensor function appropriate to a layered earth, and is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic Green's tensor functions over the scattering current About a 3-D conductive body near the earth's surface, interpretation using 2-D transverse electric modeling schemes can imply highly erratic low resistivities at depth. This is why these routines do not account for the effect of boundary charges. However, centrally located profiles across elongate 3-D prisms may be modeled accurately with a 2-D transverse magnetic algorithm, which implicitly includes boundary charges in its formulation. Multifrequency calculations show that apparent resistivity and impedance phase are really two complementary parameters. Hence, they should be treated simultaneously in broadband MT interpretation.

• Journal of Korea Water Resources Association
• /
• v.47 no.11
• /
• pp.1077-1086
• /
• 2014

The flood wave generated due to dam break is affected by initial depth upstream since it is related with hydraulic characteristics propagating downstream, and flow resistance stress has influence on the celerity, travel distance, and approaching depth of shock wave in implementing numerical simulation. In this study, a shallow water flow model employing SU/PG scheme was developed and verified by analytic solutions; propagation characteristics of dam break according to flow resistance and initial depth were analyzed. When bottom frictional stress was applied, the flow depth was relatively higher while the travel distance of shock wave was shorter. In the case of Coulomb stress, the flow velocity behind the location of dam break became lower compared with other cases, and showed values between no stress and turbulent stress at the reach of shock wave. The value of Froude number obtained by no frictional stress at the discontinuous boundary was the closest to 1.0 regardless of initial depth. The adaption of Coulomb stress gave more appropriate results compared with turbulent stress at low initial depth. However, as the initial depth became increased, the dominance of flow resistance terms was weakened and the opposite result was observed.

• 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.21 no.1
• /
• pp.1-7
• /
• 2018

In the two-dimensional (2D) modeling of electrical method, the potential in the space domain is reconstructed with the calculated potentials in the wavenumber domain using inverse Fourier transform. The inverse Fourier integral is numerically evaluated using the transformed potential at different wavenumbers. In order to improve the precision of the integration, either the logarithmic or exponential approximation has been used depending on the size of wavenumber. Two numerical methods have been generally used to evaluate the integral; interval integration and Gaussian quadrature. However, both methods do not consider the distance from the current source. Thus the resulting potential in the space domain shows some error. Especially when the distance from the current source is very small or large, the error increases abruptly and the evaluated potential becomes extremely unstable. In this study, we developed a new method to calculate the integral accurately by introducing the distance from the current source to the rescaled Gauss abscissa and weight. The numerical tests for homogeneous half-space model show that the developed method can yield the error level lower than 0.4 percent over the various distances from the current source.

• 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 Korea Water Resources Association
• /
• v.31 no.6
• /
• pp.757-768
• /
• 1998

rate of outflow is equated to the total rate of flow, both the state equation and its linearized equation yield almost identical oscillations of water levels. This shows that effects of pipe resistance on the surges are small, and justifies a free oscillation analysis. Free oscillation equation has six eigen modes of different periods, including a rigid mode which is existed when the pumped rate of outflow differs from the total rate of inflow. The six modes constitute a variety of surges dependent on different inflow and outflow conditions. Presence of the rigid mode needs sophisticated pump operations adjusted to real flood inflows.

• Magazine of the Korea Concrete Institute
• /
• v.5 no.2
• /
• pp.181-188
• /
• 1993

In this study, a numerical method for the material nonlinear analysis of reinforced concrete shell structures including the time dependent effects due to creep and shrinkage is developed. Degenerate shell elements with the layered approach are used. The perfect or strain hardening plasticity model in compression and the linearly elastic model in tension until cracking for concrete are employed. The reinforcing bars are considered as a steel layer of equivalent thickness. Each :steel layer has an uniaxial behaviour resisting only the axial force in the bar direction. A bilinear idealization is adopted to model elasto-plastic stress-strain relationships. For the nonlinear anaysis, incremental load method combined with unbalanced load iterations for each load increment is used. To include time dependent effects of concrete, time domain is divided into several time steps which may have different length. Some numerical examples are presented to study the validity and applicability of the present method. The results are compared with experimental and numerical results obtained by other investigator.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.327-336
• /
• 1997

When a structure is built on the consolidation ground, the instant elastic deflection occures according to the characteristics of the ground and the load on it. And the corresponding contact pressure is established. But, as time passes, the secondary consolidating deflection is added to the instant elastic deflection, the upper structure, due to its flexural rigidity, resist to the additional curvature. So the variation of the contact pressure occurs. And this new contact pressure exerts influence on the consolidation form again. The new consolidation form exerts influence on the contact pressure in return. This kind of interaction continues till all the consolidation of the ground is finished. So the consolidation problem can not be definded as the linear problem. This paper intends to scheme an approximate iteration method to analyse this non-linear interaction problem between the upper structure and the lower consolidation ground which supports the former.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.12 no.1
• /
• pp.1-6
• /
• 1979

본연구에서는 그물어구의 상사를 지배하는 무차원수 K를 $$K=\frac{{\nu}^n\rho_wv^{2-n}}{{d^{1+n}(\rho-\rho_w)}$$ d, p: 재료의 직경 및 밀도 $\nu,\rho_w,v$: 물의 동점성계수, 밀도 및 속도으로 정하고, 여기에서의 직경의 비를 결정하는 방법에 따라 실물과 모형과의 상사를 완전하게 그리고 근사적으로 만족시키는 조건들을 구하였다. 즉, 원전한 상사한 경우는 직경의 비를 축척비와 같게 하고, 나아가서 다른 모든 치수의 비도 축척비와 같게 함으로써 만족된다고 하였으며, 측사적 상사의 경우느 직경의 비가 축척비 $(\frac{\lambda_2}{\lambda_1})$와 같지 않아도 된다고 하여, 그물실의 직경 d, 코의 크기 $\iota$ 및 콧수 N의 비를 $$\frac{d_2}{d_1}=\frac{\iota_2}{\iota_1}=\frac{\lambda_2}{\lambda_1}{\cdot}\frac{N_1}{N_2}$$ 으로, 줄의 직경 d', 길이 $\iota'$ 및 밀도 $\rho'$의 비를 $$\frac{d_2'}{d_1'}=\sqrt{{\frac{\lambda_2}{\lambda_1}}\cdot{\frac{d_2}{d_1}}\cdot{\frac{\rho_2-\rho_{w2}}{\rho_1-\rho_{w1}}\cdot{\frac{\rho_1'-\rho_{w1}}{\rho_2'-\rho_{w2}}}}$$, $\frac{\iota_2'}{\iota_1'}=\frac{\lambda_2}{\lambda_1}$로, 부속구의 치경 $d'$, 밀도 $\rho'$ 및 수 $N'$의 비를 $$\frac{N_2'}{N_1'}=(\frac{\lambda_2}{\lambda_1})^2(\frac{d_2}{d_1})(\frac{d_1'}{d_2'})\frac{(\rho_2-\rho_{w2})}{(\rho_1-\rho_{w1})}\frac{(\rho_1'-\rho_{w1})}{(\rho_2'-\rho_{w2})}$$으로 정하였다. 이렇게 정해진 모형어구에 대해 유속 v의 비느 $K_1=K_2$로부터 $$(\frac{u_2}{u_1})^{2-n}=(\frac{\nu_2}{\nu_1})^{-n}\;(\frac{\rho_{w1}}{\rho_{w2}})\;(\frac{\rho_2-\rho_{w2}}{\rho_1-\rho_{w1}})\;(\frac{d_2}{d_1})^{1+n}$$으로 주어지므로, 이를 이용하여 어구저항 D 및 그물감의 다리에서의 장력 $\tau$의 비를 $$\frac{D_2}{D_1}=\frac{d_2(\rho_2-\rho_{w2})}{d_1(\rho_1-\rho_{w1})}(\frac{\lambda_2}{\lambda_1})^2$$ $${\frac{\tau_2}{\tau_1}=\frac{d_2\iota_2(\rho_2-\rho_{w2})}{d_1\iota_1(\rho_1-\rho_{w1})}\;{\cdot}\frac{\lambda_2}{\lambda_1}$$로 정하였다.