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

• Journal of the Korean Society of Mineral and Energy Resources Engineers
• /
• v.55 no.6
• /
• pp.614-623
• /
• 2018

The recent occurrence of consecutive large earthquakes in the southeastern part of the Korean peninsula has brought significant attention to the prevention of earthquake damage in Korea. This article aims to explore a technology-based approach for earthquake drills using state-of-the-art visual intelligence and virtual reality technologies. The technical process consists of several stages, including acquisition of image information in living spaces using a camera, recognition of objects from the acquired image information, extraction of three dimensional geometric information, simulation of virtual earthquakes using dynamic modelling techniques such as the discrete element method, and realization of the simulated earthquake in a virtual reality environment. This article provides a comprehensive analysis of the individual processes at each stage of the technical process, a survey on the current status of related technologies, and discussion of the technical challenges in its execution.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.1
• /
• pp.15-23
• /
• 2008

In this study, modeling shear characteristics of a coarse material mainly containing crushed stones were implemented using PFC2D, a commercially available code based on DEM(Discrete Element Method). Using the DEM code, this study provides the methodology considering the shear characteristics due to a irregular grain shape, GSD(Grain Size Distribution) and porosity of coarse material which are not effectively incorporated in conventional continuum numerical codes. Direct shear test was simulated for the GSD and porosity generated sample using the code and the simulated results showed very good agreement with the laboratory test results. The current modeling approach can be applied to other coarse materials having various GSD and porosities. Using such application, prediction of the strength characteristics of coarse material in field scale would be possible, which is limited in laboratory scale so far.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.

• Journal of the Korean Geosynthetics Society
• /
• v.9 no.1
• /
• pp.47-57
• /
• 2010

Granular soils having a large particle size have been used as a filling material in the construction of foundation, harbor, dam, and so on. Consequently, the shear behavior of this granular soil plays a key role in respect of stability of structures. For example, soil particle crushing occurring at the interface between structure and soil and/or within soil mass can cause a disturbance of ground characteristics and consequently induce issues in respect of stability of structures. In order to investigate the shear behavior according to an existence and nonexistence of particle crushing, numerical analyses were conducted by using the DEM (Discrete Element Method)-based software program PFC2D (Particle Flow Code). By dividing soil particle bonding model into crushing model and noncrushing model, total four particle bonding models were simulated and their results were compared. Noncrushing model included one ball model and clump model, and crushing model included cluster model and Lobo-crushing model. The combinations of soil particle followed the research results of Lobo-Guerrero and Vallejo (2005) which were composed of eight circles. The results showed that the friction angle was in order of clump model > cluster model > one ball model. The particle bonding model compared to one ball model and noncrushing model compared to crushing model showed higher shear strength. It was also concluded that the model suggested by Lobo-Guerrero and Vallejo (2005) is not appropriate to simulate the soil particle crushing.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.7
• /
• pp.265-274
• /
• 2018

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has prompted studies directed toward the development of an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have been obtained traditionally using standard finite difference or finite element methods. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under the conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method (DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as the excessive use of storage due to the conditions of complex geometry and loading. The in-plane buckling of curved beams considering the extensibility of the arch axis was analyzed under uniformly distributed radial loads using the DQM. The critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with the precise results by other methods for cases, in which they were available. The DQM, using only a limited number of grid points, provided results that agreed very well (less than 0.3%) with the exact ones. New results according to diverse variations were obtained, showing the important roles in the buckling behavior of curved beams, and can be used in comparisons with other numerical solutions or with experimental test data.

• Journal of Korean Association for Spatial Structures
• /
• v.4 no.2 s.12
• /
• pp.63-71
• /
• 2004

The structural space is gradually wide and is wanting agreeable environment by the requirement and necessity of people who lives modem stage. The building coincides with such requirements and is the high rise building actual circumstances which is doing ultra-large. The confirmed report of the technology to organize great merit is becoming currently considerably important issue in constructing a building field. Thus, this paper examine closely for nonlinear unstable taking a picture uneasiness height of prosperity considering to initial imperfection by a numerical method with a space frame structure of discrete system in large space structure. Based on previous investigation method, this paper induce nodal stiffness matrix of solid truss elements considering geometrical nonlinear using finite element method. In this paper, three types of space structure considered; i) 1-free node space structure, ii) 2-free node space structure, iii) multi-free node space structure. It apply the above examples to a nonlinear program, next, grasp the characteristic of an unstable conduct and the result was a clearing low.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.10
• /
• pp.300-309
• /
• 2016

Curved beams are increasingly used in buildings, vehicles, ships, and aircraft, which has resulted in considerable effort towards developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of many investigations. Solutions to the relevant differential equations have traditionally been obtained by the standard finite difference or finite element methods. However, these techniques require a great deal of computer time for a large number of discrete nodes with conditions of complex geometry and loading. One efficient procedure for the solution of partial differential equations is the differential quadrature method (DQM). This method has been applied to many cases to overcome the difficulties of complex algorithms and high storage requirements for complex geometry and loading conditions. Out-of-plane buckling of curved beams with rotatory inertia were analyzed using DQM under uniformly distributed radial loads. Critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with exact results from other methods for available cases. The DQM used only a limited number of grid points and shows very good agreement with the exact results (less than 0.3% error). New results according to diverse variation are also suggested, which show important roles in the buckling behavior of curved beams and can be used for comparisons with other numerical solutions or experimental test data.

• Structural Engineering and Mechanics
• /
• v.86 no.5
• /
• pp.607-619
• /
• 2023

This study, it was tried to evaluate the asphalt behavior under tensile loading conditions through indirect Brazilian and direct tensile tests, experimentally and numerically. This paper is important from two points of view. The first one, a new test method was developed for the determination of the direct tensile strength of asphalt and its difference was obtained from the indirect test method. The second one, the effects of particle size and loading rate have been cleared on the tensile fracture mechanism. The experimental direct tensile strength of the asphalt specimens was measured in the laboratory using the compression-to-tensile load converting (CTLC) device. Some special types of asphalt specimens were prepared in the form of slabs with a central hole. The CTLC device is then equipped with this specimen and placed in the universal testing machine. Then, the direct tensile strength of asphalt specimens with different sizes of ingredients can be measured at different loading rates in the laboratory. The particle flow code (PFC) was used to numerically simulate the direct tensile strength test of asphalt samples. This numerical modeling technique is based on the versatile discrete element method (DEM). Three different particle diameters were chosen and were tested under three different loading rates. The results show that when the loading rate was 0.016 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis till coalescence to the model boundary. When the loading rate was 0.032 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis. The branching occurs in these cracks. This shows that the crack propagation is under quasi-static conditions. When the loading rate was 0.064 mm/sec, mixed tensile and shear cracks were initiated below the loading walls and branching occurred in these cracks. This shows that the crack propagation is under dynamic conditions. The loading rate increases and the tensile strength increases. Because all defects mobilized under a low loading rate and this led to decreasing the tensile strength. The experimental results for the direct tensile strengths of asphalt specimens of different ingredients were in good accordance with their corresponding results approximated by DEM software.