• Proceedings of the Korean Geotechical Society Conference
• /
• 2000.11b
• /
• pp.213-220
• /
• 2000

Numerical modeling of cut slope has some limits in simulating the real slopes. In the case of 2D analysis of slope stability, it is assumed that slope is simply straight even when it is concave or convex in plan view. In this study, 3D analysis in curved shape slopes has been conducted for the comparison with 2D analysis in terms of failure mode and factor of safety. For this, 3D analysis by FLAC3D was compared with 2D analysis in plane strain condition and axi-symmetric model condition by FLAC. It was also observed how safety factors of slopes were affected by the variation of the tensile strength and cohesion, which are important variables to decide whether the slope fails or not. 2D analysis of concave slopes under plane strain condition showed much smaller safety factors by 16-40 % errors depending on the radius of curvature of slopes, compared to the more realistic values from 3D analysis. In case of convex slopes, the lower values by 7-10 % has been reported. 2D analysis of axi-symmetric model showed also smaller safety factors by 6-10 % and by 2-4 %, in case of concave and convex slopes, respectively. Such results are expected to contribute to the better understanding of failure process and could be applied for improved design of slopes.

• Journal of Coastal Disaster Prevention
• /
• v.1 no.1
• /
• pp.1-9
• /
• 2014

In order to examine the characteristics of tsunami run-up heights on impermeable/permeable slope, a numerical wave tank by upgrading LES-WASS-3D was used in this study. Then, the model were compared with existing hydraulic model test for its verification. The numerical results well reproduced experimental results of solitary wave deformation, propagation and run-up height under various conditions. Also, the numerical simulation with a slope boundary condition has been carried out to understand solitary wave run-up on impermeable/permeable slope. It is shown that the run-up heights on permeable slope is 52.64-63.2% smaller than those on the impermeable slope because of wave energy dissipation inside the porous media. In addition, it is revealed that the numerical results with slope boundary condition agreed well with experimental results in comparison with the results by using stair type boundary condition.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.245-245
• /
• 2015

In the recent years, due to long-lasting heavy rainfall events, a large number of landslides have been observed in the mountainous area of the world. Such landslides can also form a dam as it blocks the course of a river, which may burst and cause a catastrophic flood. Numerical analysis of landslide dam formation is rarely available, while laboratory experimental studies often use assumed shape to analyze the landslide dam failure and flood hydraulics in downstream. In this study, both experimental and numerical studies have been carried out to investigate the formation of landslide dam. Two case laboratory experiments were conducted in two flumes simultaneously. The first flume (2.0 m 0.6 m 0.5 m) was set at $22^{\circ}$ and $27^{\circ}$ slope to generate the landslide using rainfall intensity of 70.0 mm/hr. On the other hand, the second flume (1.5 m 0.25 m 0.3 m) was set perpendicularly at the downstream end of the first flume to receive the landslide mass forming landslide dam. The formation of landslide dam was observed at $15^{\circ}$ slope of the second flume. The whole processes including the landslide initiation and movement of the landslide mass into the second channel was captured by three digital cameras. In numerical analysis, a two-dimensional (2D) seepage flow model, a 2D slope stability model (Spencer method) and a 2D landslide dam-geometry evaluation model were coupled as a single unit. This developed model can determine the landslide occurrence time, the failure mass and the geometry of landslide dam deposited in the second channel. The data obtained from numerical simulation results has good agreement with the experimental measurements.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.18 no.5
• /
• pp.388-397
• /
• 2012

In this study, 3D-NIT(3-Dimensional Numerical Irregular wave Tank) model in which regular wave as well as stable irregular wave can be generated in 3-dimensional numerical irregular wave tank was proposed. To verify validity, the following steps need to be conducted: 1) comparative analysis between calculated waveforms and targeted waveforms at the wave generating point, 2) comparative analysis with the existing experimental values of overtopping volume estimated, targeting shore protection structures installed on a slope bed, 3) comparison with the existing numerical and hydraulic experimental results through application in the analysis on the wave deformation by structures and wave force acting on the vertical cylindrical structures. Based on the results, characteristics of the breaking wave forces according to incident waves and interval distance of structures were identified through application of 3D-NIT model in the analysis on the breaking wave forces acting on the cylindrical structures installed on a slope bed, and reflection and overtopping was reviewed through application in the special breakwaters on the domestic fields. The numerical results obtained the 3D-NIT model are in good agreement with experimental results, and its applicaion to the complex-shpaed coastal structures is verified.

• Journal of Ocean Engineering and Technology
• /
• v.28 no.4
• /
• pp.331-337
• /
• 2014

We analyzed the 3-D characteristics of the dynamic response of seabed around a submerged breakwater due to wave loading using a 3-D numerical scheme (LES-WASS-3D). Using our model, which considers the wave-structure-sandy seabed interactions in a 3-D wave field, we were able to investigate the 3-D characteristics of the pore-water pressure in the seabed around the submerged breakwater under various incident wave conditions. To verify the 3-D numerical analysis method suggested in this study, we compared the numerical results with the existing experimental results and found good agreement between them. The numerical analysis reveals that high pore-water pressure in the seabed is generated below a large wave height at the front slope of the submerged breakwater. It was also shown that the non-dimensional pore-water pressure in the seabed increases as the wave period increases because the wave energy dissipation decreases on the submerged breakwater and seabed as the wave period increases.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.4
• /
• pp.383-389
• /
• 2006

In three dimensional wave field, a direct numerical simulation model, which is able to handle free-fall and slope-fall of a waterdrop due to overtopping, is proposed to evaluate the overtopping rate on the coastal structures like an absorbing revetment. A comparison between the numerical model and existing experimental results for overtopping rate was made to validate the proposed numerical model's accuracy, and showed fairly good agreement between them. It is confirmed in numerical and hydraulic(existing) model test that the overtopping quantity on a absorbing revetment becomes larger with an increase in Ursell number. Also, the overtopping rate estimated by 3-D numerical model is compared with it obtained by 2-D numerical model.

• Journal of Korea Water Resources Association
• /
• v.54 no.11
• /
• pp.943-953
• /
• 2021

The propagation of impact wave induced by landslide and debris flow occurred on the slope of lake, reservoir and bays is a three-dimensional natural phenomenon associated with strong interaction of debris flow and water flow in complex geometrical environments. We carried out 3D numerical modeling of such impact wave in a bay using a multiphase turbulence flow model and a rheology model for non-Newtonian debris flow. Numerical results are compared with previous experimental result to evaluate the performance of present numerical approach. The results underscore that the reasonable predictions of both thickness and speed of debris flow head penetrating below the water surface are crucial to accurately reproduce the maximum peak height and free surface profiles of impact wave. Two predictions computed using different initial debris flow thicknesses become different from the instant when the peaks of impact waves fall due to the gravity. Numerical modeling using relatively thick initial debris flow thickness appears to well reproduce the water surface profile of impact wave propagating across the bay as well as wave run-up on the opposite slope. The results show that the maximum run-up height on the opposite slope is not sensitive to the initial thickness of debris flows of same total volume. Meanwhile, appropriate rheology model for debris flow consisting of inviscid particle only should be employed to more accurately reproduce the debris flow propagating along the channel bottom.

• Economic and Environmental Geology
• /
• v.33 no.6
• /
• pp.557-563
• /
• 2000

Subsidence hazard has never been considered seriously until recent yews in Korea, although its socioeconomic impact on Korea becomes more and more enormous. There have been a few studies for the application of GIS analysis technique to the prediction of subsidence hazard. For GIS analysis, several factors, which are represented by coverage, are considered and selected for building a GIS model. Numerical method was used to quantify the importance of each factor in GIS model and the result from numerical modeling using FLAC was compared with that from previous research based on empirical methods. Analysis in 3-D needs more computer resources (i.e. memory). Therefore that in 2.5-D was considered to overcome the problem. Not only maximum vertical subsidence but also maximum horizontal strain and maximum slope have been considered for the assessment of subsidence hazard. The model can be easily modified for the purpose of applications in any subsidence area, especially cavern or abandoned mines under thick soil layer.

• Geomechanics and Engineering
• /
• v.12 no.6
• /
• pp.1021-1046
• /
• 2017

The assessment of slope stability is an essential task in geotechnical engineering. In this paper, a three-dimensional (3D) finite element analysis (FEA) was employed to investigate the performance of different shear pin arrangements to increase the stability of a soil block resting on an inclined plane with a low-interface friction plane. In the numerical models, the soil block was modeled by volume elements with linear elastic perfectly plastic material in a drained condition, while the shear pins were modeled by volume elements with linear elastic material. Interface elements were used along the bedding plane (bedding interface element) and around the shear pins (shear pin interface element) to simulate the soil-structure interaction. Bedding interface elements were used to capture the shear sliding of the soil on the low-interface friction plane while shear pin interface elements were used to model the shear bonding of the soil around the pins. A failure analysis was performed by means of the gravity loading method. The results of the 3D FEA with the numerical models were compared to those with the physical models for all cases. The effects of the number of shear pins, the shear pin locations, the different shear pin arrangements, the thickness and the width of the soil block and the associated failure mechanisms were discussed.

• Journal of the Korean Geotechnical Society
• /
• v.38 no.7
• /
• pp.5-24
• /
• 2022

The current performance evaluation of slope anchors qualitatively determines the physical bonding between the anchor head and ground as well as cracks or breakage of the anchor head. However, such performance evaluation does not measure these primary factors quantitatively. Therefore, the time-dependent management of the anchors is almost impossible. This study is an evaluation of the 3D numerical model by SfM which combines UAS images with terrestrial LiDAR to collect numerical data on the damage factors. It also utilizes the data for the quantitative maintenance of the anchor system once it is installed on slopes. The UAS 3D model, which often shows relatively low precision in the z-coordinate for vertical objects such as slopes, is combined with terrestrial LiDAR scan data to improve the accuracy of the z-coordinate measurement. After validating the system, a field test is conducted with ten anchors installed on a slope with arbitrarily damaged heads. The damages (such as cracks, breakages, and rotational displacements) are detected and numerically evaluated through the orthogonal projection of the measurement system. The results show that the introduced system at the resolution of 8K can detect cracks less than 0.3 mm in any aperture with an error range of 0.05 mm. Also, the system can successfully detect the volume of the damaged part, showing that the maximum damage area of the anchor head was within 3% of the original design guideline. Originally, the ground adhesion to the anchor head, where the z-coordinate is highly relevant, was almost impossible to measure with the UAS 3D numerical model alone because of its blind spots. However, by applying the combined system, elevation differences between the anchor bottom and the irregular ground surface was identified so that the average value at 20 various locations was calculated for the ground adhesion. Additionally, rotation angle and displacement of the anchor head less than 1" were detected. From the observations, the validity of the 3D numerical model can obtain quantitative data on anchor damage. Such data collection can potentially create a database that could be used as a fundamental resource for quantitative anchor damage evaluation in the future.