• Tunnel and Underground Space
• /
• v.33 no.3
• /
• pp.189-207
• /
• 2023

In the present study, the thermoshearing experiment on a rough rock fracture were modeled using a three-dimensional grain-based distinct element model (GBDEM). The experiment was conducted by the Korea Institute of Construction Technology to investigate the progressive shear failure of fracture under the influence of thermal stress in a critical stress state. The numerical model employs an assembly of multiple polyhedral grains and their interfaces to represent the rock sample, and calculates the coupled thermo-mechanical behavior of the grains (blocks) and the interfaces (contacts) using 3DEC, a DEM code. The primary focus was on simulating the temperature evolution, generation of thermal stress, and shear and normal displacements of the fracture. Two fracture models, namely the mated fracture model and the unmated fracture model, were constructed based on the degree of surface matedness, and their respective behaviors were compared and analyzed. By leveraging the advantage of the DEM, the contact area between the fracture surfaces was continuously monitored during the simulation, enabling an examination of its influence on shear behavior. The numerical results demonstrated distinct differences depending on the degree of the surface matedness at the initial stage. In the mated fracture model, where the surfaces were in almost full contact, the characteristic stages of peak stress and residual stress commonly observed in shear behavior of natural rock joints were reasonably replicated, despite exhibiting discrepancies with the experimental results. The analysis of contact area variation over time confirmed that our numerical model effectively simulated the abrupt normal dilation and shear slip, stress softening phenomenon, and transition to the residual state that occur during the peak stress stage. The unmated fracture model, which closely resembled the experimental specimen, showed qualitative agreement with the experimental observations, including heat transfer characteristics, the progressive shear failure process induced by heating, and the increase in thermal stress. However, there were some mismatches between the numerical and experimental results regarding the onset of fracture slip and the magnitudes of fracture stress and displacement. This research was conducted as part of DECOVALEX-2023 Task G, and we expect the numerical model to be enhanced through continued collaboration with other research teams and validated in further studies.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.151-164
• /
• 1998

Changes in tides and tidal currents in Chonsu Bay caused by the construction of Seosan A and B reclamation regions have been investigated using a depth-integrated two-dimensional tidal model. Three water level recorders were deployed for about one month and used for the specification of the open boundary condition. The computed currents were found to be in a good agreement with the measurements at two stations within Chonsu Bay. Comparison of the absolute velocities computed with the conditions before and after the embankment clearly shows that the reduction in tidal current amplitude is evident throughout the bay and the magnitude of the reduction increases to the north. Calculation shows that the embankment has advanced the time of drying at the northern part of the bay by about 51 minutes and has increased the exposure time by about 23 minutes. The high water time has advanced by about 1 hour, lowering the high water level by about 15 cm.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.5
• /
• pp.323-331
• /
• 2019

In this study, a computational model approach for the modeling of hydrodynamic pressures acting on radial gates during strong earthquakes is proposed. The use of the dynamic layering method with the Arbitrary Lagrangian Eulerian (ALE) algorithm and the SIMPLE method for simulating free reservoir surface flow in addition to moving boundary interfaces between the fluid domain and a structure due to earthquake excitation are suggested. The verification and validation of the proposed approach are realized by comparisons performed using the renowned formulation derived by the experimental results for vertical and inclined dam surfaces subjected to earthquake excitation. A parameter study for the truncated lengths of the two-dimensional fluid domain demonstrates that twice the water level leads to efficient and converged computational results. Finally, numerical simulations for large radial gates with different curvatures subjected to two strong earthquakes are successfully performed using the suggested computational model.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.136-147
• /
• 2004

In this study, the characteristics of magnetotelluric (MT) responses due to a three-dimensional (3-D) body are analyzed with 3-D numerical modeling. The first model for the analysis consists of a single isolated conductive body embedded in a resistive homogeneous half-space. The second model has an additional conductive overburden while the other conditions remain the same as the first one. The analysis of apparent resistivities shows well that the 3-D effects are dominant over some frequency range for the first model. Two mechanisms, current channeling and induction, for secondary electric fields due to the conductive body are analyzed at various frequencies: at high frequencies induction is more dominant than channeling, while at low frequencies channeling is more dominant than induction. Tippers have a strong relation to the position of anomalous body and the real and imaginary parts of induction vector also indicate the position of anomalous body. off-line conductive anomaly sometimes causes severe problem in 2-D interpretation. In such case, induction vector analysis can give information on the existence and location of the anomalous body. Each parameter of the second model shows similar responses as those of the first model. The only difference is that the magnitude of all parameters is decreased and that the domain showing the 3-D effects becomes narrower. As shown in this study, the analysis of 3-D effects provides a useful and effective means to understand the 3-D subsurface structure and to interpret MT survey data.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.12 no.4
• /
• pp.267-274
• /
• 2014

The first Low- and Intermediate-Level Waste (LILW) disposal facility with 6 silos has been constructed in granite host rock saturated with groundwater in Korea. A two-dimensional numerical modeling on gas migration was carried out using TOUGH2 with EOS5 module in the disposal facility. Laboratory-scale experiments were also performed to measure the important properties of silo concrete related with gas migration. The gas entry pressure and relative gas permeability of the concrete was determined to be $0.97{\pm}0.15bar$ and $2.44{\times}10^{-17}m^2$, respectively. The results of the numerical modeling showed that hydrogen gas generated from radioactive wastes was dissolved in groundwater and migrated to biosphere as an aqueous phase. Only a small portion of hydrogen appeared as a gas phase after 1,000 years of gas generation. The results strongly suggested that hydrogen gas does not accumulate inside the disposal facility as a gas phase. Therefore, it is expected that there would be no harmful effects on the integrity of the silo concrete due to gas generation.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.11 no.3
• /
• pp.179-196
• /
• 1986

In this paper, the peformance of a statistical packet voice/data multiplexer is studied. In ths study we assume that in the packet voice/data multiplexer two separate finite queues are used for voice and data traffics, and that voice traffic gets priority over data. For the performance analysis we divide the output link of the multiplexer into a sequence of time slots. The voice signal is modeled as an (M+1) - state Markov process, M being the packet generation period in slots. As for the data traffic, it is modeled by a simple Poisson process. In our discrete time domain analysis, the queueing behavior of voice traffic is little affected by the data traffic since voice signal has priority over data. Therefore, we first analyze the queueing behavior of voice traffic, and then using the result, we study the queueing behavior of data traffic. For the packet voice multiplexer, both inpur state and voice buffer occupancy are formulated by a two-dimensional Markov chain. For the integrated voice/data multiplexer we use a three-dimensional Markov chain that represents the input voice state and the buffer occupancies of voice and data. With these models, the numerical results for the performance have been obtained by the Gauss-Seidel iteration method. The analytical results have been verified by computer simylation. From the results we have found that there exist tradeoffs among the number of voice users, output link capacity, voic queue size and overflow probability for the voice traffic, and also exist tradeoffs among traffic load, data queue size and oveflow probability for the data traffic. Also, there exists a tradeoff between the performance of voice and data traffics for given inpur traffics and link capacity. In addition, it has been found that the average queueing delay of data traffic is longer than the maximum buffer size, when the gain of time assignment speech interpolation(TASI) is more than two and the number of voice users is small.

• Journal of the Korea Concrete Institute
• /
• v.17 no.2 s.86
• /
• pp.281-290
• /
• 2005

This paper presents experimental and analytical studies on long-term behavior of square CFT columns under central axial loading. Two loading cases are considered; (1) the load applied only at the inner concrete of the column and (2) the load applied simultaneously on both the concrete and the steel tube. Four specimens of square CFT columns were tested under the two loading cases, and basic creep test for two concrete specimens was performed to find out the creep properties of the inner concrete. Three-dimensional finite element analysis models were established and verified with the experimental results. The verification shows that the prediction for the long-term behavior of actual square CFT columns is possible from the three dimensional finite element modeling considering the bond behavior between steel tube and inner concrete. Also, experimental results and numerical calculations revealed that the bond stress Induced by the confinement pressure as well as the slip between inner concrete and steel tube were increased with time In the first loading case. However, the confinement by the loading Plate was decreased with time while increasing confinement effect by the steel tube was observed over time. In contrast no confinement effects occur in the second loading case.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.12
• /
• pp.659-667
• /
• 2017

This study analyzes the changes in the flow characteristics due to the difference in inflow discharges from the main channel and tributary at the confluence of the Nakdong and Geumho Rivers. The analysis was done using a two-dimensional numerical method. The study site has complicated flow patterns because of the discharge variation from the main stream and tributary. The study section has a meandering main channel, and the hydraulic characteristics cannot be defined with simple conditions such as the confluence angle of the channels or the ratio of the channel widths. An actual flood event in 2012 was applied in the numerical simulation. The maximum velocity occurred in the meandering section after passing the confluence, where a rapid change was expected. A high velocity and large bed change in this section were observed in the simulation results. The variation of discharges from the main channel and tributary was a more dominant factor in the flow and bed changes for the normal flow conditions than the flood event. This indicates that countermeasures for channel stabilization should be considered in the meandering section downstream of the confluence section, and countermeasures for the study section should be investigated.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.29 no.4
• /
• pp.363-371
• /
• 2023

To construct a ship, blocks of various sizes must be moved and erected . In this process, lugs are used such that they match the block fastening method and various functions suitable for the characteristics of each shipyard facility. The sizes and shapes of the lugs vary depending on the weight and shape of the block structures. The structure is reinforced by welding the doubling pads to compensate for insufficient rigidity around the holes where the shackle is fastened. As for the method of designing lugs according to lifting loading conditions, a simple calculation based on the beam theory and structural analysis using numerical modeling are performed. In the case of the analytical method, a standardized evaluation method must be established because results may differ depending on the type of element and modeling method. The application of this ambiguous methodology may cause serious safety problems during the process of moving and turning-over blocks. In this study , the effects of various parameters are compared and analyzed through numerical structural analysis to determine the modeling conditions and evaluation method that can evaluate the actual structural response of the lug. The modeling technique that represents the plate part and weld bead around the lug hole provides the most realistic behavior results. The modeling results with the same conditions as those of the actual lug where only the weld bead is connected to the main body of the lug, showed a lower ulimated strength compared with the results obtained by applying the MPC load. The two-dimensional shell element is applied to reduce the modeling and analysis time, and a safety working load was verified to be predicted by reducing the thickness of the doubling pad by 85%. The results of the effects of various parameters reviewed in the study are expected to be used as good reference data for the lug design and safe working load prediction.

• Computers and Concrete
• /
• v.31 no.4
• /
• pp.337-348
• /
• 2023

Shear wall is commonly used as a lateral force resisting system of concrete mid-rise and high-rise buildings, but it brings challenges in providing relatively large space throughout the building height. For this reason, the structure system where the upper structure with bearing, non-bearing and/or shear walls that sits on top of a transfer plate system supported by widely spaced columns at the lower stories is preferred in some regions, particularly in low to moderate seismic regions in Asia. A thick reinforced concrete (RC) plate has often been used as a transfer system, along with RC transfer girders; however, the RC plate becomes very thick for tall buildings. Applying the post-tensioning (PT) technique to RC plates can effectively reduce the thickness and reinforcement as an economical design method. Currently, a simplified model is used for numerical modeling of PT transfer plate, which does not consider the interaction of the plate and the upper structure. To observe the actual behavior of PT transfer plate under seismic loads, it is necessary to model whole parts of the structure and tendons to precisely include the interaction and the secondary effect of PT tendons in the results. This research evaluated the seismic behavior of shear wall-type residential buildings with PT transfer plates for the condition that PT tendons are included or excluded in the modeling. Three-dimensional finite element models were developed, which includes prestressing tendon elements, and response spectrum analyses were carried out to evaluate seismic forces. Two buildings with flat-shape and L-shape plans were considered, and design forces of shear walls and transfer columns for a system with and without PT tendons were compared. The results showed that, in some cases, excluding PT tendons from the model leads to an unrealistic estimation of the demands for shear walls sit on transfer plate and transfer columns due to excluding the secondary effect of PT tendons. Based on the results, generally, the secondary effect reduces shear force demand and axial-flexural demands of transfer columns but increases the shear force demand of shear walls. The results of this study suggested that, in addition to the effect of PT on the resistance of transfer plate, it is necessary to include PT tendons in the modeling to consider its effect on force demand.