• 한국터널공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.348-358
• /
• 2005

Although the evaluation of the mechanical properties and behavior of discontinuous rock masses is very important for the design of underground openings, it has always been considered the most difficult problem. One of the difficulties in describing the rock mass behavior is assigning the appropriate constitutive model. This limitation may be overcome with the progress in discrete element software such as PFC, which does not need the user to prescribe a constitutive model for rock mass. Instead, the micro-scale properties of the intact rock and joints are defined and the macro-scale response results from those properties and the geometry of the problem. In this paper, a $30m{\times}30m{\times}30m$ jointed rock mass of road tunnel site was analyzed. A discrete fracture network was developed from the joint geometry obtained from core logging and surface survey. Using the discontinuities geometry from the DFN model, PFC simulations were carried out, starting with the intact rock and systematically adding the joints and the stress-strain response was recorded for each case. With the stress-strain response curves, the mechanical properties of discontinuous rock masses were determined and compared to the results of empirical methods such as RMR, Q and GSI. The values of Young's modulus, Poisson's ratio and peak strength are almost similar from PFC model and Empirical methods. As expected, the presence of joints had a pronounced effect on mechanical properties of the rock mass. More importantly, the mechanical response of the PFC model was not determined by a user specified constitutive model.

• Tunnel and Underground Space
• /
• v.26 no.4
• /
• pp.316-326
• /
• 2016

The purpose of this study is to investigate the influence of cleavages on indirect tensile strength (ITS) of the granite. Brazilian disc tests and ring tests with three different hole sizes were performed. 2D DEM (Discrete Element Method) simulation was employed to further understand the failure process during the tests and the mechanism behind. Results show that ITS decreases across hardway, grain and rift cleavage. Measured average ITS from ring tests is about 2.5 ~ 6.4 times of those measured from Brazilian disc tests, and it decreases with increasing ratio of diameters of inner hole and specimen. Failure pattern in ring tests is influenced by both hole size and relative positions of cleavages parallel and perpendicular to the loading direction.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.4
• /
• pp.465-474
• /
• 2012

In this study, a dynamic simulation model that considers many spherical particles and multibody dynamics (MBD) entities is developed. Plenteous spherical particles are solved using the Discrete Element Method (DEM) technique and simulated on a GPU board in a PC. A fast algorithm is used to calculate the Hertzian contact forces between many spherical particles, and NVIDIA CUDA is used to increase the calculation speed. The explicit integration method is applied to solve the many spheres. MBD entities are simulated by recursive formulation. Constraints are reduced by recursive formulation, and the implicit generalized alpha method is applied to solve the dynamic model. A new algorithm is developed to simulate the DEM and MBD models simultaneously. As a numerical example, a truck car model and gear model are developed. The results show that the proposed algorithm using a general-purpose GPU in a PC has many advantages.

• Steel and Composite Structures
• /
• v.10 no.2
• /
• pp.129-149
• /
• 2010

The objective of this study is to formulate a general 3D material-structural analysis framework for the thermomechanical behavior of steel-concrete structures in a fire environment. The proposed analysis framework consists of three sequential modeling parts: fire dynamics simulation, heat transfer analysis, and a thermomechanical stress analysis of the structure. The first modeling part consists of applying the NIST (National Institute of Standards and Technology) Fire Dynamics Simulator (FDS) where coupled CFD (Computational Fluid Dynamics) with thermodynamics are combined to realistically model the fire progression within the steel-concrete structure. The goal is to generate the spatial-temporal (ST) solution variables (temperature, heat flux) on the surfaces of the structure. The FDS-ST solutions are generated in a discrete form. Continuous FDS-ST approximations are then developed to represent the temperature or heat-flux at any given time or point within the structure. An extensive numerical study is carried out to examine the best ST approximation functions that strike a balance between accuracy and simplicity. The second modeling part consists of a finite-element (FE) transient heat analysis of the structure using the continuous FDS-ST surface variables as prescribed thermal boundary conditions. The third modeling part is a thermomechanical FE structural analysis using both nonlinear material and geometry. The temperature history from the second modeling part is used at all nodal points. The ABAQUS (2003) FE code is used with external user subroutines for the second and third simulation parts in order to describe the specific heat temperature nonlinear dependency that drastically affects the transient thermal solution especially for concrete materials. User subroutines are also developed to apply the continuous FDS-ST surface nodal boundary conditions in the transient heat FE analysis. The proposed modeling framework is applied to predict the temperature and deflection of the well-documented third Cardington fire test.

• Geomechanics and Engineering
• /
• v.7 no.3
• /
• pp.317-333
• /
• 2014

The study of the mechanical behavior of rockfill materials under three-dimensional loading conditions is a current research focus area. This paper presents a microscale numerical study of rockfill deformation and strength characteristics using the Combined Finite-Discrete Element Method (FDEM). Two features unique to this study are the consideration of irregular particle shapes and particle crushability. A polydisperse assembly of irregular polyhedra was prepared to reproduce the mechanical behavior of rockfill materials subjected to axial compression at a constant mean stress for a range of intermediate principal stress ratios in the interval [0, 1]. The simulation results, including the stress-strain characteristics, relationship between principal strains, and principal deviator strains are discussed. The stress-dilatancy behavior is described using a linear dilatancy equation with its material constants varying with the intermediate principal stress ratio. The failure surface in the principal stress space and its traces in the deviatoric and meridian plane are also presented. The modified Lade-Duncan criterion most closely describes the stress points at failure.

• Structural Engineering and Mechanics
• /
• v.66 no.4
• /
• pp.453-463
• /
• 2018

Three points bending flexural test was modelled numerically to study the crack propagation in the pre-cracked beams. The pre-existing double internal cracks inside the beam models were considered to investigate the crack propagation and coalescence paths within the modelled samples. Notch configuration effects on the failure stress were considered too. This numerical analysis shown that the propagation of wing cracks emanating from the tips of the pre-existing internal cracks caused the final breaking of beams specimens. It was also shown that when two notches were overlapped, they both mobilized in the failure process and the failure stress was decreased when the notches were located in centre line. However, the failure stress was increased by increasing the bridge area angle. Finally, it was shown that in all cases, there were good agreements between the discrete element method results and, the other numerical and experimental results. In this research, it is tried to improve the understanding of the crack propagation and crack coalescence phenomena in brittle materials which is of paramount importance in the stability analyses of rock and concrete structures, such as the underground openings, rock slopes and tunnel construction.

• Geomechanics and Engineering
• /
• v.22 no.1
• /
• pp.97-108
• /
• 2020

The time-consuming and less objectivity are the main problems of conventional micromechanical parameters calibration method of Particle Flow Code simulations. Thus this study aims to address these two limitation of the conventional "trial-and-error" method. A new calibration method for the linear parallel bond model (CM-LPBM) is proposed. First, numerical simulations are conducted based on the results of the uniaxial compression tests on limestone. The macroscopic response of the numerical model agrees well with the results of the uniaxial compression tests. To reduce the number of the independent micromechanical parameters, numerical simulations are then carried out. Based on the results of the orthogonal experiments and the multi-factor variance analysis, main micromechanical parameters affecting the macro parameters of rocks are proposed. The macro-micro parameter functions are ultimately established using multiple linear regression, and the iteration correction formulas of the micromechanical parameters are obtained. To further verify the validity of the proposed method, a case study is carried out. The error between the macro mechanical response and the numerical results is less than 5%. Hence the calibration method, i.e., the CM-LPBM, is reliable for obtaining the micromechanical parameters quickly and accurately, providing reference for the calibration of micromechanical parameters.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.547-552
• /
• 2003

In mold machining, there are many concave machining regions where chatter and tool deflection occur since MRR (material removal rate) increases as curvature increases even though cutting speed and depth of cut are constant. Boolean operation between stock and tool model is widely used to compute MRR in NC milling simulation. In finish cutting, the side step is reduced to about 0.3mm and tool path length is sometimes over 300m. so Boolean operation takes long computation time and includes much error if the resolution of stock and tool model is larger than the side step. In this paper, curvature of CL(cutter location) surface and side step of tool path is used to compute the feedrate for constant MRR machining. The data structure of CL surface is Z-map generated from NC tool path. The algorithm to get local curvature from discrete data was developed and applied to compute local curvature of CL surface. The side step of tool path was computed by point density map which includes cutter location point density at each grid element. The feedrate computed from curvature and side step is inserted to new tool path to regulate MRR. The resultants wire applied to feedrate optimization system which generates new tool path with feedrate from NC codes for finish cutting. The system was applied to speaker mold machining. The finishing time was reduced to 12.6%. tool wear was reduced from 2mm to 1.1mm and chatter marks and over cut on corner were removed.

• Computers and Concrete
• /
• v.20 no.1
• /
• pp.39-47
• /
• 2017

In this paper, mode I fracture toughness of rock was determined by direct and indirect methods using Particle Flow Code simulation. Direct methods are compaction tension (CT) test and hollow centre cracked quadratic sample (HCCQS). Indirect methods are notched Brazilian disk (NBD) specimen, the semi-circular bend (SCB) specimen, hollow centre cracked disc (HCCD), the single edge-notched round bar in bending (SENRBB) specimen and edge notched disk (END). It was determined that which one of indirect fracture toughness values is close to direct one. For this purpose, initially calibration of PFC was undertaken with respect to data obtained from Brazilian laboratory tests to ensure the conformity of the simulated numerical models response. Furthermore, the simulated models in five introduced indirect tests were cross checked with the results from direct tests. By using numerical testing, the failure process was visually observed. Discrete element simulations demonstrated that the macro fractures in models are caused by microscopic tensile breakages on large numbers of bonded discs. Mode I fracture toughness of rock in direct test was less than other tests results. Fracture toughness resulted from semi-circular bend specimen test was close to direct test results. Therefore semi-circular bend specimen can be a proper test for determination of Mode I fracture toughness of rock in absence of direct test.

• The Transactions of the Korean Institute of Electrical Engineers P
• /
• v.64 no.1
• /
• pp.35-40
• /
• 2015

Generators are the most complex and expensive single element in a power system. The generator protection relays should to minimize damage during fault states and must be designed for maximum reliability. A conventional CDR(Current Differential Relaying) technique based on DFT(Discrete Fourier Transform) filter have the disadvantages that the time information can lead to loss in the process of converting the signal from the time domain to the frequency domain. A WT(Wavelet transform) and WT analysis is known that it is possible with the local analysis of the fault and transient signal. In this paper, to overcome the defects in the DFT process, an application of WT for fault detection of generator is presented. This paper describes an selection of mother Wavelet to detect faults of generator. Using collected data from the fault simulation with ATPdraw, we analyzed the several mother Wavelet through the course of MLD(multi-level decomposition) using MATLAB software. Finally, it can be seen that the proposed technique using detail coefficient of Daubechies level 2 which can be fault discriminant of generator.