• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.4
• /
• pp.734-742
• /
• 2016

Global warming and the depletion of fossil fuels have been caused by decades of reckless development. Wind energy is one form of renewable energy and is considered a future energy source. The wind tower is designed with a fundamental frequency in the soft-stiff design between the 1P and 3P range to avoid resonance. Usually, to perform natural frequency analysis of a wind tower, the boundary condition is set to the Fixed-End, and soil-pile interaction is not considered. In this study, consideration of the effect of soil-pile interaction on the wind tower was included and the difference in the natural frequency was studied. The fixed boundary condition was not affected by the soil condition and depth of the pile and the coupled spring boundary condition was unaffected by the depth of pile but affected by the depth of the pile, and the Winkler spring boundary condition is affected by both the soil condition and the depth of the pile. Therefore, the coupled spring boundary condition should be used in shallow depth soil conditions because the soil condition does not take the shallow depth soil into consideration.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.12 no.3
• /
• pp.165-172
• /
• 2009

A set of equations for description of transformation of harmonic waves is proposed here. Velocity potential function and separation of variables are introduced for the derivation. The continuity equation is in a vertical plane is integrated through the water so that a horizontal one-dimensional wave equation is produced. The new equation composed of the complex velocity potential function, further be modified into. A set up of equations composed of the wave amplitude and wave phase gradient. The horizontally one-dimensional equations on the wave amplitude and wave phase gradient are the first and second-order ordinary differential equations. They are solved in a one-way marching manner starting from a side where boundary values are supplied, i.e. the wave amplitude, the wave amplitude gradient, and the wave phase gradient. Simple spatially-centered finite difference schemes are adopted for the present set of equations. The equations set is applied to three test cases, Booij's inclined plane slope profile, Massel's smooth bed profile, and Bragg's wavy bed profile. The present equations set is satisfactorily verified against existing theories including Massel's modified mild-slope equation, Berkhoff's mild-slope equation, and the full linear equation.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.1
• /
• pp.29-39
• /
• 2012

The maximum unit point resistance ($q_{max}$) of rock socketed drilled shafts subjected to axial loads was investigated by a numerical analysis. A 3D Finite Difference Method (FDM) analysis and a Distinct Element Method (DEM) analysis were performed with varying rock elastic modulus (E), discontinuity spacing ($S_j$), discontinuity dip angle ($i_j$), and pile diameter (D). Based on the results of obtained, it was found that the ultimate point resistance ($q_{max}$) increased as rock elastic modulus (E) and rock discontinuity spacing ($S_j$) increased. But, it was found that $q_{max}$ decreased as pile diameter (D) increased. As for the influence of the dip angle of rock discontinuity ($i_j$), it was shown that $q_{max}$ decreased up to 50% of maximum value within the range of $0^{\circ}$ < $i_j$ < $60^{\circ}$ due to the shear failure at rock discontinuities. Furthermore, it was found that if $20^{\circ}{\leq}i_j{\leq}40^{\circ}$, influence of $i_j$ should be taken into account because $q_{max}$ tended to approach a minimum value as $i_j$ approached a value near the friction angle of the discontinuity (${\phi}_j$).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.9 no.1
• /
• pp.9-23
• /
• 1997

A numerical study of a two-layer, stratified flow is investigated, using the implicit finite difference method in one dimension. The results of computational method have been tested and, in case of lock exchange flow, compared with the results of experimental data. The results of model experiments with various interfacial, bottom friction coefficients along with various time weighting factor of numerical scheme and dissipative interface are shown and discussed. Two-layer model experiment has been also carried out to investigate the generation and propagation characteristics of internal tidal wave over the steep bottom topography under stratified condition. The internal wave seems to well radiate through the downstream boundary under the experiments adopting radiation conditions both at two layers and only at upper layer, confirming the applicability of radiational boundary condition in stratified flows. It is also shown that the internal wave through the downstream boundary propagates more actively with increasing thickness of lower layer in the downstream. This implies that the potential tidal energy in the interface will depend upon the thickness of lower layer for the constant thickness of upper layer.

• Journal of the Korean GEO-environmental Society
• /
• v.16 no.9
• /
• pp.33-41
• /
• 2015

The vertical bearing capacity of a bucket foundation installed in sand can be calculated as sum of the skin friction and end bearing capacity. However, the current design equations are not considering the non-associated flow characteristics of sand and the reduction in the skin friction and increase in the end bearing capacity when the vertical load is applied. In this study, we perform two-dimensional axisymmetric finite element analyses following non-associated flow rule and calculate the vertical bearing capacity of circular bucket foundation of various sizes installed in sand of different friction angles. After calculating the skin friction and end bearing force at the ultimate state, design equations are derived for each. The skin friction of bucket foundation is shown significantly small compared to the end bearing capacity. Considering the difference with the available design equation for piles, it is recommended that the equation for piles is used for the bucket foundation. A new shape-depth factor ($s_q{\cdot}d_q$) for bucket foundation is recommended which also accounts for the increment of the end bearing capacity due to skin friction. Additionally, the shape and depth factor of embedded foundation proposed from the associated flow rule can overestimate the bearing capacity in sand, so it is more adequate to use the shape-depth factor proposed in this study.

• The Journal of Korean Academy of Prosthodontics
• /
• v.37 no.4
• /
• pp.433-453
• /
• 1999

There are a lot of elements affecting the success of implant prosthesis. The quality of surgical procedure is considered as one of the key factors. To avoid the excess heat generation is an important element of successful osseointegration and it can be achieved by using a gentle surgical technique with a sharp instrument in bone drilling. This study was performed to measure and analyze comparatively the heat transmitted to sur-rounding bone at a distance of 0.5mm from the periphery of the drill hole in each drilling stage. The results were as follows. In standard system, the temperature of surrounding bone tissue ranged from $29.2^{\circ}C\;to\;48.3^{\circ}C$ with irrigation and from $34.6^{\circ}C\;to\;84.3^{\circ}C$ without irrigation. And in wide system, the temperature of surrounding bone tissue ranged from $29.5^{\circ}C\;to\;52.5^{\circ}C$ with irrigation and from $34.8^{\circ}C\;to\;87.8^{\circ}C$ without irrigation. And the temperature ranges exceeded the threshold without irrigation, while showing less than the threshold by the cooling effect of irrigation. In comparing standard system with wide system, although there was no significant difference, ${\phi}4.3mm$ pilot and ${\phi}4.3mm$ twist drill of wide system showed high value and wide system showed slightly high elevation of temperature in all depth in fixture installation. In the finite element analysis, the calculated value by the Fourier's cooling law were applied to the bone drilling surface. And through analysis using different irrigation temperatures at $28^{\circ}C,\;15^{\circ}C\;and\;5^{\circ}C$, and according to the time. The result was that the cooling water at least below $15^{\circ}C$ was required to maintain the temperature of surrounding bone less than threshold in bone drilling, the cooling water below $5^{\circ}C$ was required to gain more sufficient cooling effect, and cooling over 5 seconds was needed after bone drilling for sufficient effect.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.9 no.2
• /
• pp.149-158
• /
• 1998

The SAR patterns of biological objects contacted with coaxial waveguide antennal has been investigated, in which the biological object was modeled by a homogeneous and four-layered lossy human body. We derived the finite-difference time-domain(FDTD) algorithm and equation of MUR and generalized perfectly matched layer(GPML) ABCs in cylindrical coordination. The coupling between coaxial waveguide antenna and a biological object was analyzed by use of MUR and GPML ABCs in the FDTD method to obtain the absorbed power patterns in the media. The specific absorption rates (SAR) distribution which was corresponding to the temperature distribution was calculated in each region by use of the steady-state response in the FDTD method. The SAR patterns of the FDTD method using MUR absorbing boundary conditions(ABCs) was compared with those of the FDTD method using GPML ABCs. The comparison exhibits that the penetration depth of the SAR patterns using MUR ABCs is deeper than that of the SAR patterns using GPML ABCs because of loss in free space. However, the spread in the lateral directions of the SAR patterns using GPML ABCs is smaller than of the SAR patterns using MUR ABCs.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.2
• /
• pp.87-96
• /
• 1999

To analyze the dynamic behavior of structure, direct integration and mode superposition may be utilized in time domain analysis. As finite number of frequencies can give relatively exact solutions, mode superposition is preferable in analyzing structural behavior. In non-linear analysis, however, mode superposition is seldom used since time-varying element stiffness changes stiffness matrix, and the change of stiffness matrix leads to the change of essential constants - natural frequencies and mode shapes. In spite of these difficulties, there are some attempts to adopt mode superposition because of low cost compared to direct integration, but the result is not satisfactory. In this paper, a method using mode superposition in non-linear analysis is presented by separating local element stiffness from global stiffness matrix with the difference between linear and non-linear restoring forces to the external force vectors included. Moreover, the hysteresis model changing with the relative deformation in each floor makes it possible to analyze non-linear behavior of structure. The proposed algorithm is applied to shear beam model and the maximum displacement is compared with the result using direct integration method.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.35-47
• /
• 1996

A surface panel method treating a boundary-value problem of the Dirichlet type is presented to design a three dimensional body with free surface corresponding to a prescribed pressure distribution. An integral equation is derived from Green's theorem, giving a relation between total potential of known strength and the unknown local flux. Upon discretization, a system of linear simultaneous equations is formed including free surface boundary condition and is solved for an assumed geometry. The pseudo local flux, present due to the incorrect positioning of the assumed geometry, plays a role f the geometry corrector, with which the new geometry is computed for the next iteration. Sample designs for submerged spheroids and Wigley hull and carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. For the calculation of the wave resistance, normal dipoles and Rankine sources are distributed on the body surface and Rankine sources on the free surface. The free surface boundary condition is linearized with respect to the oncoming flow. Four-points upwind finite difference scheme is used to compute the free surface boundary condition. A hyperboloidal panel is adopted to represent the hull surface, which can compensate the defects of the low-order panel method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but by experiment.

• Geophysics and Geophysical Exploration
• /
• v.9 no.2
• /
• pp.139-147
• /
• 2006

We had conducted a numerical modeling to investigate seismic properties of gas hydrate with field parameters acquired over the East sea in 1998. We used a 2-D staggered grid finite difference method to generate synthetic elastic seismograms for multi-channel seismic survey, OBC (Ocean Bottom Cable) survey and VCS (Vertical Cable Seismic) survey. The results of this study showed that the method using staggered grid yielded stable results and could be used to seismic imaging. We could find out the high amplitude anomaly and the phase reversal phenomenon of reflection wave at interface between the gas hydrate layer and free gas layer such a BSR (Bottom Simulating Reflector) which is the evidence for existence of gas hydrate in seismic reflection data. And we computed the reflection coefficients at the incident angles corresponding to offset distance with the synthetic seismograms. The reflection coefficients acquired from the numerical modeling were nearly consistent with the reflection coefficient computed by Shuey's equation.