• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.2A
• /
• pp.233-242
• /
• 2008

This study presents a non-prismatic beam element for modeling the elastic and inelastic behavior of steel beams, which have the post-Northridge(cover plate) connections in steel moment frames that are subjected to earthquake ground motions. The elastic stiffness matrix for non-prismatric members with increased beam section (IBS) connection is in the closed-form. The plasticity model is of a discrete type and is composed of a series of nonlinear hinges connected by rigid links. The hardening rules can model the inelastic behavior for monotonic and random cyclic loading, and the effects of local buckling. Moreover the determination of yield surfaces, stiffness parameters, and hardening (or softening) rule parameters for IBS beam element were described. Analytical results of the IBS beam element show good correlation with test data and FEM results.

• Economic and Environmental Geology
• /
• v.53 no.3
• /
• pp.271-285
• /
• 2020

In this study, a numerical experiment related to the stress-strain analysis was performed on 3-D discrete fracture network(DFN) systems based on the distinct element method to evaluate the effect of joint geometry on deformability of jointed rock masses. Using one or two joint sets with deterministic orientation, a total of 12 3-D DFN blocks having 10m cube domain were generated with different joint density and size distribution. Directional deformation modulus of the DFN cube blocks were estimated along the axis directions of 3-D cartesian coordinate. In addition, deviatoric stress directions were chosen at every 30° of trend and plunge in 3-D for some DFN blocks to examine the variability of directional deformation modulus with respect to joint geometry. The directional deformation modulus of the DFN block were found to reduce with the increase of joint size distribution. The increase in joint density was less likely to have a significant effect on directional deformation modulus of the DFN block in case of the effect of rock bridges was relatively large because of short joint size distribution. It, however, was evaluated that the longer the joint size, the increase in the joint density had a more significant effect on the anisotropic deformation modulus of the DFN block. The variation of the anisotropic deformation modulus according to the variations in joint density and size distribution was highly dependent on the number of joint sets and their orientation in the DFN block. Finally, this study addressed a numerical procedure for stress-strain analysis of jointed rock masses considering joint geometry and discussed a methodology for practical application at the field scale.

• The Journal of Engineering Geology
• /
• v.32 no.1
• /
• pp.59-72
• /
• 2022

This study employed a 3-D numerical analysis based on the distinct element method to estimate the strength and deformability of a Cretaceous biotite granitic rock mass at Gijang, Busan, Korea. A workflow was proposed to evaluate the scale effect and the representative elementary volume (REV) of mechanical properties for fractured rock masses. Directional strength and deformability parameters such as block strength, deformation modulus, shear modulus, and bulk modulus were estimated for a discrete fracture network (DFN) in a cubic block the size of the REV. The size of the mechanical REV for fractured rock masses in the study area was determined to be a 15 m cube. The mean block strength and mean deformation modulus of the DFN cube block were found to be 52.8% and 57.7% of the intact rock's strength and Young's modulus, respectively. A constitutive model was derived for the study area that describes the linear-elastic and orthotropic mechanical behavior of the rock mass. The model is expected to help evaluate the stability of tunnels and underground spaces through equivalent continuum analysis.

• Tunnel and Underground Space
• /
• v.10 no.2
• /
• pp.173-179
• /
• 2000

The influence of the plastic flow rules on the elasto-plastic behaviour of a discrete joint element was investigated by performing the numerical direct shear tests under both constant normal displacement and normal displacement conditions. The finite interface elements obeying Plesha’s joint constitutive law was used to allow the relative motion of the rock blocks on the joint surface. Realistic results were obtained in the tests adopting the non-associated flow rule, while the associated flow rule overestimated the joint dilation. To overcome the computational drawbacks coming from the non-symmetric element stiffness matrix in the conventional non-associated plasticity, the symmetric formulation of the tangential stiffness matrix for a non-associated joint element was proposed. The symmetric elasto-plastic matrix it derived by assuming an imaginary equivalent joint with associated flow rule which shows the same plastic response as that of original Joint with non-associated flow rule. The validity of the formulation was confirmed through the numerical direct shear tests under constant normal stress condition.

• Magazine of the Korea Concrete Institute
• /
• v.7 no.4
• /
• pp.149-156
• /
• 1995

The fracture process zone in concrete is a region ahead of a traction-free crack, in which two major mechanisms, microcracking and bridging, play important roles. The toughness due to bridging is dominant compared to toughness induced by microcracking, so that the bridging is dominani: mechanism governing the fracture process of concrete. Fracture mechanics does work for concrete provided that the fracture process zone is being considered, so that the development of model for the fracture process zone is most important to describe fracture phenomena in concrete. In this paper the bridging zone, which is a part of extended rnacrocrack with stresses transmitted by aggregates in concrete, is modelled by a Dugdale-Barenblatt type model with linear tension-softening curve. Two finite element techniques are shown for the analysis of progressive cracking in concrete based on the discrete crack approach: one with crack element, the other without crack element. The advantage of the technique with crack element is that it dees not need to update the mesh topology to follow the progressive cracking. Numerical results by the techniques are demonstrated.

• Magazine of the Korea Concrete Institute
• /
• v.6 no.3
• /
• pp.201-214
• /
• 1994

The fracture process zone in concrete is a region ahead of a traction-free crack, in which two major mechanisms, microcracking and bridging, play important roles. The toughness due to bridging is dominant compared to toughness induced by microcracking, so that the bridging is dominani: mechanism governing the fracture process of concrete. Fracture mechanics does work for concrete provided that the fracture process zone is being considered, so that the development of model for the fracture process zone is most important to describe fracture phenomena in concrete. In this paper the bridging zone, which is a part of extended rnacrocrack with stresses transmitted by aggregates in concrete, is modelled by a Dugdale-Barenblatt type model with linear tension-softening curve. Two finite element techniques are shown for the analysis of progressive cracking in concrete based on the discrete crack approach: one with crack element, the other without crack element. The advantage of the technique with crack element is that it dees not need to update the mesh topology to follow the progressive cracking. Numerical results by the techniques are demonstrated.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.22 no.3
• /
• pp.16-20
• /
• 2018

Steel Fiber reinforced self-compacting concrete (SFRSCC) has been widely used in a number of structures, such as ordinary civil infrastructures, sky scrapers, nuclear power plants, hospitals, dams, channels and etc. Thanks to its short and discrete reinforcing fibers, its performance, including tensile strength, ductility, toughness and flexural strength gets much better in comparison with ordinary self-compacting concrete (SCC) without any reinforcing fibers. Despite all these aforementioned advantages of SFRSCC, its performance highly depends on fiber's orientation. In case of short discrete fibers, the orientation of fibers is completely random and cannot be controlled during pumping process. If fibers distribution inside hardened state concrete are randomly distributed, it leads to less resistance potential of concrete element, especially in terms of flexural and tensile strength. The maximum expected strength may not be achieved. Therefore, fiber alignment has been considered as one of the important factors in SFRSCC. To address this issue, this study investigates the effects of concrete matrix's density and inlet velocity on fiber alignment during the pumping process using a finite element method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Journal of the Korean Mathematical Society
• /
• v.50 no.5
• /
• pp.1129-1163
• /
• 2013

Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.12
• /
• pp.1180-1187
• /
• 2007

In the present work, topology optimization method using adaptive inner-front level set method is presented. In the conventional level set based topology optimization method, there exists an incapability for inner-front creation during optimization process. In this regard, as a new attempt to avoid and to overcome the limitation, an inner-front creation algorithm is proposed. In the inner-front creation algorithm, the strain energy density of a structure along with volume constraint is considered. Especially, to facilitate the inner-front creation process during the optimization process, the inner-front creation map which corresponds to the discrete valued function of strain energy density is constructed. In the evolution of the level set function during the optimization process, the least-squares finite element method (LSFEM) is employed. As an application to shell structures, the lightweight design of doubly curved shell and segmented mirror is carried out.