• Journal of the Korean Geotechnical Society
• /
• v.32 no.4
• /
• pp.5-14
• /
• 2016

Numerical simulation of dynamic soil-pile-structure interaction embedded in a dry sand was carried out. 3D model of the dynamic centrifuge model tests was formulated in a time domain to consider nonlinear behavior of soil using the finite difference method program, FLAC3D. As a modeling methodology, Mohr-Coulomb criteria was adopted as soil constitutive model. Soil nonlinearity was considered by adopting the hysteretic damping model, and an interface model which can simulate separation and slip between soil and pile was adopted. Simplified continuum modeling (Kim et al., 2012) was used as boundary condition to reduce analysis time. Calibration process for numerical modeling results and test results was performed through the parametric study. Verification process was then performed by comparing numerical modeling results with another test results. Based on the calibration and validation procedure, it is identified that proposed modeling method can properly simulate dynamic behavior of soil-pile system in dry condition.

• Geomechanics and Engineering
• /
• v.11 no.5
• /
• pp.625-638
• /
• 2016

This paper presents field study and numerical modeling results for a single-cell low-fill concrete box culvert under a flexible pavement subjected to traffic loading. The culvert in the field test was instrumented with displacement transducers to capture the deformations resulting from different combinations of static and traffic loads. A low-boy truck with a known axle configuration and loads was used to apply seven static load combinations and traffic loads at different speeds. Deflections under the culvert roof were measured during loading. Soil and pavement samples were obtained by drilling operation on the test site. The properties of the soil and pavement layers were determined in the laboratory. A 3-D numerical model of the culvert was developed using a finite difference program FLAC3D. Linear elastic models were used for the pavement layers and soil. The numerical results with the material properties determined in the laboratory were compared with the field test results. The observed deflections in the field test were generally smaller under moving loads than static loads. The maximum deflections measured during the static and traffic loads were 0.6 mm and 0.41 mm respectively. The deflections computed by the numerical method were in good agreement with those observed in the field test. The deflection profiles obtained from the field test and the numerical simulation suggest that the traffic load acted more like a concentrated load distributed over a limited area on the culvert. Elastic models for culverts, pavement layers, and surrounding soil are appropriate for numerical modeling of box culverts under loading for load rating purposes.

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

As Simpson et al. (1990) emphasized the importance of the straining process in the stratification and mixing in the estuarine circulation process, various researches have investigated on the relative contribution of each process to the overall potential energy anomaly dynamics. However, many numerical works have done only for two dimensional modeling along channel or the short distance cross sectional three dimensional simulations as Burchard et al. (2008) and the estuarine channel was not simulated so far. But, in the study on the physics of shallow coastal seas, spatial dimension in the three dimensional way affects significantly on results of a particular numerical model. Therefore, the comparison of two and three dimensional models is important to understand the real physics of mixing and stratification in an estuary. Also, as Geyer and MacCready (2013) pointed out that the lateral process seems to be important in determining the periodic stratifications, to study such process the three dimensional modeling must be required. The present study uses a numerical model to show the signification roles of each term of the time-dependent dynamic equation for the potential energy anomaly (PEA) in controlling along and lateral channel flows and different stratification structures. Moreover, we present the relationships between the ${\Phi}$-advection, the depth mean straining, vertical mixing and vertical advection can explain well how water level, salinity distribution and across velocity 2D model are slightly different from 3D.

• Journal of the Korean GEO-environmental Society
• /
• v.13 no.10
• /
• pp.55-62
• /
• 2012

In this research, the algorithm for simulating specific grain size distribution(GSD) with large diameter granular material was developed using the distinct element analysis program $PFC^{3D}$(Particle Flow Code). This modeling approach can generate the initial distinct elements without clump logic or cluster logic and prevent distinct element from escaping through the confining walls during the process. Finally the proposed distinct element model is used to simulate large triaxial compression test of the rockfill material and we compared the simulation output with lab test results. Simulation results of Assembly showed very well agreement with the GSD of the test sample and numerical modeling of granular material would be possible for various stress conditions using this application through the calibration.

• Structural Engineering and Mechanics
• /
• v.57 no.6
• /
• pp.1125-1142
• /
• 2016

Tunneling is one of the challenging tasks in civil engineering because it involves a variety of decision making and engineering judgment based on knowledge and experience. One of the challenges is to construct tunnels in risky areas under shallow overburden. In order to prevent the collapse of ceilings and walls of a large tunnels, in such conditions, either a sequential excavation method (SEM) or ground reinforcing method, or a combination of both, can be utilized. This research deals with the numerical modeling of L-profiles and pipe fore-poling pre-support systems in the adit tunnel in northwestern Iran. The first part of the adit tunnel has been drilled in alluvial material with very weak geotechnical parameters. Despite applying an SEM in constructing this tunnel, analyzing the results of numerical modeling done using FLAC3D, as well as observations during drilling, indicate the tunnel instability. To improve operational safety and to prevent collapse, pre-support systems, including pipe fore-poling and L-profiles were designed and implemented. The results of the numerical modeling coupled with monitoring during operation, as well as the results of instrumentation, indicate the efficacy of both these methods in tunnel collapse prevention. Moreover, the results of modeling using FLAC3D and SECTION BUILDER suggest a double angle with equal legs ($2L100{\times}100{\times}10mm$) in both box profile and tee array as an alternative section to pipe fore-poling system while neither $L80{\times}80{\times}8mm$ nor $2L80{\times}80{\times}8mm$ can sustain the axial and shear stresses exerted on pipe fore-poling system.

• The Journal of the Korea Contents Association
• /
• v.16 no.2
• /
• pp.558-565
• /
• 2016

Recent 3D printing technology has been grafting onto various medical practices. In light of this trend, this research is intended to examine the figuration surface's accuracy of 3D images made by using DICOM images after printing by 3D printing. The medical images were obtained from animal bone objects, while the objects were printed after undergoing STL file conversion for 3D printing purposes. Ultimately, after the 3D figuration, which was obtained by the original animal bones and 3D printing, was scanned by 3D scanner, 3D modeling was merged each other and the differences were compared. The result analysis was conducted by visual figuration comparison, color comparison of modeling's scale value, and numerical figuration comparison. The shape surface was not visually distinguished; the numerical figuration comparison was made from the values measured from the four different points on the X, Y and Z coordinates. The shape surface of the merged modeling was smaller than the original object (the animal bone) by average of -0.49 mm in the 3D printed figuration. However, not all of the shape surface was uniformly reduced in size and the differences was within range of -0.83 mm on the experiment.

• Computers and Concrete
• /
• v.22 no.1
• /
• pp.1-10
• /
• 2018

In this study an implicit algorithm for modeling of propagation of macrocracks in 3D concrete structures suffering from alkali-silica reaction has been developed and implemented. The formulation of the problem prior to the onset of localized deformation is based on a chemo-elasticity approach. The localized deformation mode, involving the formation of macrocracks, is described using a simplified form of the strong discontinuity approach (SDA) that employs a volume averaging technique enhanced by a numerical procedure for tracing the propagation path in 3D space. The latter incorporates a non-local smoothening algorithm. The formulation is illustrated by a number of numerical examples that examine the crack propagation pattern in both plain and reinforced concrete under different loading scenarios.

• 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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.573-578
• /
• 2008

Finite element analysis of impact response of a motorcycle helmet is presented in this paper. The finite element LS-DYNA3D code is used to simulate the impact response of the helmet including of plastic shell, foam liner, and magnesium headform. Since the maximum accelerations at center of gravity of the headform obtained by numerical analysis and experiment agree well, the numerical simulation is proved to be valid.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.451-456
• /
• 2007

A finite element analysis and experiment study of a motorcycle helmet are presented in this paper. The finite element LS-DYNA3D code is used to analyze the helmet. The test specimen, instruments, and setup procedures are described. Since the displacements and Von-Mises stresses obtained by numerical analysis and experiment agree well, the numerical simulation is proved to be valid.