• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.105-113
• /
• 2009

The interpretation of observed waveform characteristics identified in refraction and wide-angle reflection data increases confidence in the crustal structure model obtained. When calculating traveltimes and raypaths, wavefront methods on a regular grid based on graph theory are robust even with complicated structures, but basically compute only first arrivals. In this paper, we develop new algorithms to compute traveltimes and raypaths not only for first arrivals, but also for fast and later reflection arrivals, later refraction arrivals, and converted waves between P and S, using the modified wavefront method based on slowness network nodes mapped on a multi-layer model. Using the new algorithm, we can interpret reflected arrivals, Pg-later arrivals, strong arrivals appearing behind Pn, triplicated Moho reflected arrivals (PmP) to obtain the shape of the Moho, and phases involving conversion between P and S. Using two models of an ocean-continent transition zone and an oceanic ridge or seamount, we show the usefulness of this algorithm, which is confirmed by synthetic seismograms using the 2D Finite Difference Method (2D-FDM). Characteristics of arrivals and raypaths of the two models differ from each other in that using only first-arrival traveltime data for crustal structure analysis involves risk of erroneous interpretation in the ocean-continent transition zone, or the region around a ridge or seamount.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.5
• /
• pp.612-620
• /
• 2019

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has results in considerable effort being directed toward developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have traditionally been obtained by the standard finite difference. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method(DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as excessive use of storage due to conditions of complex geometry and loading. In this study, the in-plane extensional vibration for asymmetric curved beams with linearly varying cross-section is analyzed using the DQM. Fundamental frequency parameters are calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results are compared with the result by other methods for cases in which they are available. According to the analysis of the solutions, the DQM, used only a limited number of grid points, gives results which agree very well with the exact ones.

• Journal of the Korean Geotechnical Society
• /
• v.39 no.12
• /
• pp.47-60
• /
• 2023

At present, in South Korea, there is a growing concern regarding solar power facilities installed on slopes because they are prone to damage caused by natural disasters, such as heavy rainfall and typhoons. Each year, these solar power facilities experience soil erosion due to heavy rainfall and foundation damage or detachment caused by strong wind loads. Despite these challenges, the interaction between the ground and structures is not adequately considered. Current analyses primarily focus on the structural stability under external loads; the overall facility site's stability-excluding the solar structures-in relation to its surrounding slopes is neglected. Therefore, in this study, we use finite-difference method analysis to simulate the behavior of the foundation and piles to assess changes in lateral displacement and bending stress in piles, as well as the safety factor of sloped terrains, in response to various influencing factors, such as pile diameter, spacing between piles, pile-embedding depth, wind loads, and dry and wet conditions. The analysis results indicate that pile spacing and wind loads significantly influence lateral displacement and bending stress in piles, whereas pile-embedding depth strongly influences the safety factor of sloped terrains. Moreover, we found that under certain conditions, the design criteria in domestic standards may not be met.

• Journal of the Korean Geotechnical Society
• /
• v.40 no.4
• /
• pp.49-60
• /
• 2024

Ground-penetrating radar (GPR) enables rapid data acquisition over extensive areas, but interpreting the obtained data requires specialized knowledge. Numerous studies have utilized numerical analysis methods to examine GPR signal characteristics under various conditions. To develop more realistic numerical models, the heterogeneous nature of the ground, which causes clutter, must be considered. Clutter refers to signals reflected by objects other than the target. The Peplinski material model and fractal techniques can simulate these heterogeneous characteristics, yet there is a shortage of research on the necessary input parameters. Moreover, methods for quantitatively evaluating the similarity between field and analytical data are not well established. In this study, we calculated the autocorrelation coefficient of field data and determined the correlation length using the autocorrelation function. The correlation length represented the temporal or spatial distance over which data exhibited similarity. By comparing the correlation length of field data with that of the numerical model incorporating fractal weights, we quantitatively evaluated a numerical model for heterogeneous ground. Consequently, the results of this study demonstrated a numerical modeling technique that reflected the clutter characteristics of the field through correlation length.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.1
• /
• pp.53-66
• /
• 1980

The research work is concerned with the analytical and experimental studies on the heat transfer phenomenon around the underground concrete digester used for biogas production Systems. A mathematical and computational method was developed to estimate heat losses from underground cylindrical concrete digester used for biogas production systems. To test its feasibility and to evaluate thermal parameters of materials related, the method was applied to six physical model digesters. The cylindrical concrete digester was taken as a physical model, to which the model,atical model of heat balance can be applied. The mathematical model was transformed by means of finite element method and used to analyze temperature distribution with respect to several boundary conditions and design parameters. The design parameters of experimental digesters were selected as; three different sizes 40cm by 80cm, 80cm by 160cm and l00cm by 200cm in diameter and height; two different levels of insulation materials-plain concrete and vermiculite mixing in concrete; and two different types of installation-underground and half-exposed. In order to carry out a particular aim of this study, the liquid within the digester was substituted by water, and its temperature was controlled in five levels-35。 C, 30。 C, 25。 C, 20。C and 15。C; and the ambient air temperature and ground temperature were checked out of the system under natural winter climate conditions. The following results were drawn from the study. 1.The analytical method, by which the estimated values of temperature distribution around a cylindrical digester were obtained, was able to be generally accepted from the comparison of the estimated values with the measured. However, the difference between the estimated and measured temperature had a trend to be considerably increased when the ambient temperature was relatively low. This was mainly related variations of input parameters including the thermal conductivity of soil, applied to the numerical analysis. Consequently, the improvement of these input data for the simulated operation of the numerical analysis is expected as an approach to obtain better refined estimation. 2.The difference between estimated and measured heat losses was shown to have the similar trend to that of temperature distribution discussed above. 3.It was found that a map of isothermal lines drawn from the estimated temperature distribution was very useful for a general observation of the direction and rate of heat transfer within the boundary. From this analysis, it was interpreted that most of heat losses is passed through the triangular section bounded within 45 degrees toward the wall at the bottom edge of the digesten Therefore, any effective insulation should be considered within this region. 4.It was verified by experiment that heat loss per unit volume of liquid was reduced as the size of the digester became larger For instance, at the liquid temperature of 35˚ C, the heat loss per unit volume from the 0. 1m$^3$ digester was 1, 050 Kcal/hr m$^3$, while at for 1. 57m$^3$ digester was 150 Kcal/hr m$^3$. 5.In the light of insulation, the vermiculite concrete was consistently shown to be superior to the plain concrete. At the liquid temperature ranging from 15。 C to 350 C, the reduction of heat loss was ranged from 5% to 25% for the half-exposed digester, while from 10% to 28% for the fully underground digester. 6.In the comparison of heat loss between the half-exposed and underground digesters, the heat loss from the former was fr6m 1,6 to 2, 6 times as much as that from the latter. This leads to the evidence that the underground digester takes advantage of heat conservation during winter.

• The Journal of Engineering Geology
• /
• v.27 no.3
• /
• pp.331-343
• /
• 2017

In this study to improve stability, workability and economics of the H-Pile+Earth plate or H-Pile+Earth plate+Cutoff grouting currently in use, we had developed HCS method belonging to the retaining wall which is consisting of a combination H-Pile, Plastic Sheet Pile and Steel Square Pipe for gap maintenance and reinforcement of flexible plastic Sheet Pile, and the behavior of each member composing HCS method is investigated by three-dimensional finite element analysis. To numerically analyze the behavior of the HCS method, we have performed extensive three-dimentional finite element analysis for three kinds of plastic Sheet Pile size, two kinds of H-Pile size and three kinds of H-Pile installation interval, one kinds of Steel Square Pipe and three kinds of Steel Square Pipe installation interval. After analyzing the numerical results, we found that the combinations of $P.S.P-460{\times}131.5{\times}7t$ (PS7) and H-Pile $250{\times}250{\times}9{\times}14$ (H250), $P.S.P473{\times}133.5{\times}9t$ (PS9) and H-Pile $300{\times}200{\times}9{\times}14$ (H300) is the most economical because these combinations are considered to have a stress ratio (=applied stress/allowable stress) close to that as the stiffness of H-Pile, plastic Sheet Pile and Steel Square Pipe composite increased, the horizontal displacement of the retaining wall and the vertical displacement of the upper ground decreased. Especially, due to the arching effects caused by the difference in stiffness between H-Pile and plastic Sheet Pile, a large part of the earth pressure acting on plastic Sheet Pile caused a stress transfer to H-Pile, and the stress and displacement of plastic Sheet Pile were small. Through this study, we can confirm the behavior of each member constituting the HCS method, and based on the confirmed results of this study, it can be used to apply HCS method in reasonable, stable and economical way in the future.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.183-193
• /
• 2000

A series of studies such as geological logging data analysis, detailed geological survey, rock mass evaluation, in-situ and laboratory tests, rock strength and mechanical properties of the rock were concerned. The stability of the slope were carried out inorder to design the pit slope and individual benches using the stereographic projection analysis and numerical methods in Roto Pit of Pasir coal fetid. The bedding plane was one of the major discontinuities in the Roto Pit and the dip of which is about $60^{\circ}$in the northern part and $83^{\circ}$in the southern part. The dip of bedding becomes steeper from north to south. The plane and toppling failures are presented in many slopes. In laboratory test the average uniaxial compressive strength of mudstone was 9 MPa and that of weak sandstone was 10 MPa. In-situ test showed that the rocks of Roto north mining area are mostly weak enough to be classified in grade from R2(weak) to R3(medium strong weak) and the coal is classified in grades from R1(Very weak) to R2(Weak). The detailed stability analysis were carried out on 4 areas of Roto north(east, west, south and north), and 2 areas of Roto south(east and west). In this paper, the minimum factor of safety was set to 1.2 which is a general criterion for open pit mines. Using the stereographic projection analysis and the limit equilibrium method, slope angles were calculated as 30~$36^{\circ}$for a factor of safety greater than 1.2. Then these results were re-evaluated by numerical analysis using FLAC. The final slope angles were determined by rational described abode. A final slope of 34 degrees can guarantee the stability for the eastern part of the Roto north area, 33 degrees for the western part, 35 degrees for the northern part and 35 degrees for the southern part. For the Roto south area, 36 degrees was suggested for both sides of the pit. Once the pit slope is designed based on the stability analysis and the safety measures. the stability of 니ope should be checked periodically during the mining operations. Because the slope face will be exposed long time to the rain fall, a study such aspreventive measures against weathering and erosion is highly recommended to be implemented.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.430-440
• /
• 2000

A series of studies such as geological logging data analysis, detailed geological survey, rock mass evaluation, in-situ and laboratory tests, rock strength and mechanical properties of the rock were concerned. The stability of the slope were carried out inorder to design the pit slope and individual benches using the stereographic projection analysis and numerical methods in Roto Pit of Pasir coal field. The bedding plane was one of the major discontinuities in the Roto Pit and the dip of which is about 60$^{\circ}$ in the northern part and 83$^{\circ}$ in the southern part. The dip of bedding becomes steeper from north to south. The plane and toppling failures are presented in many slopes. In laboratory test the average uniaxial compressive strength of mudstone was 9MPa and that of weak sandstone was 10MPa. In-situ test showed that the rocks of Roto north mining area are mostly weak enough to be classified in grade from R2(weak) to R3(medium strong weak) and the coal is classified in grades from R1(Very weak) to R2(Weak). The detailed stability analysis were carried out on 4 areas of Roto north (east, west, south and north), and 2 areas of Roto south(east and west). In this paper, the minimum factor of safety was set to 1.2 which is a general criterion for open pit mines. Using the stereographic projection analysis and the limit equilibrium method, slope angles were calculated as 30∼36$^{\circ}$ for a factor of safety greater than 1.2. Then these results were re-evaluated by numerical analysis using FLAC. The final slope angles were determined by rational described above. A final slope of 34 degrees can guarantee the stability for the eastern part of the Roto north area, 33 degrees for the western part, 35 degrees for the northern part and 35 degrees for the southern part. For the Roto south area, 36 degrees was suggested for both sides of the pit. Once the pit slope is designed based on the stability analysis and the safety measures, the stability of slope should be checked periodically during the mining operations. Because the slope face will be exposed long time to the rain fall, a study such aspreventive measures against weathering and erosion is highly recommended to be implemented.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.51-60
• /
• 2010

To determine subglacial topography and internal features of the Fourcade Glacier on King George Island in Antarctica, helicopter-borne and ground-towed ground-penetrating radar (GPR) data were recorded along four profiles in November 2006. Signature deconvolution, f-k migration velocity analysis, and finite-difference depth migration applied to the mixed-phase, single-channel, ground-towed data, were effective in increasing vertical resolution, obtaining the velocity function, and yielding clear depth images, respectively. For the helicopter-borne GPR, migration velocities were obtained as root-mean-squared velocities in a two-layer model of air and ice. The radar sections show rugged subglacial topography, englacial sliding surfaces, and localised scattering noise. The maximum depth to the basement is over 79m in the subglacial valley adjacent to the south-eastern slope of the divide ridge between Fourcade and Moczydlowski Glaciers. In the ground-towed profile, we interpret a complicated conduit above possible basal water and other isolated cavities, which are a few metres wide. Near the terminus, the GPR profiles image sliding surfaces, fractures, and faults that will contribute to the tidewater calving mechanism forming icebergs in Potter Cove.

• The Journal of Engineering Geology
• /
• v.24 no.2
• /
• pp.261-271
• /
• 2014

An almost permanent anchor (friction type) is resistant to ground deformation due to the friction between the soil and grout at a fixed length from the anchor body. The purpose of this study is to calculate the force of bearing resistance for a bearing anchor in enlarged boreholes. We conducted analytical and numerical analyses, along with laboratory testing, to find the quantities of bearing resistance prior to grouting in EBA (Enlarged Bearing Anchor) construction. The force of bearing resistance from the analytical method was defined as a function of general borehole diameter, expanded borehole diameter, and soil unconfined compressive strength. We also employed the Flac 3D finite difference numerical modeling code to analyze the bearing resistance of the soil conditions. We then created a laboratory experimental model to measure bearing resistance and carried out a pull-out test. The results of these three analyses are presented here, and a regression analysis was performed between bearing resistance and uniaxial compression strength. The laboratory results yield the strongest bearing resistance, with reinforcement 28.5 times greater than the uniaxial compression strength; the analytical and numerical analyses yielded values of 13.3 and 9.9, respectively. This results means that bearing resistance of laboratory test appears to be affected by skin friction resistance. To improve the reliability of these results, a comparison field study is needed to verify which results (analytical, numerical, or laboratory) best represent field observations.