• Computers and Concrete
• /
• v.18 no.1
• /
• pp.39-51
• /
• 2016

This paper considers the tensile strength of concrete samples in direct, CTT, modified tension, splitting and ring tests using both of the experimental tests and numerical simulation (particle flow code 2D). It determined that which one of indirect tensile strength is close to direct tensile strength. Initially calibration of PFC was undertaken with respect to the data obtained from Brazilian laboratory tests to ensure the conformity of the simulated numerical models response. Furthermore, validation of the simulated models in four introduced tests was also cross checked with the results from experimental tests. By using numerical testing, the failure process was visually observed and failure patterns were watched to be reasonable in accordance with experimental results. Discrete element simulations demonstrated that the macro fractures in models are caused by microscopic tensile breakages on large numbers of bonded discs. Tensile strength of concrete in direct test was less than other tests results. Tensile strength resulted from modified tension test was close to direct test results. So modified tension test can be a proper test for determination of tensile strength of concrete in absence of direct test. Other advantages shown by modified tension tests are: (1) sample preparation is easy and (2) the use of a simple conventional compression press controlled by displacement compared with complicate device in other tests.

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

• Journal of the Korean Geotechnical Society
• /
• v.32 no.12
• /
• pp.33-43
• /
• 2016

In order to analyze ground relaxation and cavity formation mechanism due to deteriorated sewer pipe, field test was carried out and a numerical assessments were compared with the field test results. An artificial underground cavity was intended using the ice block overlaying the buried pipe and confirmed that the cavity and relaxation of the surrounding ground were gradually formed as the ice block starts to melt down. Such mechanism was highly suspected to be involved with soil particle interlocking as a soil compaction was a typical process for the buried pipes. In exploring such mechanism numerically, commercially available DEM (Discrete Element Method) code PFC2D was used and the interlocking induced cavern behaviors were successfully simulated and compared with field test results by utilizing the clump logic imbedded in PFC code.

• Geomechanics and Engineering
• /
• v.22 no.2
• /
• pp.133-142
• /
• 2020

The evolution of the mining-induced fracture network formed during longwall top coal caving (LTCC) has a great influence on the gas drainage, roof control, top coal recovery ratio and engineering safety of aquifers. To reveal the evolution of the mining-induced stress and fracture network formed during LTCC, the fracture network in front of the working face was observed by borehole video experiments. A discrete element model was established by the universal discrete element code (UDEC) to explore the local stress distribution. The regression relationship between the fractal dimension of the fracture network and mining stress was established. The results revealed the following: (1) The mining disturbance had the most severe impact on the borehole depth range between approximately 10 m and 25 m. (2) The distribution of fractures was related to the lithology and its integrity. The coal seam was mainly microfractures, which formed a complex fracture network. The hard rock stratum was mainly included longitudinal cracks and separated fissures. (3) Through a numerical simulation, the stress distribution in front of the mining face and the development of the fracturing of the overlying rock were obtained. There was a quadratic relationship between the fractal dimension of the fractures and the mining stress. The results obtained herein will provide a reference for engineering projects under similar geological conditions.

• Structural Engineering and Mechanics
• /
• v.66 no.5
• /
• pp.569-582
• /
• 2018

Effects of crack separation, bridge area, on the tensile behaviour of concrete are studied experimentally and numerically through the Brazilian tensile test. The physical data obtained from the Brazilian tests are used to calibrate the two-dimensional particle flow code based on discrete element method (DEM). Then some specially designed Brazilian disc specimens containing two parallel cracks are used to perform the physical tests in the laboratory and numerically simulated to make the suitable numerical models to be tested. The experimental and numerical results of the Brazilian disc specimens are compared to conclude the validity and applicability of these models used in this research. Validation of the simulated models can be easily checked with the results of Brazilian tests performed on non-persistent cracked physical models. The Brazilian discs used in this work have a diameter of 54 mm and contain two parallel centred cracks ($90^{\circ}$ to the horizontal) loaded indirectly under the compressive line loading. The lengths of cracks are considered as; 10 mm, 20 mm, 30 mm and 40 mm, respectively. The visually observed failure process gained through numerical Brazilian tests are found to be very similar to those obtained through the experimental tests. The fracture patterns demonstrated by DEM simulations are mostly affected by the crack separation but the tensile strength of bridge area is related to the fracture pattern and failure mechanism of the testing samples. It has also been shown that when the crack lengths are less than 30 mm, the tensile cracks may initiate from the cracks tips and propagate parallel to loading direction till coalesce with the other cracks tips while when the cracks lengths are more than 30 mm, these tensile cracks may propagate through the intact concrete itself rather than that of the bridge area.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.18 no.2
• /
• pp.129-142
• /
• 2016

The underground cavity known as one of the reasons of ground surface settlement is a discontinuous character. Therefore, it is limited to analyze with continuum analysis. In this research, The spherical underground cavity affecting the ground surface settlement is studied with Discrete Element Method. Ground properties, depth and diameter of the spherical underground cavity are chosen as factors of the spherical underground cavity and the effect of the each factor variations on the ground surface settlement is analyzed. Relative depth to the diameter of the spherical underground cavity is also studied. The result of the research suggests the basis of underground cavity collapse prediction and standard of support.

• Geomechanics and Engineering
• /
• v.23 no.5
• /
• pp.455-466
• /
• 2020

This paper presents analysis of the interaction between tunnel and Qanat with a particular interest for the optimization of Qanat shape using the discrete element code, PFC2D, and the results will be compared with the FEM results of PLAXIS2D. For these concerns, using software PFC2D based on Discrete Element Method (DEM), a model with dimension of 100m ^⁎ 100 m was prepared. A circular tunnel with dimension of 9 m was situated 20 m below the ground surface. Also one Qanat was situated perpendicularly above the tunnel roof. Distance between Qanat center and ground surface was 8 m. Five different shapes for Qanat were selected i.e., square, semi-circular, vertical ellipse, circular and horizontal ellipse. Confining pressure of 5 MPa was applied to the model. The vertical displacement of balls situated in ground surface was picked up to measure the ground subsidence. Also two measuring circles were situated at the tunnel roof and at the Qanat roof to check the vertical displacements. The properties of the alluvial soil of Tehran city are: γ_dry=19 (KN/㎥), E= 750 (kg/㎠), ν=0.35, c=0.3(kg/㎠), φ=34°. In order to validate the DEM results, a comparison between the numerical results (obtained in this study) and analytical and field monitoring have been done. The PFC2D results are compared with the FEM results. The results shows that when Qanat has rectangular shape, the tensile stress concentration at the Qanat corners has maximum value while it has minimum value for vertical ellipse shape. The ground subsidence for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. The vertical displacements at the tunnel roof for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. Historical shape of Qante approved the finding of this research.

• Nuclear Engineering and Technology
• /
• v.46 no.6
• /
• pp.783-796
• /
• 2014

The coarse mesh finite difference (CMFD) method is applied to the discontinuous finite element method based discrete ordinate calculation for source convergence acceleration. The three-dimensional (3-D) DFEM-Sn code FEDONA is developed for general geometry applications as a framework for the CMFD implementation. Detailed methods for applying the CMFD acceleration are established, such as the method to acquire the coarse mesh flux and current by combining unstructured tetrahedron elements to rectangular coarse mesh geometry, and the alternating calculation method to exchange the updated flux information between the CMFD and DFEM-Sn. The partial current based CMFD (p-CMFD) is also implemented for comparison of the acceleration performance. The modified p-CMFD method is proposed to correct the weakness of the original p-CMFD formulation. The performance of CMFD acceleration is examined first for simple two-dimensional multigroup problems to investigate the effect of the problem and coarse mesh sizes. It is shown that smaller coarse meshes are more effective in the CMFD acceleration and the modified p-CMFD has similar effectiveness as the standard CMFD. The effectiveness of CMFD acceleration is then assessed for three-dimensional benchmark problems such as the IAEA (International Atomic Energy Agency) and C5G7MOX problems. It is demonstrated that a sufficiently converged solution is obtained within 7 outer iterations which would require 175 iterations with the normal DFEM-Sn calculations for the IAEA problem. It is claimed that the CMFD accelerated DFEM-Sn method can be effectively used in the practical eigenvalue calculations involving general geometries.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.23 no.6
• /
• pp.423-438
• /
• 2021

The blade of motor grader is used for scattering and leveling the aggregates on the foundation of road construction site. The paper performed a design improvement research of the blade part to enhance the working efficiency of motor graders. The scattering works of aggregates by blade driving were simulated by DEM (discrete element method) of a dynamic code. The four design parameters were selected and a specific leveling scenario for the simulation was determined. The nine blade models were numerically experimented, and the sensitivity of each factors was analyzed. Next, the design factors that influence a blade performance have been selected by ANOVA, and these key design factors were applied to the progressive quadratic response surface method (PQRSM). The optimum set of design factors of the blade was finally proposed.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2002.10a
• /
• pp.233-236
• /
• 2002

An efficient method was developed in this study to obtain optimal stacking sequences, thicknesses, and minimum weights of stiffened laminated composite shells under combined loading conditions and stiffener layouts using genetic algorithms (GAs) and finite element analyses. Among many parameters in designing composite laminates determining a optimal stacking sequence that may be formulated as an integer programming problem is a primary concern. Of many optimization algorithms, GAs are powerful methodology for the problem with discrete variables. In this paper the optimal stacking sequence was determined, which gives the maximum critical buckling load factor and the minimum weight as well. To solve this problem, both the finite element analysis by ABAQUS and the GA-based optimization procedure have been implemented together with an interface code. Throughout many parametric studies using this analysis tool, the influences of stiffener sizes and three different types of stiffener layouts on the stacking sequence changes were throughly investigated subjected to various combined loading conditions.