• Journal of the Korean Geotechnical Society
• /
• v.27 no.11
• /
• pp.39-45
• /
• 2011

This study presents the thermal behaviors of glass beads-rubber mixtures depending on the volumetric fraction of each constituent and relative size between them. The transient plane source method is used to measure the effective thermal conductivity of mixtures. The discrete element method (DEM) and the thermal network model are integrated to investigate the particle-scale mechanism of heat transfer in granular packings. Results show that 1) the effective thermal conductivity decreases as the rubber fraction increases, and 2) the relative size between two solid particles dominates the spatial configuration of inter-particle contact condition that in tum determines the majority of heat propagation path through particle contacts. For the mixtures whose volumetric fraction of rubber is identical, the less conductive materials (e.g., rubber particles) with a large size facilitate heat transfer in granular materials. The experimental results and particle-scale observation highlight that the thermal conduction behavior is dominated not only by the volumetric fraction but also the spatial configuration of each constituent.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.2
• /
• pp.145-152
• /
• 2017

Mixing of molecular grains having different characteristics is very important in many industries such as the food and pharmaceutical industries. With the development of computer simulations, it is common practice to find the optimal mixing conditions through a simulation before the actual mixing task to estimate the proper level of mixing. Accordingly, there has been an increasing need for a mixing index to measure the mix of particles in the simulation process. Mixing indices, which have been widely used so far, can largely be classified into two types: first is the statistical-based mixing index, which is prepared using the sampling method, and the second is the mixing index that is prepared using all the particles. In this paper, we calculated mixing indices in different ways for the data in the course of mixing the particles using the DEM simulation. Additionally, we compared the performance, advantages, and disadvantages of each mixing index. Therefore, I propose a standard that can be used to select an appropriate mixing index.

• Journal of the Korean Geotechnical Society
• /
• v.33 no.11
• /
• pp.83-95
• /
• 2017

This study examined the magnitude and distribution of earth pressure on the excavation wall in jointed rock mass by considering different groundwater conditions under various rock types, joint inclination angles, and earth pressure coefficients. Based on a physical model test (Son and Park, 2014), extended studies were conducted considering rock-structure interactions based on the discrete element method, which can consider the joints characteristics of rock mass. The results showed that the earth pressure was highly influenced by the groundwater condition as well as the rock type, joint inclination angle, and earth pressure coefficient. The results were also compared with Peck's earth pressure for soil ground, and clearly showed that the earth pressure in jointed rock mass can be greatly different from that in soil ground.

• Geomechanics and Engineering
• /
• v.22 no.5
• /
• pp.385-396
• /
• 2020

A uniaxial compression test was performed to analyse the mechanical properties and macroscale and mesoscale failure mechanisms of sandstone with pyrite concretions. The effect of the pyrite concretions on the evolution of macroscale cracks in the sandstone was further investigated through numerical simulations with Particle Flow Code in 2D (PFC2D). The results revealed that pyrite concretions substantially influence the mechanical properties and macroscale and mesoscale failure characteristics of sandstone. During the initial loading stage, significant stress concentrations occurred around the edges of the pyrite concretion accompanied by the preferential generation of cracks. Meanwhile, the events and cumulative energy counts of the acoustic emission (AE) signal increased rapidly because of friction sliding between the concretion and sandstone matrix. As the axial stress increased, the degree of the stress concentration remained relatively unchanged around the edges of the concretions. The cracks continued growing rapidly around the edges of the concretions and gradually expanded toward the centre of the sample. During this stage, the AE events and cumulative energy counts increased quite slowly. As the axial stress approached the peak strength of the sandstone, the cracks that developed around the edges of the concretion started to merge with cracks that propagated at the top-left and bottom-right corners of the sample. This crack evolution ultimately resulted in the shear failure of the sandstone sample around the edges of the pyrite concretions.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.330-339
• /
• 2008

Ground-based interceptors(GBI) comprise a major element of the strategic defense against hostile targets like Intercontinental Ballistic Missiles(ICBM) and reentry vehicles(RV) dispersed from them. An optimum design of the subsystems is required to increase the performance and reliability of these GBI. Propulsion subsystem design and optimization is the motivation for this effort. This paper describes an effort in which an entire GBI missile system, including a multi-stage solid rocket booster, is considered simultaneously in a Genetic Algorithm(GA) performance optimization process. Single goal, constrained optimization is performed. For specified payload and miss distance, time of flight, the most important component in the optimization process is the booster, for its takeoff weight, time of flight, or a combination of the two. The GBI is assumed to be a multistage missile that uses target location data provided by two ground based RF radar sensors and two low earth orbit(LEO) IR sensors. 3Dimensional model is developed for a multistage target with a boost phase acceleration profile that depends on total mass, propellant mass and the specific impulse in the gravity field. The monostatic radar cross section (RCS) data of a three stage ICBM is used. For preliminary design, GBI is assumed to have a fixed initial position from the target launch point and zero launch delay. GBI carries the Kill Vehicle(KV) to an optimal position in space to allow it to complete the intercept. The objective is to design and optimize the propulsion system for the GBI that will fulfill mission requirements and objectives. The KV weight and volume requirements are specified in the problem definition before the optimization is computed. We have considered only continuous design variables, while considering discrete variables as input. Though the number of stages should also be one of the design variables, however, in this paper it is fixed as three. The elite solution from GA is passed on to(Sequential Quadratic Programming) SQP as near optimal guess. The SQP then performs local convergence to identify the minimum mass of the GBI. The performance of the three staged GBI is validated using a ballistic missile intercept scenario modeled in Matlab/SIMULINK.

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

• Korean Journal of Materials Research
• /
• v.31 no.5
• /
• pp.286-295
• /
• 2021

This study attempts to find optimal conditions of the friction coefficient using a discrete element method (DEM) simulation with various friction coefficient conditions and three different grinding media with various ball sizes in a traditional ball mill (TBM). Using ball motion of the DEM simulation are obtained using the optimal friction coefficient compared with actual motion; photographs are taken by the digital camera and the snapshot images are analyzed. In the simulation, the rotation speed of the mill, the materials and velocity of the grinding media, and the friction coefficient between the balls and the wall of the pot are fixed as the actual experimental conditions. We observe the velocity according to the friction coefficient from the DEM simulation. The friction coefficient is found to increase with the velocity. Milling experiments using a traditional ball mill with the same experimental conditions as those of the DEM simulation are conducted to verify the simulated results. In addition, particle morphology change of copper powder is investigated and analyzed using scanning electron microscopy (SEM) for the milling experiment.

• Structural Engineering and Mechanics
• /
• v.69 no.1
• /
• pp.11-20
• /
• 2019

The biaxial failure mechanism of transversally bedding concrete layers was numerically simulated using a sophisticated two-dimensional discrete element method (DEM) implemented in the particle flow code (PFC2D). This numerical modelling code was first calibrated by uniaxial compression and Brazilian testing results to ensure the conformity of the simulated numerical model's response. Secondly, 21 rectangular models with dimension of $54mm{\times}108mm$ were built. Each model contains two transversely bedding layers. The first bedding layer has low mechanical properties, less than mechanical properties of intact material, and second bedding layer has high mechanical properties, more than mechanical properties of intact material. The angle of first bedding layer, with weak mechanical properties, related to loading direction was $0^{\circ}$, $15^{\circ}$, $30^{\circ}$, $45^{\circ}$, $60^{\circ}$, $75^{\circ}$ and $90^{\circ}$ while the angle of second layer, with high mechanical properties, related to loading direction was $90^{\circ}$, $105^{\circ}$, $120^{\circ}$, $135^{\circ}$, $150^{\circ}$, $160^{\circ}$ and $180^{\circ}$. Is to be note that the angle between bedding layer was $90^{\circ}$ in all bedding configurations. Also, three different pairs of the thickness were chosen in models, i.e., 5 mm/10 mm, 10 mm/10 mm and 20 mm/10 mm. The result shows that in all configurations, shear cracks develop between the weaker bedding layers. Shear cracks angel related to normal load change from $0^{\circ}$ to $90^{\circ}$ with increment of $15^{\circ}$. Numbers of shear cracks are constant by increasing the bedding thickness. It's to be noted that in some configuration, tensile cracks develop through the intact area of material model. There is not any failure in direction of bedding plane interface with higher strength.

• Wind and Structures
• /
• v.29 no.4
• /
• pp.247-270
• /
• 2019

Wind-resistant design of existing cooling tower structures overlooks the impacts of rainfall. However, rainstorm will influence aerodynamic force on the tower surface directly. Under this circumstance, the structural response of the super-large cooling tower (SLCT) will become more complicated, and then the stability and safety of SLCT will receive significant impact. In this paper, surrounding wind fields of the world highest (210 m) cooling tower in Northwest China underthree typical wind velocities were simulated based on the wind-rain two-way coupling algorithm. Next, wind-rain coupling synchronous iteration calculations were conducted under 9 different wind speed-rainfall intensity combinations by adding the discrete phase model (DPM). On this basis, the influencing laws of different wind speed-rainfall intensity combinations on wind-driving rain, adhesive force of rain drops and rain pressure coefficients were discussed. The acting mechanisms of speed line, turbulence energy strength as well as running speed and trajectory of rain drops on structural surface in the wind-rain coupling field were disclosed. Moreover, the fitting formula of wind-rain coupling equivalent pressure coefficient of the cooling tower was proposed. A systematic contrast analysis on its 3D distribution pattern was carried out. Finally, coupling model of SLCT under different working conditions was constructed by combining the finite element method. Structural response, buckling stability and local stability of SLCT under different wind velocities and wind speed-rainfall intensity combinations were compared and analyzed. Major research conclusions can provide references to determine loads of similar SLCT accurately under extremely complicated working conditions.

• Smart Structures and Systems
• /
• v.23 no.5
• /
• pp.479-493
• /
• 2019

In this paper, the interaction between notch and micro pore under uniaxial compression has been performed experimentally and numerically. Firstly calibration of PFC2D was performed using Brazilian tensile strength, uniaxial tensile strength and biaxial tensile strength. Secondly uniaxial compression test consisting internal notch and micro pore was performed experimentally and numerically. 9 models consisting notch and micro pore were built, experimentally and numerically. Dimension of these models are 10 cm*1 cm*5 cm. the length of joint is 2 cm. the angularities of joint are $30^{\circ}$, $45^{\circ}$ and $60^{\circ}$. For each joint angularity, micro pore was situated 2 cm above the lower tip of the joint, 2 cm above the middle of the joint and 2 cm above the upper of the joint, separately. Dimension of numerical models are 5.4 cm*10.8 cm. The size of the cracks was 2 cm and its orientation was $30^{\circ}$, $45^{\circ}$ and $60^{\circ}$. Diameter of pore was 1cm which situated at the upper of the notch i.e., 2 cm above the upper notch tip, 2 cm above the middle of the notch and 2 cm above the lower of the notch tip. The results show that failure pattern was affected by notch orientation and pore position while uniaxial compressive strength is affected by failure pattern.