• Advances in Computational Design
• /
• v.5 no.3
• /
• pp.277-289
• /
• 2020

This article describes the technology for constructing of a multiply-connected three-dimensional area's finite element representation. Representation of finite-element configuration of an area is described by a discrete set that consist of the number of nodes and elements of the finite-element grid, that are orderly set of nodes' coordinates and numbers of finite elements. Corresponding theorems are given, to prove the correctness of the solution method. The adequacy of multiply-connected area topology's finite element model is shown. The merging of subareas is based on the criterion of boundary nodes' coincidence by establishing a simple hierarchy of volumes, surfaces, lines and points. Renumbering nodes is carried out by the frontal method, where nodes located on the outer edges of the structure are used as the initial front.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1996.04a
• /
• pp.247-252
• /
• 1996

In the concrete structures, cracks occur in various causes and the cracks seriously affect the functions of structures. The analysis techniques of progressive crack in the concrete have been improved with the advance of numerical techniques. The discrete crack model used in finite element program for the analysis of progressive failure is very effective, but it can not be easily implemented into numerical procedures because of difficult handing of nodal points in finite element meshes for crack growth. This paper introduces one of the techniques which skips the difficulty. In this paper, the modeling of progressive failure using finite element formulation is explained for the analysis of concrete fracture. The discontinuous element using the discontinuous shape function and the dual mapping technique in the numerical integration are implemented into finite element code for this purpose. It is shown that developed finite element program can predict the quasi-brittle behavior of concrete including ultimate load. The comparisons of the analysis results with other data are also shown.

• Coupled systems mechanics
• /
• v.4 no.4
• /
• pp.365-383
• /
• 2015

Strip footing is an important type of shallow foundations and is commonly used beneath the walls. Analysis of shallow foundation involves the determination of stresses and deformations. Element free Galerkin method, one of the important mesh free methods, is used for the determination of stresses and deformations. Element free Galerkin method is an efficient and accurate method as compared to finite element method. The Element Free Galerkin method uses only a set of nodes and a description of model boundary is required to generate the discrete equation. Strip footing of width 2 m subjected to a loading intensity of 200 kPa is studied. The results obtained are agreeing with the values obtained using analytical solutions available in the literature. Parametric study is done and the effect of modulus of deformation, Poisson's ratio and scaling parameter on deformation and stresses are determined.

• Smart Structures and Systems
• /
• v.26 no.5
• /
• pp.605-617
• /
• 2020

In this study, the experimental tests for the direct tensile strength measurement of Ultra-High Performance Concrete (UHPC) were numerically modeled by using the discrete element method (circle type element) and Finite Element Method (FEM). The experimental tests used for the laboratory tensile strength measurement is the Compressive-to-Tensile Load Conversion (CTLC) device. In this paper, the failure process including the cracks initiation, propagation and coalescence studied and then the direct tensile strength of the UHPC specimens measured by the novel apparatus i.e., CTLC device. For this purpose, the UHPC member (each containing a central hole) prepared, and situated in the CTLC device which in turn placed in the universal testing machine. The direct tensile strength of the member is measured due to the direct tensile stress which is applied to this specimen by the CTLC device. This novel device transferring the applied compressive load to that of the tensile during the testing process. The UHPC beam specimen of size 150 × 60 × 190 mm and internal hole of 75 × 60 mm was used in this study. The rate of the applied compressive load to CTLC device through the universal testing machine was 0.02 MPa/s. The direct tensile strength of UHPC was found using a new formula based on the present analyses. The numerical simulation given in this study gives the tensile strength and failure behavior of the UHPC very close to those obtained experimentally by the CTLC device implemented in the universal testing machine. The percent variation between experimental results and numerical results was found as nearly 2%. PFC2D simulations of the direct tensile strength measuring specimen and ABAQUS simulation of the tested CTLC specimens both demonstrate the validity and capability of the proposed testing procedure for the direct tensile strength measurement of UHPC specimens.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.11-19
• /
• 2005

First, methods of numerical analysis of gas-particle flows is classified into micro, meso and macro scale approaches based on the concept of multi-scale mechanics. Next, the explanation moves on to discrete particle simulation where motion of individual particles is calculated numerically using the Newtonian equations of motion. The author focuses on the cases where particle-to-particle interaction has significant effects on the phenomena. Concerning the particle-to-particle interaction, two cases are considered: the one is collision-dominated flows and the other is the contact-dominated flows. To treat this interaction mathematically, techniques named DEM(Distinct Element Method) or DSMC (Direct Simulation Monte Carlo) have been developed DEM, which has been developed in the field of soil mechanics, is useful for the contact -dominated flows and DSMC method, developed in molecular gas flows, is for the collision-dominated flows. Combining DEM or DSMC with CFD (computer fluid dynamics), the discrete particle simulation becomes a more practical tool for industrial flows because not only the particle-particle interaction but particle-fluid interaction can be handled. As examples of simulations, various results are shown, such as hopper flows, particle segregation phenomena, particle mixing in a rotating drum, dense phase pneumatic conveying, spouted bed, dense phase fluidized bed, fast circulating fluidized bed and so on.

• Geomechanics and Engineering
• /
• v.35 no.5
• /
• pp.475-486
• /
• 2023

Ballast particles have an irregular shape and are discrete in nature. Due to the discrete nature of ballast, it exhibits complex mechanical behaviour under loading conditions. The discrete element method (DEM) can model the behaviour of discrete particles under a multitude of loading conditions. DEM is used in this paper to simulate a series of three-dimensional direct shear tests in order to investigate the shear behaviour of railway ballast and its interaction at the microscopic level. Particle flow code in three dimension (PFC3D) models the irregular shape of ballast particles as clump particles. To investigate the influence of particle size distribution (PSD), real PSD of Indian railway ballast specification IRS:GE:1:2004, China high-speed rail (HSR) and French rail specifications are generated. PFC3D built-in linear contact model is used to simulate the interaction of ballast particles under various normal stresses, shearing rate and shear box sizes. The results indicate how shear resistance and volumetric changes in ballast assembly are affected by normal stress, shearing rate, PSD and shear box size. In addition to macroscopic behaviour, DEM represents the microscopic behaviour of ballast particles in the form of particle displacement at different stages of the shearing process.

• Journal of the Korean Mathematical Society
• /
• v.41 no.4
• /
• pp.681-715
• /
• 2004

Mathematical formulation and numerical solutions of an optimal Dirichlet boundary control problem for the Boussinesq equations are considered. The solution of the optimal control problem is obtained by adjusting of the temperature on the boundary. We analyze finite element approximations. A gradient method for the solution of the discrete optimal control problem is presented and analyzed. Finally, the results of some computational experiments are presented.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.10a
• /
• pp.88-95
• /
• 1997

In this paper, we present an efficient method for finite element vibration analysis of constantly rotating blade disks which are deformed to some considerable extent by centrifugal force, Coriolis force and operating load, and vibrate due to several types of exciting forces. A blade disk which is a structure with cyclic symmetry is divided into substructures with the same geometry. Only one substructure is modeled and can be analysed rapidly and exactly using discrete Fourier transform by means of a computer with small memory.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.26 no.4
• /
• pp.323-332
• /
• 2022

This paper presents a numerical study on multigrid algorithms of V-cycle type for problems posed in the Hilbert space H(curl) in three dimensions. The multigrid methods are designed for discrete problems originated from the discretization using the hexahedral Nédélec edge element of the lowest-order. Our suggested methods are associated with smoothers constructed by substructuring based on domain decomposition methods of nonoverlapping type. Numerical experiments to demonstrate the robustness and the effectiveness of the suggested algorithms are also provided.

• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1129-1141
• /
• 2011

In this paper, a class of large time-stepping method based on the finite element discretization for the Cahn-Hilliard equation with the Neumann boundary conditions is developed. The equation is discretized by finite element method in space and semi-implicit schemes in time. For the first order fully discrete scheme, convergence property is investigated by using finite element analysis. Numerical experiment is presented, which demonstrates the effectiveness of the large time-stepping approaches.