• Geomechanics and Engineering
• /
• v.34 no.6
• /
• pp.683-696
• /
• 2023

The destruction and fracture of rock masses are crucial components in engineering and there is an increasing demand for the study of the influence of rock mass heterogeneity on the safety of engineering projects. The numerical manifold method (NMM) has a unified solution format for continuous and discontinuous problems. In most NMM studies, material homogeneity has been assumed and despite this simplification, fracture mechanics remain complex and simulations are inefficient because of the complicated topology updating operations that are needed after crack propagation. These operations become computationally expensive especially in the cases of heterogeneous materials. In this study, a heterogeneous model algorithm based on stochastic theory was developed and introduced into the NMM. A new fracture algorithm was developed to simulate the rupture zone. The algorithm was validated for the examples of the four-point shear beam and semi-circular bend. Results show that the algorithm can efficiently simulate the rupture zone of heterogeneous rock masses. Heterogeneity has a powerful effect on the macroscopic failure characteristics and uniaxial compressive strength of rock masses. The peak strength of homogeneous material (with heterogeneity or standard deviation of 0) is 2.4 times that of heterogeneous material (with heterogeneity of 11.0). Moreover, the local distribution of parameter values can affect the configuration of rupture zones in rock masses. The local distribution also influences the peak value on the stress-strain curve and the residual strength. The post-peak stress-strain curve envelope from 60 random calculations can be used as an estimate of the strength of engineering rock masses.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.4A
• /
• pp.529-536
• /
• 2008

The failure behavior of structures is changed under different loading rates, which might arise from the rate dependency of materials. This phenomenon has been focused in the engineering fields. However, the failure mechanism is not fully understood yet, so that it is hard to be implemented in numerical simulations. In this study, the numerical experiments to a brittle material are simulated by the Molecular Dynamics (MD) for understanding the rate dependent failure behavior. The material specimen with a notch is modeled for the compact tension test simulation. Lennard-Jones potential is used to describe the properties of a brittle material. Several dynamic failure features under 6 different loading rates are achieved from the numerical experiments, where remarkable characteristics such as crack roughness, crack recession/arrest, and crack branching are observed during the crack propagation. These observations are interpreted by the energy inflow-consumption rates. This study will provides insight about the dynamic failure mechanism under different loading rates. In addition, the applicability of the MD to the macroscopic mechanics is estimated by simulating the previous experimental research.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.6D
• /
• pp.779-784
• /
• 2011

Freeway FFS (Free-Flow Speed) under the state of LOS (Level of Service) A has been regarded as a homogeneous state to all levels of traffic volume. The features of low-density FFS is served as the basic data for a decision on limited maximal speed, accident analysis, simulation modeling etc. A few researches for the macroscopic characteristics of design and operation speed have been reported, and any study for the microscopic features of freeway free-flow speed under the state of low density has not been done. Therefore, the characteristics of low-density FFS according to the level of traffic volume (1-3 veh/30 sec) and daily hour periods (0-5, 6-8, 9-11, 12-19, 20-23) is microscopically analyzed in this study with huge volume and speed data. It was found that speed distributions are changed and show different behaviors under the conditions of traffic volume levels and daily time periods. V85s (85th percentile speed) at early morning and night periods decrease when levels of traffic volume increase, whereas V85s at day time do not show considerable differences. Especially, FFSs of levels of traffic volume at early morning and night periods, despite low density, is analyzed as heterogeneity, but homogeneity at day time.

• Journal of Korea Port Economic Association
• /
• v.26 no.3
• /
• pp.130-142
• /
• 2010

The purpose of this paper is to propose applying agent-based model(ABM) for analyzing logistics system. Logistics problems become more complex and multi-faceted. As a result, the behavior of the system becomes more and more difficult to predict, we see the limitations of the traditional top-down approach in handling complexity. The ABM, that is the bottom-up approach, provides new modeling framework in system modeling. The ABM focuses on the interactions of subsystem or agents in whole system. Then the macroscopic picture of the whole system behavior is emerged as the microscopic interactions of agents are aggregated. The ABM assumes that each agent acts based on simple rules learned from dynamic interactions among other agents or its surrounding environment. The ABM has a great advantage in understanding emergent phenomenon that cannot be explained only through considering individual attributes. The ABM is an extremely useful method to analyze complex system such as logistics system. Therefore, the great research efforts and applications on the ABM to logistics system are encouraged in future.

• Tunnel and Underground Space
• /
• v.26 no.3
• /
• pp.235-252
• /
• 2016

This study presents a methodology to reproduce the mechanical behavior of isotropic or transversely isotropic rock using the polygonal grain-based distinct element model. A numerical technique to monitor the evolution of micro-cracks during the simulation was developed in the present study, which enabled us to examine the contribution of tensile cracking and shear cracking to the progressive process of the failure. The numerical results demonstrated good agreement with general observations from rock specimens in terms of the behavior and the evolution of micro-cracks, suggesting the capability of the model to represent the mechanical behavior of rock. We also carried out a parametric study as a fundamental work to examine the relationships between the microscopic properties of the constituents and the macroscopic behavior of the model. Depending on the micro-properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics. In addition, a numerical technique to reproduce the transversely isotropic rock was suggested and applied to Asan gneiss from Korea. The behavior of the numerical model was in good agreement with the results obtained in the laboratory-scale experiments of the rock.

• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.3117-3129
• /
• 2022

Deeply understanding the phase change of thin liquid sodium film inside wick pore is very important for further studying high-temperature sodium heat pipe's heat transfer. For the first time, the evaporation and condensation of thin liquid sodium film are investigated by the Lennard-Jones potential of molecular dynamics. Based on the startup and normal operation of the sodium heat pipe, three different cases are simulated. First, the equilibrium is achieved and the Mass Accommodation Coefficients of the three cases are 0.3886, 0.2119, 0.2615 respectively. Secondly, the non-equilibrium is built. The change of liquid film thickness, the number of gas atoms, the net evaporation flux (J_net), the heat transfer coefficient (h) at the liquid-gas interface are acquired. Results indicate that the magnitude of the J_net and the h increase with the basic equilibrium temperature. In 520-600 K (the startup of the heat pipe), the h has approached 5-6 W m^-2 K^-1 while liquid film thickness is in 11-13 nm. The fact shows that during the initial startup of the sodium heat pipe, the thermal resistance at the liquid-gas interface can't be negligible. This work is the complement and extension for macroscopic investigation of heat transfer inside sodium heat pipe. It can provide a reference for further numerical simulation and optimal design of the sodium heat pipe in the future.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.42 no.5
• /
• pp.681-688
• /
• 2022

According to a recent study, most construction companies in Korea do not manage risk effectively, and it is judged that the risk management system needs to be improved. In addition, most risk-related studies deal with risks from a macroscopic perspective, and there are few studies dealing with process risks at the project construction stage. Therefore, this study tried to suggest a risk response plan through analysis and classification of risk factors that may occur in retaining work among process risks. To this end, a workshop was held for risk experts to identify and analyze risks that may occur during the construction of retaining work for apartments. As a result of the study, it was expected that savings of KRW 4.97 billion would be possible in the 95 % confidence interval, and the maximum possible cost was reduced from KRW 15 billion to about KRW 10 billion. Based on the risk reduction ratio, it was found that risks that can be reduced without any special input cost, risks with large effects in response to risks, and risks with insignificant effects were found. Therefore, using the types and risk factors presented in this study as guides, it is expected that it will be helpful in successfully operating the project if an appropriate response strategy is prepared and systematically responded to the site conditions.

• Journal of Korean Society of Transportation
• /
• v.20 no.3
• /
• pp.149-158
• /
• 2002

For the purpose of preciously describing real time traffic pattern in urban road network, dynamic network loading(DNL) models able to simulate traffic behavior are required. A number of different methods are available, including macroscopic, microscopic dynamic network models, as well as analytical model. Equivalency minimization problem and Variation inequality problem are the analytical models, which include explicit mathematical travel cost function for describing traffic behaviors on the network. While microscopic simulation models move vehicles according to behavioral car-following and cell-transmission. However, DNL models embedding such travel time function have some limitations ; analytical model has lacking of describing traffic characteristics such as relations between flow and speed, between speed and density Microscopic simulation models are the most detailed and realistic, but they are difficult to calibrate and may not be the most practical tools for large-scale networks. To cope with such problems, this paper develops a new DNL model appropriate for dynamic traffic assignment(DTA), The model is combined with vertical queue model representing vehicles as vertical queues at the end of links. In order to compare and to assess the model, we use a contrived example network. From the numerical results, we found that the DNL model presented in the paper were able to describe traffic characteristics with reasonable amount of computing time. The model also showed good relationship between travel time and traffic flow and expressed the feature of backward turn at near capacity.

• Journal of the Mineralogical Society of Korea
• /
• v.32 no.1
• /
• pp.41-49
• /
• 2019

Probing the Na environments in Na silicate and aluminosilicate glasses is essential to the macroscopic properties of melts in the Earth. In particular, exploring the atomic structure of Na silicate and aluminosilicate glasses reveals Na-O distance, which plays an important role in transport properties of melts. Here we report the local environment around Na using $^{23}Na$ magic angle spinning (MAS) NMR. We also obtain $^{23}Na$ isotropic chemical shift (${\delta}_{iso}$) of Na silicate and aluminosilicate glasses with varying composition using Dmfit program. The Q mas 1/2 model simulates the experimental results with three simulated peaks while the CzSimple model simulates with one peak. The ${\delta}_{iso}$ decreases with increasing $SiO_2$ content in Na silicate and aluminosilicate glasses. The ${\delta}_{iso}$ increases with increasing $Na_2O$ content in Na-Ca silicate and Na aluminosilicate glasses when the $SiO_2$ content is fixed. Considering the ${\delta}_{iso}$ of Na aluminosilicate glasses available in the previous studies, together with the current simulation results, we confirm that the ${\delta}_{iso}$ has positive correlation with Al / (Al + Si). Those experimental results were reproduced better using Q mas 1/2 model. The disorder of Na in Na silicate and aluminosilicate glasses can be revealed through the simulation of 1D $^{23}Na$ MAS NMR spectra using Dmfit program in a short time.

• Korean Journal of Mineralogy and Petrology
• /
• v.34 no.3
• /
• pp.157-167
• /
• 2021

Lead (Pb) is one of the key trace elements, exhibiting a peculiar partitioning behavior into silicate melts in contact with minerals. Partitioning behaviors of Pb between silicate mineral and melt have been known to depend on melt composition and thus, the atomic structures of corresponding silicate liquids. Despite the importance, detailed structural studies of Pb-bearing silicate melts are still lacking due to experimental difficulties. Here, we explored the effect of lead content on the atomic structures, particularly the evolution of silicate networks in Pb-bearing sodium metasilicate ([(PbO)_x(Na₂O)_1-x]·SiO₂) glasses as a model system for trace metal bearing natural silicate melts, using ²⁹Si solid-state nuclear magnetic resonance (NMR) spectroscopy. As the PbO content increases, the ²⁹Si peak widths increase, and the maximum peak positions shift from -76.2, -77.8, -80.3, -81.5, -84.6, to -87.7 ppm with increasing PbO contents of 0, 0.25, 0.5, 0.67, 0.86, and 1, respectively. The ²⁹Si MAS NMR spectra for the glasses were simulated with Gaussian functions for Qⁿ species (SiO₄ tetrahedra with n BOs) for providing quantitative resolution. The simulation results reveal the evolution of each Qⁿ species with varying PbO content. Na-endmember Na₂SiO₃ glass consists of predominant Q² species together with equal proportions of Q¹ and Q³. As Pb replaces Na, the fraction of Q² species tends to decrease, while those for Q¹ and Q³ species increase indicating an increase in disproportionation among Qⁿ species. Simulation results on the ²⁹Si NMR spectrum showed increases in structural disorder and chemical disorder as evidenced by an increase in disproportionation factor with an increase in average cation field strengths of the network modifying cations. Changes in the topological and configurational disorder of the model silicate melt by Pb imply an intrinsic origin of macroscopic properties such as element partitioning behavior.