• Tunnel and Underground Space
• /
• v.30 no.1
• /
• pp.39-62
• /
• 2020

In this study, Barcelona Basic Model (BBM) was implemented into TOUGH2-MP/FLAC3D for the numerical analysis of coupled thermo-hydro-mechanical (THM) behavior of unsaturated soils and the prediction of long-term behaviors. Similar to the methodology described in a previous study for the implementation of BBM into TOUGH-FLAC, the User Defined Model (UDM) of FLAC based on the Modified Cam Clay Model (MCCM) and the FISH function of FLAC3D were used to extend the existing MCCM module in FLAC3D for the implementation of BBM into TOUGH2-MP/FLAC3D. In the developed BBM module in TOUGH2-MP/FLAC3D, the plastic strains due to change in suction increase (SI) in addition to mean effective stress are calculated. In addition to loading-collapse (LC) yield surface, suction increase (SI) yield surface is changed by hardening rules in the developed BBM module. Several numerical simulations were conducted to verify and validate the implementation of BBM: using an example presented in the FLAC3D manual for the standard MCCM, simulation results using COMSOL, and experimental data presented in SKB Reports. In addition, the developed BBM analysis module was validated by simultaneously performing a series of modeling tests that were performed for the validation of the Quick tools developed for the purpose of effectively deriving BBM parameters, and by comparing the Quick tools and Code_Bright results reported in a previous study.

• Geotechnical Engineering
• /
• v.14 no.2
• /
• pp.93-106
• /
• 1998

The behavior of axially loaded pile was analyzed by two methodologies: one is the finite difference method using load transfer curves recommended by API(1993) , and the other is the numerical analysis using the FLAC program. From both analyses, load-displacement curves and load distributions along the depth were evaluated appropriately for the measured. The analysis using the FLAC could capture the nonlinearity of load-displacement curve even for unloading and reloading cases, since the unloaded stress paths of fill layer elements occurred on the failure envelop. Futhermore, the measured load transfer curves were compared with the API recommendations and with the calculations obtained front the results of the FLAC analysis for the interpretation of the transfer behavior between the soil and the pile under axial loadings. It was concluded that the atrial behavior of open ended piles at Pusan could be evaluated by both the finite difference analysis using API load transfer curves and the numerical analysis using FLAC.

• Tunnel and Underground Space
• /
• v.15 no.5 s.58
• /
• pp.344-351
• /
• 2005

Direct back analysis software is coded based on Simplex method using FLAC as a subroutine. For the verification of the software, 12 different cases are assumed combining various displacements in different measuring locations around a tunnel. The number of displacements for cases varies from 3 to 240. It is verified that the exact elasticity and the primary stress state of rock around a tunnel could be found through iterative calculation regardless of the locations and number of displacements and initial values needed for the direct back analysis.

• Tunnel and Underground Space
• /
• v.12 no.1
• /
• pp.37-42
• /
• 2002

A back analysis algorithm was developed to determine the major parameters for tunnel design; the elastic modulus(E) and the ratio of horizontal to vertical stress(K). The algorithm is based on direct search method and was coded by FISH language of FLAC, a commercial finite difference program. Developed code was applied on some models to verify the validity and estimate the efficiency of the algorithm. Verification by theoretical solutions and published results of Gens' research, was successful.

• The Journal of Engineering Geology
• /
• v.12 no.3
• /
• pp.345-360
• /
• 2002

This study was for analysing the sedimentary rock slope stability and providing the reinforcement method that can heighten stability. The study area consists of Cretaceous basalt or basaltic tuff belonging to Hak-Bong Basalt Formation in Ha-Yang Group. Nature of geological structure confirmed in this area ars bedding, joint and fault. Majority of geological structure that affect most relationship rock slope stability is bedding. It is shown that dip direction is 120~160/15~25. In other structure, joint sets are shown that dip direction of set 1 is 310~330/65~85 and set 2 is 230~250/70~85. Joint set 3 shows above 85$^{\circ}$ high angle on NE trend although do not show clear. Stability analysis about rock slope used kinematic analysis, limit equilibrium method and FLAC by numerical analysis method. FLAC is continuum model that use Fintie Defferentce Method, but could use Interfaces Module and get discrete model's analysis effect such as UDEC.

• Tunnel and Underground Space
• /
• v.30 no.3
• /
• pp.242-255
• /
• 2020

The numerical simulations of Heater Experiment-D (HE-D) at the Mont Terri rock laboratory in Switzerland were performed to investigate an applicability of FLAC3D to reproduce the coupled thermo-hydro-mechanical (THM) behaviour in Opalinus Clay, as part of the DECOVLEX-2015 project Task B. To investigate the reliability of numerical simulations of the coupled behaviour using FLAC3D code, the simulation results were compared with the observations from the in-situ experiment, such as temperature at 16 sensors, pore pressure at 6 sensors, and strain at 22 measurement points. An anisotropic heat conduction model, fluid flow model, and transversely isotropic elastic model in FLAC3D successfully represented the coupled thermo-hydraulic behaviour in terms of evolution for temperature and pore pressure, however, performance of the models for mechanical behavior is not satisfactory compared with the measured strain.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.335-358
• /
• 2020

We present a numerical model to simulate coupled hydro-mechanical behavior of fault using TOUGH-FLAC simulator. This study aims to develop a numerical method to estimate fluid injection-induced fault reactivation in low permeability rock and to access the relevant hydro-mechanical stability in rock as part of DECOVALEX-2019 Task B. A coupled fluid flow and mechanical interface model to explicitly represent a fault was suggested and validated from the applications to benchmark simulations and the field experiment at Mont Terri underground laboratory in Switzerland. The pressure build-up, hydraulic aperture evolution, displacement, and stress responses matched those obtained at the site, which indicates the capability of the model to appropriately capture the hydro-mechanical processes in rock fault.

• Explosives and Blasting
• /
• v.41 no.3
• /
• pp.38-61
• /
• 2023

An Excavation Damaged Zone(EDZ) caused by blasting impact changes rock properties, in situ stress distribution, etc., and its effects are noticeable at around a radioactive waste repository located at deep underground. In particular, the increase in permeability due to the formation of cracks may significantly increase the amount of groundwater inflow and the possibility of radioactive nuclide outflow. In this study, FLAC2D and FLAC3D were used to analyze the mechanical and thermal behaviors for three categories: a)No EDZ, b)Uniform EDZ, and c)Random EDZ. It was found that the tunnel displacement in the Random EDZ case was 423% higher than that in the No EDZ case and was 16% higher than that in the Uniform EDZ case. Tunnel inflow in the Random EDZ was also 17.3% and 10.8% higher than that in the No EDZ and the Uniform EDZ case, respectively. The permeability around the tunnel was increased by up to 10 times in the corner of the tunnel wall and roof due to the stress redistribution after excavation. From the computer simulation, it was found that the permeability around the tunnel wall was partially increased but the overall tunnel inflow was decreased with increase of stress ratio. Mechanical analysis using FLAC 3D showed similar results. Slight difference between 2D and 3D could be explained with the development of plastic zone during the advance of tunnel excavation in 3D.

• Geomechanics and Engineering
• /
• v.8 no.3
• /
• pp.415-440
• /
• 2015

In the present case study, High Strain Dynamic Testing of piles is conducted at 3 different locations of Kolkata city of India. The raw field data acquired is analyzed using Pile Driving Analyzer (PDA) and CAPWAP (Case Pile Wave Analysis Programme) computer software and load settlement curves along with variation of force and velocity with time is obtained. A finite difference based numerical software FLAC3D has been used for simulating the field conditions by simulating similar soil-pile models for each case. The net pile displacement and ultimate pile capacity determined from the field tests and estimated by using numerical analyses are compared. It is seen that the ultimate capacity of the pile computed using FLAC3D differs from the field test results by around 9%, thereby indicating the efficiency of FLAC3D as reliable numerical software for analyzing pile foundations subjected to impact loading. Moreover, various parameters like top layers of cohesive soil varying from soft to stiff consistency, pile length, pile diameter, pile impedance and critical height of fall of the hammer have been found to influence both pile displacement and net pile capacity substantially. It may, therefore, be suggested to include the test in relevant IS code of practice.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.5
• /
• pp.163-173
• /
• 2003

In this paper, a numerical comparison of predictions by limit equilibrium analysis and 3D analysis is presented for slope/pile system. Special attention is given to the coupled analysis based on the explicit finite difference code, FLAC 3D. To this end, an internal routine (FISH) was developed to calculate a factor of safety for a file reinforced slope according to shear strength reduction technique. The case of coupled analyses was performed for stabilizing piles in slope in which the pile response and slope stability are considered simultaneously. In this study, by using these methods, the failure surfaces and factors of safety were compared and analyzed in several cases, such as toe, middle and top of the slope, respectively. Furthermore, the coupled method based on shear strength reduction technique was verified by the comparison with other analysis results.