• Tunnel and Underground Space
• /
• v.13 no.6
• /
• pp.465-475
• /
• 2003

Groundwater inflow into the caverns constructed in fractured rock mass was simulated by numerical modeling, NAPSAC (DFN, discrete fracture network model) and NAMMU (CPM, continuous porous media model), a finite-element software package for groundwater flow in 3D fractured media developed by AEA Technology, UK. The input parameters for modeling were determined on surface fracture survey, core logging and single hole hydraulic test data. In order to predict the groundwater inflow more accurately, the anisotropic hydraulic conductivity was considered. The anisotropic hydraulic conductivities were calculated from the fracture network properties. With a minor adjustment during model calibration, the numerical modeling is able to reproduce reasonably groundwater inflows into cavern and the travel length and times to the ground surface along the flow paths in the normal, dry and rainy seasons.

• Journal of Korean Society of Steel Construction
• /
• v.24 no.1
• /
• pp.47-58
• /
• 2012

Dynamic behaviors of a corrugated steel plate tunnel lining system are examined under wind loads due to passing vehicles. Applied wind loads are simulated by applying the time functions as a vehicle moves through the tunnel. Wind loads are described by the pressure and suction as a vehicle arrives and leaves target positions in the tunnel. The tunnel lining is modeled using the simplified shell elements that retain the characteristics of the corrugated shapes. The displacements of the tunnel lining are evaluated under various conditions regarding wind velocity and the passing vehicles. The responses are found to increase as the vehicle velocity and wind velocity increase. A maximum displacement of 25mm occurs when two vehicles are crossing at the speed of 120km/h. A row of vehicles running consecutively minimally affects the dynamic responses with less than 2.5% of the dynamic responses enlarged and attributed to one running vehicle. It should be noted that the dynamic responses of the tunnel lining should be considered when there is no shotcrete applied.

• Tunnel and Underground Space
• /
• v.12 no.3
• /
• pp.220-228
• /
• 2002

Thermal flow through jointed rock mass was analyzed by numerical methods. The effect of a single set of joints on the heat conduction was analyzed by one-dimensional model and compared with the analytical solution. When a joint is completely dry, the joint behaves as a thermal break inducing jumps in temperature distribution even at steady state. Therefore when joints are completely dry, individual joint has to be taken into consideration to get a good result. When joints are partially or fully saturated, the thermal conductivity of the joints increases drastically and the jumps in temperature distribution become less severe. Therefore the effect of joint in heat conduction can be well absorbed by continuum anisotropic model whose thermal properties represent overall thermal properties of the intact part and the discontinuities. Since the effect of joints becomes less important as the degree of the saturation increases, the overall thermal response of the rock mass also becomes close to isotropic. Therefore it can be concluded that a great effort has to be made to obtain a precise in-situ thermal properties in order to get a good prediction of the thermal response of a jointed rock mass.

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

• Journal of the Korean Geotechnical Society
• /
• v.29 no.12
• /
• pp.35-44
• /
• 2013

In this study, a particle method without using grid was applied for analysing large deformation problems in soil flows instead of using ordinary finite element or finite difference methods. In the particle method, a continuum equation was discretized by various particle interaction models corresponding to differential operators such as gradient, divergence, and Laplacian. Soil behavior changes from solid to liquid state with increasing water content or external load. The Mohr-Coulomb failure criterion was incorporated into the particle method to analyze such three-dimensional soil behavior. The yielding and hardening behavior of soil before failure was analyzed by treating soil as a viscous liquid. First of all, a sand column test without confining pressure and strength was carried out and then a self-standing clay column test with cohesion was carried out. Large deformation from such column tests due to soil yielding or failure was used for verifying the developed particle method. The developed particle method was able to simulate the three-dimensional plastic deformation of soils due to yielding before failure and calculate the variation of normal and shear stresses both in sand and clay columns.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.6 no.3
• /
• pp.201-209
• /
• 2002

Crib wall system is one of segmental crib type wall. Crib walls are constructed from separate members with no bonds between them other than frictional. The wall units are divide into two main types termed headers and stretchers. The headers run from the front to the back of the wall, perpendicular to the wall face. The cells are created by forming a grid by stacking individual wall components known as headers and stretchers. The body of wall consists of a system of open cell which are filled with a granular material. The design of crib retaining wall is usually based on conventional design methods derived from Rankine and Coulomb theory so that is able to resist the thrust of soil behind it, because it may be assumed that the wall acts as a rigid body. However, deformation characteristics of crib walls cannot be assumed as monolithic. They consist of individual members which have been stacked to creat a three dimensional grid. Therefore, the segmental grid allows relative movement between the individual member within the wall. The three dimensional flexible grid leads to stress distribution by interaction behavior between soil and crib wall. Therefore, in this study, in order to analysis the trends of deflection of crib wall system, new numerical models based on the results of Brandl's full scale test are introduced for design concept.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.3
• /
• pp.255-263
• /
• 2017

This paper presents a beam element capable of conducting material and geometric nonlinear analysis for applications requiring the ultimate behavioral analysis of structures with composite cross-sections. The element formulation is based on co-rotational kinematics to simulate geometrically nonlinear behaviors, and it uses the fiber section method to calculate the stiffness and internal forces of the element. The proposed element was implemented using an in-house numerical program in which an arc-length method was adopted to trace severe nonlinear responses(such as snap-through or snapback), as well as ductile behavior after the peak load. To verify the proposed method of element formulation and the accuracy of the program that was used to employ the element, several numerical studies were conducted and the results from these numerical models were compared with those of three-dimensional continuum models and previous studies, to demonstrate the accuracy and computational efficiency of the element. Additionally, by evaluating an example case of a frame structure with a composite member, the effects of differences between composite material properties such as the elastic modulus ratio and strength ratio were analyzed. It was found that increasing the elastic modulus of the external layer of a composite cross-section caused quasi-brittle behavior, while similar responses of the composite structure to those of homogeneous and linear materials were shown to increase the yield strength of the external layer.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.9
• /
• pp.51-62
• /
• 2016

Parametric studies for various site conditions by using 3d numerical model were carried out in order to estimate dynamic behavior of soil-pile-structure system in dry soil deposits. Proposed model was analyzed in time domain using FLAC3D which is commercial finite difference code to properly simulate nonlinear response of soil under strong earthquake. Mohr-Coulomb criterion 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 was used as boundary condition to reduce analysis time. Also, initial shear modulus and yield depth were appropriately determined for accurate simulation of system's nonlinear behavior. Parametric study was performed by varying weight of superstructure, pile length, pile head fixity, soil relative density with proposed numerical model. From the results of parametric study, it is identified that inertial force induced by superstructure is dominant on dynamic behavior of soil-pile-structure system and effect of kinematic force induced by soil movement was relatively small. Difference in dynamic behavior according to the pile length and pile head fixity was also numerically investigated.

• Tunnel and Underground Space
• /
• v.18 no.1
• /
• pp.44-57
• /
• 2008

In this study, a 2D finite-element analysis, using the SEEP/W program, was carried out to estimate the amount of groundwater flawing into a tunnel, as well as the groundwater tables around wetland areas during and after a tunnel excavation through rock mass. Four sites along the Wonhyo-tunnel in Cheonseong Mountain (Gyeongnam, Korea) were analysed, where the model damain of the tunnel included both wetland and fault zone. The anisotropy of the hydraulic conductivities of the rock mass was calculated using the DFN model, and then used as an input parameter for the cantinuum model. Parametric study on the influencing factors was perofrmed to minimize uncertainties in the hydraulic properties. Moreover, the volumetric water content and hydraulic conductivity functions were applied ta the model to reflect the ability of a medium ta store and transport water under both saturated and unsaturated conditions. The conductivity of fault zone was assumed ta be $10^{-5}m/sec\;or\;10^{-6}m/sec$ and the conductivity of grouting zone was assumed as 1/10, 1/50 or 1/100 of the conductivity of rock mass. Totally $6{\sim}8$ cases of transient flow simulation were peformed at each site. The hydraulic conductivities of fault zone showed a significant influence on groundwater inflow when the fault zone crossed the tunnel. Also, groundwater table around wetland maintained in case that the hydraulic conductivity of grouting zone was reduced ta be less than 1/50 of the hydraulic conductivity of rock mass.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.4
• /
• pp.227-233
• /
• 1992

To study the large displacement and large rotation problems, geometrically nonlinear formulation of eccentric degenerated beam element has been developed, where the restrictions of infinitesimal rotation increments are removed and the incremental equations are derived using the Taylor series expansion of the displacement function at time t+dt. The geometrically nonlinear analyses are carried out for the cases of cantilever, square frame, shallow arch and 45-degree bend beam and all of them are compared with each of the other results published. The element developed in the present research can be efficiently utilized for analysis of the nonlinear behaviours of structures when displacements and rotations are large.